“Candidatus Thermonerobacter thiotrophicus,” A Non-phototrophic Member of the Bacteroidetes/Chlorobi With Dissimilatory Sulfur Metabolism in Hot Spring Mat Communities

In this study we present evidence for a novel, thermophilic bacterium with dissimilatory sulfur metabolism, tentatively named “Candidatus Thermonerobacter thiotrophicus,” which is affiliated with the Bacteroides/Ignavibacteria/Chlorobi and which we predict to be a sulfate reducer. Dissimilatory sulfate reduction (DSR) is an important and ancient metabolic process for energy conservation with global importance for geochemical sulfur and carbon cycling. Characterized sulfate-reducing microorganisms (SRM) are found in a limited number of bacterial and archaeal phyla. However, based on highly diverse environmental dsrAB sequences, a variety of uncultivated and unidentified SRM must exist. The recent development of high-throughput sequencing methods allows the phylogenetic identification of some of these uncultured SRM. In this study, we identified a novel putative SRM inhabiting hot spring microbial mats that is a member of the OPB56 clade (“Ca. Kapabacteria”) within the Bacteroidetes/Chlorobi superphylum. Partial genomes for this new organism were retrieved from metagenomes from three different hot springs in Yellowstone National Park, United States, and Japan. Supporting the prediction of a sulfate-reducing metabolism for this organism during period of anoxia, diel metatranscriptomic analyses indicate highest relative transcript levels in situ for all DSR-related genes at night. The presence of terminal oxidases, which are transcribed during the day, further suggests that these organisms might also perform aerobic respiration. The relative phylogenetic proximity to the sulfur-oxidizing, chlorophototrophic Chlorobi further raises new questions about the evolution of dissimilatory sulfur metabolism.


INTRODUCTION
Microbial sulfate respiration is an ancient metabolic process for energy conservation, which may have originated as early as 3.47 billion years ago (Shen et al., 2001;Wacey et al., 2011). Sulfate-reducing microorganisms (SRM) have global importance for biogeochemical sulfur and carbon cycling and are prevalent in both marine sediment (Bowles et al., 2014) and freshwater wetland ecosystems (Pester et al., 2012). It has been estimated that more than 50% of the organic carbon in marine sediments is mineralized via sulfate reduction (Jørgensen and Fenchel, 1974;Jørgensen, 1977Jørgensen, , 1982. Despite its suggested antiquity, dissimilatory sulfate reduction (DSR) is patchily distributed among Archaea and Bacteria and has been experimentally observed solely within isolated members of four bacterial (Deltaproteobacteria, Nitrospirae, Firmicutes, and Thermodesulfobacteria) and two archaeal (Euryarchaeota and Crenarchaeota) lineages (Rabus et al., 2013). However, dsrAB genes have recently been detected in metagenomic assembled genomes (MAGs) affiliated with members of the phyla Acidobacteria, Chloroflexi, and several other phylum-level taxa, suggesting the capacity to perform DSR by some members of these phyla (Anantharaman et al., 2016Hausmann et al., 2018).
In all recognized SRM, the reduction of sulfate to sulfide is mediated by three enzymes. In the first step of the process, ATP sulfurylase (Sat) activates chemically refractory sulfate to adenosine-5 -phosphosulfate (APS). The second enzyme, APS reductase (Apr, encoded by aprBA), reductively cleaves APS to AMP and bisulfite. The third enzyme, dissimilatory (bi)sulfite reductase (Dsr), reduces bisulfite to a protein-based trisulfide, in which the bisulfite-derived sulfur is bound to two cysteines of DsrC, the cellular levels of which can be used for determination of physiological sulfate reduction rates (Santos et al., 2015). The Dsr-encoding genes dsrAB are often used as the diagnostic enzyme for sulfate reduction in environmental studies (Rabus et al., 2013). Sulfate reduction is coupled with the oxidation of menaquinol in the membrane, which leads to the generation of a proton gradient via the membrane-bound redox complexes. Qmo (quinone-interacting membrane-bound oxidoreductase, encoded by qmoABC) transfers electrons to Apr, and the membrane-bound enzyme complex DsrMKJOP, which is homologous to HmeABCDE [heterodisulfide reductase (Hdr)like, menaquinol-oxidizing enzyme], transfers electrons to Dsr during reduction of (bi)sulfite to sulfide (Santos et al., 2015;Dahl, 2017).
Proteins homologous to enzymes of the DSR pathway are also present and highly conserved in anoxygenic photolithotrophic and chemolithotrophic sulfur-oxidizing microorganisms (SOM), including, the strictly anaerobic, chlorophototrophic green sulfur bacteria (GSB) from the phylum Chlorobi (Dahl, 2017). Because these genes are used in both types of dissimilatory sulfur metabolism, the oxidative or reductive directionality of the pathway cannot easily be deduced from the presence of these genes alone. However, phylogenies for Dsr-encoding dsrAB and 16S rRNA genes are largely congruent with clear and distinct clades for SRM and SOM. Observed exceptions are mostly indicative of horizontal gene transfer (HGT) of dsrAB among major SRM taxa. Thus, the dsrAB genes are suitable diagnostic, phylogenetic, and functional marker genes (Klein et al., 2001;Zverlov et al., 2005;Loy et al., 2009;Müller et al., 2015;Pelikan et al., 2016;Hausmann et al., 2018), and have frequently been used as phylogenetic marker genes in amplicon-based sequencing in environmental studies (Dhillon et al., 2003;Nakagawa et al., 2004;Leloup et al., 2006;Dillon et al., 2007;Loy et al., 2009;Moreau et al., 2010;Mori et al., 2010;Pester et al., 2010;Lenk et al., 2011). However, it was recently shown that DsrAB phylogeny may be insufficient to distinguish the pathway directionality of sulfate reduction/sulfide oxidation in SRM and SOM, respectively. For example, Desulfurivibrio alkaliphilus, a deltaproteobacterium, is a SOM, but its DsrAB sequences are phylogenetically closer to those of SRM (Thorup et al., 2017).
The phylum Chlorobi has long been considered to be synonymous with its phototrophic, sulfur-oxidizing members, the GSB. The recent isolation and characterization of two chemoheterotrophic relatives, Ignavibacterium album and Melioribacter roseus, led to the restriction of the GSB to the family Chlorobiaceae within the class Chlorobea, one of six lineages within the Chlorobi (Iino et al., 2010; also see Liu et al., 2012a,b). The taxonomic level for these lineages has not been conclusively resolved. They were initially introduced as class-level lineages (Iino et al., 2010), but phylum level lineages have subsequently been suggested for some of them (e.g., Ignavibacteriae and "Ca. Kapabacteria, " but also Chlorobaeota, Podosokorskaya et al., 2013;Kantor et al., 2015;Oren et al., 2015). Sulfur oxidation in GSB relies on DSR genes similar to the ones found in SRM and other SOB. However, with regard to their DSR genes, GSB hold a special position and seem to contain a chimeric DSR system, in which some genes are most closely related to those of SOM, while others resemble those of SRM. It thus has been proposed that GSB have acquired either their whole or partial DSR system horizontally, although it is not clear whether the genes came from a SOM, a SRM, or both (Sander et al., 2006).
Microbial mats in the effluent channels of Octopus Spring and Mushroom Spring (MS), which are geochemically similar alkaline, siliceous hot springs in the Lower Geyser Basin of Yellowstone National Park (YNP), United States, are among the most intensively studied natural microbial communities . These mat communities have been studied for over 50 years and serve as model systems for exploring principles of microbial ecology (Brock, 1967;Ward et al., 1998;Ward et al., 2008). Most studies have focused on the numerous chlorophototrophic bacterial populations that occur in these mats (including Cyanobacteria, Chloroflexi, and the newly discovered species Chloracidobacterium thermophilum and "Candidatus Thermochlorobacter aerophilum" (Bryant et al., 2007;Liu et al., 2012b;Tank et al., 2017). Next-generation sequencing and metagenomic analyses were recently applied to characterize chemotrophic members of the community (Thiel et al., 2016. Although sulfate levels in the MS microbial mat and the spring water are low (<200 µM), the mat community sustains a highly active sulfur cycle Dillon et al., 2007;Thiel et al., 2017). Previous cloning experiments targeting the dsrAB genes revealed four putative SRM in these mats: a Thermodesulfovibrio sp. (Nitrospira) and three additional unidentified phylotypes (Dillon et al., 2007).
Here, we describe a putative SRM that was initially discovered in MS and that is also found in other hot springs. We describe three partial genomes associated with non-chlorophototrophic organisms affiliated with "Chlorobi-lineage 5" (also known as "OPB56 clade" and "Ca. Kapabacteria") (Kantor et al., 2015;Hiras et al., 2016). These organisms were identified by metagenomic analysis of hot spring microbial mat communities, and they represent the first SRM (or the first sulfur-oxidizing chemolithotroph) from the Bacteroidetes/Chlorobi group. The presence and diel transcription patterns of the genes strongly argue that the corresponding gene products are functional in these organisms. These findings raise important questions concerning the evolution of dissimilatory sulfur metabolism, the timing of the acquisition of these dsr genes within ancestral members of the Chlorobi, and the physiological flexibility of these novel organisms. Finally, these newly discovered, sulfurmetabolizing members of the Bacteroidetes/Chlorobi group may help to clarify the evolutionary origin(s) of dissimilatory sulfur oxidation within the GSB.
Mushroom Spring samples were collected from an oxygenic phototrophic microbial mat in one of the effluent channels of the slightly alkaline, siliceous spring at 60 • C water temperature (GPS position 44.539/−110.798; Figure 1A) as previously described (Thiel et al., 2016. Samples for metatranscriptome analysis were taken at 19 different time points over the course of 1 day in 2014 (October 6, 13:00 to October 7, 17:00). Samples were taken every 1-3 h, and four biological replicates were collected. The microbial mat consists of a 1-2 mm upper green layer, which includes many types of chlorophototrophic bacteria (Tank et al., 2017) and an orange-colored undermat layer that is metabolically active to a depth of approximately 5-7 mm below the surface. Metagenomic, as well as partial metatranscriptomic analyses, of the chlorophototrophic community of the upper green layer have previously been published (Klatt et al., , 2013Liu et al., 2011;Thiel et al., 2016Thiel et al., , 2017. In this study, the two layers were divided and separately analyzed. Samples from Nakabusa hot springs were collected from an oxygenic, chlorophototrophic microbial mat from a small side pool that is fed by a stream of slightly alkaline spring water with a temperature of ∼55 • C ( Figure 1B). Triplicate samples, taken from the so-called "Stream Site" = "Site B" (36.393/−137.748, Nishihara et al., 2018), were collected from the region illustrated in Figure 1B over the time course of sampling. Samples for DNA extraction were collected in triplicate at the beginning and the end of a 39-h diel metatranscriptome experiment on November 3rd and 4th, 2016 using a #4 cork borer (∅ = 8 mm). Samples for RNA extraction were collected hourly during the diel experiment in 2016 (November 3, 5:00 to November 6, 19:00, three biological replicates). The uppermost 3-mm mat layer, consisting of approximately equal parts of green and orange color, was separated from the rest using a knife, and the samples were quickly frozen in a dry-ice cooled, 70% (v/v) ethanol bath in 2ml screw cap tubes. The samples were kept frozen on dry ice and were then transferred to a −80 • C freezer until used. RNA from 20 of the 39 time points was sequenced for the metatranscriptomic analysis.
Iron-silica oxide samples were collected from two locations at Chocolate Pots (GPS position 44.710/−110.741; Figure 1C), a Fe-rich circumneutral-pH thermal spring, in September 2012 and used to initiate two microbial Fe(III)-reducing enrichment cultures, hereafter referred to as "vent" and "midway." Temperature and pH at the two sample locations were 50.7 • C, pH 5.9, and 42.8 • C, pH 7.2 for the hot spring vent and midway down the flow path, respectively (Fortney et al., 2016). Fe(III)-reducing enrichment cultures were incubated in artificial spring water media following the geochemical measurements made by Parenteau and Cady (2010), at pH 6.1, 45 • C, and pH 6.7, 40 • C for the vent and midway enrichment cultures, respectively. Electron donor in the form of 2 mM acetate and 2 mM lactate was added to the media, and 10% (v/v) transfer to new media was made every 2 weeks, hereafter referred to as a generation (for details see Fortney et al., 2016).

DNA Extraction
Genomic DNA was isolated from five replicate mat samples from MS using an enzymatic cell lysis protocol followed by phenol extraction as previously described (Thiel et al., 2016). For samples from Nakabusa hot springs, DNA was extracted and pooled from six replicate samples (three taken before and three after the diel experiment) using the MoBio PowerBiofilm DNA extraction kit (QIAGEN Inc., United States), following the manufacturer's protocol. All centrifugation procedures were performed at 13,000 × g at room temperature. DNA concentration and purity was tested using dsDNA Broadrange (BR) assay on a Qubit fluorometer (Life Technologies, Grand Island, NY, United States) and Bio-Spec Nano (UV-VIS Spectrophotometer, Shimadzu, Japan).
Genomic DNA was extracted from generation 33 of the Chocolate Pots vent and midway Fe(III)-reducing enrichment cultures using the MoBio PowerSoil DNA isolation kit as previously described (Fortney et al., 2016).

RNA Extraction
RNA from MS microbial mat samples was extracted using a hot-phenol extraction protocol adapted from Steunou et al. (2006) and Aiba et al. (1981). In brief, frozen mat samples were homogenized for 10 s using 0.5 g 150-200 µm glass beads in acidic acetate-EDTA buffer (250 µl 10 mM Na-acetate, pH 4.5; 37.5 µl 500 mM Na 2 -EDTA, pH 8.0). After the addition of 375 µl of lysis buffer [16 mM Na-acetate, 2% (w/v) SDS, pH 4.5], the sample was vortexed and incubated at 65 • C for 3 min. Another 3-min incubation step followed at 65 • C in hot acidic phenol, pH 4.5 (700 µl), before the samples were quickly cooled on ice prior to the initiation of phase separation by centrifugation (2 min at 17,000 × g). The aqueous phase was collected, and nucleic acids were purified with phenol:chloroform (1:1 v/v) and chloroform steps; the RNA was then precipitated using 0.1 volume of 10 M LiCl + 2.5 volumes 100% ethanol at −20 • C for 30 min. Precipitated RNA was collected by centrifugation at 4 • C (17,000 × g, 60 min), washed with 80% (v/v) ethanol, dried in a vacuum concentrator at room temperature, and resuspended in RNase-free water. Extracted RNA was treated with DNase following the Joint Genome Institute (JGI) recommendations with DNase I (2 U per 10 µg RNA; Ambion AM222) and then purified and concentrated using the QIAGEN RNeasy R MinElute R kit.
RNA was extracted from mat samples (0.10-0.21 g wet weight) from Nakabusa hot spring after centrifugation for 1 min at 13,000 × g at room temperature to remove excess liquid. These samples were subjected to RNA extraction using the RNeasy PowerBiofilm Kit (QIAGEN, Inc., United States) following the protocol of the manufacturer. The RNA was treated with DNase I, eluted in RNAse-free water (50 µl), and its concentration and purity were evaluated using the RNA High Sensitivity (HS) assay with Qubit fluorometer (Life Technologies, Grand Island, NY, United States) and Bio-Spec Nano (UV-VIS Spectrophotometer, Shimadzu, Japan).

Sequencing
Genomic DNA extracted from MS mat samples was sequenced at the DOE JGI, United States, using Illumina HiSeq technology. Metagenome and 16S rRNA gene amplicon sequencing was conducted on Illumina HiSeq and MiSeq instruments, respectively, as described previously (Thiel et al., 2016. Genomic DNA extracted from Nakabusa hot spring microbial mat samples was sequenced using an Illumina MiSeq platform using a single run to obtain 2 × 300 paired-end reads at Fasmac (Atsugi, Japan). Purified DNA (300 ng) was used for library construction using the KAPA HyperPlus Library Prep Kit for Illumina (KAPA Biosystems) according to the protocol of the manufacturer.
Genomic DNA from the Chocolate Pots enrichment cultures was sequenced at the University of Wisconsin Biotechnology Center (UWBC 1 ) using paired-end 2 × 250 bp Illumina MiSeq shotgun metagenomic sequencing technology.
Extracted total RNA samples from both microbial mats (MS and Nakabusa) were treated to remove rRNA using the TruSeq Ribo-Zero Bacteria kit (Illumina). The mRNA-enriched RNA was sequenced at JGI (MS samples) or DNALink Inc. (Seoul, South Korea, Nakabusa hot springs) using Illumina HiSeq 2000 and Next Seq500 instruments, respectively.

Sequence Data Analyses
Metagenomic DNA sequence data from MS were analyzed as previously described (Thiel et al., 2016. Metagenomic data obtained from Nakabusa hot spring samples were analyzed as follows. Raw reads were trimmed, low-quality reads were removed using sickle tools ver. 1.33, and remaining reads were assembled using the SPAdes Genome Assembler ver. 3.7.1. after error correction using "BayesHammer" (Nikolenko et al., 2013).
RNA samples isolated from MS collected over a diel cycle were sequenced on an Illumina HiSeq instrument at the JGI, United States. The raw data were pre-processed, quality checked, and trimmed as described for metagenome sequencing reads (Thiel et al., 2016). For a selection of 20 time-point samples from the Nakabusa diel experiment, RNA was sequenced using an Illumina Next Seq500 instrument. The raw RNA reads obtained were pre-processed using FastQC 3 . Adapter and low-quality reads were trimmed by Cutadapt version 1.12 (Martin, 2011). Quality-checked reads were mapped to the assembled reference metagenome using bowtie2 version 2.3.0 (Langmead and Salzberg, 2012) using the default settings allowing no mismatches, and aligned using the EDGE-pro algorithm (Magoc et al., 2013), using the rRNA depletion option. Differential expression analysis was conducted with DESeq2 (Love et al., 2014 4 ).

Binning, Annotation, and Phylogenetic Affiliations Using Marker Genes
Metagenomic contigs of the metagenomes obtained from MS and Nakabusa were binned on the basis of tetranucleotide frequency patterns, and the bins were visualized using emergent self-organizing maps (ESOM) as described . Metagenomic contigs from the Chocolate Pots vent and midway metagenomic assemblies were binned using CONCOCT (Alneberg et al., 2014), and the taxonomic identities of the vent and midway MAGs were determined using PhyloSift (Darling et al., 2014). Unbinned vent and midway metagenomic assemblies were uploaded to Integrated Microbial Genomes (IMG) with Microbiome Expert Reviewer (IMG/M ER 5 ) for gene annotation using the default IMG gene-calling method (Mavromatis et al., 2009). The scaffolds obtained from assembly and binning of the metagenomic sequences were treated as a single genome and automatically annotated using the RAST annotation server (Rapid Annotation using SEED Technology) (Aziz et al., 2008;Overbeek et al., 2014). Completeness and purity of MAGs were assessed using CheckM (Parks et al., 2015).
Phylogenetic affiliations were evaluated using the AMPHORA2 metagenomic workflow suite, which uses 31 universal bacterial phylogenetic marker genes for taxonomic affiliation of a genome sequence. This suite also uses Phyla_ AMPHORA, which searches (partial) genomes for phylum-level, specific phylogenetic marker genes of 20 different bacterial phyla (Wu and Scott, 2012;Wang and Wu, 2013).
A set of 1544 representative, complete prokaryotic genomes was downloaded from the NCBI database 6 [Filter settings "prokaryotes" (organism group), "Representative" (RefSeq category), "Complete genome" (Assembly level); accessed on March 15, 2018]. Phylogenetic marker genes were extracted from the genome sequences using AMPHORA2 (Wu and Scott, 2012). Genomes were filtered by the presence and number of phylogenetic marker genes detected, and a selection of 1152 representative prokaryotic genomes containing one copy of each marker gene per genome was used for further analysis. Additional partial genomes closely related to the metagenomic bins analyzed in this study were selected for analyses based on BLAST results of phylogenetic marker gene sequences. A total of 1171 (partial) genomes were used for phylogenetic analysis of phylogenetic marker genes.
Phylogenetic marker genes were extracted using AMPHORA2 from the (partial) genomes, and each set of extracted amino acid sequences was aligned using MEGA7 (Kumar et al., 2016) (Linux command line version M7CC, implemented MUSCLE align, data type "Protein, " standard settings except for "Gap Open" penalty value of −2.00 and "Gap Extend" penalty of −0.50, output format "FASTA"). The aligned sequences of 30 of the 31 marker genes (smpB was omitted due to low coverage) were concatenated using "Java Sequence Matrix" 7 . The output Nexus file was converted into FASTA format using AliView (Larsson, 2014). Alignments were imported into ARB (Ludwig et al., 2004) and manually refined. The amino acid alignment was filtered to exclude unalignable, highly variable regions within and between genes and was exported for phylogenetic analysis conducted using MEGA7. The maximum-likelihood phylogenetic analysis was based on the WAG model, with Gamma distributed rates among sites (five categories) using all sites selection. All nonmentioned parameters were used in default setting. Robustness of the trees was tested using non-parametric bootstrap analysis with 100 replicates in all cases.
Phylogenetic affiliations of MAGs and reference genomes were tested based on 16S rRNA gene sequences, when available. Phylogenetic trees were conducted based on the maximumlikelihood method using the "PhyML (DNA)" tool implemented in the software package ARB 8 (Ludwig et al., 2004). Alignments of sequences obtained in this study as well as from database references were manually refined based on the SILVA/ARB SSU RefNR99 database Release 123 9 . GTR substitution model, four substitution rate categories, ML-based frequency estimates, as well as estimated ts/tv ration, gamma distribution parameter, and proportion of invariable sites settings were used. Robustness of branching was tested using 100 bootstrap replicates.

Phylogenetic Analysis of dsr Genes
Based on their amino acid sequences, all genes related to DSR or oxidation were analyzed using Basic Local Alignment Tool (BLAST; Altschul et al., 1990) and conserved domain database (CDD) (Marchler-Bauer et al., 2015) available through the National Center for Biotechnology Information (NCBI), Bethesda, MD, United States. Phylogenetic analyses were conducted on concatenated sequences for DsrAB. Reference protein sequences for phylogenetic analyses were obtained by BLAST searches against a selected number of representative genome sequences in the JGI IMG system (Markowitz et al., 2014) and by BLAST searches against the NCBI nr database. The dsrAB/DsrAB genes/proteins were analyzed using the updated publicly available dsrAB/DsrAB database (Müller et al., 2015) and the ARB software package (Ludwig et al., 2004). Phylogenetic positions were determined by using the parsimony quick-add tool to add the DsrAB sequence to an existing phylogenetic tree. Using a number of reference sequences, a full maximumlikelihood tree was calculated using the PHYML (Amino Acids) tool implemented in ARB. Substitution model WAG, estimated gamma distribution, four substitution-rate categories, as well as estimated proportion of invariable sites and 100 bootstrap replicates were used.

Transcriptomic Analyses
Using a proprietary mapping pipeline 10 , reads from the MS metatranscriptome were mapped onto the corresponding upper green mat layer and undermat metagenomes (IDs 3300010182 and 3300002493, respectively). Reads from the Nakabusa metatranscriptome were mapped onto the assembled Nakabusa metagenome using bowtie2 version 2.3.0, the default setting bowtie2 adding "-un-conc" option to collect "unmapped" reads. The read counts corresponding to mapped genes in the metagenome, which had been identified as belonging to bin MS-B_bin-24  (MS), and bin Naka2016_bin-10_OTU-51 (Nakabusa), were extracted with user-generated scripts and normalized for each time point as described (Liu et al., 2011(Liu et al., , 2012b. Briefly, read counts were first normalized for each time point by the total number of read counts retrieved for the target organism at that time point. Secondly, the relative expression of each gene during the diel cycle was calculated by normalizing against the mean of all the reads at each time point for that particular gene. This method allows one to compare the relative transcript abundance levels, rather than absolute values, for each gene across the diel cycle. Genes were clustered according to their diel gene expression patterns using the k-means algorithm in Cluster 3.0 (Eisen et al., 1998), and results were visualized with Java Treeview (Saldanha, 2004) as described (Liu et al., 2012b). Normalized expression levels for each gene were imported into Cluster, adjusted by a log transformation, centered by mean, and then clustered using the k-means algorithm with k = 6, runs = 1000, and using other default parameters.

Metagenomic Data
Metagenomic sequencing and assembly of DNA from the undermat of MS (Figure 1) produced a 232-Mb metagenome comprising 315,154 total contigs . The metagenome had a mean GC value of 54%, a maximum scaffold length of 158 kb, and a N/L50 value of 32,529/1.24 kb, which defines the number of fragments at or above the Length50 cutoff. There were 13,766 contigs >2.5 kb, which led to the recovery of 44 metagenomic bins (MAGs). Bin MS-B_bin-24 (Thiel et al., 2016 contained 2.6-Mb in a total of 109 contigs with a mean GC value of 38.5%. It contained a partial 23S rRNA gene, as well as 46 tRNA and 53 ribosomal protein genes ( Table 1). Gene IDs for this MAG refer to the scaffold and gene IDs of the whole metagenome assembly, which is available at JGI IMG/MER (Taxon ID 3300002493). For convenience, the prefix of "JGI24185J35167_1" for each of the scaffold/gene IDs are omitted. Metagenomic sequencing of DNA from a hot-spring microbial mat from Nakabusa, Japan, resulted in 32,921,964 paired-end reads, which were assembled into a 320-Mb metagenome comprising 192,524 total contigs with a mean GC value of 54%, a maximum contig length of 1,260 kb, and a N50 value of 3,718 bp. Binning of 16,846 contigs ≥2.5 kb led to 76 MAGs; Naka2016_bin-10_OTU-51 contained 2.62 Mb in 32 contigs and importantly contained a complete rRNA operon ( Table 1, and see below). Metagenomic sequencing and assembly of the DNA extracted from the Chocolate Pots vent and midway Fe(III)-reducing enrichment cultures yielded metagenomes of 42.1 Mb and 35,690 contigs, and 86.0 Mb and 93,531 contigs, respectively. The vent enrichment culture metagenomic assembly had a mean GC content of 54.1%, a maximum scaffold length of 657 kb, and an N50 value of 2,072. The midway metagenomic assembly had a mean GC content of 52.7%, a maximum scaffold length of 744 kb, and an N50 value of 1,149. The vent and midway enrichment culture metagenomic assemblies are available at JGI/MER under Taxon Object IDs 3300002900 and 3300002895, respectively. A total of 1,520 contigs >3 kb were binned into eight MAGs with an average read coverage of 16.2× in the vent enrichment culture metagenomic assembly. 2,148 contigs >3 kb were binned into 15 MAGs with an average read coverage of Relative completeness and contamination were assessed using CheckM (Parks et al., 2015).
FIGURE 2 | Taxonomic affiliations of metagenomic bins based on phylogenetic marker genes of 20 different bacterial phyla (Phyla_AMPHORA; Wang and Wu, 2013). The majority of phylogenetic marker genes detected in the partial genomes belonged to the taxa Chlorobi (100-167 of 808 phylum-specific marker genes) or Bacteroidetes (55-65 of 215 phylum-specific marker genes) (Wang and Wu, 2013). Detailed number of genes detected for each phylum can be found in Table 2. 13.4× in the midway assembly. A MAG "CP-midway-bin-8, " which was assembled from metagenomic sequences from the midway enrichment culture from Chocolate Pots thermal feature, YNP, had a total sequence length of 2.85 Mb, 33.8% GC content, and comprised 359 contigs with an average read coverage of 7.5×. Statistics for these three MAGs are provided in Table 1. Each of the three MAGs described above contained dsrAB genes, and these metagenomic bins representing putatively novel SRM/SOM organisms were analyzed further.

Phylogenetic Analysis of Marker Genes
Phylogenetic marker-gene analysis by Amphoranet and Phyla_Amphora indicated that the three MAGs were taxonomically affiliated with the Bacteroidetes/Chlorobi group (Figures 2, 3 and Table 2). About half of the phylum-level, phylogenetic marker genes in each MAG were affiliated with Chlorobi, and a further 20% of the marker genes were affiliated with Bacteroidetes. The remaining marker genes were distributed among several phyla, with Firmicutes and Aquificae being the next most abundant (Figure 2). Specifically, the individual MAGs contained 100-167 of 808 phylum-level marker genes specific for Chlorobi; furthermore, the bins contained 55-65 of 215 Bacteroidetes-specific marker genes (Figure 2). By comparison, the genomes of M. roseus and I. album contain 184 and 185 Chlorobi-specific, and 73 and 76 Bacteroidetesspecific marker genes, respectively. A phylogenetic analysis of the concatenated amino acid sequences for 30 marker-gene products further supported the affiliation of the MAGs with the Bacteroidetes/Chlorobi superphylum (Figure 3).

Phylogenetic Affiliation Based on 16S rRNA
Metagenomic bin "Naka2016_bin-10_OTU-51" from Nakabusa hot spring contained a complete rRNA operon on contig_660. However, MS-B_bin-24 only contained a partial 23S rRNA gene, and the ChocPot-Midway-bin-8 did not contain any rRNA sequences ( Table 1). The 23S rRNA gene sequences obtained from Mushroom and Nakabusa hot springs share 98% nucleotide identity, indicating a close relationship between the two organisms. They are additionally closely similar to Bacteroides and Chlorobi fosmid clones JFF029_C06 and JFF027_B02 (NCBI acc. nos. AP011722 and AP011715, 86-90% nt ID) as well as to the type strains of I. album and M. roseus (acc. nos. CP003418/CP003557, 85% nt identity). Phylogenetic analysis of the available 16S rRNA gene sequence placed the organisms represented by the hot-spring-derived MAGs into the "Chlorobilineage 5" (Hiras et al., 2016, also known as "OPB56 clade"; it has recently been suggested that the organisms in this group represent a novel phylum with the name "Candidatus Kapabacteria" Kantor et al., 2015) (Figure 4). A closely related 16S rRNA gene sequence (97% nt identity to the sequence from Nakabusa, MS-B_OTU-29, Thiel et al., 2016) was found in the metagenome produced from MS; it was too short to be assigned during binning but is also consistent with this assignment. All closely related reference genomes that cluster with the novel MAGs in the phylogenetic tree on the basis of the concatenated phylogenetic marker genes in Figure 3 were also placed within "Chlorobi lineage 5"-i.e., the OPB56 cluster on the basis of 16S rRNA sequences (Figure 4). This includes "Ca. Kapabacteria thiocyanatum" as well as organisms represented by MAGs obtained from bioreactors (Kantor et al., 2015), anammox wastewater treatment plants (Speth et al., 2016), compost (Hiras et al., 2016), and a microbial mat from a subsurface geothermal stream (Kato et al., 2018), respectively. Although the ribosomal genes indicate a close relationship between the organisms obtained from MS and Nakabusa, a genome wide comparison indicates an ANI between the two partial genomes of only 75%, which does not support the conclusion that these two organisms are the same species. Therefore, we propose the genus name "Candidatus Thermonerobacter" (Thermonerobacter, thermós, Greek for hot/thermal; nero, Greek for water; bacter, from bakterion, Greek for bacterium) for the novel organisms obtained from MS and Nakabusa hot springs microbial mats, and the species name "Candidatus Thermonerobacter thiotrophicus" (thiotrophicus, thio from Greek theion for sulfur; trophicus from Greek trophus, eater) for the organism inhabiting MS.

Sulfur Metabolism and Clustering of dsr Genes
Complete sets of genes for DSR (or sulfide oxidation) were present in two of the three partial genomes recovered from the hot-spring metagenomes. The slightly less complete metagenomic bin (ChocPot-Midway-bin-8) obtained from the Chocolate Pots enrichment culture was missing 2 (dsrP and qmoA) of the 20 required genes ( Table 3). The genes required for DSR were located on four (MS), two (Nakabusa), or three (Chocolate Pots) different contigs, respectively ( Figure 5). The overall organization of the genes required for sulfide oxidation/sulfate reduction was identical for the three hot-spring MAGs. The genes encoding ATP sulfurylase (sat) and APS reductase (aprBA) occur in an apparent operon with genes encoding the membrane-bound QmoABC complex (qmoABC) in GSB (Frigaard and Bryant, 2008;Bryant et al., 2012). In the MAGs described here, the sat genes were always found adjacent to aprBA, but in contrast to GSB, aprA was upstream of two conserved hypothetical genes before qmoC, which encodes the membrane-bound, heme b-binding subunit of the heterodisulfide-reductase-like "quinone-interacting, membrane-bound oxidoreductase" (Qmo) complex ( Figure 5). The qmoA and qmoB genes were adjacent in two of the three MAGs, but these genes were found in different chromosomal locations in the three MAGs ( Figure 5). Although qmoA was not present in the Chocolate Pots metagenome, it probably occurs upstream of qmoB, which is located at the 5 -end of a contig.
The dsrAB genes are highly conserved components of the Dsr and rDsr (i.e., reversely operating Dsr found in sulfur-oxidizing bacteria) enzymes, which are responsible for the reduction of bisulfite to sulfide or the oxidation of sulfide to bisulfite, respectively (Frigaard and Dahl, 2008;Müller et al., 2015;Santos et al., 2015;Dahl, 2017). In the three hot-spring MAGs, the sat-aprBA-hyp-hyp-qmoC operon was followed by genes encoding the Dsr/rDsr enzyme. The dsrAB genes were followed by dsrD-dsrN-hyp-dsrL-dsrC and homologs of dsrM and dsrK (dsrM-2 and dsrK-2 in Table 3 and Figure 5). The co-occurrence of dsrD and dsrL is unusual, and has only been found in a newly described and uncultured member of the phylum Acidobacteria with dissimilatory sulfur metabolism . Ubiquity of DsrD among SRM suggests an essential function in sulfate reduction, while DsrL usually has been associated with SOMs (Hittel and Voordouw, 2000;Lübbe et al., 2006;Weissgerber et al., 2014;Hausmann et al., 2018). However, since the functions of DsrL and DsrD are still unresolved, functional predictions based only on the presence of these genes are not possible. Genes encoding the presumed menaquinol-oxidizing enzyme complex (dsrTMKJOP), which encode subunits similar to those of membrane-bound heterodisulfide reductases, showed the same organizational pattern in the metagenomes of bins   recovered from the microbial mat communities of Mushroom and Nakabusa hot springs. In both cases dsrT occurs upstream from dsrMKJO, but dsrP is located elsewhere in the genome ( Figure 5). Although not originally annotated by RAST, a partial sequence for dsrO was found downstream from dsrJ on contig_2519 (Figure 5). The dsrP gene is assumed to be present elsewhere in the genome on a contig that was not included in the metagenomic bin.
dsrAB Identity and Phylogeny The dsrAB genes, which were detected in the bin MS-B_bin-24 metagenome ( Table 1) and which encode the alpha and beta subunits of Dsr/rDsr, are most closely related to dsrAB clone sequences obtained from MS in a previous study (Dillon et al., 2007; Table 4 and Figure 6). The dsrAB sequences obtained from the MS metagenome share 98-99% nucleotide and 94-98% amino acid sequence identity to clones MS3.094, MS3.117, and MS3.098 (NCBI acc. nos. EF429279, EF429284, EF429281); these were identified as "clade 2" with unknown phylogenetic affiliation in a dsrAB-targeted, cloning study of the MS microbial mat community (Dillon et al., 2007). Those clones are also the best BLAST hits to dsrAB sequences derived from the Nakabusa hot-spring metagenome (Table 4). Furthermore, phylogenetic analysis of the dsrAB sequences from the three hot-spring MAGs were closely affiliated with other organisms with "unclassified environmental sequences" within the "environmental supercluster 1" of unidentified and uncultured reductive bacterial-(=sulfate-reducing bacteria)-type dsrAB clade (Müller et al., 2015; Figure 6 and Supplementary Figure S1). The sequences found in this study do not cluster with any of the recently published DsrAB sequences found in other MAGs, such as the peat soil members of the Acidobacteria .

Additional Genes Possibly Affiliated With Sulfur Metabolism
In addition to the 20 genes required for DSR (or sulfide oxidation), several other genes related to sulfur metabolism were found in the bins. Three genes encoding a membranebound molybdopterin oxidoreductase of the Psr/Psh family (MS gene IDs 01959_007-009) were present in the MS and Nakabusa MAGs. These genes are most similar to the psrABC genes found in M. roseus [MROS_1774 to 1776] and I. album [IALB_1661 to 1663] (60-78% aa identity). The subunits of this oxidoreductase are similar to the corresponding components of polysulfide and nitrite reductases, and they show 37-52% amino acid sequence identity to polysulfide reductase-like complex 2 (PSRLC2) in GSB (Hinsley and Berks, 2002;Liu et al., 2012a;Kadnikov et al., 2013). A putative transcriptional regulator of the MerR family (MS gene ID 01959_010) precedes these three genes in the MS and Nakabusa MAGs. Transcription patterns (see below) as well as nearby genes located downstream of psrABC support a function in sulfur metabolism: (1) a putative sulfur transporter of the DUF395 family containing two sulfur transport domains (MS gene ID 01959_005/006 and Pfam family Sulf_transp PF04143) and (2) a rhodanese-related sulfurtransferase (MS gene ID 01959_004). Sulfur transporter and sulfurtransferase genes are also present in the Nakabusa MAG, notably on the same contig and only a few genes downstream of the DSR operon. Homologs of PSRLC3, an alternative sulfite-oxidizing system found in GSB, were not found in the metagenomic sequences (Frigaard and Dahl, 2008). Finally, no sulfide:quinone reductases (sqr) were detected in the partial genomes. This is in contrast to GSB but is nevertheless consistent with the idea that these members of the Bacteriodetes/Chlorobi group reduce sulfate rather than oxidize sulfide. Interestingly, the I. album genome encodes two sulfide:quinone reductase genes but lacks the genes for sulfate reduction (Liu et al., 2012a). These uncultured organisms include a complete set of genes for glycolysis and the TCA cycle which suggests that they have the ability to oxidize organic substrates completely to CO 2 . Many SRM can use the Wood-Ljungdahl pathway (Fuchs, 2011), including carbon monoxide dehydrogenase (EC 1.2.7.4)/CO-methylating acetyl CoA synthase (acetyl-CoA synthase, EC 2.3.1.169), to synthesize acetate from CO 2 . However, the genes for this key enzyme were absent from all three MAGs. Likewise, genes for other CO 2 fixation pathways were not present in the MAGs, which implies that these organisms cannot grow autotrophically. Nevertheless, genes for pyruvate:ferredoxin oxidoreductase (PFOR, pyruvate synthase; EC 1.2.7, gene ID 000011_21) and 2-oxoglutarate:ferredoxin oxidoreductase (KFOR EC 1.2.7.3; gene ID 000227_19-20), which can carboxylate acetyl-CoA and succinyl-CoA, respectively, were present in the metagenomes (see description of the section "Amino Acid Metabolism"). Thus, these organisms might incorporate some CO 2 into biomass as apparently occurs for C. thermophilum (Tank and Bryant, 2015a,b). The presence of acetyl-CoA synthetase (gene ID 000018_23, EC 6.2.1.1) further suggests these organisms can probably assimilate acetate. Finally, the presence of genes for phosphate acetyltransferase (phosphotransacetylase, EC 2.3.1.8; gene ID 000091_20) and acetate kinase (EC 2.7.2.1, gene ID 000091_21) suggests that under some conditions, these organisms may use acetyl-phosphate to produce ATP by substrate-level phosphorylation with subsequent excretion of acetate.

Electron Transfer and Respiration
Similar to the genome of I. album (Liu et al., 2012a), the three hot-spring MAGs contain two sets of genes encoding type-I NADH dehydrogenase: one set encodes 14 subunits (nuoABCDEFGHIJKLMN) that are present in different parts of the genome, while the other set lacks the nuoEFG genes and comprises a single 11-gene operon (Figure 7). The absence of nuoC in the metagenomic bin derived from Chocolate Pots is assumed to be due to incomplete genome coverage. Genes closely related to NuoEFG (or to hydrogenase subunits, HydABC) are encoded elsewhere in the genome and appear to encode a soluble, FIGURE 6 | Maximum-likelihood tree based on DsrAB amino acid sequences displaying the phylogenetic affiliation of the novel thermophilic putative SRB Bacteroidetes/Chlorobi group member within the Environmental Supercluster. Sequences obtained in this study are indicated in red bold font. Sequences shown in black bold font were obtained from MAGs with known phylogenetic affiliations from previous studies. The tree was constructed using an alignment in the publicly available DsrAB ARB database (Müller et al., 2015; http://www.microbial-ecology.net/download). Robustness of branching patterns was tested with 100 bootstrap replicates. Scale bar represents 0.1 changes per amino acid position.
bifurcating hydrogenase (Schut and Adams, 2009). Such an enzyme could reversibly oxidize hydrogen and produce NADH and reduced ferredoxin as products. These products could be re-oxidized by the two type-1 NADH dehydrogenase complexes as shown in Figure 7. If the NADH dehydrogenase complex and bifurcating Hyd act together, up to four protons could be translocated per hydrogen molecule oxidized because both NADH and reduced ferredoxin can be used to reduce quinones. As such, this mechanism would be more favorable than using membrane-bound hydrogenases to oxidize hydrogen and directly produce reduced quinones, which would account for only two protons translocated (Buckel and Thauer, 2018a,b).
Three genes encoding succinate dehydrogenase (complex II) as well as five genes encoding alternative complex III (ACIII) are present in each of the four partial genomes. The identified genes for the latter share highest sequence similarities to corresponding act genes of I. album (IALB_1386, 1388, 1389, 1390, and 1391) and M. roseus (MROS_0041-45). Similar to Dokdonia sp. MED134, Schleiferia sp., and other members of the Bacteroidetes, as well as Bacterioidetes/Chlorobi members I. album and M. roseus, genes encoding ubiquinol:cytochrome c oxidoreductase (complex III) were not detected (González et al., 2011;Liu et al., 2012a;Kadnikov et al., 2013;Thiel et al., 2014).
Four genes (MS gene ID 000332_3-6) encoding an aa 3 -type cytochrome c oxidase (complex IV) are present in each of the MAGs, which suggests these putative SRM have the ability to perform aerobic respiration for growth under some conditions or can detoxify and tolerate relatively high oxygen concentrations.
Additionally, genes for a cytochrome bd-type quinol oxidase are present; this type of quinol oxidase typically has a very high affinity for oxygen, and thus it can function efficiently at very low oxygen concentrations (Forte et al., 2017). In contrast to the genomes of I. album, M. roseus, and some GSB strains, genes encoding an alternative, high-affinity cytochrome cbb 3type oxidase were not observed in the four partial genomes. As found for I. album, a complete Na + -translocating Fd:NAD + oxidoreductase (RNF complex) was present in the partial genome from Nakabusa. RnfCDEGAB-encoding genes most closely related to genes in Ignavibacteriaceae and the Nakabusa sequences were also detected in the MS metagenome on several short contigs not included in the binning. This might suggest that the RNF complex also occurs in the organism found in MS, but at the moment this cannot be inferred conclusively (Thiel et al., 2016).
In addition to possible aerobic respiration and DSR, genes encoding other terminal reductases of anaerobic respiratory pathways are present in the analyzed MAGs. Although no genes for nitrate reductase were found in any of the MAGs, all three contain NrfAH-type nitrite reductase (MS gene ID 00001139-40), nitric oxide reductase (CytS, MS gene ID 000001125), and nitrous oxide (NosZ, MS gene ID 00000354) reductase (Figure 7). This is similar to the situation in M. roseus .

Motility
Genes necessary for flagellar motility were present in the metagenomes, but the flagellar genes do not form a single gene cluster as in the I. album genome. Furthermore, a complete set of chemotaxis proteins, CheABWRYZ were found in the MS and Nakabusa MAGs, although only some genes were initially affiliated with MS-B_bin-24 while others were found on contigs too short to be reliably binned. Currently, the identities of potential attractants and/or repellants are unknown.

Amino Acid Metabolism
On the basis of the partial genomic information, a complete pathway for lysine biosynthesis is present; this is similar to Chloroherpeton thalassium but unlike I. album, which lacks this ability. Only incomplete metabolic pathway information for biosynthesis of other amino acids exists at present. However, based on the information currently available in the partial genomes, these organisms probably are unable to synthesize two branched-chain amino acids, valine and isoleucine. In contrast, the ability to synthesize leucine is assumed based on the presence of all genes needed from the intermediate 2-isopropylmalate, although 2-isopropylmalate synthase (EC 2.3.3.13) has not been detected. Notably, they should be able to degrade all branched-chain amino acids like "Ca. T. aerophilum, " and C. thermophilum (Liu et al., 2012b;Tank and Bryant, 2015a,b). Thus, these organisms may primarily synthesize 2-oxoglutarate by carboxylation of succinyl-CoA as has been suggested to occur in C. thermophilum and "Ca. T. aerophilum" (Liu et al., 2012b;Tank and Bryant, 2015a,b;Tank et al., 2017). The currently available information from the three metagenomes suggests that these SRM may be unable to synthesize aromatic amino acids (phenylalanine, tyrosine, tryptophan), histidine, arginine, proline, and possibly cysteine, while serine, glycine, and alanine are synthesized (Figure 7).

Pigments
The partial genomes do not encode any genes for the synthesis of (bacterio)chlorophyll or other proteins required for photosynthesis. The MAGs do, however, encode enzymes for carotenoid biosynthesis, and it is predicted that these organisms should produce cyclic xanthophyll carotenoids. Genes encoding phytoene synthase (crtB, 3 copies), phytoene dehydrogenase (crtI), lycopene cyclase (crtY), β-carotene hydroxylase (crtZ), as well as genes possibly encoding pro-ζ-carotene-producing phytoene desaturase (crtP) and pro-ζ-carotene desaturase (crtQ) were present in the MS MAG. The other two MAGs contained partial gene sets for carotenoid biosynthesis (crtB, crtP in both; crtI only in Chocolate Pots).

Nitrogen Metabolism
No genes for nitrogen fixation (nif ) were found in any of the partial genomes, so it is unlikely that these organisms can reduce dinitrogen to ammonia. The presence of glutamine synthetase (GS) as well as glutamate:2-oxoglutarate amidotransferase (GOGAT or glutamate synthase) could suggest that these organisms use ammonia as an N-source. However, similar to I. album, no genes for ammonium transporters have yet been found in the partial genomes. Thus, the use of free ammonium as an N-source appears unlikely. More likely, the presence of genes for amino acid and oligopeptide transporters in (almost) all partial genomes suggests that these organisms predominantly use amino acids as nitrogen sources. This is a trait that is also shared with several other mat organisms, including C. thermophilum, "Ca. T. aerophilum, " and Synechococcus spp. (Tank and Bryant, 2015a,b;Tank et al., 2017). As mentioned above, and similar to M. roseus , all three MAGs contain genes indicating an ability to reduce nitrite, but not nitrate (CytS, NosZ, and NrfAH; Figure 7). The absence of nitrate reductase in M. roseus was correlated with the observation that nitrate did not stimulate growth, but the organism was able to use nitrite as terminal electron acceptor .

Transcription of Genes in situ
Metatranscriptomic data derived from the upper green mat layer of the MS community over a diel cycle were analyzed to identify the relative transcript abundances for all genes assigned to the MAG "MS-B_bin-24, " representing "Ca. T. thiotrophicus." The diel expression patterns for all genes were then subjected to k-means clustering, which revealed six differential expression pattern categories, A-F (Supplementary Figure S2 and Supplementary Table S1). The highest relative transcript levels for most genes occurred at night (clusters D and E) with over 900 genes (ca. 44%) found in these clusters, and many fewer genes had highest transcript levels during the day (cluster C, see below). The relative transcript abundances for the 20 genes associated with sulfate reduction ( Table 3) were among those genes with the highest absolute transcript levels and belonged to categories D and E. These genes showed highest relative transcript abundances during the night and transition periods when oxygen levels were low (Figure 8 and Supplementary Table S2). A sharp decrease in transcript levels occurred during the day when cyanobacteria were actively performing photosynthesis and when oxygen levels were maximal. A similar pattern was also observed for "Ca. Thermonerobacter" in the mats of Nakabusa; genes associated with DSR showed highest transcript levels during the night and lowest transcript levels during the day (data not shown). Interestingly, in the undermat of MS, at depths where oxygen concentrations never reach the hyperoxic levels of the uppermost part of the mat, transcript levels for the genes associated with sulfur metabolism did not decrease as greatly during the day as in the upper green mat layer (Figure 8, Supplementary  Table S2). The transcript levels for a predicted sulfate transporter (MS gene ID 000020_7) matched the expression pattern for the sulfur metabolism genes (Supplementary Table S1). Other respiration-associated genes in the night expression categories (D and E) included genes involved in additional anaerobic processes, such as the reduction of nitrite, nitrous oxide, and thiosulfate (Supplementary Figure S3A and Supplementary  Table S2).
Relative transcript levels for the two sets of nuo genes were highest at night. Similarly, relative transcript levels for genes encoding ACIII, cytochrome bd oxidase, and ATP synthase were also highest at night (Supplementary Figure S3B and Supplementary Table S2). All of these observations support the hypothesis that "Ca. T. thiotrophicus" performs anaerobic respiration at night. Interestingly, the relative transcript levels for genes encoding the aa 3 -type cytochrome c oxidase belong to expression pattern A, which are expressed at highest levels during light-dark transition periods in the evening and the morning (Supplementary Figure S4 and Supplementary Table S2). These transition periods correspond to times when the cyanobacteria are reducing dinitrogen to ammonia, but also represent periods when the oxygen levels in the mat change from oxic to anoxic (evening) and anoxic to oxic (morning).
Transcript levels for genes associated with carbohydrate metabolism, including those for glycolysis/gluconeogenesis and the TCA cycle, as well as those involved in acetaterelated metabolism (i.e., acetate kinase and acetyl-CoA acetyltransferase), also belonged to categories D and E and had highest transcript levels at night. However, two genes encoding glycogen synthase belonged to category B and were mostly unchanged throughout the day (Supplementary Figure S5 and Supplementary Tables S1, S2).
Fewer genes (n = 222) belonged to the smallest transcript category, C: genes that have highest transcript abundances during the day. Genes involved in oxygen protection, such as superoxide dismutase (MS gene ID 000129_2), peroxiredoxins (000011_73, 000020_4, 000236_9, 000026_30, 000091_19), glutathione peroxidase (000015_37), a nitroreductase (000041_8), and others (Supplementary Table S1), are included in this category (Supplementary Figure S6 and Supplementary  Table S2). Other day-genes in category C include the suf genes (000129_15-20), which encode proteins involved in iron-sulfur cluster biosynthesis under oxic conditions in many organisms (Supplementary Figure S6 and Supplementary Table S2) (Ayala-Castro et al., 2008) and, various genes coding for molecular chaperones (0054414-6) and several RpoE-like sigma factors (sigma factor 24) (MS gene ID 000006_60, 000002_32, 000100_29, 000020_22, 000036_46) possibly associated with stress conditions (Raina et al., 1995;Supplementary Table S1). The iron chelator hydroxamate (gene ID 004701_1) and a chelator transporter that might presumably be involved in harvesting iron for some of the enzymes mentioned above also showed higher relative transcription levels during the day (Supplementary Table S1, see the section "Discussion").
Genes encoding flagellar components were largely constitutively expressed, and most belong to diel expression category B with more or less unchanging relative transcription levels throughout the diel cycle. Thus, the transcription patterns did not support the hypothesis that flagella have a specialized function and are only used under specific growth conditions in the mat (Supplementary Table S1). Relative transcript levels for chemotaxis-related genes, however, were higher during the evening and were lower when conditions in the mat become oxic (category D, Supplementary Table S1).
In order to verify that the gene expression patterns observed were accurate representations of the "Ca. T. thiotrophicus" ecophysiology and not isolated events, we sampled and analyzed metatranscriptomes at selected times in the 24-h period immediately following the diel cycle examined in this study. These additional time points showed similar transcript levels to those measured in the more complete diel cycle analysis (second day data shown within all metatranscriptome figures; please see time points at 13:00, 15:00, 17:00 for both days data, and time point 14:00 for second day data only. Please note that for some processes, e.g., sulfate reduction, the data on both days are identical and time points overlap). Similarly, the metatranscriptomic analysis from Nakabusa included independent samples at time points on two different days as a control for day to day reproducibility. In the samples from Nakabusa, similar transcript levels were also found for samples taken at the same time points on consecutive days (data not shown).

DISCUSSION
"Ca. Thermonerobacter thiotrophicus": A Member of the Chlorobi-Lineage 5 (OPB56, "Ca. Kapabacteria") With Dissimilatory Sulfur Metabolism In this study we describe the initial genomic and transcriptomic characterization of "Ca. T. thiotrophicus, " a member of the Bacteroidetes/Chlorobi group from MS microbial mats with dissimilatory sulfur metabolism, and predicted sulfate-reducing directionality. "Ca. T. thiotrophicus" is the first described putatively SRM from the Bacteroidetes/Chlorobi group, among which to this point only (phototrophic) members with sulfideoxidizing, dissimilatory sulfur metabolism have been observed. SRM were first shown to occur in the microbial mats associated with MS on the basis of dsrAB cloning experiments, microelectrode measurements, and enzymatic assays (Dillon et al., 2007). Sequences from this earlier study, labeled as "unidentified clade-2, " clustered phylogenetically with the "environmental supercluster" of reductive DsrAB sequences. These sequences were also detected in this metagenomic study, and a MAG for the uncultured SRM was obtained. Additional MAGs containing highly similar dsrAB gene sequences to those in MS were identified in metagenomes from Nakabusa hot springs in Japan and the iron-rich, Chocolate Pots hot spring in YNP.
No cultured representative of the "environmental supercluster" of dsrAB sequences has yet been described, and no taxonomic classification of the dsrAB-containing organisms in MS had been possible until this study. Recently, (meta)genomics have allowed insights into the functional potentials of many uncultured organisms, and partial genomes of uncultured members of different phyla, including Acidobacteria, Nitrospira, Chloroflexi, "Ignavibacteria, " and a number of "Candidatus" phyla have been affiliated with this largely uncharacterized "environmental supercluster" Hausmann et al., 2018;Zecchin et al., 2018). Using ribosomal genes present in the Nakabusa MAG as well as other phylogenetic marker genes from the three MAGs studied, we were able to identify the uncultured, putative SRM from these hot springs as members of the "Chlorobi lineage 5, " also known as "clade OPB56, " or "Ca. Kapabacteria" (Iino et al., 2010;Kantor et al., 2015). Previous partial genomes of representatives of this clade suggested that these organisms are chemoheterotrophs that metabolize small organic molecules and are broadly distributed in thermal environments (Hiras et al., 2016). However, DSR was not indicated by those studies. We propose the name of "Candidatus Thermonerobacter" for this hot spring-associated SRM from the OPB56 clade, with the member from MS as "Candidatus Thermonerobacter thiotrophicus." Moreover, we propose that the dsrAB-containing "Ignavibacteria" genome of MAG GWB2_35_12 (GenBank WGS entry MGZV) identified by Anantharaman et al. (2018) should be re-classified as an OPB56 member, a sister clade to Ignavibacteraceae. Based on phylogenetic analyses of the 16S rRNA (Figure 4) and phylogenetic marker genes (Figure 3), this organism is very similar and shares most of the genomic properties with the organisms described in this study, although it comes from a very different environment, the groundwater aquifer system at Rifle, CO, United States. Some minor differences in the organization of the genes for sulfate reduction are observed (Figure 5), but overall this organism appears to be quite similar to those described here from hot springs.

Metabolic Prediction and Versatility/Activity in situ
Despite the low sulfate concentrations in the mat (<200 µM), substantial rates of sulfate reduction were measured in a previous study of MS (Dillon et al., 2007). The overall activity was attributed to Thermodesulfovibrio sp. as well as to organisms associated with three unidentified dsrAB clone sequences retrieved from the mats. In addition, although there is almost no measurable sulfide in the spring water (0.003 mg/l = 0.08 µM, USGS report -2002, the MS mats contain significant numbers of sulfideoxidizing chlorophototrophs (e.g., Chloroflexus sp., Roseiflexus sp., "Ca. Chloranaerofilum corporosum, " and possibly even Synechococcus spp.) (Thiel et al., 2016Tank et al., 2017), as well as sulfur-oxidizing members, such as Thermocrinis sp. , which would provide sulfate for DSR in the mats. Active sulfate reduction during evening and night hours is further supported by the observation that sulfide accumulates in the mat during the night but is undetectable during the day (E. Trampe, M. Lichtenberg, V. Thiel, M. Kühl, unpublished data;Dillon et al., 2007).
Given the predicted significance of a tightly coupled sulfur cycle, identifying organisms responsible for sulfate reduction had been an important goal of ongoing studies on the MS mats and those associated with similar hot springs. Our data indicate that "Ca. T. thiotrophicus" is involved in such sulfur cycling. Relative transcript levels for the dsr genes in this organism varied from night to day to a far greater extent than any other transcripts in these organisms, including those associated with other anaerobic processes such as nitrite and thiosulfate respirations. The patterns of night/day expression variation were consistent in the upper layers and undermats of both MS and Nakabusa hot spring. Interestingly, in the undermat of MS, at depths where oxygen concentrations never reach the hyperoxic levels of the uppermost part of the mat, transcripts levels for the genes associated with sulfur metabolism did not decrease as greatly during the day as they did in the upper green mat layer (Figure 7). This possibly indicates that these genes are regulated by environmental oxygen conditions and that active expression of the genes for sulfate reduction occurs at times when the mats are anoxic. In fact, most genes involved in energy conservation, including those involved in carbohydrate oxidation, electron transport and ATP synthase exhibited higher relative transcript levels at night than during the day.
In contrast to the large number of genes and processes with a nighttime transcription pattern, we found few genes that could be involved in energy conservation with higher relative transcript levels during the day. These include a putative 2-oxovalerate:ferredoxin oxidoreductase (VOR), an aldo/keto oxidoreductase with homology to oxidoreductases associated with the oxidation of amino acids, two genes coding for ferredoxin, and two genes with homology to cytochrome oxidase biosynthesis components (Supplementary Figure S7 and Supplementary Table S2). Given the lack of additional metabolic transcripts during the day, it is unclear whether the products of these genes participate in energy conservation, and if so, whether they contribute to a common metabolic process or act independently of one another. In contrast to the relatively small number of metabolic genes with highest transcript levels during the day, genes with products involved in stress responses, such as those associated with oxygen protection, and suf genes used to produce FeS clusters under oxygen-stress conditions (Jang and Imlay, 2010;Zang et al., 2017;Supplementary Figure S6 and Supplementary Table S2) are highly represented in the day-time expression cluster. These patterns reflect the significant oxygen levels this organism sees during the day. Lastly, the partial genomes encode two oxygen-dependent terminal oxidases (caa 3 and bd-type), of which cytochrome oxidase caa 3 shows a dual transcription peak pattern during the day, peaking in the morning and a second time in the evening (Supplementary Figure S5 and Supplementary Table S2). This might point toward active oxygen respiration in this organism, possibly during times of low oxygen concentrations. Although no previously isolated natural SRM has been demonstrated to grow by aerobic respiration, a new study has shown that this ability can evolve under oxygen-driven conditions (Schoeffler et al., 2018). This finding supports the idea that novel and uncultured hot spring organisms may also have this capability. As Desulfovibrio species contain the electron-transport complexes required for both sulfate and oxygen reduction (Figure 7), only a few mutations, presumably involving either inactivation or overexpression of the gene encoding heterodisulfide reductase, were required to permit the emergence of oxygen-respiration enabled growth (Schoeffler et al., 2018). Similarly, genes for both systems, sulfate and oxygen respiration, were detected in the "Ca. T. thiotrophicus" MAG under strongly aerobic day-time conditions, which would favor the evolution of a flexible, alternating facultative anaerobic metabolism. However, the observation that most genes, including those involved in electron transfer, are transcribed during the night might argue against energy production by aerobic respiration, at least under highly oxic conditions. Alternatively, this cytochrome oxidase might function in oxygen protection, as has been suggested to occur in some strictly anaerobic GSB (Li et al., 2009).
Sequences affiliated with "Ca. Thermonerobacter thiotrophicus, " the putative, sulfate-reducing member of OPB56 (Bacteroidetes/Chlorobi), have been detected in both the upper layer and undermat of MS. Based on amplicon data for MS-B_OTU-29 (Figure 4; Thiel et al., 2016) relative abundance of this member does not differ between the two layers. The upper layer, in particular, experiences a strong diel cycling of oxygen concentrations that can reach up to 800% air saturation. The presence and transcription of oxygen protection genes in the Bacteroidetes/Chlorobi member during the day support this hypothesis; most genes with high relative transcript levels during the day belong to the categories of oxygen and stress protection. Based on these data we predict that pathways of conserving energy in "Ca. T. thiotrophicus" are active during the anoxic periods of the diel cycle whereas, during the oxic periods, cells respond to oxygen stress by activating mechanisms that slow down or pause its metabolism. As expected from these predictions, shifts in relative transcript levels are not as pronounced in the undermat, where the oxygen levels do not reach the supersaturated levels in the upper green layer. However, an alternating aerobic/anaerobic metabolism involving sulfate as well as oxygen respiration cannot be excluded but can only be experimentally verified in future studies when pure cultures of these novel organisms are available.

SRM or SOM?
In natural systems, the oxidation state of sulfur can range from −2 to +6 and can support microbial energy metabolism as a reductant or oxidant. DSR is performed by microorganisms using sulfate as a terminal electron acceptor during anaerobic respiration ("sulfate respiration"), whereas in dissimilatory sulfur oxidation, sulfur compounds are used as electron donors in respiratory or photoautotrophic metabolism. Despite the different directionality of the reactions, a similar set of DSR genes is found in both SRM and SOM, which indicates a common origin from an early ancestor (Wagner et al., 1998;Grein et al., 2013;Müller et al., 2015). Although DSR genes have been found in both SRM and SOM, specific patterns and differences occur in the two groups.
The dsr genes encoding the Hrd-like menaquinol-oxidizing enzyme (dsrTMKJOP) as well as the genes sat-aprBA-qmoABC are not unique to SRM, but rather are also present in many SOM. Thus, the presence of these genes is not diagnostic for either SRM or SOM. The dsrAB genes, on the other hand, have long been used as phylogenetic markers to distinguish between these two metabolic modes, and to a more limited extent to identify known phylogenetic groups. Based on DsrAB phylogeny, the uncultured organisms described here were affiliated with "unclassified environmental sequences" within the "environmental supercluster 1" of unidentified and uncultured SRM (Müller et al., 2015). This affiliation is strongly suggestive that the hot-spring organisms perform DSR. Further support for this contention is found by the absence of the dsrEFH genes as well as aprM, which are usually present in SOM but not SRM. However, the recent findings of Thorup et al. (2017), who describe a novel deltaproteobacterium, D. alkaliphilus, which contains all genes phylogenetically characteristic of a SRM, but performs sulfide oxidation as its sole mode of dissimilatory sulfur metabolism, challenges the above-mentioned arguments. Nevertheless, the presence of all SRM "marker genes" as well as the absence of all SOM "marker genes, " in combination with the very high relative transcript levels for all genes for sulfate reduction during the anoxic night period, lead us to predict that these uncultured organisms perform sulfate reduction during periods when their mat environments are anoxic. Similarly, the relatively low transcript levels during oxic day times, as well as non-detectable sulfide concentrations in the microbial mats during the day (van der Meer et al., 2005; E. Trampe, M. Lichtenberg, V. Thiel, and M. Kuehl, unpublished data) do not support a metabolism based on sulfide oxidation coupled to aerobic respiration, as has been discussed for the dsrAB-carrying Acidobacteria, which carry a similar genome organization of the DSR pathway as observed here . Finally, the third alternative, using the DSR-related genes for anaerobic sulfide-oxidation coupled with dissimilatory nitrate reduction to ammonium (DNRA) as reported for D. alkaliphilus (Thorup et al., 2017), is considered to be unlikely because these organsims lack a nitrate reductase as well as any complete pathway for carbon fixation.

Implications for Our Understanding of the Evolution of Dissimilatory Sulfur Metabolism
For about a century, sulfur-oxidizing GSB were the only known members of the phylum Chlorobi, and the two terms were often used synonymously. The detection of DSR genes in putative sulfate-reducing chemoheterotrophic organisms within the same (super)phylum presents new and unexpected evidence in the evolutionary history of dissimilatory sulfur metabolism. Culture-independent studies, principally 16S rRNA gene analyses, have indicated that a far greater diversity exists within this phylum and have disclosed a substantial variety of uncultivated members affiliated within this group. However, based on the limited metabolic diversity of GSB, initially these distant relatives were incorrectly assumed to be sulfide-oxidizing, anaerobic chlorophototrophs (Iino et al., 2010;Liu et al., 2012a). The isolation and characterization of two chemoorganoheterotrophic strains, I. album (Iino et al., 2010;Liu et al., 2012a) and M. roseus Podosokorskaya et al., 2013), clearly demonstrated that the Chlorobi/Ignavibacteria-like (super)phylum encompassed much greater metabolic diversity than previously imagined. Presently, this group comprises five lineages in addition to Chlorobia. The two non-phototrophic isolates mentioned above are affiliated with "Chlorobi lineage 1, " which was proposed as a new class, Ignavibacteria (Iino et al., 2010), but it has also been suggested to represent a novel phylum, Ignavibacteriae . No isolated representatives have yet been described for any of the other four lineages, but partial genomes available for members of "Chlorobi lineage 5" (also known as OPB56 group; recently suggested to represent a novel phylum, "Ca. Kapabacteria;" Kantor et al., 2015) indicate that this lineage includes organisms predicted to be facultative anaerobic chemoheterotrophs, possibly with a putatively predatory lifestyle (Kantor et al., 2015;Hiras et al., 2016). However, no sulfatereducing or chemotrophic sulfur-oxidizing member of the Chlorobi/Ignavibacteria-like superphylum has yet been isolated and characterized.
As discussed above, GSB (and other SOM) contain a set of DSR genes that they utilize for sulfide/sulfur oxidation, but that are very similar to those of sulfate-reducing organisms. Interestingly, not all chlorophototrophic members of the Chlorobi contain DSR genes. Chloroherpeton thalassium, an early diverging chlorophototrophic member of the Chlorobiales, lacks the rDSR system (Frigaard and Bryant, 2008) and only possesses sulfide:quinone reductase, like sulfide-oxidizing members of the Chloroflexi. Furthermore, the iron-oxidizing anaerobes, Chlorobium ferrooxidans and Chlorobium phaeoferrooxidans and the photoheterotrophic aerobe, "Ca. T. aerophilum, " lack all of the genes for sulfide oxidation to sulfite. This and additional phylogenetic analyses suggest that the DSR pathway may have been acquired by the Chlorobi lineage after the divergence of Chloroherpeton-like organisms (Frigaard and Bryant, 2008;Liu et al., 2012b). On the other hand, DSR genes in sulfur-oxidizing, chlorophototrophic GSB, seem to represent a hybrid stage or an evolutionary chimera, in which some characteristics clearly are associated with SOM (such as their dsrAB phylogenetic affiliation), while other characteristics and genes are shared with SRM (such as dsrTMKJOP). This chimeric status of the DSR system in the sulfur-oxidizing, phototrophic Chlorobi suggests that this group did not inherit or acquire its genes from a typical SRM. In terms of the evolution of the DSR in the Chlorobi, three hypothetical scenarios can be considered: (1) the DSR genes were first acquired in one of the groups and were subsequently transferred horizontally to the other Chlorobi, thereby changing the utilization of the genes from oxidation to reduction or vice versa; (2) the DSR genes were acquired at an evolutionary earlier stage by a putative SRM, then loss of the genes occurred in some lineages and the genes were converted from reductive to oxidative dissimilatory metabolism in another lineage; or (3) a horizontal acquisition of rDSR genes in GSB after the loss of the sulfate-reducing set of DSR genes after the divergence from the chemotrophic Chlorobi lineage(s). It has been shown that GSB are rather susceptible to HGT, which could further support the hypothesis of a possible acquisition of DSR-related genes from two different sources (Llorens-Marès et al., 2017). It should further be noted that chlorophototrophy in the Chlorobi was probably obtained by HGT, and if this was the case in an organism that contained genes for DSR, there would have been strong selection pressure to alter the directionality of the reactions of sulfur metabolism. A detailed phylogenetic analysis of DSRrelated genes will be performed to address these questions in the future.

Concluding Comments
In this study partial genomes assembled from metagenomic data are presented, which describe putative, novel SRM, tentatively named "Ca. Thermonerobacter (thiotrophicus)" that occur in hot spring mat communities. The three MAGs described here contain complete (or nearly complete) sets of genes encoding all known enzymes required for DSR. Active sulfate reduction has previously been reported to occur in one of the three hotspring microbial mats. Metatranscriptomic studies performed on samples collected over a complete diel cycle for two of the hot springs showed that relative transcript levels for all 20 genes related to DSR were highest during the night. We hypothesize that the novel but presently uncultured organisms perform DSR at night. The presence of oxidative terminal cytochrome c oxidases in the genomes, and their transcription during the day (Supplementary Figures S5, S7C, Supplementary Table S2) might further indicate the facultative ability to also perform aerobic respiration. Future studies in situ and/or with enrichment or axenic cultures of these novel DSR-containing organisms will test this hypothesis and should conclusively demonstrate whether these "Chlorobi lineage 5" (OPB56, "Ca. Kapabacteria") members, "Ca. Thermonerobacter (thiotrophicus), " indeed are SRM or SOM.

AUTHOR CONTRIBUTIONS
VT conducted sampling, DNA and RNA extraction, metagenome binning for MS samples (YNP, WY, United States), as well as data analysis of MAGs. AGC conducted RNAseq data analysis for MS (YNP, WY, United States) and NK (Japan) samples and MAGs. NF conducted sampling, DNA extraction, and metagenome binning for CP samples (YNP, WY, United States). JM conducted sampling, DNA and RNA extraction, and binning for NK samples (Japan). MT conducted sampling for MS (YNP, WY, United States) and NK (Japan) samples, contributed to DNA and RNA extraction, data analysis and interpretation of MAGs. ER and EB planned and funded the study at CP (YNP, WY, United States). DW and DB planned and funded the study at MS (YNP, WY, United States). VT and SH planned and funded the study at NK (Japan). VT, AGC, and DB wrote the manuscript. All co-authors contributed to the manuscript with sections of their expertise, feedback, and fruitful discussions.

FUNDING
This study was partly funded by the Division of Chemical Sciences, Geosciences, and Biosciences, Office of Basic Energy Sciences of the Department of Energy through Grant DE-FG02-94ER20137. DB and DW additionally acknowledge support from the NASA Exobiology program (NNX09AM87G and NNX16SJ62G). This work was also partly supported by the US Department of Energy (DOE), Office of Biological and Environmental Research (BER), as part of BER's Genomic Science Program 395 (GSP). This contribution originates from the GSP Foundational Scientific Focus Area (FSFA) at the Pacific Northwest National Laboratory (PNNL) under a subcontract to DB. The nucleotide sequencing was performed as part of a Community Sequencing Program (Project CSP-411) and was performed by the US Department of Energy JGI, which is supported by the Office of Science of the US Department of Energy under Contract No. DE-AC02-05CH11231. Work by ER, EB and NF was supported by the NASA Astrobiology Institute, grant NNA13AA94A. EB was supported by the National Science Foundation, grant #1820658. Research by VT, MT, and SH was supported by the Institute for Fermentation, Osaka, Japan. JM was funded by the Tokyo Human Resources Fund for City Diplomacy.

ACKNOWLEDGMENTS
The authors would like to thank all of the JGI staff members who contributed to obtaining the sequence data from MS, members of the Wisconsin Biotech Center for sequencing the Chocolate Pots enrichment metagenomes, as well as the staff at DNALink for their service and generating the metagenome and metatranscriptome data for the Nakabusa hot spring mat samples. The materials from YNP used in this study were collected under permit no. YELL-SCI-0129 held by DW and YELL-2012-SCI-5972 held by ER administered under the authority of YNP. The authors especially thank Christie Hendrix and Stacey Gunther for their advice and assistance. They are grateful to Mr. Takahito Momose, the president of Nakabusa Onsen Inn Inc. (Japan), for letting them use their hot springs.

SUPPLEMENTARY MATERIAL
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb. 2018.03159/full#supplementary-material FIGURE S1 | Phylogenetic tree based on DsrAB amino acid sequences displaying the phylogenetic affiliation of the putative thermophilic SRM Bacteroidetes-Chlorobi group member (red box within the Environmental Supercluster in the reductive-type DsrAB clade). Sequences of GSB, SOM members of the Chlorobi cluster with the rDSR (sulfur oxidizing) DsrAB clade. The tree was calculated from with the publicly available DsrAB ARB database (Müller et al., 2015; http://www.microbial-ecology.net/download). New sequences were added using the parsimony method without changing the tree topology.
FIGURE S2 | Heatmaps showing expression patterns of "Ca. T. thiotrophicus" genes throughout a diel cycle in the upper green layer of MS (sampling times shown at the bottom). Red = high relative transcript abundance; green = low relative transcript abundance; black = no transcripts detected. Higher color intensity reflects relatively higher transcript levels. Genes are clustered according to their major patterns of expression as described in the text. Red FIGURE S3 | Relative transcript levels (y-axis) across a diel cycle (x-axis) for in the upper green layer of MS of genes associated with anaerobic respiration (A) and other electron transfer processes (B). White areas indicate light periods during the day and irradiance intensities in µmol photons m −2 s −1 (400-700 nm, secondary y-axis), while gray areas indicate dark periods during the night.
FIGURE S4 | Relative transcript levels (x-axis) of aa 3 -type cytochrome c oxidase across a diel cycle (x-axis) in the upper green mat layer of MS. White areas indicate light periods during the day and irradiance intensities in µmol photons m −2 s −1 (400-700 nm, secondary y-axis), while gray areas indicate dark periods during the night.
FIGURE S5 | Relative transcription levels (y-axis) of genes involved in glycolysis, the TCA cycle, and other central carbon metabolism pathways across a diel cycle (x-axis) in MS MAG (MS-B_bin-24), in the upper green mat layer of MS. White areas indicate light periods during the day and irradiance intensities in µmol photons m −2 s −1 (400-700 nm, secondary y-axis), while gray areas indicate dark periods during the night.
FIGURE S6 | Relative transcript levels (y-axis) of genes associated with oxygen stress/oxic conditions during a diel cycle (y-axis) in the upper green mat layer of MS. White areas indicate light periods during the day and irradiance intensities in µmol photons m −2 s −1 (400-700 nm, secondary y-axis), while gray areas indicate dark periods during the night.