Skip to main content

ORIGINAL RESEARCH article

Front. Microbiol., 23 August 2022
Sec. Terrestrial Microbiology

Soil metatranscriptome demonstrates a shift in C, N, and S metabolisms of a grassland ecosystem in response to elevated atmospheric CO2

  • 1Institute of Applied Microbiology, Justus Liebig University, Giessen, Germany
  • 2Faculty of Basic and Biomedical Sciences, Simón Bolívar University, Barranquilla, Colombia
  • 3Institute of Plant Ecology, Justus Liebig University, Giessen, Germany
  • 4School of Biology and Environmental Science and Earth Institute, University College Dublin, Dublin, Ireland

Soil organisms play an important role in the equilibrium and cycling of nutrients. Because elevated CO2 (eCO2) affects plant metabolism, including rhizodeposition, it directly impacts the soil microbiome and microbial processes. Therefore, eCO2 directly influences the cycling of different elements in terrestrial ecosystems. Hence, possible changes in the cycles of carbon (C), nitrogen (N), and sulfur (S) were analyzed, alongside the assessment of changes in the composition and structure of the soil microbiome through a functional metatranscriptomics approach (cDNA from mRNA) from soil samples taken at the Giessen free-air CO2 enrichment (Gi-FACE) experiment. Results showed changes in the expression of C cycle genes under eCO2 with an increase in the transcript abundance for carbohydrate and amino acid uptake, and degradation, alongside an increase in the transcript abundance for cellulose, chitin, and lignin degradation and prokaryotic carbon fixation. In addition, N cycle changes included a decrease in the transcript abundance of N2O reductase, involved in the last step of the denitrification process, which explains the increase of N2O emissions in the Gi-FACE. Also, a shift in nitrate (NO3-) metabolism occurred, with an increase in transcript abundance for the dissimilatory NO3- reduction to ammonium (NH4+) (DNRA) pathway. S metabolism showed increased transcripts for sulfate (SO42-) assimilation under eCO2 conditions. Furthermore, soil bacteriome, mycobiome, and virome significantly differed between ambient and elevated CO2 conditions. The results exhibited the effects of eCO2 on the transcript abundance of C, N, and S cycles, and the soil microbiome. This finding showed a direct connection between eCO2 and the increased greenhouse gas emission, as well as the importance of soil nutrient availability to maintain the balance of soil ecosystems.

Introduction

World atmospheric carbon dioxide (CO2) concentration has increased by about 50%, from pre-industrial levels of about 278 parts per million volume (ppmV) to the current concentration of more than 415 ppmV (IPCC, 2021; NASA, 2022), and the current anthropogenic emissions of the greenhouse gases (GHG) are the highest in history (IPCC, 2021). Because terrestrial ecosystems act as a “sink” for a significant portion of the global carbon (C), fluctuations in net C exchange between soil and atmosphere impact the CO2 concentration in the atmosphere profoundly (DOE.2020, 2020). Hence, the response of terrestrial ecosystems to increasingly higher concentrations of CO2 under a changing climate has important implications for the global carbon cycle (Vestergard et al., 2016). In this sense, it has been widely described that elevated CO2 (eCO2) concentration affects plants in such a way that it decreases transpiration (Owensby et al., 1997; Kimball, 2016) and increases growth (Idso, 1994; He et al., 1995), plant yield (Kimball, 1983), photosynthetic capacity (Habash et al., 1995; He et al., 1995; Johnson and Pregitzer, 2007), below-ground biomass (Jongen et al., 1995), and the efflux amounts of root exudates (Phillips et al., 2012; Jia et al., 2014; Dong et al., 2021).

Consequently, the supply of fresh plant-derived C to the soil matrix due to eCO2 may accelerate the decomposition of soil organic matter (SOM) and decrease soil C stocks (Fontaine et al., 2004; Blagodatskaya and Kuzyakov, 2008), a process known as “the priming effect”. This alteration in the increased decomposition of SOM has been previously reported in different ecosystems, such as grasslands (Vestergard et al., 2016; Liu et al., 2017), forests (Phillips et al., 2012; Qiao et al., 2014; Liu et al., 2017), and crop fields (Trivedi et al., 2016). Old SOM pools contain significant physically and chemically protected N stocks; consequently, soil microorganisms under plenty of C supply gain access to a reservoir of N to meet their enhanced N demand (Derrien et al., 2014; Vestergard et al., 2016; Liu et al., 2017), causing alterations in soil N balance and N cycle. This process has been described for the Gi-FACE grassland by (Müller et al. 2009), who reported that under eCO2, the mineralization of labile organic N became more important. (Müller et al. 2009) found that eCO2 caused an increase in dissimilatory NO3- reduction to NH4+ (DNRA) and immobilization of NO3- and NH4+ ions. Other alterations in the N cycle due to eCO2 have been described by Kammann et al. (2008), who indicated an increase in N2O emissions. Likewise, Moser et al. (2018) reported that N2O emissions were 1.79-fold higher for the Gi-FACE grassland under the eCO2 treatment. Also, Moser et al. (2018) described that N2O emissions from denitrification, nitrification, and heterotrophic nitrification showed a 2.09-fold, 1.64-fold, and 1.66-fold increase, respectively. More recently, Du et al. (2022), based on a meta-analysis, indicated that eCO2 significantly increased N2O emissions, NO3- content, and soil microbial biomass N by 44, 13, and 7% for agricultural soils.

Likewise, C and N cycle changes are directly related to the soil microbiome and soil microbial processes. For example, Xu et al. (2013) described that the abundance of genes involved in labile C degradation and C and N fixation, and denitrification processes significantly increased under eCO2. Similarly, He et al. (2014) and Xiong et al. (2015) have reported a shift in soil microbial communities under eCO2 in a soybean and maize agroecosystem, respectively. These changes included stimulation of key functional genes involved in carbon fixation and degradation, nitrogen fixation, denitrification, methane metabolism, and phosphorus cycling. Simonin et al. (2015) reported that shoot biomass, root biomass, and soil respiration were increased under eCO2 and N supply, and these variables were positively correlated with the abundance of ammonia-oxidizing bacteria. Le Roux et al. (2016) described that the potential nitrite oxidation rate was enhanced in soil by eCO2. Furthermore, the increase in soil microbial C and N cycling may be accompanied by microbial sulfur (S) and phosphorus (P) demand (Xiong et al., 2015; Yu et al., 2018a, 2021). Regarding S cycle alterations under eCO2, (Yu et al. 2018a,b); Yu et al. (2021) have reported an increase in S cycling in semiarid grassland soils exposed to eCO2, indicating a significant increase in the abundance of dsrA, dsrB, and sox genes. Likewise, Padhy et al. (2020) described that several genera, such as Desulfatibacillum, Desulfotomaculum, Desulfococcus, and Desulfitobacterium, were more abundant under eCO2 conditions in a lowland rice field and that several enzymes involved in S assimilation pathways showed higher counts at eCO2 concentrations as well.

Nonetheless, all the above-mentioned studies utilized a DNA-based approach to assessing the changes in the gene abundance involved in the C, N, and S cycles and the microbiome composition under eCO2 conditions. DNA-based approaches could lead to biases because DNA from dead cells or free DNA represented a significant fraction of microbial DNA in many soils (Carini et al., 2016). In addition, DNA from dead cells can remain in soils for weeks to years and may cloud DNA-based assessments of microbiome analyses (Dlott et al., 2015; Morrissey et al., 2015). Therefore, using RNA instead of DNA for metastudies provides an ideal tool to study the microbial populations that actively participate in various ecological processes (Sharma and Sharma, 2018). In this sense, some studies were done in the Giessen free-air CO2 enrichment experiment (Gi-FACE) in Giessen, Germany, which addressed this issue by performing microbiome metatranscriptomics analyses with rRNA and mRNA. Their findings were that eCO2 significantly affected the expression of 16S rRNA, transcription machinery, oxidative phosphorylation, translation and transcription factors, membrane transport, and nucleotide metabolism, among others, associated with rhizosphere microbiomes and plant roots. Likewise, the structure and composition of the rhizosphere soil microbiome were the most affected by eCO2 (Bei et al., 2019; Rosado-Porto et al., 2021). Furthermore, these reports showed that through the use of RNA instead of DNA, it was possible to assess the effects of eCO2 on the soil microbiome in the Gi-FACE, contrary to the previous studies, which reported little or no effect of it (Regan et al., 2011; de Menezes et al., 2016; Brenzinger et al., 2017).

Nevertheless, in the current literature, the use of mRNA metatranscriptomics to assess the effects of eCO2 conditions on C, N, and S cycle processes has not been described. mRNA metatranscriptomics allows addressing which genes are transcribed and to what extent, thereby enabling the demonstration of the functions of a potential range of microorganisms (Franzosa et al., 2014). From such functional data, active metabolic pathways can be identified in the microbiome and can be associated with particular environmental conditions, offering a more informative perspective, as it can reveal details about populations that are transcriptionally active (Bashiardes et al., 2016). Therefore, we hypothesized that the soil metatranscriptome was significantly affected by the eCO2 treatment, including the abundance of transcripts involved in nutrient cycling, because of the shifted microbial community under eCO2 (Bei et al., 2019; Rosado-Porto et al., 2021). The aims of the present study were: i) to assess the effect of long-term eCO2 concentrations on active soil bacteriome, mycobiome, protistome, and virome through an mRNA-based approach; ii) to evaluate the influence of eCO2 on C, N, and S cycle expressed genes in a grassland ecosystem; and iii) to propose an interaction model of C, N, and S cycle processes under eCO2 conditions.

Materials and methods

Study site description

The Gi-FACE study is located at 50°32'N and 8°41.3'E near Giessen, Germany, at an elevation of 172 m above sea level. It consists of three pairs of rings with a diameter of 8 m; each pair consists of an ambient and an elevated CO2 treatment ring (Jäger et al., 2003). From May 1998 until the present, eCO2 rings have been continuously enriched by 20% above ambient CO2 concentrations during daylight hours. Ambient and elevated CO2 rings are separated by at least 20 m, and each pair is placed at the vertices of an equilateral triangle. The presence of a slight slope within the experimental site (between 0.5 and 3.5°) places the rings on a moisture gradient, such that pair 1 has the lowest mean moisture content (38.8 ± 10.2%) and pair 2 has the highest mean moisture content (46.1 ± 13.2%), whereas pair 3 is intermediate (40.7 ± 11%) (Jäger et al., 2003; de Menezes et al., 2016). The average annual air temperature and precipitation are 9.4 °C and 580 mm, respectively.

The vegetation is an Arrhenatheretum elatioris Br.Bl. Filipendula ulmaria subcommunity, dominated by Arrhenatherum elatius, Galium album, and Geranium pratense. At least 12 grass species, 15 non-leguminous herbs, and up to 5 legumes with small biomass contributions (<5%) are present within a single plot (Andresen et al., 2018). The experimental field has not been plowed for more than 100 years. It has received N fertilization in the form of granular mineral calcium-ammonium-nitrate (40 kg N ha−1 year−1) once a year since 1995 and has been mown two times a year since 1993. The soil at the Gi-FACE site is classified as Fluvic Gleysol; its texture is a sandy clay loam over a clay layer, with pH = 6.2 and average C and N contents of 4.5 and 0.45%, respectively, as measured in 2001 (Jäger et al., 2003).

Soil sampling, total RNA extraction, and ribodepletion

Soil sampling was performed utilizing sawed-off 50 ml syringes (11 × 3 cm), and four samples were taken to a depth of ~10 cm within each ring in September 2017. Once taken, samples were refrigerated and transported to the lab for immediate processing. Upon arrival in the laboratory (10-min driving), soil cores were gently shaken by hand to remove loosely attached soil (bulk soil), while the soil that remained attached to the roots was considered rhizosphere soil. Rhizosphere soil was detached from the roots with sterile tweezers and directly sieved. If roots were still present, they were sorted out before sieving. Bulk and rhizosphere soils were sieved (<2 mm) and stored at −80 °C for further analyses. Sample processing took less than 1 h before freezing.

Total RNA extraction was performed following a modified protocol of Mettel et al. (2010), as described by Rosado-Porto et al. (2021). After extraction, samples were treated for DNA digestion with RNase-Free DNase Set (QIAGEN GmbH - Germany) according to the manufacturer's instructions. DNase reaction was stopped with 10 μl of 50 mM EDTA. With the DNA-free RNA, a PCR was carried out, using the universal 16S rRNA gene primers 27F (5'-AGAGTTTGATCMTGGATCMTGGCTCAG-3') and 1492R (5'- GGTTACCTTGTTACGACTT-3') (Lane, 1991; Weisburg et al., 1991) and checked on agarose gel electrophoresis to verify the absence of remaining DNA in the samples as described by Rosado-Porto et al. (2021). Afterward, total RNA technical replicates were pooled into a composite sample according to the ring number and rhizosphere or bulk soil. Later, total RNA samples were ribodepleted using the MICROBExpress™ Kit (Life Technologies, 5791, Carlsbad – California, USA), following the manufacturer's instructions. Finally, mRNA was reverse transcribed to produce double-stranded cDNA (LGC Genomics GmbH, Berlin, Germany).

cDNA sequencing and metatranscriptomics analysis

The cDNA products were sequenced with Illumina MiSeq V3 (2 × 300 bp), 40 M read pairs/12 Gb of raw data (LGC Genomics GmbH, Berlin, Germany). After sequencing, all libraries for each sequencing lane were demultiplexed using the Illumina bcl2fastq 2.17.1.14 software (Illumina, 2019). Later, sequencing adapter remnants were removed, and reads with a final length of <100 bases were discarded. Afterward, the sequencing outputs were analyzed using SqueezeMeta version 1.3.1 (Tamames and Puente-Sánchez, 2019). Next, sequence assembly was performed using Megahit (Li et al., 2015), and the removal of short contigs (<200 bps) was done with Prinseq (Schmieder and Edwards, 2011). Afterward, RNAs, tRNA/tmRNA, and open reading frames (ORFs) were predicted using Barrnap (Seemann, 2014), Aragorn (Laslett and Canback, 2004), and Prodigal (Hyatt et al., 2010), respectively. Subsequently, Diamond was utilized (Buchfink et al., 2015) to perform the alignment and search of similarities against GenBank (Clark et al., 2016), eggNOG (Huerta-Cepas et al., 2016), and Kyoto Encyclopedia of Genes and Genomes (KEGG) (Kanehisa and Goto, 2000) databases, using an e-value of 10 × 10−3 and minimum identity values of 40 and 30 for taxonomical and functional identities, respectively. Additionally, HMM homology searches were done by HMMER3 (Eddy, 2009) for the Pfam database (Finn et al., 2016) applying an e-value of 10 × 10−10. Moreover, additional ORFs were produced by Diamond BlastX (Buchfink et al., 2015), implementing an e-value of 10 × 10−3 and a minimum identity value of 40. Later, taxonomical classification of the mRNA transcripts was performed using the hits for each query gene utilizing the results of the Diamond search against the GenBank nr database and applying the lowest common ancestor (LCA) algorithm, allowing a 0.2% of different taxa from the LCA and the minimum number of hits per taxa of 2. The read mapping against contigs was performed using Bowtie2 (Langmead and Salzberg, 2012), and the binning was done utilizing MaxBin2 (Wu et al., 2016) and Metabat2 (Kang et al., 2019), and later bin statistics were computed using CheckM (Parks et al., 2015).

Diversity and differential abundance analyses

For the analysis of SqueezeMeta output, data were imported into R studio software 1.1.419 with package SQMtools version 0.6.1. (Puente-Sánchez et al., 2020). For diversity assessment of bacteria, archaea, fungi, viruses, and protists, frequency tables were created and analyzed with package Phyloseq 1.28.0 (McMurdie and Holmes, 2013). Core features for each of the above taxonomical groups were calculated for eCO2 and aCO2 conditions by transforming the frequency table counts to relative abundance with Microbiome package version 1.8.0 (Lahti and Shetty, 2019). Later, features with a total relative abundance ≥10 × 10−4% and present in ≥95% of samples were included as part of the core. Likewise, for KEGG and GenBank Clusters of Orthologous Groups (COG) protein outputs, features with unknown functions or unassigned names were removed from the frequency tables, and core features were calculated as described above.

For beta diversity analysis, core datasets were transformed using the centered log-ratio (clr) method (Aitchison, 1982, 1986) using the R package ALDEx2 1.22.0 (Fernandes et al., 2013, 2014). Afterward, community dissimilarity distance matrices were generated using the Aitchison distance (Aitchison, 1982, 1986) and visualized using principal components analysis (PCA) (Jolliffe and Cadima, 2016). Statistical differences between CO2 conditions were assessed by a permutational multivariate analysis of variance using the Adonis method and employing 999 permutations (Anderson, 2001).

Differential abundance analysis of core features was done with R package ALDEx2 1.22.0 (Fernandes et al., 2013, 2014) by performing the clr transformation using as denominator the geometric mean abundance of all features and 128 Monte Carlo instances. Subsequently, features with absolute ALDEx2 effect sizes of >0.8, >0.5, and >0.2 were considered to have a significantly greater, moderate, and slightly higher abundance, respectively, between aCO2 and eCO2 rings (Sawilowsky, 2009).

Pathway reconstruction analysis

Pathway prediction for KEGG (Kanehisa and Goto, 2000) and MetaCyc (Caspi et al., 2018) databases was done using MinPath (Ye and Doak, 2009). In addition, pathway reconstruction and assessment of the log2 fold change between aCO2 and eCO2 rings were performed with SQMtools version 0.6.1. (Puente-Sánchez et al., 2020) and its function “exportPathway” and analyzing feature frequencies as relative abundances.

Results

Sequencing and assembly

In total, 23,970,892,090 bases were obtained, comprising 90,534,066 raw sequences, from which 72,253,754 sequences were mapped and assembled with Megahit, with the percentage of sequences successfully mapped per sample ranging between 81.14 and 78.12%. A total of 1,396,973,823 bases from the mapped sequences were retained after short contigs were removed and assembled into 3,997,902 contigs with lengths ranging from 9,714 to 200 bases. From the contigs, there were predicted 4,063,836 ORFs, 1,199,550 rRNAs, and 2,406 tRNAs/tmRNAs, which subsequently were annotated, producing 92,698 successfully annotated taxa and 483,556 and 1,163,975 KEGG and COG annotations, respectively. Details regarding the number of total annotated transcripts per sample are provided in Supplementary material S3.

Beta diversity and microbe differential abundance

Metatranscriptome results from the Gi-FACE exhibited changes in the composition and structure of soil microbial communities due to elevated concentrations of atmospheric CO2. Our data indicated that the soil core bacteriome (p < 0.05), mycobiome (p < 0.05), and virome (p < 0.05) were the most affected by eCO2 concentrations and showed significantly different compositions between aCO2 and eCO2 rings, according to the permutational multivariate analysis (Figures 1A–C). In contrast, the general structure of the Gi-FACE soil core archaeome and protistome was not significantly affected by eCO2 (Figures 1D,E).

FIGURE 1
www.frontiersin.org

Figure 1. Beta diversity analysis of core microbial taxa and transcripts from the Gi-FACE metatranscriptome. Principal components analysis (PCA) of (A) bacteria, (B) fungi, (C) virus, (D) archaea, (E) protist, (F) other eukarya, (G) GenBank COG, and (H) KEGG functions; p-values from PERMANOVA test.

Moreover, differential abundance results from ALDEx2 (Figure 2) indicated that several taxa were significantly increased or decreased under eCO2 conditions and that these affected taxa shaped the soil microbiome of the Gi-FACE. Besides, differential abundance results showed that the number of bacterial taxa that were positively stimulated under eCO2 is greater than the number of taxa that were negatively affected. Among the bacterial taxa which were highly stimulated under eCO2 conditions are Flavobacterium, Ruminiclostridium, Gemmata, Dehalococcoides, Minicystis, Ureaplasma, Saccharopolyspora, Asaia, Nocardioides, Defluviimonas, Bacillus, Nannocystis, Glaesserella, Pedosphaera, Arenimonas, Nitrospirae bacterium, Blastopirellula, Amycolatopsis, Tatlockia, Povalibacter, Thermasporomyces, Halolactibacillus, Clostridium, Pedobacter, Aminipila, Rhodovastum, Pirellula, and Burkholderia, which showed ALDEx2 effect sizes between 1.5 and 0.8 (Figure 2A, Supplementary material S1.1). Likewise, soil mycobiome was shaped by several fungi greatly affected under eCO2 conditions, most belonging to phyla Basidiomycota, Mucoromycota, and Ascomycota, as is the case of the genus Aspergillus (phylum Ascomycota), which showed an ALDEx2 effect size of 1.15 (Figure 2B, Supplementary material S1.2). Additionally, fungi like Rhizopus, Cadophora, Gigaspora, Histoplasma, and Aplosporella were also highly stimulated in eCO2 rings with effect sizes from 0.86 to 1.33 (Figure 2B). Regarding the Gi-FACE soil virome, viruses like Brome mosaic virus, Panicovirus, and Cocksfoot mild mosaic virus decreased in eCO2 rings presenting effect sizes ranging from −0.58 to −0.74. In contrast, viruses such as Penicillium discovirus, unclassified Picornavirales, and unclassified Endornaviridae were positively affected under eCO2 conditions with effect sizes from 0.54 to 0.87 (Figure 2D, Supplementary material S1). Moreover, some viral features belonging to the families Leviridae, Siphoviridae, Bromoviridae, and Dicistroviridae were affected by eCO2 as well (Supplementary material S1.3).

FIGURE 2
www.frontiersin.org

Figure 2. Differential abundance of core microbial taxa from the Gi-FACE metatranscriptome of (A) bacteria, (B) fungi, (C) archaea, (D) virus, and (E) protist. ALDEx2 results of features with an ALDEx2 effect size > 0.5 using centered log ratio (clr) transformation and the geometric mean abundance of all features.

Although our data did not show that eCO2 significantly influenced the general structure of the soil archaeome and protistome, the differential abundance test showed that some archaea and protist taxa were either positively or negatively affected under eCO2 conditions (Figures 1D,E, 2C,E, Supplementary materials S1.4, S1.5).

Functional metatranscriptome and differential abundance

Beta diversity of transcripts analyzed against GenBank, COG, and KEGG databases showed that the functional metatranscriptome was greatly affected under eCO2 conditions in which the annotations performed to both databases were significantly different in their structure and composition between eCO2 and aCO2 conditions (Figures 1G,H). After removing unclassified and non-characterized proteins, 7,780 remained for GenBank COG and 8,880 for KEGG datasets. Furthermore, our data indicated that the sequences analyzed against both databases showed similar results regarding the number of proteins with an ALDEx2 effect size greater than 0.5. In the case of GenBank COG data, 146 transcripts were moderately or greatly stimulated under eCO2 conditions, in contrast to 161 negatively affected under these conditions. Similarly, KEGG results showed that the abundance of 147 and 156 transcripts was positively and negatively affected, respectively (Supplementary material S2.1).

Moreover, eCO2 conditions positively influenced the transcript abundance of several COG categories, such as energy production and conversion, inorganic ion transport and metabolism, cell envelope biogenesis, outer membrane intracellular trafficking carbohydrate transport and metabolism, and signal transduction mechanisms Figure 3, (Supplementary material S2.2). In contrast, categories for translation, ribosomal structure, and biogenesis; transcription; secondary metabolite biosynthesis, transport, and catabolism; nucleotide transport and metabolism; DNA replication, recombination, and repair; and coenzyme metabolism were negatively affected at eCO2 concentrations (Figure 3, Supplementary material S2.2).

FIGURE 3
www.frontiersin.org

Figure 3. Differential abundance of Gi-FACE metatranscriptome core transcripts annotated against GenBank Clusters of Orthologous Groups (COG) and grouped according to COG categories. Results expressed as relative abundance (right) and ALDEx2 effect size (left) of transcripts with an ALDEx2 effect size > 0.5 using centered log ratio (clr) transformation and the geometric mean abundance of all features.

Nitrogen cycle

The data showed that, under eCO2 conditions, an increase in the dissimilatory nitrate (NO3-) reduction to ammonium (NH4+) (DNRA) pathway and a decrease in the assimilatory NO3- reduction to NH4+ occurred (Figure 4A). Specifically, the transcripts for the DNRA enzymes nitrite reductase (NADH) (NirBD) and nitrate reductase (NarGHI) showed greater abundances under eCO2 conditions (Figures 4A,B). In contrast, the transcript abundance for the enzymes nitrate reductase (NAD(P)H) (NR), ferredoxin-nitrite reductase (NirA), and assimilatory nitrate reductase (NasAB) was negatively affected at eCO2 concentrations (Figures 4A,B, Supplementary material S2.3).

FIGURE 4
www.frontiersin.org

Figure 4. Reconstructed KEGG (Kanehisa and Goto, 2000) N pathways of NO3- assimilatory and dissimilatory reduction, denitrification, and nitrification processes. (A) N transformations expressed as log2 fold change of transcripts as relative abundance. ALDEx2 effect size: (**) >0.5, (*) >0.2. (B) Differential abundance of N cycle transcripts with ALDEx2 effect sizes > 0.2.

Similarly, the denitrification process showed changes as well. Transcript abundance for the denitrification enzymes nitrate reductase/nitrite oxidoreductase (NarGHI/NapAB) and nitric oxide reductase (NorBC) showed higher levels under eCO2 conditions with ALDEx2 effect sizes of 0.44 and 0.64. These enzymes are responsible for the transformation of NO3- to nitrite (NO2-) and the reduction of nitric oxide (NO) to nitrous oxide (N2O), respectively. On the contrary, the transcription of the gene for the enzyme nitrous oxide reductase (NosZ), which catalyzes the transformation of N2O to atmospheric nitrogen (N2), was reduced in the eCO2 rings (Figures 4A,B, Supplementary material S2.3).

Likewise, the nitrification process was also affected. The data exhibited changes in the expression patterns of the enzymes methane/ammonia monooxygenase (AmoCAB) and nitrate reductase/nitrite oxidoreductase (NrxAB), which were negatively and positively affected, respectively (Figures 4A,B, Supplementary material S2.3). Furthermore, pathway reconstruction and differential abundance analyses did not show significant changes in the abundance of N fixation enzymes under eCO2 conditions.

Sulfur cycle

The metatranscriptomics results indicated changes in the dissimilatory and assimilatory pathways of sulfate (SO42-) reduction. Transcript abundance for the enzymes sulfate adenylyltransferase (Sat) and adenylylsulfate reductase (AprAB), part of the dissimilatory SO42- reduction pathway, was highly decreased under eCO2 conditions. Furthermore, the transcript abundance for the enzyme AprAB was the one that showed the highest decrease, with an ALDEx2 effect size of −0.86. This protein catalyzes the transformation of sulfite (SO32-) to adenosine 5'-phosphosulfate (APS) (Figures 5A,B).

FIGURE 5
www.frontiersin.org

Figure 5. Reconstructed KEGG (Kanehisa and Goto, 2000) pathways of S metabolism. (A) Assimilatory SO42- reduction and dissimilatory SO42- reduction and oxidation processes expressed as log2 fold change of transcripts as relative abundance. ALDEx2 effect size: (***) >0.8, (**) >0.5, (*) >0.2. (B) Differential abundance of S cycle transcripts with absolute ALDEx2 effect sizes > 0.1 involved in assimilatory SO42- reduction, dissimilatory SO42- reduction, and oxidation, and uptake of S from organic compounds and sulfide synthesis.

Similarly, the decrease in transcript abundance of the sulfate adenylyltransferase subunit 2 (CysND) and sulfate adenylyltransferase (PAPSS), part of the assimilatory sulfate reduction pathway, suggested a decrease in the reduction of SO42- to APS at eCO2 concentrations. These two proteins showed ALDEx2 effect sizes of −0.58 and −0.41, respectively (Figures 5A,B, Supplementary material S2.4). Nonetheless, the enzymes adenylylsulfate kinase (CysC) and sulfite reductase (NADPH) (CysJI), which are involved in the reduction of APS to 3′-Phosphoadenosine-5′-Phosphosulfate (PAPS) and the reduction of SO32- to sulfide (S2−), respectively, were increased at eCO2 concentrations (Figures 5A,B, Supplementary material S2.4).

Moreover, our data showed that several enzymes belonging to pathways responsible for the transformation of organic S compounds had higher transcript abundances at eCO2, as is the case with dimethylsulfone monooxygenase, thiosulfate dehydrogenase [quinone], and taurine dioxygenase (Figure 5B, Supplementary material S2.4). Although the SOX system for the oxidation of S was generally not over-expressed under eCO2 concentrations, transcripts of the enzyme sulfane dehydrogenase subunit (SoxC) showed a slight increase at these conditions, with an ALDEx2 effect size of 0.28.

Carbon cycle

Functional metatranscriptome showed changes in the metabolism of C compounds. The main changes comprised a general increase in transcripts from the glycolytic and pentose phosphate pathways, which included the increase in the abundance of transcripts for the enzymes phosphoglucomutase, glucose-6-phosphate isomerase, phosphoenolpyruvate carboxykinase (ATP), pyruvate water dikinase, 2-oxoglutarate, gluconate 2-dehydrogenase, gluconolactonase, transketolase, and xylulose-5-phosphate/fructose-6-phosphate phosphoketolase, all with ALDEx2 effect sizes ranging from 0.79 to 0.52 (Figure 6A). Likewise, the data exhibited an increase in the transcription of genes coding for enzymes responsible for the degradation of chitin, cellulose, and aromatic compounds, for example, alpha-N-arabinofuranosidase; endo-1,4-beta-xylanase, and chitinase (Figure 6A). In contrast, the metabolism of fatty acids, starch, and sucrose was negatively affected under eCO2 conditions, with the most affected features having ALDEx2 effect sizes ranging from −0.84 to −0.51 (Figure 6A, Supplementary material S2.5).

FIGURE 6
www.frontiersin.org

Figure 6. Differential abundance of transcripts grouped by KEGG Orthology (KO) second level of (A) carbon compounds metabolism and (B) ABC transporters.

Furthermore, our results indicated a stimulation in the metabolism of aromatic, branched-chain, and sulfur amino acids. In the case of sulfur amino acid metabolism, an increase in the transcript abundance of enzymes for the metabolism of homocysteine, taurine, and thiol groups occurred (Supplementary material S2.5). Likewise, the transcript abundance of several enzymes involved in the degradation of aromatic amino acids was highly stimulated, for example, aminocarboxymuconate-semialdehyde decarboxylase case enoyl-CoA hydratase, amidase, monoamine oxidase, acylpyruvate hydrolase, and gentisate 1,2-dioxygenase, with ALDEx2 effect sizes between 1.10 and 0.54. Moreover, the transcript abundance for genes involved in the Arnon–Buchanan cycle (reductive citric acid cycle) increased under eCO2, including key enzymes like phosphoenolpyruvate carboxykinase (ATP), pyruvate water dikinase, and pyruvate ferredoxin oxidoreductase (Figure 6A, Supplementary material S2.5).

ABC membrane transporters

The metatranscriptomic data on the ABC membrane transport proteins suggested changes in the uptake and transport of different carbon compounds under eCO2 conditions. For example, membrane transporters for glucose/mannose, α-glucoside, ribose/D-xylose, and chitobiose increased at eCO2 concentrations. In contrast, there was a decrease in the expression of membrane transporters for raffinose/stachyose/melibiose, rhamnose, galactose oligomer/maltooligosaccharide, maltose, and fructose (Figure 6B, Supplementary material S2.6). Similarly, a shift in the ABC transporters for amino acids occurred, with an increase in the transcript abundance of the transporters for general L-amino acids and branched-chain amino acids and a decrease in glutamate/aspartate and oligopeptide transporters (Figure 6B, Supplementary material S2.6). Additionally, other membrane transport proteins over-expressed under eCO2 conditions, including transporters for osmoprotectants, lipopolysaccharides, and iron (II) (Figure 6B, Supplementary material S2.6).

Discussion

Soil microbiome response to eCO2

Our results on the functional metatranscriptome of the Gi-FACE confirm previous reports from Bei et al. (2019) and Rosado-Porto et al. (2021) on the changes in microbiome composition due to eCO2 concentrations. Additionally, the outcome expands our understanding of eCO2 concentrations in N, S, and C cycles. The metatranscriptome data presented in this study are a one-time snapshot of the gene expression of the microbiome. Although transcriptomes are dynamic and variable in response to diverse environmental factors (Nuccio et al., 2021; Zhao et al., 2021), it nevertheless provides valuable information on soil ecosystems' response to climate change scenarios.

Our data confirm that the structure of the Gi-FACE soil bacteriome was strongly affected under eCO2 (Figure 1A). Prior studies have already portrayed differences in the composition of the bacteriomes between aCO2 and eCO2 conditions (Bei et al., 2019; Rosado-Porto et al., 2021), who have described significant changes mainly in the rhizosphere. Additionally, several bacterial taxa found in the present study (Figure 2A, Supplementary material S1.1) have already been described as stimulated under eCO2 conditions, as is the case of genera Bacillus, Burkholderia, Mesorhizobium, Streptomyces, and Dongia (Rosado-Porto et al., 2021). Besides the soil bacteriome, the results showed that the soil mycobiome was greatly affected at eCO2 concentrations. Like Bei et al. (2019), our data indicated that the Gi-FACE mycobiome was composed mainly of phyla Basidiomycota, Mucoromycota, and Ascomycota (Supplementary material S1.2). Moreover, several highly affected fungi belonged to Ascomycota families, such as Mycosphaerellaceae, Didymosphaeriaceae, Ophiocordycipitaceae, Saccharomycetaceae, and Aspergillaceae, and Mucoromycota families, such as Cunninghamellaceae, Rhizopodaceae, and Glomeraceae (Supplementary material S1.2). Nevertheless, although our results showed a significant effect of eCO2 concentrations on the mycobiome composition, the reports of its effect on soil fungal communities vary according to different authors. Carney et al. (2007) described a decrease in fungal abundance under eCO2 conditions, whereas some others reported no significant change in the fungal communities (He et al., 2010; Hayden et al., 2012). This indicates that the response of fungal communities to eCO2 depends on other environmental factors like temperature and soil moisture and may be ecosystem specific as well.

The effects of CO2 concentration on the metagenome or metatranscriptome of soil archaeomes have not been widely studied; however, some reports described a strong influence of CO2 concentrations on soil archaeal communities (Hayden et al., 2012; Lee et al., 2015; Lee and Kang, 2016). Although the Gi-FACE archaeome did not show significant differences in its structure and composition in response to eCO2 concentrations, some taxa showed changes in their abundance (Figure 2C, Supplementary material S1), most belonging to the family Nitrosopumilaceae (phylum Thaumarchaeota). In addition, in the present study, the core archaeome was mainly composed of the phylum Euryarchaeota, contrary to that reported by Bei et al. (2019), who reported the phylum Thaumarchaeota as the most abundant.

Furthermore, our data exhibited that alongside the soil bacteriome and mycobiome, the soil core virome was affected by the eCO2 in the Gi-FACE (Figure 2D, Supplementary material S1.3). So far, in the current literature, there are no reports about the effects of eCO2 on the soil virome. Moreover, some reports indicated that, in general, the diversity of the soil virome is highly underestimated, with most of the current information focused on bacterial phages, while almost nothing is known about viruses that infect non-bacterial soil microbes, such as the archaea, fungi, and soil protozoa (Williamson et al., 2017; Pratama and van Elsas, 2018). Our results on the differential abundance of the core virome under eCO2 conditions suggested that several viral transcripts were reacting to changes in bacterial, archaeal, and fungal taxa. For example, the Leviridae and Siphoviridae families include viruses that use bacteria and archaea as hosts and in our study were increased under eCO2 conditions (Supplementary material S1.3) (Duin and Olsthoorn, 2012; Hendrix et al., 2012; Krupovic et al., 2020). Similarly, some fungal viruses, such as Mitovirus and Penicillium discovirus, have shown significant changes in their abundance under eCO2 conditions, which might be linked to the changes in some fungal features as observed in Penicillium oxalicum and members belonging to the class Ophiostomatales (Hong et al., 1999; Krishnamurthy, 2017).

Changes in C compound assimilation and priming effect

The data showed changes in transcript abundance from pathways involved in the metabolism of different C compounds, indicating that C dynamics have changed due to eCO2. It has been widely described that eCO2 increases the efflux of soluble sugars, amino acids, phenolic acids, and organic acids in the root exudates (Phillips et al., 2012; Jia et al., 2014; Dong et al., 2021), which produces the so-called “priming effect,” thus leading to an acceleration in SOM decomposition (Fontaine et al., 2004; Blagodatskaya and Kuzyakov, 2008). The metatranscriptome data showed an increase in the priming effect due to eCO2 concentrations in the Gi-FACE soil. The priming effect in our data is represented mainly by the over-expression of transcripts from the glycolysis and pentose phosphate pathways and an increase in transcript abundance for certain amino acid metabolism, alongside an increase in transcript abundance of enzymes responsible for the degradation of chitin, cellulose, and lignin (Figure 7). Similarly, He et al. (2010, 2014), Xiong et al. (2015), and (Yu et al. 2018a,b) have reported the stimulation of functional gene abundance for C compound degradation, either labile or recalcitrant under eCO2. Likewise, other authors have described the increase in the degradation of soil organic polymers as part of the decomposition of older soil C (Van Groenigen et al., 2005; Xie et al., 2005; Niklaus and Falloon, 2006; Vestergard et al., 2016). This enhancement of carbohydrate, amino acid, and SOM degradation would be reflected in a higher respiration rate and, consequently, a higher efflux of CO2. In our previous study, CO2 soil fluxes were 35% higher in eCO2 rings compared to the aCO2 ones (Rosado-Porto et al., 2021).

FIGURE 7
www.frontiersin.org

Figure 7. Model diagram of the interaction of C, N, and S cycles in the Gi-FACE.

Furthermore, the data suggested a shift in the uptake and use of C sources at eCO2 concentrations, reflected in a shift toward higher utilization of sugars and amino acids and a decrease in the metabolism of lipids, especially fatty acids (Figures 3, 6A). Additionally, the analysis of ABC membrane transporters revealed changes in the transcript abundance for saccharide uptake systems that are more often used under eCO2 conditions, indicating a shift in preference for the uptake of glucose, mannose, α-glucosides, ribose, xylose, and chitobiose instead of raffinose, stachyose, melibiose, rhamnose, galactose, maltose, and fructose. Moreover, our results revealed an increase in transcripts for prokaryotic carbon fixation at eCO2 concentrations. This included increases in enzymes, such as phosphoenolpyruvate carboxykinase, pyruvate water dikinase, and pyruvate ferredoxin oxidoreductase, accompanied by a decrease in the ribulose-bisphosphate carboxylase (Rubisco) (Figure 6A, Supplementary material S2.5). These results are opposite to the ones reported by He et al. (2010, 2014), Xu et al. (2013), Xiong et al. (2015), and (Yu et al. 2018a,b), who described a significant increase in the gene abundance of the Rubisco enzyme under eCO2 conditions. The aforementioned results suggest that in the Gi-FACE, the C fixation performed by prokaryotes at eCO2 concentrations is very likely to be done through the reverse reductive citric acid cycle, also known as the Arnon–Buchanan cycle (Evans et al., 1966; Buchanan and Arnon, 1990; Buchanan et al., 2017) instead of the Rubisco pathway. The increase in transcript abundance of the described enzymes involved in the reverse citric acid cycle could be associated with the increase of some members of the Nitrospirae and Aquificae phyla, which were significantly augmented under eCO2 (Figure 2, Supplementary material S1.1) (Berg, 2011; Alfreider et al., 2018; Mundinger et al., 2019; Steffens et al., 2021). Nonetheless, it is difficult to differentiate CO2 fixation for autotrophic growth from anapleurotic reactions of the citric acid cycle. Moreover, this autotrophic pathway might be more widespread among anaerobic and microaerobic bacteria (Berg, 2011) and possibly in aerobic bacteria (Buchanan et al., 2017).

Shift in N cycle processes

Changes in the N cycle have been previously described in the Gi-FACE (Kammann et al., 2008; Müller et al., 2009; Moser et al., 2018). However, the underlying microbiological mechanisms driving these processes were not detected until now. The metatranscriptomic results confirmed a switch in the NO3- reduction at eCO2 concentrations, from an assimilatory process to a DNRA, reflected by the increase in transcript abundance of the enzymes nitrite reductase (NADH) (NirBD) and nitrate reductase (NarGHI), responsible for the transformations of NO3- to NO2- and from NO2- to NH4+ in the DNRA process (Figures 4A,B, 7). Previously, (Müller et al. 2009), utilizing a 15N labeling approach, identified an increase in the DNRA and the immobilization of NH4+ and NO3-.

Additionally, our functional metatranscriptomic approach gives some clarity about the processes leading to the excess production of N2O under eCO2 conditions previously described by Kammann et al. (2008), Moser et al. (2018), and more recently in the meta-analysis of Du et al. (2022). The data suggest that the alteration in the denitrification process leads to an increase in N2O production. This increase results from over-expression of the N2O-producing enzyme nitric oxide reductase (NorBC) and a reduced expression of the enzyme nitrous oxide reductase (NosZ), which lead to a decrease in N2O reduction to N2. Furthermore, results indicate that the increase in N2O is caused both by an increase in N2O production and a decrease in N2O reduction (Figure 7). These results seem to denote that changes in nitrous oxide reductase (NosZ) occur at a transcriptional level. Contrary to our data, previous reports from Liu et al. (2010) and Bakken et al. (2012) indicated that the high N2O:N2 product ratio is a post-transcriptional phenomenon only due to the sensitivity of this enzyme to lower pH values that are usually found at eCO2 concentrations and therefore affect its translation/assembly.

The results also showed changes in the nitrification process, represented by a reduction in transcripts involved in the conversion of NH4+ to hydroxylamine (H3NO) performed by the enzyme methane/ammonia monooxygenase (AmoCAB), accompanied by an increase in the transcript abundance for the enzyme nitrate reductase/nitrite oxidoreductase (NrxAB), suggesting an increase in the rate of transformation from NO2- to NO3- (Figures 4A,B, 7). Additionally, a reduction in the transcript abundance for the first nitrification step was observed, which denotes that under eCO2 conditions, soil organisms obtain N from other sources instead of NH4+. Previously, (Müller et al. 2009) described that the mineralization of labile organic N became more critical at eCO2 concentrations.

The changes in transcript abundance for the nitrification process might suggest an increase in heterotrophic nitrification by fungi. The fungal nitrification comprises the oxidation of different forms of organic N, such as L-asparagine, propionamide, malonylmonohydroxamate, and 3-nitropropionate, using peroxidase enzymes (Hora and Iyengar, 1960; Marshall and Alexander, 1962; Doxtader and Alexander, 1966). More recently, Laughlin et al. (2008) and Zhu et al. (2015) have described that fungi carried out a significant part of the nitrification in soils and that they can simultaneously oxidize NH4+ and organic N. Moreover, many of the fungal taxa that are able to perform nitrification are members of the genus Aspergillus (Hora and Iyengar, 1960; Marshall and Alexander, 1962; Doxtader and Alexander, 1966), one of the most positively affected in the Gi-FACE mycobiome (Figure 2B). Therefore, the above-mentioned results might support the idea of the mining of SOM by soil microorganisms, very likely fungi, in order to fulfill their N requirements under eCO2, and could be linked to the previously reported data of higher C:N ratios in eCO2 rings (Brenzinger et al., 2017; Rosado-Porto et al., 2021).

Additionally, our data did not show any increase in the transcript abundance of N fixation enzymes under eCO2 conditions, in contrast to the reports from He et al. (2010, 2014); Xu et al. (2013); Xiong et al. (2015), and (Yu et al. 2018a,b). These results support the idea that in the Gi-FACE, the enhanced N requirements are being met through the uptake of organic sources. According to our results, these sources of organic N might have been the aromatic, sulfur, and branched-chain amino acids, since their metabolism and uptake were augmented at eCO2 concentrations (Figures 6A, 7).

S metabolism at eCO2 concentration

Most studies about the effects of eCO2 on the cycling of nutrients have focused on C and N cycles; nonetheless, the effects of eCO2 conditions have also been assessed for other elements, including S (He et al., 2010, 2014; Yu et al., 2018; Padhy et al., 2020). There are no reports about the changes in the S cycling and metabolism in the Gi-FACE. The results in the present study exhibited alterations in the metabolism of SO42-. There was a reduction in transcripts involved in the dissimilatory metabolism of SO42- reduction, evidenced by the decrease in the expression of the enzymes sulfate adenylyltransferase (Sat) and adenylylsulfate reductase (AprAB) under eCO2 conditions (Figures 5, 7). Similarly, the first step of the assimilatory SO42- reduction metabolism changed due to eCO2, comprised by the depletion in the reduction step from SO42- to APS. However, the other steps of the assimilatory SO42- reduction, from the reduction of APS up to the production of S2−, were increased under eCO2 conditions (Figure 5). This phenomenon could indicate that similar to N metabolism, due to the augmented C supply, S has also become a limiting element for the development of soil organisms. Thus, the assimilatory metabolism of S was enhanced at eCO2 concentrations as a response to this environmental pressure. Although there are not many reports about the effect of eCO2 on the S cycle, Yu et al. (2018b) have described that under eCO2, an increase in S cycling occurred. Likewise, Padhy et al. (2020) reported an increase in the genes of the assimilatory metabolism of S under eCO2 conditions.

Moreover, our data suggest that the obtention of S in the Gi-FACE is not from inorganic sources but from organic ones. This process is likely a consequence of the priming effect and the mining of S from the SOM. According to our data, one of the sources for the supply of S might be sulfur amino acids, and molecules with thiol groups, due to the metabolism of these compounds, were augmented under eCO2 conditions (Figures 5B, 7). Moreover, although our data did not show an overall increase in the transcripts for the SOX system for the acquisition of sulfur, a slight increase in transcript abundance of the enzyme sulfane dehydrogenase (SoxC) occurred. These data indicate that soil organisms have used organic molecules to supply the S requirements in the Gi-FACE under eCO2 concentrations. He et al. (2014) have reported similar results, describing an increase in the expression of sox genes under eCO2 conditions.

Conclusion

Our research showed for the first time how eCO2 simultaneously affects the gene expression in C, N, and S cycles, and potentially affects these processes at the ecosystem scale. The lower abundance of nitrogen fixation transcripts suggests that soil microorganisms degrade the SOM to fulfill their N requirements due to a higher C supply. Likewise, an increase in transcripts for carbohydrates, amino acids, chitin, lignin, and cellulose assimilation and degradation was observed. In addition, the changes in the transcript abundance of the DNRA and denitrification shed some light on the underlying mechanisms that lead to the increase of N2O emissions previously reported in the Gi-FACE. Additionally, this research presents evidence for the first time that the gene expression in the S cycle from the Gi-FACE changes at eCO2, comprised mainly by an increase in the transcript abundance of the assimilatory SO42- metabolism. Regarding soil microbiome structure, our findings confirmed previous data on the changes of bacterial and fungal communities at eCO2 concentrations. In addition, the results revealed new evidence of eCO2 effects on the virome with still unknown effects on ecosystem processes. In summary, our findings enhance our understanding of the observed changes in greenhouse gas fluxes under eCO2 that result in positive feedback by increased N2O, CH4, and CO2 emissions and reduced CH4 uptake.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/genbank/,PRJNA810964.

Author contributions

DR-P was responsible for the experiments, data curation, data analysis, and writing of the manuscript. SR, CM, and SS contributed to methodology, review, and editing. GM contributed to data curation, review, and editing. MD was responsible for the experiments, data curation, review, and editing. All authors contributed to the article and approved the submitted version.

Funding

The work was supported partly by the LOEWE excellence cluster FACE2FACE of the Hessian State Ministry of Higher Education, Research and the Arts.

Acknowledgments

We thank Rita Geissler-Plaum and Bellinda Schneider for their excellent technical support.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.937021/full#supplementary-material

References

Aitchison, J. (1982). The statistical analysis of compositional data. J. R. Stat. Soc. Ser. B 44, 139–177. doi: 10.1111/j.2517-6161.1982.tb01195.x

CrossRef Full Text | Google Scholar

Aitchison, J. (1986). Book Review. XII. London, NY: Chapman and Hall

Google Scholar

Alfreider, A., Grimus, V., Luger, M., Ekblad, A., Salcher, M. M., and Summerer, M. (2018). Autotrophic carbon fixation strategies used by nitrifying prokaryotes in freshwater lakes. FEMS Microbiol. Ecol. 94, 1–12. doi: 10.1093/femsec/fiy163

PubMed Abstract | CrossRef Full Text | Google Scholar

Anderson, M. J. (2001). A new method for non parametric multivariate analysis of variance. Austral Ecol. 26, 32–46. doi: 10.1046/j.1442-9993.2001.01070.x

CrossRef Full Text | Google Scholar

Andresen, L. C., Yuan, N., Seibert, R., Moser, G., Kammann, C. I., Luterbacher, J., et al. (2018). Biomass responses in a temperate european grassland through 17 years of elevated CO2. Glob. Chang. Biol. 24, 3875–3885. doi: 10.1111/gcb.13705

PubMed Abstract | CrossRef Full Text | Google Scholar

Bakken, L. R., Bergaust, L., Liu, B., and Frostegård, Å. (2012). Regulation of denitrification at the cellular level: a clue to the understanding of N2O emissions from soils. Philos. Trans. R. Soc. B Biol. Sci. 367, 1226–1234. doi: 10.1098/rstb.2011.0321

PubMed Abstract | CrossRef Full Text | Google Scholar

Bashiardes, S., Zilberman-Schapira, G., and Elinav, E. (2016). Use of metatranscriptomics in microbiome research. Bioinform. Biol. Insights 10, 19–25. doi: 10.4137/BBI.S34610

PubMed Abstract | CrossRef Full Text | Google Scholar

Bei, Q., Moser, G., Wu, X., Müller, C., and Liesack, W. (2019). Metatranscriptomics reveals climate change effects on the rhizosphere microbiomes in european grassland. Soil Biol. Biochem. 138, 1–10. doi: 10.1016/j.soilbio.2019.107604

CrossRef Full Text | Google Scholar

Berg, I. A. (2011). Ecological aspects of the distribution of different autotrophic CO2 fixation pathways. Appl. Environ. Microbiol. 77, 1925–1936. doi: 10.1128/AEM.02473-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Blagodatskaya, E., and Kuzyakov, Y. (2008). Mechanisms of real and apparent priming effects and their dependence on soil microbial biomass and community structure: Critical review. Biol. Fertil. Soils 45, 115–131. doi: 10.1007/s00374-008-0334-y

CrossRef Full Text | Google Scholar

Brenzinger, K., Kujala, K., Horn, M. A., Moser, G., Guillet, C., Kammann, C., et al. (2017). Soil conditions rather than long-term exposure to elevated CO2 affect soil microbial communities associated with N-cycling. Front. Microbiol. 8, 1976. doi: 10.3389/fmicb.2017.01976

PubMed Abstract | CrossRef Full Text | Google Scholar

Buchanan, B. B., and Arnon, D. I. (1990). A reverse KREBS cycle in photosynthesis: consensus at last. Photosynth. Res. 24, 47–53. doi: 10.1007/BF00032643

PubMed Abstract | CrossRef Full Text | Google Scholar

Buchanan, B. B., Sirevåg, R., Fuchs, G., Ivanovsky, R. N., Igarashi, Y., Ishii, M., et al. (2017). The Arnon–Buchanan cycle: a retrospective, 1966–2016. Photosynth. Res. 134, 117–131. doi: 10.1007/s11120-017-0429-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Buchfink, B., Xie, C., and Huson, D. H. (2015). Fast and sensitive protein alignment using DIAMOND. Nat. Methods 12, 59–60. doi: 10.1038/nmeth.3176

PubMed Abstract | CrossRef Full Text | Google Scholar

Carini, P., Marsden, P. J., Leff, J. W., Morgan, E. E., Strickland, M. S., and Fierer, N. (2016). Relic DNA is abundant in soil and obscures estimates of soil microbial diversity. Nat. Microbiol. 2, 1–6. doi: 10.1101/043372

PubMed Abstract | CrossRef Full Text | Google Scholar

Carney, K. M., Hungate, B. A., Drake, B. G., and Megonigal, J. P. (2007). Altered soil microbial community at elevated CO2 leads to loss of soil carbon. Proc. Natl. Acad. Sci. U. S. A. 104, 4990–4995. doi: 10.1073/pnas.0610045104

PubMed Abstract | CrossRef Full Text | Google Scholar

Caspi, R., Billington, R., Fulcher, C. A., Keseler, I. M., Kothari, A., Krummenacker, M., et al. (2018). The MetaCyc database of metabolic pathways and enzymes. Nucleic Acids Res. 46, D633–D639. doi: 10.1093/nar/gkx935

PubMed Abstract | CrossRef Full Text | Google Scholar

Clark, K., Karsch-Mizrachi, I., Lipman, D. J., Ostell, J., and Sayers, E. W. (2016). GenBank. Nucleic Acids Res. 44, D67–D72. doi: 10.1093/nar/gkv1276

PubMed Abstract | CrossRef Full Text | Google Scholar

de Menezes, A. B., Müller, C., Clipson, N., and Doyle, E. (2016). The soil microbiome at the Gi-FACE experiment responds to a moisture gradient but not to CO2 enrichment. Microbiology 162, 1572–1582. doi: 10.1099/mic.0.000341

PubMed Abstract | CrossRef Full Text | Google Scholar

Derrien, D., Plain, C., Courty, P. E., Gelhaye, L., Moerdijk-Poortvliet, T. C. W., Thomas, F., et al. (2014). Does the addition of labile substrate destabilise old soil organic matter? Soil Biol. Biochem. 76, 149–160. doi: 10.1016/j.soilbio.2014.04.030

CrossRef Full Text | Google Scholar

Dlott, G., Maul, J. E., Buyer, J., and Yarwood, S. (2015). Microbial rRNA: RDNA gene ratios may be unexpectedly low due to extracellular DNA preservation in soils. J. Microbiol. Methods 115, 112–120. doi: 10.1016/j.mimet.2015.05.027

PubMed Abstract | CrossRef Full Text | Google Scholar

DOE.2020 U. S. (2020). U.S. Department of Energy Free-Air CO2 enrichment experiments: FACE results, lessons, and legacy. DOE/SC−0202.

Dong, J., Hunt, J., Delhaize, E., Zheng, S. J., Jin, C. W., and Tang, C. (2021). Impacts of elevated CO2 on plant resistance to nutrient deficiency and toxic ions via root exudates: A review. Sci. Total Environ. 754, 142434. doi: 10.1016/j.scitotenv.2020.142434

PubMed Abstract | CrossRef Full Text | Google Scholar

Doxtader, K. G., and Alexander, M. (1966). Nitrification by growing and replacement cultures of Aspergillus. Can. J. Microbiol. 12, 807–815. doi: 10.1139/m66-109

PubMed Abstract | CrossRef Full Text | Google Scholar

Du, Y., Guo, X., Li, J., Liu, Y., Luo, J., Liang, Y., et al. (2022). Elevated carbon dioxide stimulates nitrous oxide emission in agricultural soils: A global meta-analysis. Pedosphere 32, 3–14. doi: 10.1016/S1002-0160(21)60057-7

CrossRef Full Text | Google Scholar

Duin, J., and Olsthoorn, R. C. L. (2012). “Family Leviviridae,” in Virus Taxonomy: Ninth Report of the International Committee on Taxonomy of Viruses, eds. A. M. Q. King, M. J. Adams, E. B. Carstens, and E. J. Lefkowitz (London: Academic Press) 1035–1043. doi: 10.1016/B978-0-12-384684-6.00089-6

CrossRef Full Text

Eddy, S. R. (2009). a New generation of homology search tools based on probabilistic inference. Genome Inform. 23, 205–211. doi: 10.1142/9781848165632_0019

PubMed Abstract | CrossRef Full Text | Google Scholar

Evans, M. C., Buchanan, B. B., and Arnon, D. I. (1966). A new ferredoxin-dependent carbon reduction cycle in a photosynthetic bacterium. Proc. Natl. Acad. Sci. 55, 928–934. doi: 10.1073/pnas.55.4.928

PubMed Abstract | CrossRef Full Text | Google Scholar

Fernandes, A. D., Macklaim, J. M., Linn, T. G., Reid, G., and Gloor, G. B. (2013). ANOVA-like differential expression (ALDEx) analysis for mixed population RNA-Seq. PLoS ONE 8, e67019. doi: 10.1371/journal.pone.0067019

PubMed Abstract | CrossRef Full Text | Google Scholar

Fernandes, A. D., Reid, J. N. S., Macklaim, J. M., McMurrough, T. A., Edgell, D. R., and Gloor, G. B. (2014). Unifying the analysis of high-throughput sequencing datasets: Characterizing RNA-seq, 16S rRNA gene sequencing and selective growth experiments by compositional data analysis. Microbiome 2, 1–13. doi: 10.1186/2049-2618-2-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Finn, R. D., Coggill, P., Eberhardt, R. Y., Eddy, S. R., Mistry, J., Mitchell, A. L., et al. (2016). The Pfam protein families database: Towards a more sustainable future. Nucleic Acids Res. 44, D279–D285. doi: 10.1093/nar/gkv1344

PubMed Abstract | CrossRef Full Text | Google Scholar

Fontaine, S., Bardoux, G., Abbadie, L., and Mariotti, A. (2004). Carbon input to soil may decrease soil carbon content. Ecol. Lett. 7, 314–320. doi: 10.1111/j.1461-0248.2004.00579.x

CrossRef Full Text | Google Scholar

Franzosa, E. A., Morgan, X. C., Segata, N., Waldron, L., Reyes, J., Earl, A. M., et al. (2014). Relating the metatranscriptome and metagenome of the human gut. Proc. Natl. Acad. Sci. U. S. A. 111, E2329–38. doi: 10.1073/pnas.1319284111

PubMed Abstract | CrossRef Full Text | Google Scholar

Habash, D. Z., Paul, M. J., Parry, M. A. J., Keys, A. J., and Lawlor, D. W. (1995). Increased capacity for photosynthesis in wheat grown at elevated CO2 - the relationship between electron transport and carbon metabolism. Planta 197, 482–489. doi: 10.1007/BF00196670

CrossRef Full Text | Google Scholar

Hayden, H. L., Mele, P. M., Bougoure, D. S., Allan, C. Y., Norng, S., Piceno, Y. M., et al. (2012). Changes in the microbial community structure of bacteria, archaea and fungi in response to elevated CO2 and warming in an Australian native grassland soil. Environ. Microbiol. 14, 3081–3096. doi: 10.1111/j.1462-2920.2012.02855.x

PubMed Abstract | CrossRef Full Text | Google Scholar

He, P., Bader, K. P., Radunz, A., and Schmid, G. H. (1995). Consequences of high CO2-concentrations in air on growth and gas-exchange rates in Tobacco mutants. Zeitschrift für Naturforsch. 50c, 781–788. doi: 10.1515/znc-1995-11-1208

CrossRef Full Text | Google Scholar

He, Z., Xiong, J., Kent, A. D., Deng, Y., Xue, K., Wang, G., et al. (2014). Distinct responses of soil microbial communities to elevated CO2 and O3 in a soybean agro-ecosystem. ISME J. 8, 714–726. doi: 10.1038/ismej.2013.177

PubMed Abstract | CrossRef Full Text | Google Scholar

He, Z., Xu, M., Deng, Y., Kang, S., Kellogg, L., Wu, L., et al. (2010). Metagenomic analysis reveals a marked divergence in the structure of belowground microbial communities at elevated CO2. Ecol. Lett. 13, 564–575. doi: 10.1111/j.1461-0248.2010.01453.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Hendrix, R. W., Casjens, S. R., and Lavigne, R. (2012). “Family Siphoviridae,” in Virus Taxonomy: Ninth Report of the International Committee on Taxonomy of Viruses, eds. A. M. Q. King, M. J. Adams, E. B. Carstens, and E. J. Lefkowitz (London: Academic Press) 86–98. doi: 10.1016/B978-0-12-384684-6.00004-5

CrossRef Full Text | Google Scholar

Hong, Y., Dover, S. L., Cole, T. E., Brasier, C. M., and Buck, K. W. (1999). Multiple mitochondrial viruses in an isolate of the Dutch elm disease fungus Ophiostoma novo-ulmi. Virology 258, 118–127. doi: 10.1006/viro.1999.9691

PubMed Abstract | CrossRef Full Text

Hora, T. S., and Iyengar, M. R. S. (1960). Nitrification by soil fungi. Arch. Mikrobiol. 35, 252–257. doi: 10.1007/BF00446826

PubMed Abstract | CrossRef Full Text | Google Scholar

Huerta-Cepas, J., Szklarczyk, D., Forslund, K., Cook, H., Heller, D., Walter, M. C., et al. (2016). eggNOG 4.5: A hierarchical orthology framework with improved functional annotations for eukaryotic, prokaryotic and viral sequences. Nucleic Acids Res. 44, D286–D293. doi: 10.1093/nar/gkv1248

PubMed Abstract | CrossRef Full Text | Google Scholar

Hyatt, D., Chen, G.-L., LoCascio, P. F., Land, M. L., Larimer, F. W., and Hauser, L. J. (2010). Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinform. 11, 1–11. doi: 10.1186/1471-2105-11-119

PubMed Abstract | CrossRef Full Text | Google Scholar

Idso, K. (1994). Plant responses to atmospheric CO2 enrichment in the face of environmental constraints: a review of the past 10 years' research. Agric. For. Meteorol. 69, 153–203. doi: 10.1016/0168-1923(94)90025-6

CrossRef Full Text | Google Scholar

Illumina (2019). bcl2fastq2 Conversion Software v2.20. Available online at: https://support.illumina.com/downloads/bcl2fastq-conversion-software-v2-20.html (accessed March 25, 2022).

IPCC (2021). “Technical Summary,” in Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change. eds. V. Masson-Delmotte, P. Zhai, A. Pirani, S. L. Connors, C. Péan, S. Berger, et al. Cambridge: Cambridge University Press.

Jäger, H. J., Schmidt, S. W., Kammann, C., Grunhage, L., Muller, C., and Hanewald, K. (2003). The university of Giessen Free-Air Carbon Dioxide Enrichment study: Description of the experimental site and of a new enrichment system. J. Appl. Bot. Bot. 77, 117–127.

Google Scholar

Jia, X., Wang, W., Chen, Z., He, Y., and Liu, J. (2014). Concentrations of secondary metabolites in tissues and root exudates of wheat seedlings changed under elevated atmospheric CO2 and cadmium-contaminated soils. Environ. Exp. Bot. 107, 134–143. doi: 10.1016/j.envexpbot.2014.06.005

CrossRef Full Text | Google Scholar

Johnson, R. M., and Pregitzer, K. S. (2007). Concentration of sugars, phenolic acids, and amino acids in forest soils exposed to elevated atmospheric CO2 and O3. Soil Biol. Biochem. 39, 3159–3166. doi: 10.1016/j.soilbio.2007.07.010

CrossRef Full Text | Google Scholar

Jolliffe, I. T., and Cadima, J. (2016). Principal component analysis: a review and recent developments Subject Areas. Phil. Trans. R. Soc. A 374, 1–16. doi: 10.1098/rsta.2015.0202

PubMed Abstract | CrossRef Full Text | Google Scholar

Jongen, M., Jones, M. B., Hebeisen, T., Blum, H., and Hendrey, G. (1995). The effects of elevated CO2 concentrations on the root growth of Lolium perenne and Trifolium repens grown in a FACE* system. Glob. Chang. Biol. 1, 361–371. doi: 10.1111/j.1365-2486.1995.tb00034.x

CrossRef Full Text | Google Scholar

Kammann, C., Müller, C., Grünhage, L., and Jäger, H. J. (2008). Elevated CO2 stimulates N2O emissions in permanent grassland. Soil Biol. Biochem. 40, 2194–2205. doi: 10.1016/j.soilbio.2008.04.012

CrossRef Full Text | Google Scholar

Kanehisa, M., and Goto, S. (2000). KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 28, 27–30. doi: 10.1093/nar/28.1.27

PubMed Abstract | CrossRef Full Text | Google Scholar

Kang, D. D., Li, F., Kirton, E., Thomas, A., Egan, R., An, H., et al. (2019). MetaBAT 2: An adaptive binning algorithm for robust and efficient genome reconstruction from metagenome assemblies. PeerJ. 2019, 1–13. doi: 10.7717/peerj.7359

PubMed Abstract | CrossRef Full Text

Kimball, B. A. (1983). Carbon dioxide and agricultural yield: An assemblage and analysis of 430 prior observations. Agron. J. 75, 779. doi: 10.2134/agronj1983.00021962007500050014x

CrossRef Full Text | Google Scholar

Kimball, B. A. (2016). Crop responses to elevated CO2 and interactions with H2O, N, and temperature. Curr. Opin. Plant Biol. 31, 36–43. doi: 10.1016/j.pbi.2016.03.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Krishnamurthy, S. R. (2017). Expansion of microbial virology by impetus of the reduction of viral dark matter. Arts Sci. Electron. Theses Diss. doi: 10.7936/K7VX0FXS Available online at: https://openschoolarship.wustl.edu/art_sci_etds/1221

CrossRef Full Text

Krupovic, M., Dolja, V. V., and Koonin, E. V. (2020). The LUCA and its complex virome. Nat. Rev. Microbiol. 18, 661–670. doi: 10.1038/s41579-020-0408-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Lahti, L., and Shetty, S. (2019). microbiome R package. Available online at: http://microbiome.github.io (accessed January 14, 2022).

Google Scholar

Lane, D. J. (1991). “16S/23S rRNA sequencing.,” in Nucleic Acid Techniques in Bacterial Systematics, eds. E. Stackebrandt and M. Goodfellow (New York: John Wiley and Sons) 115–175.

Google Scholar

Langmead, B., and Salzberg, S. L. (2012). Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359. doi: 10.1038/nmeth.1923

PubMed Abstract | CrossRef Full Text | Google Scholar

Laslett, D., and Canback, B. (2004). ARAGORN, a program to detect tRNA genes and tmRNA genes in nucleotide sequences. Nucleic Acids Res. 32, 11–16. doi: 10.1093/nar/gkh152

PubMed Abstract | CrossRef Full Text | Google Scholar

Laughlin, R. J., Stevens, R. J., Müller, C., and Watson, C. J. (2008). Evidence that fungi can oxidize NH inline math to NO inline math in a grassland soil. Eur. J. Soil Sci. 59, 285–291. doi: 10.1111/j.1365-2389.2007.00995.x

CrossRef Full Text

Le Roux, X., Bouskill, N. J., Niboyet, A., Barthes, L., Dijkstra, P., Field, C. B., et al. (2016). Predicting the responses of soil nitrite-oxidizers to multi-factorial global change: A trait-based approach. Front. Microbiol. 7, 628. doi: 10.3389/fmicb.2016.00628

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, S. H., and Kang, H. (2016). Elevated CO2 causes a change in microbial communities of rhizosphere and bulk soil of salt marsh system. Appl. Soil Ecol. 108, 307–314. doi: 10.1016/j.apsoil.2016.09.009

CrossRef Full Text | Google Scholar

Lee, S. H., Kim, S. Y., Ding, W., and Kang, H. (2015). Impact of elevated CO2 and N addition on bacteria, fungi, and archaea in a marsh ecosystem with various types of plants. Appl. Microbiol. Biotechnol. 99, 5295–5305. doi: 10.1007/s00253-015-6385-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, D., Liu, C. M., Luo, R., Sadakane, K., and Lam, T. W. (2015). MEGAHIT: An ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics 31, 1674–1676. doi: 10.1093/bioinformatics/btv033

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, B., Mørkved, P. T., Frostegård, Å., and Bakken, L. R. (2010). Denitrification gene pools, transcription and kinetics of NO, N2O and N2 production as affected by soil pH. FEMS Microbiol. Ecol. 72, 407–417. doi: 10.1111/j.1574-6941.2010.00856.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, X. J. A., Sun, J., Mau, R. L., Finley, B. K., Compson, Z. G., van Gestel, N., et al. (2017). Labile carbon input determines the direction and magnitude of the priming effect. Appl. Soil Ecol. 109, 7–13. doi: 10.1016/j.apsoil.2016.10.002

CrossRef Full Text | Google Scholar

Marshall, K. C., and Alexander, M. (1962). Nitrification by Aspergillus flavus. J. Bacteriol. 83, 572–578. doi: 10.1128/jb.83.3.572-578.1962

PubMed Abstract | CrossRef Full Text | Google Scholar

McMurdie, P. J., and Holmes, S. (2013). Phyloseq: An R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE 8, e61217. doi: 10.1371/journal.pone.0061217

PubMed Abstract | CrossRef Full Text | Google Scholar

Mettel, C., Kim, Y., Shrestha, P. M., and Liesack, W. (2010). Extraction of mRNA from soil. Appl. Environ. Microbiol. 76, 5995–6000. doi: 10.1128/AEM.03047-09

PubMed Abstract | CrossRef Full Text | Google Scholar

Morrissey, E. M., McHugh, T. A., Preteska, L., Hayer, M., Dijkstra, P., Hungate, B. A., et al. (2015). Dynamics of extracellular DNA decomposition and bacterial community composition in soil. Soil Biol. Biochem. 86, 42–49. doi: 10.1016/j.soilbio.2015.03.020

CrossRef Full Text | Google Scholar

Moser, G., Gorenflo, A., Brenzinger, K., Keidel, L., Braker, G., Marhan, S., et al. (2018). Explaining the doubling of N2O emissions under elevated CO2 in the Giessen FACE via in-field 15N tracing. Glob. Chang. Biol. 24, 3897–3910. doi: 10.1111/gcb.14136

PubMed Abstract | CrossRef Full Text | Google Scholar

Müller, C., Rütting, T., Abbasi, M. K., Laughlin, R. J., Kammann, C., Clough, T. J., et al. (2009). Effect of elevated CO2 on soil N dynamics in a temperate grassland soil. Soil Biol. Biochem. 41, 1996–2001. doi: 10.1016/j.soilbio.2009.07.003

CrossRef Full Text | Google Scholar

Mundinger, A. B., Lawson, C. E., Jetten, M. S. M., Koch, H., and Lücker, S. (2019). Cultivation and transcriptional analysis of a canonical nitrospira under stable growth conditions. Front. Microbiol. 10, 1325. doi: 10.3389/fmicb.2019.01325

PubMed Abstract | CrossRef Full Text | Google Scholar

NASA (2022). SVS: Atmospheric carbon dioxide concentrations. Available online at: https://svs.gsfc.nasa.gov/4962 (accessed February 6, 2022).

Google Scholar

Niklaus, P. A., and Falloon, P. (2006). Estimating soil carbon sequestration under elevated CO2 by combining carbon isotope labelling with soil carbon cycle modelling. Glob. Chang. Biol. 12, 1909–1921. doi: 10.1111/j.1365-2486.2006.01215.x

CrossRef Full Text | Google Scholar

Nuccio, E. E., Nguyen, N. H., Nunes da Rocha, U., Mayali, X., Bougoure, J., Weber, P. K., et al. (2021). Community RNA-Seq: multi-kingdom responses to living versus decaying roots in soil. ISME Commun. 1, 1–10. doi: 10.1038/s43705-021-00059-3

CrossRef Full Text | Google Scholar

Owensby, C. E., Ham, J. M., Knapp, A. K., Bremer, D., and Auen, L. M. (1997). Water vapour fluxes and their impact under elevated CO2 in a C4-tallgrass prairie. Glob. Chang. Biol. 3, 189–195. doi: 10.1046/j.1365-2486.1997.00084.x

CrossRef Full Text | Google Scholar

Padhy, S. R., Bhattacharyya, P., Nayak, A. K., Dash, P. K., Roy, K. S., Baig, M. J., et al. (2020). Key metabolic pathways of sulfur metabolism and bacterial diversity under elevated CO2 and temperature in lowland rice: A metagenomic approach. Geomicrobiol. J. 37, 13–21. doi: 10.1080/01490451.2019.1657992

CrossRef Full Text | Google Scholar

Parks, D. H., Imelfort, M., Skennerton, C. T., Hugenholtz, P., and Tyson, G. W. (2015). CheckM: Assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 25, 1043–1055. doi: 10.1101/gr.186072.114

PubMed Abstract | CrossRef Full Text | Google Scholar

Phillips, R. P., Meier, I. C., Bernhardt, E. S., Grandy, A. S., Wickings, K., and Finzi, A. C. (2012). Roots and fungi accelerate carbon and nitrogen cycling in forests exposed to elevated CO2. Ecol. Lett. 15, 1042–1049. doi: 10.1111/j.1461-0248.2012.01827.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Pratama, A. A., and van Elsas, J. D. (2018). The ‘Neglected' Soil virome – potential role and impact. Trends Microbiol. 26, 649–662. doi: 10.1016/j.tim.2017.12.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Puente-Sánchez, F., García-García, N., and Tamames, J. (2020). SQMtools: automated processing and visual analysis of'omics data with R and anvi'o. BMC Bioinformatics 21, 358. doi: 10.1186/s12859-020-03703-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiao, N., Schaefer, D., Blagodatskaya, E., Zou, X., Xu, X., and Kuzyakov, Y. (2014). Labile carbon retention compensates for CO2 released by priming in forest soils. Glob. Chang. Biol. 20, 1943–1954. doi: 10.1111/gcb.12458

PubMed Abstract | CrossRef Full Text | Google Scholar

Regan, K., Kammann, C., Hartung, K., Lenhart, K., Müller, C., Philippot, L., et al. (2011). Can differences in microbial abundances help explain enhanced N2O emissions in a permanent grassland under elevated atmospheric CO2? Glob. Chang. Biol. 17, 3176–3186. doi: 10.1111/j.1365-2486.2011.02470.x

CrossRef Full Text | Google Scholar

Rosado-Porto, D., Ratering, S., Cardinale, M., Maisinger, C., Moser, G., Deppe, M., et al. (2021). Elevated atmospheric co2 modifies mostly the metabolic active rhizosphere soil microbiome in the Giessen FACE experiment. Microb. Ecol. doi: 10.1007/s00248-021-01791-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Sawilowsky, S. S. (2009). New effect size rules of thumb. J. Mod. Appl. Stat. Methods 8, 597–599. doi: 10.22237/jmasm/1257035100

CrossRef Full Text | Google Scholar

Schmieder, R., and Edwards, R. (2011). Quality control and preprocessing of metagenomic datasets. Bioinformatics 27, 863–864. doi: 10.1093/bioinformatics/btr026

PubMed Abstract | CrossRef Full Text | Google Scholar

Seemann, T. (2014). Prokka: Rapid prokaryotic genome annotation. Bioinformatics 30, 2068–2069. doi: 10.1093/bioinformatics/btu153

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharma, R., and Sharma, P. K. (2018). Metatranscriptome sequencing and analysis of agriculture soil provided significant insights about the microbial community structure and function. Ecol. Genet. Genomics 6, 9–15. doi: 10.1016/j.egg.2017.10.001

CrossRef Full Text | Google Scholar

Simonin, M., Le Roux, X., Poly, F., Lerondelle, C., Hungate, B. A., Nunan, N., et al. (2015). Coupling between and among ammonia oxidizers and nitrite oxidizers in grassland mesocosms Submitted to elevated CO2 and nitrogen supply. Microb. Ecol. 70, 809–818. doi: 10.1007/s00248-015-0604-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Steffens, L., Pettinato, E., Steiner, T. M., Mall, A., König, S., Eisenreich, W., et al. (2021). High CO2 levels drive the TCA cycle backwards towards autotrophy. Nature 592, 784–788. doi: 10.1038/s41586-021-03456-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Tamames, J., and Puente-Sánchez, F. (2019). SqueezeMeta, a highly portable, fully automatic metagenomic analysis pipeline. Front. Microbiol. 10, 1–10. doi: 10.3389/fmicb.2018.03349

PubMed Abstract | CrossRef Full Text | Google Scholar

Trivedi, P., Delgado-Baquerizo, M., Trivedi, C., Hu, H., Anderson, I. C., Jeffries, T. C., et al. (2016). Microbial regulation of the soil carbon cycle: Evidence from gene-enzyme relationships. ISME J. 10, 2593–2604. doi: 10.1038/ismej.2016.65

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Groenigen, K. J., Gorissen, A., Six, J., Harris, D., Kuikman, P. J., Van Groenigen, J. W., et al. (2005). Decomposition of 14C-labeled roots in a pasture soil exposed to 10 years of elevated CO2. Soil Biol. Biochem. 37, 497–506. doi: 10.1016/j.soilbio.2004.08.013

CrossRef Full Text | Google Scholar

Vestergard, M., Reinsch, S., Bengtson, P., Ambus, P., and Christensen, S. (2016). Enhanced priming of old, not new soil carbon at elevated atmospheric CO2. Soil Biol. Biochem. 100, 140–148. doi: 10.1016/j.soilbio.2016.06.010

CrossRef Full Text | Google Scholar

Weisburg, W. G., Barns, S. M., Pelletier, D. A., and Lane, D. J. (1991). 16S ribosomal DNA amplification for phylogenetic study. J. Bacteriol. 173, 697–703. doi: 10.1128/jb.173.2.697-703.1991

PubMed Abstract | CrossRef Full Text | Google Scholar

Williamson, K. E., Fuhrmann, J. J., Wommack, K. E., and Radosevich, M. (2017). Viruses in soil ecosystems: An unknown quantity within an unexplored territory. Annu. Rev. Virol. 4, 201–219. doi: 10.1146/annurev-virology-101416-041639

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, Y. W., Simmons, B. A., and Singer, S. W. (2016). MaxBin 2.0: An automated binning algorithm to recover genomes from multiple metagenomic datasets. Bioinformatics 32, 605–607. doi: 10.1093/bioinformatics/btv638

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, Z., Cadisch, G., Edwards, G., Baggs, E. M., and Blum, H. (2005). Carbon dynamics in a temperate grassland soil after 9 years exposure to elevated CO2 (Swiss FACE). Soil Biol. Biochem. 37, 1387–1395. doi: 10.1016/j.soilbio.2004.12.010

CrossRef Full Text | Google Scholar

Xiong, J., He, Z., Shi, S., Kent, A., Deng, Y., Wu, L., et al. (2015). Elevated CO2 shifts the functional structure and metabolic potentials of soil microbial communities in a C4 agroecosystem. Sci. Rep. 5, 1–9. doi: 10.1038/srep09316

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, M., He, Z., Deng, Y., Wu, L., Van Nostrand, J. D., Hobbie, S. E., et al. (2013). Elevated CO2 influences microbial carbon and nitrogen cycling. BMC Microbiol. 13, 1–11. doi: 10.1186/1471-2180-13-124

PubMed Abstract | CrossRef Full Text | Google Scholar

Ye, Y., and Doak, T. G. (2009). A parsimony approach to biological pathway reconstruction/inference for genomes and metagenomes. PLoS Comput. Biol. 5, 1–8. doi: 10.1371/journal.pcbi.1000465

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, H., Deng, Y., He, Z., Pendall, E., Carrillo, Y., Wang, S., et al. (2021). Stimulation of soil microbial functioning by elevated CO2 may surpass effects mediated by irrigation in a semiarid grassland. Geoderma 401, 115162. doi: 10.1016/j.geoderma.2021.115162

CrossRef Full Text | Google Scholar

Yu, H., Deng, Y., He, Z., Van Nostrand, J. D., Wang, S., Jin, D., et al. (2018a). Elevated CO2 and warming altered grassland microbial communities in soil top-layers. Front. Microbiol. 9, 1790. doi: 10.3389/fmicb.2018.01790

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, H., He, Z., Wang, A., Xie, J., Wu, L., Van Nostrand, J. D., et al. (2018b). Divergent responses of forest soil microbial communities under elevated CO2 in different depths of upper soil layers. Appl. Environ. Microbiol. 84, e01694–17. doi: 10.1128/AEM.01694-17

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, K., Ma, B., Xu, Y., Stirling, E., and Xu, J. (2021). Light exposure mediates circadian rhythms of rhizosphere microbial communities. ISME J. 15, 2655–2664. doi: 10.1038/s41396-021-00957-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, T., Meng, T., Zhang, J., Zhong, W., Müller, C., and Cai, Z. (2015). Fungi-dominant heterotrophic nitrification in a subtropical forest soil of China. J. Soils Sediments 15, 705–709. doi: 10.1007/s11368-014-1048-4

CrossRef Full Text | Google Scholar

Keywords: elevated CO2, FACE, functional metatranscriptomics, carbon cycle, nitrogen cycle, sulfur cycle, soil microbiome

Citation: Rosado-Porto D, Ratering S, Moser G, Deppe M, Müller C and Schnell S (2022) Soil metatranscriptome demonstrates a shift in C, N, and S metabolisms of a grassland ecosystem in response to elevated atmospheric CO2. Front. Microbiol. 13:937021. doi: 10.3389/fmicb.2022.937021

Received: 05 May 2022; Accepted: 01 August 2022;
Published: 23 August 2022.

Edited by:

Paul Dijkstra, Northern Arizona University, United States

Reviewed by:

Mengting Maggie Yuan, University of California, Berkeley, United States
Rongxiao Che, Yunnan University, China

Copyright © 2022 Rosado-Porto, Ratering, Moser, Deppe, Müller and Schnell. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Sylvia Schnell, Sylvia.Schnell@umwelt.uni-giessen.de

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.