Oxygen minimum zones are major sites of fixed nitrogen loss in the ocean. Recent studies have highlighted the importance of anaerobic ammonium oxidation, anammox, in pelagic nitrogen removal. Sources of ammonium for the anammox reaction, however, remain controversial, as heterotrophic denitrification and alternative anaerobic pathways of organic matter remineralization cannot account for the ammonium requirements of reported anammox rates. Here, we explore the significance of microaerobic respiration as a source of ammonium during organic matter degradation in the oxygen-deficient waters off Namibia and Peru. Experiments with additions of double-labelled oxygen revealed high aerobic activity in the upper OMZs, likely controlled by surface organic matter export. Consistently observed oxygen consumption in samples retrieved throughout the lower OMZs hints at efficient exploitation of vertically and laterally advected, oxygenated waters in this zone by aerobic microorganisms. In accordance, metagenomic and metatranscriptomic analyses identified genes encoding for aerobic terminal oxidases and demonstrated their expression by diverse microbial communities, even in virtually anoxic waters. Our results suggest that microaerobic respiration is a major mode of organic matter remineralization and source of ammonium (~45-100%) in the upper oxygen minimum zones, and reconcile hitherto observed mismatches between ammonium producing and consuming processes therein.

In Eastern Boundary Upwelling Systems nutrient-rich waters are transported to the ocean surface, fuelling high photoautotrophic primary production. Subsequent heterotrophic decomposition of the produced biomass increases the oxygen-depletion at intermediate water depths, which can result in the formation of oxygen minimum zones (OMZ). OMZs can sporadically accumulate hydrogen sulfide (H₂S), which is toxic to most multicellular organisms and has been implicated in massive fish kills. During a cruise to the OMZ off Peru in January 2009 we found a sulfidic plume in continental shelf waters, covering an area >5500 km², which contained ∼2.2×10⁴ tons of H₂S. This was the first time that H₂S was measured in the Peruvian OMZ and with ∼440 km³ the largest plume ever reported for oceanic waters. We assessed the phylogenetic and functional diversity of the inhabiting microbial community by high-throughput sequencing of DNA and RNA, while its metabolic activity was determined with rate measurements of carbon fixation and nitrogen transformation processes. The waters were dominated by several distinct γ-, δ- and ε-proteobacterial taxa associated with either sulfur oxidation or sulfate reduction. Our results suggest that these chemolithoautotrophic bacteria utilized several oxidants (oxygen, nitrate, nitrite, nitric oxide and nitrous oxide) to detoxify the sulfidic waters well below the oxic surface. The chemolithoautotrophic activity at our sampling site led to high rates of dark carbon fixation. Assuming that these chemolithoautotrophic rates were maintained throughout the sulfidic waters, they could be representing as much as ∼30% of the photoautotrophic carbon fixation. Postulated changes such as eutrophication and global warming, which lead to an expansion and intensification of OMZs, might also increase the frequency of sulfidic waters. We suggest that the chemolithoautotrophically fixed carbon may be involved in a negative feedback loop that could fuel further sulfate reduction and potentially stabilize the sulfidic OMZ waters.

Marine microbes play essential roles in global energy and nutrient cycles. A primary method of determining their diversity and distribution is through sequencing of 16S ribosomal RNA genes from environmental samples. However, the perceived community composition may vary significantly based on differences in methodology, including choice of 16S variable region(s). This study investigated the influence of 16S variable region selection (V4-V5 or V6-V8) on perceived community composition and diversity for bacteria, Archaea and chloroplasts by tag-Illumina sequencing. We used 24 samples from the photic zone of the Scotian Shelf, northwest Atlantic, collected during a spring phytoplankton bloom. Taxonomic assignment and community composition varied greatly depending on the choice of variable regions while observed patterns of beta diversity were reproducible between variable regions. V4-V5 was considered the preferred variable region for future studies based on its superior recognition of Archaea, which has received little attention in bloom dynamics. The V6-V8 region captured more of the bacterial diversity, including the abundant SAR11 clades and, to a lesser extent, that of chloroplasts. However, the magnitude of difference between variable regions for bacteria and chloroplast was less than for Archaea.

Background: Profile Hidden Markov Models (HMM) are statistical representations of protein families derived from patterns of sequence conservation in multiple alignments and have been used in identifying remote homologues with considerable success. These conservation patterns arise from fold specific signals, shared across multiple families, and function specific signals unique to the families. The availability of sequences pre-classified according to their function permits the use of negative training sequences to improve the specificity of the HMM, both by optimizing the threshold cutoff and by modifying emission probabilities to minimize the influence of fold-specific signals. A protocol to generate family specific HMMs is described that first constructs a profile HMM from an alignment of the family's sequences and then uses this model to identify sequences belonging to other classes that score above the default threshold (false positives). Ten-fold cross validation is used to optimise the discrimination threshold score for the model. The advent of fast multiple alignment methods enables the use of the profile alignments to align the true and false positive sequences, and the resulting alignments are used to modify the emission probabilities in the original model. Results: The protocol, called HMM-ModE, was validated on a set of sequences belonging to six sub-families of the AGC family of kinases. These sequences have an average sequence similarity of 63% among the group though each sub-group has a different substrate specificity. The optimisation of discrimination threshold, by using negative sequences scored against the model improves specificity in test cases from an average of 21% to 98%. Further discrimination by the HMM after modifying model probabilities using negative training sequences is provided in a few cases, the average specificity rising to 99%. Similar improvements were obtained with a sample of G-Protein coupled receptors sub-classified with respect to their substrate specificity, though the average sequence identity across the sub-families is just 20.6%. The protocol is applied in a high-throughput classification exercise on protein kinases. Conclusion: The protocol has the potential to maximise the contributions of discriminating residues to classify proteins based on their molecular function, using pre-classified positive and negative sequence training data. The high specificity of the method, and increasing availability of pre-classified sequence data holds the potential for its application in sequence annotation.

International audienceBiological nitrogen fixation (BNF) supplies nutrient-depleted oceanic surface waters with new biologically available fixed nitrogen. Diazotrophs are the only organisms that can fix dinitrogen, but the factors controlling their distribution patterns in the ocean are not well understood. In this study, the relative abundances of eight diazotrophic phylotypes in the subtropical North Atlantic Ocean were determined by quantitative PCR (qPCR) of the nifH gene using TaqMan probes. A total of 152 samples were collected at 27 stations during two GEOTRACES cruises; Lisbon, Portugal to Mindelo, Cape Verde Islands (USGT10) and Woods Hole, MA, USA via the Bermuda Time Series (BATS) to Praia, Cape Verde Islands (USGT11). Seven of the eight diazotrophic phylotypes tested were detected. These included free-living and symbiotic cyanobacteria (unicellular groups (UCYN) A, B and C, Trichodesmium, the diatom-associated cyanobacteria Rhizoselinia-Richelia and Hemiaulus-Richelia) and a γ-proteobacterium (Gamma A, AY896371). The nifH gene abundances were analyzed in the context of a large set of hydrographic parameters, macronutrient and trace metal concentrations measured in parallel with DNA samples using the PRIMER-E software. The environmental variables that most influenced the abundances and distribution of the diazotrophic phylotypes were determined. We observed a geographic segregation of diazotrophic phylotypes between east and west, with UCYN A, UCYN B and UCYN C and the Rhizosolenia-Richelia symbiont associated with the eastern North Atlantic (east of 40°W), and Trichodesmium and Gamma A detected across the basin. Hemiaulus-Richelia symbionts were primarily found in temperate waters near the North American coast. The highest diazotrophic phylotype abundance and diversity were associated with temperatures greater than 22. °C in the surface mixed layer, a high supply of iron from North African aeolian mineral dust deposition and from remineralized nutrients upwelled at the edge of the oxygen minimum zone off the northwestern coast of Afric

Iron (Fe) is an essential micronutrient for many processes in all living cells. Dissolved Fe (dFe) concentrations in the ocean are of the order of a few nM, and Fe is often a factor limiting primary production. Bioavailability of Fe in aquatic environments is believed to be primarily controlled through chelation by Fe-binding ligands. Marine microbes have evolved different mechanisms to cope with the scarcity of bioavailable dFe. Gradients in dFe concentrations and diversity of the Fe-ligand pool from coastal to open ocean waters have presumably imposed selection pressures that should be reflected in the genomes of microbial communities inhabiting the pelagic realm. We applied a hidden Markov model (HMM)-based search for proteins related to cellular iron metabolism, and in particular those involved in Fe uptake mechanisms in 164 microbial genomes belonging to diverse taxa and occupying different aquatic niches. A multivariate statistical approach demonstrated that in phototrophic organisms, there is a clear influence of the ecological niche on the diversity of Fe uptake systems. Extending the analyses to the metagenome database from the Global Ocean Sampling expedition, we demonstrated that the Fe uptake and homeostasis mechanisms differed significantly across marine niches defined by temperatures and dFe concentrations, and that this difference was linked to the distribution of microbial taxa in these niches. Using the dN/dS ratios (which signify the rate of non-synonymous mutations) of the nucleotide sequences, we identified that genes encoding for TonB, Ferritin, Ferric reductase, IdiA, ZupT, and Fe²⁺ transport proteins FeoA and FeoB were evolving at a faster rate (positive selection pressure) while genes encoding ferrisiderophore, heme and Vitamin B12 uptake systems, siderophore biosynthesis, and IsiA and IsiB were under purifying selection pressure (evolving slowly).