刺猎蝽属种间转录组比较分析

Interspecific comparative analysis of transcriptome in the genus Sclomina

【目的】建立刺猎蝽属Sclomina 3种的转录组数据库,分析近缘种间转录组表达差异,探讨在转录组水平的种间趋异情况。【方法】采用Illumina HiSeq TM 4000高通量测序平台进行刺猎蝽属齿缘刺猎蝽S.erinacea、兴仁刺猎蝽S.xingrensis和广西刺猎蝽S.guangxiensis 3种转录组测序;利用Nr,Swiss-Prot,COG/KOG和KEGG数据库进行基因功能注释;对种间差异表达基因(differentially expressed genes,DEGs)进行GO功能注释和KEGG通路富集分析。【结果】齿缘刺猎蝽、兴仁刺猎蝽和广西刺猎蝽转录组测序组装获得42215个unigenes。21117个基因在上述4个数据库中得到注释(各数据库注释基因数为:Nr,20522;Swiss-Prot,15550;COG/KOG,13969;KEGG,10850)。兴仁刺猎蝽vs齿缘刺猎蝽转录组有3390个DEGs(803个上调,2587个下调),兴仁刺猎蝽vs广西刺猎蝽转录组有12543个DEGs(6639个上调,5904个下调),齿缘刺猎蝽vs广西刺猎蝽转录组有3580个DEGs(2485个上调,1095个下调)。GO功能注释结果显示,上述DEGs被注释到生物学过程、细胞组分和分子功能3个功能大类和功能亚类,兴仁刺猎蝽vs齿缘刺猎蝽转录组中DEGs分别注释到代谢过程、细胞和细胞部分等50个功能亚类中,下调基因多,上调基因少;兴仁刺猎蝽vs广西刺猎蝽转录组中DEGs被注释到生物学过程、细胞部分和细胞等54个功能亚类中;齿缘刺猎蝽vs广西刺猎蝽转录组中DEGs被分别注释到细胞部分、代谢过程等46个功能亚类中。6254个DEGs被注释到KEGG代谢通路,涉及多个不同的代谢通路,包括代谢通路、药物代谢-细胞色素P450、核糖体等;兴仁刺猎蝽vs齿缘刺猎蝽转录组702个DEGs富集到164个KEGG通路;兴仁刺猎蝽vs广西刺猎蝽转录组中2091个DEGs富集得到201个KEGG通路;齿缘刺猎蝽vs广西刺猎蝽转录组中865个DEGs富集得到129个KEGG通路;兴仁刺猎蝽vs广西刺猎蝽转录组DEGs富集到代谢通路、药物代谢-细胞色素P450和核糖体等KEGG通路的数量最多。【结论】本研究结果丰富了刺猎蝽属的分子生物学数据;在转录组水平上呈现了齿缘刺猎蝽、兴仁刺猎蝽和广西刺猎蝽的种间差异,探索了该属物种分化的实际情况。

【Aim】To establish the transcriptome databases of three species of the genus Sclomina,analyze the expression differences of transcriptomes among related species,and investigate the interspecific divergence at the transcriptome level.【Methods】The Illumina HiSeq TM 4000 high-throughput sequencing platform was used for transcriptome sequencing in the three Sclomina species,S.erinacea,S.xingrensis and S.guangxiensis.The gene function annotation was conducted by using Nr,Swiss-Prot,COG/KOG and KEGG databases.The GO functional annotation and KEGG pathway enrichment analysis were performed for the interspecific differentially expressed genes(DEGs).【Results】A total of 42215 unigenes were obtained by transcriptome sequencing of S.erinacea,S.xingrensis and S.guangxiensis.A total of 21117 genes were annotated in the above four databases(number of genes:Nr,20522;Swiss-Prot,15550;COG/KOG,13969,KEGG,10850).There were 3390 DEGs(803 up-regulated,2587 down-regulated)in the S.xingrensis vs S.erinacea transcriptomes,12543 DEGs(6639 up-regulated,5904 down-regulated)in the S.xingrensis vs S.guangxiensis transcriptomes,and 3580 DEGs(2485 up-regulated,1095 down-regulated)in the S.erinacea vs S.guangxiensis transcriptomes.GO functional annotation result showed that all of the above DEGs were annotated into three functional classes including biological process,cellular component and molecular function,and some functional subclasses.The DEGs in the S.xingrensis vs S.erinacea transcriptomes were annotated into 50 functional subclasses,such as metabolic process,cell and cell part,and more genes were down-regulated,and a few genes were up-regulated.The DEGs in the S.xingrensis vs S.guangxiensis transcriptomes were annotated into 54 functional subclasses,such as biological process,cell part and cell.The DEGs in the S.erinacea vs S.guangxiensis transcriptomes were annotated into 46 functional subclasses,such as cell part and metabolic process.A total of 6254 DEGs were annotated to KEGG metabolic pathways involving several different metabolic pathways,including metabolic pathway,drug metabolism-cytochrome P450,ribosome,etc.A total of 702 DEGs in the S.xingrensis vs S.erinacea transcriptomes were enriched in 164 KEGG pathways.A total of 2091 DEGs in the S.xingrensis vs S.guangxiensis transcriptomes were enriched to 201 KEGG pathways.A total of 865 DEGs in the S.erinacea vs S.guangxiensis transcriptomes were enriched to 129 KEGG pathways.The number of DEGs enriched in metabolic pathway,drug metabolism-cytochrome P450 and ribosome in the S.xingrensis vs S.guangxiensis transcriptomes was the highest.【Conclusion】This study results enrich the molecular biological data of the genus Sclomina.At the level of transcriptome,the interspecies difference in S.erinacea,S.xingrensis and S.guangxiensis,and the actual situation of species divergence in this genus Sclomina are explored.