疾病相关遗传变异可在组织特异性蛋白亚型中引起错义效应

Disease-associated genetic variants can cause missense effects in tissue-specific protein isoforms.

作者信息Giovanna Weykopf, Mihaly Badonyi, Elias T Friman, Jasmine Minh Hang Nguyen, Alexis Ioannou, Benjamin J Livesey, Audrey Coutts, Elizabeth F Hird, Murray Wham, Chloe M Stanton, Veronique Vitart, Jing Su, Lee Murphy, J Kenneth Baillie, Mark D Gorrell, Joseph A Marsh, Wendy A Bickmore, Simon C Biddie
PMID42303611
发布时间2026-06-16
DOI10.1038/s41467-026-74280-w

实验完整度

高

多组学数据分析结合AlphaFold建模、细胞系RNA验证、酶活与免疫共沉淀功能验证。

主要模型

GTEx等长读长RNA-seq转录本与人群变异数据 肺上皮细胞系A549、H358、HSAEC1-KT等 HEK293T细胞过表达DPP9与NLRP1 AlphaFold 3预测的DPP9蛋白结构模型

重点核对

DPP9-AFE变异rs12610495 A>G引入p.Leu8Pro,在HEK293T中共免疫沉淀验证与NLRP1相互作用 DPP9酶活检测使用H-Gly-Pro-pNA(1 mM)及H-Val-Pro-AMC等底物,设三个生物学重复 FLEXIR-seq靶向长读长RNA测序在A549、H358、HSAEC1-KT等细胞系验证fl-DPP9-AFE表达 AlphaFold 3结合FoldX计算p.Leu8Pro的ΔΔG及复合物稳定性 qRT-PCR在A549、H358和HSAEC1-KT细胞中验证DPP9-AFE转录本

摘要

遗传变异可导致蛋白质编码突变并引发疾病。变异通常基于基因的参考转录本进行解读。然而,大多数人类多外显子基因存在可变剪接亚型。我们发现,与可变亚型外显子进化约束降低相一致,其携带的人群变异多于参考亚型外显子,且这些变异更可能引起非同义突变。定位于可变转录本的常见和罕见疾病相关变异可导致相应蛋白亚型中预测结构破坏的氨基酸替换。疾病相关变异所定位的可变转录本表现出高度组织特异性,许多在参考人类基因组中未注释,仅通过长读长RNA测序揭示。作为示例,我们报道了炎症小体调节因子DPP9的一个未注释可变转录本,它特异表达于肺上皮,携带一个与重症COVID-19和肺纤维化相关的常见遗传变异。通过靶向捕获全长转录本的深度RNA测序,我们确认了该未注释DPP9亚型的表达。该DPP9亚型变异在可变第一外显子中引起p.Leu8Pro错义突变,预测会破坏编码的α螺旋,并且我们显示该变异改变了DPP9酶活性。我们的发现强调了在变异解读中考虑可变亚型、其组织特异性表达和全长转录本的重要性,对于揭示常见和罕见疾病中被低估的机制具有重要意义。

实验结论

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

研究问题
疾病相关遗传变异是否可通过编码改变作用于组织特异性可变转录本及其蛋白亚型,从而影响疾病机制。
核心机制
常见变异rs12610495 A>G位于DPP9肺上皮特异性可变第一外显子,引起p.Leu8Pro错义突变,预测破坏N端α螺旋,并改变DPP9对底物和抑制剂的亲和力。
主要证据
长读长RNA-seq和FLEXIR-seq确认fl-DPP9-AFE在肺上皮细胞表达;AlphaFold 3/FoldX预测p.Leu8Pro破坏结构;酶活实验显示对H-Val-Pro和Val-boroPro亲和力改变;共免疫沉淀显示与NLRP1相互作用。
研究意义
提示在遗传变异解读中应考虑可变亚型及其组织特异性表达,有助于罕见病变异临床分类和复杂性状生物学洞察。

研究路径

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

1

外显子分类与变异负荷分析

比较参考外显子与可变外显子的进化保守性和人群变异负荷。

基于GTEx、ENCODE等长读长RNA-seq转录本,将外显子分为ref-exon、EIC和ENIC,统计gnomAD变异频率及编码影响。

2

疾病变异映射与组织特异性评估

评估GWAS和ClinVar变异在可变外显子中的分布及对应亚型的组织特异性。

将GWAS Catalog和ClinVar变异映射至长读长RNA-seq可变转录本,计算Gini系数和表达水平。

3

变异效应与蛋白结构预测

预测错义变异对可变亚型蛋白稳定性和结构的影响。

使用ESM-1v、AlphaFold 3和FoldX对约17,000个可变亚型进行变异效应和ΔΔG预测。

4

DPP9可变亚型表达验证

验证fl-DPP9-AFE在肺上皮细胞中的表达及组织特异性。

开发FLEXIR-seq靶向长读长RNA测序,并在肺上皮细胞系和qRT-PCR中验证DPP9-AFE转录本。

5

p.Leu8Pro功能与结构验证

验证p.Leu8Pro对DPP9酶活、抑制剂敏感性和NLRP1相互作用的影响。

在A549细胞中表达DPP9亚型,检测酶活和动力学;在HEK293T中共免疫沉淀DPP9与NLRP1。

研究方法

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

产品清单

实验环节名称品牌货号
DMEM培养基Life Technologies41965039
RPMI培养基----
McCoy’s 5A培养基Gibco26600023
SABM BulletKit培养基LonzaCC-3119 and CC-4124
角质形成细胞无血清培养基Gibco17005042
RNEasy迷你试剂盒Qiagen74104
SuperScript第一链合成系统Invitrogen11904-018
LightCycler 480 SYBR Green I预混液Roche04887352001
CFX96 Touch实时PCR检测系统Bio-Rad--
ONT cDNA-PCR条形码试剂盒Oxford Nanopore TechnologiesSQK-PCS114.24
Twist机械片段化文库制备试剂盒Twist101280
Twist通用接头系统Twist101307
Twist标准杂交和洗涤试剂盒v2Twist104445
Twist通用封闭剂Twist100578
Dynabeads M-270链霉亲和素磁珠Invitrogen65305
ONT连接试剂盒Oxford Nanopore TechnologiesSQK-LSK114
AMPure XP磁珠Beckman Coulter AgencourtA63881
Lipofectamine 3000转染试剂InvitrogenL3000015
PEI转染试剂Sigma Aldrich919012
ChromoTek GFP-Trap磁性琼脂糖珠Proteintechgtma
α-GFP抗体ProteintechPABG1
α-Flag抗体SigmaF1804
α-GAPDH抗体Proteintech10494-1-AP
H-Gly-Pro-pNA底物Bachem4025614
H-Val-Pro-AMC底物Enzo Life SciencesBML-P448
H-Trp-Pro-AMC底物EnzoBML-P450
H-Lys-Pro-AMC底物MedChemExpressP4426A
西格列汀----
Val-boroPro----

关键环节

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

环节核对要点
外显子分类与变异负荷分析
外显子分类标准:ref-exon基于MANE Select,EIC基于Gencode V46,ENIC基于长读长RNA-seq;变异频率分组阈值需核对。
阅读提示:Methods中“Defining exons classes”和“Exon class variant annotation”小节。
疾病变异映射与组织特异性评估
GWAS Catalog和ClinVar版本、Gini系数计算方法、GTEx长读长RNA-seq组织数(22种)及表达过滤条件。
阅读提示:Methods中“Tissue-specific scoring”和“Variant effect prediction in isoforms”小节。
DPP9可变亚型表达验证
FLEXIR-seq探针设计(327条,120 bp)、细胞系来源、qRT-PCR引物和参考基因(GAPDH、SAFB)。
阅读提示:Methods中“Full-Length targEted capture using eXon probes for IsofoRms (FLEXIR-seq)”和“Quantitative RT-PCR”小节。
p.Leu8Pro功能验证
酶活检测底物浓度(H-Gly-Pro-pNA 1 mM,H-Val-Pro-AMC 3 µM)、Val-boroPro浓度范围(0.01–200 nM)、细胞裂解液蛋白量(2.5 µg)及重复数(三个生物学重复)。
阅读提示:Methods中“DPP9 enzymatic and inhibition assays”小节及图5图注。
免疫共沉淀验证
HEK293T转染质粒量(DPP9 1 µg,NLRP1 1.37 µg)、PEI:DNA比例(3:1)、抗体稀释比(α-GFP 1:1000,α-Flag 1:750)。
阅读提示:Methods中“Co-Immunoprecipitation assays and immunoblotting”小节及图5G图注。