柳杉维管形成层及木质部区转录组序列及基因表达分析

Transcriptome Sequence and Gene Expression in Vascular Cambium and Xylem of Cryptomeria fortunei

柳杉是中国南方重要的针叶用材树种,具有树形高大,纹理直等特点,已广泛用于板材和建筑生产中,研究柳杉木材形成过程中维管形成层及木质部区转录组特征,为木材形成主要过程的分子机理,木材形成有关基因调控,培育优良材质的林木提供理论参考。本研究以柳杉形成层组织为材料,通过Illumina Hiseq测序平台对4个不同发育阶段的柳杉维管形成层进行转录组测序;对Unigene进行了蛋白功能注释、分类及KEGG代谢通路分析等。测序数据参照无参转录组测序流程进行分析,一共产生105.9 Gb的数据,过滤后Clean reads均达到90%以上。Clean reads经Trinity软件进行组装,一共产生64 969个Unigene,对Unigene进行拼接,拼接后的Contigs序列有29 381条、Singletons有35 588条,总长度为83 003 836 bp。将Unigen与七大功能数据库(NR, NT, GO, COG, KEGG, Swissprot, Interpro)进行比对,共有42 836个Unigene得到注释,占总Unigene的66.05%。GO注释共有89 644个转录本得到注释。对GO注释转录本进行分类,其中有38 432个转录本(42.87%)注释为生物过程,有32 749个转录本(36.53%)注释为细胞组分,有18 463个转录本(20.6%)注释为分子功能。KEGG注释有31 580个Unigene得到注释。分为6大类为:细胞代谢过程、环境信息处理、遗传信息处理、人类疾病、代谢作用、生物体系统。其中代谢作用涉及的Unigene最多达18 580个,占KEGG总注释Unigene的58.83%。另有37 762个Unigene得到COG注释,被分为25类,其中注释数目最多的类型是仅预测一般功能的Unigene,占总数的16.9%;其次是转录,比例为8.8%;再次是复制、重组及修复,占总数的7.8%。本研究最后分析了可能参与木材形成重要功能基因,为探讨木材发育的分子机制及进行重要功能基因的挖掘提供了提供了依据。

Cryptomeriafortunei is an important coniferous timber tree in southern China,which has the character- istics of tall shape and straight texture.It has been widely used in plate and building material production.The characteristics of the transcriptome ofvascular cambium and xylem zone were studied and analyzed,which would lay a theoretical basis for further study of the molecular mechanism and genetic regulation during the process of wood formation,and the cultivation of fine woods.In this study,cambium tissue of C.fortune i was used as test materials,and the transcriptome sequencing of C.fortunei vascular cambium at four different development stages was conducted based on Illumina Hiseq sequencing platform.The protein functional annotation,GO classification and KEGG metabolic pathway analysis of Unigene were performed.The sequencing data were analyzed according to the procedure of no-reference transcriptome sequencing.A total of 105.9 Gb data were obtained and the filtered clean reads reached more than 90%.Clean reads were assembled by Trinity software and a total of 64 969 Unigenes (29 381 contigs and 35 588 35 588 singletons)were obtained.The total length was 83 003 836 bp.42 914 U- nigenes (66.05%of the total Unigenes)were functionally annotated after all Unigenes aligned with seven functional databases (NR,NT,GO,COG,KEGG,Swissprot,Interpro).GO annotation showed that 89 644 transcripts were divided into 3 groups,which contained Biological Process (38 432 transcripts,42.87%),Cellular Compo- nent (32 749 transcripts,36.53%)and Molecular Function (18 463 transcripts,20.6%).31 580 Unigenes were annotated in KEGG database,which were divided into six groups Cellular Process,Environmental Information Processing,Genetic Information Processing,Human Disease,Metabolism,Organismal System.18 580 Unigenes were involved in metablism at most,accounting for 58.83%of the total Unigenes annotated in K.EGG database.37 762 Unigenes were annotated in COG database,which were divided into 25 groups.The group of general function prediction only (16.9%)was the biggest,the second one was transcription (8.8%),the latter ones were replication,re- combination and repair (7.8%).With the help of high-through sequencing,In this study,the possible functional genes involved in wood formation were analyzed,which might provide molecular data for exploring the molecular mechanism of wood development and further discovery of important functional genes.