形成中的地方种群：斐济桃花心木（Swietenia macrophylla）人工林的基因组多样性及迁地保护案例

《Trees, Forests and People》：A land race in formation: genome-wide diversity in Fiji’s mahogany ( Swietenia macrophylla) plantations and the case for ex situ conservation

　　

自20世纪60年代以来，大叶桃花心木（Swietenia macrophylla）已在斐济的维提岛（Viti Levu）和瓦努阿岛（Vanua Levu）广泛种植，利用从中美洲和南美洲多个国家进口的种子建立了约40,000公顷的人工林。该林产是一个正在形成的地方种群（land race），也是《濒危野生动植物种国际贸易公约》（Convention on International Trade in Endangered Species, CITES）附录II物种的意外迁地保存地，但其遗传组成尚未被表征。研究人员量化了斐济人工林种群内和种群间的遗传多样性，描述了种群遗传结构，并检验了空间分化，先验预期自种植以来的一至两代时间太短，不足以产生距离隔离（isolation by distance）。研究人员使用DArTseq^TM技术对来自维提岛三个地点、七个按队列分层（cohort-stratified）的种群的83个个体进行了基因分型。过滤后（检出率（call rate）≥0.80，可重复性（reproducibility）≥0.96，最小等位基因频率（minor allele frequency, MAF）≥0.05），保留了25,850个单核苷酸多态性（single nucleotide polymorphism, SNP）。所有个体均携带独特的多位点基因型（multilocus genotype, MLG）。观测杂合度（observed heterozygosity, H_O）为0.168–0.202，近交系数（inbreeding coefficient, F_IS）为0.092–0.173。全局Weir–Cockerham F_ST为0.0023，分子方差分析（analysis of molecular variance, AMOVA）将99.2%的变异归因于种群内个体间的差异。Mantel检验（r=0.008, p=0.378）和基于距离的Moran特征向量图（distance-based Moran's eigenvector maps, dbMEM）分析（p=0.159）均未检测到距离隔离，与预期一致。多样性在各地点和更新队列（recruitment cohort）中得以保留，分化可忽略不计——这是广泛、多种源（multi-provenance）奠基样本所预测的模式，而非当代基因流所致。苏瓦（Suva）拥有最高的等位基因丰富度（allelic richness）和1,233个私有等位基因（private allele），是迁地保护的重点。采伐规则应权衡亲缘关系而非仅凭表型，以防止商业采伐下的遗传侵蚀。维持这一基因资源的更新也标志着一个外来物种在高度特有性的群岛上无辅助繁殖。

大叶桃花心木（Swietenia macrophylla）是全球最具商业价值的硬木之一，原产于南美和中美洲，因过度采伐被列入《濒危野生动植物种国际贸易公约》（CITES）附录II。斐济自20世纪60年代起从多个中南美国家进口种子，在维提岛和瓦努阿岛建立约40,000公顷人工林，该林产被视为一个“形成中的地方种群”（land race），同时成为CITES附录II物种的意外迁地保存库。然而其遗传组成从未被表征。此前基于传统分子标记的研究受标记密度和基因组覆盖度限制，无法提供全基因组水平的评估。因此，研究人员利用DArTseq^TM高密度SNP技术，对斐济维提岛三个地点七个队列分层种群进行首次基因组学表征，旨在评估遗传多样性、种群结构和空间分化，为人工林保护与管理提供依据。该论文发表在《Trees, Forests and People》。

主要技术方法：研究人员在维提岛Nukurua（NKR）、Nausori Highlands（NHL）和Suva（SUV）三个地点采集叶片样本，按亲本（-p）、幼树（-sa）和幼苗（-se）队列分层，共83个个体（NKR-p=10, NKR-sa=11, NKR-se=9, NHL-p=9, NHL-sa=9, NHL-se=11, SUV-p=24）。采用DArTseq^TM技术进行基因分型，过滤后保留25,850个SNP。分析方法包括遗传多样性统计、Weir–Cockerham F_ST、AMOVA、Mantel检验、dbMEM分析、PCA、sNMF、DAPC和系统发育分析。

研究结果：
1. SNP标记质量：过滤后获得25,850个高质量SNP，保留率35.4%，平均检出率0.965，可重复性0.994，平均MAF 0.233，平均PIC 0.260，Ts/Tv比1.57。基因型累积曲线在20个位点饱和，说明标记辨别力充足。
2. 遗传多样性：83个个体均为独特多位点基因型（MLG）。所有种群的观测杂合度（H_O）为0.168–0.202，无偏期望杂合度（uH_E）为0.224–0.264，近交系数（F_IS）为0.092–0.173，显示杂合子缺失。等位基因丰富度（AR，稀化至n=9）介于1.20–1.25，SUV-p最高（1.25），且拥有1,233个私有等位基因和多态位点比例97.0%，其他种群私有等位基因数为0。总体连锁不平衡显著（r?_d=0.0066, p=0.001），但NKR-se和NHL-sa不显著。
3. 种群分化：全局Weir–Cockerham F_ST=0.0023，配对F_ST绝大多数接近零，最高为NHL-sa与NKR-p（0.013）。AMOVA显示99.21%的遗传变异存在于种群内个体间，地点间仅1.68%（p=0.001），地点内队列间为负值（-0.89%，p=0.888）。Mantel检验未检测到距离隔离（r=0.008, p=0.378），dbMEM分析亦不显著（p=0.159）。这与种植后仅一至二代、不足以形成距离隔离的预期一致。
4. 种群结构与系统发育：PCA前两轴解释方差仅9.6%和5.3%，多数个体聚于中心区，SUV-p和NHL-sa偏离。sNMF最低交叉熵在K=4，但聚类与地点/队列不对应。DAPC个体重分配至先验种群的比例低。NJ树显示大部分个体构成松散中央群，仅少数节点支持率>50%；NKR与南美参考基因组最接近。

讨论总结：研究人员认为，斐济人工林保留的遗传多样性与天然种群相当（H_O 0.14–0.28范围内），且没有出现等位基因丰富度相对杂合度下降的窄奠基特征，这与多国家多种源的奠基历史一致。正F_IS可能源于Wahlund效应、双亲近交和管理历史中的选择性采伐。SUV-p独特的高多样性和私有等位基因使其成为迁地保护优先对象，但需考虑样本量影响。极低分化归因于共享奠基种源、短时间和异交生物学。管理上，广泛跨地点采种可捕获大部分多样性；采伐规则应基于亲缘关系而非仅表型。人工林的自我更新既使其成为可持续的基因资源，也使其作为外来种在特有性高的群岛上无辅助繁殖。

结论翻译：利用25,850个基因组范围DArTseq SNP，本研究首次对斐济及更广泛太平洋地区人工种植的S. macrophylla进行了遗传多样性和种群结构表征。所有83个个体均携带独特多位点基因型，多样性在三个地点和更新队列间得以保留，种群水平分化可忽略不计（全局F_ST=0.0023），且未检测到距离隔离。这两个结果均符合对种植后仅一至二代人工林的预先预期。从地方种群的角度，这预测了广泛、多种源奠基样本的结局：共享奠基种源和不足以分化的时间解释了这一结果，无需归因于当代基因流或漂移-迁移平衡。苏瓦种群在共同样本量下具有最高等位基因丰富度（AR=1.25），并拥有1,233个私有等位基因和最高比例的多态位点，应在