Oxygen microsensors were used to characterize the balance between respiration and photosynthesis in the green berry aggregates. Aggregates examined were relatively symmetric ellipsoids of similar size with an equivalent spherical diameter of 1.7 ± 0.1 mm . Photosynthesis produced supersaturated oxygen concentrations within the aggregates: 380 µM O2 with illumination at 170 µE m−2 s −1 and 520 µM O2 at 320 µE m−2 s −1 . During these experiments, bubbles were not observed on aggregate surfaces and the aggregates never floated. Oxygen production rates per aggregate were calculated from these profiles as 13 and 31 nmol O2 per hour at 170 µE m−2 s −1 and 320 µE m−2 s −1 , respectively . We did not rigorously determine the saturating light intensity; however, we observed that illumination with a third lamp failed to stimulate increased oxygen production beyond that with two lamps , suggesting that the saturating light intensity lies in the range of 170 – 320 µE m−2 s −1 . The green berries were anoxic in darkness with a dark respiration rate of 19 nmol per aggregate per hour . This rate, equivalent to a volumetric rate of 6.8 µmol cm−3 hr−1 , is well above the theoretical threshold for respiration rates capable of causing oxygen diffusional limitation from surrounding seawater . In light-dark shifts,flower display buckets we observed a rapid response where the aggregate core transitioned from steady state supersaturation to full anoxia in 5–8 min .
Volumetric gross photosynthetic rates were calculated by two methods: from depth microprofiles via the sum of net photosynthesis and dark respiration and via the light-dark shift technique performed at a single point in the aggregate center . At both light intensities examined, rates calculated via the light-dark shift method were found to be 4.5 µmol cm−3 hr−1 , lower than those from depth microprofiles. While this difference could arise from biological variability between aggregates, we suspect that the light-dark shift rates measured at the aggregate core were lower than those we might have measured closer to the aggregate surface. Future depth integrated studies of photosynthetic rates will helpto clarify this difference and allow better characterization of respiratory activity in the light. Comparing the green berries’ dark respiration and gross photosynthesis to other photosynthetic mats and aggregates, we find them similar to the high rates measured for large , filamentous aggregates of the heterocystous cyanobacterium, Nodularia spumigena from the Baltic Sea . Indeed, our estimates of carbon fixation are close to prediction of 349 ng C per aggregate per hour calculated using Ploug et al.’s regression of volume to gross photosynthesis from a 2009 Nodularia bloom. The green berry dark respiration rate was similar, though slightly higher than that observed for Nodularia aggregates of similar diameter; however, the ratio of dark respiration to gross photosynthesis of 0.38 was identical to that observed for 7 mm3 aggregates found late in the bloom . In large Nodularia aggregates, dark anoxia is associated with active dissimilatory nitrogen cycling, including denitrification, dissimilatory reduction of nitrate to ammonia , and significant rates of cryptic nitrification .
Our metagenomic data analysis predicts genes in GB-CYN1 involved in the Calvin cycle, TCA cycle, and photosystems I and II . Sequence identities for these genes in GB-CYN1 and published Chroococcales genomes were high, ranging from 75 to 99% . The observation of photosystem II genes in GB-CYN1 indicates the metabolic potential for oxygenic photosynthesis, unlike UCYN-A which lacks photosystem II, RuBisCo and the TCA cycle. These findings confirm that GB-CYN1 resembles free-living Chroococcales of the clades B and C rather than the metabolically streamlined, endosymbiotic UCYN-A clade . Related Chroococcales species are known to produce copious quantities of exopolysaccharides , and EPS produced by Crocosphaera watsonii has been linked to the formation of transparent exopolymeric particles in pelagic environments . A cassette of genes required for EPS production has been identified by comparative genomic analysis of related Crocosphaera watsonii strains ; however, we did not recover homologs of these genes in our metagenomic sequence . We conclude that the absence of these sequences is most likely due to our incomplete sampling of the GB-CYN1 genome, though it could also indicate either an unknown pathway for EPS synthesis in GB-CYN1, or an alternate source of the green berry exopolymeric matrix . Consistent with the dynamic oxygen conditions within the green berries, we found sequences in the metagenomic data that suggest the metabolic potential for aerobic and anaerobic respiration, and anaerobic fermentative pathways assigned to several different phyla of bacteria . Abundant fixed nitrogen from diazotrophy and transient anoxia presents an ecological opportunity for dissimilatory nitrogen metabolism, a process that could cause concomitant nitrogen fixation and loss over rapid spatiotemporal scales.
Indeed, previous studies of diazotrophic Nodularia spumigena aggregates measured significant rates of both nitrification and denitrification . We specifically investigated potential marker genes for dissimilatory nitrogen metabolism within the green berries. Metagenomic sequences homologous to the nitrite reductase gene nirK were found in six unassembled reads and on a corresponding 753 bp long contig . A single read homologous to nitrous oxide reductase nosZ was also found suggesting the presence of a denitrification pathway . The top database matches to these nirK and nosZ sequences belonged to marine phytoplankton epiphytes from the Alphaproteobacteriaand Flavobacterales . In the case of the nirK metagenomic sequences, the best match was to Roseibium sp. TrichSKD4, an alphaproteobacterial species isolated from a nitrogen-fixing Trichodesmium aggregate in the Atlantic Ocean . Homologs to the napA periplasmic nitrate reductase and the nirB nitrite reductase were identified from six and five unassembled reads, respectively, and were most similar to database sequences from marine heterotrophs in the Alphaproteobacteriaand Gammaproteobacteria . We did not detect sequences supporting the presence of nitrifying bacteria or archaea, though given our limited sequence depth, this could be a function of missing data. Cyanobacterial nitrogen fixation in the global oceans is frequently aggregate-associated, as with Trichodesmium sp. colonies and rafts , filamentous heterocystous cyanobacterial colonies , or Crocosphaera watsonii associated with TEP . The green berries of the Sippewissett Salt Marsh are nitrogen-fixing macroscopic consortia of unicellular cyanobacteria , diatoms and heterotrophic bacteria. While nitrogen and carbon fixation mediated by the green berries is unlikely to play a major role in the overall marsh ecosystem due to their low abundance and patchy distribution , these consortia provide an interesting comparative system to investigate the dynamics nitrogen flux within diazotrophic cyanobacterial aggregates. More broadly, studies of coastal marine estuarine sediments have indicated complex dynamics and close spatial coupling of co-occurring nitrogen fixation and denitrification processes . In other oceanic diazotrophic cyanobacterial aggregates, similarly rapid respiratory rates create transient anoxic zones within the aggregate core , and create a heterogeneous microenvironment where both aerobic and anaerobic metabolisms co-exist. This parallel is not merely conceptual: the green berry heterotrophic bacteria were similar, both phylogenetically and in their metabolic marker genes, to those observed to colonize other marine phytoplankton aggregates. The recovery of denitrification marker gene sequences in metagenome suggests that there might be heterotrophic denitrifiers in the green berries with the potential to exploit this suboxic niche. However,flower bucket further studies including rate measurements and better genome reconstructions are needed to clarify the importance of denitrification in the berries. The presence of such a pathway for nitrogen loss in the green berries is speculative, given the fragmented metagenome and absence of activity measurements, but remains an interesting direction for future work. The existence of closely coupled nitrogen cycling within diazotrophic cyanobacterial aggregates has been explored previously, though initial studies demonstrating the association of heterocystous cyanobacterial aggregates with key bacterial species and marker genes for nitrification and denitrification measured only low to negligible rates . However, more recent investigations of in vitro rates and in situ marker gene transcription indicates that denitrification within anoxic cyanobacterial aggregates could contribute significantly to nitrogen loss in hypoxic waters, where bulk oxygen concentrations would normally inhibit such activity . In the light of modeling studies showing tight spatiotemporal coupling of nitrogen fixation and denitrification processes in the global ocean , we postulate that this coupling may be occurring at the microscale within ephemeral, aggregateassociated niches. The potential for such interactions emphasizes the need to examine biogeochemical cycles from the microbial perspective in the spatially heterogeneous marine environment.Green berries were sampled in June and July of 2010 or July and August 2014 from a single intertidal pool formed in the Little Sippewissett Salt Marsh, Falmouth, MA United States .
These aggregates, however, have also been observed at other locations throughout Little and Great Sippewissett marshes. Berries were collected from the sediment-water interface by sieving and were washed three times in 0.2 µm filter sterilized marsh water. Imaging of the berries was conducted using a Zeiss Axio IMAGER MZ epifluorescence microscope equipped with a color camera and a Zeiss LSM 710 spectral confocal scanning laser microscope . The MG-RAST 3.3 pipeline was used to classify unassembled metagenomic sequence reads as ribosomal RNA and protein coding sequences using the M5RNA and M5NR databases to provide a diversity description shown in Figure 2 . Metagenomic sequence data was mined using MG-RAST. Functions of interest were mined using the hierarchical function assignment and were assigned to GB-CYN1 when the best hits were to sequenced Chroococcales genomes. Marker genes were further investigated via profile HMM and BLASTX searches . PCR amplified SSU rRNA gene sequences were aligned to the SILVA 115 database using SINA and curated using ARB . 16S rRNA phylogenies were inferred using the GTRGAMMA rate approximation. For functional genes of interest that were found to have frameshift sequencing errors , the sequences were corrected using Framebot . The metagenomic nifH amino acid sequence was aligned to references sequences from the nifH database available from the Zehr research group using the ARB software package . ProtTest3 was used to select the fixed WAG model of amino acid evolution with an inverse gamma rate approximation for the nifH phylogeny, though similar topologies were recovered using related models. Phylogenetic reconstruction for both 16S rRNA and nifH genes were conducted as follows: maximum likelihood phylogeny was constructed with RAxML 7.2.8 with 1000 rapid bootstrap inferences and Bayesian phylogeny with MrBayes 3.3 . For Bayesian tree inferences, MCMC was run with default parameters and convergence was assessed when the reported average standard deviation of split frequencies fell below 0.01. Hidden Markov model profiles of 40 phylogenetic marker genes for Bacteria and Archaeawere used to search the green berry metagenomic peptide sequence database using the trusted cutoffs. For each of the 34 markers that have green berry hits, green berry peptide sequences were aligned with all the bacterial and archaeal references sequences by hmmalign in HMMER3 . A maximum likelihood tree was built by Fastree2 for each alignments, and the alignments and trees were examined. As a result, 33 green berry peptide sequences were selected for further analysis because they branched with Cyanobacteria with good alignments. These sequences were further filtered to include only those 29 phyeco gene families with a single copy, cyanobacterial green berry hit. Single-copied reference sequences from 126 Cyanobacteria genomes of the 29 markers and the green berry cyanobacterial metagenomic sequences were retrieved from the alignments built in the previous step and were concatenated into a large alignment. Tree topology and branch lengths were optimized by the program and SH statistics was used for branch support estimation. Whole-aggregate in vivo acetylene reduction assays were conducted on two 30 mL serum bottles of five washed green berries each. Acetylene, generated by the hydration of calcium carbide in an evacuated 150 mL serum-vial, was added to 10% of the headspace in a 30 mL serum-vial containing the berries in 5 mL of anoxic 0.2 µm filtered in situ marsh water under an 90% N2-CO2 atmosphere . Experiments were incubated on a 14 h light, 10 h dark cycle with full spectrum illumination at 30◦C for 2 days. Acetylene and ethylene were quantified using a Varian 2400 series gas chromatograph with an H2 flame ionization detector, as described previously . Rates were calculated over the 2 days incubation period with the assumption that nitrogen fixation was restricted to the dark hours of incubation, as previously described for unicellular cyanobacteria . Individual aggregates were placed on an agar plate and covered with filter-sterilized seawater collected from the marsh pool. An airstream was directed at the water so that slight ripples were visible on the surface.