Gileno V. Lacerda-Júnior1,2*
Melline F. Noronha2
Lucélia Cabral2,3
Tiago P. Delforno2
Sanderson Tarciso Pereira de Sousa2
Paulo I. Fernandes-Júnior4*
Itamar S. Melo1Valéria M. Oliveira2*- 1Brazilian Agricultural Research Corporation, Embrapa Meio Ambiente, Jaguariúna, Brazil
- 2Division of Microbial Resources (DRM), Chemical, Biological and Agricultural Pluridisciplinary Research Center (CPQBA), Campinas State University (UNICAMP), Campinas, Brazil
- 3Brazilian Bioethanol Science and Technology Laboratory (CTBE), Brazilian Center for Research in Energy and Materials (CNPEM), Campinas, Brazil
- 4Brazilian Agricultural Research Corporation, Embrapa Semiárido, Petrolina, Brazil
Drylands occupy approximately 41% of the Earth’s terrestrial surface. Climate change and land use practices are expected to affect biogeochemical cycling by the soil microbiome in these ecosystems. Understanding how soil microbial community might respond to these drivers is extremely important to mitigate the processes of land degradation and desertification. The Caatinga, an exclusively Brazilian biome composed of an extensive seasonal tropical dry forest, is exposed to variable spatiotemporal rainfall patterns as well as strong human-driven pressures. Herein, an integrated analysis of shotgun metagenomics approach coupled to meteorological data was employed to unravel the impact of seasonality and land use change on soil microbiome from preserved and agriculture-affected experimental fields in Caatinga drylands. Multivariate analysis suggested that microbial communities of preserved soils under seasonal changes were shaped primarily by water deficit, with a strong increase of Actinobacteria and Proteobacteria members in the dry and rainy seasons, respectively. In contrast, nutrient availability notably played a critical role in driving the microbial community in agriculture-affected soils. The strong enrichment of bacterial genera belonging to the poorly-known phylum Acidobacteria (‘Candidatus Solibacter’ and ‘Candidatus Koribacter’) in soils from dry season affected by ferti-irrigation practices presupposes a contrasting copiotrophic lifestyle and ecological role in mitigating the impact of chemical fertilization. Functional analyses identify overrepresented genes related to osmotic stress response (synthesis of osmoprotectant compounds, accumulation of potassium ions) and preferential carbon and nitrogen utilization when comparing the microbiome of preserved soils under seasonal changes, reflecting differences in the genetic potential for nutrient cycling and C acquisition in the environment. However, the prevalence of nitrosative stress and denitrification functions in irrigation/fertilization-affected soils of the dry season clearly suggest that nutrient input and disruption of natural water regime may impact biogeochemical cycles linked to the microbial processes, with potential impacts on the ecosystem functionality. These findings help to better understand how natural seasonality and agricultural management differentially affect soil microbial ecology from dry forests, providing support for the development of more sustainable land management in dryland ecosystems.
Introduction
In terrestrial ecosystems, soil microbes play an important role in maintaining structure, decomposing organic matter and nutrient cycling, sequestering of carbon (C) and moderating of climate (Bardgett and Van Der Putten, 2014). With increasing global desertification due to climate change and human-driven practices, it is important to understand the mechanisms that enable soil microbiota to cope with external factors, and their effects on key microbial processes of biogeochemical cycles (Drenovsky et al., 2004; Andrew et al., 2012; Shade et al., 2013).
However, the assessment of soil microbial diversity and their functional contributions is a challenge since 1 g of surface soil may harbor billions of microbes encompassing millions of individual species (Hughes et al., 2001). Recently, advances in culture-independent methods and DNA sequencing technologies have provided a deep understanding of the biotic and abiotic parameters affecting soil microbial community composition; overcoming the limitations of culture-based approaches. Particularly in arid and semi-arid environments, the temperature and water availability have been shown exert over-control on soil microbial communities (Stres et al., 2008; Pasternak et al., 2013; Nielsen and Ball, 2015; Armstrong et al., 2016; Zhao et al., 2016). Short-term water status water status is also an essential factor regulating the microbial activity and microbial community composition (evaluated by PLFA markers) in paddy soils varying in pH, soil organic matter, and soil texture (Liao et al., 2018).
The effects of land use changes and management on soil microbial communities have been widely studied in forests and grassland ecosystems, frequently comparing between tilled or fertilized systems (Rasche et al., 2011; Lauber et al., 2013; Navarrete et al., 2013; Rodrigues et al., 2013; Nacke et al., 2014; Thapa et al., 2018). In dryland regions, the crop production is one of the most devastating land use type nourishing one-third of the global population, and may significantly influence soil health and quality (United Nations Environmental Management Group, 2011). The agricultural management inputs, such as fertilizer, herbicide, and irrigation impact soil microbial diversity and the ecosystem functioning (Ding et al., 2013). Due these water-limited ecosystems cover about 41.5% of the Earth’s land surface (Sorensen, 2009) and have undeniable importance in global biogeochemical cycles (Maestre et al., 2016), studies investigating the effects of environmental parameters, as climate, land use type, and management on the soil microbial communities functions have been requested (Tian et al., 2017; Lüneberg et al., 2018; Pajares et al., 2018). In soils from semi-arid region of Mezquital Valley (P/PET 0.32) the irrigation with different water quantity and quality (freshwater, untreated wastewater, and untreated wastewater) differentially impact the entire and potentially active bacterial community. RNA-based analysis revealed gene abundances involved in nitrogen, carbon and phosphorous cycles are severely affected among different land use systems and season.
These studies become extremely important as the most recent climatic projections suggest an increase of 11–23% in the extent of global drylands by the end of this century (Huang et al., 2016). There is a need to understand how land-use changes affect soil microbial ecology in order to predict ecosystem stability for the development of more sustainable land management in dryland ecosystems.
Seasonal tropical dry forests (STDF) are widely distributed (up to 40% of all tropical forest) and the most threatened and least studied of the world’s forested ecosystems (Murphy and Lugo, 1986; Gillespie et al., 2012). “Caatinga,” an exclusively Brazilian biome composed of the largest STDF in South America, is located in the semi-arid northeastern Brazil (covering about 11% of the territory) harboring more than 23 million of people, indicating high human-driven pressure in the natural environment (Beuchle et al., 2015). In this region, the rainfall is low (below 800 mm year−1) and its distribution in the time leads to two well-defined seasons, rainy and dry (Salgado et al., 2015). The Caatinga eco-physiognomy presents a heterogeneous mosaic of thorny trees and shrubs with xerophytic survival characteristics highly adapted to water shortage. For the majority of the tree species the leaves fall in the dry season and white tree trunks and shrubs remain in the landscape (Salgado et al., 2015). Despite harboring great biodiversity and a high level of endemic species, this ecosystem is still poorly studied and preserved (Santos et al., 2011; Pacca et al., 2018). Anthropogenic processes (i.e., agriculture, livestock and predatory extractive) as well as climatic changes have caused worrying environmental damage, such as accelerated desertification and impacts on C and N cycles (de Albuquerque et al., 2012; Menezes et al., 2012; Santos et al., 2014).
Increased attention has recently been paid to the cultivation, characterization and biotechnological potential of Caatinga microbes due to their unique biological features that allow them to survive under severe climate conditions (high temperature, high ultraviolet exposure, and water shortage) (Monteiro et al., 2009; Marcon et al., 2010; Castro et al., 2014; Fernandes-Júnior et al., 2015). Metagenomics-based approaches and 16S rRNA amplicon sequencing have sought to assess the great microbial diversity and shed light on the metabolic strategies that enable the Caatinga microbes to survive under harsh conditions (Lançoni et al., 2013; Nessner Kavamura et al., 2013; Pacchioni et al., 2014; Lacerda Júnior et al., 2017; Leite et al., 2017). Although some studies have shown that soil microbial communities are shaped by seasonality, it is still not clear how microbes adapt to the conversion of intact Caatinga dry forest into cultivation fields, in addition to the underlying impact on the ecosystem functioning.
In this study, the area under investigation is located along São Francisco river valley and includes one preserved area of Caatinga dry forest (Caatinga experimental field) and another one consisting of Caatinga fragments (Bebedouro experimental field) surrounded by agricultural lands and consequently affected by water and nutrients input from cropped areas. These areas represent natural environmental models for investigating seasonality and anthropogenic stressors outlining soil microbial communities under semi-arid conditions. In this context, we hypothesized that: (1) environmental changes associated with seasonal cycles are the major factor outlining the taxonomic and functional traits of the microbiome in preserved soils from Caatinga dry forest; and (2) water and nutrient input by land use practices affect the natural composition and functional capability of soil microbiome, with potential impacts on biogeochemical processes and soil functionality. To test these hypotheses, a field survey was performed combining meteorological and soil chemical parameters with deep metagenome sequencing of the Caatinga dryland soils.
Materials and Methods
Study Area and Sampling
Bulk soil samples were collected at the Caatinga experimental field (CEF) and Bebedouro experimental field (BEF), located in São Francisco River Valley (semi-arid region of northeastern Brazil) at Petrolina city (Pernambuco state, Brazil). Geographic coordinates of the sampling sites and soil characteristics are listed in Supplementary Table S1. CEF is a preserved area of the Caatinga biome managed by Brazilian Agricultural Research Corporation (Embrapa Semiárido, Brazil) widely covered by several dry tolerant species native from Caatinga biome (da Silva et al., 2017; Figure 1A). Within the facilities of Embrapa Semiárido, the CEF had no human interventions at least for the last 43 years, time of the Embrapa settled in the place. BEF encompasses Caatinga fragments surrounded by a management area with irrigation, fertilization and fertirrigation practices during the dry season. The area in BEF is located within an irrigation perimeter with 6000 hectares used for the agricultural crop grown since the 1960s (see Supplementary Figure S1).
Figure 1. (A) Pictures of Caatinga vegetation during dry and rainy season in the semi-arid region of northeastern Brazil (pictures taken by the first author, 2014 and 2015). (B) Surface chart showing the distribution of meteorological data of relative humidity (%), air temperature (°C) and solar Radiation (MJ/m2) collected in the CEF and BEF experimental fields during the years 2014 and 2015. Soil samples were collected in October 2014 and May 2015 during dry and rainy seasons, respectively. The data were obtained from automatic meteorological stations of EMBRAPA.
Soil samples were randomly collected (depth 0- to 10-cm of the topsoil layer) from three different sites (P1, P2, and P3) across the CEF and BEF areas during the spring (October 2014) and winter (May 2015), which corresponded to the peak of the dry and wet seasons, respectively. Each replicate consisted of five subsamples (arranged 50 m away) which were collected and combined in the field. Samples were stored in a sterile plastic bag, kept refrigerated and transported to the laboratory for storage at −80°C for further DNA extraction. Sampling was authorized by the Institute of Environment and Renewable Natural Resources (IBAMA), process number 02001.004527/2011-90.
Soil Chemistry and Meteorological Analysis
Chemical measurements were performed in triplicate for each sample. Soil samples were air-dried and sieved through a 2-mm mesh sieve for determination of organic matter (OM), ammonium (N-NH4+), nitrate (N-NO3-), sulphate (S-SO4), potassium (K), calcium (Ca), magnesium (Mg), aluminum (Al), iron (Fe), copper (Cu), boron (B), zinc (Zn), manganese (Mn), sodium (Na) and pH. Chemical analyses were carried out as described by Raij et al. (2001). For soil ammonium (N-NH4+) and nitrate (N-NO3−) concentrations, values were determined after extraction of 10 g of soil in 50 mL of KCl (2 M), according to the method described by Keeney and Nelson (1982). In order to obtain a detailed description of the climatic conditions during the samplings, meteorological data were collected daily by automatic stations of Agrometeorology Laboratory (Embrapa Semiárido) in the years 2014 and 2015. Temporal data were collected on air temperature (°C), relative humidity (%), global solar radiation (MJ/m2) and rainfall index (mm). Statistical analyses of soil chemical properties were performed by two-way analysis of variance (ANOVA), using the Sisvar 5.0 statistical package (Ferreira, 2011). Statistical differences between means of dry and rainy soils were assessed according to the Bonferroni test (p-value < 0.05). All assumptions required by variance analysis were checked.
Total DNA Extraction and High-Throughput Metagenomic Sequencing
Firstly, soil samples were homogenized by sieving (2 mm) to remove rocks and plant material. Total community DNA was directly extracted using PowerMax Soil DNA Isolation Kit (Mo Bio Laboratories, United States), following manufacturer’s instructions. DNA integrity was checked on 1% agarose gel electrophoresis. DNA quality was measured using NanoDrop spectrophotometer (Thermo Scientific, United States) at A260/280 nm ratio. DNA samples were used to prepare libraries with the NexteraTM DNA Sample Preparation Kit (Illumina®-compatible). Environmental DNA samples from each site were barcoded for shotgun sequencing in one single lane of Illumina HiSeq 2000 platform (2 × 100 bp) at Multi-User Laboratory of Functional Genomics (Piracicaba, Brazil), according to the manufacturer’s instructions.
In silico Bioinformatics Analysis
For taxonomic and functional annotation, raw reads of each one of the metagenomes were submitted to the MG-RAST pipeline version 3.6 (Metagenomics Analysis Server Annotation) (Meyer et al., 2008). Paired reads were merged with a minimum overlap of 8 bp and a maximum difference of 10%. Low-quality reads and artificial duplicate reads were removed using MG-RAST quality control with default quality thresholds. The taxonomic assignment of the high-quality unassembled reads was performed via BLASTP search against SEED database (Overbeek et al., 2005) using the lowest common ancestor (LCA)-based algorithm (maximum e-value cutoff of 1e−5, minimum percentage of identity cutoff of 60% and a minimum alignment length cutoff of 50 base pairs). Functional analysis was performed against subsystems category of SEED database using BLASTP with the same cutoffs.
Statistical Analysis of Metagenome Data
Statistical analyses were performed using the STAMP (Statistical Analysis of Metagenomic Profiles) software version 2.1.3 (Parks and Beiko, 2010) to identify biologically relevant differences. The input data were relative abundance assigned to the different taxa and functional categories of SEED databases for each metagenome. Statistical differences between soils of the seasonal period (dry versus rainy) from CEF (n = 3) and BEF (n = 3) were determined by two-sided Fisher’s Exact test, and Storey’s FDR (false discovery rate) method was used for multiple test corrections, as recommended by the STAMP developers. All features with a p < 0.05 were filtered. Multivariate analyses of soil chemical variables and taxonomic profiles were performed using the PAST software version 3.15 (Hammer et al., 2001).
Results and Discussion
Meteorological Data and Chemical Properties of Soil Samples
The climatic patterns of CEF and BEF in the years 2014 and 2015 were relatively consistent with the expected ones for the well-defined dry and rainy seasons of the Caatinga biome (Salgado et al., 2015). Solar radiation and air temperature were higher during August to October 2014 (corresponding to the dry season), while higher air humidity rates were obtained in May 2015 (corresponding to the rainy season) (Figure 1B). With regards to the rainfall, although the intra-seasonal precipitation is dynamics and varies slightly from year to year, the rains were concentrated from December to June (years 2014 and 2015), while a drought period was observed from August to October with very low or even zero rainfall indexes in both experimental fields (Figure 2A). In this work, the daily meteorological monitoring ensured that soil samples were collected after a long drought period (30 October of 2014) and rainfall events (5 May of 2015), during dry and rainy seasons, respectively (Figure 2B).
Figure 2. (A) Monthly rainfall distribution at the CEF and BEF stations of the Caatinga during the years 2014 and 2015. (B) Daily rainfall data (0–30) during soil sampling on 5 May and 30 October, corresponding to the rainy (main axis) and dry (secondary axis) periods, respectively. The data were obtained from automatic meteorological stations of EMBRAP9A.
Chemical characteristics were determined from rainy and dry soils of preserved (CEF) and anthropized (BEF) Caatinga fragments. The soils of BEF, surrounded by crop lands, showed higher content of macronutrients (C, P, N, Ca, S, K, and Mg) and micronutrients (Mn, Cu, Fe, Al, and B) as well as higher values of parameters related to soil fertility such as cation-exchange capacity (CEC), base saturation (V%) and sum of exchangeable bases (SB) when compared to rainy BEF and pristine CEF soils (Table 1).
Table 1. Chemical properties of soil samples from Caatinga and Bebedouro experimental fields (n = 3).
Water content (%) ratios of 1.67 and 8 were observed between rainy and dry soils from BEF and CEF, respectively (Supplementary Table S2). Likewise, the moisture was the most significantly different (p < 0.05) parameter between rainy and dry preserved soils as determined by F test (Table 1). Despite the similar rainfall index (Figure 2B), BEF soils collected in the dry season presented remarkably higher and similar moisture as compared to preserved soils of dry (average of 0.4 %) and rainy (average of 3.1 %) seasons, respectively (Table 1). These results suggest that BEF is a good experimental field for the investigation of the impacts of natural water regime alterations and land use practices on the soil microbiome.
Impacts on the Soil Microbial Community Structure by the Natural Seasonal Regime and Agriculture Practices
It is well known that several biotic (i.e., competitiveness) and abiotic factors can contribute to shaping microbial community structure over time and space (Drenovsky et al., 2004; Andrew et al., 2012; Shade et al., 2013). In arid and semi-arid soils, the temperature and rainfall regime may exert over-control over microbial communities (Stres et al., 2008; Pasternak et al., 2013; Nielsen and Ball, 2015; Armstrong et al., 2016; Zhao et al., 2016). Previous studies based on T-RFLP analysis of 16S rRNA gene have shown that temperature and drought stress caused by seasonal differences were the major forces on the modulation of the microbial community in bulk soil and rhizosphere of Cereus jamacaru (cactus) and two leguminous trees (Mimosa tenuiflora and Piptadenia stipulacea) from the Caatinga biome (Lançoni et al., 2013; Nessner Kavamura et al., 2013). Herein, approximately 400 million (about 40 Gb) high-quality reads were obtained after control quality, with an average length of 100 bp (Supplementary Table S3). Taxonomic classification of reads based on the lowest common ancestor (LCA) algorithm showed that the Bacteria dominated CEF and BEF soil samples in the dry (94.45 and 94.11 %, respectively) and rainy (94.28 and 94.75 %, respectively) seasons. The remaining reads corresponded to Archaea and Eukarya, Viruses and unassigned sequences (Supplementary Figure S2a). Principal component analysis (PCA) based on the taxonomic profiles (at genus level) showed a clear separation between soil microbial communities of preserved (CEF) and agriculture-affected (BEF) sites as well as between the dry and rainy seasons of both experimental fields (Figure 3C). Canonical correspondence analysis (CCA) was used to determine whether correlations in microbial structure were associated with soil chemical parameters. The preserved CEF soils of rainy season tended to cluster with soil humidity; whereas soil micro- and macronutrients content (N, Cu, K, CTC, SB, P, Fe, Mg, Mn, B, Ca, and O.M) contributed most to the variance in BEF soils of the dry season (Figure 4).
Figure 3. Relative abundances of major bacterial taxonomic groups at the phylum level in Caatinga and Bebedouro experimental fields in dry and rainy seasons (A) and (B) Microbial phyla that exhibited significant differences (p < 0.05) between dry and rainy seasons of CEF and BEF soils. Significance was determined by Fisher’s exact test with Story’s FDR correction for multiple comparisons. Only taxa with difference between proportions > 0.2 (i.e., considered large effect) are shown (C) Principal Component Analysis (PCA) based taxonomic profile at genus level.
Figure 4. Canonical correspondence analysis (CCA) based on microbial community composition (at phylum level) and soil properties of the soils collected from Caatinga and Bebedouro experimental fields.
These results suggested that community structure in preserve soils is shaped primarily by the natural seasonality factors, whereas nutrient supply seems to be a strongest factor outline the microbial structure in fertilization-affected soils of the dry season.
Enrichment of Bacterial Groups According to Seasonal Variation and Agriculture Impact
At the phylum level, bacterial composition changed according to the experimental field and seasonal variation (Figure 3A). A significant enrichment (p < 0.05) of the phyla Actinobacteria and Cyanobacteria was found in the dry season of CEF soils in contrast to an increase of Proteobacteria, Acidobacteria and Bacteroidetes in the rainy season (Figure 3B). Similar responses of have been described in bacterial communities of soils collected from other sites of Caatinga biome (Lançoni et al., 2013; Nessner Kavamura et al., 2013), as well as in many drylands subjected to rainfall regimes (Barnard et al., 2013; Zhao et al., 2016). Actinobacteria, Proteobacteria and Acidobacteria were also the most responsive phyla in African and Japanese forest soils exposed to drying and rewetting treatments (Zhou et al., 2016). The increase of Actinobacteria in the drought periods is perhaps unsurprising, given their known ability to grow under high temperature, salt concentrations and radiation as well as wide metabolic capacity, which also makes them interesting targets for bioprospecting purposes (Mohammadipanah and Wink, 2016). Similarly, Cyanobacteria members dwelling in arid environments have developed a broad range of strategies for ecological adaptation, including synthesis of extracellular polysaccharides for water retention (Tamaru et al., 2005), protective shield against UV radiation by pigments (Gao and Ye, 2007), and also N2-fixation and photosynthesis for nitrogen (Tashyreva and Elster, 2015) and carbon (Harel et al., 2004) supply, respectively. Indeed, the low moisture required to be photosynthetically active supports the increase in the drought period and the hypothesis of a vital role of this photoautotrophic group as primary producers in dryland environments (Harel et al., 2004; Tracy et al., 2010; Lins et al., 2016). On the other hand, Proteobacteria and Bacteroidetes have copiotrophic strategies with quick growth responses to high moisture and resource availability (Fierer et al., 2007).
Interestingly, some similar and differential responses were observed in soils influenced by ferti-irrigation crop area. The phylum Proteobacteria was significantly overrepresented (p < 0.05) in BEF soil from the rainy season, while Acidobacteria, Firmicutes and Chloroflexi were more enriched in the dry season (Figure 3B). The non-restrictive water content (Table 1) during the dry season ensured by the irrigation practices may explain the lack of statistical difference (p < 0.05) observed in the Cyanobacteria and Actinobacteria abundances among BEF soils. Interesting, Acidobacteria appeared notably more abundant in the dry period of (Figures 3A,B). Given the high nutrient availability in the dry period of BEF soils, the notable enrichment of Acidobacteria disagree with the described oligotrophic life-strategy (analog to K-strategist: slow-growing) of this phylum according to nutrition and growth based-classification (Fierer et al., 2007). However, this scenario is consistent with the different ecological classification proposed more recently by Barnard et al. (2013), where Acidobacteria subgroups displayed a fast life-strategy in response to greater water and nutrient availability in California grasslands soils.
Even with the large proportion of unclassified reads limited by the direct annotation of short reads (Supplementary Table S3), LCA-algorithm was used for phylogenetic resolution at the genus level and provides insights into specific metabolic potential and life-strategy of enriched bacterial groups. Characterized as strictly aerobic, the Conexibacter genus (belonging to Actinobacteria Phylum) was the most represented in all soil metagenomes (Supplementary Figure S2b), except in the irrigation-affected soils (BEF) that displayed a greater abundance of ‘Candidatus Solibacter’ (Supplementary Figure S2b). Although the ecological role of Conexibacter is still not understood, previous metagenomic studies have already reported these microbes as dominant in Caatinga soils (Pacchioni et al., 2014) and quite abundant in other soils (Janssen, 2006; Seki et al., 2012; Deng et al., 2015). Moreover, soil Conexibacter strains were reported as slow-growing microorganisms and with high genomic G+C-content (over 70%), key attributes for survival under stressed arid conditions (Mann and Chen, 2010; Seki et al., 2012).
Interestingly, with the exception of Mycobacterium (enriched in the rainy season), all genera belonging to the phylum Actinobacteria (Conexibacter, Rubrobacter, Streptomyces, Frankia, Thermobifida, Thermomonospora, Nocardioides and Acidothermus) were highly enriched in the drought period of preserved CEF soils (Figure 5). Most of them are described as polyextremotolerant bacteria remarkably adapted to thrive under harsh conditions (Shivlata and Satyanarayana, 2015). These findings corroborate with previous reports on Actinobacteria members variation in response to irrigation influence, suggesting that this bacterial group may be a good marker for “health” of arid soils (Frenk et al., 2014; Wafula et al., 2015). On the other hand, taxa affiliated to Proteobacteria (Bradyrhizobium, Rhodopseudomonas, Burkholderia, Nitrobacter and Pseudomonas) and Acidobacteria (‘Candidatus Koribacter’ and ‘Candidatus Solibacter’) were slightly overrepresented to the rainy season (Figure 5).
Figure 5. Significant differences in microbial communities at the genus taxonomic levels between dry and rainy seasons of CEF and BEF. Significance was determined by Fisher’s exact test with Story’s FDR correction for multiple comparisons, with corrected p < 0.05 significance filter. Only taxa with difference between proportions > 0.2 (i.e., considered large effect) are shown.
Unlike, BEF soils of the dry season impacted by ferti-irrigation practices showed strong enrichment (p < 0.05) of ‘Candidatus Koribacter’ and ‘Candidatus Solibacter’ (belonging to the Acidobacteria phylum) (Figure 5). Although widely abundant and diverse across soils (Janssen, 2006; Fierer et al., 2007; Navarrete et al., 2013). the physiology and ecological roles of the Acidobacteria phylum are poorly understood due to the difficulty to culture representative isolates (Kielak et al., 2016). Despite they were described as oligotrophic and K-strategist lifestyle (Fierer et al., 2007), integrated data derived from physiological, genomic and metagenomic studies gathered from Kielak et al. (2016) have pointed out that different Acidobacteria subdivisions may thrive in non-oligotrophic niches. Likewise, a positive correlation of Acidobacteria subgroups and carbon availability was showed in agricultural soils from the Amazon rainforest (Navarrete et al., 2013) and ultramafic soils from a tropical savanna (Pessoa-Filho et al., 2015). Wafula et al. (2015) also reported a positive correlation of total N, organic matter and moisture with Acidobacteria communities from treated wastewater-irrigated soils. In fact, bacterial isolates of this phylum are able to grow in high concentrations of C sources (Eichorst et al., 2007; Castro et al., 2014; García-Fraile et al., 2016), and genetic evidence have revealed the presence of genes required to decompose complex substrates (García-Fraile et al., 2016; Kielak et al., 2016; Lacerda Júnior et al., 2017). Genomic insights also presume a role of ‘Candidatus Koribacter’ and ‘Candidatus Solibacter’ in N2 cycling by the reduction of nitrate, nitrite, and possibly nitric oxide (Ward et al., 2009).
All these findings support and shed new insights on the role of the Acidobacteria in the nutrient and carbon cycling in agriculture-impacted soils with high N and organic matter content.
On the other hand, soil community of rainy season was most enriched by Conexibacter, Rubrobacter, Salinispora and Mycobacterium (Actinobacteria Phylum), and other genera belonging to Proteobacteria (Bradyrhizobium, Rhodopseudomonas, Methylobacterium and Nitrobacter) (Figure 5). It is interesting to note that most of the Proteobacteria-affiliated taxa (especially Alphaproteobacteria) enriched in rainy soils from both experimental fields have generally been characterized as fast-growing and free-living N-fixing and nitrifying bacteria in soil ecosystems (Tsoy et al., 2016). The abundance of Bradyrhizobium in both sampling seasons and soils was remarkable. This genus is widely spread in the soils worldwide (Delgado-Baquerizo et al., 2018), including the bulk soil in Brazilian semi-arid lands (Pacchioni et al., 2014). Bradyrhizobium also nodulates several crops and native legumes of Brazilian semi-arid (Martins et al., 2003; Leite et al., 2017; Marinho et al., 2017; Menezes et al., 2017; Santos et al., 2017; Rodrigues et al., 2018). The abundance of Bradyrhizobium in the soil evaluated, even in low water availability, indicates their drought tolerance and biotechnological potential for rhizobial prospection studies.
Anaeromyxobacter, a facultative anaerobic myxobacterium, was the only Proteobacteria-affiliated genus significantly overrepresented (p < 0.05, Figure 5) in the soils with land-use influence. Anaeromyxobacter spp. have been proposed to contribute to in situ bioremediation applications for ferric-iron and uranium reduction, herbicides degradation, and are important in N and C cycling in aquatic, sedimentary, and agricultural soils (He and Sanford, 2003; North et al., 2004; Thomas et al., 2009; Sanford et al., 2012).
The increase of Anaeromyxobacter spp. and Acidobacteria groups in agriculture-affected soils may be explained by the process called “priming effect,” generated by the nutrient input. Furthermore, their previously recognized metabolic and genetic potential presume a role in the nutrient cycling (mainly N), mitigating the impacted of the fertilization, such as soil acidification, nitrate leaching, and loss of biodiversity (Ward et al., 2009; Sanford et al., 2012).
Differences in the Functional Traits of the Microbiome From Preserved Soils Under the Seasonal Regime
The functional composition of Caatinga soil metagenomes was analyzed using high-quality reads obtained after quality control filtering of MG-RAST pipeline. Only 20–25% of total sequences could be assigned to functional categories using the SEED database (Supplementary Table S3). Most of them were assigned to one of the following level 1 SEED subsystems related to housekeeping functions: “carbohydrates,” “protein metabolism,” “amino acids and derivatives” and “RNA and DNA metabolism.” The relatively minor categories represented specific ecological functions such as nutrient cycling, stress response and dormancy/sporulation strategies (Supplementary Figure S3). STAMP software was used in the statistical comparisons at a finer functional classification (SEED level 2 and 3) in order to the understanding of how the functional capabilities encoded by the soil microbial communities are affected in response to seasonal and anthropogenic factors.
The most striking functional difference between microbial communities from CEF soils of was detected in Virulence, Disease and Defense, followed by Carbohydrates subsystem (SEED level 1, Supplementary Figure S3). The analysis of the functions contained in these subsystems indicated that toxic compounds and antibiotic resistance genes, as well as others likely associated with microbe-microbe competition, were prevalent in moist soils. Interestingly, functions within “Carbohydrates” subsystem overrepresented in dry soils were related with the metabolism of labile C sources, such as monosaccharides (i.e., L-rhamnose, D-ribose, L-arabinose) and oligosaccharides (i.e., maltose/maltodextrin); or osmoprotective sugars (i.e., mannitol and inositol), which also seem to play a role in the drought tolerance (Figure 6 and Supplementary Figure S4; Zahid and Deppenmeier, 2016). On the other hand, moist soils showed an increase of polysaccharide degradation pathways, including xyloglucans (plant polysaccharides) and glycogen. In fact, a vast repertoire of lignocellulolytic enzymes was mined by sequence-based screening from a metagenomic library of a Caatinga soil in the rainy season (Lacerda Júnior et al., 2017). This is expected since the litterfall is the major organic source deposited in the soil of Caatinga forest (Salgado et al., 2015). Also, genes related to “trehalose uptake and utilization” and “trehalose biosynthesis” were more abundant in dry and rainy season soils, respectively (Supplementary Figure S4). Trehalose is an important compatible solute accumulated in bacterial cells as a response to hyperosmotic conditions and as an energy source for long periods of carbon starvation (Argüelles, 2000; Ruhal et al., 2013). The concomitant role of trehalose in drought tolerance and thermoprotection has been established in the soil bacterium Rhizobium etli under free-living conditions (Reina-Bueno et al., 2012). It protects and stabilizes membranes and proteins against numerous kinds of stress, including dehydration, oxidation and heat conditions (Crowe, 2007). Thus, this contrasting ability of the microbiome for acquisition and utilization or for biosynthesis of trehalose may be related to the different requirement for osmoadaptation in dry and wet soils.
Figure 6. Fisher’s exact test-based comparative analysis of selected metabolic pathways of SEED subsystems (levels 2 and 3) enriched in dry and rainy seasons of CEF and BEF. The Story’s FDR correction was used at P < 0.05 filter.
Additionally, we also explored some of the specific functional categories related to nutrient cycling. Although a subtle statistical difference (p < 0.05) was observed in the SEED level 1 “Nitrogen Metabolism” subsystems (Supplementary Figure S3), selective preferences in the biological processes used to uptake and/or convert the available nutrients were showed by lower hierarchical gene annotation. Level 3 of the SEED analysis showed a significant increase in the “nitrogen fixation” and “ammonia assimilation” pathways in moist soils in comparison to a higher abundance of genes involved in “nitrite and nitrate ammonification” (direct uptake of nitrate/nitrite and subsequent reduction to ammonia), and “allantoin utilization” in dry soils (Figure 6, p < 0.05). Slight differences were found in the “Sulfur Metabolism” subsystem (Level 1) (Supplementary Figure S3). However, lower level SEED predictions (level 2–3) revealed a higher abundance of genes related to “inorganic sulfur assimilation” in the dry soils in contrast to the enrichment of genes related to “organic sulfur compounds uptake” (especially glutathione as a sulfur source) in the wet season soils (Figure 6). Likewise, potassium (K+) metabolism was more abundant in moist soils with the majority of genes involved in K+ homeostasis, contrasting with the greatest amount of hyperosmotic K+ uptake genes in the dry soils (Figure 6).
The capacity of response to osmotic stress by adaptive mechanisms is an important characteristic that enables drought-tolerant bacteria to survive and proliferate in soils with water shortage (Lebre et al., 2017). As expected, an increase in the functional categories associated with “dormancy/sporulation strategies” and “stress response” was also found in the preserved dry soils (Supplementary Figure S3), with overrepresented sequences involved in the synthesis and uptake systems of compatible solutes including choline and betaine, or encoding putative ion transporters. Comparative analyses using Level 3 “stress response subsystem” (Figure 6) also showed enrichment of several genes encoding bacterial sigma transcription factors responsible for the regulation of bacterial gene expression induced in response to a variety of environmental stresses/physiological signals in the dry soils (Feklístov et al., 2014; Paget, 2015). In contrast, aquaporin Z genes (AqpZ), encoding membrane water channel proteins presumed to play an important role in osmoregulation by accelerating transmembrane water flux in response to abrupt changes in osmotic pressure, were more abundant in rainy season soils. Other general categories, such as “Cell Wall and Capsule,” “membrane transport,” “protein metabolism,” “RNA metabolism,” “Respiration,” “Cofactor, vitamins, prosthetic groups and pigments” categories were more abundant in preserved wet soil of CEF (Supplementary Figure S3).
The enrichment of microbial groups and functional categories related to osmoregulation and nutrient and carbohydrate metabolism shows that the microbial community of preserved soils undergoes an adjustment during the rainy and dry cycles, reflecting differences in the genetic potential for nutrient cycling and carbon acquisition in the environment. This is probably a consequence of seasonal differences in water and nutrient availability since leaf fall in the drought periods is an intrinsic feature of the Caatinga forest; and adaptive processes that ensure long-term resilience over the years of rainfall regime.
Influence of Agriculture Practices on the Functional Features of Soil Microbiome
As shown, the constant supply of nutrients and water break down the natural water regime in the soils of BEF during the dry season. Hence, the influence of land-use practices outlining the functional traits of the microbial communities was evaluated. A greater relative abundance of “carbohydrates,” “nitrogen metabolism” and “dormancy and sporulation” subsystems (SEED 1) were found in irrigation-affected soils of the dry season. A divergent trend was observed with the enrichment of functions related to water stress response, such as “biosynthesis” and “uptake and utilization” pathways of trehalose, when compared to the differences observed between CEF soils (Figure 6). It has also been shown the prevalence of genes related to di- and oligosaccharides, and polysaccharides degradation (mainly glycogen metabolism).
Interestingly, the overrepresented functions included in the “Stress Response” subsystems in irrigation-affected dry soils were associated with regulation and genetic responses to oxidative stresses such as “mycothiol biosynthesis,” “protection from reactive oxygen species” and “redox-dependent processes” (Figure 6). The increase of these functions suggests a greater functional requirement for cellular protection in response to damage environmental conditions. On the other hand, soils of the rainy season showed a greater enrichment of functions related to osmotic stress, such as “choline and betaine uptake” and “betaine biosynthesis,” and “synthesis of osmoregulated periplasmic glucans” (Figure 6). This unexpected enrichment of functions associated to resistance to osmotic stress in the rainy soils may be the reflection of an adaptation to rewetting-desiccation cycles due to sparse rainfall, in opposition to the continuous water supply by irrigation practices in the soils of the dry season.
Key processes related to the terrestrial N cycle are catalyzed by reductive and oxidative enzymes, in which microbes have a predominant role (Tsoy et al., 2016). Despite the increase of reads assigned to “Nitrogen Metabolism” subsystem (SEED Level 1) in the irrigation-affectedsoils (Supplementary Figure S3), no significant difference (p > 0.05) was found in the functional subcategory of “nitrogen fixation” (Figure 6). However, genes encoding enzymes that play roles in protecting from nitrosative stress caused by overproduction of reactive nitric oxide (NO) and NO-derived compounds, nitrate and nitrite ammonification and denitrification processes were overrepresented (Figure 6). The application of N fertilizers to soils has been shown to stimulate nitrogen loss via gaseous emission derived from the denitrification process (Wolsing and Priemé, 2004; Wallenstein et al., 2006; Attard et al., 2011; Yin et al., 2015). The influence of high nutrient input and continuous supply of water coupled with the high abundance of denitrifying bacteria in BEF soils of the dry season, provide highly favorable conditions for denitrification, as shown by Rasche et al. (2011) in temperate forest soils.
The analysis of potassium metabolism at (SEED level-3) revealed a predominance of genes related to K homeostasis in BEF soils of the dry season (Figure 6), probably caused by the high K concentration in such impacted soils (Table 1). Concerning the sulfur metabolism subsystem, organic alkanesulfonate assimilation was the most prevalent route for the environmental sulfur acquisition in the BEF soils of the rainy season. In contrast to a greater abundance of inorganic sulfur assimilation-related genes and other organosulfur reduction and oxidation processes in BEF soils of the dry season (Figure 6). Although soil sulfur content has not been measured, these results indicated differential microbial community strategies to acquire sulfur in the environment in response to impacted conditions.
Conclusion
Overall, the metagenomics survey showed that natural seasonality and land-use practices differentially affect the composition and functional traits of the soil microbiome from preserved and agriculture-influenced experimental fields of the Caatinga dry forest. Multivariate analysis suggested that environmental factors associated with seasonal cycles (mainly soil moisture) and ferti-irrigation practices (micro- and macronutrients availability) were the major environmental factors driving the microbial community of preserved and agriculture-affected soils, respectively. The strong enrichment of ‘Candidatus Koribacter’ and ‘Candidatus Solibacter’ raises novel questions about the role of the poorly-known Acidobacteria phylum in fertilization-affected soils and opens the possibility to explore such taxa as bioindicators for agricultural soil management effects in this wide semiarid area. Cellular processes related to osmoregulation and carbon and nutrient cycling are the most impacted functions in preserved soils, indicating different strategies for nutrient cycling and carbon acquisition of the microbiome under different rainfall patterns. However, the nutrient input and disruption of natural water regime by ferti-irrigation practices may impact soil biogeochemical processes linked to microbial processes, such as nitrosative stress and denitrification-related functions, with potential impacts in the soil functioning. This work is an important step toward the understanding of how environment factors affect the soil microbiome in remaining fragments of a seasonal tropical dry forest, occurring exclusively in Brazil. The next step is to understand if these long-term impacts influence the capability of the soil microbial community to maintain functional resilience to perform their normal nutrient cycling and soil structure maintenance processes across the years. These keys ecological drivers may help for the development of sustainable farming systems in dryland ecosystem.
Author Contributions
GL-J and PF-J carried out the experimental design and soil collection in the experimental fields. GL-J and SdS performed the extraction and quantification of metagenomic DNA. MN performed the bioinformatic analyses of metagenomic data. GL-J, LC, TD, and PF-J performed statistical analyses. The manuscript was written by GL-J (primarily) and PF-J. IM and VO read and refined the final manuscript.
Funding
We are grateful to São Paulo Research Foundation (FAPESP, Process No. 2013/09386-9) for providing the scholarship to GL-J. PF-J and VO thank, respectively, to Brazilian Agricultural Research Corporation (Process 03.13.07.008.00.00) and Brazilian National Council for Scientific and Technological Development – CNPq (Process 308955/2016-1) for financial support.
Conflict of Interest Statement
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.
Acknowledgments
We would like to acknowledge Dr. Magna Soelma Beserra de Moura (Agrometeorology Laboratory, Embrapa Semiárido, Brazil) for providing agrometeorological data. Special thanks to João L. Silva for support on expeditions to Caatinga, and Embrapa for the partial financial support.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.00648/full#supplementary-material
References
Andrew, D. R., Fitak, R. R., Munguia-Vega, A., Racolta, A., Martinson, V. G., and Dontsova, K. (2012). Abiotic factors shape microbial diversity in Sonoran desert soils. Appl. Environ. Microbiol. 78, 7527–7537. doi: 10.1128/AEM.01459-12
Argüelles, J. C. (2000). Physiological roles of trehalose in bacteria and yeasts: a comparative analysis. Arch. Microbiol. 174, 217–224. doi: 10.1007/s002030000192
Armstrong, A., Valverde, A., Ramond, J. B., Makhalanyane, T. P., Jansson, J. K., Hopkins, D. W., et al. (2016). Temporal dynamics of hot desert microbial communities reveal structural and functional responses to water input. Sci. Rep. 6:34434. doi: 10.1038/srep34434
Attard, E., Recous, S., Chabbi, A., De Berranger, C., Guillaumaud, N., Labreuche, J., et al. (2011). Soil environmental conditions rather than denitrifier abundance and diversity drive potential denitrification after changes in land uses. Glob. Chang. Biol. 17, 1975–1989. doi: 10.1111/j.1365-2486.2010.02340.x
Bardgett, R. D., and Van Der Putten, W. H. (2014). Belowground biodiversity and ecosystem functioning. Nature 515, 505–511. doi: 10.1038/nature13855
Barnard, R. L., Osborne, C. A., and Firestone, M. K. (2013). Responses of soil bacterial and fungal communities to extreme desiccation and rewetting. ISME J. 7, 2229–2241. doi: 10.1038/ismej.2013.104
Beuchle, R., Grecchi, R. C., Shimabukuro, Y. E., Seliger, R., Eva, H. D., Sano, E., et al. (2015). Land cover changes in the Brazilian Cerrado and Caatinga biomes from 1990 to 2010 based on a systematic remote sensing sampling approach. Appl. Geogr. 58, 116–127. doi: 10.1016/j.apgeog.2015.01.017
Castro, W. D. O., Torres-Ballesteros, A. M., Nakayama, R., Melo, S., Pellizari, V. H., Silva, A., et al. (2014). Draft genome sequence of Haloferax sp. Strain ATB1, isolated from a semi-arid region in the Brazilian Caatinga. Genome Announc. 2, 4–5. doi: 10.1128/genomeA.00812-14
Crowe, J. H. (2007). Trehalose as a “chemical chaperone”: Fact and fantasy. Adv. Exp. Med. Biol. 594, 143–158. doi: 10.1007/978-0-387-39975-1_13
da Silva, A. F., de Freitas, A. D. S., Costa, T. L., Fernandes-Júnior, P. I., Martins, L. M. V., Santos, C. E. R. S., et al. (2017). Biological nitrogen fixation in tropical dry forests with different legume diversity and abundance. Nutr. Cycl. Agroecosystems 107, 321–334. doi: 10.1007/s10705-017-9834-1
de Albuquerque, U. P., de Lima Araújo, E., El-Deir, A. C. A., de Lima, A. L. A., Souto, A., Bezerra, B. M., et al. (2012). Caatinga revisited: ecology and conservation of an important seasonal dry forest. Sci. World J. 2012, 1–18. doi: 10.1100/2012/205182
Delgado-Baquerizo, M., Oliverio, A. M., Brewer, T. E., Benavent-González, A., Eldridge, D. J., Bardgett, R. D., et al. (2018). A global atlas of the dominant bacteria found in soil. Science 359, 320–325. doi: 10.1126/science.aap9516
Deng, J., Gu, Y., Zhang, J., Xue, K., Qin, Y., Yuan, M., et al. (2015). Shifts of tundra bacterial and archaeal communities along a permafrost thaw gradient in Alaska. Mol. Ecol. 24, 222–234. doi: 10.1111/mec.13015
Ding, G. C., Piceno, Y. M., Heuer, H., Weinert, N., Dohrmann, A. B., Carrillo, A., et al. (2013). Changes of soil bacterial diversity as a consequence of agricultural land use in a semi-arid ecosystem. PLoS One 8:e59497. doi: 10.1371/journal.pone.0059497
Drenovsky, R. E., Vo, D., Graham, K. J., and Scow, K. M. (2004). Soil water content and organic carbon availability are major determinants of soil microbial community composition. Microb. Ecol. 48, 424–430. doi: 10.1007/s00248-003-1063-2
Eichorst, S. A., Breznak, J. A., and Schmidt, T. M. (2007). Isolation and characterization of soil bacteria that define Terriglobus gen. nov., in the phylum Acidobacteria. Appl. Environ. Microbiol. 73, 2708–2717. doi: 10.1128/AEM.02140-06
Feklístov, A., Sharon, B. D., Darst, S. A., and Gross, C. A. (2014). Bacterial sigma factors: a historical, structural, and genomic perspective. Annu. Rev. Microbiol. 68, 357–376. doi: 10.1146/annurev-micro-092412-155737
Fernandes-Júnior, P. I., Aidar, S., de, T., Morgante, C. V., Gava, C. A. T., Zilli, J. É, et al. (2015). The resurrection plant Tripogon spicatus (Poaceae) harbors a diversity of plant growth promoting bacteria in northeastern Brazilian Caatinga. Rev. Bras. Cienc. do Solo 39, 993–1002. doi: 10.1590/01000683rbcs20140646
Ferreira, D. F. (2011). Sisvar: a computer statistical analysis system. Cienc. Agrotecnol. 35, 1039–1042. doi: 10.1590/S1413-70542011000600001
Fierer, N., Bradford, M. A., and Jackson, R. B. (2007). Toward an ecological classification of soil bacteria. Ecology 88, 1354–1364. doi: 10.1890/05-1839
Frenk, S., Hadar, Y., and Minz, D. (2014). Resilience of soil bacterial community to irrigation with water of different qualities under Mediterranean climate. Environ. Microbiol. 16, 559–569. doi: 10.1111/1462-2920.12183
Gao, K., and Ye, C. (2007). Photosynthetic insensitivity of the terrestrial cyanobacterium Nostoc flagelliforme to solar UV radiation while rehydrated or desiccated. J. Phycol. 43, 628–635. doi: 10.1111/j.1529-8817.2007.00358.x
García-Fraile, P., Benada, O., Cajthaml, T., Baldrian, P., and Lladó, S. (2016). Terracidiphilus gabretensis gen. nov., sp. nov., an abundant and active forest soil acidobacterium important in organic matter transformation. Appl. Environ. Microbiol. 82, 560–569. doi: 10.1128/AEM.03353-15
Gillespie, T. W., Lipkin, B., Sullivan, L., Benowitz, D. R., Pau, S., and Keppel, G. (2012). The rarest and least protected forests in biodiversity hotspots. Biodivers. Conserv. 21, 3597–3611. doi: 10.1007/s10531-012-0384-1
Hammer, Ø, Harper, D. A. T., and Ryan, P. D. (2001). PAST: paleontological statistics software package for education and data analysis. Palaeontol. Electron 4, 1–9. doi: 10.1016/j.bcp.2008.05.025
Harel, Y., Ohad, I., and Kaplan, A. (2004). Activation of photosynthesis and resistance to photoinhibition in cyanobacteria within biological desert crust. Plant Physiol. 136, 3070–3079. doi: 10.1104/pp.104.047712
He, Q., and Sanford, R. A. (2003). Characterization of Fe(III) reduction by chlororespiring Anaeromxyobacter dehalogenans. Appl. Environ. Microbiol. 69, 2712–2718. doi: 10.1128/AEM.69.5.2712-2718.2003
Huang, J., Yu, H., Guan, X., Wang, G., and Guo, R. (2016). Accelerated dryland expansion under climate change. Nat. Clim. Chang. 6, 166–171. doi: 10.1038/nclimate2837
Hughes, J. B., Hellmann, J. J., Ricketts, T. H., and Bohannan, B. J. M. (2001). Counting the uncountable: statistical approaches to estimating microbial diversity. Appl. Environ. Microbiol. 67, 4399–4406. doi: 10.1128/AEM.67.10.4399-4406.2001
Janssen, P. H. (2006). Identifying the dominant soil bacterial taxa in libraries of 16S rRNA and 16S rRNA genes. Appl. Environ. Microbiol. 72, 1719–1728. doi: 10.1128/AEM.72.3.1719-1728.2006
Keeney, D., and Nelson, D. (1982). “Nitrogen-Inorganic forms,” in Methods of Soil Analysis part 2: Chemical and Microbiological Properties, 2nd Edn, eds A. L. Page, R. H. Miller, and D. R. Keeney (Madison, WI: American society of Agronomy), 643–698.
Kielak, A. M., Barreto, C. C., Kowalchuk, G. A., van Veen, J. A., and Kuramae, E. E. (2016). The ecology of Acidobacteria: moving beyond genes and genomes. Front. Microbiol. 7:744. doi: 10.3389/fmicb.2016.00744
Lacerda Júnior, G. V., Noronha, M. F., de Sousa, S. T. P., Cabral, L., Domingos, D. F., Sáber, M. L., et al. (2017). Potential of semiarid soil from Caatinga biome as a novel source for mining lignocellulose-degrading enzymes. FEMS Microbiol. Ecol. 93:fiw248. doi: 10.1093/femsec/fiw248
Lançoni, M. D., Taketani, R. G., Kavamura, V. N., and de Melo, I. S. (2013). Microbial community biogeographic patterns in the rhizosphere of two Brazilian semi-arid leguminous trees. World J. Microbiol. Biotechnol. 29, 1233–1241. doi: 10.1007/s11274-013-1286-4
Lauber, C. L., Ramirez, K. S., Aanderud, Z., Lennon, J., and Fierer, N. (2013). Temporal variability in soil microbial communities across land-use types. ISME J. 7, 1641–1650. doi: 10.1038/ismej.2013.50
Lebre, P. H., De Maayer, P., and Cowan, D. A. (2017). Xerotolerant bacteria: surviving through a dry spell. Nat. Rev. Microbiol. 15, 285–296. doi: 10.1038/nrmicro.2017.16
Leite, J., Fischer, D., Rouws, L. F. M., Fernandes-Júnior, P. I., Hofmann, A., Kublik, S., et al. (2017). Cowpea nodules harbor non-rhizobial bacterial communities that are shaped by soil type rather than plant genotype. Front. Plant Sci. 7:2064. doi: 10.3389/fpls.2016.02064
Liao, H., Chapman, S. J., Li, Y., and Yao, H. (2018). Dynamics of microbial biomass and community composition after short-term water status change in Chinese paddy soils, 2018, Environmental Science and Pollution Research, 25: 2932-2941). Environ. Sci. Pollut. Res. Int. 25, 2932–2941. doi: 10.1007/s11356-017-0690-y
Lins, R. P. M., Barbosa, L. G., Minillo, A., and De Ceballos, B. S. O. (2016). Cyanobacteria in a eutrophicated reservoir in a semi-arid region in Brazil: dominance and microcystin events of blooms. Rev. Bras. Bot. 39, 583–591. doi: 10.1007/s40415-016-0267-x
Lüneberg, K., Schneider, D., Siebe, C., and Daniel, R. (2018). Drylands soil bacterial community is affected by land use change and different irrigation practices in the Mezquital Valley, Mexico. Sci. Rep. 8:1413. doi: 10.1038/s41598-018-19743-x
Maestre, F. T., Eldridge, D. J., Soliveres, S., Kéfi, S., Delgado-Baquerizo, M., Bowker, M. A., et al. (2016). Structure and functioning of dryland ecosystems in a changing world. Annu. Rev. Ecol. Evol. Syst. 47, 215–237. doi: 10.1146/annurev-ecolsys-121415-032311
Mann, S., and Chen, Y. P. P. (2010). Bacterial genomic G + C composition-eliciting environmental adaptation. Genomics 95, 7–15. doi: 10.1016/j.ygeno.2009.09.002
Marcon, A. E., Ferreira, D., de, M., Moura, M., de, F., Campos, T. F., et al. (2010). Genotoxic analysis in aquatic environment under influence of cyanobacteria, metal and radioactivity. Chemosphere 81, 773–780. doi: 10.1016/j.chemosphere.2010.07.006
Marinho, R., de, C. N., Ferreira, L., de, V. M., da Silva, A. F., Martins, L. M. V., et al. (2017). Symbiotic and agronomic efficiency of new cowpea rhizobia from Brazilian semi-arid. Bragantia 76, 273–281. doi: 10.1590/1678-4499.003
Martins, L. M. V., Xavier, G. R., Range, F. W., Ribeiro, J. R. A., Neves, M. C. P., Morgado, L. B., et al. (2003). Contribution of biological nitrogen fixation to cowpea: a strategy for improving grain yield in the semi-arid region of Brazil. Biol. Fertil. Soils 38, 333–339. doi: 10.1007/s00374-003-0668-4
Menezes, K. A. S., Escobar, I. E. C., Fraiz, A. C. R., Martins, L. M. V., and Fernandes-Júnior, P. I. (2017). Genetic variability and symbiotic efficiency of Erythrina velutina willd. Root nodule bacteria from the semi-arid region in Northeastern Brazil. Rev. Bras. Cienc. Do Solo 41:e0160302. doi: 10.1590/18069657rbcs20160302
Menezes, R., Sampaio, E., Giongo, V., and Pérez-Marin, A. (2012). Biogeochemical cycling in terrestrial ecosystems of the Caatinga Biome. Braz. J. Biol. 72, 643–653. doi: 10.1590/S1519-69842012000400004
Meyer, F., Paarmann, D., D’Souza, M., Olson, R., Glass, E. M., Kubal, M., et al. (2008). The metagenomics RAST server - a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinformatics 9:386. doi: 10.1186/1471-2105-9-386
Mohammadipanah, F., and Wink, J. (2016). Actinobacteria from arid and desert habitats: diversity and biological activity. Front. Microbiol. 6:1541. doi: 10.3389/fmicb.2015.01541
Monteiro, J. M., Vollú, R. E., Coelho, M. R. R., Alviano, C. S., Blank, A. F., and Seldin, L. (2009). Comparison of the bacterial community and characterization of plant growth-promoting rhizobacteria from different genotypes of Chrysopogon zizanioides (L.) Roberty (vetiver) rhizospheres. J. Microbiol. 47, 363–370. doi: 10.1007/s12275-009-0048-3
Murphy, P. G., and Lugo, A. E. (1986). Ecology of tropical dry forest. Annu. Rev. Ecol. Syst. 17, 67–88. doi: 10.1146/annurev.es.17.110186.000435
Nacke, H., Fischer, C., Thürmer, A., Meinicke, P., and Daniel, R. (2014). Land use type significantly affects microbial gene transcription in soil. Microb. Ecol. 67, 919–930. doi: 10.1007/s00248-014-0377-6
Navarrete, A. A., Kuramae, E. E., de Hollander, M., Pijl, A. S., van Veen, J. A., and Tsai, S. M. (2013). Acidobacterial community responses to agricultural management of soybean in Amazon forest soils. FEMS Microbiol. Ecol. 83, 607–621. doi: 10.1111/1574-6941.12018
Nessner Kavamura, V., Taketani, R. G., Lançoni, M. D., Andreote, F. D., Mendes, R., Soares, et al. (2013). Water Regime Influences Bulk Soil and Rhizosphere of Cereus jamacaru Bacterial Communities in the Brazilian Caatinga Biome. PLoS One 8:e73606. doi: 10.1371/journal.pone.0073606
Nielsen, U. N., and Ball, B. A. (2015). Impacts of altered precipitation regimes on soil communities and biogeochemistry in arid and semi-arid ecosystems. Glob. Chang. Biol. 21, 1407–1421. doi: 10.1111/gcb.12789
North, N. N., Dollhopf, S. L., Petrie, L., Istok, J. D., Balkwill, D. L., and Kostka, J. E. (2004). Change in bacterial community structure during in situ biostimulation of subsurface sediment cocontaminated with uranium and nitrate. Appl. Environ. Microbiol. 70, 4911–4920. doi: 10.1128/AEM.70.8.4911-4920.2004
Overbeek, R., Begley, T., Butler, R. M., Choudhuri, J. V., Chuang, H. Y., Cohoon, M., et al. (2005). The subsystems approach to genome annotation and its use in the project to annotate 1000 genomes. Nucl. Acids Res. 33, 5691–5702. doi: 10.1093/nar/gki866
Pacca, E., Mattar, L., Teixeira, T., Barros, V., Cunha, B. B., Souza, J. F., et al. (2018). Federal conservation units in Brazil: the situation of biomes and regions. Floresta Ambiente 25:e20150051. doi: 10.1590/2179-8087.005115
Pacchioni, R. G., Carvalho, F. M., Thompson, C. E., Faustino, A. L. F., Nicolini, F., Pereira, T. S., et al. (2014). Taxonomic and functional profiles of soil samples from Atlantic forest and Caatinga biomes in northeastern Brazil. Microbiologyopen 3, 299–315. doi: 10.1002/mbo3.169
Paget, M. S. (2015). Bacterial sigma factors and anti-sigma factors: structure, function and distribution. Biomolecules 5, 1245–1265. doi: 10.3390/biom5031245
Pajares, S., Campo, J., Bohannan, B. J. M., and Etchevers, J. D. (2018). Environmental controls on soil microbial communities in a seasonally dry tropical forest. Appl. Environ. Microbiol. 17:e00342-18. doi: 10.1128/AEM.00342-18
Parks, D. H., and Beiko, R. G. (2010). Identifying biologically relevant differences between metagenomic communities. Bioinformatics 26, 715–721. doi: 10.1093/bioinformatics/btq041
Pasternak, Z., Al-Ashhab, A., Gatica, J., Gafny, R., Avraham, S., Minz, D., et al. (2013). Spatial and temporal biogeography of soil microbial communities in arid and semiarid regions. PLoS One 8:e69705. doi: 10.1371/journal.pone.0069705
Pessoa-Filho, M., Barreto, C. C., dos Reis, Junior, F. B., Fragoso, R. R., Costa, F. S., et al. (2015). Microbiological functioning, diversity, and structure of bacterial communities in ultramafic soils from a tropical savanna. Antonie Van Leeuwenhoek 107, 935–949. doi: 10.1007/s10482-015-0386-6
Raij, B., van Andrade, J. C., Cantarella, H., and Quaggio, J. A. (2001). Análise Química Para Avaliação da Fertilidade de Solos Tropicais. Campinas: Instituto Agronômico.
Rasche, F., Knapp, D., Kaiser, C., Koranda, M., Kitzler, B., Zechmeister-Boltenstern, S., et al. (2011). Seasonality and resource availability control bacterial and archaeal communities in soils of a temperate beech forest. ISME J. 5, 389–402. doi: 10.1038/ismej.2010.138
Reina-Bueno, M., Argandoña, M., Nieto, J. J., Hidalgo-García, A., Iglesias-Guerra, F., Delgado, M. J., et al. (2012). Role of trehalose in heat and desiccation tolerance in the soil bacterium Rhizobium etli. BMC Microbiol. 12:207. doi: 10.1186/1471-2180-12-207
Rodrigues, D. R., Silva, A. F., da Cavalcanti, M. I. P., Escobar, I. E. C., Fraiz, A. C. R., Ribeiro, P. R. A., et al. (2018). Phenotypic, genetic and symbiotic characterization of Erythrina velutina rhizobia from Caatinga dry forest. Braz. J. Microbiol. 49, 503–512. doi: 10.1016/j.bjm.2017.09.007
Rodrigues, J. L., Pellizari, V. H., Mueller, R., Baek, K., Jesus, Eda C, Paula, F. S., et al. (2013). Conversion of the Amazon rainforest to agriculture results in biotic homogenization of soil bacterial communities. Proc. Natl. Acad. Sci. U.S.A. 110, 988–993. doi: 10.1073/pnas.1220608110
Ruhal, R., Kataria, R., and Choudhury, B. (2013). Trends in bacterial trehalose metabolism and significant nodes of metabolic pathway in the direction of trehalose accumulation. Microb. Biotechnol. 6, 493–502. doi: 10.1111/1751-7915.12029
Salgado, E. V., Andrade, E. M., de Hevia, J. N., Nunes, E. P., Rodrigues, M. M., and de, A. (2015). Rainfall patterns and the contribution of litter in the caatinga dry tropical forest. Rev. Ciência Agronômica 46, 299–309. doi: 10.5935/1806-6690.20150009
Sanford, R. A., Wagner, D. D., Wu, Q., Chee-Sanford, J. C., Thomas, S. H., Cruz-Garcia, C., et al. (2012). Unexpected nondenitrifier nitrous oxide reductase gene diversity and abundance in soils. Proc. Natl. Acad. Sci. U.S.A. 109, 19709–19714. doi: 10.1073/pnas.1211238109
Santos, J. C., Leal, I. R., Almeida-Cortez, J. S., Fernandes, G. W., and Tabarelli, M. (2011). Caatinga: the scientific negligence experienced by a dry tropical forest. Trop. Conserv. Sci. 4, 276–286. doi: 10.1177/194008291100400306
Santos, J. M. F., dos Casaes Alves, P. A., Silva, V. C., Kruschewsky Rhem, M. F., James, E. K., and Gross, E. (2017). Diverse genotypes of Bradyrhizobium nodulate herbaceous Chamaecrista (Moench) (Fabaceae, Caesalpinioideae) species in Brazil. Syst. Appl. Microbiol. 40, 69–79. doi: 10.1016/j.syapm.2016.12.004
Santos, M. G., Oliveira, M. T., Figueiredo, K. V., Falcão, H. M., Arruda, E. C. P., Almeida-Cortez, J., et al. (2014). Caatinga, the Brazilian dry tropical forest: Can it tolerate climate changes? Theor. Exp. Plant Physiol. 26, 83–99. doi: 10.1007/s40626-014-0008-0
Seki, T., Matsumoto, A., Shimada, R., Inahashi, Y., Omura, S., and Takahashi, Y. (2012). Conexibacter arvalis sp. nov., isolated from a cultivated field soil sample. Int. J. Syst. Evol. Microbiol. 62, 2400–2404. doi: 10.1099/ijs.0.036095-0
Shade, A., Gregory Caporaso, J., Handelsman, J., Knight, R., and Fierer, N. (2013). A meta-analysis of changes in bacterial and archaeal communities with time. ISME J. 7, 1493–1506. doi: 10.1038/ismej.2013.54
Shivlata, L., and Satyanarayana, T. (2015). Thermophilic and alkaliphilic Actinobacteria: biology and potential applications. Front. Microbiol 6:1014. doi: 10.3389/fmicb.2015.01014
Sorensen, L. (2009). A Spatial Analysis Approach to the Global Delineation of Dryland Areas of Relevance to the CBD Programme of Work on Dry and Sub-Humid Lands. Cambridge: UNEP World Conservation Monitoring Centre.
Stres, B., Danevcic, T., Pal, L., Fuka, M. M., Resman, L., Leskovec, S., et al. (2008). Influence of temperature and soil water content on bacterial, archaeal and denitrifying microbial communities in drained fen grassland soil microcosms. FEMS Microbiol. Ecol. 66, 110–122. doi: 10.1111/j.1574-6941.2008.00555.x
Tamaru, Y., Takani, Y., Yoshida, T., and Sakamoto, T. (2005). Crucial role of extracellular polysaccharides in desiccation and freezing tolerance in the terrestrial cyanobacterium Nostoc commune. Appl. Environ. Microbiol. 71, 7327–7333. doi: 10.1128/AEM.71.11.7327-7333.2005
Tashyreva, D., and Elster, J. (2015). Effect of nitrogen starvation on desiccation tolerance of Arctic Microcoleus strains (cyanobacteria). Front. Microbiol. 6:278. doi: 10.3389/fmicb.2015.00278
Thapa, V. R., Ghimire, R., Mikha, M. M., Idowu, O. J., and Marsalis, M. A. (2018). Land Use Effects on soil health in semiarid drylands. Agric. Environ. Lett. 23, 1–5. doi: 10.2134/ael2018.05.0022
Thomas, S. H., Padilla-Crespo, E., Jardine, P. M., Sanford, R. A., and Löffler, F. E. (2009). Diversity and distribution of Anaeromyxobacter strains in a uranium-contaminated subsurface environment with a nonuniform groundwater flow. Appl. Environ. Microbiol. 75, 3679–3687. doi: 10.1128/AEM.02473-08
Tian, Q., Taniguchi, T., Shi, W. Y., Li, G., Yamanaka, N., and Du, S. (2017). Land-use types and soil chemical properties influence soil microbial communities in the semiarid Loess Plateau region in China. Sci. Rep. 28:45289. doi: 10.1038/srep45289
Tracy, C. R., Streten-Joyce, C., Dalton, R., Nussear, K. E., Gibb, K. S., and Christian, K. A. (2010). Microclimate and limits to photosynthesis in a diverse community of hypolithic cyanobacteria in northern Australia. Environ. Microbiol. 12, 592–607. doi: 10.1111/j.1462-2920.2009.02098.x
Tsoy, O. V., Ravcheev, D. A., Èuklina, J., and Gelfand, M. S. (2016). Nitrogen fixation and molecular oxygen: comparative genomic reconstruction of transcription regulation in Alphaproteobacteria. Front. Microbiol. 7:1343. doi: 10.3389/fmicb.2016.01343
United Nations Environmental Management Group (2011). Global Drylands: A UN System-wide Response. Available at: https://www.zaragoza.es/contenidos/medioambiente/onu//issue07/1107-eng.pdf
Wafula, D., White, J. R., Canion, A., Jagoe, C., Pathak, A., and Chauhan, A. (2015). Impacts of long-term irrigation of domestic treated wastewater on soil biogeochemistry and bacterial community structure. Appl. Environ. Microbiol. 81, 7143–7158. doi: 10.1128/AEM.02188-15
Wallenstein, M. D., Peterjohn, W. T., and Schlesinger, W. H. (2006). N fertilization effects on denitrification and N cycling in an aggrading forest. Ecol. Appl. 16, 2168–2176. doi: 10.1890/1051-0761(2006)016[2168:NFEODA]2.0.CO;2
Ward, N. L., Challacombe, J. F., Janssen, P. H., Henrissat, B., Coutinho, P. M., Wu, M., et al. (2009). Three genomes from the phylum Acidobacteria provide insight into the lifestyles of these microorganisms in soils. Appl. Environ. Microbiol. 75, 2046–2056. doi: 10.1128/AEM.02294-08
Wolsing, M., and Priemé, A. (2004). Observation of high seasonal variation in community structure of denitrifying bacteria in arable soil receiving artificial fertilizer and cattle manure by determining T-RFLP of nir gene fragments. FEMS Microbiol. Ecol. 48, 261–271. doi: 10.1016/j.femsec.2004.02.002
Yin, C., Fan, F., Song, A., Cui, P., Li, T., and Liang, Y. (2015). Denitrification potential under different fertilization regimes is closely coupled with changes in the denitrifying community in a black soil. Appl. Microbiol. Biotechnol. 99, 5719–5729. doi: 10.1007/s00253-015-6461-0
Zahid, N., and Deppenmeier, U. (2016). Role of mannitol dehydrogenases in osmoprotection of Gluconobacter oxydans. Appl. Microbiol. Biotechnol. 100, 9967–9978. doi: 10.1007/s00253-016-7680-8
Zhao, C., Miao, Y., Yu, C., Zhu, L., Wang, F., Jiang, L., et al. (2016). Soil microbial community composition and respiration along an experimental precipitation gradient in a semiarid steppe. Sci. Rep. 6:24317. doi: 10.1038/srep24317
Keywords: Caatinga biome, tropical dry forest, soil microbial communities, metagenomic, seasonality, land use change
Citation: Lacerda-Júnior GV, Noronha MF, Cabral L, Delforno TP, de Sousa STP, Fernandes-Júnior PI, Melo IS and Oliveira VM (2019) Land Use and Seasonal Effects on the Soil Microbiome of a Brazilian Dry Forest. Front. Microbiol. 10:648. doi: 10.3389/fmicb.2019.00648
Received: 14 December 2018; Accepted: 14 March 2019;
Published: 05 April 2019.
Edited by:
Eoin L. Brodie, United States Department of Energy (DOE), United StatesReviewed by:
Damien Robert Finn, Arizona State University, United StatesHuaiying Yao, Institute of Urban Environment (CAS), China
Copyright © 2019 Lacerda-Júnior, Noronha, Cabral, Delforno, de Sousa, Fernandes-Júnior, Melo and Oliveira. 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: Gileno V. Lacerda-Júnior, Z2lsZW5vbGFjZXJkYWpyQGdtYWlsLmNvbQ== Paulo I. Fernandes-Júnior, cGF1bG8uaXZhbkBlbWJyYXBhLmJy Valéria M. Oliveira, dm1haWFAY3BxYmEudW5pY2FtcC5icg==