基于高通量测序的大花序桉顶芽转录组生物信息学分析  

Transcriptome Sequencing and Bioinformatic Analysis of Eucalyptus cloeziana Terminal Buds

作　　者：兰俊 蒋维昕[2] 张磊 梁馨元 白天道[2] Lan Jun;Jiang Weixin;Zhang Lei;Liang Xinyuan;Bai Tiandao(Dongmen Forest Farm of Guangxi Zhuang Autonomous Region,Chongzuo,532108;Forestry College of Guangxi University,Nanning,530000)

机构地区：[1]广西壮族自治区国有东门林场,崇左532108 [2]广西大学林学院,南宁530000

出　　处：《分子植物育种》2021年第16期5342-5351,共10页Molecular Plant Breeding

基　　金：广西重点研发计划项目(2018AB44025);广西大学科研项目“大花序桉人工林集成技术研究”(20180493)共同资助。

摘　　要：为了获得珍贵用材树种大花序桉顶芽转录组数据及预测关键基因功能,本研究基于Illumina HiSeq X Ten测序技术获得大花序桉顶芽转录组原始数据,经Trinity组装拼接获得高质量Unigene,并与NR、Swiss-Prot、GO、KOG、egg NOG和KEGG等生物信息数据库进行序列比对和功能注释,利用MISA软件进行SSR位点搜索和分析。从大花序桉顶芽中共获得26 587条高质量Unigene,平均长度为1 279.69 bp;共有22 099条Unigene至少在一个数据库中被成功注释,其中,11 507条Unigene被注释到KOG数据库中25个功能类别,以参与一般功能基因的数量最多;GO数据库中,所注释到的14 105条Unigene分别匹配到生物功能、细胞组分和分子功能3大类50个功能基因区,其中执行生物过程所占比例最多;KEGG功能注释共发现有7 117个Unigene参与127条代谢通路,以代谢相关的基因最丰富;共有1 021条Unigene注释到转录因子数据库,分布于65个家族,其中比例最大的是bHLH和MYB家族;3 274条Unigene注释到植物抗性基因数据库,分布于13个类别,相匹配基因数量最大的是RLP和TNL。MISA软件共检测到12 366个SSR位点,分布密度为1/2.75 kb,重复基元类型丰富,标记开发潜力大。本研究利用高通量测序获得丰富的顶芽转录组信息,可以为大花序桉分子辅助育种提供丰富的资源。In order to obtain transcriptome data and predict the key gene function of Eucalyptus cloeziana, an important timber species in China, the Illumina HiSeq X Ten sequencing technology was conducted to carry out transcriptome sequencing of E. cloeziana terminal buds. The high-quality Unigene was assembled and spliced by Trinity, and sequence alignment and function annotation were performed with biological information databases such as NR, Swiss-Prot, GO, KOG, eggNOG and KEGG etc. SSR sites search and analysis was performed by MISA software. A total of 26 587 high quality Unigenes were obtained from the terminal bud of E. cloeziana with an aver age length of 1 279.69 bp. A total of 22 099 Unigenes were successfully annotated in at least one biological database.Of these, 11 507 Unigenes were successfully annotated with 25 functions in KOG database, and the most common function was general functional genes prediction. In the GO database, the 14 105 Unigenes annotated were matched to 50 functional gene groups in 3 categories of biological function, cell component and molecular function, respectively. Among them, the biological processes accounted for the largest proportion. Through KEGG pathways analysis, 7 117 Unigenes were successfully annotated and 127 metabolic pathways were detected, with the most abundant metabolism-related genes. Moreover, 1 021 Unigenes annotations were assigned with 65 families of transcription factor(TF) database, among which the b HLH and MYB had the largest proportion. A total of3 274 Unigenes were annotated in PRG database and formed into 13 resistant gene categories, among which RLP and TNL had the largest number of matched genes. In addition, 12 366 SSR sites were detected by MISA software,with a distribution density of 1/2.75 kb and a variety of repeating primitive types. In this study, abundant transcriptome information from terminal bud was obtained by using high-throughput sequencing, which was beneficial for molecular assisted breeding of E. cloeziana.

关 键 词：大花序桉 Illumina HiSeq X Ten 转录组 基因注释 顶芽 

分 类 号：S792.39[农业科学—林木遗传育种]