《Biology》：Population Genetic Structure and Diversity of Camellia oleifera: Insights from Comparative Chloroplast Genomics and Development of Molecular Markers
编辑推荐：
本研究旨在通过比较叶绿体基因组（plastome）分析开发叶绿体序列标记，并评估其在所采样油茶（Camellia oleifera）群体中的单倍型鉴定和遗传多样性评估中的应用。新组装的3个油茶组（Camellia sect. Oleifera）叶绿体基因组长度
本研究旨在通过比较叶绿体基因组（plastome）分析开发叶绿体序列标记，并评估其在所采样油茶（Camellia oleifera）群体中的单倍型鉴定和遗传多样性评估中的应用。新组装的3个油茶组（Camellia sect. Oleifera）叶绿体基因组长度范围为156,918–157,003 bp，并保留了典型的四部分结构。对这些基因组与13个参考叶绿体基因组的比较分析鉴定出605个长重复序列和3315个简单序列重复（SSR）。三个新叶绿体基因组共享29个偏好有义密码子，其中28个以A/U结尾。滑动窗口核苷酸多样性（π）分析突出了8个候选变异区，最大π值为0.01397。在来自37个采样群体的324个个体中，1772 bp的拼接比对包含38个变异位点和55个单倍型（Hd = 0.8052；π = 0.002183）。其中36个单倍型为单例，而H003出现在35个群体的133个个体中。序列主成分分析（PCA）和探索性ADMIXTURE分析描述了广泛的变异。在K = 2–10的展示值中，K = 2的交叉验证误差最低，其两个组分被描述性解释。这些结果记录了所采样C. oleifera材料中的叶绿体序列多样性和群体结构，并为种质资源鉴定和保护采样提供了分子资源。
油茶（Camellia oleifera）是山茶科山茶属经济木本油料作物，种子油富含营养，在中国南方地区广泛栽培，已有两千多年利用历史。然而，随着油茶产业规模扩大和人工选育品种大面积推广，自然油茶资源面临严峻挑战：栖息地破碎化、过度采挖和人为干扰导致野生群体规模明显下降，遗传资源流失风险加剧；同时，现有育种对野生种群遗传背景利用不足，造成栽培品种遗传基础狭窄，制约了抗逆性和品质提升。因此，系统分析油茶群体遗传结构、遗传多样性与进化动态，对遗传资源保护和可持续育种具有重要意义。叶绿体基因组呈母系遗传，结构保守但单拷贝区变异丰富，是理想的群体遗传标记。以往使用的ISSR、RAPD等标记存在同源性判断不明、重复性差、覆盖不足等局限，而基于叶绿体基因组高变区开发的序列标记能够更可靠地表征叶绿体序列变异。本研究旨在通过比较叶绿体基因组学开发叶绿体序列标记，并评估其在油茶群体中的单倍型鉴定和遗传多样性评估中的应用。
研究人员新组装了3个油茶组（Camellia sect. Oleifera）叶绿体基因组（YC127、XYYC和HNL2-2，分别采自浙江金华国际山茶物种园和四川雅安天全县），并与13个参考叶绿体基因组进行比较分析；通过滑动窗口分析核苷酸多样性（π）识别候选高变区，设计出三个叶绿体标记（trnG-UCC、ndhF、trnL-UAA–trnF-GAA）；对来自7个省份37个采样群体的324个油茶个体进行PCR扩增和测序；拼接后计算单倍型多样性（Hd）和核苷酸多样性（π），并采用最小生成网络、序列主成分分析（PCA）、邻接法（NJ）和ADMIXTURE软件分析群体结构。
3.1 叶绿体基因组特征：三个新叶绿体基因组长度为156,918–157,003 bp，均具有典型的四部分结构，共注释132个基因位点，包括87个编码序列（CDS）、37个tRNA和8个rRNA。GC含量在反向重复区（IRs）高于大单拷贝区（LSC）和小单拷贝区（SSC），符合rRNA和tRNA富集的特征。
3.2 重复序列与IR/SC边界比较：16个叶绿体基因组中共鉴定出605个长重复序列和3315个简单序列重复（SSR），其中单核苷酸重复占74.18%，以A/T富集为主。边界分析显示ycf1基因在所有基因组中跨越SSC-IRa边界，rps19基因在15个基因组中跨越LSC-IRb边界，仅C. oleifera MF541730.2例外。
3.3 密码子使用模式：三个新基因组共享29个偏好有义密码子（相对同义密码子使用频率RSCU>1），其中28个以A/U结尾，表明它们具有相似的密码子使用偏好。
3.4 核苷酸多样性与候选可变区：滑动窗口π最大值为0.01397，位于trnS-GCU–trnG-UCC区域；共筛选出8个候选高变异区，包括trnS-GCU–trnG-UCC、psbK–psbI、petN–psbM、trnL-UAA–trnF-GAA、petB内含子、ycf15–trnL-CAA、ndhF和SSC区的ycf1拷贝。IR区的平均π低于LSC和SSC，显示IR区更保守。
3.5 遗传多样性：三个标记分别检测到9、31和8个片段单倍型，平均多态信息含量（PIC）为0.4789，其中ndhF标记（标记2）多样性最高。在324个个体的1772 bp拼接比对中，共发现38个变异位点和55个单倍型，总体Hd=0.8052，π=0.002183。36个单倍型为单例，最频单倍型H003出现于35个群体的133个个体中。
3.6 cpDNA单倍型网络：最小生成网络显示55个单倍型，常见单倍型在群体间共享，同时存在大量以单例形式存在的稀有单倍型，表明单倍型丰富度与序列差异程度并不完全对应。
3.7 基于叶绿体DNA的遗传结构：序列PCA中PC1解释64.96%方差，PC2解释8.57%方差，个体在排序空间中重叠明显；ADMIXTURE在K=2时交叉验证误差最低，两个组分被描述性解释；NJ树内部支持率较低（50.1%–73.0%），因此多数分支关系未被充分支持。
讨论指出，16个油茶组叶绿体基因组在基因顺序和整体结构上高度保守，但IR边界存在细微差异，这种保守框架下的局部变异为标记开发提供了靶点。三个标记能够有效检测油茶种内序列变异，其中编码基因ndhF反而表现出最高的片段单倍型丰富度，说明编码区也可能含有有用的多态性。本研究中广泛共享的单倍型与大量单例单倍型共存，提示单倍型丰富度并不等同于大尺度序列差异，遗传多样性指标反映的是不同方面的特征。群体结构分析显示采样群体间存在重叠，K=2模型组分仅作为探索性描述，不能等同于两个独立的分化群体，且有限的NJ支持率也限制了对具体分支关系的解释。
研究结论可翻译为：比较叶绿体基因组分析在整体保守的基因组框架中鉴定出局部序列变异，为叶绿体标记开发提供了靶点。在来自37个采样群体的324个个体中，三个标记在1772个完整位点中共产生55个单倍型（Hd=0.8052；π=0.002183）。标记2具有最高的片段单倍型丰富度和PIC。一个在35个群体中共享的常见单倍型与36个单例单倍型共存。序列PCA和描述性K=2模型概括了叶绿体变异的总体模式，而有限的NJ引导支持率限制了对具体分支关系的解释。这些标记为种质资源表征和多样性评估提供了基础，但其对更多材料的适用性仍需进一步评估。
