《Plants》：Development of Genome-Derived InDel Markers and Genetic Diversity Analysis of Caragana acanthophylla in Xinjiang, China
编辑推荐：
Caragana acanthophylla Kom.是中国新疆一种生态上重要的耐旱灌木，但用于种质表征的物种特异性分子标记仍较有限。研究人员在新疆该物种已知分布区内从11个地点采样93个个体，每地点选3个共33个个体进行全基因组重测序，获得2,873,410
Caragana acanthophylla Kom.是中国新疆一种生态上重要的耐旱灌木，但用于种质表征的物种特异性分子标记仍较有限。研究人员在新疆该物种已知分布区内从11个地点采样93个个体，每地点选3个共33个个体进行全基因组重测序，获得2,873,410个高质量SNP（single nucleotide polymorphism）和5,679,915个InDel（insertion/deletion）。基于全基因组的SNP（single nucleotide polymorphism）PCA（principal component analysis）与遗传关系分析对33个重测序个体给出独立高分辨率评估。从34对候选引物中保留8个扩增稳定、分型清晰的多态InDel标记并用于全部93个个体。SNP数据集显示明显区域分化和更细的地点关联关系；同一33个个体用8个InDel位点可恢复部分广泛模式，尤其西部YL材料分化，但细尺度分辨率较低。全部93个个体中，InDel面板显示中低标记水平遗传多样性和可检测区域分化。AMOVA（analysis of molecular variance）中11个原始采样地点间变异占67.00%，5个探索分析组间变异占68.37%。Mantel检验在全93个个体中r=0.801、p<0.001，排除YL后r=?0.032、p=0.724，表明整体空间信号主要由YL地理隔离驱动。结果表明8标记面板是初步种质表征和扩大样本筛查的实用低成本工具，而全基因组SNP数据对种群水平推断分辨率更高。
研究背景方面，Caragana属隶属Fabaceae，是干旱与半干旱区重要灌木类群。Caragana acanthophylla长期适应高温、干旱、盐碱及贫瘠生境，具有强抗逆性与环境适应性，是水土保持、防风固沙、土壤改良及区域生态修复的重要种质资源。已有研究多集中于生态分布、种子萌发、苗木培育与资源利用，而该物种在新疆已知分布范围内的种群水平遗传信息仍有限。尽管同属其他种研究表明地理隔离、环境异质性与历史过程可促成 arid 区种群分化，但C. acanthophylla缺乏物种特异性分子标记与群体遗传学证据，限制了种质评价、保护取样和种源管理。
研究人员以新疆11个地点93个天然个体为材料，其中每地点3个共33个个体进行全基因组重测序，并基于Caragana arborescens染色体级参考基因组筛选InDel（insertion/deletion，插入缺失）变异、设计PCR（polymerase chain reaction）引物，最终建立8个多态InDel标记并推广至全部93个个体；同时利用全基因组SNP（single nucleotide polymorphism，单核苷酸多态性）对33个重测序个体作独立遗传关系评估。论文发表于《Plants》。
关键技术方法上，研究人员在新疆11个地点采集93个C. acanthophylla天然个体，地点间距大于10 km、个体间距大于10 m；选取33个个体开展全基因组重测序，比对C. arborescens参考基因组并调用SNP与InDel；依InDel长度、染色体分布、侧翼保守性和片段大小差异初筛34对引物，经PCR扩增与凝胶分型保留8个多态位点；以GenAlEx进行遗传多样性、PCoA（principal coordinate analysis）、AMOVA与Mantel检验，以MEGA构建UPGMA树，以STRUCTURE推断群体结构，以GCTA和TreeBeST对33个重测序个体做SNP水平PCA与NJ（neighbor-joining）树分析。
研究结果如下。2.1 全基因组候选InDel位点发现与分布：33个个体重测序获得高质量数据，比对率、覆盖度与深度均较高；共鉴定5,679,915个InDel，总长达30,018,850 bp，占参考基因组2.4687%，删除型数量多于插入型但插入型总更长，染色体间分布不均，说明候选位点资源充足。2.2 InDel引物多态性分析：34对候选引物保留8个；等位基因数Na均值1.950，有效等位基因数Ne、Shannon指数I、观测杂合度Ho、期望杂合度He与PIC（polymorphism information content）差异明显，部分位点Fis极端，部分偏离HWE（Hardy–Weinberg equilibrium），Fst与Nm显示位点间分化与基因流估计差异大。2.3 33个重测序个体全基因组SNP遗传关系：PCA中PC1与PC2解释31.79%和5.88%变异，YL沿PC1明显分离，非YL材料也有分化；NJ树显示地点关联聚类更强，说明SNP分辨率高于8个InDel。2.4 93个个体InDel遗传多样性：五组分析中Na、Ne、I、Ho、He、uHe总体中低；AS的I最高，YL最低，表明YL遗传变异较低而AS相对较高。2.5 93个个体InDel遗传关系：PCoA第一、二坐标解释58.51%和7.41%，YL居一侧；STRUCTURE在K=2时YL与其他四组祖先成分差异最明显，但非YL组重叠与混合较强，不支持锐利离散种群。2.6 UPGMA聚类：材料按遗传距离与地理来源聚类，YL成独立大支；WS与CJ内部相似，JM介于其间并与邻组重叠，AS近CJ但不完全嵌套，故保留为独立探索组。2.7 AMOVA：11地点间变异67.00%，五组间68.37%，说明空间分化显著，可能与断续分布、山地—盆地景观及遗传漂变有关。2.8 地理与遗传距离关系：全93个体Mantel r=0.801且显著，排除YL后r=?0.032不显著，证明整体隔离—距离信号主要由YL地理隔离贡献，而非其余样本普遍IBD（isolation by distance）。
讨论部分总结：8个PCR基础InDel标记成本低、适于初步种质筛查，但标记数与基因组覆盖有限，细尺度关系不及SNP。中低He与I可能源于断续分布、有效种群小、生境破碎或连接度受限，但仅8位点下属合理推测而非确证。AS多样、YL低多样却最分化，YL隔离获SNP与InDel两尺度支持。AMOVA、PCoA、STRUCTURE与UPGMA一致突出YL，但33个体SNP每地点仅3个、93个体仅8位点，群体划定仍需更密基因组与更宽地点内采样；五组是探索框架，不是正式保护单元或种源单元。保护上应兼顾高多样AS与高分化YL，并保留产地记录。
结论部分翻译：对33个C. acanthophylla个体的全基因组重测序鉴定出2,873,410个高质量SNP和5,679,915个InDel，为遗传评估与标记开发提供基因组资源。从34对候选引物中保留8个扩增稳定、分型清晰的多态InDel位点，并用于11个原始地点93个个体。全基因组SNP分析揭示明显区域分化与更细地点关联；同一33个体中8个InDel恢复了部分广泛区域模式，尤其YL分化，但细尺度分辨率较低。全部93个体中InDel面板显示中低标记水平多样性和可检测区域分化；AS多样性较高，YL多样性最低但区分度最强。基于11地点和五组的AMOVA模式相似，多数变异在采样单元间。遗传与地理距离的整体关联主要由YL地理分离驱动，排除YL后不再显著。这些发现支持代表性收集差异化YL材料，并仔细记录AS及其他区域种源。所开发InDel面板是初步种质表征、种源筛查和扩大个体采样的实用低成本工具；但五個探索分析组不应视为正式划定的生物种群、保护单元或种源单元。未来研究应把全基因组标记与更宽地点内采样及独立技术重复结合。
