Next Article in Journal
Experimental Evaluation of δ2H, δ13C and δ15N Variability in Blood and Feathers of Wild and Captive Birds: Implications for Interspecific Food Web Studies
Next Article in Special Issue
It’s All about the Zone: Spider Assemblages in Different Ecological Zones of Levantine Caves
Previous Article in Journal
Bird Functional Diversity in Agroecosystems and Secondary Forests of the Tropical Andes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

16S rRNA Gene Amplicon Sequencing Data of the Iron Quadrangle Ferruginous Caves (Brazil) Shows the Importance of Conserving This Singular and Threatened Geosystem

by
Camila G. C. Lemes
1,
Morghana M. Villa
2,
Érica B. Felestrino
1,
Luiza O. Perucci
3,
Renata A. B. Assis
1,
Isabella F. Cordeiro
1,
Natasha P. Fonseca
1,
Lara C. C. Guerra
2,
Washington L. Caneschi
1,
Lauro Â. G. Moraes
3,
Flávio F. do Carmo
4,
Luciana H. Y. Kamino
4,
Pedro N. C. Vale
5,
Suzana E. S. Guima
6,
João C. Setubal
6,
André A. R. Salgado
5 and
Leandro M. Moreira
1,2,*
1
Programa de Pós-Graduação em Biotecnologia, Núcleo de Pesquisas em Ciências Biológicas, Universidade Federal de Ouro Preto, Ouro Preto 35400-000, MG, Brazil
2
Departamento de Ciências Biológicas, Instituto de Ciências Exatas e Biológicas, Universidade Federal de Ouro Preto, Ouro Preto 35400-000, MG, Brazil
3
Laboratórios Multiusuários de Genômica e Bioinformática, Universidade Federal de Ouro Preto, Ouro Preto 35400-000, MG, Brazil
4
Instituto Prístino, Belo Horizonte 30642-180, MG, Brazil
5
Departamento de Geografia, Instituto de Geociências, Universidade Federal de Minas Gerais, Belo Horizonte 31270-901, MG, Brazil
6
Departamento de Bioquímica, Instituto de Química, Universidade de São Paulo, São Paulo 05508-000, SP, Brazil
*
Author to whom correspondence should be addressed.
Submission received: 8 September 2021 / Revised: 4 October 2021 / Accepted: 4 October 2021 / Published: 11 October 2021
(This article belongs to the Special Issue Biodiversity in Karst Landscapes)

Abstract

:
The Iron Quadrangle (IQ) is one of the main iron ore producing regions of the world. The exploitation of its reserves jeopardizes the high biological endemism associated with this region. This work aimed to understand the diversity and bacterial potential associated with IQ caves. Floor and ceiling samples of seven ferruginous caves and one quartzite cave were collected, and their microbial relative abundance and diversity were established by 16S rRNA gene amplicon sequencing data. The results showed that ferruginous caves present higher microbial abundance and greater microbial diversity compared to the quartzite cave. Many species belonging to genera found in these caves, such as Pseudonocardia and Streptacidiphilus, are known to produce biomolecules of biotechnological interest as macrolides and polyketides. Moreover, comparative analysis of microbial diversity and metabolic potential in a biofilm in pendant microfeature revealed that the microbiota associated with this structure is more similar to the floor rather than ceiling samples, with the presence of genera that may participate in the genesis of these cavities, for instance, Ferrovum, Geobacter, and Sideroxydans. These results provide the first glimpse of the microbial life in these environments and emphasize the need of conservation programs for these areas, which are under intense anthropogenic exploration.

1. Introduction

For many researchers, karst is a type of landscape that only occurs in carbonates. However, during the last decades several karst landscapes in non-carbonate rocks have been identified in South America [1,2,3,4,5]. Although most of these landscapes are located on siliciclastic rocks, karstic geoforms, mainly caves, have also been found in iron ore, such as banded iron formation and canga (superficial weathering product). Many of these iron ore caves are in the Iron Quadrangle (IQ) in southeastern Brazil [6]. In IQ, although most caves are small, they are high in density and several have strong evidence that the dissolution processes were important to their morphogenesis [7,8,9,10]. This allows to classify the region as karstic or, at least, pseudo-karstic.
The IQ is located in the center-south region of the State of Minas Gerais, Brazil, and corresponds to an area of approximately 7200 km2 [11] (Figure 1) that share two biodiversity hotspots: the Cerrado (Brazilian Savanna) and the Atlantic Forest [12] (Supplementary Figure S1). The ironstone ranges, which are geologically formed by banded iron formation (BIF) of itabirites with very variable iron concentration of Archean and Paleoproterozoic ages, occur in this region. The iron concentration in itabirites is variable as this rock has two bands: one composed by hematite (Fe2O3) and other by quartz (SiO2) and the thickness of each one of these bands naturally varies. Duricrusts, known as cangas, formed by minerals resulting from the natural alteration of itabirite, such as limonite and goethite (FeO(OH)) [10,13], cover the upper reaches of ironstone ranges and reach up to 30 m in thickness [13,14].
The ironstone ranges are subject to abiotic stressors such as acid substrates and possess anomalous metal contents. These abiotic stressors are due to the natural weathering process of BIF’s that transforms the hematite, present in the itabirite, into limonitic capping and in laterized clay deposits, which are rich in goethite. The associated soils are extremely oligotrophic and depleted of exchangeable nutrients, such as calcium, potassium, magnesium, and total nitrogen, resulting in an environment with intense ecological filters, and harboring a high number of endemic species that are rare and endangered [15,16,17].
The IQ also shelters hundreds of ferruginous caves and is one of the regions with the highest occurrence of these geoforms in the word. In general, although they take a long time to form, ferruginous caves are small in size, as they rarely exceed 100 m in horizontal projection, with an average of 61 m2 of area and 81.3 m3 of volume and 3 m of unevenness in the cave elevation along its horizontal projection [18].
The small diameter entrances of these caves are almost always located at the edges of the canga as a result of erosive processes or, less frequently, are small vertical openings caused by collapse of the canga mantle [7]. Although small, they have a very diversified geological composition, found only in BIFs, in cangas, or in contact with both; they generally have a ceiling composed by the canga and a floor formed by the BIF [10]. Despite their distinct composition, their floors are usually filled with ferruginous sediments (rich in hematite, limonite, and goethite) of granulometries ranging from clay to boulders [18]. Its speleogenesis may be related with many processes [19], but the most recent paper defends a hypogenic genesis related to late vadose stages where the erosion only exhumes the cave [20]. Most of the IQ caves, despite not having an evident entrance, are not completely isolated from the surface. As the canga is very porous and the itabirite is fractured, flows of superficial waters enable penetration of micro particles inside the caves. The interiors of the caves present considerable aphotic extensions, and most of these have stable environmental conditions such as temperature, pH, and humidity, providing a highly selective environment for life [21]. Although the biofilm in pendant microfeature in ferruginous cavities may be rare [22,23], these structures are defined as secondary mineral deposits originating from the processes of circulation and discharge of water or from biological activity [24].
Despite its biological and speleological importance, the IQ is one of the natural areas in Brazil most endangered with a loss and degradation of natural ecosystems, due to the intense mining activity and urban expansion. The region concentrates one of the world’s largest open pit complexes [25,26] (Supplementary Figure S1). In BIF areas, there is a defined relationship between the occurrence of iron ore and high-grade ore caves, as a result of which, the full exploitation of iron resources will lead to the loss of caves [27]. The industrial waste generated during mineral processing has several impacts: pollution of groundwater, rivers, and soil by metals and trace elements (e.g., Cr, Fe, Mn, and Zn); alterations to water turbidity; as well as silting of stretches of streams located downstream of the dams. In addition, thousands of tons of chemicals used for the treatment and processing of minerals (such as caustic soda, sodium hypochlorite, ether, and hydrogen peroxide) are capable of altering the natural conditions of water bodies, soil, biota, ecosystems, and human health [28]. Additionally, in the last years, discharges of tailing dams have placed the IQ, causing disasters as well as irreparable environmental damage [29,30]. However, there are still no studies that reveal the environmental impacts and damages caused by scale mining industry. Therefore, it is likely that such influences have been source of losses of microbial biodiversity and cave invertebrates.
Studies reporting the presence and importance of the microbiota in ferruginous caves are still incipient. Parker and colleagues have made significant contributions to the importance of the microbiota in the biogenesis of ferruginous cavities [23,31]. According to these authors, the action of these bacteria when reducing iron III to iron II allows their solubilization and subsequent leaching through groundwater, which, taking advantage of the system of rock failures and fractures and the microporosities of the canga, penetrates and seeps into the interior the cave. Thus, expanding knowledge about the diversity and importance of the microbiota associated with these cavities, especially in the IQ region, provide important characterization of the ecosystems within these unique environments [32] (Supplementary Figure S1).
In this context, this work aimed to investigate the prokaryotic diversity associated with the natural subterranean cavities in ferruginous rocks by 16S rRNA gene amplicon sequencing, and to predict the characteristics and metabolic potentials of the microorganisms identified in this environment. The reported data will contribute to a broader initiative that aims to establish conservation areas within the IQ in the future. In addition, the data may support the conservation in situ of iron caves into newly protect areas or to stimulate studies in few caves located in conservation areas, such as the Gandarela National Park, the main preservation area in the entire region, and the Natural Monument of Serra da Moeda [33].

2. Materials and Methods

2.1. Collection Sites

The determining factors in the selection of caves included geographic location, variation in altitude, and degree of human interference (Table 1 and Figure 1). Selection of the quartzitic cave used as an experimental control took into account the geographic proximity and lithology similarity to the IQ. According to the literature, caves are stable environments, where changes in temperature and humidity are lesser than those in the atmosphere, and therefore the internal temperature range between 11 °C and 22 °C and humidity between 48% and 98% in the IQ ferruginous caves throughout the year, while pH ranges from 3.5 to 4.9 [34]. The samples were collected at points close to the cavity entrances, therefore in a borderline region between the photic and aphotic zones (when possible), and in different aphotic regions depending on the dimensions of the cavities.
Samples of 50 g of solid material from the floor, ceiling, and biofilm hanging from the ceiling were collected from seven caves associated with canga outcrops and a quartzite cave (experimental control), using previously sterilized materials such as spatulas and 50-mL Falcon® tubes (Life Sciences, Corning, NY, USA).

2.2. Total Genomic DNA Extraction

Total microbial DNA was extracted from 0.25 g of each sample using the PowerSoil DNA Mobio™ kit (GeneWorks, Adelaide, South Australia, Australia), following the protocol recommended by the manufacturer. After extraction, the DNA samples were quantified and their purity was analyzed using NanoDrop™ (Thermo Fisher Scientific, Waltham, MA, USA). The sample quality was verified by electrophoresis on 0.8% agarose gel. No blank samples were used as a contamination control associated with the sequencing kit reagents [35].

2.3. Partial Amplification of the 16S Ribosomal Gene

To amplify 16S ribosomal genes of the V4-V5 region, 515F and 926R oligonucleotides were used [36]. These oligonucleotides were designed to include the adapter sequences used in the Ion Torrent sequencing protocol and were linked to a 10 nt barcode for each sample [36]. The complete list of oligonucleotides used in this study can be found in Supplementary Table S1. The PCR was prepared in a final volume of 30 μL containing 0.38 mM dNTP, 2.5 mM MgCl2, and 0.15 μM of each oligonucleotide, 20 ng DNA, 0.08 U of Platinum™ Taq DNA polymerase (Invitrogen, Waltham, MA, USA), and 10× Buffer. Amplification was performed in a thermocycler Biocycler® (Applied Biosystems, Waltham, MA, USA) with an initial cycle of 5 min at 94 °C, followed by 35 cycles of 30 s at 94 °C, 1 min at 57 °C, and 1 min at 72 °C, with a final cycle of 5 min at 72 °C. The amplicons were analyzed on an 0.8% agarose gel. Once quality was checked, 22 μL of the final PCR product was purified using the GFX PCR Kit™ and Gel Band Purification™ (GE Healthcare, Chicago, IL, USA) Kit, following manufacturer’s recommendations. The quality of the purified product was again measured on 0.8% agarose gel.

2.4. 16S rRNA Amplicon Sequencing Using the Ion Torrent Platform

The DNA library was quantified using Qubit dsDNA HS™ (High Sensitivity) using the Qubit 2.0™ fluorometer (Life technologies, Corning, NY, USA) and the samples were then normalized to 70 pM. For library assembly, 20 μL of each sample was placed in a single 1.5 mL tube and homogenized. Then, 25 μL of this DNA pool was loaded onto a 318 v2™ chip (Life technologies, Corning, NY, USA), and the samples were sequenced using the Ion PGM™ platform (Life technologies, Corning, NY, USA) following the manufacturer’s recommendations. After sequencing, the polyclonal sequences were filtered using the software Torrent Suite™, version 5.0.5, and the data obtained were exported as FastQ files. The homopolymers, stretches of the same nucleotide sequence, are also detected at very high accuracy by the Ion technology [37]. A 5-mer is currently called with greater 97.5% per base accuracy according to the IonTorrent Platform, although several studies report that there should be a concern with the theme [37,38,39,40].

2.5. 16S rRNA Amplicon Sequencing Data Analyses

The data were processed using the clustering-based pipeline adapted from the Brazilian Microbiome Project (BMP) [41]. In the pre-processing step, reads were filtered using Prinseq [42]. Low quality reads of less than 100 bp, and with an average quality (Q) ≤ 20 were discarded. Filtered reads were dereplicated using the USEARCH v.8.1 script [43], performing clustering of abundant sequences and discarding singletons (single unique reads, without similarity). For rarefaction analysis, we adopted the sequencing depth of 28,252 reads for all samples, which was the number of reads of the smallest sample. After rarefaction, the reads were grouped into OTUs (Operational Taxonomic Units) following the UPARSE method [44], where the identity of sequences corresponding to the same OUT was ≥97%. For each OTU generated, a taxonomy was assigned using the RDP 11.5 (Ribosomal Database Project) classifier [45] available in QIIME 2 v.2021.2 [46]. Additionally, all sequences associated with chloroplasts were manually excluded from the analyses. A table was then obtained in the BIOM (The Biological Observation Matrix) format containing the corresponding read counts for the OTUs.
Alignment of OTUs was processed by PyNAST [47], available in QIIME 2. A phylogenetic tree was constructed using FastTree [48], also available in QIIME 2. Diversity analyses were performed from the BIOM file and the phylogenetic tree to generate the results of alpha-diversity, beta-diversity, frequency tables, and taxonomic classification charts. For the rarefaction plot we used the alpha pipeline from the diversity plugin available on QIIME 2 v.2021.2 [46]. To obtain the observed features for all samples, the following groups were compared: cave vs. cave (regardless of coming from ceiling or floor), ceiling vs. floor (regardless of the cave), and ceiling vs. floor separated by caves. The Kruskal–Wallis test [49] was used for these comparisons (alpha-group-significance visualizer from the diversity plugin).
The Phyloseq package v.1.34.0 [50] was used for alpha and beta diversity analysis. From the BIOM table imported into R software, alpha diversity metrics such as observed features, Shannon [51], Simpson [52], and Chao1 [53] were calculated and plotted using the function plot_richness from phyloseq. For ordination plots of beta diversity metrics, sampling counts from the imported BIOM were first transformed to even sampling depth with transform_sample_counts function (formula: 1 × 106 ∙ x/sum(x), where x is the number of OTU counts in a sample, and sum(x) is the sum of all OTUs counts in a sample). Ordinate function from phyloseq was then used to calculate both unweighted and weighted Unifrac distance [54] and the graphs of the Principal Coordinate Analysis (PCoA) were generated from these distances [55]. Ellipses were computed for the ordination plot with stat_ellipse function from ggplot2 v.3.3.3 [56] considering a multivariate t-distribution with 0.95 level. For beta diversity statistical analysis, we used the diversity plugin from QIIME2. Unweighted Unifrac and weighted Unifrac [54,55] distances were calculated between samples and added into distance matrices (beta and beta-phylogenetic pipeline). ANOSIM and PERMANOVA tests [57] were used with 999 permutations in order to compare the ceiling group vs. the floor group and one cave vs. another cave (beta-group-significance visualizer with pairwise option from the diversity plugin). Biofilm in pendant microfeature sample was not included for these tests, as the material associated with the biofilm structure only allowed for one sample to be collected.
All Metagenome sequences were deposited in the Sequence Read Archive under Bioproject PRJNA575049. The biosample IDs (SAMN12876970 to SAMN12876998) and their corresponding cave names are present in the Supplementary Table S1.

2.6. Functional Metabolism Prediction from 16S rRNA Amplicon Sequences

To perform functional metabolic predictions, we used the program PICRUSt2 [58]. A Nearest Sequenced Taxon Index (weighted NSTI) of score 0.03 was used as default. Functional classification was based on KEGG [59], which classifies the metabolic pathways into three possible levels, as detailed below. The analysis was performed for Level 1 categories: Metabolism, Cell Processes, Cell Process/Signaling, Environmental Information Processing, and their respective subcategories (levels 2 and 3). Statistical analyses were performed using STAMP (STatistical Analysis of Metagenomic Profiles) software 2.1.3 [60], for multiple groups ANOVA was used followed by post hoc Tukey–Kramer test (p < 0.05) and for two groups Welch’s t-test was used (p < 0.05).

3. Results

3.1. Sample Quality and the Number of Sequences Obtained

Thirty samples were collected from seven ferruginous and one quartzite caves: 16 samples from the floor, 13 from the ceiling, and one from a biofilm in pendant microfeature (Table 1 and Figure 1). Sample 18 was disregarded from the analysis for technical problems with hybridization of the specific primer-barcode. No blank samples were used as a contamination control. More than 8 million reads were generated, with fragments varying between 100 and 450 bp (featuring a bimodal distribution in the lengths of 200 bp and another in 400 bp), totaling 76% of the total capacity utilization of the chip (Supplementary Table S1). After removal of the smaller-size sequences (between 100 and 300 bp, thus eliminating the reads that made up the first peak of the read distribution), 6 million reads were considered valid (74% of the chip capacity), distributed among 29 samples with coverage varying between 84,000 and 358,000 reads per sample (Supplementary Table S1).

3.2. Alpha Diversity Characterization

Rarefaction curves were determined for each cave based on Observed_OTUs (Figure 2a). Cave 1, followed by caves 5, 6, and 4, were those that presented the highest number of observed OTUs, whereas cave 7 (quartzite), together with cave 8 and 3, presented the fewest. The three representative caves of Chapada de Canga (1, 2, and 3) presented distinct rarefaction curves. Despite the observed differences, there was no statistical difference between any of the rarefaction curves.
The various measures of alpha diversity for a given cave yielded similar results (Figure 2b and Supplementary Figure S2A,B). General observations that can be made are that floor samples tend to be more diverse than ceiling samples, although there are no statistical differences; and that there was substantial variation in alpha diversity in samples from the same cave, especially for caves 2 and 8.

3.3. Beta-Diversity Characterization

PCoA (unweighted and weighted Unifrac, Figure 2c) analysis of samples showed that floor samples could be clustered together more tightly than ceiling samples, with statistical difference between ceiling and floor (Supplementary Table S2). Outliers in the case of floor samples were 25 (quartzite cave) and 27 (cave 8). Note also that its ceiling sample (26) was placed next to its floor sample, indicating that these samples have similar composition. This was not the case for cave 8, as its ceiling samples (29 and 30) were placed far from the floor sample.

3.4. Comparison of Prokaryotic Community Structures among Caves

The composition of the communities in each cave was investigated. It was found that between 14 and 26 phyla associated with 101 to 183 genera can be found per investigated sample, with a number of non-redundant phyla varying between 23 and 27, and non-redundant genera ranging from 161 to 216 per investigated cave (Supplementary Tables S3 and S4). In the ferruginous caves, the most abundant phyla were Acidobacteria (43.6%), Chloroflexi (15.4%), Proteobacteria (14.8%), Actinobacteria (4%), Nitrospirae (2.0%), and Firmicutes (1.8%), which together represented 86% of the entire prokaryotic community (Figure 3). The most representative classes in the ferruginous caves were Alphaproteobacteria (16.4%), Betaproteobacteria (15.8%), and Actinobacteria (15.4%), whereas Gammaproteobacteria (16%), Betaproteobacteria (15%), and Actinobacteria (15%) dominated in the quartzite cave (Figure 4a). The genera with greatest relative abundance in the quartzite cave were Aciditerrimonas and Leucobacter (Actinobacteria), Bryobacter, and GP6 (Acidobacteria), Alcaligenes (Betaproteobacteria), Jahnella (Deltaproteobacteria), Coxiella, Rhodonobacter, and Xanthomonas (Gammaproteobacteria), Gemmatimonas (Gemmatimonadetes), and Nitrospira (Nitrospira). In contrast, the most abundant genera in the ferruginous caves were Blastopirellula (Planctomycetia), WPS1 (WPS), and Beijerinckia (Alphaproteobacteria) (Supplementary Tables S3 and S4 and Figure 4b).

3.5. Metabolic Analysis

A general metabolic prediction analysis identified similarities in metabolic potential between the ferruginous microbiota and the quartzite cave microbiota (Figure 5A). However, some categories of biotechnological interest, when analyzed separately showed statistical differences between the quartzite cave and ferruginous cave samples (Figure 5B). Within all 61 metabolic pathways associated with the level Metabolism (level 1), four subcategories on level 3 are highlighted: novobiocin biosynthesis (Biosynthesis and Biodegradation of secondary metabolites), nitrogen metabolism (Energy metabolism), biosynthesis of 12, 14, and 16 membered macrolides, and biosynthesis of type II polyketide products (Metabolism of Terpenoids and Polyketides). Regarding categories associated with (i) Cellular Processes/Signaling and (ii) Environmental Information and Processing (all level 1), we highlight the level 3 subcategories involved in Signal Transduction Mechanisms (Signal Transduction) and Secretion System (Membrane Transport) (Figure 5B). The samples from the ferruginous caves showed a higher proportion of sequences for the categories signal transduction mechanisms, nitrogen metabolism, novobiocin biosynthesis, and secretion system. For the pathways: biosynthesis of 12, 14, and 16 membered macrolides and biosynthesis of type II polyketide products the highest proportions were from the quartzite cave.

3.6. Ferruginous Caves Associated with Chapada de Canga

Three of the eight caves investigated in this study are located in the Chapada de Canga region (Figure 6a). Despite their geographical proximity, the observed microbiological richness differed between these three caves. Although no statistical difference was observed while cave 1 presented highest richness (Figure 2a), cave 3 presented the least richness and still showed the highest diversity among the caves in this region (with 32 specific OTUs) (Figure 6b). Considering only the samples of the cave floor, the highest diversity (at the genus level) was observed in cave 2, with 27 specific OTUs (2.1% of all OTUs) (Figure 6c). Specific OTUs are those identified only at that location when compared to the other Chapada de Canga caves in this study. The samples from cave 3, with 38 specific OTUs (22.5% of all OTUs) (Figure 6d), showed the highest diversity among the ceiling samples. Among the specific OTUs of the ceiling of each cave, cave 2 stands out due to its high relative abundance of Actinobacteria (30%). Less abundances of Actinobacteria, 0 and 19%, were observed in the ceiling samples of cave 1 and 3, respectively.
Cave 2 was analyzed in detail as it is the only one to present a biofilm hanging from the ceiling. When comparing the floor, ceiling, and biofilm samples obtained in cave 2, it was possible to identify different number of specific OTUs. A greater number of specific genera (24) was observed in the floor samples, compared with 10 and 9 in the biofilm and ceiling, respectively (Figure 7a). Interestingly, the number of genera shared exclusively between the floor and biofilm (26) and between the floor and ceiling (25) was much larger compared to the number of genera shared between the ceiling and biofilm (1), which are physically connected. In order to make this correlation more evident between the microbiota of the floor and biofilm, functional prediction of microbiota metabolism involving the three substrates was analyzed. Among the 24 categories investigated, 21 presented direct correspondence between the floor and biofilm, with emphasis on the pathways of secondary compounds and cellular processes (Figure 7b).

3.7. Inference of Ecological Relations in Ferruginous Caves

A systematic analysis of metabolic capacity of each genus identified in cave 2 (Supplementary Table S5) was performed to establish a model that indicates the importance of the microbiota for maintaining ecological relations within the ferruginous caves.

4. Discussion

Ferruginous Caves in IQ

The physical environmental heterogeneity inside the iron cavities is directly reflected in their microbiological biodiversity, and provides unique habitats with pronounced differences in resources, with some being wet and others with and without light [23,31]. In addition, the fine sediments injected into the floor of the caves through microcracks and microchannels drainage are formed in an anaerobic environment. Finally, as the dissolution of silica and iron is very slow [61] the iron-developed caves of the IQ have taken a long time to form and have been isolated from the surface [20]. The chemical and physical properties of the caves indicate that the microbiomes inhabiting their interiors have to possess mechanisms to survive in such peculiar environments. The characteristics and aspects of the environment in IQ form this singular environment, which adds even greater scientific interest as it is a key area for the conservation efforts focused on protecting microbiological diversity. These conservation efforts also favor the investigations about bacteria associated with rare plants of ferruginous rocky fields as a reservoir of the microbiota of biotechnological and ecological [62,63] and, in addition, protect the plant biodiversity, which includes hundreds of endemic [64,65,66] and rare species [65] (Supplementary Figure S1). Even with the imminent risk of loss of these environments, studies involving such cavities are still incipient, and the microbiota associated with this geosystem are understudied. Such studies can contribute to conservation efforts, as well as broaden our understanding of the rehabilitation capacity of areas degraded by mineral extraction, as demonstrated by a recent work involving the mined areas of the ferruginous fields of the Carajás region, also in Brazil [67].
To our knowledge, this work is the first report of the 16S rRNA gene amplicon sequencing data analyses associated with seven ferruginous caves located in the Brazilian IQ. It is important to highlight that, although no control sample has been incorporated in this work (blank sample), possible contamination from water and other reagents used in sequencing [35], if they have occurred, probably interfered very little and equally in the results obtained, as all samples were processed and sequenced using the same inputs. The results obtained showed that when compared to a quartzite cave, ferruginous caves presented high microbial richness (Figure 2), even though a significant statistical difference was not observed. Microorganisms have a biogeochemical function, including biogenesis processes in the iron caves themselves, active participation in iron oxi-reduction cycles, generation of sulfate and phosphate, and a role in the formation of microchannels and speleothems [68,69,70]. Energy flux within ferruginous caves is not fully understood. In general, microorganisms in cave systems can play key roles in maintaining subterranean biodiversity, allowing nutrient cycling [71] and potentially the survival of plants established above these cavities, as well as of animals. Although the roles of microorganisms and their ecological interactions need further investigation, the fact that extensive branching root systems are established within the structures of the cavities indicates complex ecological functions (Supplementary Figure S4).
Unlike all other cavities, cave 8 showed signs of human interference, such as traces of fire and litter (PET bottles, cigarette packs and plastic bags). Campfires were found in the vicinity of the cave and also about three meters inside it. A comparison between samples revealed that the microbiota in samples collected closest to the anthropogenic intervention zone in cave 8 differed from the microbiota in areas farther from it. Such a distinct spatial difference was not observed in other caves (where more than one sampling point was tested).
It has been demonstrated that caves accessible to the human visitors and the presence of exogenous organic matter affect diversity and richness when compared to undisturbed areas [72]. In this context, the slow genesis and seclusion of the IQ caves have been essential factors in defining and sustaining their microbial diversity. As observed in cave 8, even minor human interventions cause disturbances, indicating a low resilience of these subterranean systems [17]. This further highlights the mineral extraction activities of the IQ region as an important factor contributing to the loss of potential uncharted genetic heritage, as mining of the ores that constitute the cavities leads to a complete destruction of the caves.
Despite similar geomorphological composition to that of cave 8 (canga ceiling and itabirite floor), much greater richness in terms of observed OTUs was observed in caves 1, 4, 5, and 6. As cave 4 has small dimensions of its conduits, only one ceiling sample was analyzed. It presented a diversity of microorganisms very similar to other caves that had more points sampled. Although different points of cave 5 and 6 were sampled, metabolic diversity predictions were quite homogeneous among the samples investigated, demonstrating that the composition of microbial communities is quite conserved. In other words, the distance of the sample collection point from the entrance of the cave did not affect the microbial diversity, which supports our abovementioned speculation about the effect of anthropogenic interference on the microbiota in cave 8.
Beta-diversity analyses, especially using the weighted Unifrac metric, showed that the majority of floor samples clustered together; ceiling samples on the other hand tended to be more scattered. Other factors, such as altitude or location (within or outside protected areas), did not seem to affect the clustering. In contrast to ferruginous caves, floor and ceiling samples from the quartzite cave were quite similar to each other.
It was observed that the phylum Acidobacteria dominated in ferruginous caves, while the phylum Proteobacteria was the most abundant in other kinds of caves [73,74,75,76]. This observation can be explained by differences in geomorphological and geochemical characteristics, genesis and/or ecological interactions [10]. Additionally, Acidobacteria are associated with low pH of the environment, which is in agreement with the reference values (pH below 5.09) observed in ferruginous caves. The most abundant genus in the canga caves was Blastopirellula, whereas the genera that present the greatest abundance in the quartzite cave included Aciditerrimonas, Bryobacter, Rhodanobacter, and Nitrospira.
From the abundance of identified taxa, it was possible to trace a metabolic profile and to identify traits in the microbiota associated with ferruginous caves that could have use in biotechnologies. Although it is known that only a small proportion of microorganisms can be cultivated [77], thus allowing a targeted analysis of the metabolites produced, it is already demonstrated that it is possible to identify and isolate new biomolecules of wide interest without dependence on culture medium [77,78]. Such techniques could be used in the future to analyze samples from the ferruginous caves, in order to support or refute the conclusions of this study that are described below.
Out of the 328 metabolic pathways identified at level 3, six were selected based on their biological importance and biotechnological interest. Except biosynthesis of 12, 14, and 16 membered macrolides and biosynthesis of type II polyketide products categories, the quartzite cave samples contained reduced proportions of sequences associated with metabolic predictions in all other selected categories.
Macrolides and novobiocin are antibiotics naturally produced by some species belonging to the Staphylococcus and Streptomyces genera, as well as Cyanobacteria [79,80,81], and are of wide interest to the pharmaceutical industry [82,83,84]. However, polyketides are an important and structurally diverse group of secondary prokaryotic metabolites that have several biological functions, but have also been widely reported as molecules with antibiotic and immunosuppressive activities [85]. Considering that the main bacterial genera that produce these compounds belong to the phylum Actinobacteria, which was detected in higher proportions in the caves 8, 7, and 2 (averages 16.6%, 15.4 % and 14.2% of the total OTUs of the floor, ceiling, and biofilm samples from each cave, respectively), the natural cavities could become a promising source of isolates that might find use in the production of new drugs. In addition to potential biotechnological applications, Acidobacteria play important ecological roles, such as plant biomass degradation [86], and plant growth promotion [87].
Adaptive responses to the environment are dependent on signal transduction mechanisms, which can activation secretion systems in prokaryotes [88]. These systems, in turn, can secrete toxins into the environment. Therefore, these categories are related to one another, and may be related to the ecological complexity established in these environments [89].
In caves of other lithologies, the importance of nitrogen fixation and the nitrogen cycle have already been reported as fundamental to maintain the entire ecosystem [90]. As Desai and collaborators (2013) have suggested, diazotrophy might be widespread in chemosynthetic communities in dark caves. Elevated abundance of chemosynthetic organisms (including Aquisphaera, Ferrovum, Georgfuchsia, Methylocystis, and Propionibacterium) were determined by the metagenomics analysis, indicating the occurrence of fixation of atmospheric nitrogen in the ferruginous caves.
The Chapada de Canga, located in the east of the IQ, extends for about 3900 ha [91] and has a unique biological, archaeological, speleological, geological, and landscape heritage [91]. Despite these characteristics, this region is surrounded by areas of mining and tailing dams, one of which ruptured in 2015, causing irreparable environmental, economic, and social harm [30,92,93]. Three out of the seven caves selected for this study are located in the region affected by the dam’s spill. Our work has shown that each of these three caves presented completely different diversity, despite similar morphology, dimensions, and proximity to each other. It was found that, although the richness of cave 3 is the smallest of all analyzed in relation to the abundance of OTUs, its diversity is the largest among the Chapada de Canga. Another feature of the caves of Chapada de Canga is related to the abundance of specific phyla. The phylum Actinobacteria includes bacteria that are widely studied as they present potential antagonistic activity to human, animal, and plant pathogens [94] which contributes to give importance to the caves of Chapada de Canga, where abundances of this phylum were identified. Bacteria belonging to this phylum have been detected in several environments, including extreme habitats [95,96,97], and have been characterized as potential producers of biomolecules with anticancer and antimalaria action among other pharmacological interests [98]. These findings are in agreement with the functional metabolic predictions in categories shown, and indicate that cavities could serve as a source of organisms for a new biomolecule production.
A biofilm hanging from the ceiling was identified in cave 2 and sampled for metagenomic analysis. The mineralogy of speleothems is diverse, including, e.g., oxides and hydroxides of iron, phosphates, and sulfates [99], which affect the composition of microbiomes that, on the other hand, play a fundamental role in the mineral deposit formation [100]. The prokaryotic community found in the biofilm was more similar to the floor samples than the ceiling samples, and included genera unique to each of the niches.
This high correspondence between biofilm and floor samples was also seen in functional metabolic predictions. In almost all pathways of secondary compound metabolism and for most of the cell pathways investigated, there was a direct correspondence between the floor and biofilm, except for the nitrotoluene degradation pathways, novobiocin biosynthesis, and bacterial toxins in this context. This result dynamically highlights two important biological events commonly observed in cave environments: the origin of the microbiota associated with the biofilm as a product of the weathering of the soil and the canga above the cavity, and the biofilm as a physical chemistry structure for transferring the associated microbiota to the floor [101].
Results showed that the genera associated with biofilm hanging from the ceiling have a high capacity for oxidation and reduction in iron, phenomena directly associated with the genesis process of the cavities [20,68]. Among the OTUs unique to the biofilm (that were not present in the floor and ceiling samples from the same cave), the most abundant genera were Ferrovum, Geobacter, and Sideroxydans, described as iron oxidizers [102,103,104]. Ferrovum and Sideroxydans species are often associated with caves [105], being directly involved with mineral precipitation and thus formation of biofilm [106,107]. In addition, Geobacter species have been described as capable of reducing toxic and radioactive metals [108,109], and they therefore play an important role in bioremediation of contaminated sites [110,111].
Methylocapsa, Proprionibacterium, Streptacidiphilus, and Pseudonocardia were identified as unique genera in cave 2 ceiling samples. Methylocapsa species have already been described as capable of fixing atmospheric nitrogen [112], whereas bacteria belonging to the genera Proprionibacterium, Streptacidiphilus, and Pseudonocardia have been described as species with high potential for antifungal activity [113,114,115,116]. Additionally, all these genera are directly or indirectly associated with nutrient cycling (e.g., recalcitrant organic compound and organic matter; Atashgahi et al. [117]) maintenance of biogeochemical cycles (carbon; [118] and nitrogen; [119]), and other metabolic functions (fixation, reduction, and oxidation of metals). These organisms are therefore suitable for several biotechnological applications, such as rhizoremediation [120], effluent treatment [121], promotion of plant growth [120], plant pathogen control [122], and metal biotransformation [123]. In the floor samples, the genera Burkholderia, Rhodanobacter, and Variovorax dominated. Some species of Burkholderia are of broad agricultural and industrial interest [124,125]. Rhodanobacter has been reported to remove uranium from uranium-contaminated sites [126,127,128], and Variovorax has been described as a genus capable of degrading metaldehyde [129].

5. Conclusions

This work presents, for the first time, the microbial diversity associated with ferruginous caves from the Brazilian Iron Quadrangle identified by high-throughput sequencing of the 16S rRNA gene, bringing new findings that could potentially be used in biotechnology. In addition, great variations in richness and abundance were observed among ferruginous caves, including those the Chapada de Canga region that are geographically extremely close. Comparisons of prokaryotic diversity in floor and ceiling samples also showed differences in microbiomes within individual caves, which may have a direct correlation, respectively, with the chronology and geomorphological processes linked to the genesis of these cavities. Finally, the data presented in this study highlight the need for conservation of these areas exposed to intense anthropogenic pressure.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/d13100494/s1, Figure S1: Maps characterizing the iron quadrilateral region, Figure S2: Diversity indices, Figure S3: Functional metabolic prediction, Figure S4: Plant root system throughout the structure of the cavities, Table S1: Sequencing data of collected samples, Table S2: Beta-Diversity indices analysis, Table S3: The most abundant OTUs by sample and the unique OTUs by cave, Table S4: Relative abundance (%) of all OTUs identified, and Table S5: Features of the genera found on the floor, ceiling, and biofilm hanging from the ceiling of cave 2 in Chapada de Canga.

Author Contributions

C.G.C.L., M.M.V., L.M.M., L.H.Y.K. and F.F.d.C. collected the cave samples. C.G.C.L., L.O.P., É.B.F. and L.M.M. conceived and designed all experiments. C.G.C.L., L.M.M., L.Â.G.M., S.E.S.G. and J.C.S. analyzed the biological data. L.M.M., L.Â.G.M. and J.C.S. contributed with reagents, materials, and analysis tools. C.G.C.L., L.M.M., P.N.C.V. and A.A.R.S. prepared the figures and tables. C.G.C.L., M.M.V., É.B.F., L.O.P., R.A.B.A., I.F.C., N.P.F., L.C.C.G., W.L.C., L.Â.G.M., F.F.d.C., L.H.Y.K., P.N.C.V., A.A.R.S., J.C.S. and L.M.M. wrote, revised, and approved the final version of the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—Brazil (CAPES)—Finance Code 001—(the BIGA Project, CFP 51/2013, process 3385/2013), Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq Process 481226/2013-3), and Fundação de Amparo à Pesquisa do Estado de Minas Gerais—FAPEMIG (process APQ-02387-14 and process APQ-02357-17). L.M.M., J.C.S., A.A.R.S. have a research fellowship from CNPq. L.M.M. has UFOP grants. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are publicly available. This data can be found in the Supplementary Materials and here https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/sra/?term=PRJNA575049. Registration date: 30 September 2019.

Acknowledgments

We thank the members of the Laboratory of Genomics and Bacteria-Environment Interaction of the Federal University of Ouro Preto (UFOP) and members of the Bioinformatics Laboratory of the University of São Paulo (USP) for scientific support, especially Carlos Morais Piroupo by the help with the 16S rRNA gene amplicon sequencing data submission.

Conflicts of Interest

The authors declare no conflict of interest.

Ethics Approval and Consent to Participate

The cave samples acquisition and analysis were authorized by SISBIO process 54015-1 (authentication code: 13676514).

References

  1. Auler, A.; Farrant, A. A brief introduction to karst and caves in Brazil. Proc. Univ. Bristol Spelaeol. Soc. 1996, 20, 187–200. [Google Scholar]
  2. Piccini, L. Karst in siliceous rocks—Karst landforms and caves in the Auyan-Tepui Massif (Est. Bolivar, Venezuela). Int. J. Speleol. 1995, 24, 41–54. [Google Scholar] [CrossRef]
  3. Silveira, O.; Maurity, C.; Pinheiro, R.; Henriques, A.; Kern, D.; Souza, S. Estudo das cavernas da província espeleológica arenítica de Monte Alegre—PA. Cad. Geociênc. 1995, 15, 57–63. [Google Scholar]
  4. Hardt, R.; Pinto, S.d.A.F. Carste Em Litologias Não Carbonáticas. Rev. Bras. De Geomorfol. 2009, 10, 99–105. [Google Scholar] [CrossRef] [Green Version]
  5. Piccini, L.; Mecchia, M. Solution weathering rate and origin of karst landforms and caves in the quartzite of Auyan-tepui (Gran Sabana, Venezuela). Geomorphology 2009, 106, 15–25. [Google Scholar] [CrossRef]
  6. Gomes, M.; Azevedo, Ú.; Ferreira, R.; Goulart, F. Landscape Fragmentation around Ferruginous Caves of the Iron Quadrangle, Minas Gerais, Brazil. Cuad. Geogr. Rev. Colomb. Geogr. 2019, 28. [Google Scholar] [CrossRef]
  7. Simmons, G.C. Canga Caves in the Quadrilátero Ferrífero, Minas Gerais, Brazil. Bull. Natl. Speleol. Soc. 1963, 22, 66–72. [Google Scholar]
  8. Pereira, M.; Stávale, Y.; Salgado, A. Estudo da gênese das cavidades e depressões em minério de ferro—Quadrilátero Ferrífero/MG: Serras do Rola Moça e do Gandarela. Rev. Bras. Geomorfol. 2013, 13. [Google Scholar] [CrossRef] [Green Version]
  9. Pereira, M.; Rodet, J.; Salgado, A. Aspectos genéticos e morfológicos das cavidades naturais da Serra da Piedade, Quadrilatero Ferrifero/MG (Brasil). Rev. Bras. Geomorfol. 2012, 13, 465–476. [Google Scholar] [CrossRef] [Green Version]
  10. Calux, A.S.; Cassimiro, J.R.; Salgado, A.A.R. Caves in iron formations in the Quadrilátero Ferrífero, Minas Gerais, southeastern Brazil: Lithological, morphological and hydrological settings and speleogenesis. Z. Geomorphol. 2019, 62, 125–144. [Google Scholar] [CrossRef]
  11. Dorr, J.V.N. Physiographic, Stratigraphic and Structural Development of the Quadrilátero Ferrífero Minas Gerais, Brazil, 2nd ed.; United States Government Printing Office: Washington, DC, USA, 1969; 110p. [Google Scholar]
  12. Myers, N.; Mittermeier, R.A.; Mittermeier, C.G.; da Fonseca, G.A.; Kent, J. Biodiversity hotspots for conservation priorities. Nature 2000, 403, 853–858. [Google Scholar] [CrossRef]
  13. Monteiro, H.S.; Vasconcelos, P.M.; Farley, K.A.; Spier, C.A.; Mello, C.L. (U–Th)/He geochronology of goethite and the origin and evolution of cangas. Geochim. Cosmochim. Acta 2014, 131, 267–289. [Google Scholar] [CrossRef]
  14. Salgado, A.A.R.; Carmo, F.F. Quadrilátero Ferrífero: A Beautiful and Neglected Landscape Between the Gold and Iron Ore Reservoirs. In Landscapes and Landforms of Brazil; Vieira, B.C., Salgado, A.A.R., Santas, L., Eds.; Springer: London, UK, 2015; pp. 319–330. [Google Scholar]
  15. Jacobi, C.M.; Carmo, F.F.; Campos, I.C. Soaring extinction threats to endemic plants in Brazilian metal-rich regions. Ambio 2011, 40, 540–543. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Carmo, F.F.; Jacobi, C.M. Diversity and plant trait-soil relationships among rock outcrops in the Brazilian Atlantic rainforest. Plant. Soil 2016, 403, 7–20. [Google Scholar] [CrossRef]
  17. Ferreira, R.L.; Oliveira, M.P.A.; Silva, M.S. Subterranean Biodiversity in Ferruginous Landscapes: Analysis and Synthesis. In Cave Ecology; Moldovan, O., Kováč, Ľ., Halse, S., Eds.; Springer: Cham, Switzerland, 2018; pp. 435–447. [Google Scholar]
  18. Calux, A.S.; Cassimiro, J.R. Geoespeleologia das cavernas em rochas ferríferas: Aspectos dimensionais, morfológicos, hidrológicos e sedimentares. In Patrimônio Espeleológico em Rochas Ferruginosas: Propostas para Sua Conservação no Quadrilátero Ferrífero, Minas Gerais; Sociedade Brasileira de Espeleologia, 2015; 26p. [Google Scholar]
  19. Dutra, G. Síntese dos processos de gênese de cavidades em litologias de ferro. In Proceedings of the 32° Congresso Brasileiro de Espeleologia, Barreiras, Brazil, 11–14 July 2013; pp. 415–426. [Google Scholar]
  20. Auler, A.; Piló, L.; Parker, C.; Senko, J.; Sasowsky, I.D.; Barton, H. Hypogene Cave Patterns in Iron Ore Caves: Convergence of Forms or Processes? In Hypogene Cave Morphologies; Klimchouk, A., Sasowsky, I.D., Mylroie, J., Engel, S.A., Engel, A.S., Eds.; Karst Waters Institute Special Publication: San Salvador Island, Bahamas, 2014; Volume 18, pp. 15–19. [Google Scholar]
  21. Culver, D.C. Cave Life: Evolution and Ecology; Harvard University Press: Cambridge, MA, USA, 1982; p. 189. [Google Scholar]
  22. Scherer, R.S.; Piló, L.B.; Souza-Filho, P.W.M.; Barata, A.S.; Scherer, B.S. Ocorrência de espeleotemas e feições morfológicas raras em cavernas ferríferas da Serra dos Carajás, no Pará. In Proceedings of the 34° Congresso Brasileiro de Espeleologia, Ouro Preto, Brazil, 13–18 June 2017; pp. 409–416. [Google Scholar]
  23. Parker, C.W.; Auler, A.S.; Barton, M.D.; Sasowsky, I.D.; Senko, J.M.; Barton, H.A. Fe(III) Reducing Microorganisms from Iron Ore Caves Demonstrate Fermentative Fe(III) Reduction and Promote Cave Formation. Geomicrobiol. J. 2018, 35, 311–322. [Google Scholar] [CrossRef]
  24. Gonçalves, D.F.; Sousa, D.L. Aspectos morfológicos de espeleotemas em feições pseudocársticas da região de Carajás – PA. In Proceedings of the 31° Congresso Brasileiro de Espeleologia, Ponta Grossa, Brazil, 21–24 2011; pp. 141–145. [Google Scholar]
  25. Carmo, F.F.; Kamino, L.H.Y. Geossistemas Ferruginosos do Brasil—Áreas Prioritárias para Conservação da Diversidade Geológica e Biológica, Patrimônio Cultural e Serviços Ambientais, 1st ed.; 3i Editora: Belo Horizonte, Brazil, 2015; Volume 1. [Google Scholar]
  26. Salles, D.M.; Carmo, F.F.; Jacobi, C.M. Habitat Loss Challenges the Conservation of Endemic Plants in Mining-Targeted Brazilian Mountains. Environ. Conserv. 2019, 46, 140–146. [Google Scholar] [CrossRef] [Green Version]
  27. Auler, A.S.; Piló, L.B. Caves and Mining in Brazil: The Dilemma of Cave Preservation Within a Mining Context. In Hydrogeological and Environmental Investigations in Karst Systems; Andreo, B., Carrasco, F., Durán, J.J., Jiménez, P., LaMoreaux, J.W., Eds.; Springer: Berlin/Heidelberg, Germany, 2015; pp. 487–496. [Google Scholar]
  28. Kamino, L.H.Y.; Pereira, E.O.; Carmo, F.F. Conservation paradox: Large-scale mining waste in protected areas in two global hotspots, southeastern Brazil. Ambio 2020, 49, 1629–1638. [Google Scholar] [CrossRef]
  29. Hatje, V.; Pedreira, R.M.A.; Rezende, C.E.; Schettini, C.A.F.; Souza, G.C.; Marin, D.C.; Hackspacher, P.C. The environmental impacts of one of the largest tailing dam failures worldwide. Sci. Rep. 2017, 7, 10706. [Google Scholar] [CrossRef] [PubMed]
  30. Carmo, F.F.; Kamino, L.H.Y.; Junior, R.T.; Campos, I.C.; Carmo, F.F.; Silvino, G.; Castro, K.J.d.S.X.d.; Mauro, M.L.; Rodrigues, N.U.A.; Miranda, M.P.d.S.; et al. Fundão tailings dam failures: The environment tragedy of the largest technological disaster of Brazilian mining in global context. Perspect. Ecol. Conserv. 2017, 15, 145–151. [Google Scholar] [CrossRef]
  31. Parker, C.W.; Wolf, J.A.; Auler, A.S.; Barton, H.A.; Senko, J.M. Microbial Reducibility of Fe(III) Phases Associated with the Genesis of Iron Ore Caves in the Iron Quadrangle, Minas Gerais, Brazil. Minerals 2013, 3, 395. [Google Scholar] [CrossRef] [Green Version]
  32. Makhalanyane, T.P.; Valverde, A.; Gunnigle, E.; Frossard, A.; Ramond, J.B.; Cowan, D.A. Microbial ecology of hot desert edaphic systems. FEMS Microbiol. Rev. 2015, 39, 203–221. [Google Scholar] [CrossRef] [PubMed]
  33. Jacobi, C.M.; Carmo, F.F.; Vincent, R.C.; Stehmann, J.R. Plant communities on ironstone outcrops: A diverse and endangered Brazilian ecosystem. Biodivers. Conserv. 2007, 16, 2185–2200. [Google Scholar] [CrossRef]
  34. Jacobi, C.M.; Carmo, F.F.; Carmo, F.F.; Campos, I.C. Iron geosystems: Priority areas for conservation in Brazil. In Mining in Ecologically Sensitive Landscapes; Tibbett, M., Ed.; CRC Press: Boca Raton, FL, USA, 2015. [Google Scholar]
  35. Salter, S.J.; Cox, M.J.; Turek, E.M.; Calus, S.T.; Cookson, W.O.; Moffatt, M.F.; Turner, P.; Parkhill, J.; Loman, N.J.; Walker, A.W. Reagent and laboratory contamination can critically impact sequence-based microbiome analyses. BMC Biol. 2014, 12, 87. [Google Scholar] [CrossRef] [Green Version]
  36. Fouhy, F.; Clooney, A.G.; Stanton, C.; Claesson, M.J.; Cotter, P.D. 16S rRNA gene sequencing of mock microbial populations- impact of DNA extraction method, primer choice and sequencing platform. BMC Microbiol. 2016, 16, 123. [Google Scholar] [CrossRef]
  37. Ion Personal Genome Machine™ Performance Overview. Application Note. Ion. Torrent Life Technol. 2011. Available online: http://tools.thermofisher.com/content/sfs/brochures/cms_094272.pdf (accessed on 15 September 2021).
  38. Song, L.; Huang, W.; Kang, J.; Huang, Y.; Ren, H.; Ding, K. Comparison of error correction algorithms for Ion Torrent PGM data: Application to hepatitis B virus. Sci. Rep. 2017, 7, 8106. [Google Scholar] [CrossRef]
  39. Feng, W.; Zhao, S.; Xue, D.; Song, F.; Li, Z.; Chen, D.; He, B.; Hao, Y.; Wang, Y.; Liu, Y. Improving alignment accuracy on homopolymer regions for semiconductor-based sequencing technologies. BMC Genom. 2016, 17, 521. [Google Scholar] [CrossRef] [Green Version]
  40. Golan, D.; Medvedev, P. Using state machines to model the Ion Torrent sequencing process and to improve read error rates. Bioinformatics 2013, 29, i344–i351. [Google Scholar] [CrossRef] [Green Version]
  41. Pylro, V.S.; Roesch, L.F.W.; Ortega, J.M.; Amaral, A.M.; Tótola, M.R.; Hirsch, P.R.; Rosado, A.S.; Góes-Neto, A.; Costa da Silva, A.L.; Rosa, C.A.; et al. Brazilian Microbiome Project: Revealing the Unexplored Microbial Diversity—Challenges and Prospects. Microb. Ecol. 2014, 67, 237–241. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Schmieder, R.; Edwards, R. Quality control and preprocessing of metagenomic datasets. Bioinformatics 2011, 27, 863–864. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Edgar, R.C. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 2010, 26, 2460–2461. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Edgar, R.C. UPARSE: Highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 2013, 10, 996. [Google Scholar] [CrossRef] [PubMed]
  45. Wang, Q.; Garrity, G.M.; Tiedje, J.M.; Cole, J.R. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 2007, 73, 5261–5267. [Google Scholar] [CrossRef] [Green Version]
  46. Bolyen, E.; Rideout, J.R.; Dillon, M.R.; Bokulich, N.A.; Abnet, C.C.; Al-Ghalith, G.A.; Alexander, H.; Alm, E.J.; Arumugam, M.; Asnicar, F.; et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 2019, 37, 852–857. [Google Scholar] [CrossRef]
  47. Caporaso, J.G.; Bittinger, K.; Bushman, F.D.; DeSantis, T.Z.; Andersen, G.L.; Knight, R. PyNAST: A flexible tool for aligning sequences to a template alignment. Bioinformatics 2010, 26, 266–267. [Google Scholar] [CrossRef] [Green Version]
  48. Price, M.N.; Dehal, P.S.; Arkin, A.P. FastTree: Computing large minimum evolution trees with profiles instead of a distance matrix. Mol. Biol Evol. 2009, 26, 1641–1650. [Google Scholar] [CrossRef]
  49. Kruskal, W.H.; Wallis, W.A. Use of Ranks in One-Criterion Variance Analysis. J. Am. Stat. Assoc. 1952, 47, 583–621. [Google Scholar] [CrossRef]
  50. McMurdie, P.J.; Holmes, S. phyloseq: An R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE 2013, 8, e61217. [Google Scholar] [CrossRef] [Green Version]
  51. Shannon, C.E.; Weaver, W. The Mathematical Theory of Communication; University of Illinois Press: Champaign, IL, USA, 1949; 117p. [Google Scholar]
  52. Simpson, E.H. Measurement of Diversity. Nature 1949, 163, 688. [Google Scholar] [CrossRef]
  53. Chao, A. Non-parametric estimation of the classes in a population. Scand. J. Stat. 1984, 11, 265–270. [Google Scholar] [CrossRef]
  54. Hamady, M.; Lozupone, C.; Knight, R. Fast UniFrac: Facilitating high-throughput phylogenetic analyses of microbial communities including analysis of pyrosequencing and PhyloChip data. ISME J. 2010, 4, 17–27. [Google Scholar] [CrossRef]
  55. Lozupone, C.; Knight, R. UniFrac: A new phylogenetic method for comparing microbial communities. Appl. Environ. Microbiol. 2005, 71, 8228–8235. [Google Scholar] [CrossRef] [Green Version]
  56. Wickham, H. ggplot2: Elegant Graphics for Data Analysis, 2nd ed.; Springer International Publishing: New York, NY, USA, 2016; p. 260. [Google Scholar]
  57. Aitchison, J.; Barceló-Vidal, C.; Martín-Fernández, J.A.; Pawlowsky-Glahn, V. Logratio Analysis and Compositional Distance. Math. Geol. 2000, 32, 271–275. [Google Scholar] [CrossRef]
  58. Douglas, G.M.; Maffei, V.J.; Zaneveld, J.R.; Yurgel, S.N.; Brown, J.R.; Taylor, C.M.; Huttenhower, C.; Langille, M.G.I. PICRUSt2 for prediction of metagenome functions. Nat. Biotechnol. 2020. [Google Scholar] [CrossRef]
  59. Kanehisa, M.; Sato, Y.; Kawashima, M.; Furumichi, M.; Tanabe, M. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 2015, 44, D457–D462. [Google Scholar] [CrossRef] [Green Version]
  60. Parks, D.H.; Tyson, G.W.; Hugenholtz, P.; Beiko, R.G. STAMP: Statistical analysis of taxonomic and functional profiles. Bioinformatics 2014, 30, 3123–3124. [Google Scholar] [CrossRef] [Green Version]
  61. Maurity, C.; Kotschoubey, B. Evolução recente da cobertura de alteração no Platô N1—Serra dos Carajás-PA. Degradação, pseudocarstificação, espeleotemas. Bol. Mus. Para. Emilio Goeldi Sér. Ciênc. Terra 1995, 7, 331–362. [Google Scholar]
  62. Caneschi, W.L.; Felestrino, É.B.; Fonseca, N.P.; Villa, M.M.; Lemes, C.G.C.; Cordeiro, I.F.; Assis, R.d.A.B.; Sanchez, A.B.; Vieira, I.T.; Kamino, L.H.Y.; et al. Brazilian Ironstone Plant Communities as Reservoirs of Culturable Bacteria With Diverse Biotechnological Potential. Front. Microbiol. 2018, 9, 1638. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Felestrino, É.B.; Vieira, I.T.; Caneschi, W.L.; Cordeiro, I.F.; Assis, R.A.B.; Lemes, C.G.C.; Fonseca, N.P.; Sanchez, A.B.; Cepeda, J.C.C.; Ferro, J.A.; et al. Biotechnological potential of plant growth-promoting bacteria from the roots and rhizospheres of endemic plants in ironstone vegetation in southeastern Brazil. World J. Microbiol. Biotechnol. 2018, 34, 156. [Google Scholar] [CrossRef] [PubMed]
  64. Drummond, G.M.; Martins, C.S.; Machado, A.B.M.; Sebaio, F.A.; Antonini, Y. Biodiversidade em Minas Gerais: Um Atlas para Sua Conservação, 2nd ed.; Fundação Biodiversitas: Belo Horizonte, Brazil, 2005; 222p. [Google Scholar]
  65. Giuliett, A.M.; Rapini, A.; Gomes de Andrade, M.; Queiroz, L.; Da Silva, J.M. Plantas Raras do Brasil; Conservação Internacional: Belo Horizonte, Brazil, 2009; 495p. [Google Scholar]
  66. Carmo, F.F.; Jacobi, C.M. A vegetação de canga no Quadrilátero Ferrífero, Minas Gerais: Caracterização e contexto fitogeográfico. Rodriguésia 2013, 64, 527–541. [Google Scholar] [CrossRef] [Green Version]
  67. Gastauer, M.; Vera, M.P.O.; de Souza, K.P.; Pires, E.S.; Alves, R.; Caldeira, C.F.; Ramos, S.J.; Oliveira, G. A metagenomic survey of soil microbial communities along a rehabilitation chronosequence after iron ore mining. Sci. Data 2019, 6, 190008. [Google Scholar] [CrossRef] [PubMed]
  68. Levett, A.; Gagen, E.; Shuster, J.; Rintoul, L.; Tobin, M.; Vongsvivut, J.; Bambery, K.; Vasconcelos, P.; Southam, G. Evidence of biogeochemical processes in iron duricrust formation. J. S. Am. Earth Sci. 2016, 71, 131–142. [Google Scholar] [CrossRef] [Green Version]
  69. Parker, C.; Auler, A.; Senko, J.; Sasowsky, I.; Piló, L.; Smith, M.; Johnston, M.; Barton, H. Microbial Iron Cycling and Biospeleogenesis: Cave Development in the Carajás Formation, Brazil. In Proceedings of the Proceedings of the 16th International Congress of Speleology, Brno, Czech Republic, 21 July 2013. [Google Scholar]
  70. Figueira, R.L.; Coimbra Horbe, A.M.; Herrera Aragón, F.F.; Gonçalves, D.F. Exotic sulphate and phosphate speleothems in caves from eastern Amazonia (Carajás, Brazil): Crystallographic and chemical insights. J. S. Am. Earth Sci. 2019, 90, 412–422. [Google Scholar] [CrossRef]
  71. Jacoby, R.; Peukert, M.; Succurro, A.; Koprivova, A.; Kopriva, S. The Role of Soil Microorganisms in Plant Mineral Nutrition-Current Knowledge and Future Directions. Front. Plant Sci. 2017, 8, 1617. [Google Scholar] [CrossRef] [Green Version]
  72. Adetutu, E.M.; Thorpe, K.; Shahsavari, E.; Bourne, S.; Cao, X.S.; Fard, R.; Kirby, G.; Ball, A. Bacterial community survey of sediments at Naracoorte Caves, Australia. Int. J. Speleol. 2012, 41, 137–147. [Google Scholar] [CrossRef] [Green Version]
  73. Barton, H.A.; Taylor, N.M.; Lubbers, B.R.; Pemberton, A.C. DNA extraction from low-biomass carbonate rock: An improved method with reduced contamination and the low-biomass contaminant database. J. Microbiol. Methods 2006, 66, 21–31. [Google Scholar] [CrossRef]
  74. Ortiz, M.; Neilson, J.W.; Nelson, W.M.; Legatzki, A.; Byrne, A.; Yu, Y.; Wing, R.A.; Soderlund, C.A.; Pryor, B.M.; Pierson, L.S., III; et al. Profiling bacterial diversity and taxonomic composition on speleothem surfaces in Kartchner Caverns, AZ. Microb. Ecol. 2013, 65, 371–383. [Google Scholar] [CrossRef]
  75. Gulecal-Pektas, Y. Bacterial Diversity and Composition in Oylat Cave (Turkey) with Combined Sanger/Pyrosequencing Approach. Pol. J. Microbiol. 2016, 65, 69–75. [Google Scholar] [CrossRef] [Green Version]
  76. D’Auria, G.; Artacho, A.; Rojas, R.A.; Bautista, J.S.; Mendez, R.; Gamboa, M.T.; Gamboa, J.R.; Gomez-Cruz, R. Metagenomics of Bacterial Diversity in Villa Luz Caves with Sulfur Water Springs. Genes 2018, 9, 55. [Google Scholar] [CrossRef] [Green Version]
  77. Ling, L.L.; Schneider, T.; Peoples, A.J.; Spoering, A.L.; Engels, I.; Conlon, B.P.; Mueller, A.; Schaberle, T.F.; Hughes, D.E.; Epstein, S.; et al. A new antibiotic kills pathogens without detectable resistance. Nature 2015, 517, 455–459. [Google Scholar] [CrossRef]
  78. Nichols, D.; Cahoon, N.; Trakhtenberg, E.M.; Pham, L.; Mehta, A.; Belanger, A.; Kanigan, T.; Lewis, K.; Epstein, S.S. Use of ichip for high-throughput in situ cultivation of “uncultivable” microbial species. Appl. Environ. Microbiol. 2010, 76, 2445–2450. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Nemoto, P.A.; Zhang, X.; Djuric, S.; Ma, Z. Novel 12-membered ring macrolides with activity against erythromycin-resistant organisms. J. Antibiot. 2003, 56, 392–398. [Google Scholar] [CrossRef] [Green Version]
  80. Vickers, A.A.; Chopra, I.; O’Neill, A.J. Intrinsic novobiocin resistance in Staphylococcus saprophyticus. Antimicrob. Agents Chemother. 2007, 51, 4484–4485. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Wang, M.; Zhang, J.; He, S.; Yan, X. A Review Study on Macrolides Isolated from Cyanobacteria. Mar. Drugs 2017, 15, 126. [Google Scholar] [CrossRef] [Green Version]
  82. Singh, G.; Jayanarayan, K.G.; Dey, C.S. Novobiocin induces apoptosis-like cell death in topoisomerase II over-expressing arsenite resistant Leishmania donovani. Mol. Biochem. Parasitol. 2005, 141, 57–69. [Google Scholar] [CrossRef]
  83. Romero, L.; Huerfano, C.; Grillo-Ardila, C.F. Macrolides for treatment of Haemophilus ducreyi infection in sexually active adults. Cochrane Database Syst. Rev. 2017, 12, CD012492. [Google Scholar] [CrossRef] [PubMed]
  84. Kelly, C.; Chalmers, J.D.; Crossingham, I.; Relph, N.; Felix, L.M.; Evans, D.J.; Milan, S.J.; Spencer, S. Macrolide antibiotics for bronchiectasis. Cochrane Database Syst. Rev. 2018, 3, CD012406. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. O’Hagan, D. The polyketide metabolites; E. Horwood: Sawston, UK, 1991; 176p. [Google Scholar]
  86. Lewin, G.R.; Carlos, C.; Chevrette, M.G.; Horn, H.A.; McDonald, B.R.; Stankey, R.J.; Fox, B.G.; Currie, C.R. Evolution and Ecology of Actinobacteria and Their Bioenergy Applications. Annu Rev. Microbiol. 2016, 70, 235–254. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Palaniyandi, S.A.; Yang, S.H.; Zhang, L.; Suh, J.W. Effects of actinobacteria on plant disease suppression and growth promotion. Appl. Microbiol. Biotechnol. 2013, 97, 9621–9636. [Google Scholar] [CrossRef]
  88. Harshey, R.M.; Kawagishi, I.; Maddock, J.; Kenney, L.J. Function, Diversity, and Evolution of Signal Transduction in Prokaryotes. Dev. Cell 2003, 4, 459–465. [Google Scholar] [CrossRef] [Green Version]
  89. Nazir, R.; Mazurier, S.; Yang, P.; Lemanceau, P.; van Elsas, J.D. The Ecological Role of Type Three Secretion Systems in the Interaction of Bacteria with Fungi in Soil and Related Habitats Is Diverse and Context-Dependent. Front. Microbiol. 2017, 8, 38. [Google Scholar] [CrossRef] [Green Version]
  90. Desai, M.S.; Assig, K.; Dattagupta, S. Nitrogen fixation in distinct microbial niches within a chemoautotrophy-driven cave ecosystem. ISME J. 2013, 7, 2411. [Google Scholar] [CrossRef] [Green Version]
  91. Carmo, F.F.; Kamino, L.H.Y. Chapada de Canga: Uma introdução In Chapada de Canga: Patrimônio Natural e Cultural de Relevante Interesse para Conservação; Carmo, F.F., Kamino, L.H.Y., Eds.; 3i Editora: Belo Horizonte, Brazil, 2017; pp. 11–24. [Google Scholar]
  92. Lopes, L.M.N. O rompimento da barragem de Mariana e seus impactos socioambientais. Sinapse Múlt. 2016, 5, 1–14. [Google Scholar]
  93. Freitas, C.M.; Silva, M.A.; Menezes, F.C. O desastre na barragem de mineração da samarco—Fratura exposta dos limites do Brasil na redução de risco de desasastres. Ciência Cult. 2016, 68, 25–30. [Google Scholar] [CrossRef] [Green Version]
  94. Barka, E.A.; Vatsa, P.; Sanchez, L.; Gaveau-Vaillant, N.; Jacquard, C.; Meier-Kolthoff, J.P.; Klenk, H.P.; Clement, C.; Ouhdouch, Y.; van Wezel, G.P. Taxonomy, Physiology, and Natural Products of Actinobacteria. Microbiol. Mol. Biol. Rev. MMBR 2016, 80, 1–43. [Google Scholar] [CrossRef] [Green Version]
  95. Qin, S.; Chen, H.H.; Zhao, G.Z.; Li, J.; Zhu, W.Y.; Xu, L.H.; Jiang, J.H.; Li, W.J. Abundant and diverse endophytic actinobacteria associated with medicinal plant Maytenus austroyunnanensis in Xishuangbanna tropical rainforest revealed by culture-dependent and culture-independent methods. Environ. Microbiol. Rep. 2012, 4, 522–531. [Google Scholar] [CrossRef] [PubMed]
  96. Qin, S.; Li, W.J.; Dastager, S.G.; Hozzein, W.N. Actinobacteria in Special and Extreme Habitats: Diversity, Function Roles, and Environmental Adaptations. Front. Microbiol. 2016, 7, 1415. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  97. Hamedi, J.; Mohammadipanah, F.; Ventosa, A. Systematic and biotechnological aspects of halophilic and halotolerant actinomycetes. Extrem. Life Under Extrem. Cond. 2013, 17, 1–13. [Google Scholar] [CrossRef] [PubMed]
  98. Genilloud, O. Actinomycetes: Still a source of novel antibiotics. Nat. Prod. Rep. 2017, 34, 1203–1232. [Google Scholar] [CrossRef] [PubMed]
  99. Piló, L.B.; Coelho, A.; Reino, J.C.R. Geoespeleologia em rochas ferríferas: Cenário atual e conservação. In Geossistemas Ferruginosos do Brasil: Áreas Prioritárias para Conservação da Diversidade Geológica e Biológica, Patrimônio Cultural e Serviços Ambientais; Carmo, F.F.d., Kamino, L.H.Y., Eds.; 3i Editora: Belo Horizonte, Brazil, 2015; pp. 125–148. [Google Scholar]
  100. Sallstedt, T.; Ivarsson, M.; Lundberg, J.; Sjöberg, R.; Romaní, J.R.V. Speleothem and biofilm formation in a granite/dolerite cave, Northern Sweden. Int. J. Speleol. 2014, 43, 305–313. [Google Scholar] [CrossRef] [Green Version]
  101. Turrini, P.; Tescari, M.; Visaggio, D.; Pirolo, M.; Lugli, G.A.; Ventura, M.; Frangipani, E.; Visca, P. The microbial community of a biofilm lining the wall of a pristine cave in Western New Guinea. Microbiol. Res. 2020, 241, 126584. [Google Scholar] [CrossRef] [PubMed]
  102. Caccavo, F., Jr.; Lonergan, D.J.; Lovley, D.R.; Davis, M.; Stolz, J.F.; McInerney, M.J. Geobacter sulfurreducens sp. nov., a hydrogen- and acetate-oxidizing dissimilatory metal-reducing microorganism. Appl. Environ. Microbiol. 1994, 60, 3752–3759. [Google Scholar] [CrossRef] [Green Version]
  103. Emerson, D.; Moyer, C. Isolation and characterization of novel iron-oxidizing bacteria that grow at circumneutral pH. Appl. Environ. Microbiol. 1997, 63, 4784–4792. [Google Scholar] [CrossRef] [Green Version]
  104. Johnson, D.B.; Hallberg, K.B.; Hedrich, S. Uncovering a microbial enigma: Isolation and characterization of the streamer-generating, iron-oxidizing, acidophilic bacterium “Ferrovum myxofaciens”. Appl. Environ. Microbiol. 2014, 80, 672–680. [Google Scholar] [CrossRef] [Green Version]
  105. Kimura, S.; Bryan, C.; B Hallberg, K.; Johnson, D. Biodiversity and Geochemistry of an extremely acidic, low-temperature subterranean environment sustained by chemolithotrophy. Environ. Microbiol. 2011, 13, 2092–2104. [Google Scholar] [CrossRef]
  106. Northup, D.E.; Barns, S.M.; Yu, L.E.; Spilde, M.N.; Schelble, R.T.; Dano, K.E.; Crossey, L.J.; Connolly, C.A.; Boston, P.J.; Natvig, D.O.; et al. Diverse microbial communities inhabiting ferromanganese deposits in Lechuguilla and Spider Caves. Environ. Microbiol. 2003, 5, 1071–1086. [Google Scholar] [CrossRef]
  107. Pacton, M.; Breitenbach, S.; Lechleitner, F.; Vaks, A.; Rollion-Bard, C.; Gutareva, O.S.; Osintcev, A.V.; Vasconcelos, C. The role of microorganisms in the formation of a stalactite in Botovskaya Cave, Siberia—Paleoenvironmental implications. Biogeosciences 2013, 10, 6115–6130. [Google Scholar] [CrossRef] [Green Version]
  108. Malvankar, N.S.; Vargas, M.; Nevin, K.P.; Franks, A.E.; Leang, C.; Kim, B.-C.; Inoue, K.; Mester, T.; Covalla, S.F.; Johnson, J.P.; et al. Tunable metallic-like conductivity in microbial nanowire networks. Nat. Nanotechnol. 2011, 6, 573. [Google Scholar] [CrossRef]
  109. Prakash, D.; Gabani, P.; Chandel, A.K.; Ronen, Z.; Singh, O.V. Bioremediation: A genuine technology to remediate radionuclides from the environment. Microb. Biotechnol. 2013, 6, 349–360. [Google Scholar] [CrossRef] [PubMed]
  110. Lee, K.Y.; Bosch, J.; Meckenstock, R.U. Use of metal-reducing bacteria for bioremediation of soil contaminated with mixed organic and inorganic pollutants. Environ. Geochem. Health 2012, 34, 135–142. [Google Scholar] [CrossRef]
  111. Newsome, L.; Morris, K.; Lloyd, J.R. The biogeochemistry and bioremediation of uranium and other priority radionuclides. Chem. Geol. 2014, 363, 164–184. [Google Scholar] [CrossRef]
  112. Dedysh, S.N.; Khmelenina, V.N.; Suzina, N.E.; Trotsenko, Y.A.; Semrau, J.D.; Liesack, W.; Tiedje, J.M. Methylocapsa acidiphila gen. nov., sp. nov., a novel methane-oxidizing and dinitrogen-fixing acidophilic bacterium from Sphagnum bog. Int. J. Syst. Evol. Microbiol. 2002, 52, 251–261. [Google Scholar] [CrossRef] [Green Version]
  113. Lind, H.; Jonsson, H.; Schnürer, J. Antifungal effect of dairy propionibacteria—Contribution of organic acids. Int. J. Food Microbiol. 2005, 98, 157–165. [Google Scholar] [CrossRef]
  114. Le Lay, C.; Coton, E.; Le Blay, G.; Chobert, J.M.; Haertle, T.; Choiset, Y.; van Long, N.N.; Meslet-Cladiere, L.; Mounier, J. Identification and quantification of antifungal compounds produced by lactic acid bacteria and propionibacteria. Int. J. Food Microbiol. 2016, 239, 79–85. [Google Scholar] [CrossRef]
  115. Williams, S.T.; Khan, M.R. Antibiotics—A soil microbiologist’s viewpoint. Postepy Hig. Med. Dosw. 1974, 28, 395–408. [Google Scholar]
  116. Sen, R.; Ishak, H.D.; Estrada, D.; Dowd, S.E.; Hong, E.; Mueller, U.G. Generalized antifungal activity and 454-screening of Pseudonocardia and Amycolatopsis bacteria in nests of fungus-growing ants. Proc. Natl. Acad. Sci. USA 2009, 106, 17805–17810. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  117. Atashgahi, S.; Liebensteiner, M.G.; Janssen, D.B.; Smidt, H.; Stams, A.J.M.; Sipkema, D. Microbial Synthesis and Transformation of Inorganic and Organic Chlorine Compounds. Front. Microbiol. 2018, 9, 3079. [Google Scholar] [CrossRef] [PubMed]
  118. Schimel, J.; Schaeffer, S. Microbial control over carbon cycling in soil. Front. Microbiol. 2012, 3, 348. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  119. Takai, K. The Nitrogen Cycle: A Large, Fast, and Mystifying Cycle. Microbes Environ. 2019, 34, 223–225. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  120. Yu, X.; Li, Y.; Zhang, C.; Liu, H.; Liu, J.; Zheng, W.; Kang, X.; Leng, X.; Zhao, K.; Gu, Y.; et al. Culturable Heavy Metal-Resistant and Plant Growth Promoting Bacteria in V-Ti Magnetite Mine Tailing Soil from Panzhihua, China. PLoS ONE 2014, 9, e106618. [Google Scholar] [CrossRef] [PubMed]
  121. Ouyang, E.; Lu, Y.; Ouyang, J.; Wang, L.; Wang, X. Bacterial community analysis of anoxic/aeration (A/O) system in a combined process for gibberellin wastewater treatment. PLoS ONE 2017, 12, e0186743. [Google Scholar] [CrossRef] [Green Version]
  122. Daranas, N.; Rosello, G.; Cabrefiga, J.; Donati, I.; Frances, J.; Badosa, E.; Spinelli, F.; Montesinos, E.; Bonaterra, A. Biological control of bacterial plant diseases with Lactobacillus plantarum strains selected for their broad-spectrum activity. Ann. Appl. Biol. 2019, 174, 92–105. [Google Scholar] [CrossRef] [Green Version]
  123. Koppel, N.; Maini Rekdal, V.; Balskus, E.P. Chemical transformation of xenobiotics by the human gut microbiota. Science 2017, 356, eaag2770. [Google Scholar] [CrossRef] [PubMed]
  124. Mahenthiralingam, E.; Urban, T.A.; Goldberg, J.B. The multifarious, multireplicon Burkholderia cepacia complex. Nat. Rev. Microbiol. 2005, 3, 144–156. [Google Scholar] [CrossRef]
  125. Mahenthiralingam, E.; Baldwin, A.; Dowson, C.G. Burkholderia cepacia complex bacteria: Opportunistic pathogens with important natural biology. J. Appl. Microbiol. 2008, 104, 1539–1551. [Google Scholar] [CrossRef] [PubMed]
  126. Prakash, O.; Green, S.J.; Jasrotia, P.; Overholt, W.A.; Canion, A.; Watson, D.B.; Brooks, S.C.; Kostka, J.E. Rhodanobacter denitrificans sp. nov., isolated from nitrate-rich zones of a contaminated aquifer. Int. J. Syst. Evol. Microbiol. 2012, 62, 2457–2462. [Google Scholar] [CrossRef] [Green Version]
  127. Sousa, T.; Chung, A.P.; Pereira, A.; Piedade, A.P.; Morais, P.V. Aerobic uranium immobilization by Rhodanobacter A2-61 through formation of intracellular uranium-phosphate complexes. Met. Integr. Biometal Sci. 2013, 5, 390–397. [Google Scholar] [CrossRef]
  128. Chung, A.; Sousa, T.; Pereira, A.; Morais, P. Microorganisms—Tools for Bioremediation of Uranium Contaminated Environments. Procedia Earth Planet. Sci. 2014, 8, 53–58. [Google Scholar] [CrossRef] [Green Version]
  129. Satola, B.; Wubbeler, J.H.; Steinbuchel, A. Metabolic characteristics of the species Variovorax paradoxus. Appl. Microbiol. Biotechnol. 2013, 97, 541–560. [Google Scholar] [CrossRef]
Figure 1. Thematic maps of the IQ region, MG, Brazil. (a) IQ delimitation and location of the eight caves investigated in this study (Details see Table 1); (b) Lithologies that have natural cavities; (c) Potential for the formation of natural cavities. All maps were obtained from “Instituto Prístino: atlas digital geoambiental. WebGis system (Web Geographical Information System) of free access to the environmental database”. Available online: https://institutopristino.org.br/atlas accessed on 15 June 2020.
Figure 1. Thematic maps of the IQ region, MG, Brazil. (a) IQ delimitation and location of the eight caves investigated in this study (Details see Table 1); (b) Lithologies that have natural cavities; (c) Potential for the formation of natural cavities. All maps were obtained from “Instituto Prístino: atlas digital geoambiental. WebGis system (Web Geographical Information System) of free access to the environmental database”. Available online: https://institutopristino.org.br/atlas accessed on 15 June 2020.
Diversity 13 00494 g001
Figure 2. Global analysis of the 16S rRNA amplicon sequencing profile of the samples. (a) Rarefaction curve using the Observed_Otus metric (p < 0.05). (b) Diversity indices for the characterization of the prokaryotic community associated to cave samples. The expanded version of these results is shown in the Supplementary Figure S2. (c) Analysis of beta-diversity by PcoA using unweighted and weighted unifrac distance (p < 0.05). Orange line highlights the grouping of floor samples and green line the ceiling samples. Quartzitic cave samples are separated from ferruginous cave samples. Note that the floor sample from cave 8 (sample 27) is distant from all other samples (see details at Supplementary Figure S3). * quartzite cave samples.
Figure 2. Global analysis of the 16S rRNA amplicon sequencing profile of the samples. (a) Rarefaction curve using the Observed_Otus metric (p < 0.05). (b) Diversity indices for the characterization of the prokaryotic community associated to cave samples. The expanded version of these results is shown in the Supplementary Figure S2. (c) Analysis of beta-diversity by PcoA using unweighted and weighted unifrac distance (p < 0.05). Orange line highlights the grouping of floor samples and green line the ceiling samples. Quartzitic cave samples are separated from ferruginous cave samples. Note that the floor sample from cave 8 (sample 27) is distant from all other samples (see details at Supplementary Figure S3). * quartzite cave samples.
Diversity 13 00494 g002
Figure 3. Diversity of the phyla present in the ferruginous and quartzitic caves, obtained from the table of OTUs. Acidobacteria, Chloroflexi, Proteobacteria, Actinobacteria, Nitrospirae, and Firmicutes showed relative abundance > 1% in the composition of the total samples.
Figure 3. Diversity of the phyla present in the ferruginous and quartzitic caves, obtained from the table of OTUs. Acidobacteria, Chloroflexi, Proteobacteria, Actinobacteria, Nitrospirae, and Firmicutes showed relative abundance > 1% in the composition of the total samples.
Diversity 13 00494 g003
Figure 4. Relative abundance and diversity analysis. (a) Relative abundances and bacterial diversity at class level in ferruginous and quartzitic caves. The numbers 1 to 8 indicate the eight caves investigated in this study and georeferenced in Figure 1. (b) Heatmap of genus abundances (trend) in the samples collected in caves 1 to 8 (according to Table 1). The numbers in red circles indicate which cave presented greater relative abundance of the genus in comparison to other investigated caves.
Figure 4. Relative abundance and diversity analysis. (a) Relative abundances and bacterial diversity at class level in ferruginous and quartzitic caves. The numbers 1 to 8 indicate the eight caves investigated in this study and georeferenced in Figure 1. (b) Heatmap of genus abundances (trend) in the samples collected in caves 1 to 8 (according to Table 1). The numbers in red circles indicate which cave presented greater relative abundance of the genus in comparison to other investigated caves.
Diversity 13 00494 g004
Figure 5. Relative proportions of the functional metabolic profiles based on sequences for ferruginous and quartzite caves. (a) Relative abundance of sequences between the caves based on Picrust predictions for functional categories inferred at level 3. Categories highlighted in yellow are presented in detail in b. (b) Comparative analysis of proportion of sequences (%) for predicted metabolic functions on genus level between the ferruginous caves and the quartzite cave. Welch’s t-test was used (p < 0.05).
Figure 5. Relative proportions of the functional metabolic profiles based on sequences for ferruginous and quartzite caves. (a) Relative abundance of sequences between the caves based on Picrust predictions for functional categories inferred at level 3. Categories highlighted in yellow are presented in detail in b. (b) Comparative analysis of proportion of sequences (%) for predicted metabolic functions on genus level between the ferruginous caves and the quartzite cave. Welch’s t-test was used (p < 0.05).
Diversity 13 00494 g005
Figure 6. Geographical location of caves of Chapada de Canga and relations between genera determined in the caves. (a) A site map of Chapada de Canga caves. For more details see Figure 1. The map was obtained from “Instituto Prístino: atlas digital geoambiental. WebGis system (Web Geographical Information System): https://institutopristino.org.br/atlas, accessed on 15 June 2020. (b) Venn diagram representing the number of OTUs (on genus level) shared by the three caves of Chapada de Canga. (c,d) Venn diagrams representing unique and shared genera between each Chapada cave for floor and ceiling, respectively.
Figure 6. Geographical location of caves of Chapada de Canga and relations between genera determined in the caves. (a) A site map of Chapada de Canga caves. For more details see Figure 1. The map was obtained from “Instituto Prístino: atlas digital geoambiental. WebGis system (Web Geographical Information System): https://institutopristino.org.br/atlas, accessed on 15 June 2020. (b) Venn diagram representing the number of OTUs (on genus level) shared by the three caves of Chapada de Canga. (c,d) Venn diagrams representing unique and shared genera between each Chapada cave for floor and ceiling, respectively.
Diversity 13 00494 g006
Figure 7. Diversity and prediction of functional metabolism in cave 2. (a) Venn diagram showing relations between floor, ceiling, and biofilm hanging from the ceiling samples from cave 2 (Chapada de Canga), on phylum level. (b) Analysis of the functional prediction and metabolic pathways associated with the ceiling, floor, and biofilm of cave 2. Elevated abundances for each pathway are shown in green, reduced abundances in red. ANOVA followed by post hoc Tukey–Kramer test (p < 0.05), however no statistical difference was observed to the categories available.
Figure 7. Diversity and prediction of functional metabolism in cave 2. (a) Venn diagram showing relations between floor, ceiling, and biofilm hanging from the ceiling samples from cave 2 (Chapada de Canga), on phylum level. (b) Analysis of the functional prediction and metabolic pathways associated with the ceiling, floor, and biofilm of cave 2. Elevated abundances for each pathway are shown in green, reduced abundances in red. ANOVA followed by post hoc Tukey–Kramer test (p < 0.05), however no statistical difference was observed to the categories available.
Diversity 13 00494 g007
Table 1. Characterization of sampled caves.
Table 1. Characterization of sampled caves.
Cave NumberLocalizationGeomorphologyOther Features
CountyLocalityBiomeLatitude (S)Longitude (W)Altitude (m)LithologyCeilingFloorHuman ActivityProximity to Mine (up to 500 m)Protected Area
1MarianaChapadaAtlantic Forest20°9′53.78″43°24′19.49″879FerruginousCangaCangaNo~250 mNo
2MarianaChapadaAtlantic Forest20°9′54.03″43°24′29.72″865FerruginousCangaCangaNo~350 mNo
3MarianaChapadaAtlantic Forest20°9′49.15″43°24′21.85″880FerruginousCangaCangaNo~350 mNo
4CaetéGandarela’s MountainAtlantic Forest20°3′19.58″43°41′42.41″1624FerruginousCangaItabiriteNoNoYes
5Santa BárbaraGandarela’s MountainAtlantic Forest20°3′23.40″43°35′59.94″1236FerruginousCangaItabiriteNoNoYes
6Ouro PretoLavras NovasAtlantic Forest20°26′35.76″43°31′31.18″1480FerruginousCangaItabiriteNoNoNo
7Ouro PretoLavras NovasAtlantic Forest20°28′42.20″43°31′15.72″1346QuartiziteQuartiziteQuartiziteNoNoNo
8Nova LimaMoeda’s MountainAtlantic Forest20°13′18.26″43°58′38.48″1488FerruginousCangaItabiriteYesNoNo
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lemes, C.G.C.; Villa, M.M.; Felestrino, É.B.; Perucci, L.O.; Assis, R.A.B.; Cordeiro, I.F.; Fonseca, N.P.; Guerra, L.C.C.; Caneschi, W.L.; Moraes, L.Â.G.; et al. 16S rRNA Gene Amplicon Sequencing Data of the Iron Quadrangle Ferruginous Caves (Brazil) Shows the Importance of Conserving This Singular and Threatened Geosystem. Diversity 2021, 13, 494. https://0-doi-org.brum.beds.ac.uk/10.3390/d13100494

AMA Style

Lemes CGC, Villa MM, Felestrino ÉB, Perucci LO, Assis RAB, Cordeiro IF, Fonseca NP, Guerra LCC, Caneschi WL, Moraes LÂG, et al. 16S rRNA Gene Amplicon Sequencing Data of the Iron Quadrangle Ferruginous Caves (Brazil) Shows the Importance of Conserving This Singular and Threatened Geosystem. Diversity. 2021; 13(10):494. https://0-doi-org.brum.beds.ac.uk/10.3390/d13100494

Chicago/Turabian Style

Lemes, Camila G. C., Morghana M. Villa, Érica B. Felestrino, Luiza O. Perucci, Renata A. B. Assis, Isabella F. Cordeiro, Natasha P. Fonseca, Lara C. C. Guerra, Washington L. Caneschi, Lauro Â. G. Moraes, and et al. 2021. "16S rRNA Gene Amplicon Sequencing Data of the Iron Quadrangle Ferruginous Caves (Brazil) Shows the Importance of Conserving This Singular and Threatened Geosystem" Diversity 13, no. 10: 494. https://0-doi-org.brum.beds.ac.uk/10.3390/d13100494

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop