==== Front Braz J MicrobiolBraz. J. MicrobiolBrazilian Journal of Microbiology1517-83821678-4405Elsevier S1517-8382(17)30177-610.1016/j.bjm.2017.08.010Short CommunicationArchaea diversity in vegetation gradients from the Brazilian Cerrado Araujo Ademir Sergio Ferreira de asfaruaj@yahoo.com.bra⁎Mendes Lucas Wiliam bBezerra Walderly Melgaço cNunes Luis Alfredo Pinheiro Leal aLyra Maria do Carmo Catanho Pereira de dFigueiredo Marcia do Vale Barreto dMelo Vania Maria Maciel ca Federal University of Piauí, Agricultural Science Center, Soil Quality Lab., Teresina, PI, Brazilb University of Sao Paulo USP, Center for Nuclear Energy in Agriculture CENA, Cell and Molecular Biology Laboratory, Piracicaba, SP, Brazilc Federal University of Ceara, Laboratório de Ecologia Microbiana e Biotecnologia, Fortaleza, CE, Brazild Agronomic Institute of Pernambuco, Recife, PE, Brazil⁎ Corresponding author at: Soil Quality Lab., Agricultural Science Center, Federal University of Piauí, Teresina, PI, Brazil. asfaruaj@yahoo.com.br11 2 2018 Jul-Sep 2018 11 2 2018 49 3 522 528 2 3 2017 25 8 2017 © 2018 Sociedade Brasileira de Microbiologia. Published by Elsevier Editora Ltda.2018Sociedade Brasileira de MicrobiologiaThis is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).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 ecologySoil microbiology16S rRNATropical soilsAssociate Editor: Valeria Oliveira ==== Body 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.1 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.2, 3 Since the recognition of the archaeal domain in the tree of life,4 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.5 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 equitans6 and Thaumarchaeota.7 Most recently, two superphylum have been proposed, namely (i) TACK, which includes the phyla Thaumarchaeota, Aigarchaeota, Crenarchaeota and Korarchaeota,8 and (ii) Asgard, which consist in a sister group to TACK and are considered more closely related to the original eukaryote.9 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.10 Communities of autotrophic and heterotrophic Archaea catalyze iron and sulfur oxidation, which influence the release rate of metals and sulfur to the environment.11 Also, regarding the carbon cycle, all known methanogen organisms belong exclusively to Archaea domain and are generally found in oxygen-depleted environments.12 In Brazil, studies regarding the diversity of Archaea have been conducted for some important ecosystems, such as Amazonian and Atlantic forests,13, 14 as well as in mangroves.15 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.3, 16 However, it has been hypothesized that the vegetation and soil conditions in the Cerrado from Northeastern Brazil differs from those of the Central Plateau.17 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.18 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).Fig. 1 Map presenting the gradient of Cerrado at Sete Cidades National Park, Brazil. Table 1 Vegetation diversity indices in the Cerrado areas.19 Table 1 Campo graminoide Cerrado stricto sensu Cerradao Floresta decidual Plant richnesse 4.7 11 17 18 Plant diversityf 0.2 0.85 1.10 1.11 Plant densityg 4.7 27.1 35.0 51.8 Vegetationh a b c d a Andropogon fastigiatus; Aristida longifolia; Eragrostis maypurensis. b Andropogon fastigiatus; Aristida longifolia; Terminalia fagifolia; Magonia pubescens; Hymenaea courbaril; Plathymenia reticulata; Qualea grandiflora; Combretum mellifluum; Lippia origanoides; Anacardium occidentale; Simarouba versicolor; Vatairea macrocarpa. c Aspidosperma discolor; Parkia platycephala; Terminalia fagifolia; Piptadenia moniliformis; Plathymenia reticulata; Qualea parviflora; Anacardium occidentale; Copaifera coriacea; Thiloa glaucocarpa; Casearia grandiflora. d Aspidosperma multiflorum; Aspidosperma subincanum; Campomanesia aromática; Casearia lasiophylla; Casearia ulmifolia; Copaifera coriacea; Ephedranthus pisocarpus; Piptadenia moniliformis; Pterocarpus violaceus; Thiloa galucocarpa. e Species/100 m2. f H/100 m2. g Individual/100 m2. h Species of plants present. 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.20 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).21 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 follows21: 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 Project22 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).23 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)24 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.25 To determine the differences in abundance of archaeal groups among soil samples, the Statistical Analysis of Metagenomic Profiles (STAMP) software package was used.26 The OTU table was used as input, and P-values were calculated using the two-sided Welch's t-test,27 while confidence intervals were calculated using DP Welch's inverted and correction was made using Benjamini–Hochberg false discovery rate.28 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. Table 2Site Moisture (%) Temperature (°C) TOC (g kg−1) pH P (mg kg−1) K (cmolc kg−1) CEC (cmolc kg−1) Campo graminoide 7.3c 32a 4.3c 4.5c 3.9c 1.4b 2.28b Cerrado stricto sensu 10.5b 30b 8.3b 4.3b 3.9b 1.8b 2.31b Cerradao 11.9b 30b 9.1b 4.6b 4.5b 1.6b 2.35b Floresta decidual 31.8a 28c 15.2a 4.9a 5.3a 3.8a 4.91a TOC, total organic C; CEC, cation exchange capacity. Values followed by the same letter within each column are not significantly different at the 5% level, as determined by Student's t-test. 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.3 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.18 Interestingly, our results highlight the soil pH as a significant driver of archaea communities in these soils. According to Gorres,29 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.18 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.30 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.31Fig. 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.3, 16 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.32 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,33 and although ammonia oxidation occurs in both Bacteria and Archaea, the archaeal group dominates in both soil and marine environments.34 Nitrogen-related microorganisms are commonly found abundantly in forest soils,30 which may explain the higher abundance of Nitrosphaerales in our forest samples. Also, Navarrete et al.35 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.36 Members of this group are commonly found in oligotrophic soils, such as boreal forests,37 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.38 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.39 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.40 In our samples, the high soil pH found in the Floresta decidual area favored the dominance of Nitrosphaerales, which is sensitive to low pH.41 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.,41 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. Conflicts of interest The authors declare no conflicts of interest. 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. ==== Refs References 1 Schimel J. Playing scales in the methane cycle: from microbial ecology to the globe Proc Natl Acad Sci USA 101 2004 12400 12401 15314221 2 Siles J.A. 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 72 2016 207 220 26961712 3 Catao E. Castro A.P. Barreto C.C. Kruger R.H. Kyaw C.M. Diversity of archaea in Brazilian savana soils Arch Microbiol 195 2013 507 512 23515915 4 Woese C.R. Kandler O. Wheelis M.L. Towards a natural system of organisms: proposal for the domains Archaea, Bacteria, and Eucarya Proc Natl Acad Sci USA 87 1990 4576 4579 2112744 5 Schloss P.D. Girard R.A. Martin T. Status of the Archaeal and Bacterial census: an update MBio 7 2016 e00201 e216 27190214 6 Huber H. Hohn M.J. Rachel R. Fuchs T. Wimmer V.C. Stetter K.O. A new phylum of Archaea represented by a nanosized hyperthermophilic symbiont Nature 417 2002 63 67 11986665 7 Brochier-Armanet C. Boussau B. Gribaldo S. Forterre P. Mesophilic Crenarchaeota: proposal for a third archaeal phylum, the Thaumarchaeota Nat Rev Microbiol 6 2008 245 252 18274537 8 Guy L. Ettema T.J.G. The archaeal ‘TACK’ superphylum and the origin of eukaryotes Trends Microbiol 19 2011 580 587 22018741 9 Zaremba-Niedzwiedzka K. Caceres E.F. Saw J.H. Asgard archaea illuminate the origin of eukaryotic cellular complexity Nature 541 2017 353 358 28077874 10 Cabello P. Roldán M.D. Moreno-Vivián C. Nitrate reduction and the nitrogen cycle in archaea Microbiology 150 2004 3527 3546 15528644 11 Baker B.J. Banfield J.F. Microbial communities in acid mine drainage FEMS Microbiol Ecol 44 2003 139 152 19719632 12 Spang A. Ettema T.J.G. The methanogenic roots of Archaea Nat Microbiol 2 2017 17109 28741608 13 Taketani R.G. Tsai S.M. infuence of diferent land uses on the structure of archaeal communities in Amazonian anthrosols based on 16S rRNA and amoA genes Microb Ecol 59 2010 734 743 20204349 14 Etto R.M. Cruz L.M. Jesus E.C. Prokaryotic communities of acidic peatlands from the Southern Brazilian Atlantic Forest Braz J Microbiol 43 2012 661 674 24031878 15 Mendes L.W. Taketani R.G. Navarrete A.A. Tsai S.M. Shifts in phylogenetic diversity of archaeal communities in mangrove sediments at different sites and depths in southeastern Brazil Res Microbiol 163 2012 366 377 22683855 16 Castro A.P. Silva M.R.S.S. Quirino B.F. Bustamante M.M.C. Krüger R.H. Microbial diversity in Cerrado biome (Neotropical Savanna) soils PLoS One 11 2016 e0148785 26849674 17 Castro A.A.J.F. Martins F.R. Fernandes A.G. The woody flora of Cerrado vegetation in the State of Piauí, Northeastern Brazil Edinb J Bot 55 1998 455 472 18 Araujo A.S.F. Bezerra W.M. Santos V.M. Distinct bacterial communities across a gradient of vegetation from a preserved Brazilian Cerrado Anto van Leeuw 110 2017 457 469 19 Oliveira M.E.A. Martins F.R. Castro A.A.J.F. Santos J.R. 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 13 2007 1775 1783 20 Tedesco M.J. Gianello C. Bissani C.A. Analises de solos plantas e outros materiais 1995 UFRGS Porto Alegre 21 Caporaso J.G. Lauber C.L. Walters W.A. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample Proc Natl Acad Sci USA 108 2011 4516 4522 20534432 22 Pylro V.S. Roesch L.F.W. Ortega J.M. Brazilian Microbiome Project: revealing the unexplored microbial diversity – challenges and prospects Microb Ecol 67 2014 237 241 24173537 23 DeSantis T.Z. Hugenholtz P. Larsen N. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB App Environ Microbiol 72 2006 5069 5072 24 Anderson M. A new method for non-parametric multivariate analysis of variance Aust Ecol 26 2001 32 46 25 Hammer Ø. Harper D.A.T. Ryan P.D. Past: Paleontological Statistics Software package for education and data analysis Paleontol Elect 4 2001 9 26 Parks D.H. Tyson G.W. Hugenholtz P. Beiko R.G. STAMP: statistical analysis of taxonomic and functional profiles Bioinformatics 30 2014 3123 3124 25061070 27 Welch B.L. The generalization of “Student's” problem when several different population variances are involved Biometrika 34 1947 28 35 20287819 28 Benjamini Y. Hochberg Y. Controling the false discovery rate – a practical and powerful approach to multiple testing J Royal St Soc 57 1995 289 300 29 Gorres C.-M. Conrad R. Petersen S.O. Effect of soil properties and hydrology on Archaeal community composition in three temperate grasslands on peat FEMS Microbiol Ecol 85 2013 227 240 23521431 30 Mendes L.W. Tsai S.M. Navarrete A.A. Hollander M. van Veen J.A. Kuramae E.E. Soil-borne microbiome: linking diversity to function Microb Ecol 70 2015 255 265 25586384 31 Nicol G.W. Glover L.A. Prosser J.I. Molecular analysis of methanogenic archaeal communities in man-aged and natural upland pasture soils Glob Change Biol 9 2003 1451 1457 32 Nicol G.W. Tscherko D. Embley T.M. Prosser J.I. Primary succession of soil Crenarchaeota across a receding glacier foreland Environ Microbiol 7 2005 337 347 15683394 33 Stieglmeier M. Mooshammer M. Kitzler B. Aerobic nitrous oxide production through N -nitrosating hybrid formation in ammonia-oxidizing archaea ISME J 8 2014 1135 1146 24401864 34 Leininger S. Urich T. Schloter M. Archaea predominate among ammonia-oxidizing prokaryotes in soils Nature 442 2006 806 809 16915287 35 Navarrete A.A. Taketani R.G. Mendes L.W. Cannavan F.S. Moreira F.M.S. Tsai S.M. Land-use systems affect archaeal community structure and functional diversity in Western Amazon soils Rev Bras Ci Solo 35 2011 1527 1540 36 Vetriani C. Jannasch J.W. MacGregor B.J. Stahl D.A. Reysenbach A.-L. Population structure and phylogenetic characterization of marine benthic archaea in deep-sea sediments Appl Environ Microbiol 65 1999 4375 4384 10508063 37 Jurgens G. Lindström K. Saano A. Novel group within the kingdom Crenarchaeota from boreal forest soil Appl Environ Microbiol 63 1997 3090 3095 16535669 38 Iino T. Tamaki H. Tamazawa S. 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 28 2013 244 250 39 Tripathi B.M. Kim M. Lai-Hoe A. pH dominates variation in tropical soil archaeal diversity and community structure FEMS Microbiol Ecol 86 2013 303 311 23773164 40 Gubry-Rangin C. Hai B. Quince C. Niche specialization of terrestrial archaeal ammonia oxidizers Proc Natl Acad Sci USA 108 2011 21206 21211 22158986 41 Hu H.-W. Zhang L.-M. Yuan C.-L. The large-scale distribution of ammonia oxidizers in paddy soils is driven by soil pH, geographic distance, and climatic factors Front Microbiol 6 2015 938 946 26388866