Original Research ARTICLE
Comparative Metagenomics Provides Insight Into the Ecosystem Functioning of the Shark Bay Stromatolites, Western Australia
- 1Space Life Science Lab, Department of Microbiology and Cell Science, University of Florida, Gainesville, FL, United States
- 2Department of Microbiology and Cell Science, Genetics Institute, Institute for Food and Agricultural Sciences, University of Florida, Gainesville, FL, United States
- 3Genomics of Gene Expression Laboratory, Prince Felipe Research Center, Valencia, Spain
- 4EURA NOVA, Marseille, France
- 5Rosenstiel School of Marine and Atmospheric Science, University of Miami, Miami, FL, United States
Stromatolites are organosedimentary build-ups that have formed as a result of the sediment trapping, binding and precipitating activities of microbes. Today, extant systems provide an ideal platform for understanding the structure, composition, and interactions between stromatolite-forming microbial communities and their respective environments. In this study, we compared the metagenomes of three prevalent stromatolite-forming microbial mat types in the Spaven Province of Hamelin Pool, Shark Bay located in Western Australia. These stromatolite-forming mat types included an intertidal pustular mat as well as a smooth and colloform mat types located in the subtidal zone. Additionally, the metagenomes of an adjacent, non-lithifying mat located in the upper intertidal zone were also sequenced for comparative purposes. Taxonomic and functional gene analyses revealed distinctive differences between the lithifying and non-lithifying mat types, which strongly correlated with water depth. Three distinct populations emerged including the upper intertidal non-lithifying mats, the intertidal pustular mats associated with unlaminated carbonate build-ups, and the subtidal colloform and smooth mat types associated with laminated structures. Functional analysis of metagenomes revealed that amongst stromatolite-forming mats there was an enrichment of photosynthesis pathways in the pustular stromatolite-forming mats. In the colloform and smooth stromatolite-forming mats, however, there was an increase in the abundance of genes associated with those heterotrophic metabolisms typically associated with carbonate mineralization, such as sulfate reduction. The comparative metagenomic analyses suggest that stromatolites of Hamelin Pool may form by two distinctive processes that are highly dependent on water depth. These results provide key insight into the potential adaptive strategies and synergistic interactions between microbes and their environments that may lead to stromatolite formation and accretion.
Stromatolites are one of the most prevalent and recognizable components of the fossil record, dating back more than 3.7 Ga years (Grotzinger and Knoll, 1999; Nutman et al., 2016). These lithified, organosedimentary structures are a type of microbialite, formed by the sediment trapping, binding and/or carbonate precipitating activities of microorganisms in response to their local environment (Awramik et al., 1976; Burne and Moore, 1987). Stromatolites are distinguished from other microbialites by the presence of laminated microstructures, which are formed through iterative successions of sedimentation, microbial mat growth, and lithification (Reid et al., 2000; Dupraz and Visscher, 2005; Pace et al., 2018). These long-lived microbialite ecosystems have had a profound impact on the habitability of the planet, as they are attributed to changing the global redox conditions via oxygenic photosynthesis (Des Marais, 1991; Lyons et al., 2014). These properties make contemporary microbialites ideal model systems for investigating carbon cycling and the underlying processes associated with the precipitation and dissolution of calcium carbonate (Riding, 2000; Dupraz and Visscher, 2005; ?).
Modern microbialite ecosystems are widespread and are found in a diverse range of environments including, but not limited to, lacustrine ecosystems (e.g., Cuatro Ciénegas, Mexico; Lake Alchichica, Mexico; Pavilion and Kelly Lake, British Columbia) (Osborne et al., 1982; Ferris et al., 1997; Laval et al., 2000; Souza et al., 2006; Couradeau et al., 2011; Theisen et al., 2015; Chagas et al., 2016; White et al., 2016), freshwater systems (e.g., Deer Cave and Giblin River; Lundberg and McFarlane, 2011; Proemse et al., 2017); geothermal springs (e.g., Yellowstone National Park, WY, United States) (Inskeep et al., 2004; Pepe-Ranney et al., 2012), open marine environments (e.g., Exuma Sound, The Bahamas) (Dravis, 1983; Dill et al., 1986; Reid et al., 2000; Khodadad and Foster, 2012), and hypersaline waters (e.g., Kiritimati Atoll, Great Salt Lake) (Burns et al., 2004; Schneider et al., 2013; Ruvindy et al., 2016; Suosaari et al., 2016a; Lindsay et al., 2017).
Of the many extant examples of living, accreting microbialites, the largest, most extensive marine stromatolite-forming ecosystem is within the hypersaline embayment of Hamelin Pool, Shark Bay, Western Australia, a UNESCO world heritage site (Figure 1A; Playford, 1990; Reid et al., 2003; Suosaari et al., 2016a,b). Since their discovery in the 1950’s, the stromatolites of Hamelin Pool have served as important models to understand the formation of both modern and ancient stromatolite systems (Logan et al., 1974; Playford and Cockbain, 1976; Walter, 1976; Reid et al., 2003). Based on these pioneering studies, the stromatolite-forming microbial mats of Hamelin Pool are morphologically characterized into three canonical types: (1) ‘colloform’ mats, which form structures with moderately laminated carbonate fabrics; (2) ‘smooth’ mats, which are associated with highly laminated structures; and (3) ‘pustular’ stromatolite-forming mats, which exhibit no internal lamination (Logan et al., 1974; Playford and Cockbain, 1976). Historically, all carbonate build-ups within Hamelin pool have been referred to as ‘stromatolites’ regardless of the degree of lamination (Logan et al., 1974; Awramik et al., 1976; Playford, 1990; Jahnert and Collins, 2012; Suosaari et al., 2016a).
FIGURE 1. Overview of Hamelin Pool collection site and locations of stromatolite-forming and non-lithifying mat types within the pool. (A) Map of Provinces located within Hamelin Pool, Shark Bay. The stars reflect the collection sites of the colloform (yellow), smooth (red), pustular (pink), and non-lithifying mats (black) within the Spaven Province. The numbers 1, 2, 3 reflect transects associated with the collection of the mat types. (B) Cartoon depicting the position within the pool of the distinct microbial structures targeted in this study. Bar = 10 km.
Each of these distinctive stromatolite-forming mat types is found in different tidal zones within the pool (Playford and Cockbain, 1976; Jahnert and Collins, 2011; Suosaari et al., 2016b). The colloform stromatolite mats are the most seaward, located in the subtidal zone and can be found in depths up to four meters, whereas the pustular stromatolite-forming mats are the most shoreward in the intertidal zone (Figure 1B). In addition to these three prevalent stromatolite-forming mat communities, there are non-lithifying microbial mats in the upper intertidal zone including a smooth mat type and a pustular mat type (Wong et al., 2015; Suosaari et al., 2016b). These non-lithifying mats do not develop carbonate build-ups, but rather form extensive sheets, or fields of mats, along the coastline (Suosaari et al., 2016a).
The locations of the dominant stromatolite-forming mat types and morphological variations of stromatolites within Hamelin Pool have been extensively surveyed (Jahnert and Collins, 2013; Suosaari et al., 2016a,b). A recent mapping effort revealed eight geographical zones, or “Provinces,” each with distinctive and differentiated stromatolite macrostructures (Figure 1A; Suosaari et al., 2016a,b). One of the largest zones is Spaven Province, located along the western edge of the pool, spanning more than 20 km. The Spaven Province is distinguished by elongated and nested stromatolites, with the direction of the elongation being perpendicular from the shore (Suosaari et al., 2016b). All three of the canonical stromatolite-forming mat types, as well as the non-lithifying mats, are present within the Spaven Province; however, no in-depth metagenomic analyses have been completed for stromatolite-forming mats in this portion of Hamelin Pool.
Over the past decade, several studies have begun to characterize the microbial diversity associated with the stromatolite-forming and non-lithifying mats within the southern region of Hamelin Pool by targeting the SSU rRNA gene (Burns et al., 2004; Papineau et al., 2005; Goh et al., 2009; Wong et al., 2015, 2017; Suosaari et al., 2016a). These studies indicate that microbial populations within the different mat types are distinctive, with minimal heterogeneity (Wong et al., 2015; Suosaari et al., 2016a). More recently, metagenomic sequencing was used to compare the non-stromatolite forming smooth and pustular mats, with an intertidal, columnar stromatolite build-up from the Nilemah province (Ruvindy et al., 2016). These results indicated that the stromatolite-forming communities exhibit a metabolic potential that is distinct from that of non-lithifying smooth and pustular mats.
In this study, we build on this previous work by comparing the metagenomes of the three main lithifying stromatolite-forming mats (i.e., colloform, smooth, and pustular) to assess whether there are distinctive metagenomic signatures for each type. We also examine the non-lithifying pustular sheet-forming mat in the upper intertidal zone, as a control, to ascertain whether there are significant differences between the lithifying and non-lithifying mats within the Spaven Province. By employing a comparative metagenomic approach, the metabolic potential for each of the dominant mat types can be identified, thereby advancing our understanding of microbes and processes associated with stromatolite formation and accretion within Hamelin Pool.
Materials and Methods
Stromatolite-forming and non-lithifying mats were collected from the Spaven Province of Hamelin Pool, Shark Bay, Western Australia in April 2014. Although the various mat types are found throughout Hamelin Pool, relative abundance at specific locations is variable. Replicate stromatolite heads collected for this study spanned approximately 10 km within the Spaven Province (Table 1 and Figure 1A). Samples from each replicate head and non-lithifying sheet mats were collected using a sterile Harris 8.0 mm Uni-Core sampler (Ted Pella, Inc., Redding, CA, United States), immediately placed in RNAlater (Life Technologies, Inc., Grand Island, NY, United States), and stored at -20°C. Upon returning to the Space Life Science Lab, cores were stored at -20°C until DNA extraction.
DNA Extraction and Sequencing
Total genomic DNA was isolated (n = 8–10 extractions) from each replicate mat type samples using a xanthogenate bead beating method as previously described (Foster et al., 2009; Khodadad and Foster, 2012). Several modifications were made due to the high exopolysaccharide content of the mats. Each extraction replicate (60 mg) underwent three freeze-thaw cycles in liquid nitrogen, followed by immediate vortexing at max speed for 5 min with zirconia beads (Biospec Products, 2.0, 0.7, 0.1 mm diameter). All samples were visually inspected after the freeze-thaw incubations to ensure sufficient cell lysis. For those samples that did not exhibit full cell lysis due to extensive carbonate precipitate within the stromatolite-forming mat, samples were ground in liquid nitrogen prior to DNA extraction. Following extraction, DNA precipitation was conducted with 100% cold ethanol at -80°C, centrifuged, re-suspended in 70% cold ethanol, centrifuged and air dried for 2–5 min. DNA recovery was achieved with the DNeasy PowerSoil Kit (Qiagen, Germantown, MD, United States). Lastly, a second round of DNA precipitation was conducted with a final 0.3 M sodium acetate solution with cold 100% ethanol at -80°C to ensure removal of exopolysaccharides. DNA concentration was measured with Qubit® 2.0 fluorometer (Thermo Fisher Scientific, Waltham, MA, United States), normalized, and replicate extractions (n = 8 – 10) were pooled. DNA was sequenced at the University of Florida’s Interdisciplinary Center for Biotechnology Research using the NextSeq 500/550 High Output v2 Kit (paired-end, 150 cycles) on a NextSeq 500 sequencing system (Illumina, San Diego, CA, United States). All raw reads have been deposited in the NCBI Sequence Read Archive under BioProject number PRJNA429237.
Taxonomic Classification and Functional Annotation of Metagenomes
Raw sequences were quality filtered using sickle v1.33 with default parameters (e.g., minimum values of 20 bp length and Phred score of >20) for paired-end reads (Joshi and Fass, 2011). The range of alignment lengths for all 16 metagenomes ranged from 9 to 106 bp with a median of 48 bp. Additionally, the bit scores ranged from 20 to 125 with a median of 58 for all 16 metagenomes. For taxonomic classification, SSU rRNA sequences were recovered with SortMeRNA v2.1 using default parameters (Kopylova et al., 2012) against the SILVA_128 SSU Ref NR99 (Quast et al., 2013) database. Next, the SSU rRNA sequences were analyzed with QIIME version 1.9.1 (Caporaso et al., 2010) using the UCLUST method [open-reference operational taxonomic unit (OTU) picking] against the Silva1 and Greengenes v13.8 (DeSantis et al., 2006) databases. For further taxonomic analysis, the subsequent OTU table was inputted into MEGAN version 6.7.15, using the lowest common ancestor method with default parameters (Huson et al., 2016).
Quality controlled reads were annotated using BLASTx v2.2.26+ directly on the unassembled reads (Altschul et al., 1990) against the 2016 non-redundant UniProtKB/Swiss-Prot database (The UniProt Consortium, 2017) and filtered with an E-value cut-off of 10-3. An assembly free approach was used due to the variable taxonomic abundances within the stromatolite-forming mats, which has been shown to increase chimera formation frequency (Howe and Chain, 2015; Ghurye et al., 2016). Filtered BLAStx hits were programmatically linked to KEGG Orthology (KO) (Kanehisa et al., 2004) identifiers using the UniProt Retrieve/ID mapping tool. Unique KOs with their respective counts were then inputted into MEGAN to facilitate annotation of KEGG pathways (Huson et al., 2016). Pathway information for KOs not included in the MEGAN database, which was 62.3%, were identified via manual lookup on the KEGG website2.
Visualization and Statistical Analysis
To preserve the biological variability and variances within the metagenomic library sizes, the raw OTU table (absolute count) was normalized using the DESeq2 normalization technique (Love et al., 2014) with QIIME. Taxonomic diversity analysis was performed on each of the samples from the normalized OTU counts. For alpha diversity metrics, Shannon–Weaver (Shannon and Weaver, 1949) and Faith’s Phylogenetic Diversity (Faith, 1992) were used and significance of the diversity indices was performed with an adonis test (i.e., an analog of the nonparametric permutational manova) using 999 permutations. Beta-diversity metrics were generated from the unweighted UniFrac distance matrix and significance was determined using an adonis test. Phyla abundances were compared between any two mat types using a Wilcoxon test and statistical significance was et at p-value < 0.05 and visualized using the R package MetacodeR (Foster et al., 2017). Differential abundance analysis of the KEGG orthologs was conducted using the DESeq2 v1.16.1 statistical package in R (Love et al., 2014), which has been shown to be highly effective in comparing metagenomic data (McMurdie and Holmes, 2014; Rodriguez-R et al., 2015; Warden et al., 2016). As input, DESeq2 requires un-normalized count data for the statistical model to hold as it is designed to account for library size differences internally. Statistically significant differences between pairs of mat types were identified using a negative binomial Wald test; raw p-values were corrected for using the Benjamin–Hochberg adjustment and all adjusted p-values < 0.05 were considered to indicate differentially abundant KEGG orthologs between the mat types. Orthologs (i.e., KEGG level 4) were then mapped to their respective pathways (i.e., KEGG level 3) and results of the differential analysis were visualized with ggplot2 (Wickham, 2009). PCoA of the DESeq2 normalized gene counts was conducted using vegan v2.4-4 (Oksanen, 2017) with Bray–Curtis distances as input and visualized with ggplot2.
Site Description of Stromatolite-Forming Mats in Spaven Province
Within the Spaven Province, the three canonical, stromatolite-forming mat types were abundant; however, only the colloform and smooth mat types were associated with the elongated, nested stromatolite morphology characteristic of the Province (Figures 2A,B). Both the colloform and smooth stromatolites ranged from between 50 and 75 cm high and several meters in length. The pustular stromatolite-forming mats were associated with rounded, discrete heads that ranged from 30 to 50 cm in height and up to 1 m in width (Figure 2C). The non-lithifying mats formed sheets that extended for 100s of meters in the intertidal to supratidal zones (Figure 2D). The surfaces of each mat type showed distinctive features. The colloform stromatolites (Figure 2E) exhibited convoluted surfaces, whereas the smooth stromatolite-forming mats typically had flat, relatively unbroken surfaces (Figure 2F). Both the stromatolite-forming and non-lithifying pustular mats exhibited nodular surfaces, although the surface community on the stromatolite-forming pustular mats was darkly pigmented and often crusty compared to the non-lithifying mats (Figures 2G,H). In cross-section, all three of the stromatolite-forming mats exhibited extensive precipitation with a pronounced layer of cyanobacteria a few mm beneath the surface (Figures 2I–K), whereas the non-lithifying mat did not exhibit any lithification and often crumbled upon examination (Figure 2L).
FIGURE 2. Morphological features of mat types targeted in this study. (A–D) Underwater images depicting the macroscopic morphologies of the elongated nested build-ups associated with the colloform (A) and smooth (B) mat types, as well as the circular carbonate build-ups associated with the pustular (C) stromatolite-forming mats. The non-lithifying pustular mats form extensive sheets within the upper intertidal zone. (A,B) Bar = 50 cm; (C,D) Bar = 100 cm. (E–H) Surface images of the mat types showing the convoluted colloform (E), flat smooth (F), and nodular surfaces of the lithifying (G) and non-lithifying (H) pustular mats. Bar = 5 cm. (I–L) Cross-section of the mat surfaces depicting the domal structure of the colloform mat (I), flat surfaces of the smooth mats (J), heavily pigmented nodular structure of the lithifying pustular mat (K), and lighter pigmentation in the non-lithifying pustular mats (L). Bar = 0.5 cm. Arrows indicate cyanobacterial layer within the stromatolite forming mats.
Environmental data collected at the time of sampling as part of an independent study (Suosaari et al., 2016a,b) showed no significant temperature or tidal differences throughout the 20 km area of Spaven Province. The annual mean pH of the Spaven Province is 8.1 (Oehlert and Suosaari, personal communication) and the April water temperatures in Spaven Province at the collection sites ranged from 25.2 to 26°C at noon (the time of collection); salinity varied from 54.5 to 56 ppt.
Overview of Metagenomic Sequencing of Spaven Stromatolite-Forming Mats
Replicate metagenomic libraries were generated for each of the colloform, smooth, and pustular stromatolite-forming mats (Table 1). Additionally, libraries were also created for non-lithifying pustular mats to enable a direct comparison with the stromatolite-forming mats. The high-throughput sequencing effort produced a total of 16 metagenomes that contained 380,819,456 raw reads with an average of 23,801,216 per sample. For downstream analysis, an average of 2.4% of the sequences were removed to produce high-quality libraries (Phred quality score > 20; Table 1).
For taxonomic classification, initial assessment of the data revealed that the stromatolite-forming mat communities were more than 99% bacterial with archaea and eukaryotes comprising less than 1% of the total recovered reads. Few viral sequences were recovered and viruses were likely missed due to the DNA extraction and library preparation approaches used in this study. Bioinformatic mining of the SSU rRNA sequences from the metagenomes was conducted, producing an average of 65,588 SSU rRNA sequences per sample (Table 1). For functional annotation, all sequences were blasted against the UniProtKB/Swiss-Prot database, rendering the deduction of functionality from manually curated protein products and filtered with an e-value < 0.001, generating an average of 4,238,097 reads/sample (Table 1) that contained a KEGG Orthology (KO) identifier.
Comparison of Overall Microbial Diversity by Mat Type and Location
The taxa of the replicate metagenomes were compared using a principal coordinates analysis (PCoA) computed from the unweighted UniFrac distances between samples (Figure 3). When all OTUs (i.e., bacteria, archaea, eukaryotes) associated with the metagenomes were analyzed, three groupings were discernable (Figure 3A; adonis, p = 0.001; R2 = 0.27). First, the non-lithifying mats formed a cluster that was distinctive from the three stromatolite-forming mats. Second, the pustular stromatolite-forming mats formed a separate, although less cohesive, grouping, suggesting there is a higher level of heterogeneity within that mat type. Lastly, the colloform and smooth stromatolite-forming mats formed a single grouping, indicating a high degree of similarity between these two mat types. As sample collection spanned several km, the taxa of the metagenomes were also compared based on location within Spaven Province; Figure 3B). The taxonomic differences between the mat types was weakly linked to the collection area (adonis, p = 0.01; R2 = 0.16) and clustered primarily based on mat type.
FIGURE 3. Comparison of microbial diversity between lithifying stromatolite-forming mats and non-lithifying microbial mats in Spaven Province. Principal coordinates analysis plotting the unweighted Unifrac distances between samples according to mat type (A) and the three collection transects (Spaven 1, 2, 3) within Spaven Province (B).
Taxonomic Assessment of the Spaven Stromatolite-Forming Mats
Analysis of the 16 metagenomes revealed 33 bacterial, 3 archaeal and 3 eukaryotic phyla or superphyla. The overall alpha diversity within the different mat types was examined using both Shannon-Wiener Diversity and Faith’s Phylogenetic Diversity metrics (Supplementary Figures S1A,B). No statistical differences for either metric were observed (p > 0.05), indicating that the overall level of microbial diversity between the different mat types was comparable. The indices values were, however, much higher than previous reports (Allen et al., 2009; Goh et al., 2009; Garby et al., 2013; Wong et al., 2015; Suosaari et al., 2016b) likely reflecting the use of newer sequencing platforms (e.g., Illumina) that lead to an increase in the overall microbial diversity detected within all mat types compared to previous technologies.
In all the different mat types, the dominant bacterial phyla included the Proteobacteria (55–69%) and Cyanobacteria (15–29%), with the Planctomycetes (5–7%), Bacteroidetes (3–7%), Verrucomicrobia (2–3%), Chloroflexi (2–5%), and Actinobacteria (0.25–2%) phyla present in all of the mat types (Figure 4A and Supplementary Figure S1C). At much lower abundances (relative abundance < 1% across all replicates), a total of 26 additional bacterial phyla were observed (Supplementary Figure S1D). Although bacteria dominated all mat types, the recovered archaea were predominantly Thaumarchaeota (Supplementary Figure S1E) whereas the eukaryotic taxa were primarily associated with the Opisthokonta and the SAR supergroup, which includes the stramenopiles, alveolates and Rhizarias (Supplementary Figure S1F). These results are consistent with previously documented taxa within the southern region of Hamelin Pool (Burns et al., 2004; Goh et al., 2009; Edgcomb et al., 2014; Wong et al., 2015; Ruvindy et al., 2016; Suosaari et al., 2016a).
FIGURE 4. Taxonomic distribution of bacteria within the targeted mat types. (A) Overview of the top 35 bacterial taxa recovered from the 16 metagenomic libraries generated in this study and serves as orientation for (B). (B) Differential heat trees of pairwise comparisons between different mat types highlighting significant taxonomic enrichments (Wilcoxon test, p < 0.05). The color of each taxon represents the log-2 ratio of the median of the proportions observed with each mat type and corresponds to the different pairwise comparisons listed along the columns or rows. Node/edge color and size display the relative proportion for each taxon with lines in red reflecting those taxa enriched in mats along the x-axis and lines in blue enriched along the y-axis.
Although few differences were observed between mat types at the phylum-level, in-depth pairwise comparisons of the bacteria revealed significant differences at lower taxonomic levels (Wilcoxon test, p < 0.05) and are visualized as differential abundance heat trees in Figure 4B. The most pronounced differences occurred between the stromatolite-forming mat types and the non-lithifying mats. There was a pronounced overall enrichment of Proteobacteria, in particular, Alpha-, Delta- and Gammaproteobacteria taxa in the non-lithifying mats compared to the stromatolite-forming mat types (Figure 4B and Supplementary Figure S2); however, some individual proteobacterial taxa were enriched in the stromatolite-forming mats. For example, in the smooth mats, there was a statistically significant enrichment of Desulfovibrionales compared to the non-lithifying mats, whereas in the colloform mats increases in several Rhodobacterales taxa were observed. In the colloform and smooth stromatolite-forming mats, there was also an enrichment of Bacteroidetes, specifically in the Flavobacteriales. The colloform mats also were enriched in Saprospirales and Cytophagales compared to the non-lithifying mats.
Differences were also observed between the three dominant stromatolite-forming mat types. The most pronounced difference was an enrichment of the coccoid cyanobacteria Synechococcales within the pustular mats compared to the colloform and smooth mat types (Figure 4B and Supplementary Figure S3; Wilcoxon test, p < 0.05). There was also a differential abundance of two different Chroococcales lineages. The colloform and smooth mats were enriched in the genus Chroococcus, whereas the pustular mats were enriched in the Gomphosphaeriaceae family. Both the smooth and colloform mats showed an enrichment of Verrucomicrobia and Bacteroidetes compared to the pustular stromatolite-forming mats. There were very few differences observed between the colloform and smooth stromatolite-forming mats with an increased abundance of the Cytophagales in the colloform mats and a few enriched Alphaproteobacteria (Rhodobacterales and Rhizobiales) in the smooth stromatolite-forming mats (Figure 4B).
Comparison of the Functional Genes Within the Stromatolite-Forming Mats
To complement the taxonomic comparison, an ortholog-based comparison between the different mat types was performed using PCoA analysis of KEGG orthologies (Figure 5). When all annotated genes were considered, there were primarily two distinct clusters, one that included the colloform and smooth stromatolite-forming mats and a second that contained the pustular stromatolite-forming mats and non-lithifying mats (Figure 5A; adonis, p = 0.001; R2 = 0.40).
FIGURE 5. Comparison of the functional gene profiles between the stromatolite-forming and non-lithifying mat types. Principal coordinate analysis (PCoA) of all DESeq2 normalized KEGG level 4 gene counts of (A) all recovered genes and (B) differentially abundant genes (p-value < 0.001) from the different mat types.
In an effort to further characterize the metabolic differences between the mat types, pairwise comparisons of the differentially abundant genes were conducted with DESeq2 (Table 2; Love et al., 2014). There were 1232 differentially abundant genes (p < 0.05) between the non-lithifying mats and the colloform stromatolite-forming mats and 1033 between the non-lithifying mats and smooth stromatolite-forming mats, respectively (Table 2). A PCoA visualization was performed with statistically significant genes (p < 0.001), as was determined by differential abundance analysis with DESeq2 (Figure 5B). This approach more clearly revealed differences between the non-lithifying and pustular mat types, however, the colloform and smooth mats were still tightly clustered.
TABLE 2. Differential gene abundance comparison between lithifying and non-lithifying mat types abundance mat types.
Interestingly, there were only four significant differentially abundant genes (p < 0.05) between the colloform and the smooth stromatolite-forming mat types, suggesting the genetic profiles of the microbes in these communities are highly similar. For those pairwise comparisons that had >1000 differentially abundant genes, filtering criteria were employed to identify those genes with a p-value < 0.001 and an absolute log2-fold change > 1 (Table 2 and Supplementary Table S1).
In pairwise comparisons between the shoreward, non-lithifying mats and the stromatolite-forming mats, there were enrichments in numerous glycan and polysaccharide biosynthesis and degradation genes in the colloform and smooth stromatolite-forming mats, such as beta-porphyranases, fuctosidases, fructosideases, and rhamnotransferases (Supplementary Figures S4–S7). These genes were associated with a wide distribution of taxa including Chroococcales, Planctomycetales, several orders of Bacteroidetes, Firmicutes, and Proteobacteria (Supplementary Figures S8–S12). There was also a prevalence of genes associated with dissimilatory sulfur metabolism in the colloform and smooth stromatolite-forming mats compared to the non-lithifying mats, such as genes encoding thiosulfate reductases, dimethyl sulfoxide reductases (Supplementary Figures S4–S7 and Supplementary Table S1). There were fewer significant differences between the non-lithifying and pustular stromatolite-forming mats with the highest differentially abundant gene with regard to fold change associated with pigment transport in the pustular stromatolite and several uncharacterized proteins (Supplementary Figures S4–S7 and Supplementary Table S1).
In the upper intertidal, non-lithifying mats there was a higher representation of genes associated with oxidative and osmotic stress responses, including nitric oxide reductases, peroxidases, glucosylglycerol-phosphate synthase (GgpP), and glycine betaine/proline transporters (Supplementary Figures S4–S7). These genes were associated with a wide range of Proteobacteria, including several methylotrophic bacteria and purple-sulfur bacteria (Supplementary Figures S8–S12). There was also a significant enrichment of genes associated with metalloid and heavy metal cycling in the non-lithifying mats, such as arsenite transporters (e.g., arsAB), arsenite-mycothiol transferase (arsC) for removal, arsenic resistance genes (e.g., arsH), and arsenite oxidases (e.g., aoxAB) as well as genes that are part of the cobalt-zinc-cadmium efflux system (Supplementary Figures S4–S7). Genes associated with arsenate reductase were recovered from the non-lithifying mats, but they were not differentially abundant compared to the stromatolite-forming mat types. The genes associated with arsenic metabolism were widespread in the Actinobacteria, Cyanobacteria, and Firmicutes phyla, whereas the genes associated with the cobalt-zinc-cadmium system were found only in the Firmicutes and cyanobacterial Synechococcales (Supplementary Figures S8–S12).
In pairwise comparisons between the different stromatolite-forming mats, the most significant differentially abundant genes were observed in the pustular stromatolite-forming mats (Figure 6 and Supplementary Figures S4–S7). There was an increase in the abundance of genes associated with carbon-concentrating mechanisms and photosynthesis compared to the colloform and smooth stromatolite-forming mats (p < 0.001; Figure 6 and Supplementary Figures S4–S7). These differentially abundant photosynthesis genes were associated with both photosystem I and II pathways, as well as numerous antenna proteins that included phycocyanins and the blue phycocyanobilin (Supplementary Table S1). The differentially abundant photosynthesis genes were widely distributed in the coccoid Chroococcales and Synechococcales orders as well as the filamentous Nostocales cyanobacteria (Figure 7 and Supplementary Figures S8–S12). As with the non-lithifying mats, the colloform and smooth stromatolite-forming mats were enriched in genes associated with heterotrophic metabolisms, such as methane and sulfur reduction metabolisms. Differentially abundant genes associated with these pathways included co-methyltransferase, heterodisulfide reductase, thioredoxin reductase, adenylylsulfate (APS) reductases and sulfide dehydrogenases (Supplementary Table S1). The methanogenesis genes were associated with Euryarchaeota (e.g., Methanobacteriales, Methanococcales, Methanopyrales, and Methanosarcinales), whereas the dissimilatory sulfur reduction genes were derived from the Desulfovibrionales (Supplementary Figures S8–S12).
FIGURE 6. Pairwise comparison of differentially abundant genes between pustular and smooth stromatolite-forming mats. Representative differentially abundant genes (p < 0.001) colored by their respective pathways. A positive log2-fold change designates differential abundance in the smooth mat type and a negative log2-fold change designates differential abundance within the pustular mats. Lines indicate the standard error of the mean.
FIGURE 7. Taxa associated with selected functional pathways. Heatmap showing the relative contribution of each taxon with respect to the differentially abundant genes within the pustular and the smooth stromatolite-forming mats.
Hamelin Pool represents the largest known marine habitat for actively accreting stromatolites; however, little is known about the molecular pathways underlying the different stromatolite-forming mat types in this ecosystem. In this study, the metagenomes of the three canonical microbial mat types associated with carbonate build-ups were compared to each other and to an adjacent sheet-forming mat to assess differences in the metabolic signatures of lithifying and non-lithifying microbial mats. The results of this study provide evidence that: (1) the stromatolite-forming mats types have distinct taxonomic and functional gene profiles compared to non-lithifying mats; (2) colloform and smooth stromatolite-forming mats exhibit distinctive surface morphologies, yet show few taxonomic and functional gene differences; and (3) intertidal pustular stromatolite-forming mats are enriched in photosynthetic genes compared to the colloform and smooth stromatolite-forming mats, suggesting two distinctive metabolic processes driving stromatolite formation within Hamelin Pool.
The metagenomic comparisons in this study provided an important, in-depth analysis of those differentially abundant metabolic pathways between lithifying and non-lithifying mat communities. Our current study expands upon previous metagenomic analyses (Ruvindy et al., 2016) to examine all three of the prevalent lithifying mat types. Our results revealed that one of the pronounced differences within the non-lithifying mats compared to the subtidal stromatolite-forming mat types was the enrichment of genes associated with metalloid cycling (p < 0.001), in particular, genes associated with arsenic metabolism. The presence of arsenic metabolism has been previously reported in the non-lithifying mats of the southern Nilemah Province within Hamelin Pool (Ruvindy et al., 2016), hypersaline hot spring biofilms from Mono Lake (Kulp et al., 2008), as well as within the high altitude stromatolites of Socompa Lake, Argentina (Kurth et al., 2017). The increase in abundance of these genes associated with metal cycling within the mats of the upper zone may reflect the hypersaline nature of these habitats (Hamelin Pool. 66 – 88 ppt; Mono Lake 84 ppt; Socampa Lake ∼90 ppt; Farías et al., 2013). For example, evaporation can increase the abundance and availability of metalloids, such as arsenic and boron, within microbial communities (Kulp et al., 2007). Additionally, under extreme saline conditions some metabolisms, such as sulfate reduction and methanogenesis, are diminished (Oren, 1999; Kulp et al., 2007; Wong et al., 2017) and reliance on arsenic cycling, particularly under anoxic conditions, may reflect an important metabolic strategy for these mat communities exposed to the energetically taxing high salt and desiccation conditions of the upper, intertidal zone.
Another pronounced difference within the non-lithifying mats compared to the stromatolite-forming mat types in Spaven Province was the enrichment of oxidative (e.g., nitric oxide reductases, peroxidases) and osmotic stress (e.g., glycine betaine/proline synthesis and transport) genes in the upper tidal non-lithifying mats, whereas genes typically associated with UV stress were found to be evenly represented in all of the different mat types examined. The enriched oxidative and osmotic stress genes in the non-lithifying mats corresponded to several of observed osmoadaptive genes and pathways recovered from hypersaline mat environments, including Hamelin Pool (Goh et al., 2010; Goh et al., 2011; Gudhka et al., 2015; Gunde-Cimerman et al., 2018). As there is a pronounced desiccation gradient within the pool (Burne and Johnson, 2012; Suosaari et al., 2016b), the enrichment of genes associated with the production of osmoregulation and oxidative stress responses in the non-lithifying mats likely reflects the dynamic tidal extremes within the pool. Additionally, several taxonomic analyses of the non-lithifying mats in the Nilemah Province of Hamelin Pool have shown distinctive gradations of key functional groups of microbes, such as sulfate-reducing bacteria and methanogens, at different tidal regimes further suggesting that water levels are influencing the mat communities (Wong et al., 2017). The enriched genes associated with pigment transport were associated with cyanobacteria and may reflect the higher light and salt selection pressures on the cyanobacteria within this mat type, as several cyanobacteria have shown differential pigment production and transport responses to light and salt stress (Riethman et al., 1988; Sudhir and Murthy, 2004). Several comparative genomic studies of marine cyanobacteria have revealed genome expansion as the result of nutrient or light stresses (Swingley et al., 2008) and light pressures can be a major driver of genome differences in pigments and the photosynthesis apparatus between closely related cyanobacterial ecotypes (Rocap et al., 2003).
Over the course of a year, the water levels within the pool can vary by as much as two meters and are highly dependent on seasonal, meteorological and astronomical effects (Burne and Johnson, 2012; Playford et al., 2013; Suosaari et al., 2016b). Although water levels are typically higher during the austral summer and fall, which coincided with the April collection of the samples used in this study, the location of Spaven Province on the eastern margin of the pool experiences the largest meteorological effect of tidal variation throughout the year, spanning 1.67 m (Suosaari et al., 2016a), indicating the non-lithifying mats would experience frequent exposure and desiccation throughout the diel cycle. Temperatures within the pool can also fluctuate between 11 and 33°C through the year (Suosaari et al., 2016a); however, due to the shallowness of the pool, there is no vertical temperature stratification within the water column (Burling et al., 2003). Therefore, it is more likely that water depth, has a stronger influence than temperature on driving the differences in the metagenomes between the mat types.
The concept that water depth influences the morphology and community composition of the stromatolite-forming mats has long been suggested (e.g., Logan, 1961; Logan et al., 1974; Golubic and Hofmann, 1976; Playford, 1990; Jahnert and Collins, 2011, 2012; Playford et al., 2013; Suosaari et al., 2016a; Wong et al., 2017). More recent studies suggest that water level and the underlying shelf physiography may be two of the most important drivers of microbial mat zonation within Hamelin Pool (Suosaari et al., 2016b). The PCoA comparisons of both the taxonomic and functional gene diversity (Figures 3, 5) of the different mat types provide support for these concepts, as the stromatolite-forming mats formed three distinct clusters that corresponded with the upper tidal, intertidal and subtidal zones.
Interestingly, within the subtidal zone, the colloform and the smooth stromatolite-forming mats showed few taxonomic and functional gene differences. In fact, only four of the 15,378 recovered genes were significantly different between colloform and smooth stromatolites (Table 2 and Supplementary Figure S5) indicating that these two subtidal mat types share almost identical metabolic potential. Despite the similar metabolic profiles between the two subtidal stromatolite-forming mat types, each of the associated carbonate build-ups exhibited different surface morphologies, with the colloform mats exhibiting a convoluted surface and the smooth mats having a flat, even surface. These differences could reflect differential gene expression between the two mats, as other metatranscriptomic analyses of Bahamian thrombolites have shown that the metatranscriptome can vary significantly throughout diel and seasonal cycles and does not fully reflect the metagenome (Mobberley et al., 2015; Louyakis et al., 2018).
Alternatively, the differences in morphology may reflect environmental factors impacting the mat surfaces, specifically, the shear forces of winds and currents. Although water circulation within Hamelin Pool is restricted due to the presence of the Faure Sill in the northern part of the pool, it is primarily driven by the wind resulting in Langmuir circulation, which are shallow, counter-rotating gyres that align with the winds (Playford et al., 2013; Suosaari et al., 2016b). In the Spaven Province, current velocities range from as low as 0.001 m/s to a maximum of 0.54 m/s with a mean of 0.125 m/s (Suosaari et al., 2016b). Waves driven by the wind could be differentially impacting the colloform and smooth stromatolite-forming mat communities in the Spaven Province resulting in the disparate morphologies. Additional analyses, such as metatranscriptomics and proteomics, will be required to ascertain whether the colloform and smooth microbial communities exhibit differential microbial activities that account for the differences in morphology or whether these differences are the product of extrinsic factors.
One key difference amongst the stromatolite-forming mat communities was the enrichment of genes associated with photosynthesis in the pustular mats, including genes associated with photosystems I/II, as well as a range of pigments and antennae proteins. These differentially abundant genes were widely distributed within the cyanobacteria and were associated with the coccoid Chroococcales and Synechococcales, as well as the filamentous Nostocales orders. Due to the paucity of sequenced genomes for cyanobacteria, many of the genes were unable to be identified beyond the order level. However, over the past few decades, microscopic analyses of the pustular mats have shown that Entophysalis spp., which belongs to the order Chroococcales as well as the filamentous Nostocales Scytonema spp. and Dichothrix spp., dominate the pustular mat types (Logan et al., 1974; Golubic and Hofmann, 1976; Collins and Jahnert, 2014; Suosaari et al., 2016a). These taxa have long been thought to play a key role in the formation of stromatolites, both in living and in ancient systems (Golubic and Hofmann, 1976; Golubic and Campbell, 1981). Additionally, recent metagenomic and metatranscriptomic examination of the unlaminated thrombolites of Highborne Cay, The Bahamas have shown that the filamentous Dichothrix spp. and its associated coccoid cyanobacteria are the most transcriptionally active taxa within the thrombolite communities, with most transcripts associated with photosynthesis (Mobberley et al., 2015; Louyakis et al., 2017, 2018). These results coupled with stable isotope analyses of the calcium carbonate have revealed that the primary metabolism driving precipitation in the intertidal Bahamian thrombolites is photosynthesis (Louyakis et al., 2017).
The enrichment of photosynthetic genes in the pustular stromatolite-forming mats of Hamelin Pool suggests that, like the Bahamian thrombolites, photosynthesis may be the primary metabolism driving precipitation in this mat type. Additionally, as the carbonate build-ups associated within the pustular mats are typically unlaminated, and have been called stromatolites rather than thrombolites for historical reasons (Logan, 1961; Logan et al., 1974; Golubic and Hofmann, 1976; Playford, 1990; Jahnert and Collins, 2012, 2013; Playford et al., 2013), the processes associated with pustular carbonate build-ups may be distinct from the other lithifying mat structures and may be more thrombolite-like. Additional research into the biogeochemical and transcriptional activities of the filamentous Dichothrix spp. and coccoid Entophysalis spp. is needed to fully elucidate the role of these abundant taxa in the formation of the pustular-associated carbonate structures.
Whereas the pustular stromatolite-forming mats harbored an increase in photosynthesis genes, the colloform and smooth stromatolite-forming mats exhibited a differential abundance of genes associated with certain heterotrophic pathways known to promote carbonate precipitation. For example, there was a differential abundance of genes associated with dissimilatory sulfate reduction and methanogenesis in the colloform and smooth mats. Previous research on the stromatolites of Highborne Cay, The Bahamas has shown a strong correlation between sulfate-reducing activity and lithified micritic lamination (Visscher et al., 1998, 2000). Bicarbonate is a product of sulfate reduction, which can alter the local pH within the mat microenvironment, thereby promoting carbonate precipitation (Dupraz and Visscher, 2005; Visscher and Stolz, 2005; ?). Methanogenesis can also increase the local pH subsequently facilitating conditions for carbonate precipitation (Visscher and Stolz, 2005) and several genes associated with that pathway were enriched in the lithifying colloform and smooth mats compared to the pustular stromatolite-forming mats.
In addition to creating localized geochemical environments that favor carbonate precipitation (?), nucleation sites must also be available. The primary nucleation sites of the lithifying mats are within the exopolymeric substances (EPS) secreted outside the cells forming a matrix (Braissant et al., 2007; ?). The EPS matrix provides a multi-faceted role for the mat community, not only acting as nucleation sites for carbonate precipitation but also serves as a protective barrier against environmental stressors and enhances community stability in high-energy environments (Decho et al., 2005; Decho and Gutierrez, 2017). Within the EPS, negatively charged acidic groups can bind Ca2+, thereby sequestering it within the matrix (Kawaguchi and Decho, 2002) and through modification to the EPS, either by environmental or metabolic degradation, the ions can be released, thus increasing local calcium concentration available for mineralization (Visscher et al., 2000; Dupraz and Visscher, 2005; ?).
In the colloform and smooth stromatolite-forming mats, there was an increase in the abundance of genes associated with formation and degradation of EPS. Specifically, there was an increase in the abundance of metabolic enzymes associated with modification of dicarboxylic acids (e.g., maleic acid) and deoxyhexoses (e.g., fucose and rhamnose) from a wide range of taxa including Bacteroidetes, Firmicutes, Planctomycetes, and Proteobacteria compared to the pustular stromatolite-forming mats. This enrichment suggests that the EPS material of the colloform and smooth stromatolite-forming mats is highly labile and may fuel heterotrophic activity within the mats. The processing and restructuring of the EPS by heterotrophic activity may change the physiochemical properties of the EPS matrix, thereby influencing the morphology and mineralogy of the precipitates within colloform and smooth stromatolites. Together, these results suggest the processes facilitating carbonate precipitation in the laminated colloform and smooth stromatolite-forming mats may rely more on heterotrophic processes than in the unlaminated pustular stromatolite-forming mat type.
In summary, our results provide the first metagenomic comparison of the three dominant stromatolite-forming mat types within Hamelin Pool. Our results reveal key metabolic differences amongst the different stromatolite-forming mat communities, in particular between the intertidal and subtidal zones. The results suggest that carbonate precipitation in the unlaminated, pustular mat build-ups may be primarily photosynthesis driven, whereas, in the laminated colloform and smooth stromatolite structures, carbonate precipitation may be the product of a synergism between autotrophic and heterotrophic processes in the associated mats. More in-depth analysis of the microbial activities within each of these lithifying mat types will be required to more fully understand the metabolic drivers of carbonate precipitation within Hamelin Pool. Taken together, this comparative metagenomic analysis has provided confirmation of many of the prior geological and morphological-based studies in Hamelin Pool, thereby providing important insight into the feedbacks between microbial mat communities and their environments, which together drive stromatolite formation.
JF and PR conducted all the fieldwork. JB conducted the experiments, sequencing analysis, and wrote the first draft of the manuscript. AC, GC, CP, AL, and JF assisted in developing the bioinformatic pipeline and analysis of the data. All authors assisted in the writing and editing of the manuscript.
This work was funded by a NASA Exobiology and Evolutionary Biology program (NNX14AK14G) awarded to JF and PR. JB was supported by the NSF Graduate Research Fellowship Program and a NASA Earth and Space Science Fellowship.
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.
The authors would also like to thank the late Dr. Phillip Playford for discussions regarding this project as well as his discovery and pioneering work on the stromatolite systems of Shark Bay. We thank Erica Suosaari for her geological expertise and for help in navigating through Hamelin Pool. We thank Amanda Oehlert, Jessica Babilonia, Lexi Duscher, and Michael Morrison for their technical assistance. We also thank Bush Heritage Australia for logistical assistance and to the Western Australia Department of Parks and Wildlife (formerly Department of Environment and Conservation) and the federal Department of Sustainability, Environment, Population, and Communities for sampling permits. This manuscript is Hamelin Stromatolite Contribution Series #5.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.01359/full#supplementary-material
Allen, M. A., Goh, F., Burns, B. P., and Neilan, B. A. (2009). Bacterial, archaeal and eukaryotic diversity of smooth and pustular microbial mat communities in the hypersaline lagoon of Shark Bay. Geobiology 7, 82–96. doi: 10.1111/j.1472-4669.2008.00187.x
Braissant, O., Decho, A. W., Dupraz, C., Glunk, C., Przekop, K. M., and Visscher, P. T. (2007). Exopolymeric substances of sulfate-reducing bacteria: interactions with calcium at alkaline pH and implication for formation of carbonate minerals. Geobiology 5, 401–411. doi: 10.1111/j.1472-4669.2007.00117.x
Burns, B. P., Goh, F., Allen, M., and Neilan, B. A. (2004). Microbial diversity of extant stromatolites in the hypersaline marine environment of Shark Bay, Australia. Environ. Microbiol. 6, 1096–1101. doi: 10.1111/j.1462-2920.2004.00651.x
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. doi: 10.1038/nmeth.f.303
Chagas, A. A. P., Webb, G. E., Burne, R. A., and Southam, G. (2016). Modern lacustrine microbialites; towards a synthesis of aqueous and carbonate geochemistry and mineralogy. Earth Sci. Rev. 162, 338–363. doi: 10.1016/j.earscirev.2016.09.012
Couradeau, E., Benzerara, K., Moreira, D., Gerard, E., Kazmierczak, J., Tavera, R., et al. (2011). Prokaryotic and eukaryotic community structure in field and cultured microbialites from the alkaline Lake Alchichica (Mexico). PLoS One 6:e28767. doi: 10.1371/journal.pone.0028767
Decho, A. W., Visscher, P. T., and Reid, R. P. (2005). Production and cycling of natural microbial exopolymers (EPS) within a marine stromatolites. Palaeogeogr. Palaeoclimatol. Palaeoecol. 219, 71–86. doi: 10.1016/j.palaeo.2004.10.015
DeSantis, T. Z., Hugenholtz, P., Larsen, N., Rojas, M., Brodie, E. L., Keller, K., et al. (2006). Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl. Environ. Microbiol. 72, 5069–5072. doi: 10.1128/AEM.03006-05
Dupraz, C., Reid, R. P., Braissant, O., Decho, A. W., Norman, R. S., and Visscher, P. T. (2009). Processes of carbonate precipitation in modern microbial mats. Earth Sci. Rev. 96, 141–162. doi: 10.1016/j.earscirev.2008.10.005
Edgcomb, V. P., Bernhard, J. M., Summons, R. E., Orsi, W., Beaudoin, D., and Visscher, P. T. (2014). Active eukaryotes in microbialites from Highborne Cay, Bahamas, and Hamelin Pool (Shark Bay), Australia. ISME J. 8, 418–429. doi: 10.1038/ismej.2013.130
Farías, M. E., Rascovan, N., Toneatti, D. M., Albarracin, V. H., Flores, M. R., Poire, D. G., et al. (2013). The discovery of stromatolites developing at 3570 m above sea level in a high-altitude volcanic lake Socompa, Argentinean Andes. PLoS One 8:e53497. doi: 10.1371/journal.pone.0053497
Foster, J. S., Green, S. J., Ahrendt, S. R., Hetherington, K. L., Golubic, S., Reid, R. P., et al. (2009). Molecular and morphological characterization of cyanobacterial diversity in the marine stromatolites of Highborne Cay, Bahamas. ISME J. 3, 573–587. doi: 10.1038/ismej.2008.129
Foster, Z. S., Sharpton, T. J., and Grunwald, N. J. (2017). Metacoder: an R package for visualization and manipulation of community taxonomic diversity data. PLoS Comput. Biol. 13:e1005404. doi: 10.1371/journal.pcbi.1005404
Garby, T. J., Walter, M. R., Larkum, A. W., and Neilan, B. A. (2013). Diversity of cyanobacterial biomarker genes from the stromatolites of Shark Bay, Western Australia. Environ. Microbiol. 15, 1464–1475. doi: 10.1111/j.1462-2920.2012.02809.x
Goh, F., Allen, M. A., Leuko, S., Kawaguchi, T., Decho, A. W., Burns, B. P., et al. (2009). Determining the specific microbial populations and their spatial distribution within the stromatolite ecosystem of Shark Bay. ISME J. 3, 383–396. doi: 10.1038/ismej.2008.114
Goh, F., Barrow, K. D., Burns, B. P., and Neilan, B. A. (2010). Identification and regulation of novel compatible solutes from hypersaline stromatolite-associated cyanobacteria. Arch. Microbiol. 192, 1031–1038. doi: 10.1007/s00203-010-0634-0
Goh, F., Jeon, Y. J., Barrow, K., Neilan, B. A., and Burns, B. P. (2011). Osmoadaptive strategies of the archaeon Halococcus hamelinensis isolated from a hypersaline stromatolite environment. Astrobiology 11, 529–536. doi: 10.1089/ast.2010.0591
Golubic, S., and Hofmann, J. J. (1976). Comparison of Holocene and mid-Precambrian Entophysalidaceae (Cyanophyta) in stromatolitic alga mats: cell division and degradation. J. Paleontol. 50, 1074–1082.
Grotzinger, J. P., and Knoll, A. H. (1999). Stromatolites in Precambrian carbonates: evolutionary mileposts or environmental dipsticks? Annu. Rev. Earth Planet. Sci. 27, 313–358. doi: 10.1146/annurev.earth.27.1.313
Gudhka, R. K., Neilan, B. A., and Burns, B. P. (2015). Adaptation, ecology, and evolution of the halophilic stromatolite archaeon Halococcus hamelinensis inferred through genome analyses. Archaea 2015:241608. doi: 10.1155/2015/241608
Gunde-Cimerman, N., Plemenitas, A., and Oren, A. (2018). Strategies of adaptation of microorganisms of the three domains of life to high salt concentrations. FEMS Microbiol. Rev. 42, 353–375. doi: 10.1093/femsre/fuy009
Howe, A., and Chain, P. S. (2015). Challenges and opportunities in understanding microbial communities with metagenome assembly (accompanied by IPython Notebook tutorial). Front. Microbiol. 6:678. doi: 10.3389/fmicb.2015.00678
Huson, D. H., Beier, S., Flade, I., Gorska, A., El-Hadidi, M., Mitra, S., et al. (2016). MEGAN community edition - interactive exploration and analysis of large-scale microbiome sequencing data. PLoS Comput. Biol. 12:e1004957. doi: 10.1371/journal.pcbi.1004957
Inskeep, W. P., Macur, R. E., Harrison, G., Bostick, B. C., and Fendorf, S. (2004). Biomineralization of As(V)-hydrous ferric oxyhydroxide in microbial mats of an acid-sulfate-chloride geothermal spring, Yellowstone National Park. Geochim. Cosmochim. Acta 68, 3141–3155. doi: 10.1016/j.gca.2003.09.020
Jahnert, R. J., and Collins, J. B. (2012). Characteristics distribution and morphogenesis of subtidal microbial systems in Shark Bay, Australia. Mar. Geol. 30, 115–136. doi: 10.1016/j.margeo.2012.02.009
Joshi, N., and Fass, J. (2011). Sickle: A Sliding-Window, Adaptive, Quality-Based Trimming Tool for FastQ Files (Version 1.33). Available at: https://github.com/najoshi/sickle
Kawaguchi, T., and Decho, A. W. (2002). Isolation and biochemical characterization of extracellular polymeric secretions (EPS) from modern soft marine stromatolites (Bahamas) and its inhibitory effect on CaCO3 precipitation. Prep. Biochem. Biotechnol. 32, 51–63. doi: 10.1081/PB-120013161
Khodadad, C. L., and Foster, J. S. (2012). Metagenomic and metabolic profiling of nonlithifying and lithifying stromatolitic mats of Highborne Cay, The Bahamas. PLoS One 7:e38229. doi: 10.1371/journal.pone.0038229
Kulp, T. R., Han, S., Saltikov, C. W., Lanoil, B. D., Zargar, K., and Oremland, R. S. (2007). Effects of imposed salinity gradients on dissimilatory arsenate reduction, sulfate reduction, and other microbial processes in sediments from two California soda lakes. Appl. Environ. Microbiol. 73, 5130–5137. doi: 10.1128/AEM.00771-07
Kulp, T. R., Hoeft, S. E., Asao, M., Madigan, M. T., Hollibaugh, J. T., Fisher, J. C., et al. (2008). Arsenic(III) fuels anoxygenic photosynthesis in hot spring biofilms from Mono Lake, California. Science 321, 967–970. doi: 10.1126/science.1160799
Kurth, D., Amadio, A., Ordonez, O. F., Albarracin, V. H., Gartner, W., and Farias, M. E. (2017). Arsenic metabolism in high altitude modern stromatolites revealed by metagenomic analysis. Sci. Rep. 7:1024. doi: 10.1038/s41598-017-00896-0
Laval, B., Cady, S. L., Pollack, J. C., McKay, C. P., Bird, J. S., Grotzinger, J. P., et al. (2000). Modern freshwater microbialite analogues for ancient dendritic reef structures. Nature 407, 626–629. doi: 10.1038/35036579
Lindsay, M. R., Anderson, C., Fox, N., Scofield, G., Allen, J., Anderson, E., et al. (2017). Microbialite response to an anthropogenic salinity gradient in Great Salt Lake, Utah. Geobiology 15, 131–145. doi: 10.1111/gbi.12201
Logan, B. W., Hoffman, P., and Gebelein, C. D. (1974). “Alga mats, cryptalgal fabrics and structures, Hamelin Pool, Western Australia,” in Evolution and diagenesis of quaternary carbonate sequences, Shark Bay, Western Australia, ed. B. W. Logan (Tulsa, OK: American Association of Petrology), 140–194.
Louyakis, A. S., Gourle, H., Casaburi, G., Bonjawo, R. M. E., Duscher, A. A., and Foster, J. S. (2018). A year in the life of a thrombolite: comparative metatranscriptomics reveals dynamic metabolic changes over diel and seasonal cycles. Environ. Microbiol. 20, 842–861. doi: 10.1111/1462-2920.14029
Louyakis, A. S., Mobberley, J. M., Vitek, B. E., Visscher, P. T., Hagan, P. D., Reid, R. P., et al. (2017). A study of the microbial spatial heterogeneity of Bahamian thrombolites using molecular, biochemical, and stable isotope analyses. Astrobiology 17, 413–430. doi: 10.1089/ast.2016.1563
Lundberg, J., and McFarlane, D. A. (2011). Subaerial freshwater phosphatic stromatolites in Deer Cave, Sarawak—A unique geobiological cave formation. Geomorphology 128, 57–72. doi: 10.1016/j.geomorph.2010.12.022
Mobberley, J. M., Khodadad, C. L., Visscher, P. T., Reid, R. P., Hagan, P., and Foster, J. S. (2015). Inner workings of thrombolites: spatial gradients of metabolic activity as revealed by metatranscriptome profiling. Sci. Rep. 5:12601. doi: 10.1038/srep12601
Nutman, A. P., Bennett, V. C., Friend, C. R., Van Kranendonk, M. J., and Chivas, A. R. (2016). Rapid emergence of life shown by discovery of 3,700-million-year-old microbial structures. Nature 537, 535–538. doi: 10.1038/nature19355
Oksanen, J. (2017). Vegan: An Introduction to Ordination. R Package Version. Available at: https://cran.r-project.org/web/packages/vegan/vignettes/intro-vegan.pdf
Pace, A., Bourillot, R., Bouton, A., Vennin, E., Braissant, O., Dupraz, C., et al. (2018). Formation of stromatolite lamina at the interface of oxygenic-anoxygenic photosynthesis. Geobiology 16, 378–398. doi: 10.1111/gbi.12281
Papineau, D., Walker, J. J., Mojzsis, S. J., and Pace, N. R. (2005). Composition and structure of microbial communities from stromatolites of Hamelin Pool in Shark Bay, Western Australia. Appl. Environ. Microbiol. 71, 4822–4832. doi: 10.1128/AEM.71.8.4822-4832.2005
Pepe-Ranney, C., Berelson, W. M., Corsetti, F. A., Treants, M., and Spear, J. R. (2012). Cyanobacterial construction of hot spring siliceous stromatolites in Yellowstone National Park. Environ. Microbiol. 14, 1182–1197. doi: 10.1111/j.1462-2920.2012.02698.x
Playford, P. E., and Cockbain, A. E. (1976). “Modern algal stromatolites at Hamelin Pool, a hypersaline barred basin in Shark Bay, Western Australia,” in Stromatolites - Developments in Sedimentology, ed. M. R. Walter (Amsterdam: Elsevier), 389–411.
Proemse, B. C., Eberhard, R. S., Sharples, C., Bowman, J. P., Richards, K., Comfort, M., et al. (2017). Stromatolites on the rise in peat-bound karstic wetlands. Sci. Rep. 7:15384. doi: 10.1038/s41598-017-15507-1
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). 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
Reid, R. P., Visscher, P. T., Decho, A. W., Stolz, J. F., Bebout, B. M., Dupraz, C., et al. (2000). The role of microbes in accretion, lamination and early lithification of modern marine stromatolites. Nature 406, 989–992. doi: 10.1038/35023158
Riethman, H., Bullerjahn, G. S., Reddy, K. J., and Sherman, L. A. (1988). “Regulation of cyanobacterial pigment-protein composition and organization by environmental factors,” in Molecular Biology of Photosynthesis, ed. Govinedjee (Dordrecht: Springer), 229–257.
Rocap, G., Larimer, F. W., Lamerdin, J., Malfatti, S., Chain, P., Ahlgren, N. A., et al. (2003). Genome divergence in two Prochlorococcus ecotypes reflects oceanic niche differentiation. Nature 424, 1042–1047. doi: 10.1038/nature01947
Rodriguez-R, L. M., Overholt, W. A., Hagan, C., Huettel, M., Kostka, J. E., et al. (2015). Microbial community successional patterns in beach sands impacted by the Deepwater Horizon oil spill. ISME J. 9, 1928–1940. doi: 10.1038/ismej.2015.5
Ruvindy, R., White, R. A. III, Neilan, B. A., and Burns, B. P. (2016). Unravelling core microbial metabolisms in the hypersaline microbial mats of Shark Bay using high-throughput metagenomics. ISME J. 10, 183–196. doi: 10.1038/ismej.2015.87
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
Souza, V., Espinosa-Asuar, L., Escalante, A. E., Eguiarte, L. E., Farmer, J., Forney, L., et al. (2006). An endangered oasis of aquatic microbial biodiversity in the Chihuahuan desert. Proc. Natl. Acad. Sci. U.S.A. 103, 6565–6570. doi: 10.1073/pnas.0601434103
Suosaari, E. P., Reid, R. P., Abreau, T. A., Playford, P. E., Holley, D. K., McNamara, K. J., et al. (2016a). Environmental pressures influencing living stromatolites in Hamelin Pool, Shark Bay, Western Australia. PALAIOS 31, 483–496. doi: 10.2110/palo.2016.023
Suosaari, E. P., Reid, R. P., Playford, P. E., Foster, J. S., Stolz, J. F., Casaburi, G., et al. (2016b). New multi-scale perspectives on the stromatolites of Shark Bay, Western Australia. Sci. Rep. 6:20557. doi: 10.1038/srep20557
Swingley, W. D., Chen, M., Cheung, P. C., Conrad, A. L., Dejesa, L. C., Hao, J., et al. (2008). Niche adaptation and genome expansion in the chlorophyll d-producing cyanobacterium Acaryochloris marina. Proc. Natl. Acad. Sci. U.S.A. 105, 2005–2010. doi: 10.1073/pnas.0709772105
Theisen, C. H., Sumner, D. Y., Mackey, T. J., Lim, D. S., Brady, A. L., and Slater, G. F. (2015). Carbonate fabrics in the modern microbialites of Pavilion Lake: two suites of microfabrics that reflect variation in microbial community morphology, growth habit, and lithification. Geobiology 13, 357–372. doi: 10.1111/gbi.12134
Visscher, P. T., Reid, R. P., and Bebout, B. M. (2000). Microscale observations of sulfate reduction: correlation of microbial activity with lithified micritic laminae in modern marine stromatolites. Geology 28, 919–922. doi: 10.1130/0091-7613 (2000)28<919:MOOSRC>2.0.CO;2
Visscher, P. T., Reid, R. P., Bebout, B. M., Hoeft, S. E., Macintyre, I. G., and Thompson, J. A. (1998). Formation of lithified micritic laminae in modern marine stromatolites (Bahamas): the role of sulfur cycling. Am. Mineral. 8, 1482–1493. doi: 10.2138/am-1998-11-1236
Warden, J. G., Casaburi, G., Omelon, C. R., Bennett, P. C., Breecker, D. O., and Foster, J. S. (2016). Characterization of microbial mat microbiomes in the modern thrombolite ecosystem of Lake Clifton, Western Australia using shotgun metagenomics. Front. Microbiol. 7:1064. doi: 10.3389/fmicb.2016.01064
White, R. A., Chan, A. M., Gavelis, G. S., Leander, B. S., Brady, A. L., Slater, G. F., et al. (2016). Metagenomic analysis suggests modern freshwater microbialites harbor a distinct core microbial community. Front. Microbiol. 6:1531. doi: 10.3389/fmicb.2015.01531
Wong, H. L., Smith, D. L., Visscher, P. T., and Burns, B. P. (2015). Niche differentiation of bacterial communities at a millimeter scale in Shark Bay microbial mats. Sci. Rep. 5:15607. doi: 10.1038/srep15607
Keywords: stromatolite, microbial mat, comparative metagenomics, Hamelin Pool, Shark Bay
Citation: Babilonia J, Conesa A, Casaburi G, Pereira C, Louyakis AS, Reid RP and Foster JS (2018) Comparative Metagenomics Provides Insight Into the Ecosystem Functioning of the Shark Bay Stromatolites, Western Australia. Front. Microbiol. 9:1359. doi: 10.3389/fmicb.2018.01359
Received: 07 January 2018; Accepted: 05 June 2018;
Published: 25 June 2018.
Edited by:Dennis A. Bazylinski, University of Nevada, Las Vegas, United States
Reviewed by:Virginia Helena Albarracín, Center for Electron Microscopy (CIME), Argentina
Brendan Paul Burns, University of New South Wales, Australia
Copyright © 2018 Babilonia, Conesa, Casaburi, Pereira, Louyakis, Reid and Foster. 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 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.
*Correspondence: Jamie S. Foster, email@example.com