Prospecting environmental mycobacteria: combined molecular approaches reveal unprecedented diversity

Background: Environmental mycobacteria (EM) include species commonly found in various terrestrial and aquatic environments, encompassing animal and human pathogens in addition to saprophytes. Approximately 150 EM species can be separated into fast and slow growers based on sequence and copy number differences of their 16S rRNA genes. Cultivation methods are not appropriate for diversity studies; few studies have investigated EM diversity in soil despite their importance as potential reservoirs of pathogens and their hypothesized role in masking or blocking M. bovis BCG vaccine. Methods: We report here the development, optimization and validation of molecular assays targeting the 16S rRNA gene to assess diversity and prevalence of fast and slow growing EM in representative soils from semi tropical and temperate areas. New primer sets were designed also to target uniquely slow growing mycobacteria and used with PCR-DGGE, tag-encoded Titanium amplicon pyrosequencing and quantitative PCR. Results: PCR-DGGE and pyrosequencing provided a consensus of EM diversity; for example, a high abundance of pyrosequencing reads and DGGE bands corresponded to M. moriokaense, M. colombiense and M. riyadhense. As expected pyrosequencing provided more comprehensive information; additional prevalent species included M. chlorophenolicum, M. neglectum, M. gordonae, M. aemonae. Prevalence of the total Mycobacterium genus in the soil samples ranged from 2.3×107 to 2.7×108 gene targets g−1; slow growers prevalence from 2.9×105 to 1.2×107 cells g−1. Conclusions: This combined molecular approach enabled an unprecedented qualitative and quantitative assessment of EM across soil samples. Good concordance was found between methods and the bioinformatics analysis was validated by random resampling. Sequences from most pathogenic groups associated with slow growth were identified in extenso in all soils tested with a specific assay, allowing to unmask them from the Mycobacterium whole genus, in which, as minority members, they would have remained undetected

Interpreting 16S metagenomic data without clustering to achieve sub-OTU resolution

The standard approach to analyzing 16S tag sequence data, which relies on clustering reads by sequence similarity into Operational Taxonomic Units (OTUs), underexploits the accuracy of modern sequencing technology. We present a clustering-free approach to multi-sample Illumina datasets that can identify independent bacterial subpopulations regardless of the similarity of their 16S tag sequences. Using published data from a longitudinal time-series study of human tongue microbiota, we are able to resolve within standard 97% similarity OTUs up to 20 distinct subpopulations, all ecologically distinct but with 16S tags differing by as little as 1 nucleotide (99.2% similarity). A comparative analysis of oral communities of two cohabiting individuals reveals that most such subpopulations are shared between the two communities at 100% sequence identity, and that dynamical similarity between subpopulations in one host is strongly predictive of dynamical similarity between the same subpopulations in the other host. Our method can also be applied to samples collected in cross-sectional studies and can be used with the 454 sequencing platform. We discuss how the sub-OTU resolution of our approach can provide new insight into factors shaping community assembly.Comment: Updated to match the published version. 12 pages, 5 figures + supplement. Significantly revised for clarity, references added, results not change

Microbiome profiling by Illumina sequencing of combinatorial sequence-tagged PCR products

We developed a low-cost, high-throughput microbiome profiling method that uses combinatorial sequence tags attached to PCR primers that amplify the rRNA V6 region. Amplified PCR products are sequenced using an Illumina paired-end protocol to generate millions of overlapping reads. Combinatorial sequence tagging can be used to examine hundreds of samples with far fewer primers than is required when sequence tags are incorporated at only a single end. The number of reads generated permitted saturating or near-saturating analysis of samples of the vaginal microbiome. The large number of reads al- lowed an in-depth analysis of errors, and we found that PCR-induced errors composed the vast majority of non-organism derived species variants, an ob- servation that has significant implications for sequence clustering of similar high-throughput data. We show that the short reads are sufficient to assign organisms to the genus or species level in most cases. We suggest that this method will be useful for the deep sequencing of any short nucleotide region that is taxonomically informative; these include the V3, V5 regions of the bac- terial 16S rRNA genes and the eukaryotic V9 region that is gaining popularity for sampling protist diversity.Comment: 28 pages, 13 figure

Evolutionary Analysis of Mitogenomes from Parasitic and Free-Living Flatworms

Copyright: © 2015 Solà et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. The attached file is the published version of the article

Shallow water marine sediment bacterial community shifts along a natural CO2 gradient in the Mediterranean Sea off Vulcano, Italy.

The effects of increasing atmospheric CO(2) on ocean ecosystems are a major environmental concern, as rapid shoaling of the carbonate saturation horizon is exposing vast areas of marine sediments to corrosive waters worldwide. Natural CO(2) gradients off Vulcano, Italy, have revealed profound ecosystem changes along rocky shore habitats as carbonate saturation levels decrease, but no investigations have yet been made of the sedimentary habitat. Here, we sampled the upper 2 cm of volcanic sand in three zones, ambient (median pCO(2) 419 μatm, minimum Ω(arag) 3.77), moderately CO(2)-enriched (median pCO(2) 592 μatm, minimum Ω(arag) 2.96), and highly CO(2)-enriched (median pCO(2) 1611 μatm, minimum Ω(arag) 0.35). We tested the hypothesis that increasing levels of seawater pCO(2) would cause significant shifts in sediment bacterial community composition, as shown recently in epilithic biofilms at the study site. In this study, 454 pyrosequencing of the V1 to V3 region of the 16S rRNA gene revealed a shift in community composition with increasing pCO(2). The relative abundances of most of the dominant genera were unaffected by the pCO(2) gradient, although there were significant differences for some 5 % of the genera present (viz. Georgenia, Lutibacter, Photobacterium, Acinetobacter, and Paenibacillus), and Shannon Diversity was greatest in sediments subject to long-term acidification (>100 years). Overall, this supports the view that globally increased ocean pCO(2) will be associated with changes in sediment bacterial community composition but that most of these organisms are resilient. However, further work is required to assess whether these results apply to other types of coastal sediments and whether the changes in relative abundance of bacterial taxa that we observed can significantly alter the biogeochemical functions of marine sediments

Bacterial diversity and community composition from seasurface to subseafloor

© The International Society for Microbial Ecology, 2015. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in ISME Journal 10 (2016): 979–989, doi:10.1038/ismej.2015.175.We investigated compositional relationships between bacterial communities in the water column and those in deep-sea sediment at three environmentally distinct Pacific sites (two in the Equatorial Pacific and one in the North Pacific Gyre). Through pyrosequencing of the v4–v6 hypervariable regions of the 16S ribosomal RNA gene, we characterized 450 104 pyrotags representing 29 814 operational taxonomic units (OTUs, 97% similarity). Hierarchical clustering and non-metric multidimensional scaling partition the samples into four broad groups, regardless of geographic location: a photic-zone community, a subphotic community, a shallow sedimentary community and a subseafloor sedimentary community (greater than or equal to1.5 meters below seafloor). Abundance-weighted community compositions of water-column samples exhibit a similar trend with depth at all sites, with successive epipelagic, mesopelagic, bathypelagic and abyssopelagic communities. Taxonomic richness is generally highest in the water-column O2 minimum zone and lowest in the subseafloor sediment. OTUs represented by abundant tags in the subseafloor sediment are often present but represented by few tags in the water column, and represented by moderately abundant tags in the shallow sediment. In contrast, OTUs represented by abundant tags in the water are generally absent from the subseafloor sediment. These results are consistent with (i) dispersal of marine sedimentary bacteria via the ocean, and (ii) selection of the subseafloor sedimentary community from within the community present in shallow sediment.This study was funded by the Biological Oceanography Program of the US National Science Foundation (grant OCE-0752336) and by the NSF-funded Center for Dark Energy Biosphere Investigations (grant NSF-OCE-0939564)

Poly-β-hydroxybutyrate administration during early life: effects on performance, immunity and microbial community of European sea bass yolk-sac larvae

The reliable production of marine fish larvae is one of the major bottlenecks in aquaculture due to high mortalities mainly caused by infectious diseases. To evaluate if the compound poly-β-hydroxybutyrate (PHB) might be a suitable immunoprophylactic measure in fish larviculture, its capacity to improve immunity and performance in European sea bass (Dicentrarchus labrax) yolk-sac larvae was explored. PHB was applied from mouth opening onwards to stimulate the developing larval immune system at the earliest possible point in time. Larval survival, growth, microbiota composition, gene expression profiles and disease resistance were assessed. PHB administration improved larval survival and, furthermore, altered the larva-associated microbiota composition. The bacterial challenge test using pathogenic Vibrio anguillarum revealed that the larval disease resistance was not influenced by PHB. The expression profiles of 26 genes involved e.g. in the immune response showed that PHB affected the expression of the antimicrobial peptides ferritin (fer) and dicentracin (dic), however, the response to PHB was inconsistent and weaker than previously demonstrated for sea bass post-larvae. Hence, the present study highlights the need for more research focusing on the immunostimulation of different early developmental stages for gaining a more comprehensive picture and advancing a sustainable production of high quality fry

An ultra-deep sequencing strategy to detect sub-clonal TP53 mutations in presentation chronic lymphocytic leukemia cases using multiple polymerases

Chronic lymphocytic leukaemia (CLL) is the most common clonal B-cell disorder characterized by clonal diversity, a relapsing and remitting course, and in its aggressive forms remains largely incurable. Current front-line regimes include agents such as fludarabine, which act primarily via the DNA damage response pathway. Key to this is the transcription factor p53. Mutations in the TP53 gene, altering p53 functionality, are associated with genetic instability, and are present in aggressive CLL. Furthermore, the emergence of clonal TP53 mutations in relapsed CLL, refractory to DNA-damaging therapy, suggests that accurate detection of sub-clonal TP53 mutations prior to and during treatment may be indicative of early relapse. In this study, we describe a novel deep sequencing workflow using multiple polymerases to generate sequencing libraries (MuPol-Seq), facilitating accurate detection of TP53 mutations at a frequency as low as 0.3%, in presentation CLL cases tested. As these mutations were mostly clustered within the regions of TP53 encoding DNA-binding domains, essential for DNA contact and structural architecture, they are likely to be of prognostic relevance in disease progression. The workflow described here has the potential to be implemented routinely to identify rare mutations across a range of diseases

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/

Seasonal dynamics of active SAR11 ecotypes in the oligotrophic Northwest Mediterranean Sea

A seven-year oceanographic time series in NW Mediterranean surface waters was combined with pyrosequencing of ribosomal RNA (16S rRNA) and ribosomal RNA gene copies (16S rDNA) to examine the environmental controls on SAR11 ecotype dynamics and potential activity. SAR11 diversity exhibited pronounced seasonal cycles remarkably similar to total bacterial diversity. The timing of diversity maxima was similar across narrow and broad phylogenetic clades and strongly associated with deep winter mixing. Diversity minima were associated with periods of stratification that were low in nutrients and phytoplankton biomass and characterised by intense phosphate limitation (turnover time80%) by SAR11 Ia. A partial least squares (PLS) regression model was developed that could reliably predict sequence abundances of SAR11 ecotypes (Q2=0.70) from measured environmental variables, of which mixed layer depth was quantitatively the most important. Comparison of clade-level SAR11 rRNA:rDNA signals with leucine incorporation enabled us to partially validate the use of these ratios as an in-situ activity measure. However, temporal trends in the activity of SAR11 ecotypes and their relationship to environmental variables were unclear. The strong and predictable temporal patterns observed in SAR11 sequence abundance was not linked to metabolic activity of different ecotypes at the phylogenetic and temporal resolution of our study