为了探明秦岭石蝴蝶花瓣数量变异原因,该研究采用Illumina HiSeq 2500高通量测序技术对秦岭石蝴蝶两种花型发育的早期和晚期进行转录组测序,挖掘参与其花发育相关的差异基因,并探讨花器变异的可能机制.结果显示:(1)与NR数据库进行比对,共有52677个Unigene注释到NR库,占Unigene总数的46.25％,与旋蒴苣苔(Dorcoceras hygrometricum)的序列同源性最高(54.29％).(2)GO富集分析结果显示,正常2-3型和变异2-4型差异基因富集最显著的GO条目为:细胞组分类的细胞膜、分子功能类的反转运蛋白活性和酶抑制剂活性、生物过程类的跨膜转运和催化活性的负调控等;KEGG富集分析结果显示,正常2-3型和变异2-4型差异基因富集最显著的KEGG通路为:植物激素信号转导、脂肪酸延长、戊糖和葡萄糖醛酸酯的相互转化、苯丙烷生物合成、硫代谢、类黄酮生物合成、玉米素的生物合成等途径.(3)对表达差异基因进行筛选并进一步对4个比较组进行交叉比对分析,确定了6个可能与秦岭石蝴蝶花器官发育相关的基因,分别为PqMIF2、PqMYB340、PqMYB305、PqGA-TA12、PqCCD4和PqZBED;qRT-PCR验证发现,其表达趋势和转录组测序分析结果一致.该研究结果为秦岭石蝴蝶的花器官发育和系统进化及其濒危机制方面研究提供了重要的信息.
秦岭石蝴蝶 花器变异 转录组 花发育基因 KEGG分析