Groundtruthing next-gen sequencing for microbial ecology-biases and errors in community structure estimates from PCR amplicon pyrosequencing

SSUnique: Detecting Sequence Novelty in Microbiome Surveys

Extensive SSU rRNA gene sequence libraries, constructed from DNA extracts of environmental or host-associated samples, often contain many unclassified sequences, many representing organisms with novel taxonomy (taxonomic “blind spots”) and potentially unique ecology. This novelty is poorly explored in standard workflows, which narrows the breadth and discovery potential of such studies. Here we present the SSUnique analysis pipeline, which will promote the exploration of unclassified diversity in microbiome research and, importantly, enable the discovery of substantial novel taxonomic lineages through the analysis of a large variety of existing data sets.

High-throughput sequencing of small-subunit (SSU) rRNA genes has revolutionized understanding of microbial communities and facilitated investigations into ecological dynamics at unprecedented scales. Such extensive SSU rRNA gene sequence libraries, constructed from DNA extracts of environmental or host-associated samples, often contain a substantial proportion of unclassified sequences, many representing organisms with novel taxonomy (taxonomic “blind spots”) and potentially unique ecology. Indeed, these novel taxonomic lineages are associated with so-called microbial “dark matter,” which is the genomic potential of these lineages. Unfortunately, characterization beyond “unclassified” is challenging due to relatively short read lengths and large data set sizes. Here we demonstrate how mining of phylogenetically novel sequences from microbial ecosystems can be automated using SSUnique, a software pipeline that filters unclassified and/or rare operational taxonomic units (OTUs) from 16S rRNA gene sequence libraries by screening against consensus structural models for SSU rRNA. Phylogenetic position is inferred against a reference data set, and additional characterization of novel clades is also included, such as targeted probe/primer design and mining of assembled metagenomes for genomic context. We show how SSUnique reproduced a previous analysis of phylogenetic novelty from an Arctic tundra soil and demonstrate the recovery of highly novel clades from data sets associated with both the Earth Microbiome Project (EMP) and Human Microbiome Project (HMP). We anticipate that SSUnique will add to the expanding computational toolbox supporting high-throughput sequencing approaches for the study of microbial ecology and phylogeny.

IMPORTANCE Extensive SSU rRNA gene sequence libraries, constructed from DNA extracts of environmental or host-associated samples, often contain many unclassified sequences, many representing organisms with novel taxonomy (taxonomic “blind spots”) and potentially unique ecology. This novelty is poorly explored in standard workflows, which narrows the breadth and discovery potential of such studies. Here we present the SSUnique analysis pipeline, which will promote the exploration of unclassified diversity in microbiome research and, importantly, enable the discovery of substantial novel taxonomic lineages through the analysis of a large variety of existing data sets.

Imbalance of Gastrointestinal Microbiota in the Pathogenesis of Helicobacter pylori-Associated Diseases

The development of new nucleotide sequencing techniques and advanced bioinformatics tools has opened the field for studying the diversity and complexity of the gastrointestinal microbiome independent of traditional cultural methods. Owing largely to the gastric acid barrier, the human stomach was long thought to be sterile. The discovery of Helicobacter pylori, the gram-negative bacterium that infects upwards of 50% of the global population, has started a major paradigm shift in our understanding of the stomach as an ecologic niche for bacteria. Recent sequencing analysis of gastric microbiota showed that H. pylori was not alone and the interaction of H. pylori with those microorganisms might play a part in H. pylori-associated diseases such as gastric cancer. In this review, we summarize the available literature about the changes of gastrointestinal microbiota after H. pylori infection in humans and animal models, and discuss the possible underlying mechanisms including the alterations of the gastric environment, the secretion of hormones and the degree of inflammatory response. In general, information regarding the composition and function of gastrointestinal microbiome is still in its infancy, future studies are needed to elucidate whether and to what extent H. pylori infection perturbs the established microbiota. It is assumed that clarifying the role of gastrointestinal communities in H. pylori-associated diseases will provide an opportunity for translational application as a biomarker for the risk of serious H. pylori diseases and perhaps identify specific organisms for therapeutic eradication.

A framework for inferring biological communities from environmental DNA

Environmental DNA (eDNA), genetic material recovered from an environmental medium such as soil, water, or feces, reflects the membership of the ecological community present in the sampled environment. As such, eDNA is a potentially rich source of data for basic ecology, conservation, and management, because it offers the prospect of quantitatively reconstructing whole ecological communities from easily obtained samples. However, like all sampling methods, eDNA sequencing is subject to methodological limitations that can generate biased descriptions of ecological communities. Here, we demonstrate parallels between eDNA sampling and traditional sampling techniques, and use these parallels to offer a statistical structure for framing the challenges faced by eDNA and for illuminating the gaps in our current knowledge. Although the current state of knowledge on some of these steps precludes a full estimate of biomass for each taxon in a sampled eDNA community, we provide a map that illustrates potential methods for bridging these gaps. Additionally, we use an original data set to estimate the relative abundances of taxon-specific template DNA prior to PCR, given the abundance of DNA sequences recovered post-PCR-and-sequencing, a critical step in the chain of eDNA inference. While we focus on the use of eDNA samples to determine the relative abundance of taxa within a community, our approach also applies to single-taxon applications (including applications using qPCR), studies of diversity, and studies focused on occurrence. By grounding inferences about eDNA community composition in a rigorous statistical framework, and by making these inferences explicit, we hope to improve the inferential potential for the emerging field of community-level eDNA analysis.

Direct 16S rRNA-seq from bacterial communities: a PCR-independent approach to simultaneously assess microbial diversity and functional activity potential of each taxon

The analysis of environmental microbial communities has largely relied on a PCR-dependent amplification of genes entailing species identity as 16S rRNA. This approach is susceptible to biases depending on the level of primer matching in different species. Moreover, possible yet-to-discover taxa whose rRNA could differ enough from known ones would not be revealed. DNA-based methods moreover do not provide information on the actual physiological relevance of each taxon within an environment and are affected by the variable number of rRNA operons in different genomes. To overcome these drawbacks we propose an approach of direct sequencing of 16S ribosomal RNA without any primer- or PCR-dependent step. The method was tested on a microbial community developing in an anammox bioreactor sampled at different time-points. A conventional PCR-based amplicon pyrosequencing was run in parallel. The community resulting from direct rRNA sequencing was highly consistent with the known biochemical processes operative in the reactor. As direct rRNA-seq is based not only on taxon abundance but also on physiological activity, no comparison between its results and those from PCR-based approaches can be applied. The novel principle is in this respect proposed not as an alternative but rather as a complementary methodology in microbial community studies.

Activation of bile acid signaling improves metabolic phenotypes in high-fat diet-induced obese mice

The metabolic benefits induced by gastric bypass, currently the most effective treatment for morbid obesity, are associated with bile acid (BA) delivery to the distal intestine. However, mechanistic insights into BA signaling in the mediation of metabolic benefits remain an area of study. The bile diversion () mouse model, in which the gallbladder is anastomosed to the distal jejunum, was used to test the specific role of BA in the regulation of glucose and lipid homeostasis. Metabolic phenotype, including body weight and composition, glucose tolerance, energy expenditure, thermogenesis genes, total BA and BA composition in the circulation and portal vein, and gut microbiota were examined. BD improves the metabolic phenotype, which is in accord with increased circulating primary BAs and regulation of enterohormones. BD-induced hypertrophy of the proximal intestine in the absence of BA was reversed by BA oral gavage, but without influencing BD metabolic benefits. BD-enhanced energy expenditure was associated with elevated TGR5, D2, and thermogenic genes, including UCP1, PRDM16, PGC-1α, PGC-1β, and PDGFRα in epididymal white adipose tissue (WAT) and inguinal WAT, but not in brown adipose tissue. BD resulted in an altered gut microbiota profile (i.e., Firmicutes bacteria were decreased, Bacteroidetes were increased, and Akkermansia was positively correlated with higher levels of circulating primary BAs). Our study demonstrates that enhancement of BA signaling regulates glucose and lipid homeostasis, promotes thermogenesis, and modulates the gut microbiota, which collectively resulted in an improved metabolic phenotype.

Taxon interactions control the distributions of cryoconite bacteria colonizing a High Arctic ice cap

Microbial colonization of glacial ice surfaces incurs feedbacks which affect the melting rate of the ice surface. Ecosystems formed as microbe-mineral aggregates termed cryoconite locally reduce ice surface albedo and represent foci of biodiversity and biogeochemical cycling. Consequently, greater understanding the ecological processes in the formation of functional cryoconite ecosystems upon glacier surfaces is sought. Here, we present the first bacterial biogeography of an ice cap, evaluating the respective roles of dispersal, environmental and biotic filtration occurring at local scales in the assembly of cryoconite microbiota. 16S rRNA gene amplicon semiconductor sequencing of cryoconite colonizing a Svalbard ice cap coupled with digital elevation modelling of physical parameters reveals the bacterial community is dominated by a ubiquitous core of generalist taxa, with evidence for a moderate pairwise distance-decay relationship. While geographic position and melt season duration are prominent among environmental predictors of community structure, the core population of taxa appears highly influential in structuring the bacterial community. Taxon co-occurrence network analysis reveals a highly modular community structured by positive interactions with bottleneck taxa, predominantly Actinobacteria affiliated to isolates from soil humus. In contrast, the filamentous cyanobacterial taxon (assigned to Leptolyngbya/Phormidesmis pristleyi) which dominates the community and binds together granular cryoconite are poorly connected to other taxa. While our study targeted one ice cap, the prominent role of generalist core taxa with close environmental relatives across the global cryosphere indicate discrete roles for cosmopolitan Actinobacteria and Cyanobacteria as respective keystone taxa and ecosystem engineers of cryoconite ecosystems colonizing ice caps.

Reducing Salt in Raw Pork Sausages Increases Spoilage and Correlates with Reduced Bacterial Diversity

Raw sausages are perishable foodstuffs; reducing their salt content raises questions about a possible increased spoilage of these products. In this study, we evaluated the influence of salt reduction (from 2.0% to 1.5% [wt/wt]), in combination with two types of packaging (modified atmosphere [50% mix of CO2-N2] and vacuum packaging), on the onset of spoilage and on the diversity of spoilage-associated bacteria. After 21 days of storage at 8°C, spoilage was easily observed, characterized by noticeable graying of the products and the production of gas and off-odors defined as rancid, sulfurous, or sour. At least one of these types of spoilage occurred in each sample, and the global spoilage intensity was more pronounced in samples stored under modified atmosphere than under vacuum packaging and in samples with the lower salt content. Metagenetic 16S rRNA pyrosequencing revealed that vacuum-packaged samples contained a higher total bacterial richness (n = 69 operational taxonomic units [OTUs]) than samples under the other packaging condition (n = 46 OTUs). The core community was composed of 6 OTUs (Lactobacillus sakei, Lactococcus piscium, Carnobacterium divergens, Carnobacterium maltaromaticum, Serratia proteamaculans, and Brochothrix thermosphacta), whereas 13 OTUs taxonomically assigned to the Enterobacteriaceae, Enterococcaceae, and Leuconostocaceae families comprised a less-abundant subpopulation. This subdominant community was significantly more abundant when 2.0% salt and vacuum packaging were used, and this correlated with a lower degree of spoilage. Our results demonstrate that salt reduction, particularly when it is combined with CO2-enriched packaging, promotes faster spoilage of raw sausages by lowering the overall bacterial diversity (both richness and evenness).

IMPORTANCE

Our study takes place in the context of raw meat product manufacturing and is linked to a requirement for salt reduction. Health guidelines are calling for a reduction in dietary salt intake. However, salt has been used for a very long time as a hurdle technology, and salt reduction in meat products raises the question of spoilage and waste of food. The study was conceived to assess the role of sodium chloride reduction in meat products, both at the level of spoilage development and at the level of bacterial diversity, using 16S rRNA amplicon sequencing and raw pork sausage as a meat model.

On the use of high-throughput sequencing for the study of cyanobacterial diversity in Antarctic aquatic mats

The study of Antarctic cyanobacterial diversity has been mostly limited to morphological identification and traditional molecular techniques. High-throughput sequencing (HTS) allows a much better understanding of microbial distribution in the environment, but its application is hampered by several methodological and analytical challenges. In this work, we explored the use of HTS as a tool for the study of cyanobacterial diversity in Antarctic aquatic mats. Our results highlight the importance of using artificial communities to validate the parameters of the bioinformatics procedure used to analyze natural communities, since pipeline-dependent biases had a strong effect on the observed community structures. Analysis of microbial mats from five Antarctic lakes and an aquatic biofilm from the Sub-Antarctic showed that HTS is a valuable tool for the assessment of cyanobacterial diversity. The majority of the operational taxonomic units retrieved were related to filamentous taxa such as Leptolyngbya and Phormidium, which are common genera in Antarctic lacustrine microbial mats. However, other phylotypes related to different taxa such as Geitlerinema, Pseudanabaena, Synechococcus, Chamaesiphon, Calothrix, and Coleodesmium were also found. Results revealed a much higher diversity than what had been reported using traditional methods and also highlighted remarkable differences between the cyanobacterial communities of the studied lakes. The aquatic biofilm from the Sub-Antarctic had a distinct cyanobacterial community from the Antarctic lakes, which in turn displayed a salinity-dependent community structure at the phylotype level.

DAMe: a toolkit for the initial processing of datasets with PCR replicates of double-tagged amplicons for DNA metabarcoding analyses

Background

DNA metabarcoding is an approach for identifying multiple taxa in an environmental sample using specific genetic loci and taxa-specific primers. When combined with high-throughput sequencing it enables the taxonomic characterization of large numbers of samples in a relatively time- and cost-efficient manner. One recent laboratory development is the addition of 5′-nucleotide tags to both primers producing double-tagged amplicons and the use of multiple PCR replicates to filter erroneous sequences. However, there is currently no available toolkit for the straightforward analysis of datasets produced in this way.

Results

We present DAMe, a toolkit for the processing of datasets generated by double-tagged amplicons from multiple PCR replicates derived from an unlimited number of samples. Specifically, DAMe can be used to (i) sort amplicons by tag combination, (ii) evaluate PCR replicates dissimilarity, and (iii) filter sequences derived from sequencing/PCR errors, chimeras, and contamination. This is attained by calculating the following parameters: (i) sequence content similarity between the PCR replicates from each sample, (ii) reproducibility of each unique sequence across the PCR replicates, and (iii) copy number of the unique sequences in each PCR replicate. We showcase the insights that can be obtained using DAMe prior to taxonomic assignment, by applying it to two real datasets that vary in their complexity regarding number of samples, sequencing libraries, PCR replicates, and used tag combinations. Finally, we use a third mock dataset to demonstrate the impact and importance of filtering the sequences with DAMe.

Conclusions

DAMe allows the user-friendly manipulation of amplicons derived from multiple samples with PCR replicates built in a single or multiple sequencing libraries. It allows the user to: (i) collapse amplicons into unique sequences and sort them by tag combination while retaining the sample identifier and copy number information, (ii) identify sequences carrying unused tag combinations, (iii) evaluate the comparability of PCR replicates of the same sample, and (iv) filter tagged amplicons from a number of PCR replicates using parameters of minimum length, copy number, and reproducibility across the PCR replicates. This enables an efficient analysis of complex datasets, and ultimately increases the ease of handling datasets from large-scale studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s13104-016-2064-9) contains supplementary material, which is available to authorized users.

Intricacies of assessing the human microbiome in epidemiologic studies

PURPOSE

In the past decade, remarkable relationships have been documented between dysbiosis of the human microbiota and adverse health outcomes. This review seeks to highlight some of the challenges and pitfalls that may be encountered during all stages of microbiota research, from study design and sample collection, to nucleic acid extraction and sequencing, and bioinformatic and statistical analysis.

METHODS

Literature focused on human microbiota research was reviewed and summarized.

RESULTS

Although most studies have focused on surveying the composition of the microbiota, fewer have explored the causal roles of these bacteria, archaea, viruses, and fungi in affecting disease states. Microbiome research is in its relatively early years and many aspects remain challenging, including the complexity and personalized aspects of microbial communities, the influence of exogenous and often confounding factors, the need to apply fundamental principles of ecology and epidemiology, the necessity for new software tools, and the rapidly evolving genomic, technological, and analytical landscapes.

CONCLUSIONS

Incorporating human microbiome research in large epidemiologic studies will soon help us unravel the intricate relationships that we have with our microbial partners and provide interventional opportunities to improve human health.

The role of the intestinal microbiota in type 1 diabetes mellitus

Type 1 diabetes mellitus (T1DM) is a chronic immune-mediated disease with a subclinical prodromal period, characterized by selective loss of insulin-producing-β cells in the pancreatic islets of genetically susceptible individuals. The incidence of T1DM has increased several fold in most developed countries since World War II, in conjunction with other immune-mediated diseases. Rapid environmental changes and modern lifestyles are probably the driving factors that underlie this increase. These effects might be mediated by changes in the human microbiota, particularly the intestinal microbiota. Research on the gut microbiome of individuals at risk of developing T1DM and in patients with established disease is still in its infancy, but initial findings indicate that the intestinal microbiome of individuals with prediabetes or diabetes mellitus is different to that of healthy individuals. The gut microbiota in individuals with preclinical T1DM is characterized by Bacteroidetes dominating at the phylum level, a dearth of butyrate-producing bacteria, reduced bacterial and functional diversity and low community stability. However, these changes seem to emerge after the appearance of autoantibodies that are predictive of T1DM, which suggests that the intestinal microbiota might be involved in the progression from β-cell autoimmunity to clinical disease rather than in the initiation of the disease process.

Microbial eukaryotic distributions and diversity patterns in a deep-sea methane seep ecosystem

Although chemosynthetic ecosystems are known to support diverse assemblages of microorganisms, the ecological and environmental factors that structure microbial eukaryotes (heterotrophic protists and fungi) are poorly characterized. In this study, we examined the geographic, geochemical and ecological factors that influence microbial eukaryotic composition and distribution patterns within Hydrate Ridge, a methane seep ecosystem off the coast of Oregon using a combination of high-throughput 18S rRNA tag sequencing, terminal restriction fragment length polymorphism fingerprinting, and cloning and sequencing of full-length 18S rRNA genes. Microbial eukaryotic composition and diversity varied as a function of substrate (carbonate versus sediment), activity (low activity versus active seep sites), sulfide concentration, and region (North versus South Hydrate Ridge). Sulfide concentration was correlated with changes in microbial eukaryotic composition and richness. This work also revealed the influence of oxygen content in the overlying water column and water depth on microbial eukaryotic composition and diversity, and identified distinct patterns from those previously observed for bacteria, archaea and macrofauna in methane seep ecosystems. Characterizing the structure of microbial eukaryotic communities in response to environmental variability is a key step towards understanding if and how microbial eukaryotes influence seep ecosystem structure and function.

Exploring the Gastrointestinal "Nemabiome": Deep Amplicon Sequencing to Quantify the Species Composition of Parasitic Nematode Communities

Parasitic helminth infections have a considerable impact on global human health as well as animal welfare and production. Although co-infection with multiple parasite species within a host is common, there is a dearth of tools with which to study the composition of these complex parasite communities. Helminth species vary in their pathogenicity, epidemiology and drug sensitivity and the interactions that occur between co-infecting species and their hosts are poorly understood. We describe the first application of deep amplicon sequencing to study parasitic nematode communities as well as introduce the concept of the gastro-intestinal “nemabiome”. The approach is analogous to 16S rDNA deep sequencing used to explore microbial communities, but utilizes the nematode ITS-2 rDNA locus instead. Gastro-intestinal parasites of cattle were used to develop the concept, as this host has many well-defined gastro-intestinal nematode species that commonly occur as complex co-infections. Further, the availability of pure mono-parasite populations from experimentally infected cattle allowed us to prepare mock parasite communities to determine, and correct for, species representation biases in the sequence data. We demonstrate that, once these biases have been corrected, accurate relative quantitation of gastro-intestinal parasitic nematode communities in cattle fecal samples can be achieved. We have validated the accuracy of the method applied to field-samples by comparing the results of detailed morphological examination of L3 larvae populations with those of the sequencing assay. The results illustrate the insights that can be gained into the species composition of parasite communities, using grazing cattle in the mid-west USA as an example. However, both the technical approach and the concept of the ‘nemabiome’ have a wide range of potential applications in human and veterinary medicine. These include investigations of host-parasite and parasite-parasite interactions during co-infection, parasite epidemiology, parasite ecology and the response of parasite populations to both drug treatments and control programs.

Cracking the Code of Human Diseases Using Next-Generation Sequencing: Applications, Challenges, and Perspectives

Next-generation sequencing (NGS) technologies have greatly impacted on every field of molecular research mainly because they reduce costs and increase throughput of DNA sequencing. These features, together with the technology's flexibility, have opened the way to a variety of applications including the study of the molecular basis of human diseases. Several analytical approaches have been developed to selectively enrich regions of interest from the whole genome in order to identify germinal and/or somatic sequence variants and to study DNA methylation. These approaches are now widely used in research, and they are already being used in routine molecular diagnostics. However, some issues are still controversial, namely, standardization of methods, data analysis and storage, and ethical aspects. Besides providing an overview of the NGS-based approaches most frequently used to study the molecular basis of human diseases at DNA level, we discuss the principal challenges and applications of NGS in the field of human genomics.

AmpliconDuo: A Split-Sample Filtering Protocol for High-Throughput Amplicon Sequencing of Microbial Communities

High throughput sequencing (HTSeq) of small ribosomal subunit amplicons has the potential for a comprehensive characterization of microbial community compositions, down to rare species. However, the error-prone nature of the multi-step experimental process requires that the resulting raw sequences are subjected to quality control procedures. These procedures often involve an abundance cutoff for rare sequences or clustering of sequences, both of which limit genetic resolution. Here we propose a simple experimental protocol that retains the high genetic resolution granted by HTSeq methods while effectively removing many low abundance sequences that are likely due to PCR and sequencing errors. According to this protocol, we split samples and submit both halves to independent PCR and sequencing runs. The resulting sequence data is graphically and quantitatively characterized by the discordance between the two experimental branches, allowing for a quick identification of problematic samples. Further, we discard sequences that are not found in both branches (“AmpliconDuo filter”). We show that the majority of sequences removed in this way, mostly low abundance but also some higher abundance sequences, show features expected from random modifications of true sequences as introduced by PCR and sequencing errors. On the other hand, the filter retains many low abundance sequences observed in both branches and thus provides a more reliable census of the rare biosphere. We find that the AmpliconDuo filter increases biological resolution as it increases apparent community similarity between biologically similar communities, while it does not affect apparent community similarities between biologically dissimilar communities. The filter does not distort overall apparent community compositions. Finally, we quantitatively explain the effect of the AmpliconDuo filter by a simple mathematical model.

Synthesis of phylogeny and taxonomy into a comprehensive tree of life

Reconstructing the phylogenetic relationships that unite all lineages (the tree of life) is a grand challenge. The paucity of homologous character data across disparately related lineages currently renders direct phylogenetic inference untenable. To reconstruct a comprehensive tree of life, we therefore synthesized published phylogenies, together with taxonomic classifications for taxa never incorporated into a phylogeny. We present a draft tree containing 2.3 million tips-the Open Tree of Life. Realization of this tree required the assembly of two additional community resources: (i) a comprehensive global reference taxonomy and (ii) a database of published phylogenetic trees mapped to this taxonomy. Our open source framework facilitates community comment and contribution, enabling the tree to be continuously updated when new phylogenetic and taxonomic data become digitally available. Although data coverage and phylogenetic conflict across the Open Tree of Life illuminate gaps in both the underlying data available for phylogenetic reconstruction and the publication of trees as digital objects, the tree provides a compelling starting point for community contribution. This comprehensive tree will fuel fundamental research on the nature of biological diversity, ultimately providing up-to-date phylogenies for downstream applications in comparative biology, ecology, conservation biology, climate change, agriculture, and genomics.

Comparison of Major and Minor Viral SNPs Identified through Single Template Sequencing and Pyrosequencing in Acute HIV-1 Infection

Massively parallel sequencing (MPS) technologies, such as 454-pyrosequencing, allow for the identification of variants in sequence populations at lower levels than consensus sequencing and most single-template Sanger sequencing experiments. We sought to determine if the greater depth of population sampling attainable using MPS technology would allow detection of minor variants in HIV founder virus populations very early in infection in instances where Sanger sequencing detects only a single variant. We compared single nucleotide polymorphisms (SNPs) during acute HIV-1 infection from 32 subjects using both single template Sanger and 454-pyrosequencing. Pyrosequences from a median of 2400 viral templates per subject and encompassing 40% of the HIV-1 genome, were compared to a median of five individually amplified near full-length viral genomes sequenced using Sanger technology. There was no difference in the consensus nucleotide sequences over the 3.6kb compared in 84% of the subjects infected with single founders and 33% of subjects infected with multiple founder variants: among the subjects with disagreements, mismatches were found in less than 1% of the sites evaluated (of a total of nearly 117,000 sites across all subjects). The majority of the SNPs observed only in pyrosequences were present at less than 2% of the subject’s viral sequence population. These results demonstrate the utility of the Sanger approach for study of early HIV infection and provide guidance regarding the design, utility and limitations of population sequencing from variable template sources, and emphasize parameters for improving the interpretation of massively parallel sequencing data to address important questions regarding target sequence evolution.

Primer and platform effects on 16S rRNA tag sequencing

Sequencing of 16S rRNA gene tags is a popular method for profiling and comparing microbial communities. The protocols and methods used, however, vary considerably with regard to amplification primers, sequencing primers, sequencing technologies; as well as quality filtering and clustering. How results are affected by these choices, and whether data produced with different protocols can be meaningfully compared, is often unknown. Here we compare results obtained using three different amplification primer sets (targeting V4, V6–V8, and V7–V8) and two sequencing technologies (454 pyrosequencing and Illumina MiSeq) using DNA from a mock community containing a known number of species as well as complex environmental samples whose PCR-independent profiles were estimated using shotgun sequencing. We find that paired-end MiSeq reads produce higher quality data and enabled the use of more aggressive quality control parameters over 454, resulting in a higher retention rate of high quality reads for downstream data analysis. While primer choice considerably influences quantitative abundance estimations, sequencing platform has relatively minor effects when matched primers are used. Beta diversity metrics are surprisingly robust to both primer and sequencing platform biases.

A flexible and economical barcoding approach for highly multiplexed amplicon sequencing of diverse target genes

High throughput sequencing of phylogenetic and functional gene amplicons provides tremendous insight into the structure and functional potential of complex microbial communities. Here, we introduce a highly adaptable and economical PCR approach to barcoding and pooling libraries of numerous target genes. In this approach, we replace gene- and sequencing platform-specific fusion primers with general, interchangeable barcoding primers, enabling nearly limitless customized barcode-primer combinations. Compared to barcoding with long fusion primers, our multiple-target gene approach is more economical because it overall requires lower number of primers and is based on short primers with generally lower synthesis and purification costs. To highlight our approach, we pooled over 900 different small-subunit rRNA and functional gene amplicon libraries obtained from various environmental or host-associated microbial community samples into a single, paired-end Illumina MiSeq run. Although the amplicon regions ranged in size from approximately 290 to 720 bp, we found no significant systematic sequencing bias related to amplicon length or gene target. Our results indicate that this flexible multiplexing approach produces large, diverse, and high quality sets of amplicon sequence data for modern studies in microbial ecology.

Environmental Monitoring: Inferring the Diatom Index from Next-Generation Sequencing Data

Diatoms are widely used as bioindicators for the assessment of water quality in rivers and streams. Classically, the diatom biotic indices are based on the relative abundance of morphologically identified species weighted by their autoecological value. Obtaining such indices is time-consuming, costly, and requires excellent taxonomic expertise, which is not always available. Here we tested the possibility to overcome these limitations using a next-generation sequencing (NGS) approach to identify and quantify diatoms found in environmental DNA and RNA samples. We analyzed 27 river sites in the Geneva area (Switzerland), in order to compare the values of the Swiss Diatom Index (DI-CH) computed either by microscopic quantification of diatom species or directly from NGS data. Despite gaps in the reference database and variations in relative abundance of analyzed species, the diatom index shows a significant correlation between morphological and molecular data indicating similar biological quality status for the majority of sites. This proof-of-concept study demonstrates the potential of the NGS approach for identification and quantification of diatoms in environmental samples, opening new avenues toward the routine application of genetic tools for bioassessment and biomonitoring of aquatic ecosystems.

Deep sequencing of the viral phoH gene reveals temporal variation, depth-specific composition, and persistent dominance of the same viral phoH genes in the Sargasso Sea

Deep sequencing of the viral phoH gene, a host-derived auxiliary metabolic gene, was used to track viral diversity throughout the water column at the Bermuda Atlantic Time-series Study (BATS) site in the summer (September) and winter (March) of three years. Viral phoH sequences reveal differences in the viral communities throughout a depth profile and between seasons in the same year. Variation was also detected between the same seasons in subsequent years, though these differences were not as great as the summer/winter distinctions. Over 3,600 phoH operational taxonomic units (OTUs; 97% sequence identity) were identified. Despite high richness, most phoH sequences belong to a few large, common OTUs whereas the majority of the OTUs are small and rare. While many OTUs make sporadic appearances at just a few times or depths, a small number of OTUs dominate the community throughout the seasons, depths, and years.

Do Bacterial Symbionts Govern Aphid's Dropping Behavior?

Defensive symbiosis is amongst nature's most important interactions shaping the ecology and evolution of all partners involved. The pea aphid, Acyrthosiphon pisum Harris (Hemiptera: Aphididae), harbors one obligatory bacterial symbiont and up to seven different facultative symbionts, some of which are known to protect the aphid from pathogens, natural enemies, and other mortality factors. Pea aphids typically drop off the plant when a mammalian herbivore approaches it to avoid incidental predation. Here, we examined whether bacterial symbionts govern the pea aphid dropping behavior by comparing the bacterial fauna in dropping and nondropping aphids of two A. pisum populations, using two molecular techniques: high-throughput profiling of community structure using 16 S reads sequenced on the Illumina platform, and diagnostic polymerase chain reaction (PCR). We found that in addition to the obligatory symbiont, Buchnera aphidicola, the tested colonies of A. pisum harbored the facultative symbionts Serratia symbiotica, Regiella insecticola and Rickettsia, with no significant differences in infection proportions between dropping and nondropping aphids. While S. symbiotica was detected by both techniques, R. insecticola and Rickettsia could be detected only by diagnostic PCR. We therefore conclude that A. pisum's dropping behavior is not affected by its bacterial symbionts and is possibly affected by other factors.

Deconstructing the polymerase chain reaction: understanding and correcting bias associated with primer degeneracies and primer-template mismatches

The polymerase chain reaction (PCR) is sensitive to mismatches between primer and template, and mismatches can lead to inefficient amplification of targeted regions of DNA template. In PCRs in which a degenerate primer pool is employed, each primer can behave differently. Therefore, inefficiencies due to different primer melting temperatures within a degenerate primer pool, in addition to mismatches between primer binding sites and primers, can lead to a distortion of the true relative abundance of targets in the original DNA pool. A theoretical analysis indicated that a combination of primer-template and primer-amplicon interactions during PCR cycles 3–12 is potentially responsible for this distortion. To test this hypothesis, we developed a novel amplification strategy, entitled “Polymerase-exonuclease (PEX) PCR”, in which primer-template interactions and primer-amplicon interactions are separated. The PEX PCR method substantially and significantly improved the evenness of recovery of sequences from a mock community of known composition, and allowed for amplification of templates with introduced mismatches near the 3’ end of the primer annealing sites. When the PEX PCR method was applied to genomic DNA extracted from complex environmental samples, a significant shift in the observed microbial community was detected. Furthermore, the PEX PCR method provides a mechanism to identify which primers in a primer pool are annealing to target gDNA. Primer utilization patterns revealed that at high annealing temperatures in the PEX PCR method, perfect match annealing predominates, while at lower annealing temperatures, primers with up to four mismatches with templates can contribute substantially to amplification. The PEX PCR method is simple to perform, is limited to PCR mixes and a single exonuclease step which can be performed without reaction cleanup, and is recommended for reactions in which degenerate primer pools are used or when mismatches between primers and template are possible.

From benchtop to desktop: important considerations when designing amplicon sequencing workflows

Amplicon sequencing has been the method of choice in many high-throughput DNA sequencing (HTS) applications. To date there has been a heavy focus on the means by which to analyse the burgeoning amount of data afforded by HTS. In contrast, there has been a distinct lack of attention paid to considerations surrounding the importance of sample preparation and the fidelity of library generation. No amount of high-end bioinformatics can compensate for poorly prepared samples and it is therefore imperative that careful attention is given to sample preparation and library generation within workflows, especially those involving multiple PCR steps. This paper redresses this imbalance by focusing on aspects pertaining to the benchtop within typical amplicon workflows: sample screening, the target region, and library generation. Empirical data is provided to illustrate the scope of the problem. Lastly, the impact of various data analysis parameters is also investigated in the context of how the data was initially generated. It is hoped this paper may serve to highlight the importance of pre-analysis workflows in achieving meaningful, future-proof data that can be analysed appropriately. As amplicon sequencing gains traction in a variety of diagnostic applications from forensics to environmental DNA (eDNA) it is paramount workflows and analytics are both fit for purpose.

Pyrosequencing based profiling of the bacterial community in the Chilika Lake, the largest lagoon of India

Brackish water lake is the most extraordinary reservoir for bacterial community with an adaptability of tolerance to saline stress. In the present study, metagenomic approach was implemented utilising 454-pyrosequencing platform to gain deeper insights into the bacterial diversity profile of the soil sediment of Chilika Lake, Odisha, India. Metagenome contained 68,150 sequences with 31,896,430 bp and 56.79% G + C content. Metagenome sequences data are now available at NCBI under the Sequence Read Archive (SRA) database with accession no. SRX753382. Bacterial community metagenome sequences were analysed by MG-RAST server representing the presence of 16,212 species belonging to 45 different phyla. The dominating phyla were Proteobacteria, Chloroflexi, Firmicutes, Acidobacteria, Actinobacteria, Bacteroidetes and Planctomycetes. The analysis of bacterial community datasets obtained from two different saline soil sediments revealed significant differences in bacterial community composition and diversity value providing better understanding of the ecosystem dynamics of Chilika Lake.

RNA preservation agents and nucleic acid extraction method bias perceived bacterial community composition

Bias is a pervasive problem when characterizing microbial communities. An important source is the difference in lysis efficiencies of different populations, which vary depending on the extraction protocol used. To avoid such biases impacting comparisons between gene and transcript abundances in the environment, the use of one protocol that simultaneously extracts both types of nucleic acids from microbial community samples has gained popularity. However, knowledge regarding tradeoffs to combined nucleic acid extraction protocols is limited, particularly regarding yield and biases in the observed community composition. Here, we evaluated a commercially available protocol for simultaneous extraction of DNA and RNA, which we adapted for freshwater microbial community samples that were collected on filters. DNA and RNA yields were comparable to other commonly used, but independent DNA and RNA extraction protocols. RNA protection agents benefited RNA quality, but decreased DNA yields significantly. Choice of extraction protocol influenced the perceived bacterial community composition, with strong method-dependent biases observed for specific phyla such as the Verrucomicrobia. The combined DNA/RNA extraction protocol detected significantly higher levels of Verrucomicrobia than the other protocols, and those higher numbers were confirmed by microscopic analysis. Use of RNA protection agents as well as independent sequencing runs caused a significant shift in community composition as well, albeit smaller than the shift caused by using different extraction protocols. Despite methodological biases, sample origin was the strongest determinant of community composition. However, when the abundance of specific phylogenetic groups is of interest, researchers need to be aware of the biases their methods introduce. This is particularly relevant if different methods are used for DNA and RNA extraction, in addition to using RNA protection agents only for RNA samples.

The truth about metagenomics: quantifying and counteracting bias in 16S rRNA studies

BACKGROUND

Characterizing microbial communities via next-generation sequencing is subject to a number of pitfalls involving sample processing. The observed community composition can be a severe distortion of the quantities of bacteria actually present in the microbiome, hampering analysis and threatening the validity of conclusions from metagenomic studies. We introduce an experimental protocol using mock communities for quantifying and characterizing bias introduced in the sample processing pipeline. We used 80 bacterial mock communities comprised of prescribed proportions of cells from seven vaginally-relevant bacterial strains to assess the bias introduced in the sample processing pipeline. We created two additional sets of 80 mock communities by mixing prescribed quantities of DNA and PCR product to quantify the relative contribution to bias of (1) DNA extraction, (2) PCR amplification, and (3) sequencing and taxonomic classification for particular choices of protocols for each step. We developed models to predict the "true" composition of environmental samples based on the observed proportions, and applied them to a set of clinical vaginal samples from a single subject during four visits.

RESULTS

We observed that using different DNA extraction kits can produce dramatically different results but bias is introduced regardless of the choice of kit. We observed error rates from bias of over 85% in some samples, while technical variation was very low at less than 5% for most bacteria. The effects of DNA extraction and PCR amplification for our protocols were much larger than those due to sequencing and classification. The processing steps affected different bacteria in different ways, resulting in amplified and suppressed observed proportions of a community. When predictive models were applied to clinical samples from a subject, the predicted microbiome profiles were better reflections of the physiology and diagnosis of the subject at the visits than the observed community compositions.

CONCLUSIONS

Bias in 16S studies due to DNA extraction and PCR amplification will continue to require attention despite further advances in sequencing technology. Analysis of mock communities can help assess bias and facilitate the interpretation of results from environmental samples.

Ecology and exploration of the rare biosphere

The profound influence of microorganisms on human life and global biogeochemical cycles underlines the value of studying the biogeography of microorganisms, exploring microbial genomes and expanding our understanding of most microbial species on Earth: that is, those present at low relative abundance. The detection and subsequent analysis of low-abundance microbial populations—the 'rare biosphere'—have demonstrated the persistence, population dynamics, dispersion and predation of these microbial species. We discuss the ecology of rare microbial populations, and highlight molecular and computational methods for targeting taxonomic 'blind spots' within the rare biosphere of complex microbial communities.

The gut microbiota and inflammatory noncommunicable diseases: associations and potentials for gut microbiota therapies

Rapid environmental transition and modern lifestyles are likely driving changes in the biodiversity of the human gut microbiota. With clear effects on physiologic, immunologic, and metabolic processes in human health, aberrations in the gut microbiome and intestinal homeostasis have the capacity for multisystem effects. Changes in microbial composition are implicated in the increasing propensity for a broad range of inflammatory diseases, such as allergic disease, asthma, inflammatory bowel disease (IBD), obesity, and associated noncommunicable diseases (NCDs). There are also suggestive implications for neurodevelopment and mental health. These diverse multisystem influences have sparked interest in strategies that might favorably modulate the gut microbiota to reduce the risk of many NCDs. For example, specific prebiotics promote favorable intestinal colonization, and their fermented products have anti-inflammatory properties. Specific probiotics also have immunomodulatory and metabolic effects. However, when evaluated in clinical trials, the effects are variable, preliminary, or limited in magnitude. Fecal microbiota transplantation is another emerging therapy that regulates inflammation in experimental models. In human subjects it has been successfully used in cases of Clostridium difficile infection and IBD, although controlled trials are lacking for IBD. Here we discuss relationships between gut colonization and inflammatory NCDs and gut microbiota modulation strategies for their treatment and prevention.

Complementary amplicon-based genomic approaches for the study of fungal communities in humans

Recent studies highlight the importance of intestinal fungal microbiota in the development of human disease. Infants, in particular, are an important population in which to study intestinal microbiomes because microbial community structure and dynamics during this formative window of life have the potential to influence host immunity and metabolism. When compared to bacteria, much less is known about the early development of human fungal communities, owing partly to their lower abundance and the relative lack of established molecular and taxonomic tools for their study. Herein, we describe the development, validation, and use of complementary amplicon-based genomic strategies to characterize infant fungal communities and provide quantitative information about Candida, an important fungal genus with respect to intestinal colonization and human disease. Fungal communities were characterized from 11 infant fecal samples using primers that target the internal transcribed spacer (ITS) 2 locus, a region that provides taxonomic discrimination of medically relevant fungi. Each sample yielded an average of 27,553 fungal sequences and Candida albicans was the most abundant species identified by sequencing and quantitative PCR (qPCR). Low numbers of Candida krusei and Candida parapsilosis sequences were observed in several samples, but their presence was detected by species-specific qPCR in only one sample, highlighting a challenge inherent in the study of low-abundance organisms. Overall, the sequencing results revealed that infant fecal samples had fungal diversity comparable to that of bacterial communities in similar-aged infants, which correlated with the relative abundance of C. albicans. We conclude that targeted sequencing of fungal ITS2 amplicons in conjunction with qPCR analyses of specific fungi provides an informative picture of fungal community structure in the human intestinal tract. Our data suggests that the infant intestine harbors diverse fungal species and is consistent with prior culture-based analyses showing that the predominant fungus in the infant intestine is C. albicans.

Accurate multiplexing and filtering for high-throughput amplicon-sequencing

Tagging amplicons with tag sequences appended to PCR primers allow the multiplexing of numerous samples for high-throughput sequencing (HTS). This approach is routinely used in HTS-based diversity analyses, especially in microbial ecology and biomedical diagnostics. However, amplicon library preparation is subject to pervasive sample sequence cross-contaminations as a result of tag switching events referred to as mistagging. Here, we sequenced seven amplicon libraries prepared using various multiplexing designs in order to measure the magnitude of this phenomenon and its impact on diversity analyses. Up to 28.2% of the unique sequences correspond to undetectable (critical) mistags in single- or saturated double-tagging libraries. We show the advantage of multiplexing samples following Latin Square Designs in order to optimize the detection of mistags and maximize the information on their distribution across samples. We use this information in designs incorporating PCR replicates to filter the critical mistags and to recover the exact composition of mock community samples. Being parameter-free and data-driven, our approach can provide more accurate and reproducible HTS data sets, improving the reliability of their interpretations.

Changes in bacterial and eukaryotic communities during sewage decomposition in Mississippi river water

Microbial decay processes are one of the mechanisms whereby sewage contamination is reduced in the environment. This decomposition process involves a highly complex array of bacterial and eukaryotic communities from both sewage and ambient waters. However, relatively little is known about how these communities change due to mixing and subsequent decomposition of the sewage contaminant. We investigated decay of sewage in upper Mississippi River using Illumina sequencing of 16S and 18S rRNA gene hypervariable regions and qPCR for human-associated and general fecal Bacteroidales indicators. Mixtures of primary treated sewage and river water were placed in dialysis bags and incubated in situ under ambient conditions for seven days. We assessed changes in microbial community composition under two treatments in a replicated factorial design: sunlight exposure versus shaded and presence versus absence of native river microbiota. Initial diversity was higher in sewage compared to river water for 16S sequences, but the reverse was observed for 18S sequences. Both treatments significantly shifted community composition for eukaryotes and bacteria (P < 0.05). Data indicated that the presence of native river microbiota, rather than exposure to sunlight, accounted for the majority of variation between treatments for both 16S (R = 0.50; P > 0.001) and 18S (R = 0.91; P = 0.001) communities. A comparison of 16S sequence data and fecal indicator qPCR measurements indicated that the latter was a good predictor of overall bacterial community change over time (rho: 0.804-0.814, P = 0.001). These findings suggest that biotic interactions, such as predation by bacterivorous protozoa, can be critical factors in the decomposition of sewage in freshwater habitats and support the use of Bacteroidales genetic markers as indicators of fecal pollution.

Life history correlates of fecal bacterial species richness in a wild population of the blue tit Cyanistes caeruleus

Very little is known about the normal gastrointestinal flora of wild birds, or how it might affect or reflect the host's life-history traits. The aim of this study was to survey the species richness of bacteria in the feces of a wild population of blue tits Cyanistes caeruleus and to explore the relationships between bacterial species richness and various life-history traits, such as age, sex, and reproductive success. Using PCR-TGGE, 55 operational taxonomic units (OTUs) were identified in blue tit feces. DNA sequencing revealed that the 16S rRNA gene was amplified from a diverse range of bacteria, including those that shared closest homology with Bacillus licheniformis, Campylobacter lari, Pseudomonas spp., and Salmonella spp. For adults, there was a significant negative relationship between bacterial species richness and the likelihood of being detected alive the following breeding season; bacterial richness was consistent across years but declined through the breeding season; and breeding pairs had significantly more similar bacterial richness than expected by chance alone. Reduced adult survival was correlated with the presence of an OTU most closely resembling C. lari; enhanced adult survival was associated with an OTU most similar to Arthrobacter spp. For nestlings, there was no significant change in bacterial species richness between the first and second week after hatching, and nestlings sharing the same nest had significantly more similar bacterial richness. Collectively, these results provide compelling evidence that bacterial species richness was associated with several aspects of the life history of their hosts.

The role of the gut microbiome in the healthy adult status

The gut microbiome, which hosts up to 1000 bacterial species that encode about 5 million genes, perform many of the functions required for host physiology and survival. Consequently, it is also known as "our forgotten organ". The recent development of next-generation sequencing technologies has greatly improved metagenomic research. In particular, it has increased our knowledge about the microbiome and its mutually beneficial relationships with the human host. Microbial colonization begins immediately at birth. Although influenced by a variety of stimuli, namely, diet, physical activity, travel, illness, hormonal cycles and therapies, the microbiome is practically stable in healthy adults. This suggests that the microbiome plays a role in the maintenance of a healthy state in adulthood. Quantitative and qualitative alterations in the composition of the gut microbiome could lead to pathological dysbiosis, and have been related to an increasing number of intestinal and extra-intestinal diseases. With the increase in knowledge about gut microbiome functions, it is becoming increasingly more possible to develop novel diagnostic, prognostic and, most important, therapeutic strategies based on microbiome manipulation.

High temporal and spatial diversity in marine RNA viruses implies that they have an important role in mortality and structuring plankton communities

Viruses in the order Picornavirales infect eukaryotes, and are widely distributed in coastal waters. Amplicon deep-sequencing of the RNA dependent RNA polymerase (RdRp) revealed diverse and highly uneven communities of picorna-like viruses in the coastal waters of British Columbia (BC), Canada. Almost 300 000 pyrosequence reads revealed 145 operational taxonomic units (OTUs) based on 95% sequence similarity at the amino-acid level. Each sample had between 24 and 71 OTUs and there was little overlap among samples. Phylogenetic analysis revealed that some clades of OTUs were only found at one site; whereas, other clades included OTUs from all sites. Since most of these OTUs are likely from viruses that infect eukaryotic phytoplankton, and viral isolates infecting phytoplankton are strain-specific; each OTU probably arose from the lysis of a specific phytoplankton taxon. Moreover, the patchiness in OTU distribution, and the high turnover of viruses in the mixed layer, implies continuous infection and lysis by RNA viruses of a diverse array of eukaryotic phytoplankton taxa. Hence, these viruses are likely important elements structuring the phytoplankton community, and play a significant role in nutrient cycling and energy transfer.

Séance: reference-based phylogenetic analysis for 18S rRNA studies

Background

Marker gene studies often use short amplicons spanning one or more hypervariable regions from an rRNA gene to interrogate the community structure of uncultured environmental samples. Target regions are chosen for their discriminatory power, but the limited phylogenetic signal of short high-throughput sequencing reads precludes accurate phylogenetic analysis. This is particularly unfortunate in the study of microscopic eukaryotes where horizontal gene flow is limited and the rRNA gene is expected to accurately reflect the species phylogeny. A promising alternative to full phylogenetic analysis is phylogenetic placement, where a reference phylogeny is inferred using the complete marker gene and iteratively extended with the short sequences from a metagenetic sample under study.

Results

Based on the phylogenetic placement approach we built Séance, a community analysis pipeline focused on the analysis of 18S marker gene data. Séance combines the alignment extension and phylogenetic placement capabilities of the Pagan multiple sequence alignment program with a suite of tools to preprocess, cluster and visualise datasets composed of many samples. We showcase Séance by analysing 454 data from a longitudinal study of intestinal parasite communities in wild rufous mouse lemurs (Microcebus rufus) as well as in simulation. We demonstrate both improved OTU picking at higher levels of sequence similarity for 454 data and show the accuracy of phylogenetic placement to be comparable to maximum likelihood methods for lower numbers of taxa.

Conclusions

Séance is an open source community analysis pipeline that provides reference-based phylogenetic analysis for rRNA marker gene studies. Whilst in this article we focus on studying nematodes using the 18S marker gene, the concepts are generic and reference data for alternative marker genes can be easily created. Séance can be downloaded from http://wasabiapp.org/software/seance/.

Electronic supplementary material

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

Faecal microbiota of forage-fed horses in New Zealand and the population dynamics of microbial communities following dietary change

The effects of abrupt dietary transition on the faecal microbiota of forage-fed horses over a 3-week period were investigated. Yearling Thoroughbred fillies reared as a cohort were exclusively fed on either an ensiled conserved forage-grain diet ("Group A"; n = 6) or pasture ("Group B"; n = 6) for three weeks prior to the study. After the Day 0 faecal samples were collected, horses of Group A were abruptly transitioned to pasture. Both groups continued to graze similar pasture for three weeks, with faecal samples collected at 4-day intervals. DNA was isolated from the faeces and microbial 16S and 18S rRNA gene amplicons were generated and analysed by pyrosequencing. The faecal bacterial communities of both groups of horses were highly diverse (Simpson's index of diversity > 0.8), with differences between the two groups on Day 0 (P < 0.017 adjusted for multiple comparisons). There were differences between Groups A and B in the relative abundances of four genera, BF311 (family Bacteroidaceae; P = 0.003), CF231 (family Paraprevotellaceae; P = 0.004), and currently unclassified members within the order Clostridiales (P = 0.003) and within the family Lachnospiraceae (P = 0.006). The bacterial community of Group A horses became similar to Group B within four days of feeding on pasture, whereas the structure of the archaeal community remained constant pre- and post-dietary change. The community structure of the faecal microbiota (bacteria, archaea and ciliate protozoa) of pasture-fed horses was also identified. The initial differences observed appeared to be linked to recent dietary history, with the bacterial community of the forage-fed horses responding rapidly to abrupt dietary change.

16S rRNA gene survey of microbial communities in Winogradsky columns

A Winogradsky column is a clear glass or plastic column filled with enriched sediment. Over time, microbial communities in the sediment grow in a stratified ecosystem with an oxic top layer and anoxic sub-surface layers. Winogradsky columns have been used extensively to demonstrate microbial nutrient cycling and metabolic diversity in undergraduate microbiology labs. In this study, we used high-throughput 16s rRNA gene sequencing to investigate the microbial diversity of Winogradsky columns. Specifically, we tested the impact of sediment source, supplemental cellulose source, and depth within the column, on microbial community structure. We found that the Winogradsky columns were highly diverse communities but are dominated by three phyla: Proteobacteria, Bacteroidetes, and Firmicutes. The community is structured by a founding population dependent on the source of sediment used to prepare the columns and is differentiated by depth within the column. Numerous biomarkers were identified distinguishing sample depth, including Cyanobacteria, Alphaproteobacteria, and Betaproteobacteria as biomarkers of the soil-water interface, and Clostridia as a biomarker of the deepest depth. Supplemental cellulose source impacted community structure but less strongly than depth and sediment source. In columns dominated by Firmicutes, the family Peptococcaceae was the most abundant sulfate reducer, while in columns abundant in Proteobacteria, several Deltaproteobacteria families, including Desulfobacteraceae, were found, showing that different taxonomic groups carry out sulfur cycling in different columns. This study brings this historical method for enrichment culture of chemolithotrophs and other soil bacteria into the modern era of microbiology and demonstrates the potential of the Winogradsky column as a model system for investigating the effect of environmental variables on soil microbial communities.

Intracellular diversity of the V4 and V9 regions of the 18S rRNA in marine protists (radiolarians) assessed by high-throughput sequencing

Metabarcoding is a powerful tool for exploring microbial diversity in the environment, but its accurate interpretation is impeded by diverse technical (e.g. PCR and sequencing errors) and biological biases (e.g. intra-individual polymorphism) that remain poorly understood. To help interpret environmental metabarcoding datasets, we investigated the intracellular diversity of the V4 and V9 regions of the 18S rRNA gene from Acantharia and Nassellaria (radiolarians) using 454 pyrosequencing. Individual cells of radiolarians were isolated, and PCRs were performed with generalist primers to amplify the V4 and V9 regions. Different denoising procedures were employed to filter the pyrosequenced raw amplicons (Acacia, AmpliconNoise, Linkage method). For each of the six isolated cells, an average of 541 V4 and 562 V9 amplicons assigned to radiolarians were obtained, from which one numerically dominant sequence and several minor variants were found. At the 97% identity, a diversity metrics commonly used in environmental surveys, up to 5 distinct OTUs were detected in a single cell. However, most amplicons grouped within a single OTU whereas other OTUs contained very few amplicons. Different analytical methods provided evidence that most minor variants forming different OTUs correspond to PCR and sequencing artifacts. Duplicate PCR and sequencing from the same DNA extract of a single cell had only 9 to 16% of unique amplicons in common, and alignment visualization of V4 and V9 amplicons showed that most minor variants contained substitutions in highly-conserved regions. We conclude that intracellular variability of the 18S rRNA in radiolarians is very limited despite its multi-copy nature and the existence of multiple nuclei in these protists. Our study recommends some technical guidelines to conservatively discard artificial amplicons from metabarcoding datasets, and thus properly assess the diversity and richness of protists in the environment.

Richness and diversity in dust stormborne biomes at the southeast mediterranean

Dust storms include particulate matter that is transported over land and sea with biota that could impact downwind ecosystems. In addition to the physico-chemical compositions, organismal diversities of dust from two storm events in southern Israel, December 2012 (Ev12) and January 2013 (Ev13), were determined by pyro-sequencing using primers universal to 16S and 18S rRNA genes and compared. The bio-assemblages in the collected dust samples were affiliated with scores of different taxa. Distinct patterns of richness and diversity of the two events were influenced by the origins of the air masses: Ev13 was rich with reads affiliated to Betaproteobacteria and Embryophyta, consistent with a European origin. Ev12, originated in north-Africa, contained significantly more of the Actinobacteria and fungi, without conifers. The abundance of bacterial and eukaryotic reads demonstrates dissemination of biological material in dust that may impose health hazards of pathogens and allergens, and influence vegetation migration throughout the world.

The bias associated with amplicon sequencing does not affect the quantitative assessment of bacterial community dynamics

The performance of two sets of primers targeting variable regions of the 16S rRNA gene V1–V3 and V4 was compared in their ability to describe changes of bacterial diversity and temporal turnover in full-scale activated sludge. Duplicate sets of high-throughput amplicon sequencing data of the two 16S rRNA regions shared a collection of core taxa that were observed across a series of twelve monthly samples, although the relative abundance of each taxon was substantially different between regions. A case in point was the changes in the relative abundance of filamentous bacteria Thiothrix, which caused a large effect on diversity indices, but only in the V1–V3 data set. Yet the relative abundance of Thiothrix in the amplicon sequencing data from both regions correlated with the estimation of its abundance determined using fluorescence in situ hybridization. In nonmetric multidimensional analysis samples were distributed along the first ordination axis according to the sequenced region rather than according to sample identities. The dynamics of microbial communities indicated that V1–V3 and the V4 regions of the 16S rRNA gene yielded comparable patterns of: 1) the changes occurring within the communities along fixed time intervals, 2) the slow turnover of activated sludge communities and 3) the rate of species replacement calculated from the taxa–time relationships. The temperature was the only operational variable that showed significant correlation with the composition of bacterial communities over time for the sets of data obtained with both pairs of primers. In conclusion, we show that despite the bias introduced by amplicon sequencing, the variable regions V1–V3 and V4 can be confidently used for the quantitative assessment of bacterial community dynamics, and provide a proper qualitative account of general taxa in the community, especially when the data are obtained over a convenient time window rather than at a single time point.

Bacterial diversity assessment in Antarctic terrestrial and aquatic microbial mats: a comparison between bidirectional pyrosequencing and cultivation

The application of high-throughput sequencing of the 16S rRNA gene has increased the size of microbial diversity datasets by several orders of magnitude, providing improved access to the rare biosphere compared with cultivation-based approaches and more established cultivation-independent techniques. By contrast, cultivation-based approaches allow the retrieval of both common and uncommon bacteria that can grow in the conditions used and provide access to strains for biotechnological applications. We performed bidirectional pyrosequencing of the bacterial 16S rRNA gene diversity in two terrestrial and seven aquatic Antarctic microbial mat samples previously studied by heterotrophic cultivation. While, not unexpectedly, 77.5% of genera recovered by pyrosequencing were not among the isolates, 25.6% of the genera picked up by cultivation were not detected by pyrosequencing. To allow comparison between both techniques, we focused on the five phyla (Proteobacteria, Actinobacteria, Bacteroidetes, Firmicutes and Deinococcus-Thermus) recovered by heterotrophic cultivation. Four of these phyla were among the most abundantly recovered by pyrosequencing. Strikingly, there was relatively little overlap between cultivation and the forward and reverse pyrosequencing-based datasets at the genus (17.1–22.2%) and OTU (3.5–3.6%) level (defined on a 97% similarity cut-off level). Comparison of the V1–V2 and V3–V2 datasets of the 16S rRNA gene revealed remarkable differences in number of OTUs and genera recovered. The forward dataset missed 33% of the genera from the reverse dataset despite comprising 50% more OTUs, while the reverse dataset did not contain 40% of the genera of the forward dataset. Similar observations were evident when comparing the forward and reverse cultivation datasets. Our results indicate that the region under consideration can have a large impact on perceived diversity, and should be considered when comparing different datasets. Finally, a high number of OTUs could not be classified using the RDP reference database, suggesting the presence of a large amount of novel diversity.

Biofilm community dynamics in bench-scale annular reactors simulating arrestment of chloraminated drinking water nitrification

Annular reactors (ARs) were used to study biofilm community succession and provide ecological insight during nitrification arrestment through simultaneously increasing monochloramine (NH2Cl) and chlorine to nitrogen mass ratios, resulting in four operational periods (I-IV). Analysis of 16S rRNA-encoding gene sequence reads (454-pyrosequencing) examined viable and total biofilm communities and found total samples were representative of the underlying viable community. Bacterial community structure showed dynamic changes corresponding with AR operational parameters. Period I (complete nitrification and no NH2Cl residual) was dominated by Bradyrhizobium (total cumulative distribution: 38%), while environmental Legionella-like phylotypes peaked (19%) during Period II (complete nitrification and minimal NH2Cl residual). Nitrospira moscoviensis (nitrite-oxidizing bacteria) was detected in early periods (2%) but decreased to <0.02% in later periods, corresponding to nitrite accumulation. Methylobacterium (19%) and members of Nitrosomonadaceae (42%) dominated Period III (complete ammonia and partial nitrite oxidation and low NH2Cl residual). An increase in Afipia (haloacetic acid-degrading bacteria) relative abundance (<2% to 42%) occurred during Period IV (minimal nitrification and moderate to high NH2Cl residual). Microbial community and operational data provided no evidence of taxa-time relationship, but rapid community transitions indicated that the system had experienced ecological regime shifts to alternative stable states.

Meta'omic analytic techniques for studying the intestinal microbiome

Nucleotide sequencing has become increasingly common and affordable, and is now a vital tool for studies of the human microbiome. Comprehensive microbial community surveys such as MetaHit and the Human Microbiome Project have described the composition and molecular functional profile of the healthy (normal) intestinal microbiome. This knowledge will increase our ability to analyze host and microbial DNA (genome) and RNA (transcriptome) sequences. Bioinformatic and statistical tools then can be used to identify dysbioses that might cause disease, and potential treatments. Analyses that identify perturbations in specific molecules can leverage thousands of culture-based isolate genomes to contextualize culture-independent sequences, or may integrate sequence data with whole-community functional assays such as metaproteomic or metabolomic analyses. We review the state of available systems-level models for studies of the intestinal microbiome, along with analytic techniques and tools that can be used to determine its functional capabilities in healthy and unhealthy individuals.

Towards molecular biomarkers for biogas production from lignocellulose-rich substrates

Biogas production from lignocellulose-rich agricultural residues is gaining increasingly importance in sustainable energy production. Hydrolysis/acidogenesis (H/A) of lignocellulose as the initial rate-limiting step deserves particular optimization. A mixture of straw/hay was methanized applying two-phase digester systems with an initial H/A reactor and a one-stage system at different, meso- and thermophilic temperatures. H/A was intensified with increasing pH values and increasing temperature. H/A fermenters, however, were prone to switch to methanogenic systems at these conditions. Substrate turnover was accelerated in the bi-phasic process but did not reach the methanation efficiency of the single-stage digestion. There was no indication that two different cellulolytic inocula could establish in the given process. Bacterial communities were analyzed applying conventional amplicon clone sequencing targeting the hypervariable 16S rRNA gene region V6-V8 and by metagenome analyses applying direct DNA pyrosequencing without a PCR step. Corresponding results suggested that PCR did not introduce a bias but offered better phylogenetic resolution. Certain Clostridium IV and Prevotella members were most abundant in the H/A system operated at 38 °C, certain Clostridium III and Lachnospiraceae bacteria in the 45 °C, and certain Clostridium IV and Thermohydrogenium/Thermoanaerobacterium members in the 55 °C H/A system. Clostridium III representatives, Lachnospiraceae and Thermotogae dominated in the thermophilic single-stage system, in which also a higher portion of known syntrophic acetate oxidizers was found. Specific (RT-)qPCR systems were designed and applied for the most significant and abundant populations to assess their activity in the different digestion systems. The RT-qPCR results agreed with the DNA based community profiles obtained at the different temperatures. Up to 10(12) 16S rRNA copies mL(-1) were determined in H/A fermenters with prevalence of rRNA of a Ruminococcaceae subgroup. Besides, Thermohydrogenium/Thermoanaerobacterium rRNA prevailed at thermophilic and Prevotellaceae rRNA at mesophilic conditions. The developed (RT)-qPCR systems can be used as biomarkers to optimize biogas production from straw/hay and possibly other lignocellulosic substrates.

Diversity and transport of microorganisms in intertidal sands of the California coast

Forced by tides and waves, large volumes of seawater are flushed through the beach daily. Organic material and nutrients in seawater are remineralized and cycled as they pass through the beach. Microorganisms are responsible for most of the biogeochemical cycling in the beach; however, few studies have characterized their diversity in intertidal sands, and little work has characterized the extent to which microbes are transported between different compartments of the beach. The present study uses next-generation massively parallel sequencing to characterize the microbial community present at 49 beaches along the coast of California. In addition, we characterize the transport of microorganisms within intertidal sands using laboratory column experiments. We identified extensive diversity in the beach sands. Nearly 1,000 unique taxa were identified in sands from 10 or more unique beaches, suggesting the existence of a group of "cosmopolitan" sand microorganisms. A biogeographical analysis identified a taxon-distance relationship among the beaches. In addition, sands with similar grain size, organic carbon content, exposed to a similar wave climate, and having the same degree of anthropogenic influence tended to have similar microbial communities. Column experiments identified microbes readily mobilized by seawater infiltrating through unsaturated intertidal sands. The ease with which microbes were mobilized suggests that intertidal sands may represent a reservoir of bacteria that seed the beach aquifer where they may partake in biogeochemical cycling.