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, u planting gutter 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 Alphaproteobacteria and 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 Alphaproteobacteria and 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 genesequences in metagenome suggests that there might be heterotrophic denitrifiers in the green berries with the potential to exploit this suboxic niche. However, 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, aggregate associated 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, planting gutter 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 Archaea were 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. A maximum likelihood tree was built using PHYML3.0 with the LG substitution model . 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 head space 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. The water had a salinity of 3.5 psu and the temperature during measurements was either 24 or 27◦C, which was taken into account for calibration. We used Clark-type oxygen microsensors , and calibrated before and after measurements in air-saturated and anaerobic 3.5% NaCl. The micro-profiling apparatus and software Sensor Trace Suite was also provided from Unisense. Starting at the water surface, microprofiles through individual aggregates were measured in the light, and light-dark shifts. Time-lapse light dark-shift recordings included sensor tips placed directly above the aggregate ; and inserting the tip into the central core of the aggregate. Recordings of the oxygen signal were taken every second. Light sources were 65 W halogen lamps and experiments were conducted with either 170 µmol photons m−2 s −1 , or 320 µmol photons m−2 s −1 . Dark conditions were realized by switching off the lamps, removing them from the table, and carefully placing a carton box over the entire profiling setup to avoid residual light from the room. Background light intensities under the box were < 1 µmol m−2 s −1 . Theoretical limits of oxygen and DIC flux and whole aggregate O2 flux calculations were calculated from depth concentration profiles according to Ploug et al. , with a diffusion coefficient for O2 in 3.5% saline water of 2.175 × 10−5 cm2 s −1 at 24◦C and 2.3535 × 10−5 cm2 s −1 at 27◦C. Inside the aggregate, the apparent diffusivity of O2 was assumed to be 0.95 . Carbon fixation was estimated based on a photosynthetic quotient of 1.2 . The diffusion of oxygen in agar was not found to be different than in water over a wide range of salinities . Biometrics is a well-valued approach to subject identification for protecting access to a source of information and providing identification without the need for a password or other physical methods. It captures the physiological characteristics of a subject for identification and classification. Biometrics requires measurable physical or behavioral characteristics capable of satisfying the following considerations: universality, distinctiveness, collectability, circumvention, permanence, acceptability, and performance. There are different types of biometrics depending on the physiological characteristic, such as fingerprint, facial structure, iris, and voice, to the less popular ones such as EEG, ECG, deoxyribonucleic acid , gait, and ear shape.