A method for studying protistan diversity using massively parallel sequencing of V9 hypervariable regions of small-subunit ribosomal RNA genes

© 2009 The Authors. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in PLoS ONE 4 (2009): e6372, doi:10.1371/journal.pone.0006372.Massively parallel pyrosequencing of amplicons from the V6 hypervariable regions of small-subunit (SSU) ribosomal RNA (rRNA) genes is commonly used to assess diversity and richness in bacterial and archaeal populations. Recent advances in pyrosequencing technology provide read lengths of up to 240 nucleotides. Amplicon pyrosequencing can now be applied to longer variable regions of the SSU rRNA gene including the V9 region in eukaryotes. We present a protocol for the amplicon pyrosequencing of V9 regions for eukaryotic environmental samples for biodiversity inventories and species richness estimation. The International Census of Marine Microbes (ICoMM) and the Microbial Inventory Research Across Diverse Aquatic Long Term Ecological Research Sites (MIRADA-LTERs) projects are already employing this protocol for tag sequencing of eukaryotic samples in a wide diversity of both marine and freshwater environments. Massively parallel pyrosequencing of eukaryotic V9 hypervariable regions of SSU rRNA genes provides a means of estimating species richness from deeply-sampled populations and for discovering novel species from the environment.This work was supported by grants from the W.M. Keck Foundation and the Woods Hole Center for Oceans and Human Health from the National Institutes of Health and National Science Foundation (NIH/NIEHS 1 P50 ES012742-01 and NSF/OCE 0430724-J) (LAZ and SH)

Differences in osmotolerance in freshwater and brackish water populations of Theodoxus fluviatilis (Gastropoda: Neritidae) are associated with differential protein expression

The euryhaline gastropod Theodoxus fluviatilis is found in northern Germany in freshwater or in brackish water habitats in the Baltic Sea. Previous studies have revealed that individuals from both habitats are not distinguishable by morphological characters or by sequence comparison of DNA encoding 16S RNA or cytochrome C. As reported in this study, animals collected in the two habitats differ substantially in their physiological ability to adapt to different salinities. Comparison of accumulation rates of ninhydrin-positive substances (NPS) in foot muscle upon transfer of animals to higher medium salinities revealed that brackish water animals were perfectly able to mobilize NPS, while freshwater animals had only limited ability to do so. In an attempt to explore whether this difference in physiology may be caused by genetic differentiation, we compared protein expression patterns of soluble foot muscle proteins using 2D gel electrophoresis and silver staining. Of the 40 consistently detected protein spots, 27 showed similar levels in protein expression in animals collected from freshwater or brackish water habitats, respectively. In 12 spots, however, protein concentration was higher in brackish water than in freshwater animals. In four of these spots, expression levels followed increases or decreases in medium salinities. In a different set of 4 of these 12 spots, protein levels were always higher in brackish water as compared to freshwater animals, regardless of their physiological situation (14 days in artificial pond water or in medium with a salinity of 16‰). The remaining 4 of the 12 spots had complex expression patterns. Protein levels of the remaining single spot were generally higher in freshwater animals than in brackish water animals. These expression patterns may indicate that freshwater and brackish water animals of T. fluviatilis belong to different locally adapted populations with subtle genetic differentiation

jMOTU and Taxonerator: Turning DNA Barcode Sequences into Annotated Operational Taxonomic Units

BACKGROUND: DNA barcoding and other DNA sequence-based techniques for investigating and estimating biodiversity require explicit methods for associating individual sequences with taxa, as it is at the taxon level that biodiversity is assessed. For many projects, the bioinformatic analyses required pose problems for laboratories whose prime expertise is not in bioinformatics. User-friendly tools are required for both clustering sequences into molecular operational taxonomic units (MOTU) and for associating these MOTU with known organismal taxonomies. RESULTS: Here we present jMOTU, a Java program for the analysis of DNA barcode datasets that uses an explicit, determinate algorithm to define MOTU. We demonstrate its usefulness for both individual specimen-based Sanger sequencing surveys and bulk-environment metagenetic surveys using long-read next-generation sequencing data. jMOTU is driven through a graphical user interface, and can analyse tens of thousands of sequences in a short time on a desktop computer. A companion program, Taxonerator, that adds traditional taxonomic annotation to MOTU, is also presented. Clustering and taxonomic annotation data are stored in a relational database, and are thus amenable to subsequent data mining and web presentation. CONCLUSIONS: jMOTU efficiently and robustly identifies the molecular taxa present in survey datasets, and Taxonerator decorates the MOTU with putative identifications. jMOTU and Taxonerator are freely available from http://www.nematodes.org/

Eukaryotic richness in the abyss: insights from pyrotag sequencing

Background: The deep sea floor is considered one of the most diverse ecosystems on Earth. Recent environmental DNA surveys based on clone libraries of rRNA genes confirm this observation and reveal a high diversity of eukaryotes present in deep-sea sediment samples. However, environmental clone-library surveys yield only a modest number of sequences with which to evaluate the diversity of abyssal eukaryotes. Methodology/Principal Findings: Here, we examined the richness of eukaryotic DNA in deep Arctic and Southern Ocean samples using massively parallel sequencing of the 18S ribosomal RNA (rRNA) V9 hypervariable region. In very small volumes of sediments, ranging from 0.35 to 0.7 g, we recovered up to 7,499 unique sequences per sample. By clustering sequences having up to 3 differences, we observed from 942 to 1756 Operational Taxonomic Units (OTUs) per sample. Taxonomic analyses of these OTUs showed that DNA of all major groups of eukaryotes is represented at the deep-sea floor. The dinoflagellates, cercozoans, ciliates, and euglenozoans predominate, contributing to 17%, 16%, 10%, and 8% of all assigned OTUs, respectively. Interestingly, many sequences represent photosynthetic taxa or are similar to those reported from the environmental surveys of surface waters. Moreover, each sample contained from 31 to 71 different metazoan OTUs despite the small sample volume collected. This indicates that a significant faction of the eukaryotic DNA sequences likely do not belong to living organisms, but represent either free, extracellular DNA or remains and resting stages of planktonic species. Conclusions/Significance: In view of our study, the deep-sea floor appears as a global DNA repository, which preserves genetic information about organisms living in the sediment, as well as in the water column above it. This information can be used for future monitoring of past and present environmental changes.French ANR Aquaparadox; ANR DeepOases; Swiss National Science Foundation [31003A-125372]; WM Keck foundationinfo:eu-repo/semantics/publishedVersio

Contrasting Microbial Community Assembly Hypotheses: A Reconciling Tale from the Río Tinto

The Río Tinto (RT) is distinguished from other acid mine drainage systems by its natural and ancient origins. Microbial life from all three domains flourishes in this ecosystem, but bacteria dominate metabolic processes that perpetuate environmental extremes. While the patchy geochemistry of the RT likely influences the dynamics of bacterial populations, demonstrating which environmental variables shape microbial diversity and unveiling the mechanisms underlying observed patterns, remain major challenges in microbial ecology whose answers rely upon detailed assessments of community structures coupled with fine-scale measurements of physico-chemical parameters.By using high-throughput environmental tag sequencing we achieved saturation of richness estimators for the first time in the RT. We found that environmental factors dictate the distribution of the most abundant taxa in this system, but stochastic niche differentiation processes, such as mutation and dispersal, also contribute to observed diversity patterns.We predict that studies providing clues to the evolutionary and ecological processes underlying microbial distributions will reconcile the ongoing debate between the Baas Becking vs. Hubbell community assembly hypotheses

PCR Primers for Metazoan Mitochondrial 12S Ribosomal DNA Sequences

Background: Assessment of the biodiversity of communities of small organisms is most readily done using PCR-based analysis of environmental samples consisting of mixtures of individuals. Known as metagenetics, this approach has transformed understanding of microbial communities and is beginning to be applied to metazoans as well. Unlike microbial studies, where analysis of the 16S ribosomal DNA sequence is standard, the best gene for metazoan metagenetics is less clear. In this study we designed a set of PCR primers for the mitochondrial 12S ribosomal DNA sequence based on 64 complete mitochondrial genomes and then tested their efficacy. Methodology/Principal Findings: A total of the 64 complete mitochondrial genome sequences representing all metazoan classes available in GenBank were downloaded using the NCBI Taxonomy Browser. Alignment of sequences was performed for the excised mitochondrial 12S ribosomal DNA sequences, and conserved regions were identified for all 64 mitochondrial genomes. These regions were used to design a primer pair that flanks a more variable region in the gene. Then all of the complete metazoan mitochondrial genomes available in NCBI’s Organelle Genome Resources database were used to determine the percentage of taxa that would likely be amplified using these primers. Results suggest that these primers will amplify target sequences for many metazoans. Conclusions/Significance: Newly designed 12S ribosomal DNA primers have considerable potential for metazoa

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

Differences in stress tolerance and brood size between a non-indigenous and an indigenous gammarid in the northern Baltic Sea

Differences in stress tolerance and reproductive traits may drive the competitive hierarchy between nonindigenous and indigenous species and turn the former ones into successful invaders. In the northern Baltic Sea, the non-indigenous Gammarus tigrinus is a recent invader of littoral ecosystems and now occupies comparable ecological niches as the indigenous G. zaddachi. In laboratory experiments on specimens collected between June and August 2009 around Tva¨rminne in southern Finland (59°500N/23°150E), the tolerances towards heat stress and hypoxia were determined for the two species using lethal time, LT50, as response variable. The brood size of the two species was also studied and some observations were made on maturation of juveniles. Gammarus tigrinus was more resistant to hypoxia and survived at higher temperatures than G. zaddachi. Brood size was also greater in G. tigrinus than in G. zaddachi and G. tigrinus matured at a smaller size and earlier than G. zaddachi. Hence, there are clear competitive advantages for the non-indigenous G. tigrinus compared to the indigenous G. zaddachi, and these may be further strengthened through ongoing environmental changes related to increased eutrophication and a warming climate in the Baltic Sea region

Comparative analysis of bones, mites, soil chemistry, nematodes and soil micro-Eukaryotes from a suspected homicide to estimate the post-mortem interval

Criminal investigations of suspected murder cases require estimating the post-mortem interval (PMI, or time after death) which is challenging for longer periods. Here we present the case of human remains found in a Swiss forest. We have used a multidisciplinary approach involving the analysis of bones, soil chemical characteristics, mites and nematodes (by microscopy) and micro-Eukaryotes (by Illumina high throughput sequencing). We analysed soil samples collected beneath the remains of the head, upper and lower body and “control” samples taken a few meters away. The PMI estimated on hair 14C-data via bomb peak radiocarbon dating gave a time range of 1 to 2 years before the finding of the remains on site. Cluster analyses for chemical constituents, nematodes, mites and micro- Eukaryotes revealed two clusters 1) head and upper body and 2) lower body and controls. From mite evidence, we conclude that the body was likely to have been brought to the site after death. However, chemical analyses, nematode community analyses and the analyses of micro-Eukaryotes indicate that decomposition took place at least partly on site. This study illustrates the usefulness of combining several lines of evidence for the study of homicide cases to better calibrate PMI inference tools