Differential transcriptome analysis of Sporocytophaga sp. CX11 and identification of candidate genes involved in lignocellulose degradation

Abstract

Cellulose is the most abundant renewable bioresources on earth, and the biodegradation and utilization of cellulose would contribute to the sustainable development of global environment. Sporocytophaga species are common aerobic cellulose-degrading bacteria in soil, which can adhere to the surface of cellulose matrix and motile by gliding. In this study, a differential transcriptome analysis of Sporocytophaga sp. CX11 was performed and a total of 4,217 differentially expressed genes (DEGs) were identified. Gene Ontology enrichment results showed that there are three GO categories related to cellulose degradation function among the annotated DEGs. A total of 177 DEGs were identified as genes encoding carbohydrate-active enzymes (CAZymes), among which 54 significantly upregulated CAZymes were mainly cellulases, hemicellulases, pectinases, etc. 39 DEGs were screened to associate with gliding function. In order to explore unannotated genes potentially related to cellulose metabolism, cluster analysis was performed using the Short-Time Series Expression Miner algorithm (STEM). 281 unannotated genes were predicted to be associated with the initial-middle stage of cellulose degradation and 289 unannotated genes might function in the middle-last stage of cellulose degradation. Sporocytophaga sp. CX11 could produce extracellular endo-xylanase, endo-glucanase, FPase and β-glucosidase, respectively, according to different carbon source conditions. Altogether, this study provides valuable insights into the transcriptome information of Sporocytophaga sp. CX11, which would be useful to explore its application in biodegradation and utilization of cellulose resources.

Keywords

Cellulose biodegradation / Sporocytophaga / Differential transcriptome analysis / Carbohydrate-active enzymes / Gliding motility

References

[1] Anderson K, Salyers A. Biochemical evidence that starch breakdown by Bacteroides thetaiotaomicron involves outer membrane starch-binding sites and periplasmic starch-degrading enzymes. J Bacteriol, 1989, 171(6): 3192-3198.

[2] Bayer EA, Belaich JP, Shoham Y, Lamed R. The cellulosomes: multienzyme machines for degradation of plant cell wall polysaccharides. Annu Rev Microbiol, 2004, 58: 521-554.

[3] Biely P. Microbial carbohydrate esterases deacetylating plant polysaccharides. Biotechnol Adv, 2012, 30(6): 1575-1588.

[4] Boraston AB, Bolam DN, Gilbert HJ, Davies GJ. Carbohydrate-binding modules: fine-tuning polysaccharide recognition. Biochem J, 2004, 382: 769-781.

[5] Burchard PR. Gliding motility of prokaryotes: ultrastructure, physiology, and genetics. Annu Rev Microbiol, 1981, 35(35): 497-529.

[6] Chen X, Mao X, Huang J, Yang D, . KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res, 2011, 39: 316-322.

[7] Drula E, Garron M, Dogan S, Lombard V, Henrissat B, Terrapon N. The carbohydrate-active enzyme database: functions and literature. Nucleic Acids Res, 2022, 50: D571-D577.

[8] Ernst J, Nau GJ, Bar-Joseph Z. Clustering short time series gene expression data. Bioinformatics, 2005, 21: i159-i168.

[9] Fang W, Song R, Zhang X, Zhang X, Zhang X, Wang X, Fang Z, Xiao Y. Characterization of a novel beta-glucosidase from Gongronella sp. W5 and its application in the hydrolysis of soybean isoflavone glycosides. J Agric Food Chem, 2014, 62(48): 11688-11695.

[10] Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, Heger A, Hetherington K, Holm L, Mistry J, Sonnhammer EL, Tate J, Punta M. Pfam: the protein families database. Nucleic Acids Res, 2014, 42: D222-230.

[11] Gaskell J, Blanchette RA, Stewart PE, BonDurant SS, Adams M, Sabat G, Kersten P, Cullen D. Transcriptome and secretome analyses of the wood decay fungus Wolfiporia cocos support alternative mechanisms of lignocellulose conversion. Appl Environ Microbiol, 2016, 82(13): 3979-3987.

[12] Gouvêa PD, Bernardi AV, Gerolamo LE, Emerson D, Riano-Pachón D, Uyemura SA, Dinamarco TM. Transcriptome and secretome analysis of Aspergillus fumigatus in the presence of sugarcane bagasse. BMC Genomics, 2018, 19(1): 232.

[13] Guerrero EB, de Villegas RMD, Soria MA, Santangelo MP, Campos E, Talia PM. Characterization of two GH5 endoglucanases from termite microbiome using synthetic metagenomics. Appl Microbiol Biotechnol, 2020, 104(19): 8351-8366.

[14] Holt SC, Leadbetter ER. Fine structure of Sporocytophaga myxococcoides. Arch Mikrobiol, 1967, 57(3): 199-213.

[15] Hutchinson HB, Clayton J. On the decomposition of cellulose by an aerobic organism (Spirochaeta cytophaga, n. sp.). J Agric Sci, 1919, 9(02): 143.

[16] Lane S, Dong J, Jin YS. Value-added biotransformation of cellulosic sugars by engineered Saccharomyces cerevisiae. Bioresour Technol, 2018, 260: 380-394.

[17] Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods, 2012, 9(4): 357-359.

[18] Liu L, Gao P, Chen G, Wang L. Draft genome sequence of cellulose-digesting bacterium Sporocytophaga myxococcoides PG-01. Genome Announc, 2014, 2(6): e01154-e1214.

[19] Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods, 2001, 25(4): 402-408.

[20] Macconaill LE, Burns RT, Anwesha N, Coleman HA, Slevin MK, Kristina G, Madelyn L, Kevin L, Mirna J, Mcneill MS. Unique, dual-indexed sequencing adapters with UMIs effectively eliminate index cross-talk and significantly improve sensitivity of massively parallel sequencing. BMC Genomics, 2018, 19(1): 30.

[21] Mcbride MJ, Xie G, Martens EC, Lapidus A, Henrissat B, Rhodes RG, Goltsman E, Wang W, Xu J, Hunnicutt DW. Novel features of the polysaccharide-digesting gliding bacterium Flavobacterium johnsoniae as revealed by genome sequence analysis. Appl Environ Microbiol, 2009, 75(21): 6864-6875.

[22] Miller GL. Use of dinitrosalicylic acid reagent for determination of reducing sugar. Anal Chem, 1959, 31(3): 420-428.

[23] Rhodes RG, Samarasam MN, Groll EV, Mcbride MJ. Mutations in Flavobacterium johnsoniae sprE result in defects in gliding motility and protein secretion. J Bacteriol, 2011, 193(19): 5322.

[24] Sato K, Naya M, Hatano Y, Kondo Y, Sato M, Narita Y, Nagano K, Naito M, Nakayama K, Sato C. Colony spreading of the gliding bacterium Flavobacterium johnsoniae in the absence of the motility adhesin SprB. Sci Rep, 2021, 11(1): 967.

[25] Stanier RY. The cytophaga group: a contribution to the biology of mycobacteria. Bacteriol Rev, 1942, 6(3): 143-196.

[26] Suen G, Weimer PJ, Stevenson DM, Aylward FO, Boyum J, Deneke J, Drinkwater C, Ivanova NN, Mikhailova N, Chertkov O, Goodwin LA, Currie CR, Mead D, Brumm PJ. The complete genome sequence of Fibrobacter succinogenes S85 reveals a cellulolytic and metabolic specialist. PLoS ONE, 2011, 6(4): e18814.

[27] Sukumaran RK, Christopher M, Kooloth-Valappil P, Sreeja-Raju A, Mathew RM, Sankar M, Puthiyamadam A, Adarsh VP, Aswathi A, Rebinro V, Abraham A, Pandey A. Addressing challenges in production of cellulases for biomass hydrolysis: targeted interventions into the genetics of cellulase producing fungi. Bioresour Technol, 2021

[28] Taillefer M, Arntzen M, Henrissat B, Pope PB, Larsbrink J. Proteomic dissection of the cellulolytic machineries used by soil-dwelling bacteroidetes. mSystems, 2018, 3: e00240.

[29] Wang L, Feng Z, Wang X, Wang X, Zhang X. DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics, 2010, 26(1): 136-138.

[30] Wilson DB. Evidence for a novel mechanism of microbial cellulose degradation. Cellulose, 2009, 16(4): 723-727.

[31] Xie G, Bruce DC, Challacombe JF, Chertkov O, Detter JC, Gilna P, Han CS, Lucas S, Misra M, Myers GL, Richardson P, Tapia R, Thayer N, Thompson LS, Brettin TS, Henrissat B, Wilson DB, McBride MJ. Genome sequence of the cellulolytic gliding bacterium Cytophaga hutchinsonii. Appl Environ Microbiol, 2007, 73(11): 3536-3546.

[32] Yang Q, Guo B, Sun H, Zhang J, Liu S, Hexige S, Yu X, Wang X. Identification of the key genes implicated in the transformation of OLP to OSCC using RNA-sequencing. Oncol Rep, 2017, 37(4): 2355-2365.

[33] Yang M, Zhao J, Yuan Y, Chen X, Yang F, Li X. Comparative metagenomic discovery of the dynamic cellulose-degrading process from a synergistic cellulolytic microbiota. Cellulose, 2021, 28(4): 2105-2123.

[34] Ye J, Fang L, Zheng H, Zhang Y, Chen J, Zhang Z, Wang J, Li S, Li R, Bolund L, Wang J. WEGO: a web tool for plotting GO annotations. Nucleic Acids Res, 2006, 34: W293-297.

[35] Young M, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol, 2010, 11(2): R14.

[36] Zhang HP, Himmel ME, Mielenz JR. Research review paper outlook for cellulase improvement: screening and selection strategies. Biotechnol Adv, 2006, 24(5): 452-481.

[37] Zhao D, Wang Y, Wang S, Zhang W, Qi Q, Lu X. A disulfide oxidoreductase (CHU_1165) is essential for cellulose degradation by affecting outer membrane proteins in cytophaga hutchinsonii. Appl Environ Microbiol, 2020

[38] Zhao Y, Li M-C, Konaté MM, Chen L, Das B, Karlovich C. TPM, FPKM, or normalized counts? A comparative study of quantification measures for the analysis of RNA-seq data from the NCI patient-derived models repository. J Transl Med, 2021, 19(1): 269.

[39] Zheng Q, Wang X. GOEAST: a web-based software toolkit for Gene Ontology enrichment analysis. Nucleic Acids Res, 2008, 36(S2): W358-W363.

[40] Zhu Y, McBride MJ. Deletion of the Cytophaga hutchinsonii type IX secretion system gene sprP results in defects in gliding motility and cellulose utilization. Appl Microbiol Biotechnol, 2014, 98(2): 763-775.

[41] Zhu Y, McBride MJ. The unusual cellulose utilization system of the aerobic soil bacterium Cytophaga hutchinsonii. Appl Microbiol Biotechnol, 2017, 101(19): 7113-7127.

Funding

National Natural Science Foundation of China(32072160)

Science and Technology Department of Liaoning(J2020041)

Key Research Projects of The Educational Department of Liaoning Province(LJKZZ20220061)