基于DGRP的黑腹果蝇基础心率遗传结构及其膳食盐调控：基因型-环境互作与心脏特异性RNAi验证

《PLOS Genetics》：Genetic architecture of basal heart rate and its modulation by dietary salt in Drosophila melanogaster

　　

心率（heart rate, HR）是一种复杂的数量表型，受遗传和环境因素共同影响，并且是心脏生理调控中的关键因素。尽管在人类和其他生物中已发现众多与HR变异相关的基因，但这些基因如何相互作用以影响自然环境中的HR变化及其对环境胁迫的响应，仍未完全阐明。在本研究中，研究人员利用黑腹果蝇遗传参考面板（Drosophila Genetic Reference Panel, DGRP）检查了基础HR及其对盐胁迫响应的遗传基础。154个DGRP品系在标准培养基和含0.1% NaCl培养基条件下，分别对雌性和雄性第三龄幼虫进行HR测定，表现出中等广义遗传力（broad-sense heritability, H²）和显著遗传变异。混合效应模型（mixed-effects model）显示，基因型与环境之间的互作显著影响HR。全基因组关联分析（genome-wide association analysis, GWAS）鉴定出与离子转运、心脏发育、蛋白质稳态、代谢和应激响应等关键生物学过程相关的基因变异。在标准条件和盐水补充条件下鉴定出不同基因，提示环境胁迫释放了隐秘遗传变异（cryptic genetic variation, CGV）。研究人员利用心脏特异性GAL4驱动子在RNA干扰（RNA interference, RNAi）实验中检测候选基因，发现相当数量这些基因依据性别和环境显著影响HR。这些发现揭示，HR变异受到基因型-环境互作（genotype-by-environment interaction, G×E）和性别特异性遗传效应的强烈塑造，并且在基础和胁迫条件下有不同的基因集合发挥作用。环境依赖性遗传结构的鉴定表明，盐胁迫可暴露影响心脏功能的隐秘遗传变异。通过整合全基因组关联分析与基于RNAi的靶向功能验证，该研究提供了一个将自然遗传变异与生理表型以情境依赖方式联系起来的框架。

心血管疾病是全球主要死亡原因之一，HR作为心脏性能、自主平衡、离子通道活动和血管张力的综合指标，是反映机体生理状态的早期生物标志物。静息HR升高会增加心肌氧需求并缩短舒张充盈时间，与心血管疾病和全因死亡风险升高相关，但其临床效应可因遗传背景和疾病情境而异。HR是多基因性状，尽管全基因组研究已发现大量相关基因，但其全基因组遗传结构、环境胁迫下的调控方式以及何时成为可靠心血管风险预测指标仍不清楚。果蝇心脏发育、收缩机制和离子通道功能与人类直系同源物在结构和功能上相似，约75%人类疾病相关基因在果蝇中有同源物，UAS/GAL4和RNAi工具使其适合研究复杂表型的遗传基础。高盐饮食与血压升高及离子平衡紊乱相关，但膳食盐对HR遗传决定因素的影响，尤其在自然遗传变异panel中进行全基因组筛选，尚未见报道。因此，研究人员使用154个DGRP自交系，在标准培养基和0.1% NaCl培养基中饲养第三龄幼虫（third larval stage, L3），雌雄分开，研究基础HR及其盐响应的遗传结构，并结合GWAS和心脏特异性RNAi进行功能验证。该研究揭示HR受G×E和性别特异性遗传效应强烈塑造，盐胁迫可暴露CGV，并提示离子稳态、蛋白质稳态、代谢和应激信号等进化保守网络共同调控心脏功能。

关键技术方法概括：研究以154个DGRP自交系为样本队列（DGRP lines由Trudy F. C. Mackay提供），设置标准培养基和0.1% NaCl培养基，取L3幼虫，雌雄分开，每品系×处理×性别20头，以20秒视频中收缩数计数HR。统计分析采用方差分析、混合效应模型、Welch’s t检验和Benjamini–Hochberg错误发现率（false discovery rate, FDR）校正，估计H²和G×E。GWAS基于品系均值，采用PLINK初筛和FaST-LMM线性混合模型（linear mixed model, LMM），纳入基因组亲缘关系矩阵（genomic relatedness matrix, GRM），过滤条件包括次等位基因频率（minor allele frequency, MAF）≥0.05、call rate≥0.80、最小等位基因计数≥4，阈值p<10^-5；使用PANGEA进行基因集富集分析（gene set enrichment analysis, GSEA）。功能验证采用UAS/GAL4系统，RNAi品系来自维也纳果蝇资源中心（Vienna Drosophila Resource Center, VDRC）KK和GD文库，使用Mef2-GAL4和Hand-GAL4驱动，对照为VDRC60100和VDRC60000，并用VDRC60101和独立GD RNAi评估40D插入相关效应。

研究结果：在“Quantitative genetic variation in heart rate”中，研究人员对154个DGRP品系在标准与盐补充条件下的HR进行分析。标准条件下雌性最低和最高平均HR为51.4和143.7次/20秒，雄性为50.2和140.9次/20秒；盐条件下雌性为50.2和185.0次/20秒，雄性为52.1和188.1次/20秒。各处理和性别的方差分析显示line项均高度显著（p<0.001），表明HR存在显著遗传变异。标准与盐条件下幼虫HR具有中等广义遗传力，雌性H²=0.45，雄性H²=0.47；盐胁迫下遗传组分相对增加至H²=0.51。混合效应模型显示treatment和sex固定效应高度显著（p<0.001），treatment×sex互作较弱但显著（p=0.046），line×treatment互作高度显著并解释相当部分方差。跨环境相关分析显示两性r=0.65（p<0.001），偏离1，支持G×E。Welch’s t检验显示雌性60个品系HR在两种处理间显著差异，其中48个升高、12个降低；雄性46个品系显著，其中37个升高、9个降低。结论是幼虫HR为具中等广义遗传力的数量性状，盐响应具有品系特异性异质性。

在“GWA analysis of heart rate and gene set enrichment analyses”中，研究人员对雄性、雌性、性别平均值及性别差异进行GWAS。控制培养基中鉴定44个候选变异和35个候选基因，盐补充培养基中鉴定56个候选变异和39个基因，两种培养条件HR差异分析鉴定34个候选变异和30个基因。三项GWAS共鉴定134个候选变异，其中101个独特关联到基因区域，包括101个单核苷酸多态性（single nucleotide polymorphism, SNP）、8个缺失、2个插入和1个多核苷酸多态性（multiple nucleotide polymorphism, MNP）。变异主要位于内含子（91个）和下游（22个），另有上游11个、3’UTR 7个、5’UTR 2个、外显子6个；功能上6个非同义编码、11个同义编码，提示非编码和调控区可能贡献HR变异。Q-Q图提示部分关联变异可能为生物信号。GSEA显示名义富集于心脏和肌肉发育、离子跨膜转运及调控心脏稳态的信号通路，基因-疾病关联提示高血压、心力衰竭和QT间期异常。结论是标准与盐条件下候选基因不同，提示环境胁迫释放CGV，且候选基因需结合统计证据和功能背景优先验证。

在“Functional assessment”中，研究人员用UAS/GAL4系统在体内验证候选基因，RNAi品系来自VDRC KK和GD文库，分别使用Mef2-GAL4和Hand-GAL4。Mef2驱动下，标准培养基雌性有20个候选基因经两种统计检验显著，占分析29个基因的69%，其中16个敲低后HR升高、4个降低；盐培养基雌性也有20个基因显著。环境特异效应表现为sns在标准培养基升高HR而盐下无显著变化，MSI在标准无影响而盐下显著升高HR。Mef2驱动雄性标准培养基有15个基因显著，其中11个升高、4个降低，CG15144、CG5375、RluA-1、ringer、sns显示性别特异效应；盐下雄性有19个基因显著，CG15144、RluA-1、ringer、MSI显示环境依赖效应。14个在两种性别和两种统计检验中经Mef2验证的候选基因进一步用Hand驱动检测。雌性标准培养基中7个基因显著：Adgf-A、CG42541、CG6656、EMC3、Mgat4a、bma和nAChRα3，其中5个升高，CG6656和EMC3降低；雄性中7个基因显著，基因集与雌性基本相同，除一个候选不同，CG6656和EMC3降低，其余5个升高。SLO2和shn与Mef2驱动杂交不产生可存活后代，但与Hand驱动杂交在两种培养基中均显著改变HR。部分RNAi-驱动组合出现不能羽化或无存活后代，但未作为正式生存实验分析。结论是大量候选基因功能影响HR，且效应依赖性别、环境和细胞类型。

在“Assessment of potential 40D insertion-associated effects”中，鉴于部分VDRC KK RNAi品系携带40D插入可能异位表达邻近基因tiptop (tio)，研究人员用VDRC60101对照和独立GD RNAi重测相关候选基因。Mef2驱动背景下，GD文库中所有重测候选相对于对照仍显著，但仅EMC3保持相同效应方向，其余候选方向相反。Hand驱动背景下，相对于独立KK对照60101仅部分候选保留显著性，而GD RNAi支持大多数候选显著。结论是候选基因效应重现性依赖RNAi遗传背景和GAL4驱动，解释时需同时考虑显著性和表型方向。

讨论部分总结：研究支持HR为多基因复杂性状，自然遗传变异和G×E共同塑造基础HR及盐响应。标准与盐条件下GWAS候选不重叠、盐下H²略增，提示CGV释放。RNAi验证将候选基因归入离子稳态、蛋白质稳态、代谢与应激响应等互连模块。离子调节部分，nAChRα3、CG42541、Mgat4a、bma和SLO2与离子平衡和电活动相关，敲低后多升高HR，其中CG42541效应最强，增加52次/20秒，人类同源基因涉及钙通道调控和长QT综合征等。蛋白质稳态部分，EMC3敲低降低约50次/20秒，CG6656在Mef2和Hand驱动中方向相反，Prosalpha6T仅Mef2驱动显著，涉及内质网膜复合体、溶酶体质量和26S蛋白酶体。代谢与应激部分，ktub经Mef2敲低降低HR并关联GPCR信号，Adgf-A经两种驱动敲低升高HR并关联腺苷代谢，Adgf-A2（MSI）仅在盐和Mef2条件下显著，CG8635和Egfr关联NF-κB应激信号和EGFR/RTK心脏维持。研究强调这些模块并非独立，而是形成系统级网络。研究结论部分翻译如下：综上，研究人员的发现支持心脏调控的系统级模型，其中多个遗传模块，包括离子稳态（ion homeostasis）、蛋白质稳态（proteostasis）、代谢（metabolism）和应激信号（stress signaling），相互作用以决定心脏性能。这些模块并非独立作用，而是形成调控心肌细胞功能和生理输出的整合网络。在此背景下，HR呈现为由这些互连通路协调互作而非孤立基因效应塑造的复杂表型。研究人员的发现表明，心脏调控由果蝇和人类共有的进化保守遗传网络编排。该整合框架与物种间核心心脏调控机制的进化保守性一致。这种保守性在兴奋-收缩偶联（excitation–contraction coupling）、钙信号（calcium signaling）和应激响应通路等核心过程中尤为明显，这些过程通过果蝇和人类心肌细胞中的同源基因和调控网络运作。因此，即使这些网络发生细微改变，也可导致深刻的功能损害，有效重现人类心脏病理的遗传基础。本研究中基于RNAi的方法提供了功能证据，支持多个候选基因参与HR调控。未来整合蛋白质互作网络、组织特异性表达分析和细胞类型特异性方法的研究，将进一步增强对观察到的表型效应分子机制的理解。应开展额外研究，在脊椎动物模型中研究这些遗传因素调控心律的作用，并结合钙动态分析以及人源心肌细胞的转化研究。本研究强调，通过包含自然变异和基因型×环境互作（G×E）的广泛遗传背景分析来理解HR遗传结构的重要性。在此框架内，果蝇模型是理解这一复杂系统多个层面的优秀实验工具。