Metagenomic insights into seasonal variations in the soil microbial community and function in a Larix gmelinii forest of Mohe, China

Abstract

The effect of seasons on the soil microbiome in a Larix gmelinii forest of Mohe, China, where winter temperatures are generally below − 40 °C, was evaluated with metagenomics analysis. Taxonomic profiling using sequencing information revealed that Proteobacteria, Actinobacteria, Acidobacteria and Verrucomicrobia were the dominant phyla in spring, summer, and fall, as were Bradyrhizobium, Chthoniobacter, Streptomyces, Acid Candidatus Koribacter at the genus level. Some species that were abundant in spring and fall greatly diminished in abundance in summer. Clusters of orthologous groups (COG) of proteins, carbohydrate-active enzymes (CAZy), Kyoto Encyclopedia of Genes and Genomes (KEGG) and NCBI databases were used to elucidate the function of diverse proteins and metabolites of the microbial community of L. gmelinii forest. COG analysis showed that fewer genes were detected in spring than in fall and summer, indicating that many soil microbes in the L. gmelinii forest were not tolerant to cold. Based on KEGG analysis, some pathways in the soil microbes were activated in spring and autumn and deactivated in summer. CAZy analysis revealed that most CAZy were more active in summer than in spring or autumn and were severely inhibited in the spring. Many functional pathways, proteins, and CAZy involved in the community changes were concerned with cold or heat resistance. Therefore, the soil in the L. gmelinii forest can be a valuable resource for further research on heat and cold tolerance of soil microbes.

Keywords

Soil microbial community / Microbial function / Metagenomic / Seasonal variation

References

[1] Abanda-Nkpwatt D, Musch M, Tschiersch J, Boettner M, Schwab W. Molecular interaction between Methylobacterium extorquens and seedlings: growth promotion, methanol consumption, and localization of the methanol emission site. J Exp Bot, 2006, 57(15): 4025-4032.

[2] Albertsen M, Hansen LB, Saunders AM, Nielsen PH, Nielsen KL. A metagenome of a full-scale microbial community carrying out enhanced biological phosphorus removal. ISME J, 2012, 6(6): 1094-1106.

[3] Altschul SF. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res, 1997, 25(17): 3389-3402.

[4] Apostel C, Bore EK, Halicki S, Kuzyakov Y, Dippold M (2015) Carbon metabolism of soil microorganisms at low temperatures: position-specific ¹³C labeled glucose reveals the story. AGU fall meeting Ast B33D-0724

[5] Bardgett RD, Wh VDP. Belowground biodiversity and ecosystem functioning. Nature, 2014, 515(7528): 505-511.

[6] Batut J, Andersson SGE, O’Callaghan D. The evolution of chronic infection strategies in the Alphaproteobacteria. Nat Rev Microbiol, 2004, 2(12): 933-945.

[7] Bhardwaj AK, Vinothkumar K, Rajpara N. Bacterial quorum sensing inhibitors: attractive alternatives for control of infectious pathogens showing multiple drug resistance. Recent Pat Anti-Infect Drug Disc, 2013, 8(1): 68-83.

[8] Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics, 2014, 30(15): 2114-2120.

[9] Buchfink B, Xie C, Huson D. Fast and sensitive protein alignment using DIAMOND. Nat Methods, 2015, 12(1): 59-60.

[10] Chen HB, Boutros PC. VennDiagram: a package for the generation of highly-customizable Venn and Euler diagrams in R. BMC Bioinform, 2011, 12: 35.

[11] Chen S, Zhang CZ, Liu DB, Zhang ZY, Yang JC, Wang ZX. Seasonal variation in the biomass of soil decomposer microbes and its relationship to the soil habitat in the Leymus Chinensis grasslands in northeast China. Acta Ecol Sin, 1993, 15(1): 91-94. (in Chinese with English abstract)

[12] Chun JY, Sexton GL, Roth LE, Stacey G. Identification and characterization of a novel Bradyrhizobium japonicum gene involved in host-specific nitrogen fixation. J Bacteriol, 1994, 176(21): 6717-6729.

[13] Dijkstra P, Dalder JJ, Selmants PC, Hart SC, Koch GW, Schwartz E, Hungate BA. Modeling soil metabolic processes using isotopologue pairs of position-specific ¹³C-labeled glucose and pyruvate. Soil Biol Biochem, 2011, 43(9): 1848-1857.

[14] Dong LH, Yang JS, Yuan HL. Research advances in molecular ecology of ammonia oxidizing bacteria. Chin J Appl Ecol, 2008, 19(06): 1381-1388.

[15] Eddy SR. A new generation of homology search tools based on probabilistic inference. Gen Inf, 2009, 23(1): 205-211.

[16] Fan YZ, Sun HQ. The cold-adapted molecular mechanism of cold-adapted microorganism. J Hengshui Univ, 2008, 10(1): 70-72. (in Chinese with English abstract)

[17] Feng BP, Gao RH, Zhang QL, Zhao QM. Study of dynamic variation of soil microorganisms under different management modes in Larix gmelinii forest. J Inner Mong Agric Univ, 2009, 30(04): 74-79. (in Chinese with English abstract)

[18] Feng G, Xie T, Wang X, Bai J, Zhao Y. Metagenomic analysis of microbial community and function involved in cd-contaminated soil. BMC Mirobiol, 2018 18 1 11

[19] Gao X, Zhang DY. Tree-ring based reconstruction of Mohe precipitation in January for nearly a hundred years. J Northwest Univ (Nat Sci Ed), 2018, 48(05): 749-754 + 760. (in Chinese with English abstract)

[20] Gi KN. To analyze the nutritional types of rhodopillum aureus. Biol Teach Middle Sch, 2004 01 Z1 28

[21] Huerta-Cepas J, Szklarczyk D, Forslund K, Cook H, Heller D, Walter MC, Rattei T, Mende DR, Sunagawa S, Kuhn M, Jensen LJ, Cv Mering, Bork P. eggNOG 4.5: a hierarchical orthology framework with improved functional annotations for eukaryotic, prokaryotic and viral sequences. Nucleic Acids Res, 2016, 44(D1): D286-D293.

[22] Janssen PH. Identifying the dominant soil bacterial taxa in libraries of 16S rRNA and 16S rRNA genes. Appl Environ Micrbiol, 2006, 72(3): 1719-1728.

[23] Kanehisa M, Goto S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res, 2000, 28(1): 27-30.

[24] Kanehisa M, Sato Y, Morishima K. BlastKOALA and GhostKOALA: KEGG tools for functional characterization of genome and metagenome sequences. J Mol Biol, 2016, 428(4): 726-731.

[25] Kikani BA, Sharma AK, Singh SP. Metagenomic and culture-dependent analysis of the bacterial diversity of a hot spring reservoir as a function of the seasonal variation. Int J Environ Res, 2017, 11(1): 25-38.

[26] Klein HL, Kreuzer KN. Replication, recombination, and repair. Mol Cell, 2002, 9(3): 471-480.

[27] Lange M, Eisenhauer N, Sierra CA, Holger B, Christoph E, Griffiths RI, Mellado-Vázquez PG, Malik AA, Jacques R, Stefan S, Steinbeiss S, Thomson BC, Trumbore SE, Gleixner G. Plant diversity increases soil microbial activity and soil carbon storage. Nat Commun, 2015, 6: 6707.

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

[29] Lee KJ, Jung YC, Park SJ, Lee KH. Role of heat shock proteases in quorum-sensing-mediated regulation of biofilm formation by vibrio species. mBio, 2018 9 1 e02086-17

[30] Li W, Godzik A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Boinformatics, 2006, 22(13): 1658-1659.

[31] Li H, Handsaker B, Wysoker A, Fennell T, Ruan J. The sequence alignment/map format and SAMtools. Boinformatics, 2009, 25(16): 2078-2079.

[32] Liu Y, Guo J, Hu G, Zhu H. Gene prediction in metagenomic fragments based on the SVM algorithm. BMC Bioinform, 2013 14 Suppl 5 S12

[33] Lombard V, Golaconda RH, Drula E, Coutinho PM, Bernard H. The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res, 2014, 42(Database issue): D490-D495.

[34] Lu XM, Chen C, Zheng TL. Metagenomic insights into effects of chemical pollutants on microbial community composition and function in estuarine sediments receiving polluted river water. Microb Ecol, 2017, 73(4): 791-800.

[35] Ma DL, Li SS, Wang LL, Shuying Z. Seasonal discrepancy of metabolic function of soil microbial communities in Zhalong wetland. China Sciencepap, 2016, 11(09): 1046-1050.

[36] Miah S, Dey S, Haque SMS. Shifting cultivation effects on soil fungi and bacterial population in Chittagong Hill Tracts, Bangladesh. J For Res, 2010, 21(3): 311-318.

[37] Miller JR, Koren S, Sutton G. Assembly algorithms for next-generation sequencing data. Genomics, 2010, 95(6): 315-327.

[38] Peng Y, Leung HCM, Yiu SM, Chin FYL. IDBA-UD: a de novo assembler for single-cell and metagenomic sequencing data with highly uneven depth. Bioinformatics, 2012, 28(11): 1420-1428.

[39] Rampelotto PH. Biotechnology of extremophiles, 2016, Basel: Springer 2

[40] Rivkina EM, Friedmann EI, McKay CP, Gilichinsky DA. Metabolic activity of permafrost bacteria below the freezing point. Appl Environ Microbiol, 2000, 66(8): 3230-3233.

[41] Rothschild LJ, Mancinelli R. Life in extreme environment. Nature, 2001, 409(6823): 1092-1101.

[42] Rovira P, Vallejo VR. Labile and recalcitrant pools of carbon and nitrogen in organic matter decomposing at different depths in soil: an acid hydrolysis approach. Geoderma, 2002, 107(1–2): 109-141.

[43] Salam LB, Obayori SO, Nwaokorie FO, Suleiman A, Mustapha R. Metagenomic insights into effects of spent engine oil perturbation on the microbial community composition and function in a tropical agricultural soil. Environ Sci Pollut Res, 2017, 24(8): 7139-7159.

[44] Sangwan P, Chen X, Hugenholtz P, Janssen PH. Chthoniobacter flavus gen. nov. sp. nov. the first pure-culture representative of subdivision two, Spartobacteria classis nov. of the Phylum Verucomicrobia. Appl Environ Microbiol, 2004, 70(10): 5875-5881.

[45] Shannon CE. A mathematical theory of communication. Bell Syst Tech J, 1948, 27(3): 379-423.

[46] Shi N, Zhang LP, Lv ZT, Runkei Y, Zhiheng L. Study on the herbicide produced by a strain of Streptomycetes and systematics of the strain. Acta Microbiol Sin, 2004, 44(2): 148-152.

[47] Smart CC, Fleming AJ. Hormonal and environmental regulation of a plant PDR5 like ABC transporter. J Biol Chem, 1996, 271(32): 19351-19357.

[48] Soderberg MA, Rossier O, Cianciotto NP. The type II protein secretion system of legionella pneumophila promotes growth at low temperatures. J Bacteriol, 2004, 186(12): 3712-3720.

[49] Steven B, Léveillé R, Pollard WH, Whyte LG. Microbial ecology and biodiversity in permafrost. Extremophiles, 2006, 10(4): 259-267.

[50] Tatusov RL, Galperin MY, Natale DA, Koonin EV. The COG database: a tool for genome-scale analysis of protein functions and evolution. Nucleic Acids Res, 2000 28 1 33

[51] Wang KQ, Bai XB. Classification of wood surface texture based on Gauss-MRF model. J For Res, 2006, 17(1): 57-61.

[52] Wang YF, Yu SX, Liu WQ. A new species diversity index and ITS fractal analysis. Acta Phytoecol Sin, 2002, 26(4): 391-395.

[53] Wang NF, Tao Z, Fang Z, Wang ET, He JF. Diversity and structure of soil bacterial communities in the Fildes region (maritime Antarctica) as revealed by 454 pyrosequencing. Front Microbiol, 2015, 6: 1188.

[54] Williams KP, Sobral BW, Dickerman AW. A robust species tree for the Alphaproeobacteria. J Bacteriol, 2007, 189(13): 4578-4586.

[55] Wu B, Feng K, Zhi XY, He Q, He Z. Progresses in environmental microbiome diversity and function research. Acta Sci Nat Univ Sunyatseni, 2017, 56(5): 1-11. (in Chinese with English abstract)

[56] Xing P, Kong FX, Gao G. Phylogenetic diversity and seasonal variation of bacterioplankton in Lake Taihu. J Lake Sci, 2007, 19(04): 373-381.

[57] Xu GH, Zheng HY, Zhang DS, Lu YB, Li YK, Zhang SX, Liu ZZ, Wu WF. Seasonal dynamics of soil microbial ecological distribution under three main forest types in Changbai Mounain. J For Res, 1985, 21(3): 286-291.

[58] Yan FJ, Tian XL, Dong SL, Gang Y. Seasonal variation of functional diversity of aquatic microbial community in Apostichopus japoncus cultural pond. J Appl Ecol, 2014, 25(05): 1499-1505. (in Chinese with English abstract)

[59] Yang T, Li YY. Study of urease activity in the soil of alpine meadow. Acta Ecol Sin, 1985, 8(3): 28-33. (in Chinese with English abstract)

[60] Yu WS, Li GH, Wang BS, Wu HL, Zhao YW, Meng FH, Lu CA. Component characteristics of soil labile and recalcitrant carbon under long-term different fertilization systems in eastern China. J Plant Nutr Fertil, 2015, 21(03): 675-683. (in Chinese with English abstract)

[61] Zhang B. Application of metagenomics for understanding of environmental microbiology. Heilongjiang Med J, 2014, 27(02): 267-271. (in Chinese with English abstract)

[62] Zhao Z, Wang D, Wang Y, Mu Z, Zhu J. Wet deposition flux and runoff output flux of mercury in a typical small agricultural watershed in three Gorges reservoir areas. Environ Sci Pollut Res, 2015, 22(7): 5538-5551.

[63] Žifčáková L, Větrovský T, Howe A, Baldrian P. Microbial activity in forest soil reflects the changes in ecosystem properties between summer and winter. Environ Microbiol, 2016, 18(1): 288-301.

[64] Ziuzina D, Boehm D, Patil S, Cullen PJ, Bourke P. Cold plasma inactivation of bacterial biofilms and reduction of quorum sensing regulated virulence factors. PLoS ONE, 2015 10 9 e0138209