Metagenomics and metatranscriptomics reveal broadly distributed, active, novel methanotrophs in the Gulf of Mexico hypoxic zone and in the marine water column.

The northern Gulf of Mexico (nGOM) hypoxic zone is a shallow water environment where methane, a potent greenhouse gas, fluxes from sediments to bottom water and remains trapped due to summertime stratification. When the water column is destratified, an active planktonic methanotrophic community could mitigate the efflux of methane, which accumulates to high concentrations, to the atmosphere. To investigate the possibility of such a biofilter in the nGOM hypoxic zone we performed metagenome assembly, and metagenomic and metatranscriptomic read mapping. Methane monooxygenase (pmoA) was an abundant transcript, yet few canonical methanotrophs have been reported in this environment, suggesting a role for non-canonical methanotrophs. To determine the identity of these methanotrophs, we reconstructed six novel metagenome-assembled genomes (MAGs) in the Planctomycetota, Verrucomicrobiota and one putative Latescibacterota, each with at least one pmoA gene copy. Based on ribosomal protein phylogeny, closely related microbes (mostly from Tara Oceans) and isolate genomes were selected and co-analyzed with the nGOM MAGs. Gene annotation and read mapping suggested that there is a large, diverse and unrecognized community of active aerobic methanotrophs in the nGOM hypoxic zone and in the global ocean that could mitigate methane flux to the atmosphere.

Meta-analysis of salt marsh vegetation impacts and recovery: a synthesis following the Deepwater Horizon oil spill.

Marine oil spills continue to be a global issue, heightened by spill events such as the 2010 Deepwater Horizon spill in the Gulf of Mexico, the largest marine oil spill in US waters and among the largest worldwide, affecting over 1,000 km of sensitive wetland shorelines, primarily salt marshes supporting numerous ecosystem functions. To synthesize the effects of the oil spill on foundational vegetation species in the salt marsh ecosystem, Spartina alterniflora and Juncus roemerianus, we performed a meta-analysis using data from 10 studies and 255 sampling sites over seven years post-spill. We examined the hypotheses that the oil spill reduced plant cover, stem density, vegetation height, aboveground biomass, and belowground biomass, and tracked the degree of effects temporally to estimate recovery time frames. All plant metrics indicated impacts from oiling, with 20-100% maximum reductions depending on oiling level and marsh zone. Peak reductions of ~70-90% in total plant cover, total aboveground biomass, and belowground biomass were observed for heavily oiled sites at the marsh edge. Both Spartina and Juncus were impacted, with Juncus affected to a greater degree. Most plant metrics had recovery time frames of three years or longer, including multiple metrics with incomplete recovery over the duration of our data, at least seven years post-spill. Belowground biomass was particularly concerning, because it declined over time in contrast with recovery trends in most aboveground metrics, serving as a strong indicator of ongoing impact, limited recovery, and impaired resilience. We conclude that the Deepwater Horizon spill had multiyear impacts on salt marsh vegetation, with full recovery likely to exceed 10 years, particularly in heavily oiled marshes, where erosion may preclude full recovery. Vegetation impacts and delayed recovery is likely to have exerted substantial influences on ecosystem processes and associated species, especially along heavily oiled shorelines. Our synthesis affords a greater understanding of ecosystem impacts and recovery following the Deepwater Horizon oil spill, and informs environmental impact analysis, contingency planning, emergency response, damage assessment, and restoration efforts related to oil spills.

Comment on "Legacy nitrogen may prevent achievement of water quality goals in the Gulf of Mexico".

Van Meter et al (Reports, 27 April 2018, p. 427) warn that achieving nitrogen reduction goals in the Gulf of Mexico will take decades as a result of legacy nitrogen effects. We discuss limitations of the modeling approach and demonstrate that legacy effects ranging from a few years to decades are equally consistent with observations. The presented time scales for system recovery are therefore highly uncertain.

Oiling of the continental shelf and coastal marshes over eight years after the 2010 Deepwater Horizon oil spill.

We measured the temporal and spatial trajectory of oiling from the April, 2010, Deepwater Horizon oil spill in water from Louisiana's continental shelf, the estuarine waters of Barataria Bay, and in coastal marsh sediments. The concentrations of 28 target alkanes and 43 target polycyclic aromatic hydrocarbons were determined in water samples collected on 10 offshore cruises, in 19 water samples collected monthly one km offshore at 13 inshore stations in 2010 and 2013, and in 16-60 surficial marsh sediment samples collected on each of 26 trips. The concentration of total aromatics in offshore waters peaked in late summer, 2010, at 100 times above the May, 2010 values, which were already slightly contaminated. There were no differences in surface or bottom water samples. The concentration of total aromatics declined at a rate of 73% y-1 to 1/1000th of the May 2010 values by summer 2016. The concentrations inside the estuary were proportional to those one km offshore, but were 10-30% lower. The oil concentrations in sediments were initially different at 1 and 10 m distance into the marsh, but became equal after 2 years. Thus, the distinction between oiled and unoiled sites became blurred, if not non-existent then, and oiling had spread over an area wider than was visible initially. The concentrations of oil in sediments were 100-1000 times above the May 2010 values, and dropped to 10 times higher after 8 years, thereafter, demonstrating a long-term contamination by oil or oil residues that will remain for decades. The chemical signature of the oil residues offshore compared to in the marsh reflects the more aerobic offshore conditions and water-soluble tendencies of the dissolved components, whereas the anaerobic marsh sediments will retain the heavier molecular components for a long time, and have a consequential effect on the ecosystems.

Extent of the annual Gulf of Mexico hypoxic zone influences microbial community structure.

Rich geochemical datasets generated over the past 30 years have provided fine-scale resolution on the northern Gulf of Mexico (nGOM) coastal hypoxic (≤ 2 mg of O2 L-1) zone. In contrast, little is known about microbial community structure and activity in the hypoxic zone despite the implication that microbial respiration is responsible for forming low dissolved oxygen (DO) conditions. Here, we hypothesized that the extent of the hypoxic zone is a driver in determining microbial community structure, and in particular, the abundance of ammonia-oxidizing archaea (AOA). Samples collected across the shelf for two consecutive hypoxic seasons in July 2013 and 2014 were analyzed using 16S rRNA gene sequencing, oligotyping, microbial co-occurrence analysis, and quantification of thaumarchaeal 16S rRNA and archaeal ammonia-monooxygenase (amoA) genes. In 2014 Thaumarchaeota were enriched and inversely correlated with DO while Cyanobacteria, Acidimicrobiia, and Proteobacteria where more abundant in oxic samples compared to hypoxic. Oligotyping analysis of Nitrosopumilus 16S rRNA gene sequences revealed that one oligotype was significantly inversely correlated with DO in both years. Oligotyping analysis revealed single nucleotide variation among all Nitrosopumilaceae, including Nitrosopumilus 16S rRNA gene sequences, with one oligotype possibly being better adapted to hypoxia. We further provide evidence that in the hypoxic zone of both year 2013 and 2014, low DO concentrations and high Thaumarchaeota abundances influenced microbial co-occurrence patterns. Taken together, the data demonstrated that the extent of hypoxic conditions could potentially drive patterns in microbial community structure, with two years of data revealing the annual nGOM hypoxic zone to be emerging as a low DO adapted AOA hotspot.

Metagenomic Assembly and Prokaryotic Metagenome-Assembled Genome Sequences from the Northern Gulf of Mexico "Dead Zone".

Coastal regions experiencing declining dissolved oxygen are increasing in number and severity around the world. However, despite the importance of microbial metabolism in coastal hypoxia, few metagenomic surveys exist. Our data set from within the second largest human-caused hypoxic region provides opportunities to more deeply explore the microbiology of these systems.

The polychaete, Paraprionospio pinnata, is a likely vector of domoic acid to the benthic food web in the northern Gulf of Mexico.

A somewhat disparate, yet temporally cohesive, set of phytoplankton abundance, microphytobenthos, including the diatom Pseudo-nitzschia, benthic infauna, and sediment toxin data were used to develop a theory for the transfer of domoic acid (DA) from the toxic diatom to the benthos in the highly productive waters of the northern Gulf of Mexico near the Mississippi River plume. Archived samples and new data were used to test the theory that DA is likely to be incorporated into benthic consumers. High spring abundances of potentially toxic Pseudo-nitzschia diatoms were simultaneously present in the surface waters, bottom waters and on the seafloor. Examination of the gut contents of a typical deposit-feeding and suspension-feeding polychaete, Paraprionospio pinnata, during similar periods of high Pseudo-nitzschia abundance in surface water indicated consumption of the diatoms. Demersal fishes, particularly Atlantic croaker, are known to consume these polychaetes, with a potential for transfer of DA to even higher trophic levels. These findings warrant a theory to be tested with further studies about the trophic linkage of a phytoplankton toxin into the benthic food web.

Space-Time Geostatistical Assessment of Hypoxia in the Northern Gulf of Mexico.

Nearly every summer, a large hypoxic zone forms in the northern Gulf of Mexico. Research on the causes and consequences of hypoxia requires reliable estimates of hypoxic extent, which can vary at submonthly time scales due to hydro-meteorological variability. Here, we use an innovative space-time geostatistical model and data collected by multiple research organizations to estimate bottom-water dissolved oxygen (BWDO) concentrations and hypoxic area across summers from 1985 to 2016. We find that 27% of variability in BWDO is explained by deterministic trends with location, depth, and date, while correlated stochasticity accounts for 62% of observational variance within a range of 185 km and 28 days. Space-time modeling reduces uncertainty in estimated hypoxic area by 30% when compared to a spatial-only model, and results provide new insights into the temporal variability of hypoxia. For years with shelf-wide cruises in multiple months, hypoxia is most severe in July in 59% of years, 29% in August, and 12% in June. Also, midsummer cruise estimates of hypoxic area are only modestly correlated with summer-wide (June-August) average estimates ( r2 = 0.5), suggesting midsummer cruises are not necessarily reflective of seasonal hypoxic severity. Furthermore, summer-wide estimates are more strongly correlated with nutrient loading than midsummer estimates.

Declining oxygen in the global ocean and coastal waters.

Oxygen is fundamental to life. Not only is it essential for the survival of individual animals, but it regulates global cycles of major nutrients and carbon. The oxygen content of the open ocean and coastal waters has been declining for at least the past half-century, largely because of human activities that have increased global temperatures and nutrients discharged to coastal waters. These changes have accelerated consumption of oxygen by microbial respiration, reduced solubility of oxygen in water, and reduced the rate of oxygen resupply from the atmosphere to the ocean interior, with a wide range of biological and ecological consequences. Further research is needed to understand and predict long-term, global- and regional-scale oxygen changes and their effects on marine and estuarine fisheries and ecosystems.

Metabolic Roles of Uncultivated Bacterioplankton Lineages in the Northern Gulf of Mexico "Dead Zone".

Marine regions that have seasonal to long-term low dissolved oxygen (DO) concentrations, sometimes called "dead zones," are increasing in number and severity around the globe with deleterious effects on ecology and economics. One of the largest of these coastal dead zones occurs on the continental shelf of the northern Gulf of Mexico (nGOM), which results from eutrophication-enhanced bacterioplankton respiration and strong seasonal stratification. Previous research in this dead zone revealed the presence of multiple cosmopolitan bacterioplankton lineages that have eluded cultivation, and thus their metabolic roles in this ecosystem remain unknown. We used a coupled shotgun metagenomic and metatranscriptomic approach to determine the metabolic potential of Marine Group II Euryarchaeota, SAR406, and SAR202. We recovered multiple high-quality, nearly complete genomes from all three groups as well as candidate phyla usually associated with anoxic environments-Parcubacteria (OD1) and Peregrinibacteria Two additional groups with putative assignments to ACD39 and PAUC34f supplement the metabolic contributions by uncultivated taxa. Our results indicate active metabolism in all groups, including prevalent aerobic respiration, with concurrent expression of genes for nitrate reduction in SAR406 and SAR202, and dissimilatory nitrite reduction to ammonia and sulfur reduction by SAR406. We also report a variety of active heterotrophic carbon processing mechanisms, including degradation of complex carbohydrate compounds by SAR406, SAR202, ACD39, and PAUC34f. Together, these data help constrain the metabolic contributions from uncultivated groups in the nGOM during periods of low DO and suggest roles for these organisms in the breakdown of complex organic matter.IMPORTANCE Dead zones receive their name primarily from the reduction of eukaryotic macrobiota (demersal fish, shrimp, etc.) that are also key coastal fisheries. Excess nutrients contributed from anthropogenic activity such as fertilizer runoff result in algal blooms and therefore ample new carbon for aerobic microbial metabolism. Combined with strong stratification, microbial respiration reduces oxygen in shelf bottom waters to levels unfit for many animals (termed hypoxia). The nGOM shelf remains one of the largest eutrophication-driven hypoxic zones in the world, yet despite its potential as a model study system, the microbial metabolisms underlying and resulting from this phenomenon-many of which occur in bacterioplankton from poorly understood lineages-have received only preliminary study. Our work details the metabolic potential and gene expression activity for uncultivated lineages across several low DO sites in the nGOM, improving our understanding of the active biogeochemical cycling mediated by these "microbial dark matter" taxa during hypoxia.

Trends in summer bottom-water temperatures on the northern Gulf of Mexico continental shelf from 1985 to 2015.

We quantified trends in the 1985 to 2015 summer bottom-water temperature on the northern Gulf of Mexico (nGOM) continental shelf for data collected at 88 stations with depths ranging from 3 to 63 m. The analysis was supplemented with monthly data collected from 1963 to 1965 in the same area. The seasonal summer peak in average bottom-water temperature varied concurrently with air temperature, but with a 2- to 5-month lag. The summer bottom-water temperature declined gradually with depth from 30 oC at stations closest to the shore, to 20 oC at the offshore edge of the study area, and increased an average 0.051 oC y-1 between1963 and 2015. The bottom-water warming in summer for all stations was 1.9 times faster compared to the rise in local summer air temperatures, and 6.4 times faster than the concurrent increase in annual global ocean sea surface temperatures. The annual rise in average summer bottom-water temperatures on the subtropical nGOM continental shelf is comparable to the few published temperature trend estimates from colder environments. These recent changes in the heat storage on the nGOM continental shelf will affect oxygen and carbon cycling, spatial distribution of fish and shrimp, and overall species diversity.

Morphology, molecular phylogeny and azaspiracid profile of Azadinium poporum (Dinophyceae) from the Gulf of Mexico.

Azadinium poporum produces a variety of azaspiracids and consists of several ribotypes, but information on its biogeography is limited. A strain of A. poporum (GM29) was incubated from a Gulf of Mexico sediment sample. Strain GM29 was characterized by a plate pattern of po, cp, x, 4', 3a, 6″, 6C, 5S, 6‴, 2⁗, a distinct ventral pore at the junction of po and the first two apical plates, and a lack of an antapical spine, thus fitting the original description of A. poporum. The genus Azadinium has not been reported in waters of the United States of America before this study. Molecular phylogeny, based on large subunit ribosomal DNA (LSU rDNA) and internal transcribed spacer (ITS) sequences, reveals that strain GM29 is nested within the well-resolved A. poporum complex, but forms a sister clade either to ribotype B (ITS) or ribotype C (LSU). It is, therefore, designated as a new ribotype, termed as ribotype D. LSU and ITS sequences similarity among different ribotypes of A. poporum ranges from 95.4% to 98.2%, and from 97.1% to 99.2% respectively, suggesting that the LSU fragment is a better candidate for molecular discrimination. Azaspiracid profiles were analyzed using LC-MS/MS and demonstrate that strain GM29 produces predominantly AZA-2 with an amount of 45fg/cell. The results suggest that A. poporum has a wide distribution and highlights the risk potential of azaspiracid intoxication in the United States.

Archaeal enrichment in the hypoxic zone in the northern Gulf of Mexico.

Areas of low oxygen have spread exponentially over the past 40 years, and are cited as a key stressor on coastal ecosystems. The world's second largest coastal hypoxic (≤ 2 mg of O2 l(-1)) zone occurs annually in the northern Gulf of Mexico. The net effect of hypoxia is the diversion of energy flow away from higher trophic levels to microorganisms. This energy shunt is consequential to the overall productivity of hypoxic water masses and the ecosystem as a whole. In this study, water column samples were collected at 39 sites in the nGOM, 21 of which were hypoxic. Analysis of the microbial community along a hypoxic to oxic dissolved oxygen gradient revealed that the relative abundance (iTag) of Thaumarchaeota species 16S rRNA genes (> 40% of the microbial community in some hypoxic samples), the absolute abundance (quantitative polymerase chain reaction; qPCR) of Thaumarchaeota 16S rRNA genes and archaeal ammonia-monooxygenase gene copy number (qPCR) were significantly higher in hypoxic samples. Spatial interpolation of the microbial and chemical data revealed a continuous, shelfwide band of low dissolved oxygen waters that were dominated by Thaumarchaeota (and Euryarchaeota), amoA genes and high concentrations of phosphate in the nGOM, thus implicating physicochemical forcing on microbial abundance.

Cyst-motile stage relationship, morphology, ultrastructure, and molecular phylogeny of the gymnodinioid dinoflagellate Barrufeta resplendens comb. nov., formerly known as Gyrodinium resplendens, isolated from the Gulf of Mexico.

In the present study, we redescribed Gyrodinium resplendens through incubation of process bearing cysts extracted from sediment collected in the northern Gulf of Mexico. The morphology and ultrastructure of the motile stage and cyst stage were examined using light microscopy, scanning electron microscopy, and transmission electron microscopy and this revealed that the species should be transferred to the genus Barrufeta. This genus differs from other gymnodinioid genera in possessing a Smurf-cap apical structure complex (ASC) and currently encompasses only one species, Barrufeta bravensis. B. resplendens shows a Smurf-cap ASC that consists of three rows of elongated vesicles with small knobs in the middle one. B. resplendens is very similar to B. bravensis in cell morphology, but can be separated using the ultrastructure such as the shape and location of nucleus and pyrenoids, which highlights the importance of ultrastructure at inter-specific level in the genus Barrufeta. The unique cysts of B. resplendens are brown and process bearing, and have a tremic archeopyle with a zigzag margin on the dorsal side of the epicyst, and not polar as in cysts of Polykrikos. The cysts do not survive the palynological treatment used here and probably have a wide distribution. Maximum-likelihood and Bayesian inference were carried out based on partial large subunit ribosomal DNA (LSU rDNA) sequences. Molecular phylogeny supports that the genus Barrufeta is monophyletic, and that the genus Gymnodinium is polyphyletic. Our results suggest that details of the ASC together with ultrastructure are potential features to subdivide the genus Gymnodinium.

Cost-effective targeting of conservation investments to reduce the northern Gulf of Mexico hypoxic zone.

A seasonally occurring summer hypoxic (low oxygen) zone in the northern Gulf of Mexico is the second largest in the world. Reductions in nutrients from agricultural cropland in its watershed are needed to reduce the hypoxic zone size to the national policy goal of 5,000 km(2) (as a 5-y running average) set by the national Gulf of Mexico Task Force's Action Plan. We develop an integrated assessment model linking the water quality effects of cropland conservation investment decisions on the more than 550 agricultural subwatersheds that deliver nutrients into the Gulf with a hypoxic zone model. We use this integrated assessment model to identify the most cost-effective subwatersheds to target for cropland conservation investments. We consider targeting of the location (which subwatersheds to treat) and the extent of conservation investment to undertake (how much cropland within a subwatershed to treat). We use process models to simulate the dynamics of the effects of cropland conservation investments on nutrient delivery to the Gulf and use an evolutionary algorithm to solve the optimization problem. Model results suggest that by targeting cropland conservation investments to the most cost-effective location and extent of coverage, the Action Plan goal of 5,000 km(2) can be achieved at a cost of $2.7 billion annually. A large set of cost-hypoxia tradeoffs is developed, ranging from the baseline to the nontargeted adoption of the most aggressive cropland conservation investments in all subwatersheds (estimated to reduce the hypoxic zone to less than 3,000 km(2) at a cost of $5.6 billion annually).

Retrospective analysis of midsummer hypoxic area and volume in the northern Gulf of Mexico, 1985-2011.

Robust estimates of hypoxic extent (both area and volume) are important for assessing the impacts of low dissolved oxygen on aquatic ecosystems at large spatial scales. Such estimates are also important for calibrating models linking hypoxia to causal factors, such as nutrient loading and stratification, and for informing management decisions. In this study, we develop a rigorous geostatistical modeling framework to estimate the hypoxic extent in the northern Gulf of Mexico from data collected during midsummer, quasi-synoptic monitoring cruises (1985-2011). Instead of a traditional interpolation-based approach, we use a simulation-based approach that yields more robust extent estimates and quantified uncertainty. The modeling framework also makes use of covariate information (i.e., trend variables such as depth and spatial position), to reduce estimation uncertainty. Furthermore, adjustments are made to account for observational bias resulting from the use of different sampling instruments in different years. Our results suggest an increasing trend in hypoxic layer thickness (p = 0.05) from 1985 to 2011, but less than significant increases in volume (p = 0.12) and area (p = 0.42). The uncertainties in the extent estimates vary with sampling network coverage and instrument type, and generally decrease over the study period.

Climate change impacts on marine ecosystems.

In marine ecosystems, rising atmospheric CO2 and climate change are associated with concurrent shifts in temperature, circulation, stratification, nutrient input, oxygen content, and ocean acidification, with potentially wide-ranging biological effects. Population-level shifts are occurring because of physiological intolerance to new environments, altered dispersal patterns, and changes in species interactions. Together with local climate-driven invasion and extinction, these processes result in altered community structure and diversity, including possible emergence of novel ecosystems. Impacts are particularly striking for the poles and the tropics, because of the sensitivity of polar ecosystems to sea-ice retreat and poleward species migrations as well as the sensitivity of coral-algal symbiosis to minor increases in temperature. Midlatitude upwelling systems, like the California Current, exhibit strong linkages between climate and species distributions, phenology, and demography. Aggregated effects may modify energy and material flows as well as biogeochemical cycles, eventually impacting the overall ecosystem functioning and services upon which people and societies depend.

A Bayesian hierarchical modeling approach for analyzing observational data from marine ecological studies.

We introduce the Bayesian hierarchical modeling approach for analyzing observational data from marine ecological studies using a data set intended for inference on the effects of bottom-water hypoxia on macrobenthic communities in the northern Gulf of Mexico off the coast of Louisiana, USA. We illustrate (1) the process of developing a model, (2) the use of the hierarchical model results for statistical inference through innovative graphical presentation, and (3) a comparison to the conventional linear modeling approach (ANOVA). Our results indicate that the Bayesian hierarchical approach is better able to detect a "treatment" effect than classical ANOVA while avoiding several arbitrary assumptions necessary for linear models, and is also more easily interpreted when presented graphically. These results suggest that the hierarchical modeling approach is a better alternative than conventional linear models and should be considered for the analysis of observational field data from marine systems.