课题组此前的工作揭示了喜马拉雅山东部作为枯叶蛱蝶属的分化中心,存在同域分布的多个枯叶蛱蝶属物种。利用各地采集样本的测序数据,本研究发现基因组上多个区域的系统发生关系与全基因组水平的系统发生关系存在不一致,提示渗入和不完全谱系分选的存在。
图1 枯叶蛱蝶属地理分布、系统发生关系和遗传混合格局
研究利用D和fd统计量来评估物种间的渗入水平,发现同域分布的指斑枯叶蛱蝶和隐种枯叶蛱蝶存在大量的基因流。
这些渗入信号在Z染色体上富集。对渗入窗口内的基因进行功能富集分析,找到一个代表性的候选基因pale。
在中华枯叶蝶中使用CRISPR/cas9技术对pale进行基因敲除,发现该基因影响了翅鳞片的色素沉着和鳞片结构发育。
通过转录组数据进行表达网络分析,发现pale基因存在雌雄表达差异,且与其它一些渗入基因存在表达上的关联。
这种表达差异反映了pale可能与翅图案的性二型有关。由于pale位于Z染色体上,而Z染色体存在渗入信号的富集。对其进行进一步分析,研究发现这些频繁的渗入已经改变了Z染色体的系统发生关系,通过模型计算推测这些渗入能在Z染色体上富集更可能是具有适应性意义而非单纯的遗传漂变。
图6 翅图案的反射光谱测量和Z染色体上的渗入信号富集