FlowClus: efficiently filtering and denoising pyrosequenced amplicons

Amplicon sequence variant-based meiofaunal community composition revealed by DADA2 tool is compatible with species composition

The present study is aimed at implementing the morphological identification-free amplicon sequence variant (ASV) concept for describing meiofaunal species composition, while strongly indicating reasonable compatibility with the underlying species. A primer pair was constructed and demonstrated to PCR amplify a 470-490 bp 18S barcode from a variety of meiofaunal taxa, high throughput sequenced using the Illumina 300 × 2 bps platform. Sixteen 18S multi-species HTS assemblies were created from meiofaunal samples and merged to one assembly of ~2,150,000 reads. Five quality scores (q = 35, 30, 25, 20, 15) were implemented to filter five 18S barcode assemblies, which served as inputs for the DADA2 software, ending with five reference ASV libraries. Each of these libraries was clustered, applying 3% dissimilarity threshold, revealed an average number of 1.38 ± 0.078 ASVs / cluster. Hence, demonstrating high level of ASV uniqueness. The libraries which were based on q ≤ 25 reached a near-asymptote number of ASVs which together with the low average number of ASVs / cluster, strongly indicated fair representation of the actual number of the underlying species. Hence, the q = 25 library was selected to be used as metabarcoding reference library. It contained 461 ASVs and 342-3% clusters with average number of 1.34 ± 1.036 ASV / cluster and their BLASTN annotation elucidated a variety of expected meiofaunal taxa. The sixteen assemblies of sample-specific paired reads were mapped to this reference library and sample ASV profiles, namely the list of ASVs and their proportional copy numbers were created and clustered.

Ecological impact of a rare sugar on grapevine phyllosphere microbial communities

Plants host a complex microbiota inside or outside their tissues, and phyllosphere microorganisms can be influenced by environmental, nutritional and agronomic factors. Rare sugars are defined as monosaccharides with limited availability in nature and they are metabolised by only few certain microbial taxa. Among rare sugars, tagatose (TAG) is a low-calories sweetener that stimulates and inhibits beneficial and pathogenic bacteria in the human gut microbiota, respectively. Based on this differential effect on human-associated microorganisms, we investigated the effect of TAG treatments on the grapevine phyllosphere microorganisms to evaluate whether it can engineer the microbiota and modify the ratio between beneficial and pathogenic plant-associated microorganisms. TAG treatments changed the structure of the leaf microbiota and they successfully reduced leaf infections of downy mildew (caused by Plasmopara viticola) and powdery mildew (caused by Erysiphe necator) under field conditions. TAG increased the relative abundance of indigenous beneficial microorganisms, such as some potential biocontrol agents, which could partially contribute to disease control. The taxonomic composition of fungal and bacterial leaf populations differed according to grapevine locations, therefore TAG effects on the microbial structure were influenced by the composition of the originally residing microbiota. TAG is a promising biopesticide that could shift the balance of pathogenic and beneficial plant-associated microorganisms, suggesting selective nutritional/anti-nutritional properties for some specific taxa. More specifically, TAG displayed possible plant prebiotic effects on the phyllosphere microbiota and this mechanism of action could represent a novel strategy that can be further developed for sustainable plant protection.

Linking Soil Fungal Generality to Tree Richness in Young Subtropical Chinese Forests

Soil fungi are a highly diverse group of microorganisms that provide many ecosystem services. The mechanisms of soil fungal community assembly must therefore be understood to reliably predict how global changes such as climate warming and biodiversity loss will affect ecosystem functioning. To this end, we assessed fungal communities in experimental subtropical forests by pyrosequencing of the internal transcribed spacer 2 (ITS2) region, and constructed tree-fungal bipartite networks based on the co-occurrence of fungal operational taxonomic units (OTUs) and tree species. The characteristics of the networks and the observed degree of fungal specialization were then analyzed in relation to the level of tree species diversity. Unexpectedly, plots containing two tree species had higher network connectance and fungal generality values than those with higher tree diversity. Most of the frequent fungal OTUs were saprotrophs. The degree of fungal specialization was highest in tree monocultures. Ectomycorrhizal fungi had higher specialization coefficients than saprotrophic, arbuscular mycorrhizal, and plant pathogenic fungi. High tree species diversity plots with 4 to 16 different tree species sustained the greatest number of fungal species, which is assumed to be beneficial for ecosystem services because it leads to more effective resource exploitation and greater resilience due to functional redundancy.

Long-read amplicon denoising

Long-read next-generation amplicon sequencing shows promise for studying complete genes or genomes from complex and diverse populations. Current long-read sequencing technologies have challenging error profiles, hindering data processing and incorporation into downstream analyses. Here we consider the problem of how to reconstruct, free of sequencing error, the true sequence variants and their associated frequencies from PacBio reads. Called 'amplicon denoising', this problem has been extensively studied for short-read sequencing technologies, but current solutions do not always successfully generalize to long reads with high indel error rates. We introduce two methods: one that runs nearly instantly and is very accurate for medium length reads and high template coverage, and another, slower method that is more robust when reads are very long or coverage is lower. On two Mock Virus Community datasets with ground truth, each sequenced on a different PacBio instrument, and on a number of simulated datasets, we compare our two approaches to each other and to existing algorithms. We outperform all tested methods in accuracy, with competitive run times even for our slower method, successfully discriminating templates that differ by a just single nucleotide. Julia implementations of Fast Amplicon Denoising (FAD) and Robust Amplicon Denoising (RAD), and a webserver interface, are freely available.

Amelioration of obesity-related biomarkers by Lactobacillus sakei CJLS03 in a high-fat diet-induced obese murine model

Recent progresses in clinical diagnostic analyses have demonstrated the decisive influence of host gut microbiota on the status of metabolic disorders. Short chain fatty acids (SCFAs) produced by gut microbiota, in particular, are considered as a key biomarker, both of communication between gut microbiota and the host, and of impact on host metabolic homeostasis. Microbiota modulation and concomitant anti-obesity effects of probiotics have been reported by different researchers. However, the underlying modulatory functions of probiotics on gut microbiota towards host metabolic homeostasis are still not fully understood. In this study, the impact of Lactobacillus sakei CJLS03 (isolated from Korean kimchi) on obesity-related biomarkers was investigated using a diet-induced obese mouse model. Body weight increase, SCFAs, the gut microbiota and various obesity-associated biomarkers were significantly and beneficially influenced by L. sakei CJLS03 administration compared to the control groups. Analytical data on faecal samples support the role of the colonic microbial population in SCFA production. The composition of the latter may be influenced by modulation of the distal gastro-intestinal microbiota by putative probiotics such as L. sakei CJLS03.

NGmerge: merging paired-end reads via novel empirically-derived models of sequencing errors

BACKGROUND

Advances in Illumina DNA sequencing technology have produced longer paired-end reads that increasingly have sequence overlaps. These reads can be merged into a single read that spans the full length of the original DNA fragment, allowing for error correction and accurate determination of read coverage. Extant merging programs utilize simplistic or unverified models for the selection of bases and quality scores for the overlapping region of merged reads.

RESULTS

We first examined the baseline quality score - error rate relationship using sequence reads derived from PhiX. In contrast to numerous published reports, we found that the quality scores produced by Illumina were not substantially inflated above the theoretical values, once the reference genome was corrected for unreported sequence variants. The PhiX reads were then used to create empirical models of sequencing errors in overlapping regions of paired-end reads, and these models were incorporated into a novel merging program, NGmerge. We demonstrate that NGmerge corrects errors and ambiguous bases better than other merging programs, and that it assigns quality scores for merged bases that accurately reflect the error rates. Our results also show that, contrary to published analyses, the sequencing errors of paired-end reads are not independent.

CONCLUSIONS

We provide a free and open-source program, NGmerge, that performs better than existing read merging programs. NGmerge is available on GitHub ( https://github.com/harvardinformatics/NGmerge ) under the MIT License; it is written in C and supported on Linux.

Experimental Evidence of Functional Group-Dependent Effects of Tree Diversity on Soil Fungi in Subtropical Forests

Deconvoluting the relative contributions made by specific biotic and abiotic drivers to soil fungal community compositions facilitates predictions about the functional responses of ecosystems to environmental changes, such as losses of plant diversity, but it is hindered by the complex interactions involved. Experimental assembly of tree species allows separation of the respective effects of plant community composition (biotic components) and soil properties (abiotic components), enabling much greater statistical power than can be achieved in observational studies. We therefore analyzed these contributions by assessing, via pyrotag sequencing of the internal transcribed spacer (ITS2) rDNA region, fungal communities in young subtropical forest plots included in a large experiment on the effects of tree species richness. Spatial variables and soil properties were the main drivers of soil fungal alpha and beta-diversity, implying strong early-stage environmental filtering and dispersal limitation. Tree related variables, such as tree community composition, significantly affected arbuscular mycorrhizal and pathogen fungal community structure, while differences in tree host species and host abundance affected ectomycorrhizal fungal community composition. At this early stage of the experiment, only a limited amount of carbon inputs (rhizodeposits and leaf litter) was being provided to the ecosystem due to the size of the tree saplings, and persisting legacy effects were observed. We thus expect to find increasing tree related effects on fungal community composition as forest development proceeds.

A New Perspective on Polyploid Fragaria (Strawberry) Genome Composition Based on Large-Scale, Multi-Locus Phylogenetic Analysis

The subgenomic compositions of the octoploid (2n = 8× = 56) strawberry (Fragaria) species, including the economically important cultivated species Fragaria x ananassa, have been a topic of long-standing interest. Phylogenomic approaches utilizing next-generation sequencing technologies offer a new window into species relationships and the subgenomic compositions of polyploids. We have conducted a large-scale phylogenetic analysis of Fragaria (strawberry) species using the Fluidigm Access Array system and 454 sequencing platform. About 24 single-copy or low-copy nuclear genes distributed across the genome were amplified and sequenced from 96 genomic DNA samples representing 16 Fragaria species from diploid (2×) to decaploid (10×), including the most extensive sampling of octoploid taxa yet reported. Individual gene trees were constructed by different tree-building methods. Mosaic genomic structures of diploid Fragaria species consisting of sequences at different phylogenetic positions were observed. Our findings support the presence in octoploid species of genetic signatures from at least five diploid ancestors (F. vesca, F. iinumae, F. bucharica, F. viridis, and at least one additional allele contributor of unknown identity), and questions the extent to which distinct subgenomes are preserved over evolutionary time in the allopolyploid Fragaria species. In addition, our data support divergence between the two wild octoploid species, F. virginiana and F. chiloensis.

Modulation of Active Gut Microbiota by Lactobacillus rhamnosus GG in a Diet Induced Obesity Murine Model

Gut microbiota play a key role in the development of metabolic disorders. Defining and correlating structural shifts in gut microbial assemblages with conditions related to metabolic syndrome have, however, been proven difficult. Results from 16S genomic DNA and 16S ribosomal RNA analyses of fecal samples may differ widely, leading to controversial information on the whole microbial community and metabolically active microbiota. Using a C57BL/6J murine model, we compared data from 16S genomic DNA and ribosomal RNA of the fecal microbiota. The study included three groups of experimental animals comprising two groups with high fat diet induced obesity (DIO) while a third group (control) received a low fat diet. One of the DIO groups was treated with the probiotic Lactobacillus rhamnosus GG (LGG). Compared to the data obtained by DNA analysis, a significantly higher abundance of OTUs was accounted for by RNA analysis. Moreover, rRNA based analysis showed a modulation of the active gut microbial population in the DIO group receiving LGG, thus reflecting a change in the induced obesity status of the host. As one of the most widely studied probiotics the functionality of LGG has been linked to the alleviation of metabolic syndrome, and, in some cases, to an impact on the microbiome. Yet, it appears that no study has reported thus far on modulation of the active microbiota by LGG treatment. It is postulated that the resulting impact on calorie consumption affects weight gain concomitantly with modulation of the functional structure of the gut microbial population. Using the 16S rRNA based approach therefore decisively increased the precision of gut microbiota metagenome analysis.

Analytical parameters and validation of homopolymer detection in a pyrosequencing-based next generation sequencing system

Background

Current technologies in next-generation sequencing are offering high throughput reads at low costs, but still suffer from various sequencing errors. Although pyro- and ion semiconductor sequencing both have the advantage of delivering long and high quality reads, problems might occur when sequencing homopolymer-containing regions, since the repeating identical bases are going to incorporate during the same synthesis cycle, which leads to uncertainty in base calling. The aim of this study was to evaluate the analytical performance of a pyrosequencing-based next-generation sequencing system in detecting homopolymer sequences using homopolymer-preintegrated plasmid constructs and human DNA samples originating from patients with cystic fibrosis.

Results

In the plasmid system average correct genotyping was 95.8% in 4-mers, 87.4% in 5-mers and 72.1% in 6-mers. Despite the experienced low genotyping accuracy in 5- and 6-mers, it was possible to generate amplicons with more than a 90% adequate detection rate in every homopolymer tract. When homopolymers in the CFTR gene were sequenced average accuracy was 89.3%, but varied in a wide range (52.2 – 99.1%). In all but one case, an optimal amplicon-sequencing primer combination could be identified. In that single case (7A tract in exon 14 (c.2046_2052)), none of the tested primer sets produced the required analytical performance.

Conclusions

Our results show that pyrosequencing is the most reliable in case of 4-mers and as homopolymer length gradually increases, accuracy deteriorates. With careful primer selection, the NGS system was able to correctly genotype all but one of the homopolymers in the CFTR gene. In conclusion, we configured a plasmid test system that can be used to assess genotyping accuracy of NGS devices and developed an accurate NGS assay for the molecular diagnosis of CF using self-designed primers for amplification and sequencing.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4544-x) contains supplementary material, which is available to authorized users.

Revealing higher than expected meiofaunal diversity in Antarctic sediments: a metabarcoding approach

An increasing number of studies are showing that Antarctic mega- and macrofauna are highly diverse, however, little is known about meiofaunal biodiversity in sediment communities, which are a vital part of a healthy and functional ecosystem. This is the first study to analyse community DNA (targeting meiofauna) using metabarcoding to investigate biodiversity levels in sediment communities of the Antarctic Peninsula. The results show that almost all of the meiofaunal biodiversity in the benthic habitat has yet to be characterised, levels of biodiversity were higher than expected and similar to temperate regions, albeit with the existence of potentially new and locally adapted species never described before at the molecular level. The Rothera meiofaunal sample sites showed four dominant eukaryotic groups, the nematodes, arthropods, platyhelminthes, and the annelids; some of which could comprise species complexes. Comparisons with deep-sea data from the same region suggest little exchange of Operational Taxonomic Units (OTUs) between depths with the nematodes prevalent at all depths, but sharing the shallow water benthos with the copepods. This study provides a preliminary analysis of benthic Antarctic Peninsula meiofauna using high throughput sequencing which substantiates how little is known on the biodiversity of one of the most diverse, yet underexplored communities of the Antarctic: the benthos.

Lactobacillus plantarum HAC01 regulates gut microbiota and adipose tissue accumulation in a diet-induced obesity murine model

The functional features of Lactobacillus plantarum HAC01 (HAC01), isolated from fermented Korean kimchi, were studied with regard to the fat mass, immunometabolic biomarkers and dysbiosis in a diet-induced obesity (DIO) murine model. L. rhamnosus GG (LGG) served as reference strain and a PBS-treated group as control. The administration of L. plantarum HAC01 resulted in reduction of the mesenteric adipose depot, the conjunctive tissue closely associated with the gastrointestinal tract, where lipid oxidative gene expression was upregulated compared to the control group. Metagenome analysis of intestinal microbiota showed that both strains HAC01 and LGG influenced specific bacterial families such as the Lachnospiraceae and Ruminococcaceae rather than the phyla Firmicutes and Bacteroidetes as a whole. The relative abundance of the Lachnospiraceae (phylum Firmicutes) was significantly higher in both LAB-treated groups than in the control. Comparing the impact of the two Lactobacillus strains on microbial composition in the gut also suggests strain-specific effects. The study emphasises the need for deeper studies into functional specificity of a probiotic organism at the strain level. Alleviation of obesity-associated dysbiosis by modulation of the gut microbiota appears to be associated with "indicator" bacterial taxa such as the family Lachnospiraceae. This may provide further insight into mechanisms basic to the mode of probiotic action against obesity and associated dysbiosis.

Bacterial Communities Associated with Atherosclerotic Plaques from Russian Individuals with Atherosclerosis

Atherosclerosis is considered a chronic disease of the arterial wall and is the major cause of severe disease and death among individuals all over the world. Some recent studies have established the presence of bacteria in atherosclerotic plaque samples and suggested their possible contribution to the development of cardiovascular disease. The main objective of this preliminary pilot study was to better understand the bacterial diversity and abundance in human atherosclerotic plaques derived from common carotid arteries of individuals with atherosclerosis (Russian nationwide group) and contribute towards the further identification of a main group of atherosclerotic plaque bacteria by 454 pyrosequencing their 16S ribosomal RNA (16S rRNA) genes. The applied approach enabled the detection of bacterial DNA in all atherosclerotic plaques. We found that distinct members of the order Burkholderiales were present at high levels in all atherosclerotic plaques obtained from patients with atherosclerosis with the genus Curvibacter being predominant in all plaque samples. Moreover, unclassified Burkholderiales as well as members of the genera Propionibacterium and Ralstonia were typically the most significant taxa for all atherosclerotic plaques. Other genera such as Burkholderia, Corynebacterium and Sediminibacterium as well as unclassified Comamonadaceae, Oxalobacteraceae, Rhodospirillaceae, Bradyrhizobiaceae and Burkholderiaceae were always found but at low relative abundances of the total 16S rRNA gene population derived from all samples. Also, we found that some bacteria found in plaque samples correlated with some clinical parameters, including total cholesterol, alanine aminotransferase and fibrinogen levels. Finally, our study indicates that some bacterial agents at least partially may be involved in affecting the development of cardiovascular disease through different mechanisms.

Characterizing root-associated fungal communities and soils of Douglas-fir (Pseudotsuga menziesii) stands that naturally produce Oregon white truffles (Tuber oregonense and Tuber gibbosum)

Many truffle species in the genus Tuber are endemic to North America. Some of these have commercial value such as Tuber oregonense and Tuber gibbosum, commonly known as Oregon white truffles. Most of what is known about the ecology of these truffles comes from observational data. These truffle species form ectomycorrhizas with Douglas-fir (Pseudotsuga menziesii) and sometimes fruit abundantly in early successional forest regrowth. The goal of this study was to characterize fungal communities and soils associated with truffle-producing Douglas-fir sites. We extracted DNA from roots of five trees at four different truffle-producing Douglas-fir sites (n = 20). We amplified the internal transcribed spacer (ITS) region of the nuclear ribosomal DNA (nrDNA) and sequenced amplicons with 454 pyrosequencing. After quality filtering, we assembled 15,713 sequences into 150 fungal operational taxonomic units (OTUs). Pezizomycetes (Tuber and Pyronemataceae) were the most abundant taxa detected followed by Helotiales. Agaricomycetes represented most by Thelephoraceae, Russulaceae, and Inocybaceae were also abundant. A total of five Tuber species were detected. T. oregonense was the most abundant OTU, followed by T. gibbosum and Wilcoxina mikolae. Fungal root endophytes were also detected and well represented by Chalara and Phialocephala spp. Fungal community structure and soil chemistry differed between sites. This study represents the first characterization of the fungal communities in Douglas-fir stands producing Oregon white truffles. We found that Tuber species can be dominant ectomycorrhizal symbionts of Douglas-fir. Truffle fungi are also important in forest health, food webs, and as a non-timber forest resource that can contribute to rural economies.

Protective effects of Lactobacillus rhamnosus GG against dyslipidemia in high-fat diet-induced obese mice

Recent reports suggest that gut microbiota can be a major determinant of dyslipidemia and non-alcoholic fatty liver disease (NAFLD) and its modulation by treating probiotics is a valid strategy to exert a protective effect. In this study, high-fat diet (HFD)-fed mice were orally administrated with Lactobacillus rhamnosus GG (LGG) for 13 weeks. Significant reductions in the weights of the liver, mesenteric and subcutaneous adipose tissues were observed in LGG-treated HFD-fed mice compared to LGG-non-treated controls. The serum levels of triglyceride and cholesterol were also significantly reduced in LGG-treated mice. Gut microbial composition analysis showed that shifts in the diversity of dominant gut bacteria were caused by HFD and restored by LGG treatment. A remarkable decrease of hepatic fat content was also observed in LGG-treated mice, accompanied by downregulated expressions of lipogenic and pro-inflammatory genes in the liver. LGG-treated mice had lower expression levels of genes involved in cholesterol synthesis, but conversely, higher expression levels of cholesterol efflux-related genes compared to LGG-non-treated controls. The cholesterol-lowering effect of LGG was also found to be mediated by suppression of FXR and FGF15 signaling, resulting in the upregulation of hepatic CYP7A1. Our findings confirm a therapeutic potential of probiotics for ameliorating dyslipidemia and NAFLD.

Multicenter quality assessment of 16S ribosomal DNA-sequencing for microbiome analyses reveals high inter-center variability

The composition of human as well as animal microbiota has increasingly gained in interest since metabolites and structural components of endogenous microorganisms fundamentally influence all aspects of host physiology. Since many of the bacteria are still unculturable, molecular techniques such as high-throughput sequencing have dramatically increased our knowledge of microbial communities. The majority of microbiome studies published thus far are based on bacterial 16S ribosomal RNA (rRNA) gene sequencing, so that they can, at least in principle, be compared to determine the role of the microbiome composition for host metabolism and physiology, developmental processes, as well as different diseases. However, differences in DNA preparation and purification, 16S rDNA PCR amplification, sequencing procedures and platforms, as well as bioinformatic analysis and quality control measures may strongly affect the microbiome composition results obtained in different laboratories. To systematically evaluate the comparability of results and identify the most influential methodological factors affecting these differences, identical human stool sample replicates spiked with quantified marker bacteria, and their subsequent DNA sequences were analyzed by nine different centers in an external quality assessment (EQA). While high intra-center reproducibility was observed in repetitive tests, significant inter-center differences of reported microbiota composition were obtained. All steps of the complex analysis workflow significantly influenced microbiome profiles, but the magnitude of variation caused by PCR primers for 16S rDNA amplification was clearly the largest. In order to advance microbiome research to a more standardized and routine medical diagnostic procedure, it is essential to establish uniform standard operating procedures throughout laboratories and to initiate regular proficiency testing.