长角血蜱饥饿雌成蜱转录组测序及分析

Sequencing and analysis on transcriptomes of unfed female Haemaphysalis longicornis

目的利用IlluminaHiSeq高通量技术对长角血蜱饥饿雌成蜱进行转录组测序。方法将测序得到的数据进行拼接组装,并利用生物信息学方法对所得序列进行基因功能注释、功能分类、代谢途径分析和简单重复序列标记等分析。结果共获得181 246 184个clean reads数据,组装后得到107 428个unigenes序列,平均长度1 246.29。以BLAST软件将所有的unigenes序列与Nr,Nt,Pfam,KOG,Swiss-prot,KEGG,GO数据库进行比对。与Nr数据库比对发现长角血蜱基因序列与同科肩突硬蜱(Ixodes scapularis)具有较高的同源性,为55.3%;根据GO数据库的注释,将其分为生物过程、细胞组分和分子功能3大类共56分支;通过与KOG数据库的注释划分25类;根据与KEGG数据库的分析发现含有代谢通路相关的基因有32类,其中较多的是信号转导(signal transduction)类,占12.13%。通过SSR位点查找共鉴定出45 863个SSR;SNP分析发现发生转换的SNP数为195 369个,发生颠换的SNP数为96 780个。结论通过对长角血蜱饥饿雌成蜱转录组水平分析,可为后续长角血蜱功能基因的挖掘及表达谱等方面的研究奠定基础。

Objective To sequence transcriptomes of unfed female Haemaphysalis longicornis using IlluminaHiSeq high-throughput technology.Methods The data on sequences the transcriptomes were spliced and assembled,and the obtained sequences were analyzed with functional annotation,functional classification,metabolic pathway analysis and simple repeated sequence markers using bioinformatics methods.Results A total of 181 246 184 clean reads data were obtained and 107 428 unigene sequences were obtained after assembly,with an average length of 1 246.29.All unigene sequences were aligned with the Nr,Nt,Pfam,KOG,Swiss-prot,KEGG,GO databases using BLAST software.Compared with the Nr database,the long-horned blood scorpion gene sequence has a high homology(55.3%)with that of Ixodes scapularis.According to the annotation of the GO database,the functions of the all unigene sequences were divided into 3 categories(biological process,cellular component and molecular function)covering 56 functional groups;based on to the annotations of the KOG database,all the unigene sequences were assigned into 25 categories;while,according to the analysis of the KEGG database,there are 32 groups of genes involved in metabolic pathways and a major part of them(12.13%)are involved in signal transduction.A total of 45 863 simple sequence repeats(SSRs)were identified with SSR locus search.Single nucleotide polymorphism(SNP)analysis indicated that the number of SNPs for base transition was 195 369 and that for base transversion was 96 780.Conclusion The analysis on transcriptomes of unfed female Haemaphysalis longicornis lays a foundation for subsequent researches on gene expression and expression of the tick.