结合全长基因剪接分析和SpliceAI预测解读SPINK1基因所有可能编码区单核苷酸变异对剪接的影响

Combining full-length gene assay and SpliceAI to interpret the splicing impact of all possible SPINK1 coding variants.

作者信息Hao Wu, Jin-Huan Lin, Xin-Ying Tang, Gaëlle Marenne, Wen-Bin Zou, Sacha Schutz, Emmanuelle Masson, Emmanuelle Génin, Yann Fichou, Gerald Le Gac, Claude Férec, Zhuan Liao, Jian-Min Chen
PMID38414044
发布时间2024-02-27
DOI10.1186/s40246-024-00586-9

实验完整度

研究包含回顾性与前瞻性FLGSA实验验证、SpliceAI预测、数据外推及进一步验证,且结果包含体外剪接分析和功能验证。

主要模型

HEK293T细胞 全长SPINK1基因表达载体(pcDNA3.1/V5-His-TOPO)转染模型

重点核对

SPINK1基因NM_001379610.1参考序列 SpliceAI预测的Δ评分(回顾性:Illumina预计算评分;前瞻性:SpliceAI Lookup) FLGSA方法使用pCDNA3.1/V5-His-TOPO载体 HEK293T细胞培养与转染条件(DMEM基础培养基,10%胎牛血清,3.5×10^5细胞/孔,2.5μg质粒) RT-PCR产物Sanger测序验证

摘要

背景:基因编码序列中的单核苷酸变异(SNVs)可显著影响前体mRNA剪接,对致病机制和精准医学具有重要意义。本研究旨在利用成熟的全长基因剪接分析(FLGSA)结合SpliceAI预测,前瞻性解读与慢性胰腺炎相关的SPINK1基因(含四个外显子)所有潜在编码SNVs的剪接效应。结果:本研究首先对先前已通过FLGSA评估的27个SPINK1编码SNVs进行回顾性分析,随后对35个新进行FLGSA检测的SPINK1编码SNVs进行前瞻性分析,接着进行数据外推,最后进行进一步验证。我们共分析了67个SPINK1编码SNVs,占720个可能编码SNVs的9.3%。在这67个经FLGSA分析的SNVs中,发现12个影响剪接。通过详细比较FLGSA结果和SpliceAI预测,我们推断SPINK1基因中其余653个未检测的编码SNVs不太可能显著影响剪接。在12个影响剪接的事件中,九个同时产生正常剪接和异常剪接的转录本,而其余三个仅产生异常剪接的转录本。这些影响剪接的SNVs仅出现在外显子1和2中,特别是在这些外显子的第一个和/或最后一个编码核苷酸上。在12个影响剪接的事件中,11个为错义变异(占506个潜在错义变异的2.17%),一个为同义变异(占164个潜在同义变异的0.61%)。值得注意的是,将SpliceAI的截断值调整为0.30而非常规的0.20,可在不降低敏感性的情况下提高特异性。结论:通过整合FLGSA和SpliceAI,我们确定SPINK1中所有可能编码SNVs中只有不到2%(1.67%)显著影响剪接结果。我们的发现强调了在研究基因的更广泛基因组序列背景下进行剪接分析的重要性,并突显了中间SpliceAI评分(0.20至0.80)固有的不确定性。本研究首次以高准确度前瞻性解读疾病相关基因中所有可能编码SNVs,代表了在外显子组和基因组测序时代从回顾性向前瞻性变异分析转变的有意义尝试。补充信息:在线版本包含补充材料,可在10.1186/s40246-024-00586-9获取。关键词:慢性胰腺炎,全长基因剪接分析(FLGSA),遗传学精准医学,前体mRNA剪接,单核苷酸变异(SNVs),SpliceAI,剪接预测算法,剪接位点,SPINK1基因,变异解读。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
SPINK1基因所有可能的编码SNVs中有多少会影响剪接,以及FLGSA和SpliceAI预测的关联性如何?
核心机制
SPINK1编码SNVs通过影响剪接位点(如供体或受体位点)影响前体mRNA剪接,导致异常剪接转录本的产生。
主要证据
通过FLGSA分析67个SPINK1编码SNVs,发现12个影响剪接,其中9个产生正常与异常剪接混合转录本,3个仅产生异常剪接;结合SpliceAI预测,推断剩余未分析变异不显著影响剪接。
研究意义
本研究首次前瞻性解读所有潜在编码SNVs,强调在基因组背景下进行剪接分析的重要性,为从回顾性向前瞻性变异分析转变提供了潜在框架。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

回顾性关联分析

分析已知SPINK1编码SNVs的FLGSA结果与SpliceAI预测的一致性。

收集27个已知SPINK1编码SNVs的FLGSA数据,与Illumina预计算SpliceAI Δ评分进行比较。

2

选择新SPINK1编码SNVs进行FLGSA

基于SpliceAI预测评分选择潜在影响剪接的变异进行验证。

根据Illumina预计算SpliceAI评分,针对每个外显子选择候选变异,包括所有外显子起始和终止位置的SNVs,以及Δ评分≥0.20的SNVs,最终选择35个变异进行FLGSA分析。

3

FLGSA剪接分析

检测新选择的35个SPINK1编码SNVs对剪接的影响。

将变异引入全长SPINK1基因表达载体,经转染HEK293T细胞,提取RNA,RT-PCR扩增,Sanger测序分析剪接产物。

4

FLGSA数据与SpliceAI预测的关联分析

比较新分析的35个SNVs的FLGSA结果与SpliceAI预测评分,并推断未分析变异的影响。

比较FLGSA结果与对应SpliceAI评分,发现Δ评分>0.80的变异产生异常剪接,Δ评分<0.30的变异仅产生正常剪接;基于外显子内比较,推断中间Δ评分(0.20-0.80)变异的影响不确定,但通过外推认为未分析变异不太可能影响剪接。

5

数据外推

推断所有未分析的SPINK1编码SNVs的剪接影响。

基于67个变异(27个已知+35个新+5个验证)的FLGSA数据与SpliceAI预测的关联,推断剩余的653个未分析变异不太可能显著影响剪接。

6

进一步验证

验证外推结论的正确性。

选择5个未分析但预测评分较高的变异进行FLGSA分析,结果显示均仅产生正常剪接转录本,验证了外推。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

验证样本或模型特征
验证机制链路

产品清单

实验环节名称品牌货号
DMEM基础培养基Gibco--
胎牛血清Procell--
HieffTrans通用转染试剂Yeasen--
FastPure细胞/组织总RNA提取试剂盒V2(含gDNA清除剂)Vazyme--
HiScript III第一链cDNA合成试剂盒Vazyme--
2×Taq预混液Vazyme--
凝胶回收试剂盒Omega Bio-Tek--
BigDye Terminator v3.1循环测序试剂盒Applied Biosystems--
pcDNA3.1/V5-His-TOPO载体----

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
细胞培养与转染
HEK293T细胞系、培养基(DMEM)、血清(10%胎牛血清)、细胞密度(3.5×10^5细胞/孔)、转染试剂(HieffTrans)、质粒量(2.5μg)
阅读提示:Methods: Cell culture, transfection, RNA extraction, and reverse transcription (RT)
RNA提取与逆转录
转染时间(48小时)、RNA提取试剂盒、RT试剂盒、总RNA量(1μg)
阅读提示:Methods: Cell culture, transfection, RNA extraction, and reverse transcription (RT)
RT-PCR及测序
PCR反应体系(12.5μL 2×Taq Master Mix, 1μL cDNA, 0.4μM引物)、引物序列、PCR循环参数(94°C 5min; 35 cycles 94°C 30s, 55°C 30s, 72°C 5min; 72°C 7min)
阅读提示:Methods: RT-PCR and sequencing of the resulting products
FLGSA变异选择
SPINK1参考序列(NM_001379610.1)、SpliceAI评分阈值(0.20和0.80)、外显子起始和终止位置SNVs的选择
阅读提示:Methods: Selection of potential SPINK1 coding variants for FLGSA & Results: Selection of potential SPINK1 coding SNVs for FLGSA
数据分析与推断
SpliceAI评分类型(AG, AL, DG, DL)、最高Δ评分的使用、阈值调整(0.30)
阅读提示:Results: Overall correlation between FLGSA findings and SpliceAI predicted scores across all exons