A Metagenomic Assembly-Based Approach to Decoding Taxa in the Dead Zone

Coastal regions of eutrophication-driven oxygen depletion are widespread and increasing in number. Also known as dead zones, these regions take their name from the deleterious effects of hypoxia (dissolved oxygen less than 2 mg/L) on shrimp, demersal fish, and other animal life. Dead zones result from nutrient enrichment of primary production, concomitant consumption by chemoorganotrophic aerobic microorganisms, and strong stratification that prevents ventilation of bottom water. One of the largest dead zones in the world occurs seasonally in the northern Gulf of Mexico (nGOM), where hypoxia can reach up to 22,000 square kilometers. To explore the underlying genomic variation and metabolic potential of microorganisms in hypoxia, we performed metagenomic and metatranscriptomic sequencing on six samples from the 2013 nGOM dead zone from both hypoxic and oxic bottom waters. Over 217 Mb of sequence was assembled into contigs of at least 3 kb with IDBA-UD, with 72 greater than 100 kb, and the largest 495 kb in length. Annotation by IMG recovered over 224 thousand genes in these contigs. Binning with tetra-ESOM and quality filtering based on relative coverage of sample-specific reads led to the recovery of 83 partial to near complete (31 over 70%) high-quality genomes. These metagenomes represent key microbial taxa previously determined to be numerically abundant from 16S rRNA data, such as Thaumarcheaota, Marine Group II Euryarchaeota, SAR406, Synechococcus spp., Actinobacteria, and Planctomycetes. Ongoing work includes the recruitment of metatranscriptomic data to binned contigs for evaluation of relative gene expression, metabolic reconstruction, and comparative genomics with related organisms elsewhere in the global oceans. These data will provide us with detailed information regarding the metabolic potential and activity of many of the key players in the nGOM dead zone.