《The Plant Genome》：Development and deployment of a public mid-density genotyping panel for crimson clover (Trifolium incarnatum L.) based on whole-genome re-sequencing of diverse germplasm
绛三叶草（Trifolium incarnatum L.）是一种专性异交的冷季一年生豆科植物，因饲草和覆盖作物价值而受重视，但支持系统改良的基因组资源有限。研究人员开展了全球绛三叶草种质的首次全基因组重测序（whole-genome re-sequencing, WGR），以：（i）开发用于育种应用的公共中密度富集捕获面板；（ii）表征采样基因型间的遗传多样性，以验证种质对面板设计的代表性。研究人员从45份种质材料中对194个基因型进行了测序：45个以～50X测序，产生5.84 million个变异；149个以～2.54X测序，产生17.05 million个变异。过滤后，研究人员分别从高覆盖度数据集保留542,790个高置信度单核苷酸多态性（single-nucleotide polymorphisms, SNPs），从低深度队列保留约2.4 million个SNP。基于这些资源，研究人员设计了一个28,867-SNP的Twist杂交捕获面板，其富集基因区并均匀分布于7条染色体。该面板正部署于Auburn University的绛三叶草育种项目中，以支持群体改良和品种开发。遗传多样性分析（主成分分析 [principal component analysis, PCA]、ADMIXTURE）揭示栽培品种紧密聚类、野生和状态不确定种质材料分散更广，总体分化低（固定指数 [fixation index, FST] = 0.0105），并存在杂合度过剩（近交系数 [inbreeding coefficient, FIS] = ?0.0592），与高水平异交一致。由此产生的基因组资源为性状发现、预测育种和多样性监测提供了可重复的中密度基因分型平台。总之，这些进展使绛三叶草基因组资源与其他豆科植物如大豆（Glycine max L. Merr.）和紫花苜蓿（Medicago sativa L.）相当，为美国可持续农业中这一关键覆盖作物和饲草作物的基因组学辅助改良奠定了坚实基础。
论文解读：
绛三叶草（Trifolium incarnatum L.）是专性异交的冷季一年生豆科植物，作为饲草和覆盖作物具有高生物量、固氮和土壤培肥价值。美国2022—2023年全国覆盖作物调查显示，绛三叶草在豆科中排名第一、所有覆盖作物中排名第四，仅次于黑麦（Secale cereale L.）、萝卜（Raphanus sativus L.）和冬小麦（Triticum aestivum L.）。种子年产量约2.6–4.4 million kg，农场销售额$3.9–$6.2 million，生产集中在俄勒冈Willamette Valley；2018—2020年收获面积2315–3650 ha，平均产量约1120 kg ha
?1。然而，与大豆、紫花苜蓿和红三叶草（Trifolium pratense L.）相比，绛三叶草基因组资源匮乏；在Trifolium属中，多数研究集中于红三叶草和白三叶草（Trifolium repens L.），而绛三叶草仅知C值基因组大小约652 Mb、2n=14，以及早期扩增片段长度多态性（amplified fragment length polymorphism, AFLP）多样性调查。简单序列重复（simple sequence repeats, SSR）和测序基因分型（genotyping-by-sequencing, GBS）等简化代表方法仅捕获小部分基因组，常有缺失数据、位点丢失和批次效应。WGR虽成本高，但可提供发现层，用于设计高效简化代表基因分型工具；杂交捕获富集面板可选择性测序数千个预定区域，提供一致标记集、降低缺失数据并改善跨研究可比性。因此，研究人员开展首次绛三叶草全球种质WGR，目标为生成综合SNP资源、开发公共中密度全基因组杂交捕获面板（约28K SNP）、表征全球面板遗传多样性以验证种质代表性。结果形成首个全基因组变异资源，将变异发现与应用育种桥接。
研究人员对45份种质、194个基因型进行两级WGR，45个～50X和149个～2.54X；高覆盖识别5.84 million变异，低覆盖17.05 million变异；过滤后分别保留542,790和约2.4 million高置信SNP；设计28,867-SNP Twist面板；遗传多样性显示栽培品种紧密聚类、野生和不确定材料分散、总体F
ST=0.0105、F
IS=?0.0592；面板已部署于Auburn育种项目。其意义是提供可重复中密度平台，用于性状发现、预测育种和多样性监测，使绛三叶草与大豆和紫花苜蓿等豆科资源相当。
关键技术方法（不超过250字）：研究人员从45份绛三叶草种质中获得194个基因型，其中40份来自USDA NPGS，2份来自Auburn University（AU），2份来自University of Florida（UF），1份商业品种Dixie作为重复对照；田间存活194株，包括154个PI、18个AU、5个UF和17个Dixie。采用Illumina NovaSeq 6000/X进行150 bp双端测序，形成45个～50X和149个～2X的WGR队列。读段经BBTools过滤后以BWA-MEM比对至红三叶草参考基因组GCF_020283565.1 ARS_RC_1.1；SAMtools、bcftools和PLINK用于变异检测与过滤，Beagle v5.4用于定相和插补。面板设计整合JGI 50X、JGI 2X及6803个Khufu SNP，经质量过滤、链歧义去除、LD修剪、基因区并集、±60 bp扩展和染色体居中分箱选择，由Twist Bioscience合成。验证使用2024—2025年1039份育种样本经Khufu低通测序和参考引导插补；另以50X创始人面板引导2X插补。
3.1 WGR与变异及SNP发现：研究人员通过45个～50X WGR基因型获得平均映射效率78.16%，覆盖39.48X–101.06X，均值65.44X；比对至红三叶草参考基因组423–437 Mb；识别约5.8 million总变异，其中约4.96 million SNP；去除>30%缺失、MAF≥0.05并插补后保留542,790高置信SNP；染色体2最多93,614，染色体5最少55,178；平均密度约1347 SNP/Mb，部分>2500 SNP/Mb，染色体1、2、3最高。2X队列149个基因型平均映射效率82.19%，覆盖0.67X–6.98X，均值2.54X；总变异17,050,990，其中SNP 16,051,990；过滤后3,913,665，Beagle插补后2,411,468高置信SNP。结论是两级数据提供全面基因组基础。
3.2 开发SNP富集面板：研究人员整合JGI 50X、JGI 2X和6803个Khufu SNP，共17.8 million raw SNP，经多步过滤：QC后约365,000，去除链歧义后约300,000，LD修剪后约119,000，基因区并集后±60 bp扩展产生约66,000非冗余窗口，约84%与基因区相交；bin-center选择24,000，与Khufu合并为29,191，top-up 809个基因区间，最终约30,000捕获探针；最终面板含28,867 unique probes，染色体分布Chr1=3168、Chr2=4400、Chr3=3213、Chr4=3574、Chr5=3421、Chr6=2725、Chr7=3499；多数染色体基因区>75%，最高79.1%，最低66.0%；平均间距约16.9 kb；Twist验证>99.9%探针回收、>98%映射效率、归一化信息熵>0.92、平均覆盖>30X。结论是设计均匀、特异、可制造。
3.3 基因分型完整性与等位基因频率特征：研究人员将面板应用于1039份育种样本，经Khufu变异检测和插补获得72,181 SNP；所有样本和几乎所有位点call rate>95%，多数>99%；69.3%位点MAF≥0.05，即50,021个信息性SNP；MAF均值0.239、中位0.236，富集于0.3–0.5中间频率。结论是面板在育种规模上缺失率低、标记信息量高。
3.4 参考基因组映射率与富集面板覆盖：研究人员分析45个50X基因型映射率均值79.51±3.73%，范围62.49%–93.16%；在410个非重叠1-Mb窗口中，平均映射深度与探针数弱正相关（Spearman ρ=0.333，p<0.001；置换p<0.0001），ρ
2=0.111；每染色体平均探针密度45.2–56.3 probes/Mb；无1-Mb窗口无探针，最小10、中位53、IQR 46–57；去除18个离群窗口后ρ=0.374。结论是映射质量影响较小，bin-center算法和基因区标准保障均匀覆盖。
3.5 地理分布与遗传多样性：45份种质来自北美、欧洲、亚洲和非洲，其中28份栽培包括Dixie、AU Sunrise和UF Early White，19份欧洲野生；基于Euclidean距离的UPGMA显示中等地理分组，但区域间无清晰分离，反映异交和历史种质交换；PI655002（Hatif à Fleur Rouge）、PI418900（279）和PI591666（93-64）形成独特分支。结论是欧洲野生种质存在局部遗传分化。
3.6 采样基因型间的遗传多样性模式：PCA在50X中PC1解释15.72%、PC4解释9.4%，按地理和改良状态聚类；2X中PC1 20.27%、PC2 12.93%，美国与欧洲基因型分离。ADMIXTURE K=2交叉验证误差最低（0.52792±0.00003），K=3为0.53386±0.00005；K=3第三成分低频，主要见于少数亚洲和非洲基因型。F
ST显示美国与欧洲0.0066、美国与非洲0.0117、欧洲与非洲0.0172；亚洲与美国0.062、与欧洲0.0706、与非洲0.1537；改良类别中野生与育种0.0019、栽培与育种0.0088、野生与栽培0.0027、野生与不确定0.0051；总体F
ST=0.0105，F
IS=?0.0592。50X Venn分析显示美国独有SNP 2778、亚洲1061、非洲1012、欧洲197，四区域共享37,892；改良类别共享核心460,066，独有SNP为育种材料3215、栽培675、不确定58、野生92；栽培各染色体独有SNP最多，其次野生，育种较少。结论是弱群体分化、高异交、存在驯化相关遗传整合。
讨论部分总结：4.1 研究人员认为本研究是首个全基因组评估，50X产生约5.8 million变异和约4.96 million SNP，过滤后542,790高置信SNP，优于以往GBS和SSR；2X增加149个基因型，产生约2.4 million高置信SNP，平均PIC=0.202；2X变异数更多主要因联合变异检测。4.2 资源可转移至开放授粉专性异交物种育种；2022年启动轮回选择，包括SPT、PYT和ALT，未来在SPT阶段利用PYT/AYT产量数据进行基因组选择。28,867 Twist面板含28,931个探针、24,106个合并区域、21,051个成功平铺，覆盖84.3%（3.14 Mb），>75%靶向基因区，间距约16–17 kb；已用于Auburn轮回选择家系；与大豆SoySNP50K、苜蓿、花生Axiom_Arachis2等平台逻辑一致。Twist确认合成和杂交质量；但T. incarnatum特有结构或片段重复可能造成参考评估未检出的脱靶，未来de novo组装可评估。还指出不能直接混合50X和2X进行多样性分析，覆盖偏差会扭曲AF估计，PCA可见分离，因此多样性、ADMIXTURE和F
ST应仅用2X。4.3 中密度面板验证中，1039份育种材料经Khufu低通和插补后call rate高；MAF 0.3–0.5富集，同时保留低频等位基因；与固定SNP阵列相比，杂交捕获模块可更新，减少确定偏倚；牛Twist捕获与SNP阵列一致约98%、R
2=0.98；较GBS位点稳定、减少丢失；50,021个信息性常见标记支撑育种。4.4 两级框架以50X founders（C0）为参考，未来C1–C3低通～2X；跨队列协调后488,078个候选位点，Beagle插补平均DR
2=0.95，下游用DR
2≥0.8和MAF≥0.01；可扩展至<1X（约0.8X）并保持常见变异MAF≥0.05高准确度。4.5 遗传多样性与Steiner et al.早期AFLP结果一致，美国栽培窄基，欧洲和野生分散；WGR将解析度从161个AFLP扩展到194个基因型的2.95 million高置信SNP；US与欧洲共享背景，总体F
ST=0.0105、F
IS=?0.0592，支持广泛随机交配、弱地理分化和驯化相关遗传整合。4.6 限制是缺乏T. incarnatum参考，使用红三叶草ARS_RC_1.1；平均映射率79.51±3.73%，约20%读段未比对；面板覆盖在可映射区域均匀，但T. incarnatum特有区域未代表；已启动单倍型解析染色体级组装。
研究结论翻译：本研究提供了绛三叶草首个全基因组基因组资源，并为现代育种和基因库支撑的作物改良建立了实用基础。通过对全球代表性面板进行WGR，研究人员生成了高置信度的全基因组变异目录，其分辨率远超此前基于标记的研究。这些数据表征了采样基因型间的遗传多样性，捕获了稀有和私有等位基因，并提供了该物种现有遗传多样性的最完整视图。从育种角度看，这些资源使绛三叶草从历史上资源不足的物种转变为能够充分参与预测性和数据驱动改良管道的物种。公共可用的中密度（28,867个SNP）杂交捕获富集面板的开发，将全基因组变异发现与实际育种桥接起来。该面板已合成、验证并部署在Auburn绛三叶草育种项目中，使全基因组标记信息能够直接整合到轮回选择管道中。除支持基因组预测外，该平台还能追踪选择周期中的等位基因频率变化，并更有效管理遗传多样性，这是开放授粉、高度杂合物种持续遗传增益的关键要求。同样重要的是这些资源对基因库和种质管理的价值。全基因组SNP数据集能够稳健地表征种质材料、鉴定遗传上独特或代表性不足的材料，并检测收藏中的冗余。标准化和可转移标记的可用性进一步支持基因库与育种项目之间的互操作性，促进对保护多样性的更有效利用。这些进展使绛三叶草与大豆、苜蓿和花生中已建立的基因分型平台和SNP面板资源保持一致，使该物种与其他资源充足的豆科植物达到同等水平，并实现长期基因组学驱动的改良。
