Acessibilidade / Reportar erro

Archaea diversity in vegetation gradients from the Brazilian Cerrado

Abstract

We used 16S rRNA sequencing to assess the archaeal communities across a gradient of Cerrado. The archaeal communities differed across the gradient. Crenarcheota was the most abundant phyla, with Nitrosphaerales and NRPJ as the predominant classes. Euryachaeota was also found across the Cerrado gradient, including the classes Metanocellales and Methanomassiliicoccaceae.

Keywords:
Microbial ecology; Soil microbiology; 16S rRNA; Tropical soils

Soil microorganisms constitute a third of the Earth's biomass, playing important ecological and biogeochemical roles in terrestrial ecology. Among the microorganisms that inhabit soils, the domain Archaea is particularly important because it is abundant and plays important roles in C and N cycling.11 Schimel J. Playing scales in the methane cycle: from microbial ecology to the globe. Proc Natl Acad Sci USA. 2004;101:12400-12401. This domain constitutes one of the three major evolutionary lineages of life on Earth, and although it has long been believed that these organisms mainly inhabit extreme environments, previous studies demonstrated that they are widely distributed around the world in different types of soils.22 Siles JA, Margesin R. Abundance and diversity of bacterial, archaeal, and fungal communities along an altitudinal gradient in Alpine Forest Soils: what are the driving factors?. Microb Ecol. 2016;72:207-220.,33 Catao E, Castro AP, Barreto CC, Kruger RH, Kyaw CM. Diversity of archaea in Brazilian savana soils. Arch Microbiol. 2013;195:507-512.

Since the recognition of the archaeal domain in the tree of life,44 Woese CR, Kandler O, Wheelis ML. Towards a natural system of organisms: proposal for the domains Archaea, Bacteria, and Eucarya. Proc Natl Acad Sci USA. 1990;87:4576-4579. studies have been focused to explore their presence and function in a wide range of environments. With the advance of DNA-based molecular tools and genomic analysis, the knowledge of this domain has evolved and the current archaeal classification is rapidly moving. Presently, 20 archaeal phyla are recognized by small subunit RNA databases; however, it has been reported that 14 phyla have no cultured representative.55 Schloss PD, Girard RA, Martin T. Status of the Archaeal and Bacterial census: an update. MBio. 2016;7:e00201-e216. Most of culturable members of Archaea belong to the two main phyla Euryarchaeota and Crenarchaeota. Other new phyla have been proposed in the last years, for example Nanoarchaeota, based on the isolation of a hypertehermophilic symbiont Nanoarchaeum equitans66 Huber H, Hohn MJ, Rachel R, Fuchs T, Wimmer VC, Stetter KO. A new phylum of Archaea represented by a nanosized hyperthermophilic symbiont. Nature. 2002;417:63-67. and Thaumarchaeota.77 Brochier-Armanet C, Boussau B, Gribaldo S, Forterre P. Mesophilic Crenarchaeota: proposal for a third archaeal phylum, the Thaumarchaeota. Nat Rev Microbiol. 2008;6:245-252. Most recently, two superphylum have been proposed, namely (i) TACK, which includes the phyla Thaumarchaeota, Aigarchaeota, Crenarchaeota and Korarchaeota,88 Guy L, Ettema TJG. The archaeal ‘TACK' superphylum and the origin of eukaryotes. Trends Microbiol. 2011;19:580-587. and (ii) Asgard, which consist in a sister group to TACK and are considered more closely related to the original eukaryote.99 Zaremba-Niedzwiedzka K, Caceres EF, Saw JH, et al. Asgard archaea illuminate the origin of eukaryotic cellular complexity. Nature. 2017;541:353-358. Together with the advance of Archaea classification, their functional role in the ecosystem has also been the focus of recent studies. Members of Archaea carry out many steps in the nitrogen cycle, such as nitrate-based respiration and denitrification.1010 Cabello P, Roldán MD, Moreno-Vivián C. Nitrate reduction and the nitrogen cycle in archaea. Microbiology. 2004;150:3527-3546. Communities of autotrophic and heterotrophic Archaea catalyze iron and sulfur oxidation, which influence the release rate of metals and sulfur to the environment.1111 Baker BJ, Banfield JF. Microbial communities in acid mine drainage. FEMS Microbiol Ecol. 2003;44:139-152. Also, regarding the carbon cycle, all known methanogen organisms belong exclusively to Archaea domain and are generally found in oxygen-depleted environments.1212 Spang A, Ettema TJG. The methanogenic roots of Archaea. Nat Microbiol. 2017;2:17109.

In Brazil, studies regarding the diversity of Archaea have been conducted for some important ecosystems, such as Amazonian and Atlantic forests,1313 Taketani RG, Tsai SM. infuence of diferent land uses on the structure of archaeal communities in Amazonian anthrosols based on 16S rRNA and amoA genes. Microb Ecol. 2010;59:734-743.,1414 Etto RM, Cruz LM, Jesus EC, et al. Prokaryotic communities of acidic peatlands from the Southern Brazilian Atlantic Forest. Braz J Microbiol. 2012;43:661-674. as well as in mangroves.1515 Mendes LW, Taketani RG, Navarrete AA, Tsai SM. Shifts in phylogenetic diversity of archaeal communities in mangrove sediments at different sites and depths in southeastern Brazil. Res Microbiol. 2012;163:366-377. In the case of the Brazilian Cerrado, previous studies investigated the diversity of archaea in soils from the Central Plateau and found that Crenarcheota was the most abundant archaeal group in the ecosystem.33 Catao E, Castro AP, Barreto CC, Kruger RH, Kyaw CM. Diversity of archaea in Brazilian savana soils. Arch Microbiol. 2013;195:507-512.,1616 Castro AP, Silva MRSS, Quirino BF, Bustamante MMC, Krüger RH. Microbial diversity in Cerrado biome (Neotropical Savanna) soils. PLoS One. 2016;11:e0148785. However, it has been hypothesized that the vegetation and soil conditions in the Cerrado from Northeastern Brazil differs from those of the Central Plateau.1717 Castro AAJF, Martins FR, Fernandes AG. The woody flora of Cerrado vegetation in the State of Piauí, Northeastern Brazil. Edinb J Bot. 1998;55:455-472. Therefore, the pattern of Archaea communities would be different because the vegetation and soil conditions drive soil microorganisms. Previous studies have shown differences between bacterial diversity in the soils obtained from Cerrado in the Northeast when compared to those from the Central Plateau.1818 Araujo ASF, Bezerra WM, Santos VM, et al. Distinct bacterial communities across a gradient of vegetation from a preserved Brazilian Cerrado. Anto van Leeuw. 2017;110:457-469. Therefore, we used next-generation sequencing of the 16S rRNA gene in soil samples from vegetation gradients obtained from Campo graminoide, Cerrado stricto sensu, Cerradao and Floresta decidual to understand the diversity of Archaea in the soils of Cerrado in Northeastern Brazil.

The study was conducted in a preserved Brazilian Cerrado within Sete Cidades National Park (PNSC) (04°02′-08′S and 41°40′-45′W), Northeast, Brazil. This region presents two distinct seasons (wet and dry) with an annual average temperature at 25 °C and rainfall of 1558 mm distributed in February through April. We evaluated four different areas (with 1000 m2 each one) in a gradient of vegetation (‘Campo graminoide', ‘Cerrado stricto sensu', ‘Cerradao' and ‘Floresta decidual') (Table 1; Fig. 1).

Table 1
Vegetation diversity indices in the Cerrado areas.1919 Oliveira MEA, Martins FR, Castro AAJF, Santos JR. Classes de cobertura vegetal do parque nacional de sete cidades (transicao campo-floresta) utilizando imagens TM/Landsat, NE do Brasil. Proceedings XIII Simposio Brasileiro de Sensoriamento Remoto Anais. 2007;13:1775-1783.

Fig. 1
Map presenting the gradient of Cerrado at Sete Cidades National Park, Brazil.

Each area was separated in three transects (replication) where soil samples were collected at top soil layer (0-20 cm depth; three points per transect) in March (wet season), 2014. All soil samples were immediately stored in sealed plastic bags and transported in an icebox to the laboratory. A portion of the soil samples was stored in bags and kept at -20 °C for DNA analysis and another portion was air-dried, sieved through a 2-mm screen and homogenized for chemical analyses using standard laboratory protocols.2020 Tedesco MJ, Gianello C, Bissani CA. Analises de solos plantas e outros materiais. Porto Alegre: UFRGS; 1995.

Soil DNA was extracted from 0.5 g (total humid weight) of soil using the PowerLyzer PowerSoil DNA Isolation Kit (MoBIO Laboratories, Carlsbad, CA, USA), according to the manufacturer's instructions. The quality of the extracted DNA was checked with the Nanodrop ND-1000 spectrophotometer (Thermo Scientific, Waltham, MA, USA), and it was quantified by Qubit® 2.0 fluorometer using the dsDNA BR Assay kit (Invitrogen™). The integrity of the DNA was also confirmed by electrophoresis in a 0.8% agarose gel with 1× TAE buffer. The amplicon library of the archaeal 16S rRNA genes were amplified using the region-specific primers (515F/806R).2121 Caporaso JG, Lauber CL, Walters WA, et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc Natl Acad Sci USA. 2011;108:4516-4522. The first step amplification comprised 25 µL reaction containing the following: 14.8 µL of nuclease-free water (Certified Nuclease-free, Promega, Madison, WI, USA), 2.5 µL of 10× High Fidelity PCR Buffer (Invitrogen, Carlsbad, CA, USA), 1.0 µL of 50 mM MgSO4, 0.5 µL of each primer (10 µM concentration, 200 pM final concentration), 1.0 unit of Platinum Taq polymerase High Fidelity (Invitrogen, Carlsbad, CA, USA), and 4.0 µL of template DNA (10 ng). The conditions for PCR were as follows2121 Caporaso JG, Lauber CL, Walters WA, et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc Natl Acad Sci USA. 2011;108:4516-4522.: 94 °C for 4 min to denature the DNA, with 35 cycles at 94 °C for 45 s, 50 °C for 60 s, and 72 °C for 2 min, with a final extension of 10 min at 72 °C. In the second step, a unique pair of Illumina Nextera XT indexes (Illumina, San Diego, CA) was added to both ends of the amplified products. Each 50 µL reaction contained the following: 23.5 µL of nuclease-free water (Certified Nuclease-free, Promega, Madison, WI, USA), 5.0 µL of 10× High Fidelity PCR Buffer (Invitrogen, Carlsbad, CA, USA), 4.8 µL of 25 mM MgSO4, 1.5 µL of dNTP (10 mM each), 5.0 µL of each Nextera XT index (Illumina, San Diego, CA, USA), 1.0 unit of Platinum Taq polymerase High Fidelity (Invitrogen, Carlsbad, CA, USA), and 5.0 µL of each product from previous PCR. The conditions for this second round PCR were as follows: 95 °C for 3 min to denature the DNA, with 8 cycles at 95 °C for 30 s, 55 °C for 30 s, and 72 °C for 30 min, with a final extension of 5 min at 72 °C.

After indexing, the PCR products were cleaned up using Agencourt AMPure XP - PCR purification beads (Beckman Coulter, Brea, CA, USA), according to the manufacturer's manual, and quantified using the dsDNA BR assay Kit (Invitrogen, Carlsbad, CA, USA) on a Qubit 2.0 fluorometer (Invitrogen, Carlsbad, CA, USA). Once quantified, different volumes of each library were pooled into a single tube in equimolar concentration. After quantification, the molarity of the pool was determined and diluted to 2 nM, denatured, and then diluted to a final concentration of 8.0 pM with a 20% PhiX (Illumina, San Diego, CA, USA) spike for loading into the Illumina MiSeq sequencing machine (Illumina, San Diego, CA, USA).

Sequence data were processed using QIIME following the UPARSE standard pipeline according to Brazilian Microbiome Project2222 Pylro VS, Roesch LFW, Ortega JM, et al. Brazilian Microbiome Project: revealing the unexplored microbial diversity - challenges and prospects. Microb Ecol. 2014;67:237-241. to produce an OTU table and a set of representative sequences. Briefly, the reads were truncated at 240 bp and quality-filtered using a maximum expected error value of 0.5. Pre-filtered reads were dereplicated and singletons were removed and filtered for additional chimeras using the RDP_gold database using USEARCH 7.0. These sequences were clustered into OTUs at a 97% similarity cutoff following the UPARSE pipeline. After clustering, the sequences were aligned and taxonomically classified against the Greengenes database (version 13.8).2323 DeSantis TZ, Hugenholtz P, Larsen N, et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. App Environ Microbiol. 2006;72:5069-5072. The sequences were submitted to the NCBI Sequence Read Archive under the number SRP091586.

Redundancy Analysis (RDA) was used to determine the correlation between archaeal community structure and soil physicochemical properties. The OTU table was initially analyzed using Detrended Correspondence Analysis (DCA) to evaluate the gradient size of the species distribution, which indicated linearly distributed data (length of gradient <3), revealing that the best-fit model for the data was RDA. Forward selection (FS) and the Monte Carlo permutation test were applied with 1000 random permutations to verify the significance of soil chemical properties upon a microbial community. We used permutational multivariate analysis of variance (PERMANOVA)2424 Anderson M. A new method for non-parametric multivariate analysis of variance. Aust Ecol. 2001;26:32-46. to test whether sample categories harbored significantly different archaeal community structure. RDA plots were generated using the software Canoco 4.5 (Biometrics, Wageningen, The Netherlands), and PERMANOVA and alpha diversity index were calculated with the software PAST.2525 Hammer Ø, Harper DAT, Ryan PD. Past: Paleontological Statistics Software package for education and data analysis. Paleontol Elect. 2001;4:9. To determine the differences in abundance of archaeal groups among soil samples, the Statistical Analysis of Metagenomic Profiles (STAMP) software package was used.2626 Parks DH, Tyson GW, Hugenholtz P, Beiko RG. STAMP: statistical analysis of taxonomic and functional profiles. Bioinformatics. 2014;30:3123-3124. The OTU table was used as input, and P-values were calculated using the two-sided Welch's t-test,2727 Welch BL. The generalization of "Student's" problem when several different population variances are involved. Biometrika. 1947;34:28-35. while confidence intervals were calculated using DP Welch's inverted and correction was made using Benjamini-Hochberg false discovery rate.2828 Benjamini Y, Hochberg Y. Controling the false discovery rate - a practical and powerful approach to multiple testing. J Royal St Soc. 1995;57:289-300.

In this study, we examined the archaeal structure and composition in the soils from vegetation gradients from the Cerrado. We used 16S rRNA sequencing to assess the archaeal communities and correlated them with soil chemical parameters. Our results showed that soil parameters differed between the sampling sites. Areas under Campo graminoide, Cerrado stricto sensu and Cerradao presented the highest soil temperature and the lowest N, P, K, and CEC compared with Floresta decidual (Table 2). Soil moisture, pH, and TOC content were similar between Cerrado stricto sensu and Cerradao, and they were the lowest and highest in Campo graminoide and Floresta decidual, respectively.

Table 2
Soil physicochemical properties at different sites across the gradient of Cerrado.

These markedly different soil conditions, mainly between Campo graminoide and Floresta decidual, exist because of the differences in plant cover and the communities, both of which can influence soil properties. For example, Campo graminoide is covered mostly by grasses, whereas Floresta decidual is covered by trees. RDA showed the differences between the archaeal communities across the gradients (Fig. 2) and indicated that the Floresta decidual was the most distinct compared to the other areas.

Fig. 2
Redundancy analysis (RDA) of archaeal community patterns and soil characteristics from samples of Campo graminoide, Cerrado stricto sensu, Cerradao and Floresta decidual. Arrows indicate correlation between environmental parameters and archaeal profile. The significance of these correlations were evaluated via the Monte Carlo permutation test and is indicated by * (P < 0.05).

Additionally, RDA showed that the composition of archaea communities was strongly correlated with specific soil properties across the Cerrado gradients. Our results showed that the archaeal community of Floresta decidual correlated with soil moisture, TOC, CEC, P, K, and N. On the other hand, soil temperature and pH correlated with Campo graminoide, Cerrado stricto sensu and Cerradao. These results indicate that the archaea community differs across the gradients of the Cerrado and is influenced by soil conditions. Because Campo graminoide, Cerrado stricto sensu and Cerradao had similar soil conditions, their archaea communities were more similar to each other when compared to Floresta decidual. According to Catão et al.33 Catao E, Castro AP, Barreto CC, Kruger RH, Kyaw CM. Diversity of archaea in Brazilian savana soils. Arch Microbiol. 2013;195:507-512. the soil conditions drive the differences in the diversity of Archaea in the Brazilian Cerrado. These authors compared and contrasted Cerrado denso and Mata de Galeria from the Brazilian Cerrado and found differences in the community structures of Archaea. The varied soil conditions found across the gradients of vegetation have influenced the bacterial diversity.1818 Araujo ASF, Bezerra WM, Santos VM, et al. Distinct bacterial communities across a gradient of vegetation from a preserved Brazilian Cerrado. Anto van Leeuw. 2017;110:457-469. Interestingly, our results highlight the soil pH as a significant driver of archaea communities in these soils. According to Gorres,2929 Gorres C-M, Conrad R, Petersen SO. Effect of soil properties and hydrology on Archaeal community composition in three temperate grasslands on peat. FEMS Microbiol Ecol. 2013;85:227-240. soil pH is the main driver of archaeal community structure in soils. This strong correlation between soil pH and the microbial community structure may be because pH and the other physicochemical soil parameters are so closely related. Most microorganisms have intracellular pH that varies in the range of pH 7 ± 1; therefore, a small variation of the environmental pH could expose microorganisms to stress. Archaea were detected in all analyzed samples and showed very low diversity compared to the bacterial community found in the same samples.1818 Araujo ASF, Bezerra WM, Santos VM, et al. Distinct bacterial communities across a gradient of vegetation from a preserved Brazilian Cerrado. Anto van Leeuw. 2017;110:457-469. Interestingly, Floresta decidual presented the lowest alpha diversity when compared to the other areas (Fig. 3). It has been suggested that an environment in equilibrium, such as forest, the ecosystem functioning is maintained based on lower levels of diversity, but high abundance of microorganisms; in contrast, environments under stress would present an increased diversity, which leads to a higher functional diversity and, consequently, the maintenance of essential ecosystem functions.3030 Mendes LW, Tsai SM, Navarrete AA, Hollander M, van Veen JA, Kuramae EE. Soil-borne microbiome: linking diversity to function. Microb Ecol. 2015;70:255-265. Considering that Cerrado areas are more prone to environmental stress, such as high temperature, low moisture and oligotrophy, the higher archaeal diversity would be expected.

Fig. 3
Diversity measurement based on Shannon's index of archaeal communities in soils from Campo graminoide, Cerrado stricto sensu, Cerradao and Floresta decidual.

A total of 3078 reads were clustered in 33 OTUs according to the Greengenes taxonomical classification. We found Crenarchaeota to be the most abundant phyla across the Cerrado gradients, among which Nitrosphaerales and NRPJ were the predominant classes (Fig. 4). Organisms belonging to the largely non-thermophilic Crenarchaeota appear to be ubiquitous in soil systems and are dominant over euryarchaeal populations in grassland soils.3131 Nicol GW, Glover LA, Prosser JI. Molecular analysis of methanogenic archaeal communities in man-aged and natural upland pasture soils. Glob Change Biol. 2003;9:1451-1457.

Fig. 4
Average relative abundance of the most abundant archaeal groups in soils from Campo graminoide, Cerrado stricto sensu, Cerradao and Floresta decidual as revealed by the 16S rRNA gene ribotyping. For each sample type, the number of replicates is n = 9. Different lower letters indicate significant difference (FDR, P < 0.05) among samples.

These findings agree with previous studies in the Brazilian Cerrado, which have reported Crenarcheota as most dominant phylum and Nitrosphaerales as an important group.33 Catao E, Castro AP, Barreto CC, Kruger RH, Kyaw CM. Diversity of archaea in Brazilian savana soils. Arch Microbiol. 2013;195:507-512.,1616 Castro AP, Silva MRSS, Quirino BF, Bustamante MMC, Krüger RH. Microbial diversity in Cerrado biome (Neotropical Savanna) soils. PLoS One. 2016;11:e0148785. The relative abundance of Nitrosphaerales and NRPJ differed across the regions. Floresta decidual showed the highest abundance of Nitrosphaerales, whereas Campo graminoide, Cerrado stricto sensu and Cerradao showed the highest abundance of NRPJ. Previously, Nicol et al.3232 Nicol GW, Tscherko D, Embley TM, Prosser JI. Primary succession of soil Crenarchaeota across a receding glacier foreland. Environ Microbiol. 2005;7:337-347. have indicated that Crenarchaeota from soils may have specific association with plant roots, playing an important role in the rhizosphere. Nitrosphaerales are important ammonia-oxidizer organisms,3333 Stieglmeier M, Mooshammer M, Kitzler B, et al. Aerobic nitrous oxide production through N-nitrosating hybrid formation in ammonia-oxidizing archaea. ISME J. 2014;8:1135-1146. and although ammonia oxidation occurs in both Bacteria and Archaea, the archaeal group dominates in both soil and marine environments.3434 Leininger S, Urich T, Schloter M, et al. Archaea predominate among ammonia-oxidizing prokaryotes in soils. Nature. 2006;442:806-809. Nitrogen-related microorganisms are commonly found abundantly in forest soils,3030 Mendes LW, Tsai SM, Navarrete AA, Hollander M, van Veen JA, Kuramae EE. Soil-borne microbiome: linking diversity to function. Microb Ecol. 2015;70:255-265. which may explain the higher abundance of Nitrosphaerales in our forest samples. Also, Navarrete et al.3535 Navarrete AA, Taketani RG, Mendes LW, Cannavan FS, Moreira FMS, Tsai SM. Land-use systems affect archaeal community structure and functional diversity in Western Amazon soils. Rev Bras Ci Solo. 2011;35:1527-1540. showed that the diversity of ammonia-oxidizing archaeal communities tend to be higher in primary forests. On the other hand, the NRPJ group was more abundant in Cerrado sites. This group belongs to the Marine Benthic Group A, an uncultured archaeal cluster initially described from cold continental slope and deep-sea sediments.3636 Vetriani C, Jannasch JW, MacGregor BJ, Stahl DA, Reysenbach A-L. Population structure and phylogenetic characterization of marine benthic archaea in deep-sea sediments. Appl Environ Microbiol. 1999;65:4375-4384. Members of this group are commonly found in oligotrophic soils, such as boreal forests,3737 Jurgens G, Lindström K, Saano A. Novel group within the kingdom Crenarchaeota from boreal forest soil. Appl Environ Microbiol. 1997;63:3090-3095. which may explain its higher abundance in our Cerrado samples. In summary, our results showed a high abundance of Crenarchaeota groups in our samples, with differential abundance between sites, suggesting that members of this group play distinct and important roles in the soil ecosystem.

The phylum Euryarchaeota is considered the most physiologically diverse and includes methanogens, halophiles and thermophile organisms. Members of this phylum were also found across the Cerrado gradients, including the classes Metanocellales and Methanomassiliicoccaceae. Interestingly, the Metanocellales class was detected only at Campo graminoide, which also had a high abundance of Methanomassiliicoccaceae. Similarly, Cerrado stricto sensu showed high abundance of Methanomassiliicoccaceae. The family Methanomassiliicoccaceae comprises of a methanogenic lineage of the class Thermoplasmata, and members of the Methanocellales family are methanogenic organisms, contributing to the global methane cycle.3838 Iino T, Tamaki H, Tamazawa S, et al. Candidatus Methanogranum caenicola: a novel methanogen from the anaerobic digested sludge, and proposal of Methanomassiliicoccaceae fam. nov. and Methanomassiliicoccales ord. nov., for a methanogenic lineage of the class Thermoplasmata. Microb Environ. 2013;28:244-250. In this sense, our results show that important groups related to carbon cycle are abundant in Campo graminoide and Cerrado stricto sensu.

As shown in the RDA, pH was the main driver of the archaea communities across the gradients. It has been described that pH dominates the variation of archaeal diversity and community composition in tropical soils.3939 Tripathi BM, Kim M, Lai-Hoe A, et al. pH dominates variation in tropical soil archaeal diversity and community structure. FEMS Microbiol Ecol. 2013;86:303-311. In a large-scale study, the authors showed a niche specialization of terrestrial archaeal ammonia oxidizers based on pH, showing that amoA abundance and diversity increased with soil pH.4040 Gubry-Rangin C, Hai B, Quince C, et al. Niche specialization of terrestrial archaeal ammonia oxidizers. Proc Natl Acad Sci USA. 2011;108:21206-21211. In our samples, the high soil pH found in the Floresta decidual area favored the dominance of Nitrosphaerales, which is sensitive to low pH.4141 Hu H-W, Zhang L-M, Yuan C-L, et al. The large-scale distribution of ammonia oxidizers in paddy soils is driven by soil pH, geographic distance, and climatic factors. Front Microbiol. 2015;6:938-946. On the other hand, a high abundance of Euryachaeota, specifically Metanocellales, in the Campo graminoide region may be related to the low soil pH found in this area. According to Hu et al.,4141 Hu H-W, Zhang L-M, Yuan C-L, et al. The large-scale distribution of ammonia oxidizers in paddy soils is driven by soil pH, geographic distance, and climatic factors. Front Microbiol. 2015;6:938-946. soil pH is a driver of methanogenesis and subsequently may promote shifts in the archaeal community composition from Crenarchaeota to methanogenic Euryarchaeota.

In conclusion, our data showed that although the most of detected taxa are shared between all four areas, the archaeal structure and abundance differed across these gradients and was strongly driven by soil physicochemical properties. The soil pH was the main driver of archaeal communities in the studied soils, and this can be explained by the difference between the sites, where Floresta decidual presented higher values while Campo graminoide was more acid. Results showed that Crenarcheota was the most dominant phylum, followed by Euryarchaeota. Additionally, differential soil parameters lead to distinct functional potential of the communities, where Floresta decidual presented high abundance of groups related to the nitrogen metabolism while Campo graminoide and Cerrado stricto sensu presented high abundance of groups related to carbon metabolism. Also, Nitrosphaerales was the most abundant order found across the gradients, and further studies are needed to evaluate their potential in N cycling.

Acknowledgments

The authors thank "Fundação de Amparo a Pesquisa no Estado do Piauí" (FAPEPI) for financial support to this project through of PRONEX (FAPEPI/CNPq 004/2012). We thank the Unidade Multiusuário do Núcleo de Pesquisa e Desenvolvimento de Medicamentos (UM-NPDM) from Federal University of Ceará for using MiSeq Illumina sequencing facilities.

References

  • 1
    Schimel J. Playing scales in the methane cycle: from microbial ecology to the globe. Proc Natl Acad Sci USA 2004;101:12400-12401.
  • 2
    Siles JA, Margesin R. Abundance and diversity of bacterial, archaeal, and fungal communities along an altitudinal gradient in Alpine Forest Soils: what are the driving factors?. Microb Ecol 2016;72:207-220.
  • 3
    Catao E, Castro AP, Barreto CC, Kruger RH, Kyaw CM. Diversity of archaea in Brazilian savana soils. Arch Microbiol 2013;195:507-512.
  • 4
    Woese CR, Kandler O, Wheelis ML. Towards a natural system of organisms: proposal for the domains Archaea, Bacteria, and Eucarya. Proc Natl Acad Sci USA. 1990;87:4576-4579.
  • 5
    Schloss PD, Girard RA, Martin T. Status of the Archaeal and Bacterial census: an update. MBio 2016;7:e00201-e216.
  • 6
    Huber H, Hohn MJ, Rachel R, Fuchs T, Wimmer VC, Stetter KO. A new phylum of Archaea represented by a nanosized hyperthermophilic symbiont. Nature 2002;417:63-67.
  • 7
    Brochier-Armanet C, Boussau B, Gribaldo S, Forterre P. Mesophilic Crenarchaeota: proposal for a third archaeal phylum, the Thaumarchaeota. Nat Rev Microbiol 2008;6:245-252.
  • 8
    Guy L, Ettema TJG. The archaeal ‘TACK' superphylum and the origin of eukaryotes. Trends Microbiol. 2011;19:580-587.
  • 9
    Zaremba-Niedzwiedzka K, Caceres EF, Saw JH, et al. Asgard archaea illuminate the origin of eukaryotic cellular complexity. Nature 2017;541:353-358.
  • 10
    Cabello P, Roldán MD, Moreno-Vivián C. Nitrate reduction and the nitrogen cycle in archaea. Microbiology 2004;150:3527-3546.
  • 11
    Baker BJ, Banfield JF. Microbial communities in acid mine drainage. FEMS Microbiol Ecol 2003;44:139-152.
  • 12
    Spang A, Ettema TJG. The methanogenic roots of Archaea. Nat Microbiol. 2017;2:17109.
  • 13
    Taketani RG, Tsai SM. infuence of diferent land uses on the structure of archaeal communities in Amazonian anthrosols based on 16S rRNA and amoA genes. Microb Ecol. 2010;59:734-743.
  • 14
    Etto RM, Cruz LM, Jesus EC, et al. Prokaryotic communities of acidic peatlands from the Southern Brazilian Atlantic Forest. Braz J Microbiol. 2012;43:661-674.
  • 15
    Mendes LW, Taketani RG, Navarrete AA, Tsai SM. Shifts in phylogenetic diversity of archaeal communities in mangrove sediments at different sites and depths in southeastern Brazil. Res Microbiol. 2012;163:366-377.
  • 16
    Castro AP, Silva MRSS, Quirino BF, Bustamante MMC, Krüger RH. Microbial diversity in Cerrado biome (Neotropical Savanna) soils. PLoS One. 2016;11:e0148785.
  • 17
    Castro AAJF, Martins FR, Fernandes AG. The woody flora of Cerrado vegetation in the State of Piauí, Northeastern Brazil. Edinb J Bot. 1998;55:455-472.
  • 18
    Araujo ASF, Bezerra WM, Santos VM, et al. Distinct bacterial communities across a gradient of vegetation from a preserved Brazilian Cerrado. Anto van Leeuw. 2017;110:457-469.
  • 19
    Oliveira MEA, Martins FR, Castro AAJF, Santos JR. Classes de cobertura vegetal do parque nacional de sete cidades (transicao campo-floresta) utilizando imagens TM/Landsat, NE do Brasil. Proceedings XIII Simposio Brasileiro de Sensoriamento Remoto Anais 2007;13:1775-1783.
  • 20
    Tedesco MJ, Gianello C, Bissani CA. Analises de solos plantas e outros materiais Porto Alegre: UFRGS; 1995.
  • 21
    Caporaso JG, Lauber CL, Walters WA, et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc Natl Acad Sci USA. 2011;108:4516-4522.
  • 22
    Pylro VS, Roesch LFW, Ortega JM, et al. Brazilian Microbiome Project: revealing the unexplored microbial diversity - challenges and prospects. Microb Ecol. 2014;67:237-241.
  • 23
    DeSantis TZ, Hugenholtz P, Larsen N, et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. App Environ Microbiol. 2006;72:5069-5072.
  • 24
    Anderson M. A new method for non-parametric multivariate analysis of variance. Aust Ecol. 2001;26:32-46.
  • 25
    Hammer Ø, Harper DAT, Ryan PD. Past: Paleontological Statistics Software package for education and data analysis. Paleontol Elect 2001;4:9.
  • 26
    Parks DH, Tyson GW, Hugenholtz P, Beiko RG. STAMP: statistical analysis of taxonomic and functional profiles. Bioinformatics 2014;30:3123-3124.
  • 27
    Welch BL. The generalization of "Student's" problem when several different population variances are involved. Biometrika 1947;34:28-35.
  • 28
    Benjamini Y, Hochberg Y. Controling the false discovery rate - a practical and powerful approach to multiple testing. J Royal St Soc. 1995;57:289-300.
  • 29
    Gorres C-M, Conrad R, Petersen SO. Effect of soil properties and hydrology on Archaeal community composition in three temperate grasslands on peat. FEMS Microbiol Ecol. 2013;85:227-240.
  • 30
    Mendes LW, Tsai SM, Navarrete AA, Hollander M, van Veen JA, Kuramae EE. Soil-borne microbiome: linking diversity to function. Microb Ecol. 2015;70:255-265.
  • 31
    Nicol GW, Glover LA, Prosser JI. Molecular analysis of methanogenic archaeal communities in man-aged and natural upland pasture soils. Glob Change Biol. 2003;9:1451-1457.
  • 32
    Nicol GW, Tscherko D, Embley TM, Prosser JI. Primary succession of soil Crenarchaeota across a receding glacier foreland. Environ Microbiol. 2005;7:337-347.
  • 33
    Stieglmeier M, Mooshammer M, Kitzler B, et al. Aerobic nitrous oxide production through N-nitrosating hybrid formation in ammonia-oxidizing archaea. ISME J. 2014;8:1135-1146.
  • 34
    Leininger S, Urich T, Schloter M, et al. Archaea predominate among ammonia-oxidizing prokaryotes in soils. Nature 2006;442:806-809.
  • 35
    Navarrete AA, Taketani RG, Mendes LW, Cannavan FS, Moreira FMS, Tsai SM. Land-use systems affect archaeal community structure and functional diversity in Western Amazon soils. Rev Bras Ci Solo. 2011;35:1527-1540.
  • 36
    Vetriani C, Jannasch JW, MacGregor BJ, Stahl DA, Reysenbach A-L. Population structure and phylogenetic characterization of marine benthic archaea in deep-sea sediments. Appl Environ Microbiol. 1999;65:4375-4384.
  • 37
    Jurgens G, Lindström K, Saano A. Novel group within the kingdom Crenarchaeota from boreal forest soil. Appl Environ Microbiol. 1997;63:3090-3095.
  • 38
    Iino T, Tamaki H, Tamazawa S, et al. Candidatus Methanogranum caenicola: a novel methanogen from the anaerobic digested sludge, and proposal of Methanomassiliicoccaceae fam. nov. and Methanomassiliicoccales ord. nov., for a methanogenic lineage of the class Thermoplasmata. Microb Environ. 2013;28:244-250.
  • 39
    Tripathi BM, Kim M, Lai-Hoe A, et al. pH dominates variation in tropical soil archaeal diversity and community structure. FEMS Microbiol Ecol. 2013;86:303-311.
  • 40
    Gubry-Rangin C, Hai B, Quince C, et al. Niche specialization of terrestrial archaeal ammonia oxidizers. Proc Natl Acad Sci USA. 2011;108:21206-21211.
  • 41
    Hu H-W, Zhang L-M, Yuan C-L, et al. The large-scale distribution of ammonia oxidizers in paddy soils is driven by soil pH, geographic distance, and climatic factors. Front Microbiol. 2015;6:938-946.

Edited by

Associate Editor: Valeria Oliveira

Publication Dates

  • Publication in this collection
    Jul-Sep 2018

History

  • Received
    2 Mar 2017
  • Accepted
    25 Aug 2017
Sociedade Brasileira de Microbiologia USP - ICB III - Dep. de Microbiologia, Sociedade Brasileira de Microbiologia, Av. Prof. Lineu Prestes, 2415, Cidade Universitária, 05508-900 São Paulo, SP - Brasil, Ramal USP 7979, Tel. / Fax: (55 11) 3813-9647 ou 3037-7095 - São Paulo - SP - Brazil
E-mail: bjm@sbmicrobiologia.org.br