Horizontal gene transfer in an acid mine drainage microbial community

Background

Horizontal gene transfer (HGT) has been widely identified in complete prokaryotic genomes. However, the roles of HGT among members of a microbial community and in evolution remain largely unknown. With the emergence of metagenomics, it is nontrivial to investigate such horizontal flow of genetic materials among members in a microbial community from the natural environment. Because of the lack of suitable methods for metagenomics gene transfer detection, microorganisms from a low-complexity community acid mine drainage (AMD) with near-complete genomes were used to detect possible gene transfer events and suggest the biological significance.

Results

Using the annotation of coding regions by the current tools, a phylogenetic approach, and an approximately unbiased test, we found that HGTs in AMD organisms are not rare, and we predicted 119 putative transferred genes. Among them, 14 HGT events were determined to be transfer events among the AMD members. Further analysis of the 14 transferred genes revealed that the HGT events affected the functional evolution of archaea or bacteria in AMD, and it probably shaped the community structure, such as the dominance of G-plasma in archaea in AMD through HGT.

Conclusions

Our study provides a novel insight into HGT events among microorganisms in natural communities. The interconnectedness between HGT and community evolution is essential to understand microbial community formation and development.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1720-0) contains supplementary material, which is available to authorized users.     

RIEMS: a software pipeline for sensitive and comprehensive taxonomic classification of reads from metagenomics datasets

Background

Fuelled by the advent and subsequent development of next generation sequencing technologies, metagenomics became a powerful tool for the analysis of microbial communities both scientifically and diagnostically. The biggest challenge is the extraction of relevant information from the huge sequence datasets generated for metagenomics studies. Although a plethora of tools are available, data analysis is still a bottleneck.

Results

To overcome the bottleneck of data analysis, we developed an automated computational workflow called RIEMS – Reliable Information Extraction from Metagenomic Sequence datasets. RIEMS assigns every individual read sequence within a dataset taxonomically by cascading different sequence analyses with decreasing stringency of the assignments using various software applications. After completion of the analyses, the results are summarised in a clearly structured result protocol organised taxonomically. The high accuracy and performance of RIEMS analyses were proven in comparison with other tools for metagenomics data analysis using simulated sequencing read datasets.

Conclusions

RIEMS has the potential to fill the gap that still exists with regard to data analysis for metagenomics studies. The usefulness and power of RIEMS for the analysis of genuine sequencing datasets was demonstrated with an early version of RIEMS in 2011 when it was used to detect the orthobunyavirus sequences leading to the discovery of Schmallenberg virus.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0503-6) contains supplementary material, which is available to authorized users.     

Distribution and functions of TonB-dependent transporters in marine bacteria and environments: implications for dissolved organic matter utilization

Background

Bacteria play critical roles in marine nutrient cycles by incorporating and redistributing dissolved organic matter (DOM) and inorganic nutrients in the ocean. TonB-dependent transporter (TBDT) proteins allow Gram-negative bacteria to take up scarce resources from nutrient-limiting environments as well as siderophores, heme, vitamin B12, and recently identified carbohydrates. Thus, the characterization of TBDT distribution and functions is essential to better understand the contribution TBDT to DOM assimilation and its consequences on nutrient cycling in the environment.

Methodology/Principal Findings

This study presents the distribution of encoded known and putative TBDT proteins in the genomes of microorganisms and from the Global Ocean Survey data. Using a Lek clustering algorithm and substrate specificities, the TBDT sequences were mainly classified into the following three groups: (1) DOM transporters; (2) Siderophores/Vitamins transporters; and (3) Heme/Hemophores/Iron(heme)-binding protein transporters. Diverse TBDTs were found in the genomes of oligotroph Citromicrobium bathyomarinum JL354 and Citromicrobium sp JLT1363 and were highly expressed in the stationary phase of bacterial growth. The results show that the Gammaproteobacteria and the Cytophaga-Flavobacterium-Bacteroides (CFB) group bacteria accounted for the majority of the TBDT gene pool in marine surface waters.

Conclusions/Significance

The results of this study confirm the ecological importance of TBDTs in DOM assimilation for bacteria in marine environments owing to a wide range of substrate utilization potential in the ubiquitous Gammaproteobacteria and CFB group bacteria.     

Global patterns of bacterial beta-diversity in seafloor and seawater ecosystems

Background

Marine microbial communities have been essential contributors to global biomass, nutrient cycling, and biodiversity since the early history of Earth, but so far their community distribution patterns remain unknown in most marine ecosystems.

Methodology/Principal Findings

The synthesis of 9.6 million bacterial V6-rRNA amplicons for 509 samples that span the global ocean''s surface to the deep-sea floor shows that pelagic and benthic communities greatly differ, at all taxonomic levels, and share <10% bacterial types defined at 3% sequence similarity level. Surface and deep water, coastal and open ocean, and anoxic and oxic ecosystems host distinct communities that reflect productivity, land influences and other environmental constraints such as oxygen availability. The high variability of bacterial community composition specific to vent and coastal ecosystems reflects the heterogeneity and dynamic nature of these habitats. Both pelagic and benthic bacterial community distributions correlate with surface water productivity, reflecting the coupling between both realms by particle export. Also, differences in physical mixing may play a fundamental role in the distribution patterns of marine bacteria, as benthic communities showed a higher dissimilarity with increasing distance than pelagic communities.

Conclusions/Significance

This first synthesis of global bacterial distribution across different ecosystems of the World''s oceans shows remarkable horizontal and vertical large-scale patterns in bacterial communities. This opens interesting perspectives for the definition of biogeographical biomes for bacteria of ocean waters and the seabed.     

Key genes as stress indicators in the ubiquitous diatom Skeletonema marinoi

Background

The dense phytoplankton blooms that characterize productive regions and seasons in the oceans are dominated, from high to low latitudes and from coast line to open ocean, by comparatively few, often cosmopolitan species of diatoms. These key dominant species may undergo dramatic changes due to global climate change.

Results

In order to identify molecular stress-indicators for the ubiquitous diatom species Skeletonema marinoi, we tested stress-related genes in different environmental conditions (i.e. nutrient starvation/depletion, CO₂-enrichment and combined effects of these stressors) using RT-qPCR. The data show that these stressors impact algal growth rate, inducing early aging and profound changes in expression levels of the genes of interest.

Conclusions

Most analyzed genes (e.g. antioxidant-related and aldehyde dehydrogenases) were strongly down-regulated which may indicate a strategy to avoid unnecessary over-investment in their respective proteins. By contrast, key genes were activated (e.g. HSPs, GOX) which may allow the diatom species to better cope with adverse conditions. We propose the use of this panel of genes as early bio-indicators of environmental stress factors in a changing ocean.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1574-5) contains supplementary material, which is available to authorized users.     

Distribution analysis of hydrogenases in surface waters of marine and freshwater environments

Background

Surface waters of aquatic environments have been shown to both evolve and consume hydrogen and the ocean is estimated to be the principal natural source. In some marine habitats, H₂ evolution and uptake are clearly due to biological activity, while contributions of abiotic sources must be considered in others. Until now the only known biological process involved in H₂ metabolism in marine environments is nitrogen fixation.

Principal Findings

We analyzed marine and freshwater environments for the presence and distribution of genes of all known hydrogenases, the enzymes involved in biological hydrogen turnover. The total genomes and the available marine metagenome datasets were searched for hydrogenase sequences. Furthermore, we isolated DNA from samples from the North Atlantic, Mediterranean Sea, North Sea, Baltic Sea, and two fresh water lakes and amplified and sequenced part of the gene encoding the bidirectional NAD(P)-linked hydrogenase. In 21% of all marine heterotrophic bacterial genomes from surface waters, one or several hydrogenase genes were found, with the membrane-bound H₂ uptake hydrogenase being the most widespread. A clear bias of hydrogenases to environments with terrestrial influence was found. This is exemplified by the cyanobacterial bidirectional NAD(P)-linked hydrogenase that was found in freshwater and coastal areas but not in the open ocean.

Significance

This study shows that hydrogenases are surprisingly abundant in marine environments. Due to its ecological distribution the primary function of the bidirectional NAD(P)-linked hydrogenase seems to be fermentative hydrogen evolution. Moreover, our data suggests that marine surface waters could be an interesting source of oxygen-resistant uptake hydrogenases. The respective genes occur in coastal as well as open ocean habitats and we presume that they are used as additional energy scavenging devices in otherwise nutrient limited environments. The membrane-bound H₂-evolving hydrogenases might be useful as marker for bacteria living inside of marine snow particles.     

Construction of a dairy microbial genome catalog opens new perspectives for the metagenomic analysis of dairy fermented products

Background

Microbial communities of traditional cheeses are complex and insufficiently characterized. The origin, safety and functional role in cheese making of these microbial communities are still not well understood. Metagenomic analysis of these communities by high throughput shotgun sequencing is a promising approach to characterize their genomic and functional profiles. Such analyses, however, critically depend on the availability of appropriate reference genome databases against which the sequencing reads can be aligned.

Results

We built a reference genome catalog suitable for short read metagenomic analysis using a low-cost sequencing strategy. We selected 142 bacteria isolated from dairy products belonging to 137 different species and 67 genera, and succeeded to reconstruct the draft genome of 117 of them at a standard or high quality level, including isolates from the genera Kluyvera, Luteococcus and Marinilactibacillus, still missing from public database. To demonstrate the potential of this catalog, we analysed the microbial composition of the surface of two smear cheeses and one blue-veined cheese, and showed that a significant part of the microbiota of these traditional cheeses was composed of microorganisms newly sequenced in our study.

Conclusions

Our study provides data, which combined with publicly available genome references, represents the most expansive catalog to date of cheese-associated bacteria. Using this extended dairy catalog, we revealed the presence in traditional cheese of dominant microorganisms not deliberately inoculated, mainly Gram-negative genera such as Pseudoalteromonas haloplanktis or Psychrobacter immobilis, that may contribute to the characteristics of cheese produced through traditional methods.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1101) contains supplementary material, which is available to authorized users.     

Biodiversity's Big Wet Secret: The Global Distribution of Marine Biological Records Reveals Chronic Under-Exploration of the Deep Pelagic Ocean

Background

Understanding the distribution of marine biodiversity is a crucial first step towards the effective and sustainable management of marine ecosystems. Recent efforts to collate location records from marine surveys enable us to assemble a global picture of recorded marine biodiversity. They also effectively highlight gaps in our knowledge of particular marine regions. In particular, the deep pelagic ocean – the largest biome on Earth – is chronically under-represented in global databases of marine biodiversity.

Methodology/Principal Findings

We use data from the Ocean Biogeographic Information System to plot the position in the water column of ca 7 million records of marine species occurrences. Records from relatively shallow waters dominate this global picture of recorded marine biodiversity. In addition, standardising the number of records from regions of the ocean differing in depth reveals that regardless of ocean depth, most records come either from surface waters or the sea bed. Midwater biodiversity is drastically under-represented.

Conclusions/Significance

The deep pelagic ocean is the largest habitat by volume on Earth, yet it remains biodiversity''s big wet secret, as it is hugely under-represented in global databases of marine biological records. Given both its value in the provision of a range of ecosystem services, and its vulnerability to threats including overfishing and climate change, there is a pressing need to increase our knowledge of Earth''s largest ecosystem.     

Climate Change and the Potential Spreading of Marine Mucilage and Microbial Pathogens in the Mediterranean Sea

Background

Marine snow (small amorphous aggregates with colloidal properties) is present in all oceans of the world. Surface water warming and the consequent increase of water column stability can favour the coalescence of marine snow into marine mucilage, large marine aggregates representing an ephemeral and extreme habitat. Marine mucilage characterize aquatic systems with altered environmental conditions.

Methodology/Principal Findings

We investigated, by means of molecular techniques, viruses and prokaryotes within the mucilage and in surrounding seawater to examine the potential of mucilage to host new microbial diversity and/or spread marine diseases. We found that marine mucilage contained a large and unexpectedly exclusive microbial biodiversity and hosted pathogenic species that were absent in surrounding seawater. We also investigated the relationship between climate change and the frequency of mucilage in the Mediterranean Sea over the last 200 years and found that the number of mucilage outbreaks increased almost exponentially in the last 20 years. The increasing frequency of mucilage outbreaks is closely associated with the temperature anomalies.

Conclusions/Significance

We conclude that the spreading of mucilage in the Mediterranean Sea is linked to climate-driven sea surface warming. The mucilage can act as a controlling factor of microbial diversity across wide oceanic regions and could have the potential to act as a carrier of specific microorganisms, thereby increasing the spread of pathogenic bacteria.     

Genome diversification within a clonal population of pandemic Vibrio parahaemolyticus seems to depend on the life circumstances of each individual bacteria

Background

New strains of Vibrio parahaemolyticus that cause diarrhea in humans by seafood ingestion periodically emerge through continuous evolution in the ocean. Influx and expansion in the Southern Chilean ocean of a highly clonal V. parahaemolyticus (serotype O3:K6) population from South East Asia caused one of the largest seafood-related diarrhea outbreaks in the world. Here, genomics analyses of isolates from this rapidly expanding clonal population offered an opportunity to observe the molecular evolutionary changes often obscured in more diverse populations.

Results

Whole genome sequence comparison of eight independent isolates of this population from mussels or clinical cases (from different years) was performed. Differences of 1366 to 217,729 bp genome length and 13 to 164 bp single nucleotide variants (SNVs) were found. Most genomic differences corresponded to the presence of regions unique to only one or two isolates, and were probably acquired by horizontal gene transfer (HGT). Some DNA gain was chromosomal but most was in plasmids. One isolate had a large region (8,644 bp) missing, which was probably caused by excision of a prophage. Genome innovation by the presence of unique DNA, attributable to HGT from related bacteria, varied greatly among the isolates, with values of 1,366 (ten times the number of highest number of SNVs) to 217,729 (a thousand times more than the number of highest number of SNVs).

Conclusions

The evolutionary forces (SNVs, HGT) acting on each isolate of the same population were found to differ to an extent that probably depended on the ecological scenario and life circumstances of each bacterium.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1385-8) contains supplementary material, which is available to authorized users.     

A low-density SNP array for analyzing differential selection in freshwater and marine populations of threespine stickleback (Gasterosteus aculeatus)

Background

The threespine stickleback (Gasterosteus aculeatus) has become an important model species for studying both contemporary and parallel evolution. In particular, differential adaptation to freshwater and marine environments has led to high differentiation between freshwater and marine stickleback populations at the phenotypic trait of lateral plate morphology and the underlying candidate gene Ectodysplacin (EDA). Many studies have focused on this trait and candidate gene, although other genes involved in marine-freshwater adaptation may be equally important. In order to develop a resource for rapid and cost efficient analysis of genetic divergence between freshwater and marine sticklebacks, we generated a low-density SNP (Single Nucleotide Polymorphism) array encompassing markers of chromosome regions under putative directional selection, along with neutral markers for background.

Results

RAD (Restriction site Associated DNA) sequencing of sixty individuals representing two freshwater and one marine population led to the identification of 33,993 SNP markers. Ninety-six of these were chosen for the low-density SNP array, among which 70 represented SNPs under putatively directional selection in freshwater vs. marine environments, whereas 26 SNPs were assumed to be neutral. Annotation of these regions revealed several genes that are candidates for affecting stickleback phenotypic variation, some of which have been observed in previous studies whereas others are new.

Conclusions

We have developed a cost-efficient low-density SNP array that allows for rapid screening of polymorphisms in threespine stickleback. The array provides a valuable tool for analyzing adaptive divergence between freshwater and marine stickleback populations beyond the well-established candidate gene Ectodysplacin (EDA).

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-867) contains supplementary material, which is available to authorized users.     

Genome mining for ribosomally synthesized and post-translationally modified peptides (RiPPs) in anaerobic bacteria

Background

Ribosomally synthesized and post-translationally modified peptides (RiPPs) are a diverse group of biologically active bacterial molecules. Due to the conserved genomic arrangement of many of the genes involved in their synthesis, these secondary metabolite biosynthetic pathways can be predicted from genome sequence data. To date, however, despite the myriad of sequenced genomes covering many branches of the bacterial phylogenetic tree, such an analysis for a broader group of bacteria like anaerobes has not been attempted.

Results

We investigated a collection of 211 complete and published genomes, focusing on anaerobic bacteria, whose potential to encode RiPPs is relatively unknown. We showed that the presence of RiPP-genes is widespread among anaerobic representatives of the phyla Actinobacteria, Proteobacteria and Firmicutes and that, collectively, anaerobes possess the ability to synthesize a broad variety of different RiPP classes. More than 25% of anaerobes are capable of producing RiPPs either alone or in conjunction with other secondary metabolites, such as polyketides or non-ribosomal peptides.

Conclusion

Amongst the analyzed genomes, several gene clusters encode uncharacterized RiPPs, whilst others show similarity with known RiPPs. These include a number of potential class II lanthipeptides; head-to-tail cyclized peptides and lactococcin 972-like RiPP. This study presents further evidence in support of anaerobic bacteria as an untapped natural products reservoir.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-983) contains supplementary material, which is available to authorized users.     

ANGSD: Analysis of Next Generation Sequencing Data

Background

High-throughput DNA sequencing technologies are generating vast amounts of data. Fast, flexible and memory efficient implementations are needed in order to facilitate analyses of thousands of samples simultaneously.

Results

We present a multithreaded program suite called ANGSD. This program can calculate various summary statistics, and perform association mapping and population genetic analyses utilizing the full information in next generation sequencing data by working directly on the raw sequencing data or by using genotype likelihoods.

Conclusions

The open source c/c++ program ANGSD is available at http://www.popgen.dk/angsd. The program is tested and validated on GNU/Linux systems. The program facilitates multiple input formats including BAM and imputed beagle genotype probability files. The program allow the user to choose between combinations of existing methods and can perform analysis that is not implemented elsewhere.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0356-4) contains supplementary material, which is available to authorized users.     

Prokaryotic assemblages and metagenomes in pelagic zones of the South China Sea

Background

Prokaryotic microbes, the most abundant organisms in the ocean, are remarkably diverse. Despite numerous studies of marine prokaryotes, the zonation of their communities in pelagic zones has been poorly delineated. By exploiting the persistent stratification of the South China Sea (SCS), we performed a 2-year, large spatial scale (10, 100, 1000, and 3000 m) survey, which included a pilot study in 2006 and comprehensive sampling in 2007, to investigate the biological zonation of bacteria and archaea using 16S rRNA tag and shotgun metagenome sequencing.

Results

Alphaproteobacteria dominated the bacterial community in the surface SCS, where the abundance of Betaproteobacteria was seemingly associated with climatic activity. Gammaproteobacteria thrived in the deep SCS, where a noticeable amount of Cyanobacteria were also detected. Marine Groups II and III Euryarchaeota were predominant in the archaeal communities in the surface and deep SCS, respectively. Bacterial diversity was higher than archaeal diversity at all sampling depths in the SCS, and peaked at mid-depths, agreeing with the diversity pattern found in global water columns. Metagenomic analysis not only showed differential %GC values and genome sizes between the surface and deep SCS, but also demonstrated depth-dependent metabolic potentials, such as cobalamin biosynthesis at 10 m, osmoregulation at 100 m, signal transduction at 1000 m, and plasmid and phage replication at 3000 m. When compared with other oceans, urease at 10 m and both exonuclease and permease at 3000 m were more abundant in the SCS. Finally, enriched genes associated with nutrient assimilation in the sea surface and transposase in the deep-sea metagenomes exemplified the functional zonation in global oceans.

Conclusions

Prokaryotic communities in the SCS stratified with depth, with maximal bacterial diversity at mid-depth, in accordance with global water columns. The SCS had functional zonation among depths and endemically enriched metabolic potentials at the study site, in contrast to other oceans.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1434-3) contains supplementary material, which is available to authorized users.     

Metabolic reprogramming by viruses in the sunlit and dark ocean

Background

Marine ecosystem function is largely determined by matter and energy transformations mediated by microbial community interaction networks. Viral infection modulates network properties through mortality, gene transfer and metabolic reprogramming.

Results

Here we explore the nature and extent of viral metabolic reprogramming throughout the Pacific Ocean depth continuum. We describe 35 marine viral gene families with potential to reprogram metabolic flux through central metabolic pathways recovered from Pacific Ocean waters. Four of these families have been previously reported but 31 are novel. These known and new carbon pathway auxiliary metabolic genes were recovered from a total of 22 viral metagenomes in which viral auxiliary metabolic genes were differentiated from low-level cellular DNA inputs based on small subunit ribosomal RNA gene content, taxonomy, fragment recruitment and genomic context information. Auxiliary metabolic gene distribution patterns reveal that marine viruses target overlapping, but relatively distinct pathways in sunlit and dark ocean waters to redirect host carbon flux towards energy production and viral genome replication under low nutrient, niche-differentiated conditions throughout the depth continuum.

Conclusions

Given half of ocean microbes are infected by viruses at any given time, these findings of broad viral metabolic reprogramming suggest the need for renewed consideration of viruses in global ocean carbon models.     

Acidobacteria appear to dominate the microbiome of two sympatric Caribbean Sponges and one Zoanthid

Background

Marine invertebrate-associated microbial communities are interesting examples of complex symbiotic systems and are a potential source of biotechnological products.

Results

In this work, pyrosequencing-based assessment from bacterial community structures of sediments, two sponges, and one zoanthid collected in the Mexican Caribbean was performed. The results suggest that the bacterial diversity at the species level is higher in the sediments than in the animal samples. Analysis of bacterial communities’ structure showed that about two thirds of the bacterial diversity in all the samples belongs to the phyla Acidobacteria and Proteobacteria. The genus Acidobacterium appears to dominate the bacterial community in all the samples, reaching almost 80% in the sponge Hyrtios.

Conclusions

Our evidence suggests that the sympatric location of these benthonic species may lead to common bacterial structure features among their bacterial communities. The results may serve as a first insight to formulate hypotheses that lead to more extensive studies of sessile marine organisms’ microbiomes from the Mexican Caribbean.

Electronic supplementary material

The online version of this article (doi:10.1186/0717-6287-47-67) contains supplementary material, which is available to authorized users.