The Complete Genome and Physiological Analysis of the Eurythermal Firmicute Exiguobacterium chiriqhucha Strain RW2 Isolated From a Freshwater Microbialite, Widely Adaptable to Broad Thermal, pH, and Salinity Ranges

Members of the genus Exiguobacterium are found in diverse environments from marine, freshwaters, permafrost to hot springs. Exiguobacterium can grow in a wide range of temperature, pH, salinity, and heavy-metal concentrations. We characterized Exiguobacterium chiriqhucha strain RW2 isolated from a permanently cold freshwater microbialite in Pavilion Lake, British Columbia using metabolic assays, genomics, comparative genomics, phylogenetics, and fatty acid composition. Strain RW2 has the most extensive growth range for temperature (4–50°C) and pH (5–11) of known Exiguobacterium isolates. Strain RW2 genome predicts pathways for wide differential thermal, cold and osmotic stress using cold and heat shock cascades (e.g., csp and dnaK), choline and betaine uptake/biosynthesis (e.g., opu and proU), antiporters (e.g., arcD and nhaC Na+/K+), membrane fatty acid unsaturation and saturation. Here, we provide the first complete genome from Exiguobacterium chiriqhucha strain RW2, which was isolated from a freshwater microbialite. Its genome consists of a single 3,019,018 bp circular chromosome encoding over 3,000 predicted proteins, with a GC% content of 52.1%, and no plasmids. In addition to growing at a wide range of temperatures and salinities, our findings indicate that RW2 is resistant to sulfisoxazole and has the genomic potential for detoxification of heavy metals (via mercuric reductases, arsenic resistance pumps, chromate transporters, and cadmium-cobalt-zinc resistance genes), which may contribute to the metabolic potential of Pavilion Lake microbialites. Strain RW2 could also contribute to microbialite formation, as it is a robust biofilm former and encodes genes involved in the deamination of amino acids to ammonia (i.e., L-asparaginase/urease), which could potentially boost carbonate precipitation by lowering the local pH and increasing alkalinity. We also used comparative genomic analysis to predict the pathway for orange pigmentation that is conserved across the entire Exiguobacterium genus, specifically, a C30 carotenoid biosynthesis pathway is predicted to yield diaponeurosporene-4-oic acid as its final product. Carotenoids have been found to protect against ultraviolet radiation by quenching reactive oxygen, releasing excessive light energy, radical scavenging, and sunscreening. Together these results provide further insight into the potential of Exiguobacterium to exploit a wide range of environmental conditions, its potential roles in ecosystems (e.g., microbialites/microbial mats), and a blueprint model for diverse metabolic processes.

Our model microbialite ecosystem is Pavilion Lake a cold, oligotrophic ecosystem in southeastern British Columbia, Canada with dimictic circumneutral waters (median pH 8.3; mean calcium carbonate, 182 mg L −1 ), . Characterization of the limnology of Pavilion Lake is described in detail Lim et al. (2009). Resident microbialites are calcium carbonate-based thrombolites whose morphological features vary with depth , though all have thin (∼5 mm) microbial mats dominated by cyanobacteria.
While both heterotrophs and photoautotrophs (i.e., cyanobacteria) have been described and isolated from a range of microbial mats including microbialites, little work has been done on pigmented heterotrophic bacteria within microbialites. Firmicutes have been identified in the pigmented layers in microbial mats (Bottos et al., 2008;Lionard et al., 2012), and it is thought that carotenoids are responsible for their characteristic coloration (Nübel et al., 1999;Mueller et al., 2005;Klassen, 2010;Fisher et al., 2019). Fundamental questions about heterotrophic non-phototrophic arise which include (1) what the function of pigments in these microbialite associated heterotrophic bacteria is? and (2) what is the role of pigmented heterotrophs in ecosystem functioning within microbialites? To address these questions, we enriched and isolated >100 pigmented bacteria from Pavilion Lake freshwater microbialites, and grew them in the dark to select for heterotrophic or mixotrophic strains.
Among these pigmented heterotrophic enrichments included Exiguobacterium "strain RW2, " a Gram-positive member of the Firmicutes. Strain RW2 was isolated and enriched from a thrombolite at 20 m depth, where the water temperature remains around 4-10 • C throughout the year (Lim et al., 2009) (Figures 1A,B). Bacteria at this depth should be adapted to cold temperature, low phosphorus, and alkaline conditions. Recently, we have described another heterotrophic pigmented isolate Agrococcus pavilionensis strain RW1 (White et al., 2013b(White et al., , 2018a, that was co-isolated with strain RW2 from the same microbialite and enrichment. Members of this genus have a cosmopolitan distribution due to their highly adaptable physiology. Isolates of Exiguobacterium represent two major clades, clade I, which is comprised of cold-adapted strains including E. sibiricum 255-15 T (Rodrigues et al., 2008), and clade II, whose members range from temperate (e.g., E. aurantiacum DSM6208 T ) to hot environments (e.g., E. sp. AT1b) (Vishnivetskaya et al., 2009;Gutiérrez-Preciado et al., 2017), and include strain RW2. Due to the global distribution of Exiguobacterium, microbialites and/microbial mats appear to be another environment where they can adapt and colonize. Two Exiguobacterium strains have been directly isolated from microbialites: the present strain RW2 (formally pavilionensis) isolated from cold freshwater thrombolites (White et al., 2013a;Gutiérrez-Preciado et al., 2017), and strain S17, which was isolated from a warm (i.e., ∼20-24 • C) freshwater stromatolite in Lake Socompa (Ordoñez et al., 2013(Ordoñez et al., , 2015. Based on average nucleotide identity (ANI), RW2 appears to be the same species as two other strains; strain N139 and GIC31 (formally pavilionensis, Gutiérrez-Preciado et al., 2017), and it represents the type strain. Strain N139 was isolated from the water column of Leguna Negra which contains microbialites/microbial mats with extreme temperature variation ranging from −5 to 42 • C (Gomez-Javier et al., 2018;Mlewski et al., 2018). Remarkably, the third strain-GIC31 was isolated from glacial ice (Vishnivetskaya et al., 2009), potentially from a biofilm growing on top of the ice, but this is unknown. Another question arises (3) are strains RW2, N139, and GIC31 the same species? We use comparative genomics to determine the species placement of these three strains (i.e., RW2, N139, GIC31).  At the time of this writing, the genus Exiguobacterium has sixteen named species (Table 1), which were isolated from diverse environments including permafrost (Rodrigues et al., 2006), deep-sea vents (Crapart et al., 2007) and hot springs (Vishnivetskaya et al., 2009(Vishnivetskaya et al., , 2011. These widespread bacteria can tolerate of a wide range of pH (5-11), salinity (NaCl: 0-16%) and temperature (−12 • C to 55 • C) (White et al., 2013a) (Table 1). The first described member of the genus, E. aurantiacum DSM6208 T , was isolated from a potato processing plant (Collins et al., 1983). Consequently, isolates of Exiguobacterium spp. have been used as model organisms for understanding thermal adaptation to cold (i.e., −5 • C) (Rodrigues et al., 2008;Vishnivetskaya et al., 2009), hot temperatures (i.e., 55 • C) and heavy metal stress (Ordoñez et al., 2013(Ordoñez et al., , 2015. Currently, 60 genomes of Exiguobacterium are listed on NCBI, 17 of which are complete (Supplementary  Table S1). However, no complete genomes are available from microbialites, from which only a draft-genome is available (White et al., 2013a).
Here, we present the first such complete genome of Exiguobacterium isolated from modern microbialites. We also characterize strain RW2 using standard bacteriological and physiological testing, comparative genomics, and phylogenetics (i.e., 16S rRNA gene phylogeny) to describe its placement in the genus, and its broad adaptation to ranges of salinity, pH, and temperature. Lastly, from a genomic standpoint, we assess the potential of strain RW2 to contribute important metabolic functions to the microbialite community in Pavilion Lake.
We also assessed photoautotrophic and photoheterotrophic growth on MM9 minimal medium [20% glucose, 0.5% (w/v) tryptone, 0.25% (w/v) yeast extract], and on RCV medium (Beatty and Gest, 1981) in liquid medium and agar plates (1.5% agar w/v) grown both anaerobically, microaerophilically and aerobically with and without ambient sunlight on 8 h light and 16 dark cycle. No photoautotrophy or photoheterotrophy was observed.
Strain characteristics, including colony and cell morphologies, were determined by standard methods (Murray et al., 1994). Oxidase tests and biochemical enzyme assays and carbohydrate use were conducted in triplicate using API20E (BioMérieux) test strips following manufacturers instructions. Single isolated colonies from M-agar plates were washed via pelleting at 3,250 × g for 10 min three times in sterile distilled water then inoculated to API20E (BioMérieux) test strips, which includes a motility assay. Antibiotic susceptibility was determined by the Kirby-Bauer method using antibiotic disks on M-agar plates in triplicate (Collee et al., 1996).  Farrow et al., 1994 Levine andSoppeland, 1926 classified "Flavobacterium acetylicum" renamed "Brevibacterium acetylicum" then reclassified by Farrow et al. (1994) E. antarcticum Microbial mat/biofilm Fruhling et al., 2002;Carneiro et al., 2012 Lake Fryxell mat  Triplicate experiments in liquid culture were used to induce biofilm formation and measure timing (M-broth 30 • C at pH 7 and 1% NaCl w/v) at low shaking (100 rpms) and without shaking. A positive biofilm in strain RW2 was a top film appeared over medium within the flask.

Light and Scanning Electron Microscopy (SEM)
Exponentially growing cells were harvested at 26 h after being transferred to liquid M-medium and were viewed by light microscopy under oil immersion at 100×. Flagella straining was completed as per Kodaka et al. (1982) then viewed by light microscopy after 26 h of growth. For SEM, a culture in stationary phase ∼72 h after being transferred to liquid M-medium was pelleted at 3,250 × g for 10 min, then the M medium was exchanged and cells were fixed in 2.5% glutaraldehyde in PBS (137 mM NaCl, 2.7 mM KCL, 10 mM Na 2 HPO 4 2H 2 O, 2 mM KH 2 PO 4 , pH 7.4) solution for 30 min on ice. The fixed culture was filtered onto a 0.2 µm pore-size Supor polycarbonate membrane (Pall Port Washington, NY, United States). Cells on the filter were washed with PBS and post-fixed in 1% OsO 4 for 1 h. The cells and filter were passed through a graded ethanol series (25,50,70,95, 100%) at 10 min intervals, and in 100% ethanol were critical-point dried with CO 2 . A sputter coater applied 5 nm of gold/palladium alloy onto the cells before imaging by SEM using a Hitachi S4700 microscope.

Phylogenetic Analysis of 16S rRNA Gene Sequences
Thirty-two reference sequences were downloaded from NCBI representing all named typed ( T ) strains of genus Exiguobacterium, with three typed strains of Bacillus spp. as outgroups. The whole-genome assembled 16S rRNA gene sequence was used for phylogenetic analysis. Alignments of 16S rRNA gene sequences were completed in MAFFT (v. 7.310) using the options (-localpair -maxinterate 1000) which is iterative refinement method incorporating local pairwise alignment information (L-INS-i) providing the most accurate alignment (Katoh et al., 2002). A maximum likelihood phylogenetic tree was constructed using IQ-TREE (v. 1.6.1) (Hoang et al., 2017) with a total of 1,000 bootstrap replicates using UFBoot2 (Trifinopoulos et al., 2016), and visualized with SeaView (v. 4) (Gouy et al., 2010). Strain RW2 16S rRNA gene sequence and alignments are available on github.com/strain_RW2.

Whole Genome Assembly and Annotation
DNA extraction, Illumina library preparation and sequencing, data cleaning, phiX spike-in removal, and draft genome assembly are described by White et al. (2013a). The genome was closed into a single circular contig (i.e., chromosome) without plasmids using comparative genome ordering, alignment, mapping, and manual editing. For comparison, scaffolding and contig ordering of the genome of strain AT1b was compared to strain RW2 which was downloaded from NCBI. Strain RW2 contigs were ordered and aligned in progressiveMauve to strain AT1b to visualize genome order (Darling et al., 2010;Vishnivetskaya et al., 2011). The remaining gaps between contigs were closed by recursive alignments using Mauve. The ordered contigs with overlaps were merged into a single circular contig using the EMBOSS union script (Rice et al., 2000). To estimate sequencing depth represented in the assembly for coverage estimation, reads used in the assembly were mapped back to the final circular genome using Bowtie2 with the very sensitive local option (Langmead and Salzberg, 2012). The Bowtie2 read mapping output file (Sam file) was visually inspected by the Tablet program (Milne et al., 2013). The contigs that were screened for overlaps and read mapping depth (>10×) were then merged manually, based on the reference genome of AT1b.
Annotation was conducted on the RAST annotation server using the Glimmer-3 option and standard RAST against the FIGfam database release-70 (Aziz et al., 2008). Metabolic pathways were predicted in strain RW2 using MetaPathways which uses pathway tools (Paley and Karp, 2006;Konwar et al., 2013;Caspi et al., 2014), against MetaCyc/BioCyc databases via a built-in LAST (Local Alignment Search Tool) (Kiełbasa et al., 2011) for homology matches of ≥180 bp ORFs (protein-coding open reading frames) and ≥50 alignment score. A metabolic model for strain RW2 was predicted using ModelSEED (Devoid et al., 2013). The Prokaryotic Genome Annotation Pipeline (PGAP) at NCBI was also used to compare to RAST annotation for strain RW2 (Tatusova et al., 2016). We manually curated the annotations produced by RAST and NCBI (PGAP)-retaining only annotations that were a consensus between both programs.

Comparative Genomic Analysis
All genomes available for Exiguobacterium were uploaded to RAST, for synteny, AAI via RAST annotations, and metabolic model comparisons in modelSEED. Metapathways was also used to predict pathways across the genomes using BioCyc/MetaCyc. To display genomic features a circular genome plot of strain RW2 was constructed using Cgviewer server 1 (Grant and Stothard, 2008, version V 1.0, date accessed Jan 1st, 2018). Genomes of strains RW2, S17, GIC31, and N139 were compared in CGViewer using tBLASTx (e-value of 1e −5 , with 70% identity cutoff, and minimum 50 bp overlap) and then displayed in the CGviewer genome plot server (Grant and Stothard, 2008). Synteny plots were completed in the RAST server module using a BLASTbased dot plot format. Average amino-acid identity (AAI) and Reciprocal Orthology Score Average (ROSA) via functional gene similarities from RAST server annotations were calculated using with a web-based tool 2 (Krebs et al., 2013, date accessed January 1st, 2018.

Metagenomic Read Recruitment
The presence of strain RW2 within Pavilion Lake sediments, microbialites, and water column was previously described in a metagenomic study by White et al. (2016). We compared our genomic data to their published metagenomic data, by loading their predicted protein dataset onto MG-RAST (Meyer et al., 2008), and using RefSeq (protein level) annotations (using BLAT parameters, with matches of ≥60% similarity cutoff, ≥15 bp overlap, and minimum E-value of 10 −5 ) to search for signatures of strain RW2 in their Pavilion Lake metagenomes. Reads assigned to Exiguobacterium from the metagenome were mapped to the strain RW2 using Bowtie2 with the very sensitive local option (Langmead and Salzberg, 2012). These results were subject to ANOVA statistical testing using the Statistical Analysis of Metagenomic Profiles package (STAMP) (Parks and Beiko, 2010). The program FR-hit was then used to assign metagenomic read recruitment of the 7.5 million Metagenomic sequences (250 bp paired-end MiSeq reads) to Exiguobacterium genomes from strain RW2, S17, and AT1b, using default parameters with a minimum identity > 70% and an E-value > 1e −5 (Niu et al., 2011). The recruitments were then visualized with the R library ggplot2 (Wickham, 2016).

Data Availability
Strain RW2 is listed at NCBI under bioproject accession PRJNA208114 and biosample accession SAMN02471612. The Ray assembly for the draft genome is listed under DDBJ/EMBL/ GenBank under the accession number ATCL00000000, and RefSeq under NZ_ATCL00000000.1, as well as NCBI under GCA_000416965.1. Lastly, the complete genome, scripts, and supplementary data are posted on github.com/raw937/ Strain_RW2/.

Classification of RW2 Related Strains Using Phylogenetics and Whole-Genome Analysis
Using 16S rRNA gene-based phylogeny, whole-genome ANI/ AAI and digital DNA-DNA hybridization we will determine whether RW2, GIC31, and N139 are the same species. Currently, strains GIC31 and N139 are missing data required for classical bacteriological species determination via the IJSEM which include experimental fatty acid analysis, quinone characterization, peptidoglycan structure analysis, and DNA-DNA hybridization assays. However, a recent proposal from IJSEM suggests that genome-based ANI/AAI and 16S rRNA gene-based phylogeny can replace the experimental DNA-DNA hybridization and be used for naming/typing strains in the future (Chun et al., 2018).
We generated a maximum-likelihood phylogeny of 16S rRNA gene sequences from all fourteen typed strains of Exiguobacterium spp., fifteen sequences representing strains with complete genomes, and three typed strains of Bacillus spp. as outgroups. This positioned strain RW2 in clade II, along with E. sp. AT1b and E. aurantiacum DSM6208 T (Figure 2). E. aurantiacum DSM6208 T represents the first isolate of Exiguobacterium, from potato effluent (which contains soil) (Collins et al., 1983). Within clade II are two subclades; clade IIa and clade IIb, together containing isolates from a wide range of environments (Figure 2). Clade IIa stem from hot springs and both shallow and deep-sea marine ecosystems, whereas isolates of clade IIb stem terrestrial habitats including soil, microbialites (including strain RW2), and freshwater ecosystems (Figure 2). Clade IIb is well-supported, suggesting that it represents a true radiation into terrestrial habitats including lakes (i.e., strains N139, SH31, S17, RW2), glacial ice (i.e., strain GIC31), microbialites (i.e., strains S17, RW2), canals (i.e., strain HUD), and even the microbiomes of brine shrimp, which dwell in hypersaline vernal pools (i.e., strain DSM16483 T ) and canals (i.e., strain HUD) (Figure 2). However, the subclade that contains strains RW2, N139, GIC31, and HUD is not well resolved (Figure 2). Strain HUD is currently classified by NCBI as a strain of E. mexicanum, but this is challenged by our 16S rRNA gene phylogeny, given that it does not branch with the type strain of E. mexicanum (strain DSM16483) (Figure 2). The genome of the type strain of E. mexicanum DSM16483 T is currently not available to compare to strain HUD or other strains related to strain RW2. There is high bootstrap support for a clade containing strain RW2 (as the basal member) along with strains SH31 and S17. By contrast, strains HUD, GIC31, and N139 branch earlier than this clade, where the tree topology is poorly resolved (i.e., the branches are polytomic) (Figure 2). The 16S rRNA gene of RW2 is 99% similar to GIC31 and N139 based on BLASTN, whereas N139 is 100% similar to GIC31. More representative 16S rRNA gene sequences are needed to resolve the polytomy within clade IIb (Figure 2).

Morphology, Growth, and Characteristics of Strain RW2
Strain RW2 is a Gram-positive aerobe that is facultatively anaerobic. It subsists heterotrophically with no evidence of photoheterotrophy. After 48 h of growth at 30 • C on M-agar, it forms bright orange, smooth, circular colonies that are typically 3 to 4 mm in diameter ( Figure 1C and Table 3). In liquid M-medium at 30 • C, the cells were coccoid after 24-48 h in FIGURE 4 | Genome plot of strain RW2 (bp) using CGviewer. Genome key (left corner): starts with the innermost ring which is a genome ruler followed by GC skew (purple/green) and ends with two outer rings which contain Protein coding ORFs (blue), tRNAs and rRNAs (red). tBLASTx was used compare RW2 with an e-value 1e -5 , alignment cutoff -50 bp, identity cutoff -70%. Blast 1 is strain N139 vs. RW2, Blast 2 is strain S17 vs. RW2, and Blast 3 is strain GIC31 vs. RW2. logarithmic growth and irregularly shaped after 72 h in stationary phase, coincident with the beginning of biofilm formation under shaking (100 rpm) (Figures 1D,E). Biofilms first appeared within ∼72 h, and their growth culminated at 8 days for 30 • C under shaking (100 rpm).
No flagella were observed in the electron microscopy (SEM) or under standard light microscopy (Figures 1D,E). Strain RW2 had a positive motility assay provided by API20E kit. The genome annotation of strain RW2 predicts a complete flagellum biosynthesis pathway (including the fih, fli, fig gene operons). Flagella staining and motility testing was positive, meaning that strain RW2 is motile with peritrichous flagella.
Strain RW2 grows at a wide range of temperatures, added salinity, and pHs. Ranging from 4 to 50 • C, strain RW2 grows at the broadest range of temperatures currently reported for any isolate of Exiguobacterium genus) ( Table 3) though we were not able to test its growth at subzero temperatures. Likewise, strain RW2 exhibits the broadest range of pH tolerance reported within the genus, from a pH of 5 to pH 11 (Table 3). We attempted to grow strain RW2 on minimal medium to test the absolute requirement of NaCl required for growth, but we were unsuccessful. We tested overall NaCl tolerance for growth by directly adding NaCl from 0 to 16% (w/vol) to M-agar/broth. Strain RW2 only grew at 0 and 7% (Table 3). Nevertheless, strain RW2 has a higher salt tolerance than other non-marine isolates, which do not grow above 7% added NaCl (Table 3). Interestingly, strain RW2 grows beyond the range of the temperature, pH, and salinity expected at 20 m depth in Pavilion Lake microbialites, where temperatures range from 4 to 10 • C, and pH is very stable at 8.1-9.1 (Lim et al., 2009). This degree of abiotic tolerance makes RW2 unique among isolates of the genus. Cut-offs for same species with >80% probability based on logistic regression with >70% for species determination.  pH for growth  50  30  41  40  41  30 40 37 All growth temperature measurements for strain RW2 were taken after 2 days, except at 4 • C which took 10 days to grow. Non-strain RW2 data are from Chaturvedi et al. (2008). no growth (−), strong growth (+), weak growth (W), data not available (NA

Genome Properties and Comparative Genome Analysis
The initial draft genome of strain RW2 was not complete and fragmented on 23 contigs (White et al., 2013a) (Table 4). With progressiveMauve, mapping, and manual finishing, we closed the strain RW2 genome into a single circular chromosome. This chromosome is 3,019,504 bp in length, with a GC content 52.05%, and no plasmids ( Table 4). This is similar to the genome assembly of strain GIC31, which is likewise a single circular contig, with high coverage that suggests no evidence of plasmids ( Table 4).
The draft and completed genomes of Exiguobacterium spp. were similar regarding genome size, GC content, and the number of protein-coding genes. Within the genus, genome sizes range from 2.82 to 3.16 Mb, with GC content ranging from 47.5 to 53%, containing from 2,941 to 3,323 genes that encode 2,772 to 3,332 predicted proteins (Table 4). Strain RW2 had lower numbers of rRNAs, tRNAs, and had 17 pseudogenes predicted ( Table 4). Both S17 and 8-11-1 have lower synteny, possibly as a result of the fact that these are less complete genome assemblies (>25 contigs) ( Table 4 and Supplementary Figure S1). Among members of clade IIa, strain RW2 had the highest genome synteny to strain AT1b (Supplementary Figure S1). MetaCyc/BioCyc pathway predictions using Metapathways suggest that >90% of the pathways are shared even among genomes with <65% AAI (Supplementary Figure S2). Strain RW2 has the high AAI identity to strain S17 (87.76%, Supplementary Figure S2). It also  shared over 50% of RAST functions with more distant relatives in the phylum Firmicutes (e.g., Planococcus halocryophilus and Bacillus halodurans, (Supplementary Figure S1). The rarity of unique genes in strain RW2 suggests that it has a restricted accessory genome (i.e., pan-genome) while a large core genome is highly conserved across strains of Exiguobacterium.

Potential Genomic Signatures of Physiological Adaptability
Strain RW2 has physiological adaptability to a range of thermal, salinity and pHs which we further inquired the genomic basis for such growth parameters. Cold shock proteins encoded by cspA to cspI are molecular chaperones that help to regulate protein-folding under cold conditions and prevent cryodamage (Kawahara, 2008). Interestingly, strain RW2 and the other clade II members seem to lack cspA and cspB genes, while members of clade I have multiple copies of cspA and cspB (e.g., E. sibiricum strain 255-15 and E. antarcticum B7). We compared the strain RW2 genome to a distantly related member of the firmicutes known as Planococcus halocryophilus which can grow at −15 • C (Mykytczuk et al., 2013). Planococcus halocryophilus has only one copy of cspA, unlike E. sibiricum strain 255-15 and E. antarcticum B7 which have multiple copies of cspA. This suggests growth below −5 • C mechanisms are required. By contrast, strain RW2 is surprisingly thermotolerant and grows at >45 • C (White et al., 2013a) ( Table 3). Strain RW2 genome contains a complete heat shock gene cluster (dnaJ, dnaK, and GrpE), which encodes for chaperones that prevent denaturation and aggregation of proteins at high temperatures (Feder and Hofmann, 1999).
In regards to osmoregulation, choline and betaine uptake/biosynthesis pathways are effective systems for regulating osmolyte concentrations, including that of NaCl. Uptake pathways include the genes bet, opu, and proU, which encode ABC transport proteins (Sleator and Hill, 2002). We again compared the strain RW2 genome to P. halocryophilus as the latter can tolerate 18% added NaCl (Mykytczuk et al., 2013). Strain RW2 has no bet genes, one copy of opuA (type AA, AB, AC only), no opuC, two copies of opuD, and one copy of proU (proX only); whereas P. halocryophilus has one copy of betT, has two copies of opuA (type AA, AB, AC), one copy opuC, five copies of opuD, multiple copies of proU (two copies proV/X, and one proW) based on RAST annotations.
In regards to pH regulation, strain RW2 genomic evidence of regulation includes antiporters-including the arginineornithine antiporter (arcD)-can facilitate survival under acidic conditions (Fulde et al., 2014). Strain RW2 has four copies of the arcD antiporter-the same number found in E. marinum DSM16307 T and E. acetylicum DSM20416 T , which dwell in saline marine environments. By contrast, close relatives of strain RW2, such as strain AT1b and S17 have only two copies of arcD; while P. halocryophilus has zero copies based on RAST annotation. Strains of the alkaliphilic Bacillus genus commonly regulate high pH with the nhaC Na + /K + antiporter (Ito et al., 1997). Strain RW2 has three copies of nhaC gene, which is the same in AT1b (clade II) and P. halocryophilus, whereas E. sibiricum 255-15 (clade I) which only has two copies based on NCBI annotation.
The mrp and tet(L) gene cluster antiporters which also regulate high pH (Padan et al., 2005), were not found in strain RW2 or other members of the genus.
Phospholipid fatty acid profile shifted across added NaCl and pH growth conditions (Table 6 and Figure 5). Similar to the temperature profiles, the notable differences were found in i-iso-C16:1 5, iso-C17: 1 5, brC17:1 5, and iso-C18:1 5 ( Table 6 and Figure 5). Under the acidic pH 5, iso-C17:1 5 was found to comprise 7.0 ± 0.9 mol % of the total PLFAs; this value increased to 15.5 ± 1.6 mol % by pH 11 (Table 6). This shift was the major influence on the degree of unsaturation, which increased between pH 7 and pH 11 but displayed little change between pH 5 and pH 7 ( Table 6). The monoenoic PLFAs were also found to increase between pH 7 and pH 11, while there were no differences observed in these PLFAs between pH 5 and pH 7 ( Table 6).
Salinity (via added NaCl) had limited effects on the composition of the PLFA profile. The total mol % of branched monoenoic PLFAs displayed a minor decrease between 0 NaCl and 1% NaCl but remained constant after that. These values decreased to 0.12 and 2.1, respectively, at 1% NaCl and remained similar at 7% NaCl (Table 6 and Figure 5). C14:0 and C17:0 were only detected at 7% NaCl ( Table 6). On the contrary, the mol % of all branched PLFAs was the lowest at 7% NaCl. The average fatty acid chain length displayed no large changes with altered salinity ( Table 6 and Figure 5).

Resistance of Heavy Metals and Antibiotics Including Community Presence
After sequencing the metagenomes of Pavilion Lake , we questioned whether strain RW2 was part of the Pavilion Lake microbialite microbiome. With metagenomic recruitment analysis, we uncovered the presence of strain AT1b and unclassified Exiguobacterium spp., with the latter being significantly (p-value < 0.01) enriched within the microbialite, compared to the surrounding sediment and water column (Supplementary Figure S3A). We found no Exiguobacterium spp. reads present in the sediment metagenomes of Pavilion Lake (Supplementary Figure S3A). This taxonomic resolution of MG-RAST was not sufficient to detect whether the hits were from strain RW2 or another Exiguobacterium spp. To calculate the presence and abundance of strain RW2, we recruited metagenomic reads (from the 20 m depth microbialite sample against genomic databases from strains RW2, AT1b, and S17. Strain RW2 recruited more reads (0.46% of reads) than S17 (0.33%) or AT1b (0.25%-with its genome achieving nearly complete coverage-as well as more reads with 100% identity hits (Supplementary Figure S3B). We further checked whether any of the 100-95% identity hits recruited from the metagenome were related to heavy metal resistance genes present in the genome of strain RW2, using tBLASTx (1e −3 ). We found positive heavy metal gene recruitment to mercuric reductases, arsenic resistance pumps, multidrug and toxin efflux (MATE) resistance pumps, chromate transporters and cadmium-cobalt-zinc resistance genes.
By contrast, we could not determine which taxa were responsible for contributing antibiotic resistance genes to the Pavilion Lake microbialite metagenomes, other than Proteobacteria , which were affiliated with gyrB antibiotic resistance genes, known to encode multifunctional DNA repair proteins. Nevertheless, the assembled genome of strain RW2 encodes four copies of MATE resistant pumps, vancomycin B-type resistance proteins, tetracycline resistance MFS pumps, as well as four antibiotic biosynthesis monooxygenases, which are involved in many diverse processes including the biosynthesis of antibiotics.
Little is known about antibiotic sensitivity and metabolism in Exiguobacterium spp. Other than E. soli, E. indicum and E. acetylicum, isolates of Exiguobacterium spp. are sensitive to clindamycin, as is strain RW2 (Supplementary Table S2).
Additionally, strain RW2 was tested against eleven other antibiotics (Supplementary Table S2) and found only to possess resistance against sulfisoxazole (300 µg). The sul and dhf genes are found in organisms resistant to sulfisoxazole (Barman et al., 2010). This resistance may stem from the genomic presence of well-known sulfisoxazole resistance genes sul and dhf genes (encoding dihydropteroate synthase and dihydrofolate reductase), within the genomes of strain RW2, E. indicum, and E. acetylicum. While genes for vancomycin and tetracycline resistance were present in the genome of strain RW2, testing revealed the strain was sensitive to both antibiotics (Supplementary Table S2).

Carbohydrate, Amino Acid, and Nitrogen Metabolism
Biochemical tests indicated that strain RW2 uses a variety of monosaccharides, polysaccharides, and amino acids as carbon sources ( Table 7). Pathways for lactose uptake were present, but only a lacZ (β-galactosidase) for lactose utilization appeared to be present (Table 7). Strain RW2 could not use D-mannitol, Dinositol, D-sorbitol, L-rhamnose, D-melibiose, and L-arabinose as sole carbon sources (Table 7). D-mannitol utilization was predicted by the genome, yet it was not assimilated as a sole carbon source, suggesting that only fermentation occurs ( Table 7). Amygdalin was utilized by strain RW2 (Table 7), as is the case for several members of the genus. Although amygdalin metabolism was not specifically predicted in our RAST analysis of the genome, many general function betaglucosidases were present in the genome, potentially lending this function. The genome of strain RW2 also appears to possess metabolic potential to perform various forms of fermentation including mixed acid, lactate, and acetyl-CoA fermentation to butyrate. Lastly, our analysis suggests diverse substrates including chitin, N-acetylglucosamine, maltose, maltodextrin, trehalose, glycerol, glycerol-3-phosphate, glycogen, fructose, gluconate, ribose, fructooligosaccharides (FOS), raffinose and deoxyribose can be used, although none were tested as growth substrates here.
Our genome analysis suggests that strain RW2-like other members of the genus-is auxotrophic for branched-chain amino acids, as biosynthetic pathways for isoleucine, leucine, and valine were predicted to be absent, while complete degradation pathways were present. Moreover, arginine dihydrolase was All growth measurements for strain RW2 were taken after 24 h in triplicate. Non-strain RW2 data are from Chaturvedi et al. (2008). All are listed by strain name. All strain publications are further described in Table 1. functional in strain RW2; whereas, tryptophan deaminase, lysine, and ornithine decarboxylases were not (Table 7). Further growth substrates should be tested to elucidate carbohydrate and amino acid utilization in strain RW2 further. Strain RW2 genome is predicted to encode genes that generate ammonium from amino acids. An L-asparagine biosynthetic pathway is predicted to encode genes for aspartate aminotransferase, aspartate-ammonia ligase, and L-asparaginase. The enzyme L-asparaginase is known to hydrolyze L-asparagine to L-aspartate and ammonia, though we did not chemically validate its function in strain RW2. By contrast, urease-another enzyme capable of ammonification, via the catabolism of urea into ammonium and CO 2 -was not found in the genome of strain RW2. Interestingly, cells of RW2 nonetheless measured as positive on a biochemical urease test, conducted in triplicate ( Table 7). It is possible that a poorly annotated gene encodes this urease function; since there are five copies of unspecified aminohydrolases within the genome of strain RW2, and given that aminohydrolases represent the enzymatic superfamily to which urease belongs (EC 3.5.1.5).
is commonly found in members of the Firmicutes (Klassen, 2010). Our comparative genomic analysis indicates that complete pathways for carotenoid C 30 biosynthesis were present in all the published Exiguobacterium spp. genomes, including that of strain RW2 (Figure 6). Relevant enzymes include diapophytoene synthase (crtM, annotated as phytoene synthase), diapophytoene desaturase (crtN, annotated as phytoene desaturase), and diapophytoene desaturase (crtNb, annotated as the second copy of phytoene desaturase). Similar to the case in Halobacillus halophilus (Köcher et al., 2009) (Figure 6), all genes for the pathway are located on a single operon, specifically on a span of 5,673 bp, stretching from position 39,293 to 44,966 bp on the genome.

Exiguobacterium Coloration Is Based on C 30 Carotenoids With Potential Functions
Based on our comparative genomic analysis of draft and completed genomes from NCBI, we predict that a 4-4' diapocarotenoids C 30 carotenoid-based biosynthetic pathway is responsible for the orange pigmentation of Exiguobacterium colonies, with diaponeurosporene-4-oic as the final product. As reported in Halobacillus halophilus, genes encoding this pathway are arranged as an operon, and genetic ablation of the terminal gene (for diapophytoene desaturase) results in the accumulation of diaponeurosporene-4-oic acid and consequent orange pigmentation (Köcher et al., 2009).
To our knowledge, a hypothesis has not been advanced to explain the function of pigments in the genus Exiguobacterium. Low concentrations of dissolved organic carbon (DOC) in Pavilion Lake allow for high UV penetration, and carotenoids might provide a "sunscreen" function to cope with these conditions (Lim et al., 2009;Lionard et al., 2012). Three tested strains of Exiguobacterium are known to use photoenzymatic repair (PER) and nucleotide excision repair (NER) to recover from UV-B photodamage, but all had some form of UV-B resistance beyond PER and NER (Gutiérrez-Preciado et al., 2017). We hypothesize that 4-4' diapocarotenoids C 30 carotenoids protect members of Exiguobacterium from UV-B damage, given the well-known role of carotenoids quenching and scavenging reactive oxygen species, releasing excessive light energy, and sunscreening (Wada et al., 2013). Further experimentation is needed to confirm if carotenoids indeed protect firmicutes against UV radiation.

Potential Contribution to Microbialite Formation and Maintenance
While strain RW2 has many of the physiological features of an environmental generalist, it was specifically enriched in the microbialite community at 20 m depth in Pavilion Lake. Specifically, its reads were absent in the sediment metagenome sampled nearby and are poorly represented in the metagenome sampled from the water column overhead. Moreover, reads from strain RW2 comprised a higher proportion of total microbialite reads than those of strain S17 in Socompa Lake microbialites and strain AT1b in a hyperthermal microbial mat (Ordoñez et al., 2013(Ordoñez et al., , 2015. Together these data suggest that strain RW2 is a numerically significant component of the Pavilion Lake microbialites. Strain RW2 was found in lower amounts on the filters (i.e., water column metagenomes), and wasn't as prevalent as in the microbialite metagenomes. Filters can contain small amounts of microbialite microbial mat as they were sampled near the microbialites at depth; however, the water column metagenomes are distinctly different via PCA analysis from the microbialite metagenome .
Since we do not know how long strain RW2 has lived on freshwater microbialites, it is unclear how and to what extent it has co-evolved with this microbial community. Even in the absence of evolutionary specialization, there are several ways in which strain RW2 might contribute to microbialite formation and maintenance. For instance, biofilm formation has previously been nominated as a survival mechanism in nutrient-poor conditions (Shao et al., 2013), and this might apply to the cold, oligotrophic conditions from which strain RW2 was isolated (Lim et al., 2009). By adhering to sediment and precipitated carbonates, Firmicutes such as Exiguobacterium could, in turn, contribute to the coherence and spatial structure of microbialites like those found in Pavilion Lake and Laguna Negra Gomez-Javier et al., 2018). But biofilm formation is not necessarily a consequence of coevolution with microbialites since other strains form biofilms in the absence of microbialites, such as Exiguobacterium strain GIC31, which may grow directly on glacial ice. A few strains of Exiguobacterium isolates have been found to be endemic to Cuatro Ciénegas Basin which microbialites are found (Rebollar et al., 2012).
Several metabolic pathways predicted in strain RW2 might also contribute to microbialite formation by promoting inducing carbonate precipitation via urease or deamination of amino acids. For instance, the enzyme L-asparaginase generates ammonium as a by-product of the deamination of amino acids, thereby increasing alkalinity, which favors the precipitation of carbonate minerals (Castanier et al., 1999). Interestingly, our biochemical assay indicated that urease is also active in strain RW2; though we found no evidence of urease genes within its genome, nor that of any other Exiguobacterium strain. This suggests either (A) the presence of an uncharacterized "orphan enzyme" (perhaps any of the five amidohydrolases predicted in the genome of strain RW2), or (B) that a known enzyme is capable of an undescribed side reaction (Sorokina et al., 2014;Montgomery et al., 2015). Amidohydrolases belong to the same protein superfamily as urease; which is a future target for further characterization and investigation. In any case, urease-like asparaginase-can induce calcium carbonate precipitation, and it is found in carbonate precipitating members of Bacillus spp. (Hammes et al., 2003), though its actual involvement in biomineralization remains to be demonstrated.
More compelling is the possibility that strain RW2 contributes other metabolic functions-namely heavy-metal metabolismsince these Exiguobacterium genes are highly represented in the microbialite metagenome of Pavilion Lake. For reasons that are still unclear, heavy-metal metabolism appears to be an accessory feature of microbialite communities (Ruvindy et al., 2016;White et al., 2016;Kurth et al., 2017;Gomez-Javier et al., 2018). Generally, heavy metals (such as arsenic) limit microbial growth. However, we find three examples of thriving modern microbialite ecosystems in the presence of high arsenic; the microbialites of Laguna Brava (Sancho-Tomás et al., 2018), Socompa Lake (Kurth et al., 2017), andLaguna Negra (Gomez-Javier et al., 2018). The extracellular polymeric substances (EPS) within cyanobacterial microbialite mats and biofilms bind, concentrate, and remove heavy metals from the water column (Arp et al., 1999). Any heterotroph that eliminates and detoxifies heavy metals as a byproduct of its metabolism would benefit the entire microbialite microbiome and would presumably be rewarded with more substrates for growth. Strangely, Pavilion Lake has low levels of zinc (0.01-0.03 mg L −1 ) and undetectable levels of cobalt, copper, chromium, arsenic, and cadmium (Lim et al., 2009).
Genes predicted in the resistance of mercury, arsenic, chromate, cadmium, cobalt, zinc resistance were all found in Pavilion Lake microbialite microbiome based on recruitment of their metagenomic reads to the genome of strain RW2. While strain RW2 hasn't been tested for arsenic resistance; strain N139 can grow at 100 mM of arsenate (As [V]), and 2.5 mM of arsenite (As [III]) (Ordoñez et al., 2015). All three strains (RW2, GIC31, N139) have no acr3 based on genome predictions; whereas S17 does and can grow at 10 mM As [II] (Ordoñez et al., 2015). We would predict that strain RW2 would have similar arsenic resistance as strain N139, high As [V] resistance but low As [III] due to the lack of the acr3 gene.
Moreover, the genome of Agrococcus pavilionensis strain RW1-which was isolated from the same microbialite as strain RW2-also carried genes for heavy-metal resistance genes (White et al., 2013a). One possible explanation is that these genomic signatures are relicts from a time when Pavilion Lake had higher levels of dissolved heavy metals, as likely occurred in the early 20th century when mining operations were active in the region (Stevenson, 1940).
Another possibility is that metal resistance genes have other secondary alternative functions, as in Rhodobacter sphaeroides is known to have arsenic resistance genes expressed under high salt stress (Tsuzuki et al., 2011). For instance, the heavy metal resistance proteins genes in strain RW2 are retained in a low heavy metal environment because of the act as can also detoxification of other substrates (e.g., antibiotics). As such, heavy metal metabolism can drive co-selection of antibiotic resistance when aquatic systems that are impacted by agriculture or antibiotics introduced through wastewater and agricultural runoff anthropomorphic means (Seiler and Berendonk, 2012). Resistance in heavy metals has conferred resistance to antibiotics in a complex microbiome (e.g., chicken guts) (Nisanian et al., 2014). It is not clear whether this is the case for microbialites in Pavilion Lake, though strain RW2 also contributes antibiotic resistance genes to the metagenome of these microbialites, including those for sulfisoxazole resistance .

Classification of the Strains as the Same Species
Using 16S rRNA gene sequence phylogeny, whole-genome ANI/AAI, and digital DNA-DNA hybridization; we confirmed that strain RW2, GIC31, and N139 belong to same bacterial species-the recently described Exiguobacterium chiriqhucha (Gutiérrez-Preciado et al., 2017). We also provide the first complete genomes for this species-represented by our assembly of strain RW2, as well as our the complete genome of strain GIC31. Each genome consists of a single circular chromosome, with strain GIC31 also possessing one plasmid.

Adaptation to a Wide Range of Abiotic Conditions
While strain RW2 has a modest tolerance for added salts, it exhibits both the highest thermal growth range (4-50 • C) and broadest pH growth range (5-11 pH) ever reported for the genus Exiguobacterium. Strain RW2 genome encodes molecular chaperons (csp and dnaK), ABC transporters, antiporters, for thermal, freezing and osmotic stress adaptation. This is unusual, a few of these extremes would be expected in the permanently cold and alkaline habitat from which this isolate was obtained. Growth at such high temperatures and low salinities is not a requirement for life in Pavilion Lake. Though there are numerous hot springs in southwestern British Columbia (Lim et al., 2009), and geothermally heated microhabitats have been known to provide refugia for certain flora during glacial periods Alternately, if strain RW2 is a recent colonist of microbialite habitats, then these physiological tolerances might be holdovers from a more generalist lifestyle. But speculation aside, in the present environment of our isolate, cold and alkaline conditions are the most likely sources of abiotic stress.
Its thermotolerance is unique since, until now, Exiguobacterium spp., strains that grow at >45 • C have only been isolated from such hyperthermal environments (Crapart et al., 2007). Conversely, the percentage of unsaturated fatty acids decreased at higher temperatures. Indeed, the largest change in the PLFA composition of strain RW2 was prompted by increasing the temperature to 50 • C, as branched monoenoic (unsaturated) PLFAs decreased as much as 93% between 4 and 50 • C. This changing membrane profile is a characteristic cellular response to heat stress, as high levels of saturated fats decrease membrane fluidity and thereby prevent phase transitions at high temperature (Hazel, 1995). As with its lipid profile, the genome of strain RW2 seems to reflect its broad thermal and pH tolerance. It possesses a full heat-shock protein cascade (i.e., dnaJ, dnaK, and GrpE), as does its hot-spring dwelling relative, strain AT1b.
Strain RW2 considerable tolerance of acidity is surprising, as similar acidic pH tolerance has only been reported from marine isolates of Exiguobacterium (i.e., E. profundum, E. aestuarii, and E. marinum) which can grow at a pH of 5.5 (but not the upper limit of pH 11 as in strain RW2) (Kim et al., 2005;Crapart et al., 2007). Our genomic annotations indicate that strain RW2 possesses the antiporter encoded by arcD (Fulde et al., 2014), and a Na + /K + antiporter encoded by nhaC, whichin other members of Exiguobacterium-are known to regulate responses to acidity and alkalinity, respectively (Ito et al., 1997). Under decreasing pH, PLFAs increased in saturation, presumably resulting in a more ordered and rigid membrane. Such increases in saturation are thought to help resist H + influxes under acidic conditions (Mykytczuk et al., 2013). Compared to its responses to pH and temperature, PLFA composition in strain RW2 was less affected by salinity, though the percent of saturated PLFAs increased slightly at low salinities as would be expected to increase membrane rigidity and reduce the risk of lysis (Ponder et al., 2005). It is also conceivable that strain RW2 could adapt to high pH via cardiolipin biosynthesis since alkaliphilic strains of Bacillus are known to upregulate this pathway under high pH (Clejan et al., 1986). Further experimental evidence is required to understand how strain RW2 regulates its responses to changes in pH.
Strain RW2 can grow from 0 to 7% added NaCl, and could potentially respond to regulate low salt stress higher osmotic stress by increasing intracellular NaCl concentration using via the choline and betaine uptake and biosynthesis pathways encoded in its genome. Halotolerance of strain RW2 is considerably lower than marine isolates found in the clade IIa, such as Exiguobacterium aestuarii, which can grow at 17% added NaCl (Kim et al., 2005). Genome prediction suggest strain RW2 potentially regulates NaCl stress via choline and betaine uptake and biosynthesis pathway. ABC transporter genes required for higher salinity growth were present in the genome of strain RW2 (i.e., opuC, proW). However, strain RW2 was missing the bet genes thought to be required for high salinity growth (>10%) (Mykytczuk et al., 2013). P. halocryophilus has bet genes required for high salinity growth and can grow in 18% added NaCl (Mykytczuk et al., 2013). The heat shock dnaK gene cluster also can also protect cells from hyperosmotic shock induced by salt (i.e., NaCl/KCl) and carbon substrates (e.g., sugars) (Bianchi and Baneyx, 1999). Overall, the temperature, pH, and salinity-based changes in its PLFA composition seem to reflect strain RW2 adaptability likely due to membrane restructuring.

CONCLUSION
Exiguobacterium chiriqhucha strain RW2 was isolated from microbialite cyanobacterial mat within a Pavilion Lake in southeastern BC. Currently, strain RW2 has the widest pH and temperature growth range for any isolate of Exiguobacterium. We described the physiological, genomic, and fatty acid composition under different growth conditions for strain RW2 widely adaptability to various temperatures, pHs, and salinities. We revealed the putative pathway for the orange colony formation in Exiguobacterium as 4-4' diapocarotenoids C 30 carotenoidbased biosynthetic pathway using comparative genomic analysis. We suggest the potential role of these carotenoid pigments as protective potentially against photodamage or act as quenchers for reactive oxygen.
Our results suggest that strain RW2 is present within microbialite microbiome of Pavilion Lake. Metagenomic recruitment suggest genes related to the detoxification of heavy metals and antibiotics are supplied by strain RW2 genome. Strain RW2 is a robust biofilm former, and putatively increases alkalinity leading to carbonate precipitation is deamination of amino acids to ammonia (i.e., L-asparaginase/urease).
The remarkable range of possible growth conditions and metabolic potential of strain RW2 may make this species relevant to industries that are exploring the use of Exiguobacterium spp. For example, its predicted heavy-metal tolerance makes it an excellent candidate for bioremediation applications (Park and Chon, 2016), and its ability to produce carotenoids may have applications in natural pigment biosynthesis. Strain RW2 is also highly amenable to cultivation because it can grow on a diverse set of nutrients. Lastly, the ability to catabolize urea despite the apparent lack of genes encoding urease may lead to the discovery of a hitherto unknown protein or enzymatic pathway that could function in the transformation of nitrogen.

AUTHOR CONTRIBUTIONS
RW designed the study, and collected and plated the isolate, performed growth studies, extracted DNA, prepared libraries, assembled and annotated the genome, performed comparative genomic and phylogenetic analysis. SS completed PLFA, with financial support from GS. EG performed culturing experiments. RW preserved cells for scanning electron microscopy, which were imaged by GG. RW and GG wrote the manuscript. All authors participated in the manuscript drafting process.

FUNDING
Financial support for sample collection was provided by the MARSLIFE Project (9F052-10-0176) funded by the Canadian Space Agency. The laboratory work was supported by a Discovery Grant from the Natural Science and Engineering Council of Canada (CAS), and grants to CS from the Tula Foundation and the Canadian Institute for Advanced Research. Infrastructure support for laboratory work was provided to CS from the Canadian Foundation for Innovation and the British Columbia Knowledge Development Fund.

ACKNOWLEDGMENTS
We are grateful to the entire Pavilion Lake Research Project (PLRP) research team, especially Donnie Reid for his exceptional leadership of field logistics and diving operations, as well as the PLRP science divers and the DeepWorker pilots for their help recovering samples from the lake. We also thank Danielle Winget, Kynan Suttle, Jan Finke, Marli Vlok, and Tyler Nelson for assistance with sample processing in the field. Microbialite morphology photographs in Figure 1 are used with permission and copyrighted by Donnie Reid, Tyler Mackey, Amy M. Chan, Kynan Suttle, and RA. We thank the following individuals and their teams for providing the sequence data: Sugandha Dandekar (Uma) and Hemani Wijesuriya for the MiSeq data (UCLA Genotyping and Sequencing Core), Frederick Robidoux for the HiSeq data (McGill University and Genome Quebec Innovation Centre), and Sergio Pereira for the 454 FLX data (The Centre for Applied Genomics, SickKids Hospital, Toronto, ON, Canada). Special thanks to Niels W. Hanson and Kishori M. Konwar for R scripts for parsing metapathways data. Also, thank you, Brendan P. Burns, for your excellent comments on draft manuscripts. We are also grateful to Linda and Mickey Macri for hosting the PLRP and MARSLIFE projects and to the Ts'Kw'aylaxw First Nation and British Columbia Parks for their continued support of our research.