基于RNA-Seq技术分析在金属矿上部生长的芒萁的差异表达基因

以生长于广西大厂锡多金属矿上部(重金属胁迫区)和未受矿化或污染影响的矿区外围(对照区)的芒萁〔Dicranopteris pedata(Houtt.)Nakaike〕为实验材料,对芒萁叶片进行转录组高通量测序,并对组装得到的unigenes经NCBI官方非冗余蛋白质序列数据库(Nr)、NCBI官方非冗余核苷酸序列数据库(Nt)、KEGG直系同源数据库(KO)、Swiss-Prot数据库(Swiss-Prot)、蛋白质家族数据库(Pfam)、基因功能分类体系数据库(GO)和真核生物直系同源序列数据库(KOG)进行注释,同时分析重金属胁迫区和对照区芒萁叶片间的差异表达unigenes.结果显示:测序获得19.56 Gb clean data,其中,重金属胁迫区和对照区芒萁叶片分别含10.14和9.42 Gb clean data.组装得到的250582个unigenes中有120097个unigenes得到注释,占unigenes总数的47.93%.与对照区相比较,重金属胁迫区芒萁叶片中上调和下调差异表达unigenes分别有208和620个,其中120个上调差异表达unigenes注释为代谢过程,占所有上调差异表达unigenes的57.69%;285个下调差异表达unigenes注释为催化活性,占所有下调差异表达unigenes的45.97%.重金属胁迫区芒萁叶片中15个unigenes与重金属转运和耐受相关,其中c44988 g1和c84121 g1的相对表达量分别极显著和显著高于对照区.研究结果显示:芒萁响应自然金属矿化或矿山重金属污染的基因可以用于生物地球化学找矿和土壤重金属污染检测.