哈茨木霉Th-33全基因组序列测定与分析

doi:10.16409/j.cnki.2095-039x.2016.02.011

摘要/Abstract

摘要： 利用Hiseq2500高通量测序平台对哈茨木霉Th-33全基因组进行序列测定,获得196个scaffolds,共预测了10849个基因,平均长度为1776 bp(GenBank登录号:PRJNA272949)。以GO(gene ontology)数据库对预测出的基因做基因注释,共注释基因6238个;以KEGG(kyoto encyclopedia of genes and genomespathway database)数据库对预测出的基因做基因注释,有6789个基因注释到279条KEGG代谢途径。KEGG富集分析显示,对氨基苯甲酸甲酯降解代谢通路涉及基因最多,有232个基因;其次是双酚降解代谢通路,有206个基因。利用Rfam数据库对基因组序列进行RNA分类预测,共分为25个类别,包含7123个基因,其中涉及基因最多的为转录后修饰、蛋白质翻转和分子伴侣一类。比较了哈茨木霉、深绿木霉、绿木霉以及里氏木霉基因组中重寄生相关的碳水化合物活性酶、蛋白酶及次生代谢相关基因。研究结果有助于深入了解木霉菌的生防机制、推动木霉菌功能基因的挖掘和利用。

关键词: 哈茨木霉, 基因组, 基因注释, 代谢通路, 重寄生相关基因

Abstract: The whole genome of Trichoderma harzianum Th-33 was sequenced on a Hiseq2500 instrument. A total of 196 scaffolds were assembled and 10849 genes were predicted with an average length of 1776 bp (GenBank number:PRJNA272949). 6238 predicted genes were annotated in gene ontology (GO) and 6789 predicted genes were mapped to 279 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. KEGG enrichment analysis showed that the significantly enriched pathways were amino benzoate degradation pathway involving 232 genes and bisphenol degradation pathway involving 206 genes. Rfam was used for the RNA classification and predication of the genome. A total of 7123 genes were assigned to 25 categories, and the largest number of genes was involved in posttranslational modification, protein turnover, chaperones category. The mycoparasitism-related carbohydrate active enzymes, proteases and secondary metabolites-related genes in the T. harzianum, T. atroviride, T. virens and T. reesei genome were compared. The data in this study offer a better understanding of the biocontrol mechanism of Trichoderma and explore the functional genes in Trichoderma.

Key words: Trichoderma harzianum, whole genome, gene annotation, metabolic pathway, mycoparasitismrelated genes

中图分类号:

S476

孙青, 蒋细良, 庞莉, 王丽荣, 李梅. 哈茨木霉Th-33全基因组序列测定与分析[J]. 中国生物防治学报, DOI: 10.16409/j.cnki.2095-039x.2016.02.011.

SUN Qing, JIANG Xiliang, PANG Li, WANG Lirong, LI Mei. The Genome Sequence of Trichoderma harzianum Th-33[J]. journal1, DOI: 10.16409/j.cnki.2095-039x.2016.02.011.

参考文献

[1] Weindling R. Trichoderma lignorum as a parasite of other soil fungi[J]. Phytopathology, 1932, 22(8):837-845.
[2] 冀颖. 哈茨木霉Th-33ThChsC基因的克隆及功能初步分析[D]. 北京:中国农业科学院, 2011.
[3] Magill C. Current trends in fungal genomics and biology[J]. Fungal Genomics and Biology, 2012, 2(1):e103.
[4] Goffeau A, Barrell B G, Bussey H, et al. Life with 6000 genes[J]. Science, 1996, 274(5287):546, 563-567.
[5] Wood V, Gwilliam R, Rajandream M A, et al. The genome sequence of Schizosaccharomyces pombe[J]. Nature, 2002, 415(6874):871-880.
[6] Galagan J E, Calvo S E, Borkovich K A, et al. The genome sequence of the filamentous fungus Neurospora[J]. Nature, 2003, 422(6934):859-868.
[7] Martinez D, Berka R M, Henrissat B, et al. Genome sequencing and analysis of the biomass-degrading fungus Trichoderma reesei[J]. Nature Biotechnology, 2008, 26(5):553-560.
[8] Mukherjee P K, Horwitz B A, Herrera-Estrella A, et al. Trichoderma research in the genome Era[J]. Annual Review of Phytopathology, 2013, 51:105-129.
[9] Kubicek C P, Herrera-Estrella A, Seidl-Seiboth V, et al. Comparative genome sequence analysis underscores mycoparasitism as the ancestral life style of Trichoderma[J]. Genome Biology, 2011, 12(4):R40.
[10] Luo R, Liu B, Xie Y, et al. SOAPdenovo2:an empirically improved memory-efficient short-read de novo assembler[J]. GigaScience, 2012, 1(1):18.
[11] Smith A F A, Hubley R, Green P. Repeat Masker at http://repeatmasker.org.
[12] Stanke M, Steinkamp R, Waack S, et al. AUGUSTUS:a web server for gene finding in eukaryotes[J]. Nucleic Acids Research, 2004, 32(S2):309-312.
[13] Majoros W H, Pertea M, Salzberg S L. TigrScan and GlimmerHMM:two open-source ab initio eukaryotic gene-finders[J]. Bioinformatics, 2004, 20(16):2878-2879.
[14] Lukashin A, Borodovsky M. GeneMark.hmm:new solutions for gene finding[J]. Nucleic Acids Research, 1998, 26(4):1107-1115.
[15] Slater G S, Birney E. Automated generation of heuristics for biological sequence comparison[J]. BMC Bioinformatics, 2005, 6(2):31.
[16] Haas B J, Salzberg S L, Zhu W, et al. Automated eukaryotic gene structure annotation using evidence modeler and the program to assemble spliced alignments[J]. Genome Biology, 2008, 9(1):R7.
[17] Griffiths-Jones S, Bateman A, Marshall M, et al. Rfam:an RNA family database[J]. Nucleic Acids Research, 2003, 31(1):439-441.
[18] Nierman W C, Pain A, Anderson M J, et al. Genomic sequence of the pathogenic and allergenic filamentous fungus Aspergillus fumigatus[J]. Nature, 2005, 438(7071):1151-1156.
[19] Pel H J, de Winde J H, Archer D B, et al. Genome sequencing and analysis of the versatile cell factory Aspergillus niger CBS 513.88[J]. Nature Biotechnology, 2007, 25(2):221-231.
[20] Berka R M, Grigoriev I V, Otillar R, et al. Comparative genomic analysis of the thermophilic biomass-degrading fungi Myceliophthora thermophila and Thielavia terrestris[J]. Nature Biotechnology, 2011, 29(10):922-927.
[21] Martinez D, Larrondo L F, Putnam N, et al. Genome sequence of the lignocellulose degrading fungus Phanerochaete chrysosporium strain RP78[J]. Nature Biotechnology, 2004, 22(6):695-700.
[22] Espagne E, Lespinet O, Malagnac F, et al. The genome sequence of the model ascomycete fungus Podospora anserina[J]. Genome Biology, 2008, 9(5):R77.
[23] Zhou P, Zhang G, Chen S, et al. Genome sequence and transcriptome analyses of the thermophilic zygomycete fungus Rhizomucor miehei[J]. BMC Genomics, 2014, 15(4):294.
[24] Ma L J, Ibrahim A S, Skory C, et al. Genomic analysis of the basal lineage fungus Rhizopus oryzae reveals a whole-genome duplication[J]. PLoS Genetics, 2009, 5(7):e1000549.
[25] McHunu N P, Permaul K, Abdul Rahman A Y, et al. Xylanase superproducer:genome sequence of a compost-loving thermophilic fungus, Thermomyces lanuginosus strain SSBP[J]. Genome Announcements, 2013, 1(3):e00388-13.
[26] Berka R M, Grigoriev I V, Otillar R, et al. Comparative genomic analysis of the thermophilic biomass-degrading fungi Myceliophthora thermophila and Thielavia terrestris[J]. Nature Biotechnology, 2011, 29(10):922-927.
[27] Zhao Z, Liu H, Wang C, et al. Correction:Comparative analysis of fungal genomes reveals different plant cell wall degrading capacity in fungi[J]. BMC Genomics, 2014, 15(1):6.
[28] Cuomo C A, Güldener U, Xu J R, et al. The Fusarium graminearum genome reveals a link between localized polymorphism and pathogen specialization[J]. Science, 2007, 317(5843):1400-1402.
[29] Črešnar B, Petrič Š. Cytochrome P450 enzymes in the fungal kingdom[J]. Biochim Biophys Acta, 2011, 1814(1):29-35.
[30] 宋晓妍, 张玉忠, 王元秀. 木霉peptaibols抗菌肽的研究进展[J]. 微生物学报, 2011, 51(4):438-444.
[31] Mueller O, Kahmann R, Aquilar G, et al. The secretome of the maize pathogen Ustilago maydis[J]. Fungal Genetics Biology, 2008, 45(S1):63-70.
[32] Doehlemann G, van der Linde K, Assmann D, et al. Pep1, a secreted effector protein of Ustilago maydis, is required for successful invasion of plant cells[J]. PLoS Pathogens, 2009, 5(2):e1000290.
[33] Seidl V, Song L, Lindquist E, et al. Transcriptomic response of the mycoparasitic fungus Trichoderma atroviride to the presence of a fungal prey[J]. BMC Genomics, 2009, 10(4):567.
[34] Schoberle T J, Nguyen-Coleman C K, Herold J, et al. A novel C₂H₂ transcription factor that regulates gliA expression interdependently with GliZ in Aspergillus fumigates[J]. PLoS Genetics, 2014, 10(5):e1004336.
[35] Nitta M, Furukawa T, Shida Y, et al. A new Zn(Ⅱ)(2)Cys(6)-type transcription factor BglR regulates β-glucosidase expression in Trichoderma reesei[J]. Fungal Genetics Biology, 2012, 49(5):388-397.