Archaea
Volume 2019 (2019), Article ID 6957210, 11 pages
https://doi.org/10.1155/2019/6957210

Research Article

Long-Term Effects of Periodical Fires on Archaeal Communities from Brazilian Cerrado Soils

Aline Belmok,¹ Thiago Rodrigues-Oliveira,¹ Fabyano A. C. Lopes,1,2 Heloisa S. Miranda,³ Ricardo H. Krüger,¹ and Cynthia M. Kyaw¹

¹Institute of Biological Sciences, Department of Cell Biology, University of Brasilia, 70910-900 Brasília, DF, Brazil
²Institute of Biological Sciences, Department of Biochemistry and Molecular Biology, Federal University of Goiás, 74690-900 Goiânia, GO, Brazil
³Institute of Biological Sciences, Department of Ecology, University of Brasilia, 70910-900 Brasília, DF, Brazil

Correspondence should be addressed to Cynthia M. Kyaw; malta@unb.br

Received 31 August 2018; Revised 11 October 2018; Accepted 18 November 2018; Published 27 January 2019

Academic Editor: William B. Whitman

Copyright © 2019 Aline Belmok et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The Cerrado biome corresponds to an extensive area of Brazil and is considered a biodiversity hotspot. Frequent fires are a natural feature in this biome and have influences on vegetation structure and composition. However, continuous anthropogenic actions are promoting changes in fire frequency and seasonality. Despite the high biodiversity of the Cerrado, little is known about its microbiome, with few publications describing some aspects of the bacterial and fungal communities found on this biome and almost no references about archaea. In this study, we describe the archaeal diversity in Cerrado sensu stricto soils, comparing the archaeal communities from soils of an area long protected from fires to one exposed to biennial fires, using both 16S rRNA and amoA genes as molecular markers. Almost all 16S rRNA sequences from both studied areas were affiliated with I.1b and 1.1c Thaumarchaeota, groups commonly detected in terrestrial environments. A higher relative abundance of I.1b thaumarchaeal subgroup was detected in the frequently burned area even though no statistically significant differences were observed in archaeal 16S rRNA richness and diversity between the investigated areas. Many ammonia-oxidizing archaea (AOA) are affiliated with this group, which is consistent with the higher amoA diversity and OTU numbers detected in the area periodically burned. Taken together, our results suggest that, although total archaeal community richness and diversity do not seem to greatly differ between the investigated conditions, alterations in wood cover and vegetation structure caused by frequent fires likely cause long-term effects in AOA diversity in Cerrado soils.

1. Introduction

It is well established that organisms of the Archaea domain are ubiquitously distributed and represent a significant fraction of the prokaryotic cells found in aquatic and terrestrial ecosystems [1, 2]. In general, archaea can account for 1-2% of total prokaryotes in moderate aerobic soils [3], reaching up to 5% in sandy soils [4] and even 38% in acidic forest soils [5]. Furthermore, the roles of soil archaea in the nitrogen cycle have been increasingly debated [6–8], with the frequent detection of ammonia-oxidizing archaea (AOA) in these environments [9, 10]. Although most archaea found in soils (and other mesophilic environments) were initially classified as Group I Crenarchaeota [11], they are now classified as members of the Thaumarchaeota phylum [12], which contains all AOA known so far [13].

Brazil is a country widely known for its continental proportions, with many different terrestrial landscapes and biomes. However, in spite of its notable biodiversity, studies focusing on archaeal diversity in natural Brazilian environments are still scarce, with only about 50 studies published in the last 20 years according to a recent review by Rodrigues et al. [14]. Most of these studies investigated archaeal communities in soils from the Amazon or Atlantic Forest biomes, but very few have focused on archaeal diversity of the Brazilian savanna, also known as Cerrado [15–19]. This biome is extensive and extremely diverse, corresponding originally to approximately 22% of the country’s territory. Due to the large number of endemic species and the continuous degradation of its habitats, the Cerrado is considered a hotspot for conservation priorities [20]. This biome presents different physiognomic forms, ranging from grassland to forest-like vegetation [21]. The most common physiognomy is the Cerrado sensu stricto, a savanna woodland with a continuous herbaceous layer of scattered trees and shrubs and canopy closure ranging from 20 to 50% [21]. The seasonality of rains, poor nutrient soils, and fires are considered the determinants of this biome vegetation [22, 23].

Natural fires occurred in the Cerrado for millennia [24] and are still present mostly during the transition from dry to wet season [25, 26]. However, continuous occupation and conversion of large areas of Cerrado to agriculture and pasture have increased fire frequency and altered the fire season. Although many Cerrado species have adaptations to fire [22], changes in fire regimes strongly influence the structure and composition of Cerrado’s vegetation, mostly due to mortality [27, 28], topkill and resprouting from underground organs [29, 30], alteration in population dynamics resulting from the loss of reproduction investment and seedling mortality [31, 32], and significant increase in grass cover [33]. The large reduction in canopy cover results in alterations in nutrient cycling mostly due to the reduction of litter formation [34], changes in temperature due to gaps in the herbaceous layer vegetation [35], and consequently, changes in soil humidity [36, 37], all of which may alter the composition of archaeal communities [38–40].

In this work, we present a case study to investigate if the archaeal communities from typical Cerrado sensu stricto areas with long protection from fires differ from those of areas submitted to biennial fires in mid dry season. Soil samples were collected from areas of well-known fire history. Twelve DNA libraries were obtained and sequenced, six using the 16S rRNA Archaea-specific gene and six using the archaeal amoA gene. The resulting sequences were analyzed to evaluate the differences between archaeal communities in Cerrado sensu stricto soils from frequently burned and fire-protected areas.

2. Material and Methods

2.1. Study Site and Soil Sampling

The study areas are located in the Ecological Reserve of the Brazilian Institute of Geography and Statistics (RECOR-IBGE), 35 km south of Brasilia, DF (15°55S, 47°51W). The Reserve has an area of 1360 ha, encompassing the most common Cerrado physiognomies. The climate in the area is Cwa in Köppen’s classification with two well-defined seasons: a dry season lasting from June to October and a rainy season when most of the average precipitation of 1436 mm occurs. Mean temperatures are around 21°C, and the altitude varies from 1048 to 1150 m. The soil is typically yellow-red latosol (Acrustox according to the American classification) [41].

In September 2013, two Cerrado sensu stricto areas with contrasting and well-known fire history were selected for this work. The first area was protected from fire for 20 years (Control area = C), and the second, 4 m apart, was submitted to a biennial fire regime in the mid dry season (August) from 1992 to 2008 (Burned area = Q, from the Portuguese word for burned—Queimada) (Figure 1). In September 2011, a wildfire burned Q without affecting C. Therefore, Q was protected from fire for 2 years prior to our study. Canopy cover in C and Q varied from 32% to 84% and from 28% to 64%. The aboveground biomass of the herbaceous layer in C was , with the grasses and forbs representing 3% and 12% of the total. The dry mass of litter was 5.4 Mg/ha. Total herbaceous biomass in Q was , with the grasses and forbs representing 35% and 26% of total, and the dry mass of litter was 3.0 Mg/ha.

Figure 1: Location of Cerrado sensu stricto sampling sites (indicated in red), at the Reserva Ecológica do IBGE, Brasilia, Brazil. The distribution of the Cerrado biome is indicated in green on the map. The Control area (C, on the left), that has been protected from fires for 20 years, and the Burned area, submitted to a biennial fire regime since 1992 (Q, on the right), are highlighted.

In each area, three soil samples were randomly collected, within a 10 m × 10 m plot, being at least 1 m apart from each other. The soil was collected to a depth of 10 cm with a PVC tube, 10 cm in diameter, and stored at −20°C until used. Physicochemical properties for the soil samples of each area were performed by SoloQuímica Análises de Solo Ltda. Principal component analyses (PCA) were performed on abundance measures using a correlation matrix in Factoextra v.1.0.3 package implemented by R v.3.2.2 [42]. Statistical differences in soil physical-chemical properties were compared by Welch’s -test () using R statistical software v.3.2.2 [42].

2.2. DNA Extraction and 16S rRNA and amoA Gene Library Construction

Total soil DNA was extracted from 0.5 g of each sample with PowerSoil DNA Isolation Kit (MO Bio Laboratories Inc.) according to the manufacturer’s instructions. PCR assays were conducted with primer pairs 21f/958r [43] and Arch amoAf/Arch amoAr [44], for archaeal 16S rRNA and amoA genes, respectively. All PCR assays were performed in 50 μL reaction mixtures, containing 1 to 100 ng of template DNA, 1X reaction buffer (Invitrogen®), 1.5 mM MgCl₂ (Invitrogen®), 400 ng/μL bovine serum albumin, 0.5 μM of each primer, 200 μM dNTPs (Invitrogen®), and 1.25 U Taq DNA Polymerase (Invitrogen®). The PCR experiments were performed in a MJ PTC-100® (Peltier Thermal Cycles) thermocycler, with cycling conditions previously described [43, 44]. Amplified DNA was visualized on 1% agarose gel electrophoresis stained with ethidium bromide (10 mg/mL). Amplicons were purified using GeneJET PCR Purification kit (Thermo Scientific®), cloned into the pGEM-T Easy® (Promega) vector, according to manufacturer’s instructions, and transformed into Escherichia coli DH5α competent cells by heat shock treatment. Plasmidial DNA of the recombinant clones was extracted by phenol-chloroform-isoamyl alcohol at 25 : 24 : 1 (vol/vol/vol) and was sequenced by the Sanger method at Macrogen Inc. (Korea).

2.3. Phylogenetic and Statistical Analyses

All sequences were trimmed according to Phred quality superior to 20 in more than 400 nucleotides. 16S rRNA gene sequences were used for comparative analyses with the Greengenes taxonomical database [45] using Mothur v.1.24.1 [46] with an identity threshold of 90% or higher. Although Greengenes does not consider Thaumarchaeota as a separate phylum (it is considered as a crenarchaeotal class), in this work we considered it in the phylum taxonomical level as suggested by Brochier-Armanet et al. [12]. Multiple alignments of 16S rRNA and amoA gene sequences were performed with Clustal X v. 2.1 [47]. Gap columns generated by the alignment processes were filtered using Mothur v.1.24., also employed for rarefaction curve estimations. For calculations of coverage, richness, and diversity indexes—Ace, Chao, and Shannon [48–50]—libraries were normalized for the number of sequences of the smallest dataset for each gene ( for 16S rRNA and for amoA) and significant differences in indicators were evaluated by Welch’s -test (). Analysis of both 16S rRNA and amoA gene sequences was OTU-based and the identity threshold adopted was 97% nucleotide sequence similarity [13, 51].

Relative abundances of 16S rRNA OTUs from each area were evaluated using STAMP v.2.1.3 [52]. The areas were compared by Welch’s -test (), and nonmetric multidimensional scaling (NMDS) for each profile was performed using the Bray-Curtis similarity index in R software [42]. To construct phylogenetic trees with both 16S rRNA and amoA genes, multiple alignments of the deduced sequences were performed in MUSCLE v. 3.8.31 [53] using default parameters and manually edited. The phylogenetic trees were constructed in FastTree v.2.1 with default parameters and 1,000 bootstrap test [54].

3. Results

3.1. Soil Physicochemical Characteristics

The physicochemical properties of Cerrado sensu stricto soils (Supplementary Table 1) revealed an acidic pH (around 4.9) in both Control and Burned areas. Principal component analysis (PCA) of soil physicochemical data (Figure 2(a)) indicated separation of soil samples from Control and Burned areas in the first dimension axis, suggesting higher influence of Mn in soils from the Burned area and Fe and ammonium in samples from the Control area. However, Welch’s -test revealed that, among all physicochemical parameters analyzed, only NH₄⁺-N was significantly different in the studied sites, with higher concentrations in the Control area (Figure 2(b)).

(a)

(b)

Figure 2: Principal component analysis (PCA) of soil physicochemical properties based on a correlation matrix performed in Factoextra R package v.1.0.3 (a). OM: organic matter; CEC: cation exchange capacity. (b) NH₄⁺-N concentrations in soils from Control and Burned areas, with asterisk indicating significant difference detected in Welch’s -test.

3.2. Archaeal 16S rRNA Gene Analysis

DNA amplifications of the archaeal 16S rRNA gene resulted in six clone libraries: three replicates from the Control area (Ca, Cb, and Cc) and three from the biennially Burned area (Qa, Qb, and Qc). In total, 466 sequences with high Phred quality (>20) over 400 bp were obtained from these six libraries (Supplementary Table 2), all classified as Archaea according to the Greengenes taxonomical database, resulting in 32 OTUs with a 97% sequence similarity cutoff.

Rarefaction analysis (Supplementary Figure 1) and high coverage estimators (Table 1) indicate that archaeal communities in typical Cerrado sensu stricto soils from both study sites were well covered, even at OTUs with 97% sequence similarity. A comparable number of OTUs were observed in each area with no statistically significant differences detected in Ace, Chao1, and Shannon values, suggesting similar archaeal richness and diversity in soils from both areas (Table 1). Furthermore, the majority of OTUs (15) identified in soils from the Control area (representing most of the sequences obtained) were also observed in soils from the Burned area, reinforcing the similarity of the archaeal community inferred by 16S rRNA genes in soils from both conditions analyzed (Figure 3(a)). These results are in accordance with those obtained in nonmetric multidimensional scaling (NMDS) analysis performed with this marker gene, where no clear clustering pattern among replicates from the Control and Burned soils could be detected (Supplementary Figure 2).

Table 1: α-Diversity analysis of archaeal 16S rRNA gene sequences from Cerrado sensu stricto soils in both conditions analyzed in this study. C = Control area, protected from fire for 20 years, and Q = Burned Area, biennially burned. Measures were calculated independently for each replicate and the mean and standard deviation (SD) of the indicators are shown.

Area
OTU observed
Chao1
ACE
Shannon
Coverage (%)

C

Q

(a)

(b)

Figure 3: Venn diagrams showing unique and shared archaeal 16S rRNA (a) and amoA (b) OTUs (97% sequence similarity) between Cerrado sensu stricto soils from the area protected from fire for 20 years (Control) and the area submitted to biennial fire regime (Burned).

Almost all 16S rRNA gene sequences (99.36%) were affiliated to either I.1b or I.1c subgroups of the phylum Thaumarchaeota and only three sequences (0.64%), all from the replicate Cb, were identified as belonging to the Bathyarchaeota phylum (Figure 4). It is noteworthy that, despite the similarity in richness indexes (Table 1) and numbers of observed OTUs shared between both studied conditions (Figure 3(a)), a higher relative abundance of I.1c thaumarchaeal subgroup sequences (77.68%) was retrieved from soils sampled in the Control area. In contrast, most of the sequences detected in soils from the frequently burned area were affiliated with I.1b Thaumarchaeota (58.80%), as evidenced in Figure 4. Furthermore, with the exception of Cb 54, all OTUs composed exclusively by sequences obtained from Control soils clustered within the I.1c thaumarchaeal clade (Figure 4).

Figure 4: Phylogenetic tree of archaeal 16S rRNA gene OTUs (97% sequence similarity) obtained from typical Cerrado sensu stricto soils. The number of sequences represented by each OTU is shown in parentheses, with those from the Control area (C) indicated by blue squares and from the Burned area (Q) by red circles. Singleton sequences were not included. The tree was constructed using FastTree v. 2.1 [54], with 1000 bootstrap tests. Values below 70% are not shown. HWCG III: hot water crenarchaeotic group III.

3.3. Archaeal amoA Gene Analysis

Archaeal amoA gene libraries were obtained from the same soil samples, resulting in 251 sequences with high Phred quality (>20) over 400 bp (Supplementary Table 2). A total of 12 OTUs were identified, using a 97% sequence similarity cutoff. Coverage estimators indicated that sampling efforts were also satisfactory when using this gene (Table 2). Unlike the results obtained for 16S rRNA sequences, significant differences were detected in Shannon index values for this marker gene (Table 2), suggesting higher amoA diversity in soils from the biennially burned area. In addition, despite the fact that most of the archaeal amoA sequences grouped in three highly representative OTUs could be detected in both Control and Burned soils (Figure 3(b)), a number of OTUs were observed only in soils from the Burned area, indicating a possible occurrence of unique AOA in soils submitted to a biennial fire regime.

Table 2: α-Diversity analysis of archaeal amoA gene sequences from Cerrado sensu stricto soils in both conditions analyzed in this study. C = Control area, protected from fire for more than 20 years; Q = Burned Area, biennially burned. Measures were calculated independently for each replicate and the mean and standard deviation (SD) of the indicators are shown.

Area
OTU observed
Chao1
ACE
Shannon
Coverage (%)

C

Q

Significant differences detected in Welch’s -test ().

The phylogenetic tree constructed with amoA OTUs revealed that most sequences clustered with AOA from the I.1b subgroup, with the exception of Cb 48 (Figure 5). Curiously, although no 16S rRNA gene sequences were classified as I.1a thaumarchaeota, this amoA OTU was affiliated to this subgroup. Excluding Qb 23 and Qc 19, two OTUs representing only sequences from the Burned area that were closely related to Nitrososphaera spp., all other OTUs clustered in clades composed exclusively by uncultured AOA. It is worth highlighting that among these OTUs, Ca 01 contains 85.03% of Control area sequences, being extremely representative of this area’s AOA diversity.

Figure 5: Phylogenetic tree of archaeal amoA gene OTUs (97% gene sequence similarity) obtained from typical Cerrado sensu stricto soils. The number of sequences represented by each OTU is shown in parentheses, with those from the Control area (C) indicated by blue squares and from the Burned area (Q) by red circles. The tree was constructed using FastTree v. 2.1 [54], with 1000 bootstrap tests. Values below 70% are not shown.

4. Discussion

In this study, we evaluated the long-term effects of frequent fires on soil archaeal communities from the Brazilian Cerrado. Fire is a determinant in this biome, with records of significant wildfires and fire utilization by native populations dating back thousands of years [55]. However, increasing human occupation and agricultural activities in the past decades have led to changes in fire frequency and seasonality, and although many alterations in Cerrado’s vegetation structure and composition have been reported [27, 33, 56], information about the long-term effects of periodical fires in microbial communities from soils are still scarce.

The soils investigated in this study were sampled in areas of well-known fire history, where a large experiment on the effects of prescribed burnings has been conducted since 1990. With the exception of ammonium concentrations, no significant differences were detected in most physicochemical properties of soils from the Control and Burned areas (Supplementary Table 1), which is in accordance with studies previously conducted in similar sites of this ecological reserve [33]. Although short-term alterations in soil nutrient availability after fires have been described [22, 57], Cerrado soils are very resilient [22]. Since the aim of this work was to investigate long-term effects of frequent fires, soils from the Burned area were sampled two years after the last burning event, which could explain the similarity observed in most soil properties from both analyzed areas.

Almost all 16S rRNA gene sequences obtained in this study were affiliated to I.1b and I.1c subgroups of Thaumarchaeota. I.1b thaumarchaeotes are dominant in numerous terrestrial habitats [1, 39], and the ability to oxidize ammonia has been identified in many of its members, suggesting an important role of this group in nitrogen cycling in these environments [8, 10, 58]. The I.1c thaumarchaeal lineage was originally described in boreal Finnish forest soils [59], and later 16S rRNA sequences associated to this group have been detected in many other acidic soils [60–62], a characteristic feature of Cerrado [63]. So far, there are no cultured representatives of this group, and although in situ microcosm experiments indicated that these thaumarchaeotes are not associated to ammonia oxidation [64], their ecological roles are still obscure. It is worth mentioning that sequences affiliated to I.1b and I.1c groups have been detected in soils from other Cerrado’s vegetation physiognomies, indicating that these organisms could be common in these soils [15].

A single OTU detected in the Control area, representing a small fraction of 16S rRNA sequences (0.64%), was closely affiliated to a bathyarchaeote clone retrieved from Cerrado lake sediments [65]. Bathyarchaeota (previously known as Miscellaneous Crenarchaeotic Group: MCG) sequences are abundant in anoxic terrestrial and aquatic environments [66, 67], and despite the current lack of cultured representatives, metagenomic data have suggested wide metabolic capabilities in members of this phylum [68], including possible involvement in aromatic compound degradation [69], protein remineralization [70], and methane metabolism [71]. Although bathyarchaeote 16S rRNA gene sequences have been identified in other Brazilian biome soils, such as the Atlantic Forest peatlands [16] and native Amazonian forest [19], members of this phylum have not been previously detected in Cerrado soils. This fact allied to the small number of sequences of this group detected in the present study may suggest a low abundance of Bathyarchaeota in these soils. However, it is important to highlight that sampling depth and period (dry or rainy seasons) may influence these results, since members of this phylum seem to predominate in deeper or highly humid soils [66].

It is worth mentioning that, despite the similarity observed in archaeal richness and diversity inferred by 16S rRNA gene analyses (Table 1), the distribution of thaumarchaeal sequences in Control and Burned areas suggested a higher relative abundance of the I.1c group in the site long protected from fire. A predominance of I.1c thaumarcheota has been commonly reported in more acidic soils [60, 72], but as pH is highly similar in both investigated sites (Supplementary Table 1), it does not seem likely that this factor is correlated to the higher number of I.1c sequences detected in Control soils. Lanzén et al. [73] reported higher abundances of I.1c thaumarchaeota in soils covered by denser vegetation in comparison to recently cleared sites, and considering that frequent fires drastically change Cerrado’s vegetation structure and composition [33], it is possible that denser wood coverage may lead to higher I.1c subgroup abundance in soils from the site long protected from fires. However, quantitative studies are necessary to better understand the possible roles of fires in different Thaumarchaeota subgroup abundance in Cerrado soils.

Interestingly, amoA comparisons showed higher AOA diversity (Table 2) and a greater number of unique OTUs (Figure 3(b)) in soils from the frequently burned site even though analysis of 16S rRNA gene sequences did not reveal such differences (Table 1, Figure 3(a)). Furthermore, although all thaumarchaeal 16S rRNA sequences clustered within I.1b and I.1c subgroups, an amoA OTU identified in Control samples was strongly affiliated with Ca. Nitrosotalea devanaterra, an organism related to the I.1a lineage. These results highlight the importance of using different molecular markers to better describe environmental microbial communities, an aspect that is being increasingly discussed [12, 74, 75]. The detection of many amoA sequences not grouped with any previously cultured AOA, especially in soils from the Control area (Figure 5), suggests the presence of yet unknown AOA in Cerrado soils.

The estimations of greater amoA diversity and OTU numbers in soils from the Burned area are consistent with the detection of a higher number of I.1b 16S rRNA thaumarchaeal sequences in the same condition, given that this group is extensively associated to ammonia oxidation in terrestrial habitats [58, 76, 77]. Moreover, lower ammonium concentrations were detected in soils from the biennially burned area, and considering that lower ammonium availability tends to favor AOA growth in soils [78, 79], this factor could have contributed to the higher amoA diversity found in this site.

Composition differences in soil AOA and AOB (ammonia-oxidizing bacteria) communities were detected in Australian wet sclerophyll forest areas protected from fires or submitted to biennial burning treatments [80], with shifts in archaeal amoA abundance and genotypes, as well as changes in AOB communities, which could result in nitrogen cycling disturbances. However, the differences observed were mainly attributed to alterations detected in soil key parameters between the evaluated sites. Thus, frequent fires in different ecosystems may lead to long-term effects in soil AOA communities through different mechanisms.

There are few studies describing AOA activity and diversity in Cerrado habitats. It has been shown that nitrification rates are generally low in soils from this biome [81] and low AOA and AOB abundances have been reported in soils from a native field of a Cerrado phytophysiognomy denominated Campo sujo [82]. However, it is worth pointing out that, despite the identification of short-term modifications of ammonia oxidizer abundance and soil N dynamics in a managed system within the Cerrado biome [82], there are currently no studies exploring the effects of natural ecological events on nitrifier communities in Cerrado soils, especially concerning AOA. Thus, the present study suggests a possible effect of changes in vegetation cover caused by repeated fires on these communities. More studies are necessary to assert the roles of natural ecological events such as seasonality and different fire regimes on archaeal communities in this biome, as well as their ecological roles.

5. Conclusion

Our study revealed that I.1b and I.1c Thaumarchaeota are dominant in typical Cerrado sensu stricto soils from areas long protected from fire or submitted to frequent fire regime. Although total archaeal community richness and diversity inferred by 16S rRNA genes do not seem to greatly differ between the investigated conditions, amoA gene analyses revealed higher diversity and unique OTU numbers in soils from the biennially burned area. Possible differences in relative abundance of thaumarchaeal subgroups I.1b and I.1c were also observed, with higher numbers of I.1c sequences detected in soils long protected from fires and I.1b sequences in the area frequently burned. Our results suggest that alterations in vegetation structure caused by frequent fires could likely cause long-term effects in the composition of archaeal communities in Cerrado soils, especially regarding potential ammonia oxidizers. Further studies are required to determine the extent to which frequent fires could influence archaeal abundance and potential roles in nitrogen cycling in soils from this unique Brazilian biome.

Data Availability

Sequences obtained in this study were deposited in the GenBank dataset under accession numbers KR828099-KR828564 (16S rRNA gene) and KR828565-KR828812 (amoA gene).

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was supported by grants from the National Council for Scientific and Technological Development (CNPq), Brazilian Federal Agency for Support and Evaluation of Graduate Education (CAPES), and Foundation for Research Support of Distrito Federal (FAPDF). A. Belmok and F. A. C. Lopes acknowledge a fellowship from CAPES.

Supplementary Materials

Supplementary Table 1: physicochemical parameters of the Cerrado sensu stricto soils from the Control area (C) and the Burned area (Q) at the Reserva Ecológica do IBGE, Brasilia, Brazil. Asterisk indicates statistical differences identified in Welch’s -test analysis. Supplementary Table 2: number of high-quality sequences for each 16S rRNA and amoA gene library evaluated in this study. Supplementary Figure 1: rarefaction curves (with 97% sequence similarity) for archaeal 16S rRNA gene sequences from Cerrado sensu stricto soils sampled in the Control (C) and Burned (Q) areas at the Reserva Ecológica do IBGE, Brasilia, Brazil. Supplementary Figure 2: nonmetric multidimensional scaling (NMDS) plot of archaeal communities (inferred by 16S rRNA genes) from soils of Cerrado sensu stricto protected from fires (C) and submitted to biennial fires (Q) at the Reserva Ecológica do IBGE, Brasilia, Brazil. (Supplementary Materials)

References

J. C. Auguet, A. Barberan, and E. O. Casamayor, “Global ecological patterns in uncultured archaea,” The ISME Journal, vol. 4, no. 2, pp. 182–190, 2010.

C. Schleper, G. Jurgens, and M. Jonuscheit, “Genomic studies of uncultivated archaea,” Nature Reviews Microbiology, vol. 3, no. 6, pp. 479–488, 2005.

D. H. Buckley, J. R. Graber, and T. M. Schmidt, “Phylogenetic analysis of nonthermophilic members of the kingdom crenarchaeota and their diversity and abundance in soils,” Applied and Environmental Microbiology, vol. 64, no. 11, pp. 4333–4339, 1998.

T. Ochsenreiter, D. Selezi, A. Quaiser, L. Bonch-Osmolovskaya, and C. Schleper, “Diversity and abundance of Crenarchaeota in terrestrial habitats studied by 16S RNA surveys and real time PCR,” Environmental Microbiology, vol. 5, no. 9, pp. 787–797, 2003.

D. Kemnitz, S. Kolb, and R. Conrad, “High abundance of Crenarchaeota in a temperate acidic forest soil,” FEMS Microbiology Ecology, vol. 60, no. 3, pp. 442–448, 2007.

G. W. Nicol, S. Leininger, C. Schleper, and J. I. Prosser, “The influence of soil pH on the diversity, abundance and transcriptional activity of ammonia oxidizing archaea and bacteria,” Environmental Microbiology, vol. 10, no. 11, pp. 2966–2978, 2008.

K. Zhalnina, P. D. de Quadros, F. A. O. Camargo, and E. W. Triplett, “Drivers of archaeal ammonia-oxidizing communities in soil,” Frontiers in Microbiology, vol. 3, 2012.

L. M. Zhang, P. R. Offre, J. Z. He, D. T. Verhamme, G. W. Nicol, and J. I. Prosser, “Autotrophic ammonia oxidation by soil thaumarchaea,” Proceedings of the National Academy of Sciences of the United States of America, vol. 107, no. 40, pp. 17240–17245, 2010.

R. Bartossek, A. Spang, G. Weidler, A. Lanzen, and C. Schleper, “Metagenomic analysis of ammonia-oxidizing archaea affiliated with the soil group,” Frontiers in Microbiology, vol. 3, 2012.

M. Pester, T. Rattei, S. Flechl et al., “amoA-based consensus phylogeny of ammonia-oxidizing archaea and deep sequencing of amoA genes from soils of four different geographic regions,” Environmental Microbiology, vol. 14, no. 2, pp. 525–539, 2012.

E. F. DeLong, “Everything in moderation: archaea as ‘non-extremophiles’,” Current Opinion in Genetics & Development, vol. 8, no. 6, pp. 649–654, 1998.

C. Brochier-Armanet, B. Boussau, S. Gribaldo, and P. Forterre, “Mesophilic Crenarchaeota: proposal for a third archaeal phylum, the Thaumarchaeota,” Nature Reviews Microbiology, vol. 6, no. 3, pp. 245–252, 2008.

R. J. E. Alves, B. Q. Minh, T. Urich, A. von Haeseler, and C. Schleper, “Unifying the global phylogeny and environmental distribution of ammonia-oxidising archaea based on amoA genes,” Nature Communications, vol. 9, no. 1, article 1517, 2018.

T. Rodrigues, A. Belmok, E. Catão, and C. M. Kyaw, “Archaea in natural and impacted Brazilian environments,” Archaea, vol. 2016, Article ID 1259608, 14 pages, 2016.

E. Catão, A. P. Castro, C. C. Barreto, R. H. Krüger, and C. M. Kyaw, “Diversity of archaea in Brazilian savanna soils,” Archives of Microbiology, vol. 195, no. 7, pp. 507–512, 2013.

R. M. Etto, L. M. Cruz, E. C. Jesus et al., “Prokaryotic communities of acidic peatlands from the Southern Brazilian Atlantic Forest,” Brazilian Journal of Microbiology, vol. 43, no. 2, pp. 661–674, 2012.

J. E. Lima-Perim, E. M. Romagnoli, F. Dini-Andreote, A. Durrer, A. C. F. Dias, and F. D. Andreote, “Linking the composition of bacterial and archaeal communities to characteristics of soil and flora composition in the Atlantic Rainforest,” PLoS One, vol. 11, no. 1, article e0146566, 2016.

R. G. Taketani and S. M. Tsai, “The influence of different land uses on the structure of archaeal communities in Amazonian Anthrosols based on 16S rRNA and amoA genes,” Microbial Ecology, vol. 59, no. 4, pp. 734–743, 2010.

D. D. Tupinambá, M. E. Cantão, O. Y. A. Costa et al., “Archaeal community changes associated with cultivation of Amazon Forest soil with oil palm,” Archaea, vol. 2016, Article ID 3762159, 14 pages, 2016.

N. Myers, R. A. Mittermeier, C. G. Mittermeier, G. A. B. da Fonseca, and J. Kent, “Biodiversity hotspots for conservation priorities,” Nature, vol. 403, no. 6772, pp. 853–858, 2000.

J. F. Ribeiro and B. M. T. Walter, “Fitofisionomias do bioma Cerrado,” in Cerrado ambiente e flora, S. M. Sano and S. P. Almeida, Eds., pp. 89–168, DF: Embrapa-CPAC, Planaltina, Brazil, 1998.

L. M. Coutinho, “Fire in the ecology of the Brazilian Cerrado,” in Fire in the Tropical Biota - Ecosystem Processes and Global Challenges, J. G. Goldammer, Ed., pp. 82–105, Springer, 1990.

G. Eiten, “Vegetação do Cerrado,” in Cerrado: caracterização, ocupação e perspectivas, M. Novaes-Pinto, Ed., pp. 17–73, University of Brasilia Press, Brasília, 1994.

M. F. Simon, R. Grether, L. P. de Queiroz, C. Skema, R. T. Pennington, and C. E. Hughes, “Recent assembly of the Cerrado, a neotropical plant diversity hotspot, by in situ evolution of adaptations to fire,” Proceedings of the National Academy of Sciences of the United States of America, vol. 106, no. 48, pp. 20359–20364, 2009.

M. B. Ramos-Neto and V. R. Pivello, “Lightning fires in a Brazilian savanna national park: rethinking management strategies,” Environmental Management, vol. 26, no. 6, pp. 675–684, 2000.

H. França, M. B. Ramos-Neto, and A. Setzer, “O fogo no Parque Nacional das Emas,” in Série Biodiversidade, vol. 27, pp. 1–140, Ministério do Meio Ambiente, Brasília, Brazil, 2007.

H. S. Miranda, M. N. Sato, W. N. Neto, and F. S. Aires, “Fires in the Cerrado, the Brazilian savanna,” in Tropical Fire Ecology: Climate Change, Land Use and Ecosystem Dynamics, pp. 427–450, Springer-Praxis, 2009.

M. N. Sato and H. S. Miranda, “Mortalidade de plantas lenhosas do Cerrado sensu stricto submetidas a diferentes regimes de queima,” in Impactos de queimadas em áreas de cerrado e restinga, H. S. Miranda, C. H. Saito, and B. F. S. Dias, Eds., pp. 102–111, ECL/University of Brasília, Brasília, Brazil, 1996.

W. A. Hoffmann and O. T. Solbrig, “The role of topkill in the differential response of savanna woody species to fire,” Forest Ecology and Management, vol. 180, no. 1-3, pp. 273–286, 2003.

M. B. Medeiros and H. S. Miranda, “Post-fire resprouting and mortality in Cerrado woody plant species over a three-year period,” Edinburgh Journal of Botany, vol. 65, no. 1, pp. 53–68, 2008.

P. Cirne and H. S. Miranda, “Effects of prescribed fires on the survival and release of seeds of Kielmeyera coriacea (Spr.) Mart. (Clusiaceae) in savannas of Central Brazil,” Brazilian Journal of Plant Physiology, vol. 20, no. 3, pp. 197–204, 2008.

W. A. Hoffmann, “Post-burn reproduction of woody plants in a neotropical savanna: the relative importance of sexual and vegetative reproduction,” Journal of Applied Ecology, vol. 35, no. 3, pp. 422–433, 1998.

A. G. Moreira, “Effects of fire protection on savanna structure in Central Brazil,” Journal of Biogeography, vol. 27, no. 4, pp. 1021–1029, 2000.

G. B. Nardoto, M. M. da Cunha Bustamante, A. S. Pinto, and C. A. Klink, “Nutrient use efficiency at ecosystem and species level in savanna areas of Central Brazil and impacts of fire,” Journal of Tropical Ecology, vol. 22, no. 2, pp. 191–201, 2006.

L. F. Daibes, T. Zupo, F. A. O. Silveira, and A. Fidelis, “A field perspective on effects of fire and temperature fluctuation on Cerrado legume seeds,” Seed Science Research, vol. 27, no. 2, pp. 74–83, 2017.

D. C. Close, N. J. Davidson, D. W. Johnson et al., “Premature decline of Eucalyptus and altered ecosystem processes in the absence of fire in some Australian forests,” The Botanical Review, vol. 75, no. 2, pp. 191–202, 2009.

S. C. Hart, T. H. DeLuca, G. S. Newman, M. D. MacKenzie, and S. I. Boyle, “Post-fire vegetative dynamics as drivers of microbial community structure and function in forest soils,” Forest Ecology and Management, vol. 220, no. 1-3, pp. 166–184, 2005.

B. Stres, T. Danevčič, L. Pal et al., “Influence of temperature and soil water content on bacterial, archaeal and denitrifying microbial communities in drained fen grassland soil microcosms,” FEMS Microbiology Ecology, vol. 66, no. 1, pp. 110–122, 2008.

S. T. Bates, D. Berg-Lyons, J. G. Caporaso, W. A. Walters, R. Knight, and N. Fierer, “Examining the global distribution of dominant archaeal populations in soil,” The ISME Journal, vol. 5, no. 5, pp. 908–917, 2011.

Y. Shi, J. M. Adams, Y. Ni et al., “The biogeography of soil archaeal communities on the eastern Tibetan Plateau,” Scientific Reports, vol. 6, no. 1, article 38893, 2016.

IBGE, Reserva Ecológica do IBGE–Ambiente e Plantas Vasculares, Instituto Brasileiro de Geografia e Estatística, Estudos e Pesquisas, Informação Geográfica 3, Rio de Janeiro, Brazil, 2004.

R Core Team, R: A Language and Environment for Statistical Computing, R Foundation for Statistical Computing, Vienna, Austria, 2018.

E. F. DeLong, “Archaea in coastal marine environments,” Proceedings of the National Academy of Sciences of the United States of America, vol. 89, no. 12, pp. 5685–5689, 1992.

C. A. Francis, K. J. Roberts, J. M. Beman, A. E. Santoro, and B. B. Oakley, “Ubiquity and diversity of ammonia-oxidizing archaea in water columns and sediments of the ocean,” Proceedings of the National Academy of Sciences of the United States of America, vol. 102, no. 41, pp. 14683–14688, 2005.

T. Z. DeSantis, P. Hugenholtz, N. Larsen et al., “Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB,” Applied and Environmental Microbiology, vol. 72, no. 7, pp. 5069–5072, 2006.

P. D. Schloss, S. L. Westcott, T. Ryabin et al., “Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities,” Applied and Environmental Microbiology, vol. 75, no. 23, pp. 7537–7541, 2009.

M. A. Larkin, G. Blackshields, N. P. Brown et al., “Clustal W and Clustal X version 2.0,” Bioinformatics, vol. 23, no. 21, pp. 2947-2948, 2007.

A. Chao, “Nonparametric estimation of the number of classes in a population,” Scandinavian Journal of Statistics, vol. 11, pp. 265–270, 1984.

A. Chao and S. M. Lee, “Estimating the number of classes via sample coverage,” Journal of the American Statistical Association, vol. 87, no. 417, pp. 210–217, 1992.

C. E. Shannon and W. Weaver, The Mathematical Theory of Communication, University of Illinois Press, Urbana, IL, USA, 1949.

P. D. Schloss and J. Handelsman, “Status of the microbial census,” Microbiology and Molecular Biology Reviews, vol. 68, no. 4, pp. 686–691, 2004.

D. H. Parks, G. W. Tyson, P. Hugenholtz, and R. G. Beiko, “STAMP: statistical analysis of taxonomic and functional profiles,” Bioinformatics, vol. 30, no. 21, pp. 3123-3124, 2014.

R. C. Edgar, “MUSCLE: multiple sequence alignment with high accuracy and high throughput,” Nucleic Acids Research, vol. 32, no. 5, pp. 1792–1797, 2004.

M. N. Price, P. S. Dehal, and A. P. Arkin, “FastTree 2 – approximately maximum-likelihood trees for large alignments,” PLoS One, vol. 5, no. 3, article e9490, 2010.

J. Mistry, “Fire in the cerrado (savannas) of Brazil: an ecological review,” Progress in Physical Geography, vol. 22, no. 4, pp. 425–448, 1998.

V. R. Pivello, “The use of fire in the Cerrado and Amazonian rainforests of Brazil: past and present,” Fire Ecology, vol. 7, no. 1, pp. 24–39, 2011.

J. B. Kauffman, D. L. Cummings, and D. E. Ward, “Relationships of fire, biomass and nutrient dynamics along a vegetation gradient in the Brazilian Cerrado,” Journal of Ecology, vol. 82, no. 3, pp. 519–531, 1994.

Y. Wu and R. Conrad, “Ammonia oxidation-dependent growth of group I.1b Thaumarchaeota in acidic red soil microcosms,” FEMS Microbiology Ecology, vol. 89, no. 1, pp. 127–134, 2014.

G. Jurgens, K. Lindström, and A. Saano, “Novel group within the kingdom Crenarchaeota from boreal forest soil,” Applied and Environmental Microbiology, vol. 63, no. 2, pp. 803–805, 1997.

L. E. Lehtovirta, J. I. Prosser, and G. W. Nicol, “Soil pH regulates the abundance and diversity of group 1.1c Crenarchaeota,” FEMS Microbiology Ecology, vol. 70, no. 3, pp. 367–376, 2009.

N. Stopnišek, C. Gubry-Rangin, S. Höfferle, G. W. Nicol, I. Mandič-Mulec, and J. I. Prosser, “Thaumarchaeal ammonia oxidation in an acidic forest peat soil is not influenced by ammonium amendment,” Applied and Environmental Microbiology, vol. 76, no. 22, pp. 7626–7634, 2010.

S. Timonen and M. Bomberg, “Archaea in dry soil environments,” Phytochemistry Reviews, vol. 8, no. 3, pp. 505–518, 2009.

M. Haridasan, “Nutritional adaptations of native plants of the cerrado biome in acid soils,” Brazilian Journal of Plant Physiology, vol. 20, no. 3, pp. 183–195, 2008.

E. B. Weber, L. E. Lehtovirta-Morley, J. I. Prosser, and C. Gubry-Rangin, “Ammonia oxidation is not required for growth of group 1.1c soil Thaumarchaeota,” FEMS Microbiology Ecology, vol. 91, no. 3, 2015.

T. Rodrigues, E. Catão, M. M. C. Bustamante, B. F. Quirino, R. H. Kruger, and C. M. Kyaw, “Seasonal effects in a lake sediment archaeal community of the Brazilian savanna,” Archaea, vol. 2014, Article ID 957145, 9 pages, 2014.

X. Xiang, R. Wang, H. Wang, L. Gong, B. Man, and Y. Xu, “Distribution of Bathyarchaeota communities across different terrestrial settings and their potential ecological functions,” Scientific Reports, vol. 7, no. 1, article 45028, 2017.

T. Yu, Q. Liang, M. Niu, and F. Wang, “High occurrence of Bathyarchaeota (MCG) in the deep-sea sediments of South China Sea quantified using newly designed PCR primers,” Environmental Microbiology Reports, vol. 9, no. 4, pp. 374–382, 2017.

Z. Zhou, J. Pan, F. Wang, J. D. Gu, and M. Li, “Bathyarchaeota: globally distributed metabolic generalists in anoxic environments,” FEMS Microbiology Reviews, vol. 42, no. 5, pp. 639–655, 2018.

J. Meng, J. Xu, D. Qin, Y. He, X. Xiao, and F. Wang, “Genetic and functional properties of uncultivated MCG archaea assessed by metagenome and gene expression analyses,” The ISME Journal, vol. 8, no. 3, pp. 650–659, 2014.

K. G. Lloyd, L. Schreiber, D. G. Petersen et al., “Predominant archaea in marine sediments degrade detrital proteins,” Nature, vol. 496, no. 7444, pp. 215–218, 2013.

P. N. Evans, D. H. Parks, G. L. Chadwick et al., “Methane metabolism in the archaeal phylum Bathyarchaeota revealed by genome-centric metagenomics,” Science, vol. 350, no. 6259, pp. 434–438, 2015.

M. Bomberg, “The elusive boreal forest Thaumarchaeota,” Agronomy, vol. 6, no. 2, p. 36, 2016.

A. Lanzén, L. Epelde, C. Garbisu et al., “The community structures of prokaryotes and fungi in mountain pasture soils are highly correlated and primarily influenced by pH,” Frontiers in Microbiology, vol. 6, article 1321, 2015.

J. R. Freeland, “The importance of molecular markers and primer design when characterizing biodiversity from environmental DNA,” Genome, vol. 60, no. 4, pp. 358–374, 2017.

A. Patwardhan, S. Ray, and A. Roy, “Molecular markers in phylogenetic studies-a review,” Journal of Phylogenetics & Evolutionary Biology, vol. 2, no. 2, p. 131, 2014.

R. Hatzenpichler, “Diversity, physiology, and niche differentiation of ammonia-oxidizing archaea,” Applied and Environmental Microbiology, vol. 78, no. 21, pp. 7501–7510, 2012.

J.-G. Kim, M.-Y. Jung, S.-J. Park et al., “Cultivation of a highly enriched ammonia-oxidizing archaeon of thaumarchaeotal group I.1b from an agricultural soil,” Environmental Microbiology, vol. 14, no. 6, pp. 1528–1543, 2012.

A. Daebeler, P. L. E. Bodelier, M. M. Hefting, and H. J. Laanbroek, “Ammonia-limited conditions cause of Thaumarchaeal dominance in volcanic grassland soil,” FEMS Microbiology Ecology, vol. 91, no. 3, 2015.

L. Zheng, X. Zhao, G. Zhu, W. Yang, C. Xia, and T. Xu, “Occurrence and abundance of ammonia-oxidizing archaea and bacteria from the surface to below the water table, in deep soil, and their contributions to nitrification,” Microbiology, vol. 6, no. 4, p. e00488, 2017.

X. E. Long, C. Chen, Z. Xu, and J. Z. He, “Shifts in the abundance and community structure of soil ammonia oxidizers in a wet sclerophyll forest under long-term prescribed burning,” Science of the Total Environment, vol. 470-471, pp. 578–586, 2014.

G. B. Nardoto and M. M. da Cunha Bustamante, “Effects of fire on soil nitrogen dynamics and microbial biomass in savannas of Central Brazil,” Pesquisa Agropecuária Brasileira, vol. 38, no. 8, pp. 955–962, 2003.

E. C. P. Catão, F. A. C. Lopes, M. R. Rubini, G. B. Nardoto, J. I. Prosser, and R. H. Krüger, “Short-term impact of soybean management on ammonia oxidizers in a Brazilian savanna under restoration as revealed by coupling different techniques,” Biology and Fertility of Soils, vol. 52, no. 3, pp. 401–412, 2016.