地中海无花果遗传多样性研究揭示果实性状相关基因

《Frontiers in Plant Science》：Exploring the genetic diversity of Mediterranean fig trees highlights genes associated with fruit traits

本研究通过对地中海地区286份无花果种质进行全基因组重测序（WGR），系统鉴定了超过137万个高质量单核苷酸多态性（SNP）及结构变异，结合11项果实性状的表型数据，通过全基因组关联分析（GWAS）成功定位到481个显著关联位点，发现FMO1、MYB转录因子、ABC转运蛋白等关键基因调控果实重量、大小及品质形成，为无花果种质资源管理、分子标记辅助育种及基因组编辑提供了重要理论依据。

　　

引言

无花果（Ficus carica L.）作为地中海地区具有重要历史和经济价值的多年生作物，其果实、叶片和乳胶在食品、医药和化妆品领域应用广泛。尽管无花果具有显著的农艺和文化意义，但现代育种项目仍较为有限，且主要农艺性状的遗传基础尚未得到充分解析。为填补这一知识空白，本研究对来自西班牙、土耳其和突尼斯种质库的286份无花果基因型进行了全基因组重测序（WGR），旨在全面探索其遗传多样性，并鉴定与重要果实性状相关的遗传位点。

材料与方法

研究材料包括来自西班牙（61份）、土耳其（115份）和突尼斯（110份）种质库的286份无花果品种，其中包含18份雄株和268份雌株。利用Illumina NovaSeq 6000平台对幼叶提取的DNA进行双端150 bp测序，平均测序深度为18.31×。使用GATK等工具进行变异检测，共鉴定出1,374,111个高质量SNP、2,448,766个小插入/缺失（InDel）、218个拷贝数变异（CNV）和1363个结构变异（SV）。通过主成分分析（PCA）、STRUCTURE分析和系统发育树构建等方法解析种群结构，并利用GEMMA软件对11项果实性状进行全基因组关联分析（GWAS）。

结果

基因组变异特征

WGR数据分析显示，SNP密度为每252个碱基1个SNP，其中47.11%位于基因间区，20.85%位于内含子区。编码区共发现96,726个SNP，包括51,808个非同义突变和43,723个同义突变，非同义/同义突变比率（K_a/K_s）为1.2。InDel中以单碱基InDel为主（53.64%）。连锁不平衡（LD）分析表明，r2≤0.2的衰减距离约为18 kb。CNV和SV分别富集于离子转运、蛋白酶活性等生物学过程，提示这些变异可能参与环境适应性调控。

种群结构分析

PCA和STRUCTURE分析均将286份基因型划分为三个主要类群，分别对应西班牙、土耳其和突尼斯的种质来源，但存在部分基因型混杂现象，反映了历史上种质交流的影响。遗传分化指数（FST）显示群体间遗传差异适中（0.0399–0.0658），土耳其种质的核苷酸多样性最高（π=0.00178）。亲缘关系网络进一步揭示了种质间的潜在重复命名或近缘关系，如西班牙的"Clon 300"与"Granito"、土耳其的"Sarilop-1029"与西班牙的"Smyrna"均呈现高度遗传相似性。

性状关联分析

GWAS共鉴定到481个与果实性状显著关联的SNP。果实重量（WE）相关位点中，染色体2上的FMO1基因和染色体6上的MYB转录因子基因分别贡献+5.61 g和+11.36 g的加性效应。果实大小性状中，染色体13的ABC转运蛋白基因与长度（LG）相关，染色体6的WAK激酶基因与宽度（WD）相关。可溶性固形物（TSS）和可滴定酸（TA）分别与CYP94C1基因和PPR蛋白基因关联。果实硬度（FM）相关位点位于CML22钙结合蛋白基因和糖转运蛋白ERD6基因附近。此外，染色体9的PP2C63基因与结果习性（PT）显著相关。

讨论

本研究通过高密度SNP标记和精细表型鉴定，构建了目前最全面的无花果基因组资源库。种群结构分析证实了地中海无花果种质的地理分化特征，同时揭示了种质交换导致的遗传混杂。GWAS发现的候选基因多数在果实发育、糖代谢和激素信号通路中具有已知功能，如FMO1参与生长素合成、WAK激酶调控细胞扩张等，为无花果品质改良提供了关键靶点。然而，部分性状关联位点未直接落在基因编码区，提示可能存在新型调控元件，需进一步通过功能验证明确其机制。

结论

本研究系统解析了地中海无花果的遗传多样性与果实性状的遗传基础，为种质资源管理、分子标记辅助育种和基因组编辑提供了重要数据支撑。发现的候选基因和优异等位变异有望应用于无花果产量、品质和适应性的协同提升，对促进该作物的可持续利用具有重要实践价值。