《Molecular Ecology》：On the Causes of Correlated Genomic Ancestry Across Contrasting Hybridization Histories in a Monkeyflower Species Pair
编辑推荐：
杂交是塑造物种进化轨迹的强大力量，但其结果在类群间以及物种对内均高度可变。在本研究中，研究人员探究了塑造杂交祖先程度变异的进程，这种变异既存在于群体间，也存在于基因组范围内。研究人员利用低覆盖度测序数据，对来自Mimulus guttatus和Mimulus
杂交是塑造物种进化轨迹的强大力量，但其结果在类群间以及物种对内均高度可变。在本研究中，研究人员探究了塑造杂交祖先程度变异的进程，这种变异既存在于群体间，也存在于基因组范围内。研究人员利用低覆盖度测序数据，对来自Mimulus guttatus和Mimulus nasutus广泛重叠分布区内两个地理区域多个群体的782个个体进行了全基因组局部祖先推断。研究人员发现杂交祖先程度在不同群体间存在差异，支持了不同的历史杂交模式。然而，杂交祖先的基因组模式在各群体间具有相关性，表明它们是由平行进程所塑造的。地理上邻近的群体间相关性最高，包括同域和异域分布地点之间，这为渐渗并非受局部限制，而是通过迁移在整个景观中传播提供了证据。研究人员发现基因组特征可以预测杂交祖先及其在群体间的相关性。然而，与该领域其他物种的研究发现相反，这些模式不太可能由针对杂交祖先的简单连锁选择所驱动。高杂交祖先的基因组离群值通常在群体间共享，提示平行正选择对祖先起到作用。然而，已知与生殖隔离相关的位点对群体间祖先变异的预测能力较弱，表明作用于自然杂交群体的选择是高度多基因的，且其潜在遗传结构在空间中存在差异。总体而言，本研究展示了生态、人口统计和基因组特征如何相互作用以塑造杂交结果。
**论文解读：猴面花物种对中相关基因组祖先的成因**
**研究背景与科学问题**
杂交是生物进化中的重要驱动力，其结果是高度可变的，既受人口统计和随机过程影响，也受自然选择作用。基因组中的杂交祖先（渐渗）模式可揭示塑造杂交的历史和选择性过程。当前领域仍存在几个关键问题：不同杂交群体间的祖先模式相似性如何？共享模式是由基因组景观特征还是特定位点效应驱动？邻近群体间迁移对流布渐渗的影响有多大？以及已知的生殖隔离候选位点在自然群体中是否表现出可预测的平行祖先模式？Mimulus guttatus和M. nasutus系统是研究这些问题的理想模型，因其在多个独立地理区域均有杂交记录，且已有多个生殖隔离位点被定位。然而，缺乏跨地理区域的系统比较研究，限制了对杂交进程普遍性和特异性的理解。
**研究内容与结论**
研究人员利用两个独立杂交区域的数据：美国华盛顿州的Catherine Creek（CAC）和Little Maui（LM）群体，以及加利福尼亚州北部Sierra Nevada山麓和山地的新采样群体。通过低覆盖度全基因组测序，对782个个体进行了局部祖先推断。研究发现：（A）渐渗量、对称性和基因组分布在地理和基因组空间中广泛变化；（B）群体间共享祖先模式由基因组景观特征和位点特异性模式共同驱动；（C）群体间基因流将渐渗等位基因传播到异域地区；（D）定位研究中的候选位点仅能部分且不一致地预测自然渐渗模式。研究发表后，为理解杂交后选择过程提供了重要视角，表明生殖隔离的遗传基础具有高度多基因性和空间变异性。
**主要技术方法**
研究使用低覆盖度全基因组测序数据（平均约125万paired reads/个体），样本来自加利福尼亚州北部（382个个体，其中101个来自主要地点Don Pedro Reservoir）和华盛顿州（503个已发表个体）。采用ancestry_HMM方法进行局部祖先推断，基于33个自交系的参考panel（包括M. nasutus、北方和南方进化枝的M. guttatus以及CAC渐渗个体）筛选了201,387个祖先信息标记。使用PCAngsd进行群体结构分析，并通过线性模型和方差分析评估基因组特征（染色体位置、基因密度、重组率）和缺失数据对祖先频率的影响。采用置换检验评估离群值重叠的显著性，并使用卡方检验评估QTL区间内离群值的富集情况。
**研究结果**
*3.1 多个接触区域导致多种杂交结果*
通过基因组PCA分析，PC1分离南北M. guttatus，PC2分离物种。所有五个地理分组均存在可检测的杂交祖先变异，甚至异域群体也有少量渐渗。北方群体呈不对称模式，hybrid群体偏向M. guttatus祖先；南方山麓同域群体中hybrid个体较少，但guttatus个体的HI值常大于0.01；南方山地同域群体则呈现较对称的分布，包括一些携带少量M. guttatus祖先的M. nasutus个体。
*3.2 基因组间祖先频率高度变异且在不同地区群体间强相关*
所有群体中祖先频率在基因组间变异显著。北方hybrid群体的祖先频率分布呈不对称性，而南方山地hybrid群体中大量窗口处于50%频率，提示早期世代杂交的存在。任意两个guttatus或hybrid群体间的祖先频率可解释对方大量基因组变异，表明存在显著的基因组平行性，且该效应甚至延伸到南北地区间的比较。
*3.3 基因组结构特征和技术因素解释共享方差，而非连锁选择*
染色体相对位置平方、基因密度和重组率均与祖先频率显著相关，但基因密度和染色体位置是比重组率更好的预测因子。值得注意的是，基因密度/重组率比值（连锁选择强度的更直接代理）与祖先频率仅弱相关（最大R2=0.007），表明这些关联并非直接由连锁选择驱动。缺失数据比例也解释了0.4%至30.6%的祖先频率变异，提示技术变异可能影响祖先推断。
*3.4 地理邻近性是杂交祖先平行模式的主要来源*
校正基因组特征和缺失数据后，群体间残余相关显著降低但仍具显著性。群体内不同cohort间（8%-35%）和邻近群体间（5%-15%）保留最强的残余相关。南方异域M. guttatus群体的祖先频率最好由南方同域M. guttatus群体解释（5%变异），且杂交祖先在更靠近最近同域地点的异域群体中显著更高（adj.R2=0.73），支持基因流从同域向异域传播渐渗等位基因。
*3.5 共享变异模式提示正选择塑造渐渗*
在考虑基因组结构特征和缺失数据后，南方同域M. guttatus中2.1%-3.6%的变异仍可由北方群体祖先频率解释。顶部5%高minor-parent祖先的离群窗口在几乎所有比较中均比随机预期更频繁地重叠，即使在南北地区间也存在显著重叠。校正基因组特征后重叠程度降低但未消除，提示平行正选择在塑造共享祖先模式中起作用。
*3.6 生殖隔离QTL和候选基因的祖先频率偏差不一致*
对于预测M. nasutus等位基因有害的N? QTL，仅北方CAC hybrids的祖先频率显著低于预测值。对于预测M. guttatus等位基因有害的N+ QTL，所有南方M. guttatus和hybrid群体以及北方CAC M. guttatus的祖先频率显著高于预测值。四个已知生殖隔离候选位点（两个光周期位点、两个幼苗致死互作位点）周围在部分群体中观察到M. nasutus祖先频率降低，但这些位点远非其所在染色体上最极端的祖先偏差，且模式在群体间常不一致。
**讨论与结论**
研究揭示了杂交结果的广泛空间变异，支持了基因流存在下的物种形成模型。地理邻近群体的祖先相关性提示基因流在景观中传播渐渗等位基因。与许多其他系统不同，未发现连锁选择塑造基因组祖先模式的证据，这可能与植物物种中较弱的连锁选择有关。已知生殖隔离位点仅能部分且不一致地预测自然渐渗模式，表明这些位点的多态性、群体特异性选择模式以及高度多基因的生殖隔离架构使实验室定位结果难以直接应用于自然群体。总体而言，该研究强调杂交基因组学领域需要结合跨系统的比较工作、稳健理论和建模，以及更深入理解作用于杂交个体的多维选择景观，以更好地理解杂交祖先基因组变异的根本原因。
