《Current Issues in Molecular Biology》：Uncovering Functional Genetic Variation in the Indigenous Greek Eghoria Goat: An Integrative Genome-Wide Discovery and Validation Approach
编辑推荐：
艾戈利亚山羊是希腊主要的本土山羊群体，代表了一种重要的遗传资源。尽管其具有重要性，但基因组信息仍然匮乏。本研究旨在利用全基因组测序（WGS）调查艾戈利亚山羊的编码区变异，特别关注错义单核苷酸多态性（SNPs）和编码插入缺失（indels）。研究人员对艾戈利亚山
艾戈利亚山羊是希腊主要的本土山羊群体，代表了一种重要的遗传资源。尽管其具有重要性，但基因组信息仍然匮乏。本研究旨在利用全基因组测序（WGS）调查艾戈利亚山羊的编码区变异，特别关注错义单核苷酸多态性（SNPs）和编码插入缺失（indels）。研究人员对艾戈利亚山羊（n=6）进行了全基因组测序，随后进行生物信息学处理和变异过滤。通过基因本体论（GO）和京都基因与基因组百科全书（KEGG）通路分析（FDR<0.05）进行功能表征。归一化变异密度指标识别出高多态性基因。利用MassARRAY基因分型技术在一个独立群体（n=54）中对选定的错义SNPs进行了验证。全基因组分析鉴定出10,796,211个SNPs和1,022,779个indels。经过优先化筛选后，保留了15,949个错义SNPs和861个indels。功能富集分析突出了与代谢、离子转运、钙信号传导、免疫功能、转录调控和环境适应相关的通路。表现出升高的多态性密度的基因包括嗅觉受体家族成员以及与代谢调控相关的位点。验证分析确认了选定的变异在艾戈利亚山羊群体中的存在。本研究扩展了当前对本土艾戈利亚山羊基因组多样性的认知，并为未来的遗传改良、保护和育种计划提供了宝贵的资源。
**论文解读：揭示希腊本土艾戈利亚山羊功能性遗传变异——整合全基因组发现与验证**
**一、研究背景与目的**
全球畜牧业中，山羊因其对恶劣环境的适应能力及肉奶贡献而占据重要地位。希腊拥有欧盟最大的山羊群体，其中本土艾戈利亚（Eghoria）山羊占全国羊群的约90%，呈现广泛的表型多样性（如毛色、长毛），能在低投入、半粗放系统下高效利用低质草场并耐受炎热干燥气候，是重要的遗传资源。然而，本土山羊正面临种群下降和遗传侵蚀，主因是生产集约化以及与高产品种（如萨能、阿尔卑斯）的无序杂交。尽管基因组研究对保护与可持续利用至关重要，但针对艾戈利亚山羊的高分辨率全基因组研究仍然稀缺，以往研究多基于低深度测序或缺乏系统功能优先化与实验验证。因此，本研究旨在利用全基因组测序（WGS）全面表征艾戈利亚山羊的编码序列变异，重点关注错义单核苷酸多态性（SNPs）和外显子插入缺失（indels），并通过功能注释、富集分析及独立群体验证，揭示其适应性与经济性状相关的遗传基础。该论文发表在《Current Issues in Molecular Biology》。
**二、主要关键技术方法**
研究人员采用以下主要技术方法：1. 对来自希腊色萨利（Thessaly）地区的6只艾戈利亚山羊进行全基因组测序（Illumina NovaSeq X Plus平台，PE150，约22×覆盖度），采用等摩尔混合池化策略；2. 使用SAMtools/BCFtools流程进行变异检测，并用ANNOVAR进行功能注释，优先筛选出外显子区错义SNPs和编码indels；3. 利用ShinyGO v0.82进行基因本体论（GO）和KEGG通路富集分析（FDR<0.05）；4. 计算归一化变异密度（SNPs/kb和indels/kb）鉴定高多态性基因；5. 在独立群体54只艾戈利亚山羊中，使用MassARRAY系统（MALDI-TOF质谱）对24个候选错义SNPs进行个体水平基因分型验证，并计算等位基因频率、期望/观测杂合度（He/Ho）、多态信息含量（PIC）及哈代-温伯格平衡（HWE）检验。
**三、研究结果**
**3.1 全基因组测序与数据质量**：测序获得约60.8 Gb原始数据，经质控后保留60.0 Gb高质量数据，Q20>97%，Q30>93%，测序错误率约0.01%，GC含量约44%。与参考基因组（ARS1）比对率达99.53%，为可靠变异检测提供了充足深度。
**3.2 变异鉴定、优先化与基因组分布**：共鉴定出10,796,211个SNPs和1,022,779个indels，大部分位于基因间区和内含子区。经优先化筛选后，保留15,949个错义SNPs和861个外显子indels。归一化变异密度分析显示，SNP密度最高的基因包括SPACA4、OR10C1、OR5M3、OR1Q1、OR6K2等，其中多数属于嗅觉受体（OR）基因家族；indel密度最高的基因包括TMEM250、LDOC1、SPRN、ELFN1、OR5M11等，功能涉及调控、神经和免疫相关过程。
**3.3 功能富集分析**：对携带错义SNPs的基因进行富集分析，在生物学过程（BP）层面主要富集于刺激响应、发育调控和细胞组织化；分子功能（MF）层面以离子通道和跨膜转运活性为主；细胞组分（CC）层面涉及膜区和细胞外基质。KEGG通路显著富集于ECM-受体相互作用、钙信号通路、ABC转运体、味觉转导、补体与凝血级联等。对外显子indels的富集分析显示，BP层面富集于发育和系统过程、负调控及细胞组分组织；MF层面以钙依赖性磷酸酶活性、GTP酶调节活性和转录因子结合为主；CC层面涉及核质和膜结构，但未发现显著富集的KEGG通路。
**3.4 候选错义SNPs的MassARRAY验证**：在54只独立群体中成功对24个候选错义SNPs进行基因分型。样本水平基因分型成功率达90–100%，SNP水平多数达到≥95%的检出率。等位基因频率范围为0.54–0.88，均属中等至高频率，未发现稀有变异。期望杂合度（He）为0.211–0.497，观测杂合度（Ho）为0.184–0.482，多态信息含量（PIC）为0.189–0.373，所有位点均未显著偏离哈代-温伯格平衡（p>0.05）。SIFT预测显示三个变异（位于OR1Q1、OR6F1和GPR35）为有害，其余为耐受；GERP保守性评分普遍较低。
**四、讨论与结论**
**讨论总结**：本研究发现艾戈利亚山羊的编码区变异主要富集于信号转导、发育调控、细胞组织化和代谢等相关通路，与本土品种的环境适应能力相吻合。错义SNPs相关的富集通路（如钙信号、ABC转运体、ECM-受体相互作用、补体级联、味觉转导、赖氨酸降解）提示其在离子稳态、应激响应、营养利用和免疫功能中的潜在作用。高变异性基因如嗅觉受体基因家族（OR10C1等）可能与觅食行为和环境感知有关。此外，高通量变异密度也出现在免疫相关基因（S100A12、ACKR2）、代谢调控基因（FGFBP1、HTR1D）中。外显子indels则富集于转录调控、发育和细胞组织化过程，提示其对基因表达调控的潜在影响。验证的候选基因（如SIRT1、GYS2、SI、LIPC、CACNA1C、MYLK）与代谢、钙信号、生长和产奶性状相关，进一步支持了变异的功能相关性。局限性包括：采用混合池测序无法获得个体基因型、难以检测低频变异；发现阶段样本量较小；缺乏基因型-表型关联分析；未进行与其他品种的比较基因组分析。因此，所鉴定的变异应视为候选标记而非确认的功能变异。
**结论翻译**：总之，本研究提供了对希腊本土艾戈利亚山羊的首批全基因组特征描述之一，整合了变异发现、功能注释和实验验证。所鉴定的编码变异扩展了现有该本土品种的基因组资源，并证明了所应用的全基因组测序和验证策略的可靠性。这些发现为本土山羊群体的保护、遗传改良和可持续利用提供了有价值的基因组基础。未来整合更大群体、详细表型记录和功能验证的研究将有助于鉴定与适应性和经济重要性状相关的遗传标记，从而支持基因组辅助育种和保护计划的实施。
