Diel Patterns in Marine Microbial Metatranscriptomes Reflect Differences in Community Metabolic Activity Over Depth on the Continental Shelf of the North Atlantic

Microbial community dynamics are influenced not only by biological but also physical and chemical phenomena (e.g., temperature, sunlight, pH, wave energy) that vary on both short and long-time scales. In this study, samples of continental shelf waters of the northwest Atlantic Ocean were periodically collected from pre-sunrise to post-sunset and at multiple depths over summers of 2016 and 2017. Metatranscriptomic analyses revealed expression of photosynthetic genes in surface water samples corresponding to a diel relationship with sunlight. Photosynthetic genes originated from known phototrophs including Aureococcus, Ostreococcous, Synechocococus, and Prochlorococcus. Photosynthetic gene expression occurred pre-sunrise, suggesting the community initiates transcription before sunlight exposure, ostensibly to harvest energy more efficiently when the anticipated increase in light occurs. Transcripts from photoheterotrophic members of the SAR11 clade were also documented in surface samples, with rhodopsin expression being more abundant pre-sunrise and post-sunrise. Conversely, samples taken from the aphotic layer exhibited expression of transcripts related to nitrification that did not vary over the diel cycle. Nitrification gene transcripts, specifically amoA, nirK, hao, and norAB, were taxonomically related to well-known genera of ammonia oxidizers, such as Nitrospira, Candidatus Nitrosomarinus, Nitrosospira, and Nitrosopumilus. Overall, this study documents the role of light (varying with time and depth) in shaping the photosynthetic microbial community activity in the surface ocean, and further demonstrates that this diel regulation of photosynthesis is decoupled from the activity of the nitrifying microbial community in deeper and darker waters.


INTRODUCTION
Microbial communities change over short-term time scales to occupy ecological niches within marine ecosystems. Examining shifts in community gene expression over short-time periods can allow for insight into tightly coordinated ecosystem-level dynamics that orchestrate daily biogeochemical fluxes and into how they may respond in the face of rapid natural and anthropogenic environmental change. Changes in marine plankton community function are largely driven by variations in environmental conditions (e.g., physical, chemical) with implications for food web dynamics and global biogeochemical cycles (Barton et al., 2013). In the surface ocean, the relatively short time scales of diel changes in light intensity affect a multitude of different physical variables in the environment.
The function and activity of the photosynthetic community in the ocean surface has been shown to respond temporally with light exposure. As light intensity varies over the course of a day, the dynamics of photosynthesis, DOC production and utilization, nutrient uptake and regeneration and community metabolism should all be reflected in the corresponding dynamics of the metatranscriptome. For example, Aylward et al. (2015) compared daily metatranscriptome profiles of related microbial populations inhabiting the surface waters of a costal California upwelling region and oligotrophic North Pacific Subtropical Gyre (NPSG) and found that the dominant photoautotrophic microbes in each environment were central determinants of overall community transcriptome dynamics. The study also found that heterotrophic bacterial clades common to both ecosystems displayed conserved, genome-wide inter-and intrataxon transcriptional patterns and diel cycles (Aylward et al., 2015). Changes in light also have been shown to impact amino acid metabolism, DNA processing, translation, transporters, and vitamin metabolism (Ottesen et al., 2013)-tightly regulated metabolic processes responding to intense resource competition. Viruses have also been shown to have diel oscillations in surface ocean microbial communities, resulting additional recourse competition within this ecosystem . Overall, nutrient availability and the timing of growth and DNA replication are key factors playing into regulation of temporal dynamics in surface planktonic assemblages (Aylward et al., 2015). Similarly, Vislova et al. (2019) investigated depthdependent changes in diel patterns for primary and secondary production, metabolic activity, and gene expression in marine microbial communities in the NPSG and found diel transcript oscillations in microbial communities decreased in abundance and regularity with depth and light levels.
Changes in light availability can also impact nitrifying microbial communities throughout the ocean water column. The first step of nitrification is the oxidation of ammonia to hydroxylamine catalyzed by amoA (Ahlgren et al., 2019) followed by the oxidation of hydroxylamine by hydroxylamine oxidoreductase (hao) to produce nitrite in ammonia oxidizing bacteria (Poret-Peterson et al., 2008). While bacterial nirK and norAB are involved in the conventional denitrification pathway, a number of nitrifying ammonia-oxidizing archaea (AOA) and ammonia-oxidizing bacteria (AOB) also possess these genes, which are believed to be involved in so-called nitrifierdenitrification (reduction of NO 2 − to NO and N 2 O) (Wijma et al., 2004;Kozlowski et al., 2014;Lehtovirta-Morley, 2018). As nitrogen availability often limits primary productivity in both marine and terrestrial ecosystems (Canfield et al., 2010), understanding of the nitrogen cycle is critical to understanding its role in regulating primary production in many regions. Diel periodicity of nitrogen cycling in the upper ocean has been shown in many previous studies (Wheeler and Kokkinakis, 1990;Cochlan et al., 1991;Bebout et al., 1994;Brandes et al., 2007;Ward, 2008;Devol, 2015;Wilson et al., 2017), however, much less is known about possible diel periodicity of nitrogen cycling in dark layers just beneath the euphotic ocean.
The coordinated activity of microbial communities operating over daily timescales results in temporal dynamics in the activity of specific microbial taxa and metabolic pathways on diel cycles. In terms of geographic range, metatranscriptomic studies have been conducted in different locations in both the Atlantic Ocean (Gilbert et al., 2010;Hewson et al., 2010), and Pacific Ocean (e.g., Poretsky et al., 2009;Shi et al., 2011;Ottesen et al., 2013Ottesen et al., , 2014. Further, many metatranscriptomic studies have focused on either one sample taken at one timepoint and depth (Frias-Lopez et al., 2008;Wu et al., 2013;Ji et al., 2018). To date, only one study has focused on spatial and temporal changes in microbial community gene expression (Vislova et al., 2019). While these studies used different sampling designs, many focused on the functionality of the microbial community and have examined essential metabolic functions that impact biogeochemical cycles, including carbon and nitrogen cycling (Frias-Lopez et al., 2008;Poretsky et al., 2009;Hewson et al., 2010;Shi et al., 2011;Ottesen et al., 2013Ottesen et al., , 2014Wu et al., 2013;Aylward et al., 2015;Wilson et al., 2017;Ji et al., 2018;Vislova et al., 2019). Here, we expand upon these studies by examining how microbial community gene expression relating to photosynthesis and nitrification changes as a function of light as regulated over both space (i.e., depth in the water column) and time (i.e., diel cycle) at a single coastal site in the Continental shelf waters off New England (United States).

Sample Site and Collections
Water samples were collected from the North Atlantic Ocean, near the coast of New Jersey (United States) from the R/V Endeavor in August 2016 (cruise EN588, 39.417, -73.258) and August 2017 (cruise EN604, 39.188, -73.561). Samples were collected using a paint-sealed Sea-Bird rosette, outfitted with a Sea-Bird SBE 19 Plus V2 recording CTD (conductivity, temperature, depth) and 8-L X-Niskins (General Oceanics). The Niskins were acid-cleaned in the laboratory prior to use and "pre-soaked" for 12 h with local seawater collected at a shakedown station prior to deployment. This "clean" rosette was lowered on a non-metallic line (1/4 O.D. vinyl-jacketed Vectran) and bottle closure at programmed depths controlled by a Sea-Bird Auto-firing Module. Sample depths were chosen based on the real-time water column profile data obtained from the ship's CTD rosette system (Sea-Bird) that was deployed before the "clean rosette" and which also provided dissolved oxygen, salinity, photosynthetically active radiation (PAR), beam transmission, and fluorescence profiles. Samples were collected at three different depths for the 2016 cruise (5, 16, and 30 m) and 2017 cruise (3, 30, and 40 m) based on the PAR and fluorescence profiles. The 5 m samples of the 2016 cruise and the 3 m samples of the 2017 cruise were within surface sunlit layers ("photic"), the 16 m sample of the 2016 cruise was at the peak of the deep chlorophyll max ("DCM") layer (as defined by the WETlabs ECO Chlorophyll Fluorometer), and the 30 m samples from both the 2016 and 2017 cruises and the 40 m samples from the 2017 cruise were collected at depths where the PAR was 1% or lower (hereafter "aphotic"). To examine changes at these depths as a function of relative sunlight intensity (sunlight terms are only relative as PAR was not collected at each time point), samples from the 2016 cruise were collected at 04:00 (pre-sunrise), 06:30 (post-sunrise), 14:00 (peak sunlight), 19:15 (pre-sunset), and 21:30 (post-sunset). For the 2017 cruise, samples were taken at 05:00 (pre-sunrise), 09:30 (post-sunrise), 14:00 (peak sunlight), 18:45 (pre-sunset) and 22:00 (post-sunset). For both cruises, samples were transferred from Niskin bottles to a 20 L carboy and immediately filtered as described below.
Samples for dissolved organic carbon (DOC) were passed through combusted 0.2 µm GF/F filters into acidified combusted vials and stored at 4 • C until they were analyzed by a Shimadzu TOC-V analyzer (Kyoto, Japan). Samples (250 mL) for chlorophyll a quantification were collected on a 25 mm GF/F filters and stored at -20 • C. Chlorophyll a was extracted from the filter in 90% ethanol in the dark at 4 • C (9 h) measured with a 10AU fluorometer (Turner) with appropriate standards (Sigma) with the data being corrected for phaeopigments (Parsons and Strickland, 1972). Water samples for inorganic nitrogen speciation were filtered through a 0.2 µm filter, with aliquots stored frozen or analyzed immediately. Subsamples of filtered water were used for immediate shipboard measurement of ammonium (NH 4 + ) by fluorescence using the OPA method (Holmes et al., 1999) and nitrite (NO 2 − ) by spectrophotometry using the Griess reaction method (Pai et al., 2021), with detection limits of 0.02 and 0.1 µm, respectively. Frozen aliquots of nitrate (NO 3 ) were measured in the lab by reduction to NO and quantification by chemiluminescence (Braman and Hendrix, 1989), with a detection limit of ∼0.7 µm.
Samples for RNA analysis were first filtered through a 2.7 µm filter (Whatman) and then a 0.2 µm Sterivex TM filter (EMD Millipore). The Sterivex filter was immediately treated with RNA later (Qiagen) and the samples were flash frozen in liquid nitrogen. Samples were then shipped on dry ice to Central Michigan University and stored at -80 • C until processing.

RNA Extraction and Sequencing
RNA was extracted from each Sterivex filter using a MirVana miRNA isolation kits (Ambion) following the manufacturer's instructions. The extracted RNA was treated with DNase (Qiagen, Valencia, CA, United States) and then cleaned with RNeasy MinElute kit (Qiagen, Valencia, CA, United States). Total RNA was quantified using a Qubit 2.0 fluorometer. For sequencing, RNA samples were sent to the Michigan State University's Research Technology Support Facility Genomics Core (MSU RTSF). At the sequencing core, RNA was assessed with an Agilent 2100 Bioanalyzer RNA 6000 Pico Chip and RNA-Seq libraries were generated using the Ovation Complete Prokaryotic RNA-Seq System (NuGEN), which included a rRNA depletion step. Shotgun sequencing of the cDNA was performed on one lane of an Illuminia HiSeq 4000 with 2 × 150 bp chemistry. The resulting raw reads are deposited in the National Center for Biotechnology Information's (NCBI) Short Read Archives (SRA), BioProject ID: PRJNA573842 (accession numbers SAMN12827068-SAMN12827096).

Metatranscriptomic de novo Assembly and Annotations
The resulting sequencing reads were quality checked with FastQC (Holmes et al., 1999) and trimmed using Trimmomatic v0.33, with the default parameters (Bolger et al., 2014). Paired-end reads were interleaved into a single file using Khmer v2.1.1 (Crusoe et al., 2015), and rRNA sequences were removed with SortMeRNA v. 2.0 (Kopylova et al., 2012). The resulting FASTA files were combined by depth and assembled using rnaSPAdes v. 3.13.0 (Bushmanova et al., 2019). QUAST v. 5.0 (Mikheenko et al., 2016) was used to assess the quality of the assembled files (Supplementary Table 1). The depth-specific transcriptome assemblies generated from rnaSPAdes v. 3.13.0 (Bushmanova et al., 2019) were then submitted to the Department of Energy's Joint Genome Genomes Online Database (GOLD) (Mukherjee et al., 2019) for annotation (IMG Genome ID's: 3300033603, 3300033604, 3300033605, 3300034405, and 3300034345). The three annotated datasets were then compared to each other using the JGI Institute Integrated Microbial Genomes (IMG) web interface. Functional groups (e.g., photosynthesis, carbon metabolism) from each layer were identified by the KEGG orthology (KO) database or the COG database (translation, ribosome structure, and biogenesis), both within IMG. Salmon v. 0.11.3 (Patro et al., 2017) was used to quantify the number of raw reads that matched to the annotated contigs in each of the layers and normalized the output as transcripts per million (TPM). The merged counts output file from Salmon for each assembly was then merged by locus tag to the genes identified and extracted in IMG using the statistical program R v. 3.5.1 (R Core Team, 2020). Taxonomic identifications (to their best match within the database) were performed by submitting the gene FASTA amino acid (.faa) files for each annotated gene to the GhostKOALA server (Kanehisa et al., 2016). Expression levels for each taxonomic group were found by aggregating the GhostKOALA output with the Salmon expression TPM counts to show total gene expression levels within each listed taxonomic group.

Chemical and Physical Data
In both cruise years, surface water temperatures ranged between 27 and 25 • C and decreased to 8 • C at the deepest sampling  Table 2). In addition, both ammonium and nitrite were low in the surface water samples but increased with depth in both years (Supplementary  Table 3). Of the protein coding genes, 36.86% had a function prediction associated with the annotation and 62.20% had no function prediction (Supplementary Table 4).
The 2017 cruise transcripts contained a total of 4,101,045 annotated genes (Supplementary Table 3). Gene annotations were 1,733,133 (surface transcripts), 2,367,912 (aphotic transcripts), and 1,629,343 (deep aphotic transcripts) (Supplementary Table 3). Of the protein coding genes, 29.91% had a function prediction associated with the annotation and 69.44% had no function prediction associated with the annotations (Supplementary Table 4).

Photosynthesis Gene Expression
During the 2016 cruise, gene transcripts annotated as Photosystem I and II in the surface layer (5 m) had measurable expression levels at each time point, including pre-sunrise and post-sunset. Within the surface waters, an increase in expression occurred post-sunrise and peaked at the time of peak sunlight (Figure 2). Similarly, expression of ATPases and ribosomal proteins had peak expression post sunrise with relatively high expression at peak sunlight (Supplementary Figure 1). The DCM layer (16 m) also had peak expression of both Photosystem I and II occurring post-sunrise (Figure 2).
Transcripts of psbA, which encodes the D1 protein of photosystem II, in the surface water were allied with Stramenopile (Aureococcus), cyanobacteria (Synechocococus), and green algae (Ostreococcus). While expression of psbA associated with Synechococcus and Ostreococcus was highest during times of peak sunlight, Aureococcus psbA transcripts had highest levels before peak sunlight (Figure 3). Transcripts of psbA collected in the DCM were dominated by Aureococcus, with peak expression occurring post-sunrise (Figure 3). The taxa most highly expressing psaA, integral membrane protein of Photosystem I, in the surface layer was Aureococcus with the most expression occurring post-sunrise (Figure 3). Synechococccus had high psaA expression during times of sunlight, while Prochlorococcus had peak expression pre-sunset (Figure 3).
For the 2017 cruise, transcripts associated with Photosystem I and II surface layer were highest in the early morning (presunrise through to peak sunlight) and then relatively decreased to post sunset (Figure 2). While Photosystem II transcript levels within surface water in 2016 were highest during peak sunlight, the highest transcript levels were observed postsunrise in 2017 (Figure 2). Furthermore, in the 2017 surface samples the photosynthetic cytochrome B6/F complex showed peak expression pre-and post-sunrise with expression gradually declining throughout the day (Figure 2). ATPase and ribosomal protein expression also had relatively high expression pre-and post-sunrise (Supplementary Figure 1).
The surface water samples in the 2017 cruise showed the most dominant taxa expressing the psbA gene belonged to Aureococcus, Synechococcus, and Ostreococcus (Figure 3). Similarly, psaA surface transcripts were nearly consistent throughout the day for Aureococcus and Ostreococcus, but highest at peak sunlight for Synechococcus (Figure 3). Unlike the 2016 cruise photic samples, psbA and psaA were not expressed by Prochlorococcus in 2017 surface samples (Figure 3).

Photoheterotrophy Gene Expression
Photoheterotrophy was examined by assessing microbial rhodopsin expression. Specifically, rhodopsin was found in the photic layers in both years (Figure 2). Relative to the expression of genes annotated in PS I and II, rhodopsin expression in 2016 was lower, but peak expression was post sunrise in the surface layer and pre-sunrise in the DCM layer (Figure 2). The abundant rhodopsins were taxonomically identified as Alphaproteobacteria HIMB5, Alphaproteobacteria HIMB59, and Candidatus Peligibacter sp. (Figure 3). In the surface layer of the 2016 cruise, rhodopsin transcripts from Alphaproteobacteria HIMB5 showed expression throughout the day, with the highest peak expression occurring at post sunrise (Figure 3). There were no transcripts documented for other microbial rhodopsin (such as proteorhodopsin or halorhodopsin).
Rhodopsin expression patterns in 2017 were less obvious than 2016, and total rhodopsin transcripts within surface samples showed peak expression pre-sunrise (Figure 2). Similar to 2016, rhodopsins in the 2017 surface transcripts were taxonomically identified as Alphaproteobacteria HIMB5, Alphaproteobacteria HIMB59, and Candidatus Peligibacter (Figure 3) with peak expression either pre-or post-sunrise (Figure 3). Similar to 2016, there were no transcripts identified for other microbial rhodopsin (such as proteorhodopsin or halorhodopsin).

Ammonia Oxidation
Samples taken from depths deeper then 30 m from both years clearly documented the presence of genes and organisms associated with nitrification (Figure 4). In 2016, transcripts annotated as amoABC, norAB, and hao were highly expressed (Figure 4). Nitrosospira had the highest gene expression of amoA followed by Candidatus Nitrosomarinus (Figure 5). The most highly expressed transcripts of hao in the aphotic samples were taxonomically identified as Nitrosomonadaceae, Nitrosospira and Nitrosococcus (Figure 5). The major difference between samples collected in 2017 compared to 2016 was a lack of hao expression ( Figure 5). Nitrosopumilus only expressed amoA at the 30

DISCUSSION
It is well known that the availability of light regulates primary production by photosynthesis in the surface ocean (Poretsky et al., 2009;Hewson et al., 2010;Aylward et al., 2015). The energy and biomass generated via these processes varies temporally and spatially. As phototrophic organisms generate organic matter, exuded metabolites and detrital cellular material support a diverse heterotrophic community (Ottesen et al., 2014). Heterotrophic remineralization of organic matter produces reduced nitrogen compounds (especially NH 4 + , but also urea, cyanate) that stimulates and supports a chemolithoautotrophic community of nitrifying organisms, including AOA and bacteria (AOB), and nitrite-oxidizing bacteria (NOB) (Bebout et al., 1994;Zehr and Ward, 2002). The metabolic activity of each of these groups is reflected through expression of different metabolic pathways. Here, we present our examination of metatranscriptomic patterns gene expression of photosynthesis and nitrification over scales of space (depth) and time (over the course of a day and a yearly comparison).

Photosynthesis
In this coastal study, community transcripts related to photosynthesis (photosystem II and photosystem I) displayed peak gene expression occurring times of relatively higher sunlight (post-sunrise and peak-sunlight) and with lowest gene expression near sunset (Figure 2 and Supplementary Figure 1). Transcripts from photosynthetic organisms in the open ocean's surface have been shown to vary with light intensity along a diel cycle (Frias-Lopez et al., 2008;Ottesen et al., 2013Ottesen et al., , 2014Aylward et al., 2015), however, the current study did not always observe peak expression with peak sunlight. This difference could be due to variation in sample collection time, sample number, sequencing depth, or that the current study was conducted in a coastal marine ecosystem. While transcript peaks occurred near peak sunlight in this study, expression of photosynthetic genes (PS I and II) have also been documented to be expressed under conditions such as pre-sunrise (Hewson et al., 2010;Becker et al., 2019). Such dark expression likely reflects the physiological preparation of the community to efficiently harvest light as it becomes available after sunrise.
Expression of photosynthetic genes within surface waters was genetically similar to four known phototrophic microbes (Figure 3). The dominant (highest relative abundance) psbA (photosystem I) and psaA (photosystem II) transcripts in the surface and DCM layers belonged to the genera of Aureococcus, Ostreococcous, Synechocococus, and Prochlorococcus. The abundance of photosystem transcripts from certain taxa, such as Aureococcus in the DCM and Synechocococus in the surface, exhibited peak expression during peak sunlight. Similar to community transcript levels, Aureococcus and Ostreococcous within surface samples had relatively high expression of psbA and psaA pre-sunrise to peak sunlight. Both Ostreococcus (Ottesen et al., 2014;Aylward et al., 2015) and Prochlorococcus (Ottesen et al., 2013;Aylward et al., 2015;Vislova et al., 2019) have been shown to be dominant community members in surface waters of the Pacific Ocean with their expression following a diel pattern (synchronous with sunlight), but photosynthetic genes were still expressed pre-sunrise, as also observed here in the Atlantic Ocean.

Photoheterotrophy
Transcripts for production of rhodopsin, involved in harnessing light for energy (Spudich et al., 2000;Neutze et al., 2002;Pinhassi et al., 2016), maintained relatively stable expression throughout the day in photic zone samples of the 2016 cruise, with surface samples exhibiting highest levels of expression pre-sunset/post sunrise and the DCM samples showing high expression pre-sunrise and post sunset (Figure 2). Similar to photosynthic genes, the lack of peak expression during peak sunlight as seen in other studies could be due to differences in experimental design or ecosystem variability, however, two other studies have also documented peak transcript expression from phototrophs during darker periods (Ottesen et al., 2013;Becker et al., 2019). The abundant transcripts for the gene that produces rhodopsin observed in waters from the photic zone (surface and DCM) were taxonomically aligned with Alphaproteobacteria HIMB5, Alphaproteobacteria HIMB59, and Candidatus Peligibacter (Figure 3). All three of these taxa are members of the photoheterotrophic SAR11 clade (Grote et al., 2012;Zhao et al., 2017).

Ammonia Oxidation
While high levels of nitrification are generally not observed in the sunlit surface ocean, our data document transcripts suggesting an active nitrifying community at depths just below these sunlit layers (Figure 4). Taxonomic data from functional genes demonstrate that both AOA and AOB are active in the aphotic samples from the 2016 cruise (Figure 5). Both amoA and nirK allied to marine ammonia-oxidizing Thaumarchaeota Nitrosopumilus and Candidatus Nitrosomarinus (Herber et al., 2020). Transcripts of amoA and hao related to known AOB, Nitrosospira and Nitrosococcus (Hollibaugh et al., 2002;Klotz et al., 2006), and Nitrosomonadaceae (Prosser et al., 2014), were also identified ( Figure 5). Vislova et al. (2019) also noted that activity of nitrifying taxa did not vary with diel periodicity.

Comparison of Year-to-Year Variability
While back-to-back cruises in the Fall of 2016 and 2017 provide an opportunity to compare yearly variability in diel activity at this coastal site, differences in physical (mixing) and chemical (oxygen) conditions were also evident and relevant (see Figure 1 and Supplementary Table 2). In particular, the waters in 2016 were warmer, more productive, and therefore had shallower light penetration than in 2017. As such, the transition from photic to aphotic depths (as defined here) began at 12 m in 2016 but was at 28 m in 2017. Further, the surface waters in 2017 were less mixed as observed in the defined layers in the temperature and salinity plots (Figure 1).
Similar expression trends in photosynthetic transcripts documented in 2016 were also observed in 2017. The 2017 cruise collected water from the surface and two deeper aphotic layers, and not the DCM as in 2016. Thus, comparisons between the years are solely focused on the photic and deeper, aphotic samples. Similar to the 2016 cruise, the dominant psbA (photosystem I) and psaA (photosystem II) transcripts in the surface layer of the 2017 cruise belonged to the genera of Aureococcus, Ostreococcous, and Synechocococus with expression being still relatively high pre-sunrise. Further, surface transcripts collected in 2017 documented the presence of photoheterotrophs (Alphaproteobacteria HIMB5, Alphaproteobacteria HIMB59, and Candidatus Peligibacter), with higher expression still occurring during darker times of day. The presence of both phototrophs and photoheterotrophs suggests a coupling of metabolic processes, which enable the community to both fix and metabolize carbon.
Over a 2-year period both AOA, Candidatus Nitrosomarinus, and AOB, Nitrosospira, were found in these coastal euphotic waters. AOB are thought to thrive in more nutrient rich ecosystems and AOA could thrive in more oligotrophic ecosystems (Santoro et al., 2008;Martens-Habbena et al., 2009) as AOA amoA genes have higher affinities for ammonia (Martens-Habbena et al., 2009;Straka et al., 2019). In the current study, primary controlling factors (e.g., temperature, pH, substrate concentrations) did not vary year to year. Thus, it is possible that the shelf waters examined here simply represent a midpoint along the continuum between the nutrient replete coastal waters (where AOB frequently outnumber AOA) and nutrient depleted waters further offshore, where AOA are typically the dominant nitrifiers, with physical conditions not allowing either group to establish dominance.
Periodic and daily transcriptional trends in microbial metabolism occur in marine communities as a function of light attenuation. As a result, microbial communities in the ocean's surface follows a diel pattern of expression of gene transcripts related to metabolism. Interesting, photosynthetic genes described here also were expressed in pre-sunrise, ramping up for future times of more sunlight. However, below the photic zone, ocean microbial communities show more stable gene expression of nitrification. Given that both cruise (sampling separated by 1 year) transcripts showed similar expression trends, this suggests that in the photic zone microbial community function follows a consistent pattern through time. This work also focused on the coastal Atlantic Ocean, with other studies focusing on the Pacific Ocean (Ottesen et al., 2013(Ottesen et al., , 2014Aylward et al., 2015;Vislova et al., 2019), which suggests these expressional patterns occur globally.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih. gov/, PRJNA573842.

AUTHOR CONTRIBUTIONS
JZ conducted the bioinformatics, analyzed the data, and drafted the figures. DH performed the laboratory work to extraction RNA from the samples. CH, SW, CL, and DL collected the samples. CS provided support and insight for bioinformatic analysis. CH and DL designed and planned the study. All authors contributed to writing the manuscript.

FUNDING
This work was supported by the National Science Foundation grant # 1355720 (CH and CL) and from Central Michigan University Faculty Research and Creative Endeavors Committee (DL).