Original Research ARTICLE
Archaeal Communities in Deep Terrestrial Subsurface Underneath the Deccan Traps, India
- 1Environmental Microbiology and Genomics Laboratory, Department of Biotechnology, Indian Institute of Technology Kharagpur, Kharagpur, India
- 2School of Bioscience, Indian Institute of Technology Kharagpur, Kharagpur, India
- 3School of Environmental Science and Engineering, Indian Institute of Technology Kharagpur, Kharagpur, India
- 4Department of Geology and Geophysics, Indian Institute of Technology Kharagpur, Kharagpur, India
- 5Ministry of Earth Sciences, Borehole Geophysics Research Laboratory, Karad, India
- 6CSIR-National Geophysical Research Institute, Hyderabad, India
Archaeal community structure and potential functions within the deep, aphotic, oligotrophic, hot, igneous provinces of ∼65 Myr old basalt and its Archean granitic basement was explored through archaeal 16S rRNA gene amplicon sequencing from extracted environmental DNA of rocks. Rock core samples from three distinct horizons, basaltic (BS), transition (weathered granites) (TZ) and granitic (GR) showed limited organic carbon (4–48 mg/kg) and varied concentrations (<1.0–5000 mg/kg) of sulfate, nitrate, nitrite, iron and metal oxides. Quantitative PCR estimated the presence of nearly 103–104 archaeal cells per gram of rock. Archaeal communities within BS and GR horizons were distinct. The absence of any common OTU across the samples indicated restricted dispersal of archaeal cells. Younger, relatively organic carbon- and Fe2O3-rich BS rocks harbor Euryarchaeota, along with varied proportions of Thaumarchaeota and Crenarchaeota. Extreme acid loving, thermotolerant sulfur respiring Thermoplasmataceae, heterotrophic, ferrous-/H-sulfide oxidizing Ferroplasmaceae and Halobacteriaceae were more abundant and closely interrelated within BS rocks. Samples from the GR horizon represent a unique composition with higher proportions of Thaumarchaeota and uneven distribution of Euryarchaeota and Bathyarchaeota affiliated to Methanomicrobia, SAGMCG-1, FHMa11 terrestrial group, AK59 and unclassified taxa. Acetoclastic methanogenic Methanomicrobia, autotrophic SAGMCG-1 and MCG of Thaumarcheaota could be identified as the signature groups within the organic carbon lean GR horizon. Sulfur-oxidizing Sulfolobaceae was relatively more abundant in sulfate-rich amygdaloidal basalt and migmatitic gneiss samples. Methane-oxidizing ANME-3 populations were found to be ubiquitous, but their abundance varied greatly between the analyzed samples. Changes in diversity pattern among the BS and GR horizons highlighted the significance of local rock geochemistry, particularly the availability of organic carbon, Fe2O3 and other nutrients as well as physical constraints (temperature and pressure) in a niche-specific colonization of extremophilic archaeal communities. The study provided the first deep sequencing-based illustration of an intricate association between diverse extremophilic groups (acidophile-halophile-methanogenic), capable of sulfur/iron/methane metabolism and thus shed new light on their potential role in biogeochemical cycles and energy flow in deep biosphere hosted by hot, oligotrophic igneous crust.
The Earth’s deep continental biosphere is hot, oligotrophic, and limited by several physical and chemical constraints for life; yet it harbors up to 19% of the Earth’s biomass (Whitman et al., 1998; McMahon and Parnell, 2014; Kieft, 2016). Exploration of microbial life in such extreme habitat represents one of the most intriguing sectors with enormous opportunities for reconnoitering the evolution and adaptation of life on our planet and beyond (Colwell and D’Hondt, 2013). In spite of considerable research attention, the deep biosphere is still counted among the “least known frontiers of Earth” and our effort to generate a global census on deep microbial life remains limited by habitat accessibility (Breuker et al., 2011; Purkamo et al., 2018).
Despite the harsh conditions in terms of lack of readily utilizable energy and nutrient sources, extreme pressure and temperature, limited space and fluid mobility, microbial life has been observed across diverse deep crystalline biosphere settings (Fredrickson and Balkwill, 2006; Edwards et al., 2012; Miettinen et al., 2015; Kieft, 2016). It is interesting to note that the inhabitant microorganisms of such extreme habitats not only evolve striking adaptive strategies to sustain their living, but may also require those conditions for their survival (Rampelotto, 2013; Kieft, 2016). These organisms are often referred as extremophiles and include members of all three domains of life, i.e., Bacteria, Archaea, and Eukarya, of which archaea represents a small but important proportion. Nevertheless, they are proficient in adapting to different extreme conditions including survival under chronic energy stress and holding different “extremophily records” (Swan and Valentine, 2001; Rampelotto, 2013) and are the important contributors to subterranean biogeochemical cycles (Biddle et al., 2006; Silver et al., 2012; Offre et al., 2013; Magnabosco et al., 2016). Archaea are among the most hyperthermophilic, acidophilic, alkaliphilic, and halophilic microorganisms known. Previous studies have provided valuable information on the distribution and composition of archaea within deep igneous terrestrial subsurfaces and their role in carbon, nitrogen and sulfur geocycling (Offre et al., 2013; Ino et al., 2018). Archaea mediated subsurface H2 cycling, methanogenesis, anaerobic methane oxidation, iron and sulfur redox transformation, and nitrate reduction are crucial components of subsurface lithoautotrophic microbial ecosystem (SLiME) and/or mixotrophic ecosystems (Stevens and McKinley, 2000; Nealson et al., 2005; Gregory et al., 2019; Smith et al., 2019). Because of their extremophilic nature, significance in different biogeochemical cycles and syntrophic relations with bacterial members, archaea are considered as suitable residents of deep subterranean environments.
However, compared to several investigations highlighting bacterial diversity, focused studies on archaea in deep terrestrial subsurface are very limited (Takai et al., 2001; O’Connell et al., 2003; Waldron et al., 2007; Kieft, 2016; Lazar et al., 2017; Purkamo et al., 2018). One of the first studies which focused solely on archaea in deep terrestrial subsurface was carried out in South African gold mines (Takai et al., 2001). The study revealed “great” phylogenetic diversity and highlighted novel lineages of the phyla Crenarchaeota (SAGMCG-1 and SAGMCG-2) and Euryarchaeota (SAGMEG-1 and SAGMEG-2). More recently, SAGMCG was reclassified into the phylum Thaumarchaeota (Brochier-Armanet et al., 2008). The presence of niche-specific archaeal communities was observed within subsurface organic rich sedimentary rocks of Antrim Shale, United States (Waldron et al., 2007; Wuchter et al., 2013). The prevalence of lithoautotrophic archaea within carbonate-rock/siliciclastic-rock aquifer systems and their restricted dispersal from the surface to the subsurface environment was reported by Lazar et al. (2017). The widespread occurrence of methanogenic and anaerobic methane-oxidizing archaea (ANME); members of more versatile Thermoplasmatales are reported from deep granitic environment of the Witwatersrand deep subsurface, Chelungpu fault, Gifu prefecture, and from Precambrian bedrock fluids in Fennoscandian Shield (Takai et al., 2001; Moser et al., 2005; Gihring et al., 2006; Bomberg et al., 2015; Wu et al., 2015; Ino et al., 2016; Purkamo et al., 2018 and references within). Methanogenesis and methane cycling processes are a major contribution of archaea in deep crystalline rock biosphere, which are considered to be closely interconnected with redox transformation of electron acceptors (nitrate, sulfate, metal oxides, and humics) or electron donors (hydrogen, ammonia, nitrite, and sulfide) (Kietäväinen and Purkamo, 2015; Ino et al., 2018; Smith et al., 2019) and thus play significant roles in overall community function. A few studies have profiled the archaeal diversity in the deep terrestrial biosphere, but their distribution and metabolic significance in biogeochemical cycles within deep subsurface igneous provinces remains poorly explored. The hard rock terrestrial deep biosphere in basalt, metabasalt and granitic horizons have so far mainly been explored by analyzing the groundwater; studies of rock samples recovered through scientific deep drilling are rare (Breuker et al., 2011 and references within; Hubalek et al., 2016; Kieft, 2016; Ino et al., 2018).
Here, we aim to explore the archaeal communities hosted by ∼65 Ma old basalts and underlying Archean basement granitoids comprising the deep terrestrial subsurface igneous province of the Deccan traps province using the cores recovered through scientific deep drilling (Roy et al., 2013; Gupta et al., 2015). The Deccan traps are a massive continental flood basalt province (total area > 0.5 million km2) resting on ∼2.5 Ga Archean crystalline basement (Sen, 2001; Schoene et al., 2015; Bhaskar Rao et al., 2017). The Koyna–Warna region of Deccan traps is known for recurrent reservoir triggered seismicity following the impoundment of the Koyna dam in 1962 (Gupta, 1992, 2017). Exploratory deep drilling at multiple sites in the Koyna seismogenic zone allowed us to access the deep subsurface rock cores covering the entire thickness (varying between 500 and 1250 m) of Deccan basalt (designated as BS), a few 100 m of the underlying granitic basement (GR) down to 1400 m depth and an intermediate region of weathered granites (referred as transition zone or TZ) between BS and GR horizons. The goal of this study is to shed light on the potential archaeal communities across the different depths of deep basaltic-granitic subsurface. This work provides a detailed information on archaeal life in deep igneous terrestrial subsurface where extreme temperature and pressure regimes (on average 25°C increase per 1000 m in basalt and 15°C increase per 1000 m in granite; 26.7 MPa increase in lithostatic pressure per 1000 m) and geochemical conditions have shaped the environment in a unique way (Roy and Rao, 2000; Dutta et al., 2018).
Materials and Methods
Subsurface core samples were collected from different depths of two exploratory drill holes at Ukhalu (N 17° 07.552′, E 073° 52.148′) and Phansavale (N17° 09.017′, E073° 40.058′) of Koyna–Warna region of Deccan Traps, Maharashtra during July 2014 and May 2015. Eleven samples were obtained to cover the three major horizons: Deccan basalt (BS) (∼65 Ma), the underlying granitic basement rock (GR) (∼2500 Ma), and the short section between the two formations representing the weathered surface of granitic basement affected by the first lava flows (referred here as a transition zone, TZ) (see Supplementary Table S1). Samples were checked for the presence of sodium fluorescein (500 mg m–3) used during drilling (Nyyssönen et al., 2014). Only interior pieces of rock were selected for microbiological study to avoid potential drilling mud contamination as recommended elsewhere (Lever et al., 2013). Samples were collected following aseptic techniques and stored in sterile containers at 4°C for shipment. In the laboratory, the samples were stored at −20°C to limit the chance of microbial contamination until further analysis.
Rock cores were washed thoroughly under sterile condition with autoclaved, DNA free water (Thermo Fisher Scientific®) and sub-coring of the rock was done to get rock samples devoid of any possible contamination that might have occurred during drilling. At first, the uneven ends of the cores were removed using granite cutter fitted with sterilized cutting blade. Subsequently, with either a mechanical drilling instrument fitted with sterilized drill bit or with sterilized chisel and hammer, rock powder was obtained from the interior core. Sub-coring was done in different layers and sodium fluorescein was checked in the sub-cored rock powder. Rock powders were dissolved in water and aqueous fraction was used for detection of fluorescein using spectrophotometric analysis (Mota et al., 1991). The rock powders free from sodium fluorescein were stored in sterilized DNA free containers at −20°C for further geochemical and microbial analysis.
Rock powders were sieved through 2 mm mesh. Acid digestion of these rock powders was conducted using HNO3 and HF (3:1) in a microwave digester (Milestone SK12) for efficient extraction of the elements from rock powder to aqueous media followed by the measurement of major elements using Inductively Coupled Plasma Mass Spectrometry (iCAP Q, Thermo Fisher Scientific). Rock powders were ultrasonicated for major anions (Cl–, NO2–, SO42–, NO3–, and PO43–) using Dionex ICS2100 (Thermo Fisher Scientific). Total organic carbon (TOC), total inorganic carbon (TIC) and total carbon (TC) were measured using an OI analytical TOC analyzer. Before analyzing TOC, rock powders were fumigated with 1N HCl for 48 h under a fume hood. Quantification of the elemental oxides was done using a PANalytical Epsilon3 XRF instrument. Alkalinity of the rocks was measured using the USEPA method 310.2. The pH of the rock samples (in the presence of water) was measured by incubating 1 g of rock powder in 0.1 M CaCl2 at 1:10 ratio (w/v), and highly sensitive probes (Orion) fitted with an Orion multi-meter (Thermo Electron Corporation, Beverly, MA) were used (Islam et al., 2014). Four out of the eleven rock core samples used in this study (two each from borehole; U9, U8, PV4, and PV2) were used previously for a geomicrobiological investigation (Dutta et al., 2018). Geochemical properties of these four samples were included in this study for a better comparison (along with seven other samples).
Extraction of Environmental DNA
Community DNA was extracted from 11 samples using the MoBio power soil DNA isolation Kit (MoBio) following the manufacturer’s protocol. Total DNA from reagent control was also extracted using the same procedure as mentioned earlier, and were used subsequently to check any possible contamination. The quality (purity) and quantity of the extracted environmental DNA was first determined using a NanoDrop 2000 spectrophotometer, followed by precise fluorometric quantitation using Qubit (Thermo Fisher Scientific).
Quantitative Polymerase Chain Reaction (qPCR)
Quantification of the archaeal populations in the rock samples were performed by estimating the copy number of archaeal specific 16S rRNA gene. The abundance of methanogenic populations was estimated from mcrA copy numbers using quantitative PCR (qPCR). Details of qPCR primers used are provided in Supplementary Table S2. Two μl of extracted environmental DNA were added to the PCR mastermix with a total volume of 10 μl. All the reactions were set in triplicates. Quant Studio 5 was used to perform qPCR with Power SYBR green PCR mastermix (Invitrogen), with primer concentration of 5 pM and the following amplification conditions: 95°C for 10 min, 40 cycles of 95°C for 15 s, 55°C for 30 s and 72°C for 30 s. Melting curve analysis was run after each assay to check PCR specificity. Genes encoding archaeal 16S rRNA and mcrA were PCR amplified from the extracted environmental DNAs, cloned in TA cloning vector and 10-fold standard dilution series of plasmid DNAs with 102–1010 copy numbers of 16S rRNA and mcrA were used to quantify both genes in the environmental samples.
16S rRNA Gene Amplification and Sequencing
Following quality checks, the extracted environmental DNAs were sent to the Marine Biological Laboratory, Woods Hole, United States for Illumina MiSeq sequencing using the procedures outlined on the VAMP website1. Archaeal V4-V5 hypervariable region of 16S rRNA gene was sequenced on Illumina MiSeq platform. The primers used for amplification of archaeal V4-V5 region are 517F and 958R (Supplementary Table S2). The sequence reads were submitted to short read archive under SRA accession: SRP126374.
Paired-end reads obtained from Illumina MiSeq sequencing were merged into single end reads using FLASH (Magoč and Salzberg, 2011) with a minimum and maximum overlap of 10 and 100, respectively and a mismatch density of 0.25. This step was followed by quality filtering using the quantitative insights into microbial ecology (QIIME 1.9.1) (Caporaso et al., 2010). Sequences having lengths outside bounds of 250 and 500, mean quality scores below a minimum of 25 and maximum homopolymer runs exceeding a limit of 6 were filtered out. De novo-based clustering of reads to form OTU was performed using UCLUST under QIIME workflow. Sequences having greater than 97% similarity were assigned to the same OTU. Representative reads from each OTU were assigned taxonomy using UCLUST trained SILVA 128 database (Quast et al., 2012). OTUs which were present in the reagent controls were removed from OTU pool of the samples. Alpha diversity parameters were calculated using alpha_diversity.py under QIIME workflow.
Selected geochemical parameters across 11 samples were used to perform principal component analysis (PCA) using PAST software (Hammer et al., 2001). Geochemical parameters were normalized on the basis of feature scaling before PCA was plotted. Archaeal families from across the samples were selected for Spearman correlation analysis and correlation heat-map was constructed using METAGENassist (Arndt et al., 2012). The abundant archaeal classes were used to construct PCA biplot using PAST software. Non-metric multidimentional scaling (NMDS) analysis of Bray–Curtis metric based on all the archaeal families was performed using PAST software, and 21 different physico-chemical parameters were declared as environmental variables. Similarity percentage (SIMPER) analysis was also done using PAST software for 11 samples to identify the microbial families responsible for the grouping of different rock types. OTU overlap among the samples from similar rock type and OTU pool of different rock types was elucidated using InteractiVenn (Heberle et al., 2015).
Co-occurrence network analysis was performed to ascertain major connections among archaeal populations and different geochemical parameters of the samples. Correlation among different archaeal genera with geochemical parameters was calculated using the otu.association command in mothur (Schloss et al., 2009). Pairwise Spearman correlations with correlation significance value <0.05 was used for the construction of co-occurrence networks using Cytoscape 3.4.0 (Shannon et al., 2003). Two separate networks were constructed for positive and negative correlations.
Prediction of Functionality of Archaeal Communities
Metagenomic inventories of the archaeal communities were predicted using PICRUSt (Phylogenetic Investigation of Communities by Reconstruction of Unobserved States) (Langille et al., 2013). Closed referenced based OTU picking using Greengene 13.5 as a reference database was performed in QIIME, and PICRUSt was used to allocate the functional attributes by comparing the identified 16S rRNA gene sequence with the nearest match of the known genome sequence. A weighted nearest sequenced taxon index (NSTI) was calculated to evaluate the accuracy of prediction of the metagenomes. A PCA biplot was constructed on the basis of relative abundance of different metabolisms to understand the relatedness among different archaeal communities on the basis of their community functions.
Geochemical Characteristics of the Samples
Rock core samples obtained from varying depths across the basaltic layers of Deccan traps and underlying granitic basement were analyzed for major geochemical properties. Geochemical parameters pertinent to microbial activities were plotted (Figure 1A and Supplementary Table S3). The samples showed a characteristic trend in their geochemical properties with increasing depth. Based on these parameters, the distinct geochemical nature of granitic and basaltic horizons was evident, whereas samples from the transition zone showed a mixed nature. Concentrations of major anions like NO3–, and SO42–; metals like Mg, Fe, Na, Cr, Cd, and oxides of Si and K increased with depth. Samples from shallower levels were distinctive based on slightly acidic pH, higher TOC, TIC, TC, CaO, TiO2, and Fe2O3. Although certain unexpected variations in geochemical parameters across different horizons were observed, PCA showed that the samples from two different (BS and GR-TZ) zones were geochemically distinctive with respect to the measured geochemical parameters (Figure 1B). Concentrations of TOC, TIC, and TC decreased with increasing depth and were generally lower in GR horizon compared to BS horizon samples. NO3– and NO2– were overall on the lower side when compared to SO42– whereas PO43– was below detection limit for all the samples except sample PV8.
Figure 1. Distribution of different geochemical parameters and archaeal population across different samples of subterranean Deccan traps. (A) Variations of different geochemical parameters across samples. (B) Principal component analysis (PCA) of geochemical parameters of subsurface rock cores of Koyna–Warna region; Green dots represent samples from the granitic horizon, red circles represent samples from the basaltic horizon and blue dot represents sample from the transition zone. mbs – meters below surface.
Archaeal Abundance in Deccan Subsurface
Archaeal 16S rRNA gene and mcrA gene copy numbers in the analyzed rock samples were quantified using real time PCR. 16S rRNA- and mcrA-gene copy numbers could be successfully obtained for four samples (three BS and one GR samples), while for the remainder of the samples, it was not successful even after repeated attempts (Table 1). The 16S rRNA gene copies varied between 2.14 × 103 and 4.75 × 104 per gram of rock (with an average copy number of 1.6 × 104/gram of rock). Assuming an average occurrence of 1.7 copies of the 16S rRNA gene in archaeal genomes (Stoddard et al., 2015), we estimated an archaeal cell abundance ranging from 1.25 × 103 to 2.78 × 104 per gram of rock. Gene fragment of mcrA was observed in the same four samples and its copy number varied from 4.1 × 101 to 1.6 × 103 (with average copy number of 5.18 × 102). Percentage of mcrA copy number with respect to archaeal cell abundance varied between 3 and 8% and was observed to be higher in GR sample (8%) over the BS samples (average 4%; Standard deviation 1%).
Sequence information and diversity indices are summarized in Table 1. Illumina MiSeq based sequencing of the hypervariable V4-V5 region of archaeal 16S rRNA gene produced 32509 numbers of total merged reads for 11 samples. After filtering and removing potential erroneous reads, 26739 reads were used for OTU picking. A maximum of 6079 and a minimum of 693 filtered reads were obtained for PV5 and U3, respectively. OTUs containing reads from reagent (DNA extraction) control were removed from the OTU pool. A total of 20570 reads grouped into 3148 OTUs were finally obtained with a maximum of 549 (PV5) and a minimum of 140 (PV8) OTUs. Out of these total OTUs, only 6.5% were unassigned to any taxa. Noticeably, 2833 unique OTUs (∼90% of the total OTUs) were found across the samples and these OTUs represent 31% of the total reads. Interestingly, the number of observed OTUs increased with depth in basalts followed by a reduced copy number in the transition zone. Similarly, the number of OTUs increased with depth in granites except for sample U4 (Figure 2A). A comparable trend was observed in the Chao1 index used to estimate the minimal number of OTUs present. The Chao1 index suggests the presence of a maximum of 3056 OTUs and a minimum of 379 OTUs in PV5 and PV8, respectively (Figure 2B). Simpson’s (measures species dominance) and Shannon’s (evaluates species abundance and evenness) indices, ranged from 0.79 to 0.99 and 3.75 to 7.22 and both increased with depth (Figures 2C,D). Overall, the average Shannon diversity index was much higher in granites than in basalts indicating higher species richness in granites.
Figure 2. Alpha diversity parameters of different samples. (A) Observed OTUs, (B) Chao1 estimator, (C) Simpson Index, and (D) Shannon Index.
Archaeal Community Analysis
16S rRNA gene sequences assigned to archaea covered 9 phyla, 12 classes, 7 orders, 20 families, and 40 genera. Euryarchaeota was the most abundant phylum (average abundance of 59%) followed by Thaumarchaeota (30%) and Bathyarchaeota (5%) (Figure 3A). Crenarchaeota, Lokiarchaeota and WSA2 were sporadically observed across different samples. Aigarchaeota and pMC2A209 were only found in sample PV2 (BS) whereas Marine Hydrothermal Vent Group (MHVG) was observed in sample U14 (GR). A comparison between BS and GR samples revealed that the former had considerably higher proportion of Euryarchaeota (average 80%), but lower proportion of Thaumarchaeota (average 12%) and Crenarchaeota (average 2%). In GR, although the members of Euryarchaeota were detected at moderately high abundance (31%), Thaumarchaeota outnumbered and represented the most abundant taxa (54%). Bathyarchaeota present as a minor group in BS (2.6%), constituted 8.7% in GR.
Figure 3. Distribution of archaeal population across different samples of subterranean Deccan traps. Archaeal community composition of Deccan subsurface at (A) phylum level and (B) class level.
On the basis of cumulative abundance, Thermoplasmata, Methanomicrobia, South African Gold Mine Gp 1(SAGMCG-1), FHMa11 terrestrial group and AK59 were identified as the five major archaeal classes (Figure 3B). Thermoplasmata (Euryarchaeota) was much more abundant in BS- (62%) than GR-cores (13%). In contrast, average abundance of three Thaumarchaeota members namely the SAGMCG-1, FHMa11 terrestrial group and AK59; as well as Methanomicrobia (Euryarchaeota) were more abundant in the GR cores (24, 16, 13, and 18%, respectively) than in BS samples (6, 6, 0, and 10%, respectively) (Figure 4). Other archaeal classes detected in the cores recovered from the Deccan subsurface were Methanobacteria, Thermoprotei, Halobacteria, Soil Crenarchaeotic Group (SCG), and Terrestrial Hot Spring Gp (THSCG). Among these, Methanobacteria (Euryarchaeota) and Thermoprotei (Crenarchaeota) were relatively more abundant in BS (average abundance of 6 and 3%, respectively) compared to GR (0.4 and 0.7%). Other archaeal classes detected in this study were present only in the BS samples. Significant presence of Methanomicrobia (53% abundance), Methanobacteria (22%) and Thermoplasmata (13%) was observed in TZ core. Halobacteria (Euryarchaeota) and SCG (Thaumarchaeota) furthermore comprised more than 1% of the total reads in the TZ core. Methanogenic archaea were present across the samples and horizons with slightly higher average abundance in GR (18%, BS: 16%). Four out of seven known methanogenic archaeal orders were present in the order of highest to lowest abundance in GR: Methanosarcinales, Methanomicrobiales, Methanobacteriales, and Methanocellales. This order of abundance was different in the BS horizon where Methanobacteriales was the most abundant followed by Methanomicrobiales, Methanosarcinales, and Methanocellale. In general, samples from the Phansavale site harbored more methanogenic taxa (particularly for Methanomicrobia and Methanobacteria).
Figure 4. Average archaeal community composition of (A) basalt, (B) transition, and (C) granite at class level. Since only one sample from the transition zone was selected, percentage abundance of the archaeal community is shown.
Taxonomic affiliation at the lowest level (genus) was possible only for 62% of total classified reads (44% of the OTUs). Thermoplasma followed by Methanosaeta, Ferroplasma, Methanobacterium, and Methanolacinia represented the five most abundant genera across the samples. Interestingly all these major genera belonged to the phylum Euryarchaeota. Other archaeal genera detected in greater than 1% cumulative abundance were Methanolobus, ANME-3, Methanolinea, Methanocalculus, Candidatus Nitrososphaera, Natronorubrum, Methanosarcina, Methanocorpusculum, Methanospirillum, Methanogenium, Methanoregula, and Halovivax. The distribution of the major genera showed a distinct preference for either of the horizons, although a few, like ANME-3 and Terrestrial Miscellaneous Gp (TMEG) were present in a greater number of samples with varying abundance (0–4.5 and 0–1.6%, respectively). Higher average abundance of Thermoplasma (35%), Ferroplasma (21%), Methanobacterium (6%) and Methanolacinia (4%) was observed in BS. In GR horizon, Methanosaeta was more dominant (14%, compared to 3% in BS) along with several unclassified taxa affiliated to Bathyarchaeota, AK59, FHMa11 terrestrial group and SAGMAG-1 (average abundance: 9–24%).
Endemic Archaeal Population in Different Rock Types
Out of a total of 3148 OTUs only four core OTUs (OTUs present in all samples) were observed across all the BS samples whereas only one core OTU was observed for the GR horizon (Figure 5). Core OTUs of BS and GR represent 12.59 and 0.13% of total archaeal reads, respectively. Of the four core OTUs of BS samples, three were affiliated to Ferroplasma and one to the genus Thermoplasma. The single core OTU detected in GR horizon was affiliated to SAGMCG-1. No shared OTUs were observed across all the samples from BS and GR horizons. From a total pool of 2900 OTUs observed across all the BS and GR samples, only 3.5% of the OTUs were present in both the horizons whereas 55% and 41% of OTUs remained unique to BS and GR horizons (Supplementary Figure S1).
Figure 5. Venn diagram of OTUs present in (A) basaltic and (B) granitic zones. OTUs were based on 97% sequence similarity. The numbers inside the diagram are the number of OTUs.
Correlation Studies and Deciphering the Niche-Specific Colonization of Archaeal Taxa
A set of statistical tools were used to investigate the nature of correlation and the niche-specific colonization of archaeal populations within the BS and GR horizons. A PCA biplot was drawn on the basis of abundance of different archaeal classes which showed distinct assemblages of archaea across BS and GR horizons. As evident in Figure 6, regardless of the borehole sites, subsurface archaeal communities partitioned into (i) BS zone guild and (ii) GR zone guilds. Close association among Thermoplasmata, Thermoprotei, Methanobacteria and Halobacteria members in BS cores, and SAGMCG-1, FHMa11 terrestrial group, AK59 and an unclassified class of Bathyarchaeota (which represented the dominant archaeal assemblage within the GR horizon) in GR cores was observed. Sample from transition zone (PV7) did not group with any of the other groups and remained distinctly placed. One of the BS samples, PV4, clustered with the rest of the BS samples, but showed relatively higher similarity with the GR communities.
Figure 6. PCA analysis on the basis of the abundance of archaeal classes present among different samples. Red symbols represent basaltic samples; green symbol represents transition zone sample; blue symbol represents granitic samples.
Correlation among the archaeal populations was identified through the Spearman correlation method. Based on correlation values, two distinct groups/clusters each consisting of 9 and 10 archaeal families were observed (Figure 7A). Cluster 1 was exclusively represented by Euryarchaeota members and could be divided in to two distinct sub-clusters. One of the sub-clusters consisted of Methanosaetaceae, Methanoregulaceae, Marine Benthic Group D and DHVEG-1, Methanomicrobiaceae and Methanosarcinaceae whereas the other sub-cluster comprised of CCA47, Methanocorpusculaceae, Methanospirillaceae and TMEG. Most of the families in this cluster are known methanogens (both hydrogenotrophic and acetoclastic types). Interestingly, these families were more abundant in GR rock cores with relatively low contributions in BS samples (Figure 7B). Cluster 2 consisted of 10 families, of which nine of the members were affiliated to Euryarchaeaota. This cluster comprised of two sub-clusters of which one consisted of Methanomicrobiales Incertae Sedis, Halobacteriaceae, Methanobacteriaceae, KTK 4A, Methanocellaceae, AMOS1A-4113-D04, and Ferroplasmaceae. The second sub-cluster represented a strong relation among Sulfolobaceae, Thermoplasmataceae, and BSLdp215. This group consists of archaeal families known to be thermophilic/thermotolerant, methanogens, and participate in the sulfur cycle. Their presence was mainly observed in BS cores (Figure 7B).
Figure 7. Correlation among different microbial classes. (A) Heatmap of Spearman correlation across archaeal families. (B) Heatmap of relative percentage abundance of corresponding archaeal classes across all the samples.
In order to assess the collective role of environmental variables on archaeal populations inhabiting the BS and GR horizons NMDS (using archaeal families and set of environmental parameters) based analysis was performed; which separated the samples into two distinct abundance-weighted categories (Figure 8). Vectors in NMDS displayed the major environmental parameters which might govern the unique archaeal assemblages in BS and GR horizons. SIMPER analysis was used to determine the archaeal taxa responsible for the differences observed between the BS and GR archaeal microbiomes (Supplementary Table S4). The topmost “discriminating” taxa were Thermoplasmataceae, Ferroplamaceae, and Methanosaetaceae. Average taxonomic distribution patterns demonstrated that BS microbiomes contained a higher relative abundance of Thermoplasmataceae, Ferroplasmaceae, Methanobacteriaceae, Methanomicrobiaceae, Sulfolobaceaea, AMOS1A-4113-D04, and Halobacteriaceae compared to GR. On the other hand, GR microbiome displayed greater relative abundance of Methanosaetaceae, Marine Benthic Group D and DHVEG-1 and unassigned families of SAGMCG-1, AK59, FHMa11 terrestrial group and Bathyarchaeota.
Figure 8. Non-metric multidimensional scaling (NMDS) analyses of microbial communities and its association with environmental factors of subsurface rock cores of Koyna–Warna region (Bray–Curtis, Stress = 0.12). Red symbols represent basaltic samples; green symbol represents transition zone samples; blue symbols represent granitic samples.
Archaeal Co-occurrence and Effect of Geochemical Parameters
A network analysis was done to better understand the co-occurrence among different environmental parameters and archaeal taxa in Deccan subsurface. Networks for positively and negatively co-related archaeal genera with different environmental parameters were separately constructed. The positively correlated network consisted of 58 nodes and 178 edges (Figure 9A) whereas the negatively correlated network had fewer nodes (33) and edges (62) (Figure 9B). The positively correlated network consisted of five sub-networks. The main sub-network consists of a hotspot (hotspots are region where large numbers of connections are observed among nodes) constituted by 14 members. This hotspot was mainly dominated by halotolerant archaea viz. Halovivax, Haloterrigena, Natrorubrum, and others. This cluster also includes members of acidophilic Acidiplasma, methanogenic Methanobacterium and Methanocalculus and ammonia-oxidizing Candidatus Nitrososphaera. Though most of these archaeal members are present in very low abundances, it is interesting to note that these organisms were mainly present in higher depth basalts (PV4, PV2, and PV5) and transition zone of the Phansavale borehole. The peripheral nodes of the hotspot consisted of archaeal genera observed mainly in the BS cores. Of these peripheral nodes and edges, one of the significant connections was between acidophilic iron-oxidizing Ferroplasma and Fe2O3. Sub-network consisting of Methanosaeta, Methanolinea, and ANME-3 also included other geological parameters, which had an increasing trend with depth. Connections between Methanolinea, methanotrophic ANME-3 and other geochemical parameters was established by Methanosaeta, which harbors versatile capability of acetoclastic methanogenesis and is also known for its methylotrophic nature. Sub-networks comprising of Methanoregula, Methanogenium, and Methanosarcina were correlated with elevated concentrations of SO42– and NO3–. We also observed two micro-networks consisting of one connection each. One of the micro-networks consisted of sulfur-oxidizing Metallosphaera and Acidianus which were mainly found in amygdaloidal basalt (U8), which displayed elevated SO42– concentration while compared to other basalts. The second micro-network consisted of two hydrogenotrophic methanogens Methanocella and Methanocorpusculum whose co-presence was observed in the deepest recovered granitic core (U14).
Figure 9. Co-occurrence network among archaeal genera and geochemical parameters in subsurface rock cores of Koyna–Warna region. The color of the nodes represent the number of degrees it has. (A) Network showing positive correlation. (B) Network showing negative correlation.
The network showing negative correlation consisted of 33 nodes (Figure 9B). Alkalinity was the node having highest degree (13) and was connected to the groups which were mainly present in the hotspot of the positively correlated network. Ferrolpasma was negatively correlated with the geochemical parameters (Mg, SiO2, K2O, and NO2–) which were present in higher abundances in most of the granites. Interestingly, Methanosaeta, ANME-3 and Methanolinea were found to be negatively correlated with TOC.
Predicted Metabolic Properties
The metabolic properties of the archaeal communities hosted by the subsurface rocks from different horizons were predicted from 16S rRNA gene sequences using PICRUSt (Table 2). The NSTI was calculated for all the samples to evaluate the accuracy of prediction of the metabolic properties (Supplementary Table S5). The NSTI values varied between 0.05 (predicted microbiome for PV2, BS core) and 0.32 (predicted microbiome for U4, GR core). The most abundant metabolism displayed by the majority of the archaeal groups were amino acid metabolism (19–25% of all metabolism genes), carbohydrate metabolism (20–24%) and energy metabolism (13–18%). Genes allocated for amino acid metabolism and carbohydrate metabolism were elevated in GR samples (average 25 and 23%, respectively) compared to BS samples (average 22% for both) (Table 2). The PCA biplot was constructed on the basis of the relative abundance of different metabolisms to understand the grouping of the archaeal community function across different rock samples on the basis of predicted community function (Figure 10). Except for sample PV4 (a BS core recovered from a depth of 106.4 m below surface), all the basalts and granites formed distinct clusters. Predicted functions of the archaeal community hosted by sample PV4 were similar to the predicted metabolisms of the GR which was evident from the PCA. PV7 (TZ) did not cluster with BS or GS group showing distinct community metabolism. Similar grouping was observed in previous PCA and NMDS analysis which was carried out on the basis of archaeal classes and families, respectively. A similar trend obtained from different statistical analysis (on the basis of taxonomy and predicted community functions) (Figures 6, 8, 10) implies that taxonomy of different archaeal groups in particular microbiomes reflects the archaeal community function.
Figure 10. PCA biplot on the basis of the most abundant metabolism related to archaeal genes in the predicted metagenomes; red symbol depicts basaltic samples; blue symbol depicts granitic samples and green symbol depicts transition zone sample.
Deep subsurface of the Deccan traps province may be considered as an extreme environment for microbial life with several physical as well as chemical stressors. On average, temperature increases with depth at a rate of ∼25°C per km in the Deccan basalt formation and ∼15°C per km in the granitic bedrock (Roy and Rao, 1999, 2000). Typically, temperatures at a depth of 1 km are in the range 40–45°C and at 1.5 km, 55–60°C in the study region. Lithostatic pressure increases at a rate of ∼26.7 MPa per km. Nevertheless, the archaeal abundance in this extreme habitat is estimated to be in the range of 103–104 cells per gram of rock, well within the range (102–105 archaeal 16S rRNA gene copies per gram) reported earlier from several deep terrestrial subsurface ecosystems (Breuker et al., 2011; Bomberg et al., 2015; Purkamo et al., 2016). The broad geochemical spectrum of the rock cores portrays critical details of the abiotic components of this unique ecosystem and provides the perspective to understand its microbial ecology more meaningfully.
The biological, geochemical and hydrological processes at work within the deep terrestrial habitats in the Deccan traps and underlying basement rocks have likely influenced the archaeal community structure. Together with their gneissic/migmatitic nature, these rocks exhibit structural deformation and associated fractures, which could be connected to the seismicity of this region (Misra et al., 2017). Based on their geochemical properties, the three classes of rock samples could be distinctly separated; particularly the basaltic and granitic provinces. Although the overall Deccan subsurface was low in organic carbon content, relatively higher levels detected in the BS horizon could be linked to intrusion of surface organic matter, or remnants of the ancient organic matter derived from diverse past life forms present within the different layers of lava flows. The deeper GR horizon is more depleted in organic carbon, and this is in line with the extremely low porosity, lack of fluid mobility and perhaps low biological activities that might produce any “in house” organic matter (Goswami et al., 2017). Granitic rocks are formed deep within the continental crust by the cooling of intruding magma, which are often devoid of organic carbon. Unlike the scarcity of organic carbon, minerals present in deep igneous rocks are an important source of electron acceptors (iron, sulfate, nitrate/nitrite) available to the microorganisms present (Sen, 2014). Availability of sulfate in the deeper horizons (>1000 m) is widely reported through geochemical analysis of deep granitic aquifers around the globe (Ino et al., 2018 and references therein) and this corroborates well with the present study. Sulfate dependent oxidation of methane and other reduced substrates (e.g., hydrogen, ammonia, nitrite) have been identified as a potent metabolic driver in oligotrophic deep continental subsurface (Nealson et al., 2005). Scarcity of organic carbon but relative abundance of iron (as Fe2O3), oxidized nitrogen compounds and sulfate in the deep subsurface indicate the broad chemical landscape of this extreme igneous province that could potentially support a versatile microbial catabolic repertoire impacting the deep subsurface element cycling and biomass production.
A marked difference in species diversity and OTU distribution as observed with respect to the sample’s depth may be linked primarily to the local rock geochemistry and the in situ temperature and pressure conditions of the horizons. The dominance of fewer archaeal species in the shallower depth BS than in GR possibly signifies the role of specialist organisms capable of thriving in relatively nutrient wealthy niches. On the other hand, a diverse but more even community structure as detected within the GR highlights the existence of more generalist archaeal groups which may work in synergy for the optimal use of limited available nutrients. A previous study has indicated that the diversity and resilience of archaeal communities in low nutrient deep biospheres could be promoted by mutualistic interactions among community members or emergent properties resulting from their slow reproduction rates (Nyyssönen et al., 2014). Spearman correlation (inter-taxa), PCA and abundance weighted NMDS illustrate the community assemblages and the role of subsurface geochemistry (including Fe2O3, organic carbon, nitrate and sulfate) in shaping such assemblages. The significance of local geochemical and geophysical factors on microbial community structure and function in various crystalline deep biosphere has previously been reported (Brazelton et al., 2013; Magnabosco et al., 2014; Ino et al., 2018).
Distribution patterns of archaeal taxa, correlation and network analyses clearly indicate the mutual relationship among the taxa highlighting the niche-specific colonization of specific populations within the BS and GR horizons. In general, a basalt centric archaeal community in younger basaltic rocks with relatively higher concentrations of organic carbon, iron oxide and sulfate could be delineated. This was distinct from the community assemblage within the organic carbon lean, deeper, Archean granitic rocks. At lower taxonomic level these basalt and granitic centric community composition illustrated the possible road maps of subsurface biogeochemistry, carbon flux and microbial (archaeal) ecology.
The basalt centric archaeal community of acidophilic, halophilic and methanogenic populations each with their varying abundances portrays the major archaeal guild and provides a possible clue to energy and carbon flows within the younger and relatively organic carbon rich basaltic rocks. Thermoacidophilic Thermoplasma, Ferroplasma, and Sulfolobaceae members abundant here are well known resident of Fe, S rich deep subsurface extreme environments with or without tectonic, volcanic and thermal activities (Golyshina, 2011, 2014; Reysenbach and Brileya, 2014; Heberle et al., 2015; Ijiri et al., 2018; Arce-Rodríguez et al., 2019). Thermoplasma and Ferroplasma both lack cell wall, and are heterotrophic, facultative aerobic, extreme (often obligate) acidophiles (pH 0.0–1.0) with optimum growth temperature at around 60°C. Thermoplasma gains energy anaerobically from sulfur respiration and produce H2S, while Ferroplasma has been recognized with its “ubiquitous capacity” for ferrous iron oxidation and a strict dependence on low concentration of yeast extract (Macalady et al., 2004; Golyshina, 2011, 2014; Reysenbach and Brileya, 2014). The simultaneous presence of both Thermoplasma and iron oxidizing Ferroplasma has been reported recently from acidic, iron-rich deep mine environment of Pyhäsalmi Mine, Finland (Bomberg et al., 2019). Significant connections between Ferroplasma and Fe2O3 in the major network hotspot corroborate their niche-specific localization. Together with Fe-oxidizing and S-respiring populations, the presence of Sulfolobaceae is also noticed in basalts. Sulfolobaceae members are known to adapt simultaneously with high temperature fluctuations (>60°C to <100°C) and low pH (<4.0), and they are equipped with diverse chemolithotrophic, heterotrophic, and mixotrophic metabolism (utilizing broad range of carbon sources including amino acids, simple to complex substrates) (Albers and Siebers, 2014). Based on the known physiological and metabolic properties of these archaeal taxa, we could infer that in iron-sulfur rich BS rocks of elevated temperature, thermoacidophiles allow extraction of energy through their chemolithotrophic and/or dissimilatory metabolism (coupled to heterotrophy) and thus generate biomass and metabolic intermediates for other members of this ecosystem. Together with these thermophilic taxa, presence of sparsely populated alkaliphilic Halobacteriaceae members detected exclusively in the BS samples presented an interesting aspect of this guild. Taxa affiliated to Halobacteriaceae are known mesophilic or moderate thermophilic halophile, able to utilize various organic substrates, can withstand slightly acidic to alkaline pH and reduce sulfur to sulfides (Robinson et al., 2005; Oren, 2014a; Thombre et al., 2016). Relative abundance of organic carbon in BS than in GR along with presence of methanogens capable of producing local alkalinity may be attributed to the presence of Halobacteriaceae in BS rocks. Co-occurrence network indicating the closeness of anaerobic hydrogenotrophic methanogens (Methanobacterium and Methanocalculus) with Halobacteriaceae members constituting the network hotspot possibly suggest a metabolic interrelation allowing co-occurrence of these two archaeal taxa. These hydrogenotrophic methanogens are known to have a lower threshold for hydrogen (geogenic and/or biogenic) concentration (compared to their acetoclastic counterparts) and harbor enzymes responsible for carbonate precipitation, thus allowing an increase in alkalinity of the local environment. Such a local rise in pH could possibly allow the alkaliphilic members to thrive (Smith and Ferry, 1999; Ussler and Paull, 2008). Absence of these methanogens along with lower levels of organic carbon in GR horizon might be one of the major factors limiting the presence of haloarchaeal members in deeper granitic horizons. Although most of these archaeal taxa have been reported previously in various deep biosphere investigations (Bomberg et al., 2015; Purkamo et al., 2015; Lau et al., 2016; Ijiri et al., 2018; Leandro et al., 2018; Hoshino and Inagaki, 2019), their abundance and association are important findings shading light on the potential roles of these extremophiles on deep biosphere biogeochemical carbon, sulfur cycles and energy flows (Purkamo et al., 2015; Jungbluth et al., 2016; Ijiri et al., 2018; Hoshino and Inagaki, 2019).
Archaeal microbiome within the organic carbon lean, basement Archaean granite rocks is distinct and constitute predominantly by acetoclastic methanogens (Methanomicrobia), ammonia oxidizing autotrophic archaea (SAGMCG-1, AK59, FHMa11 terrestrial group) and other taxa previously reported from diverse organic carbon deprived, deep crystalline environment with low oxygen (Takai et al., 2001; Kendall and Boone, 2006; Sørensen and Teske, 2006; Schneider et al., 2013; Oren, 2014b; Lauer et al., 2016; Wong et al., 2017; Purkamo et al., 2018; Smith et al., 2019). The presence of acetoclastic methanogenic Methanomicrobia as a dominant archaea in GR could be explained by relative higher concentration of acetate in these rocks, compared to BS (data not shown) as well as a lack of hydrogen generating minerals (e.g., olivine and pyroxene). The network consisting of these acetoclastic taxa (Methanosaeta, Methanolinea) shows their close proximity with anaerobic methane oxidizers (ANME-3) and major geological parameters, which had an increasing trend with depth. Though Methanosaeta and Methanolinea are both methanogenic in nature, the close phylogenetic relationship of Methanosaeta with anaerobic-methane oxidizers (ANME-1 and ANME-2) and the presence of the genes for the CO2 reduction pathway raises the possibility that Methanosaeta could have the capacity for the oxidation of methane (Smith and Ingram-Smith, 2007; Knittel and Boetius, 2009). It was initially thought that Methanosaeta was obligatory acetoclastic methanogen, but it was later reported that Methanosaeta can make direct electrical connections with Geobacter species, accepting electrons for the reduction of carbon dioxide to methane (Rotaru et al., 2014; Zhao et al., 2015). Archaeal genus ANME-3, which is also a part of this cluster, is not only capable of anaerobic oxidation of methane but also uses methane as one of the main carbon sources (Niemann et al., 2006). ANME-3 was earlier observed in close association with bacterial genus Desulfobulbus in the sulfide rich zones and there are evidences that confirm that Desulfobulbus can thrive on methane derived carbon (Niemann et al., 2006; Cui et al., 2015; Winkel et al., 2018). Interestingly, these three archaeal genera (Methanosaeta, Methanolinea, and ANME-3) also show negative correlation with TOC, which further supports the possibility of autotrophic lifestyle in the deeper rock cores (Knittel and Boetius, 2009). This sub-network consisting of both methanotrophs and methanogens not only indicates that methane might be an important driver of life in the deeper horizons of Deccan but also supports the possibilities of SLiME (Stevens and McKinley, 2000; Nealson et al., 2005) in the deep terrestrial subsurface of Deccan traps. Metagenome function is further predicted from 16S rRNA gene sequences using PICRUSt. PCA biplot constructed on the basis of the relative abundance of different predicted metabolisms corroborate the clustering discussed previously. Comparable trends obtained from different statistical analysis (on the basis of taxonomy and predicted community functions) determine that taxonomy of different archaeal groups in a particular metagenome reflects the archaeal community function.
In conclusion, this is the first detailed study on archaeal diversity and their niche-specific diversity shift in subterranean igneous provinces of the seismically active Koyna–Warna region of Deccan traps. Our results demonstrate a distinct partitioning of archaeal communities across the basaltic and granitic horizons corroborating the characteristic geochemical/physical conditions of these two horizons. Basaltic horizon harbors thermoacidophilic, chemolithotrophic, and organotrophic archaeal groups known to occupy environments with geothermal/volcanic/tectonic activities and with significant contributions in iron and sulfur cycles. A close association of halophilic taxa and hydrogenotrophic methanogenic archaea is noted. The deeper, more oligotrophic granitic horizons, in contrast, is dominated by chemoautotrophic archaea present in close association with acetoclastic methogenic taxa. The study provides a detailed report of archaeal communities from deep crystalline crustal environment highlighting the potential role of acidophilic, iron-oxidizing, sulfur-respiring, halophilic, methanogenic, methane/ammonia oxidizing archaea in carbon and energy flows in deep terrestrial crust.
AD carried out the microbiological analysis, statistical analysis, data organization, and manuscript preparation. SDG carried out the geochemical analysis and assisted in statistical analysis. AD and AG performed sub-coring of rocks, metagenomic DNA extraction, and bioinformatics pipeline optimization. AD, AG, and JS carried out qPCR analysis. HB assisted in data interpretation and manuscript preparation. SR organized sample collection and gave inputs about geophysical/geochemical characterization of Deccan lava flow and Koyna–Warna seismic zone. AM gave inputs about conducting different geochemical analysis. PS conceived the study, designed experiments, compiled and analyzed the data, wrote the manuscript, and done overall mentoring.
This study was supported by the Ministry of Earth Sciences (MoES), Government of India, Project ID: MoES/P.O.(Seismo)/1(181)/2013 dated April 7, 2014. AG is a recipient of research fellowship awarded by the Department of Biotechnology (DBT), Government of India (DBT/2014/IITKH/13). HB is a recipient of national postdoctoral fellowship awarded by the Department of Science and Technology (DST) (PDF/2017/000657). The scientific drilling investigations in the Koyna region was supported by the Secretary of Ministry of Earth Sciences, Government of India.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
AD, SDG, and JS gratefully acknowledge the fellowship provided by IIT Kharagpur. We thank the Director and all the participants/investigators from CSIR-National Geophysical Research Institute (NGRI), Hyderabad engaged in the exploratory scientific drilling operations at Ukhalu and Phansavale in the Koyna–Warna region of Deccan traps. We would like to acknowledge the Deep Carbon Observatory, Census of Deep Life Sequencing program for amplicon sequencing of all the extracted environmental DNAs (DCO_SAR). The authors thank Harsh Gupta (NGRI), Shailesh Nayak (MoES), and B. K. Bansal (MoES) for their unstinted motivation toward this work. We also thank Saibal Gupta of Department of Geology and Geophysics, IIT Kharagpur for his critical comments and suggestions on the manuscript. We gratefully acknowledge the computational support provided by the Department of Computer Science and Engineering, IIT Kharagpur.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.01362/full#supplementary-material
Arce-Rodríguez, A., Puente-Sánchez, F., Avendaño, R., Martínez-Cruz, M., de Moor, J. M., Pieper, D. H., et al. (2019). Thermoplasmatales and sulfur-oxidizing bacteria dominate the microbial community at the surface water of a CO 2-rich hydrothermal spring located in Tenorio Volcano National Park, Costa Rica. Extremophiles 23, 177–187. doi: 10.1007/s00792-018-01072-6
Arndt, D., Xia, J., Liu, Y., Zhou, Y., Guo, A. C., Cruz, J. A., et al. (2012). METAGENassist: a comprehensive web server for comparative metagenomics. Nucleic Acids Res. 40, W88–W95. doi: 10.1093/nar/gks497
Bhaskar Rao, Y. J., Sreenivas, B., Vijaya Kumar, T., Namrata, K., Kesava, A. K., and Babu, E. V. S. S. K. (2017). Evidence for Neoarchean Basement for the Deccan Volcanic Flows around Koyna-Warna Region, Western India: Zircon U-Pb Age and Hf-isotopic Results. J. Geol. Soc. India 90, 752–760. doi: 10.1007/s12594-017-0787-4
Biddle, J. F., Lipp, J. S., Lever, M. A., Lloyd, K. G., Sørensen, K. B., Anderson, R., et al. (2006). Heterotrophic Archaea dominate sedimentary subsurface ecosystems off Peru. Proc. Natl. Acad. Sci. U.S.A. 103, 3846–3851. doi: 10.1073/pnas.0600035103
Bomberg, M., Mäkinen, J., Salo, M., and Kinnunen, P. (2019). High diversity in iron cycling microbial communities in acidic, iron-rich water of the Pyhäsalmi Mine, Finland. Geofluids 2019, 1–17. doi: 10.1155/2019/7401304
Bomberg, M., Nyyssönen, M., Pitkänen, P., Lehtinen, A., and Itävaara, M. (2015). Active microbial communities inhabit sulphate-methane interphase in deep bedrock fracture fluids in Olkiluoto, Finland. Biomed Res. Int. 2015:979530. doi: 10.1155/2015/979530
Brazelton, W. J., Morrill, P. L., Szponar, N., and Schrenk, M. O. (2013). Bacterial communities associated with subsurface geochemical processes in continental serpentinite springs. Appl. Environ. Microbiol. 79, 3906–3916. doi: 10.1128/AEM.00330-13
Breuker, A., Köweker, G., Blazejak, A., and Schippers, A. (2011). The deep biosphere in terrestrial sediments in the Chesapeake Bay area, Virginia, USA. Front. Microbiol. 2:156. doi: 10.3389/fmicb.2011.00156
Brochier-Armanet, C., Boussau, B., Gribaldo, S., and Forterre, P. (2008). Mesophilic Crenarchaeota: proposal for a third archaeal phylum, the Thaumarchaeota. Nat. Rev. Microbiol. 6, 245–252. doi: 10.1038/nrmicro1852
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336.
Dutta, A., Dutta Gupta, S., Gupta, A., Sarkar, J., Roy, S., Mukherjee, A., et al. (2018). Exploration of deep terrestrial subsurface microbiome in Late Cretaceous Deccan traps and underlying Archean basement, India. Sci. Rep. 8:17459. doi: 10.1038/s41598-018-35940-0
Edwards, K. J., Becker, K., and Colwell, F. (2012). The deep, dark energy biosphere: intraterrestrial life on earth. Annu. Rev. Earth Planet. Sci. 40, 551–568. doi: 10.1146/annurev-earth-042711-105500
Gihring, T. M., Moser, D. P., Lin, L.-H., Davidson, M., Onstott, T. C., Morgan, L., et al. (2006). The distribution of microbial taxa in the subsurface water of the Kalahari Shield. S. Afr. Geomicrobiol. J. 23, 415–430. doi: 10.1080/01490450600875696
Golyshina, O. V. (2014). “The family ferroplasmaceae,” in Prokaryotes Other Major Lineages Bact. Archaea. eds E. Rosenberg, E. F. DeLong, S. Lory, E. Stackebrandt, and F. Thompson (Berlin: Springer), 29–34. doi: 10.1007/978-3-642-38954-2_325
Goswami, D., Akkiraju, V. V., Misra, S., Roy, S., Singh, S. K., Sinha, A., et al. (2017). Rock strength measurements on Archaean basement granitoids recovered from scientific drilling in the active Koyna seismogenic zone, western India. Tectonophysics 712–713, 182–192. doi: 10.1016/j.tecto.2017.05.029
Gregory, S. P., Barnett, M. J., Field, L. P., and Milodowski, A. E. (2019). Subsurface microbial hydrogen cycling: natural occurrence and implications for industry. Microorganisms 7:53. doi: 10.3390/microorganisms7020053
Gupta, H., Rao, N. P., Roy, S., Arora, K., Tiwari, V. M., Patro, P. K., et al. (2015). Investigations related to scientific deep drilling to study reservoir-triggered earthquakes at Koyna, India. Int. J. Earth Sci. 104, 1511–1522. doi: 10.1007/s00531-014-1128-0
Heberle, H., Meirelles, G. V., da Silva, F. R., Telles, G. P., and Minghim, R. (2015). InteractiVenn: a web-based tool for the analysis of sets through Venn diagrams. BMC Bioinform. 16:169. doi: 10.1186/s12859-015-0611-3
Hubalek, V., Wu, X., Eiler, A., Buck, M., Heim, C., Dopson, M., et al. (2016). Connectivity to the surface determines diversity patterns in subsurface aquifers of the Fennoscandian shield. ISME J. 10, 2447–2458. doi: 10.1038/ismej.2016.36
Ijiri, A., Inagaki, F., Kubo, Y., Adhikari, R. R., Hattori, S., Hoshino, T., et al. (2018). Deep-biosphere methane production stimulated by geofluids in the Nankai accretionary complex. Sci. Adv. 4:eaao4631. doi: 10.1126/sciadv.aao4631
Ino, K., Hernsdorf, A. W., Konno, U., Kouduka, M., Yanagawa, K., Kato, S., et al. (2018). Ecological and genomic profiling of anaerobic methane-oxidizing archaea in a deep granitic environment. ISME J. 12, 31–47. doi: 10.1038/ismej.2017.140
Ino, K., Konno, U., Kouduka, M., Hirota, A., Togo, Y. S., Fukuda, A., et al. (2016). Deep microbial life in high-quality granitic groundwater from geochemically and geographically distinct underground boreholes. Environ. Microbiol. Rep. 8, 285–294. doi: 10.1111/1758-2229.12379
Islam, E., Paul, D., and Sar, P. (2014). Microbial diversity in Uranium deposits from Jaduguda and Bagjata Uranium mines, India as revealed by clone library and denaturing gradient gel electrophoresis analyses. Geomicrobiol. J. 31, 862–874. doi: 10.1080/01490451.2014.907375
Jungbluth, S. P., Bowers, R. M., Lin, H.-T., Cowen, J. P., and Rappé, M. S. (2016). Novel microbial assemblages inhabiting crustal fluids within mid-ocean ridge flank subsurface basalt. ISME J. 10, 2033–2047. doi: 10.1038/ismej.2015.248
Kendall, M. M., and Boone, D. R. (2006). “The order methanosarcinales,” in The Prokaryotes, eds M. Dworkin, S. Falkow, E. Rosenberg, K. H. Schleifer, and E. Stackebrandt (New York, NY: Springer), 244–256. doi: 10.1007/0-387-30743-5_12
Langille, M. G. I., Zaneveld, J., Caporaso, J. G., McDonald, D., Knights, D., Reyes, J. A., et al. (2013). Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat. Biotechnol. 31, 814–821. doi: 10.1038/nbt.2676
Lau, M. C. Y., Kieft, T. L., Kuloyo, O., Linage-Alvarez, B., van Heerden, E., Lindsay, M. R., et al. (2016). An oligotrophic deep-subsurface community dependent on syntrophy is dominated by sulfur-driven autotrophic denitrifiers. Proc. Natl. Acad. Sci. U.S.A. 113, E7927–E7936. doi: 10.1073/pnas.1612244113
Lauer, A., Sørensen, K. B., and Teske, A. (2016). Phylogenetic characterization of marine Benthic Archaea in organic-poor sediments of the Eastern Equatorial Pacific Ocean (ODP Site 1225). Microorganisms 4:E32. doi: 10.3390/microorganisms4030032
Lazar, C. S., Stoll, W., Lehmann, R., Herrmann, M., Schwab, V. F., Akob, D. M., et al. (2017). Archaeal diversity and CO2 fixers in carbonate-/siliciclastic-rock groundwater ecosystems. Archaea 2017:2136287. doi: 10.1155/2017/2136287
Leandro, T., Rodriguez, N., Rojas, P., Sanz, J. L., da Costa, M. S., and Amils, R. (2018). Study of methanogenic enrichment cultures of rock cores from the deep subsurface of the Iberian Pyritic Belt. Heliyon 4:e00605. doi: 10.1016/j.heliyon.2018.e00605
Lever, M. A., Rouxel, O., Alt, J. C., Shimizu, N., Ono, S., Coggon, R. M., et al. (2013). Evidence for microbial carbon and sulfur cycling in deeply buried ridge flank basalt. Science 339, 1305–1308. doi: 10.1126/science.1229240
Macalady, J. L., Vestling, M. M., Baumler, D., Boekelheide, N., Kaspar, C. W., and Banfield, J. F. (2004). Tetraether-linked membrane monolayers in Ferroplasma spp: a key to survival in acid. Extremophiles 8, 411–419. doi: 10.1007/s00792-004-0404-5
Magnabosco, C., Ryan, K., Lau, M. C. Y., Kuloyo, O., Lollar, B. S., Kieft, T. L., et al. (2016). A metagenomic window into carbon metabolism at 3 km depth in Precambrian continental crust. ISME J. 10:730. doi: 10.1038/ismej.2015.150
Magnabosco, C., Tekere, M., Lau, M. C. Y., Linage, B., Kuloyo, O., Erasmus, M., et al. (2014). Comparisons of the composition and biogeographic distribution of the bacterial communities occupying South African thermal springs with those inhabiting deep subsurface fracture water. Front. Microbiol. 5:679. doi: 10.3389/fmicb.2014.00679
Miettinen, H., Kietäväinen, R., Sohlberg, E., Numminen, M., Ahonen, L., and Itävaara, M. (2015). Microbiome composition and geochemical characteristics of deep subsurface high-pressure environment, Pyhäsalmi mine Finland. Front. Microbiol. 6:1203. doi: 10.3389/fmicb.2015.01203
Misra, S., Roy, S., Bartakke, V., Athavale, G., and Gupta, H. (2017). Fissures and fractures in the Koyna seismogenic zone, western India. J. Geol. Soc. India 90, 131–137. doi: 10.1007/s12594-017-0690-z
Moser, D. P., Gihring, T. M., Brockman, F. J., Fredrickson, J. K., Balkwill, D. L., Dollhopf, M. E., et al. (2005). Desulfotomaculum and Methanobacterium spp. dominate a 4-to 5-kilometer-deep fault. Appl. Environ. Microbiol. 71, 8773–8783. doi: 10.1128/aem.71.12.8773-8783.2005
Mota, M. C., Carvalho, P., Ramalho, J., and Leite, E. (1991). Spectrophotometric analysis of sodium fluorescein aqueous solutions. Determination of molar absorption coefficient. Int. Ophthalmol. 15, 321–326. doi: 10.1007/bf00128951
Nealson, K. H., Inagaki, F., and Takai, K. (2005). Hydrogen-driven subsurface lithoautotrophic microbial ecosystems (SLiMEs): do they exist and why should we care? Trends Microbiol. 13, 405–410. doi: 10.1016/j.tim.2005.07.010
Niemann, H., Lösekann, T., De Beer, D., Elvert, M., Nadalig, T., Knittel, K., et al. (2006). Novel microbial communities of the Haakon Mosby mud volcano and their role as a methane sink. Nature 443, 854–858. doi: 10.1038/nature05227
Nyyssönen, M., Hultman, J., Ahonen, L., Kukkonen, I., Paulin, L., Laine, P., et al. (2014). Taxonomically and functionally diverse microbial communities in deep crystalline rocks of the Fennoscandian shield. ISME J. 8, 126–138. doi: 10.1038/ismej.2013.125
O’Connell, S. P., Lehman, R. M., Snoeyenbos-West, O., Winston, V. D., Cummings, D. E., Watwood, M. E., et al. (2003). Detection of Euryarchaeota and Crenarchaeota in an oxic basalt aquifer. FEMS Microbiol. Ecol. 44, 165–173. doi: 10.1016/S0168-6496(02)00465-8
Oren, A. (2014a). “The family Halobacteriaceae,” in The Prokaryotes: Other Major Lineages of Bacteria and The Archaea, eds E. Rosenberg, E. F. DeLong, S. Lory, E. Stackebrandt, and F. Thompson (Berlin: Springer), 41–121. doi: 10.1007/978-3-642-38954-2_313
Oren, A. (2014b). “The family Methanosarcinaceae,” in The Prokaryotes: Other Major Lineages of Bacteria and The Archaea, eds E. Rosenberg, E. F. DeLong, S. Lory, E. Stackebrandt, and F. Thompson (Berlin: Springer), 259–281. doi: 10.1007/978-3-642-38954-2_408
Purkamo, L., Bomberg, M., Kietäväinen, R., Salavirta, H., Nyyssönen, M., Nuppunen-Puputti, M., et al. (2016). Microbial co-occurrence patterns in deep Precambrian bedrock fracture fluids. Biogeosciences 13, 3091–3108. doi: 10.5194/bg-13-3091-2016
Purkamo, L., Bomberg, M., Nyyssönen, M., Kukkonen, I., Ahonen, L., and Itävaara, M. (2015). Heterotrophic communities supplied by ancient organic carbon predominate in deep fennoscandian bedrock fluids. Microb. Ecol. 69, 319–332. doi: 10.1007/s00248-014-0490-6
Purkamo, L., Kietäväinen, R., Miettinen, H., Sohlberg, E., Kukkonen, I., Itävaara, M., et al. (2018). Diversity and functionality of archaeal, bacterial and fungal communities in deep Archaean bedrock groundwater. FEMS Microbiol. Ecol. 94:fiy116. doi: 10.1093/femsec/fiy116
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2012). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219
Reysenbach, A.-L., and Brileya, K. (2014). “The family Thermoplasmataceae,” in The Prokaryotes, eds E. Rosenberg, E. F. DeLong, S. Lory, E. Stackebrandt, and F. Thompson (Berlin: Springer), 385–387. doi: 10.1007/978-3-642-38954-2_327
Robinson, J. L., Pyzyna, B., Atrasz, R. G., Henderson, C. A., Morrill, K. L., Burd, A. M., et al. (2005). Growth kinetics of extremely halophilic Archaea (family Halobacteriaceae) as revealed by Arrhenius plots. J. Bacteriol. 187, 923–929. doi: 10.1128/jb.187.3.923-929.2005
Rotaru, A.-E., Shrestha, P. M., Liu, F., Shrestha, M., Shrestha, D., Embree, M., et al. (2014). A new model for electron flow during anaerobic digestion: direct interspecies electron transfer to Methanosaeta for the reduction of carbon dioxide to methane. Energy Environ. Sci. 7, 408–415. doi: 10.1039/c3ee42189a
Roy, S., Rao, N. P., Akkiraju, V. V., Goswami, D., Sen, M., Bansal, B. K., et al. (2013). Granitic basement below Deccan traps unearthed by drilling in the Koyna Seismic zone, western India. J. Geol. Soc. India 81, 289–290. doi: 10.1007/s12594-013-0034-6
Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541. doi: 10.1128/AEM.01541-09
Schneider, D., Arp, G., Reimer, A., Reitner, J., and Daniel, R. (2013). Phylogenetic analysis of a microbialite-forming microbial mat from a hypersaline lake of the Kiritimati Atoll, Central Pacific. PLoS One 8:e66662. doi: 10.1371/journal.pone.0066662
Schoene, B., Samperton, K. M., Eddy, M. P., Keller, G., Adatte, T., Bowring, S. A., et al. (2015). U-Pb geochronology of the Deccan Traps and relation to the end-Cretaceous mass extinction. Science 347, 182–184. doi: 10.1126/science.aaa0118
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504. doi: 10.1101/gr.1239303
Silver, B. J., Raymond, R., Sigman, D. M., Prokopeko, M., Lollar, B. S., Lacrampe-Couloume, G., et al. (2012). The origin of NO3- and N2 in deep subsurface fracture water of South Africa. Chem. Geol. 294, 51–62. doi: 10.1016/j.chemgeo.2011.11.017
Smith, A. R., Kieft, B., Mueller, R., Fisk, M. R., Mason, O. U., Popa, R., et al. (2019). Carbon fixation and energy metabolisms of a subseafloor olivine biofilm. ISME J. 1, doi: 10.1038/s41396-019-0385-0
Stevens, T. O., and McKinley, J. P. (2000). Abiotic controls on H2 production from Basalt - Water reactions and implications for aquifer biogeochemistry. Environ. Sci. Technol. 34, 826–831. doi: 10.1021/es990583g
Stoddard, S. F., Smith, B. J., Hein, R., Roller, B. R. K., and Schmidt, T. M. (2015). rrn DB: improved tools for interpreting rRNA gene abundance in bacteria and archaea and a new foundation for future development. Nucleic Acids Res. 43, D593–D598. doi: 10.1093/nar/gku1201
Takai, K. E. N., Moser, D. P., DeFlaun, M., Onstott, T. C., and Fredrickson, J. K. (2001). Archaeal diversity in waters from deep South African gold mines. Appl. Environ. Microbiol. 67, 5750–5760. doi: 10.1128/aem.67.21.5750-5760.2001
Thombre, R. S., Shinde, V. D., Oke, R. S., Dhar, S. K., and Shouche, Y. S. (2016). Biology and survival of extremely halophilic archaeon Haloarcula marismortui RR12 isolated from Mumbai salterns, India in response to salinity stress. Sci. Rep. 6:25642. doi: 10.1038/srep25642
Ussler, W., and Paull, C. K. (2008). Rates of anaerobic oxidation of methane and authigenic carbonate mineralization in methane-rich deep-sea sediments inferred from models and geochemical profiles. Earth Planet. Sci. Lett. 266, 271–287. doi: 10.1016/j.epsl.2007.10.056
Waldron, P. J., Petsch, S. T., Martini, A. M., and Nüsslein, K. (2007). Salinity constraints on subsurface archaeal diversity and methanogenesis in sedimentary rock rich in organic matter. Appl. Environ. Microbiol. 73, 4171–4179. doi: 10.1128/aem.02810-06
Winkel, M., Mitzscherling, J., Overduin, P. P., Horn, F., Winterfeld, M., Rijkers, R., et al. (2018). Anaerobic methanotrophic communities thrive in deep submarine permafrost. Sci. Rep. 8:1291. doi: 10.1038/s41598-018-19505-9
Wong, H. L., Visscher, P. T., White, R. A. III, Smith, D.-L., Patterson, M. M., and Burns, B. P. (2017). Dynamics of archaea at fine spatial scales in Shark Bay mat microbiomes. Sci. Rep. 7:46160. doi: 10.1038/srep46160
Wu, X., Holmfeldt, K., Hubalek, V., Lundin, D., Åström, M., Bertilsson, S., et al. (2015). Microbial metagenomes from three aquifers in the Fennoscandian shield terrestrial deep biosphere reveal metabolic partitioning among populations. ISME J. 10:1192. doi: 10.1038/ismej.2015.185
Wuchter, C., Banning, E., Mincer, T., Drenzek, N. J., and Coolen, M. J. L. (2013). Microbial diversity and methanogenic activity of Antrim Shale formation waters from recently fractured wells. Front. Microbiol. 4:367. doi: 10.3389/fmicb.2013.00367
Keywords: deep biosphere, archaeal diversity, Deccan traps, Illumina sequencing, metabolism
Citation: Dutta A, Sar P, Sarkar J, Dutta Gupta S, Gupta A, Bose H, Mukherjee A and Roy S (2019) Archaeal Communities in Deep Terrestrial Subsurface Underneath the Deccan Traps, India. Front. Microbiol. 10:1362. doi: 10.3389/fmicb.2019.01362
Received: 30 July 2018; Accepted: 31 May 2019;
Published: 16 July 2019.
Edited by:Dawn Cardace, University of Rhode Island, United States
Reviewed by:Wei Xie, Sun Yat-sen University, China
Marco J. L. Coolen, Curtin University, Australia
Copyright © 2019 Dutta, Sar, Sarkar, Dutta Gupta, Gupta, Bose, Mukherjee and Roy. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.