《Diversity》：The Genetic Diversity of Old and Modern Accessions of Glycine max Revealed by Genotyping-by-Sequencing (GBS)
编辑推荐：
大豆（*Glycine max* (L.) Merr.）是全球广泛种植的重要豆类作物之一。对全球大豆样本进行遗传分析对于有效的育种计划以及基因组和标记辅助选择的发展至关重要。本研究调查了来自VIR种质库（世界上最大的种子库之一）的大豆种质资源的遗传多样性。研究
大豆（*Glycine max* (L.) Merr.）是全球广泛种植的重要豆类作物之一。对全球大豆样本进行遗传分析对于有效的育种计划以及基因组和标记辅助选择的发展至关重要。本研究调查了来自VIR种质库（世界上最大的种子库之一）的大豆种质资源的遗传多样性。研究样本包括163个大豆（*G. max*）品种。研究采用基因分型测序（Genotyping-by-Sequencing，GBS）技术进行。研究目的是评估栽培大豆遗传多样性的时间变化、群体结构和连锁不平衡。研究人员根据两个标准对材料进行划分：（1）三个时间组群，即第I组（1959年之前）、第II组（1960–1989年）和第III组（1990年之后）；（2）地域特征。连锁不平衡（Linkage disequilibrium，LD）衰减分析揭示了时间组群之间的差异，其中1959年前开发的早期种质具有最长的LD衰减距离，而较晚组群的LD衰减距离较短。Nei基因多样性在各时间组群间仅显示出微小差异，较新的种质呈现中度下降趋势。关键位点被定义为在时间组群间等位基因频率差异最大且位于注释基因内部或附近的SNP。研究结果凸显了全面基因组分析对于遗传资源的保护和有效利用以及未来大豆育种策略制定具有重要意义。
**基于GBS技术的栽培大豆遗传多样性研究解读**
**1. 研究背景与问题**
大豆（*Glycine max* (L.) Merr.）起源于中国，经过数千年的驯化和传播，已成为全球广泛种植的重要经济作物。野生大豆（*Glycine soja*）约在27万年前起源于东亚，而栽培大豆的驯化发生在大约5500至6000年前的中国。随着大豆向全球传播，不同地区反复引入和人工选择导致遗传多样性的变化。传统育种方法主要依赖形态和农艺性状，难以全面评估遗传多样性，可能导致部分优良等位基因的丢失。自20世纪80年代起，分子标记（如RFLP、SSR、SNP）开始应用于大豆研究，随后高通量测序技术（NGS）使得分析数十万个SNP成为可能。然而，先前研究多集中于现代品种或国家种质库，较少利用历史种质资源系统评估不同育种时期大豆遗传多样性的时间动态。世界主要大豆种质库包括美国USDA、巴西EMBRAPA、中国国家大豆基因库、日本NIAS以及俄罗斯VIR。VIR大豆种质库是世界上最丰富的种质库之一，包含超过7800份来自72个国家的样本，其中大部分为栽培大豆。鉴于不同历史时期育种方法和优先目标的变化，利用VIR的历史种质资源研究大豆遗传多样性随时间的变化具有重要意义。本研究旨在通过GBS技术分析来自VIR种质库的163份大豆品种，评估不同育种时期和地理来源的遗传多样性、群体结构及连锁不平衡（LD）的动态变化，并鉴定与育种历史相关的基因组区域。
**2. 主要技术方法**
研究人员从VIR种质库中选取163份栽培大豆（*G. max*）品种，覆盖22个国家，并根据育种年份分为三个时间组群：第I组（1959年前，41份）、第II组（1960–1989年，68份）、第III组（1990年至今，54份）。采用改良的Adapterama III方案构建GBS文库，使用*MspI*和*PstI*限制性内切酶双酶切基因组DNA，连接双索引接头后，通过PCR扩增和纯化获得200–600 bp片段文库，并在Illumina NovaSeq 6000平台进行双端150 bp测序。原始数据经FastQC和MultiQC质控，用fastp去除接头和低质量序列，使用Stacks工具包的process_radtags进行样本拆分。高质量reads通过Bowtie2比对至大豆参考基因组（GCA_000004515.5，Williams82），平均测序深度为27.68x。使用Stacks的gstacks模块和bcftools进行变异检测，过滤后获得17,848个双等位SNP（MAF>1%，缺失率<50%）。群体结构通过ADMIXTURE软件（K=2–4）和主成分分析（PCA）推断。遗传多样性按Nei公式计算，并使用95%置信区间估计。LD衰减使用PopLDdecay计算，基于r2值在1 Kb窗口内平均。等位基因频率变化通过计算第I组和第III组之间的绝对差异，并筛选差异最大的SNP，结合参考基因组注释确定候选基因。
**3. 研究结果**
**3.1 群体结构**
ADMIXTURE交叉验证表明最可能的遗传簇数为K=4。在K=2时，来自日本、德国、波兰、瑞典及部分俄罗斯品种聚为一簇，反映其遗传相似性。K=3时，中国、俄罗斯、摩尔多瓦和乌克兰的品种形成一组，而加拿大和美国品种结构明显不同。K=4时，额外区分出包含部分中国、摩尔多瓦和乌克兰品种的特定组群。PCA分析未显示按国家严格分离，大部分地理群体重叠，但沿PC1和PC2存在部分分离，表明中等程度的遗传分化，与ADMIXTURE结果部分对应。
**3.2 遗传多样性分析**
基于17,848个SNP位点计算的Nei基因多样性平均值为0.202，三个时间组群间差异很小（第I组0.207，第II组0.196，第III组0.190），呈现轻微下降趋势。等位基因频率分析显示，第I组与第III组之间，25%的标记频率变化超过0.1，约7%超过0.2，其中18个标记变化达50%以上。变化最大的SNP位于多个染色体区域，其中NC_038253.2:33499993位点频率从0.77降至0.13（差异0.64），位于基因间区。其余17个SNP位于或邻近注释基因，如*Glyma.17G205700*、*Glyma.13G045400*、*LOC100777354*等。15个SNP在现代化品种中频率下降，而3个位于4号染色体（*Glyma.04G181200*、*Glyma.04G183600*、*Glyma.04G182900*）的SNP频率上升。区域遗传多样性分析显示，加拿大品种最低（0.172），中国品种最高（0.203），俄罗斯、美国、摩尔多瓦介于其间。
**3.3 连锁不平衡**
LD衰减分析显示，第I组LD衰减距离最长（66,186 bp），第II组最短（48,170 bp），第III组为51,174 bp。早期品种LD距离较大表明重组较少，而现代品种LD距离缩短可能反映遗传基础拓宽。
**4. 讨论与结论**
**讨论**：研究利用VIR历史种质库，揭示了不同育种时期大豆遗传结构的变化。ADMIXTURE聚类模式与种质引入历史一致，例如K=2聚类可能反映日本种质对北欧和俄罗斯育种的影响，K=3和K=4聚类可能与中国大豆种质大规模引入苏联有关。15个SNP频率下降可能反映选择压力变化，而3个SNP频率上升可能与高产改良相关。候选基因的功能注释表明其可能参与离子稳态、胁迫耐受（*Glyma.09G052000*）、碳代谢和类黄酮生物合成（*Glyma.13G045400*）、膜转运（*Glyma.14G193300*）以及光调节和种子成分形成（*Glyma.17G205700*、*Glyma.20G099500*）。LD衰减趋势与历史育种实践相符，早期品种LD较长，中期缩短，现代品种略有回升但低于早期。遗传多样性水平与先前研究相当，但中国品种多样性最高，反映了驯化中心和丰富原始材料。研究局限性包括样本量不均衡、GBS固有的缺失数据以及候选区域需进一步验证。
**结论**：通过GBS对VIR种质库163份栽培大豆的分析，研究揭示了遗传多样性和基因库结构的动态特征。早期品种LD较高，表明遗传重组较少和育种材料来源较窄；现代品种遗传基础更广。遗传多样性整体稳定，但存在轻微下降趋势。中国品种多样性最高。等位基因频率变化的位点可能为驯化和育种改良的候选区域，需功能验证。遗传聚类与地理来源无严格对应，反映了复杂的种质引入和交换历史。研究结果强调了综合基因组分析对大豆遗传资源保护、监测和有效利用的价值。未来计划扩大样本量并整合表型数据，以深入理解大豆适应和改良的机制。
