Citations

The genomic DNA input for stool samples utilized for microbiome composition has not been determined. In this study, we determined the reliable threshold level under which conclusions drawn from the data may be compromised. We also determined the type of microbial bias introduced by less-than-ideal genomic input.

Normalization is the first critical step in microbiome sequencing data analysis used to account for variable library sizes. Current RNA-Seq based normalization methods that have been adapted for microbiome data fail to consider the unique characteristics of microbiome data, which contain a vast number of zeros due to the physical absence or under-sampling of the microbes. Normalization methods that specifically address the zero inflation remain largely undeveloped. Here we propose GMPR - a simple but effective normalization method - for zero-inflated sequencing data such as microbiome data. Simulation studies and real datasets analyses demonstrate that the proposed method is more robust than competing methods, leading to more powerful detection of differentially abundant taxa and higher reproducibility of the relative abundances of taxa.

Normalization is the first critical step in microbiome sequencing data analysis used to account for variable library sizes. Current RNA-Seq based normalization methods that have been adapted for microbiome data fail to consider the unique characteristics of microbiome data, which contain a vast number of zeros due to the physical absence or under-sampling of the microbes. Normalization methods that specifically address the zero inflation remain largely undeveloped. Here we propose GMPR - a simple but effective normalization method - for zero-inflated sequencing data such as microbiome data. Simulation studies and real datasets analyses demonstrate that the proposed method is more robust than competing methods, leading to more powerful detection of differentially abundant taxa and higher reproducibility of the relative abundances of taxa.

Normalization is the first critical step in microbiome sequencing data analysis used to account for variable library sizes. Current RNA-Seq based normalization methods that have been adapted for microbiome data fail to consider the unique characteristics of microbiome data, which contain a vast number of zeros due to the physical absence or under-sampling of the microbes. Normalization methods that specifically address the zero inflation remain largely undeveloped. Here we propose GMPR - a simple but effective normalization method - for zero-inflated sequencing data such as microbiome data. Simulation studies and real datasets analyses demonstrate that the proposed method is more robust than competing methods, leading to more powerful detection of differentially abundant taxa and higher reproducibility of the relative abundances of taxa.

The clinical utility of microbiome biomarkers depends on the reliable and reproducible nature of comparative results. Underappreciation of the variation associated with common demographic, health, and behavioral factors may confound associations of interest and generate false positives. Here, we present the Midwestern Reference Panel (MWRP), a resource for comparative gut microbiome studies conducted in the Midwestern United States. We analyzed the relationships between demographic and health behavior–related factors and the microbiota in this cohort, and estimated their effect sizes. Most variables investigated were associated with the gut microbiota. Specifically, body mass index (BMI), race, sex, and alcohol use were significantly associated with microbial β-diversity (P < 0.05, unweighted UniFrac). BMI, race and alcohol use were also significantly associated with microbial α-diversity (P < 0.05, species richness). Tobacco use showed a trend toward association with the microbiota (P < 0.1, unweighted UniFrac). The effect sizes of the associations, as quantified by adjusted R 2 values based on unweighted UniFrac distances, were small (< 1% for all variables), indicating that these factors explain only a small percentage of overall microbiota variability. Nevertheless, the significant associations between these variables and the gut microbiota suggest that they could still be potential confounders in comparative studies and that controlling for these variables in study design, which is the main objective of the MWRP, is important for increasing reproducibility in comparative microbiome studies.

Abstract Motivation Recent studies have revealed a complex interplay between environment, the human microbiome and health and disease. Mediation analysis of the human microbiome in these complex relationships could potentially provide insights into the role of the microbiome in the etiology of disease and, more importantly, lead to novel clinical interventions by modulating the microbiome. However, due to the high dimensionality, sparsity, non-normality and phylogenetic structure of microbiome data, none of the existing methods are suitable for testing such clinically important mediation effect. Results We propose a distance-based approach for testing the mediation effect of the human microbiome. In the framework, the nonlinear relationship between the human microbiome and independent/dependent variables is captured implicitly through the use of sample-wise ecological distances, and the phylogenetic tree information is conveniently incorporated by using phylogeny-based distance metrics. Multiple distance metrics are utilized to maximize the power to detect various types of mediation effect. Simulation studies demonstrate that our method has correct Type I error control, and is robust and powerful under various mediation models. Application to a real gut microbiome dataset revealed that the association between the dietary fiber intake and body mass index was mediated by the gut microbiome. Availability and implementation An R package ‘MedTest’ is freely available at https://github.com/jchen1981/MedTest. Supplementary information Supplementary data are available at Bioinformatics online.

Abstract Background Illumina paired-end sequencing has been increasingly popular for 16S rRNA gene-based microbiota profiling. It provides higher phylogenetic resolution than single-end reads due to a longer read length. However, the reverse read (R2) often has significant low base quality, and a large proportion of R2s will be discarded after quality control, resulting in a mixture of paired-end and single-end reads. A typical 16S analysis pipeline usually processes either paired-end or single-end reads but not a mixture. Thus, the quantification accuracy and statistical power will be reduced due to the loss of a large amount of reads. As a result, rare taxa may not be detectable with the paired-end approach, or low taxonomic resolution will result in a single-end approach. Results To have both the higher phylogenetic resolution provided by paired-end reads and the higher sequence coverage by single-end reads, we propose a novel OTU-picking pipeline, hybrid-denovo, that can process a hybrid of single-end and paired-end reads. Using high-quality paired-end reads as a gold standard, we show that hybrid-denovo achieved the highest correlation with the gold standard and performed better than the approaches based on paired-end or single-end reads in terms of quantifying the microbial diversity and taxonomic abundances. By applying our method to a rheumatoid arthritis (RA) data set, we demonstrated that hybrid-denovo captured more microbial diversity and identified more RA-associated taxa than a paired-end or single-end approach. Conclusions Hybrid-denovo utilizes both paired-end and single-end 16S sequencing reads and is recommended for 16S rRNA gene targeted paired-end sequencing data.

AbstractMotivationOne objective of human microbiome studies is to identify differentially abundant microbes across biological conditions. Previous statistical methods focus on detecting the shift in the abundance and/or prevalence of the microbes and treat the dispersion (spread of the data) as a nuisance. These methods also assume that the dispersion is the same across conditions, an assumption which may not hold in presence of sample heterogeneity. Moreover, the widespread outliers in the microbiome sequencing data make existing parametric models not overly robust. Therefore, a robust and powerful method that allows covariate-dependent dispersion and addresses outliers is still needed for differential abundance analysis.ResultsWe introduce a novel test for differential distribution analysis of microbiome sequencing data by jointly testing the abundance, prevalence and dispersion. The test is built on a zero-inflated negative binomial regression model and winsorized count data to account for zero-inflation and outliers. Using simulated data and real microbiome sequencing datasets, we show that our test is robust across various biological conditions and overall more powerful than previous methods.Availability and implementationR package is available at https://github.com/jchen1981/MicrobiomeDDA.Supplementary informationSupplementary data are available at Bioinformatics online.

PurposeTo assess the success rate, image quality, and the ability to stage liver fibrosis of a standard 2D gradient‐recalled echo (GRE) and four different spin‐echo (SE) magnetic resonance elastography (MRE) sequences in patients with different liver iron concentrations.Materials and MethodsA total of 332 patients who underwent 3T MRE examinations that included liver fat and iron quantification were enrolled, including 136 patients with all five MRE techniques. Thirty‐four patients had biopsy results for fibrosis staging. The liver stiffness, region of interest area, image quality, and success rate of the five sequences were compared in 115/136 patients. The area under the receiver operating characteristic curves (AUCs) and the accuracies for diagnosing early‐stage fibrosis and advanced fibrosis were compared. The effect of BMI (body mass index), the R2* relaxation time, and fat fraction on the image quality and liver stiffness measurements were analyzed.ResultsThe success rates were significantly higher in the four SE sequences (99.1–100%) compared with GRE MRE (85.3%) (all P < 0.001). There were significant differences of the mean ROI area between every pair of sequences (all P < 0.0001). There were no significant differences in the AUC of the five MRE sequences for discriminating advanced fibrosis (10 P‐values ranging from 0.2410–0.9171). R2* had a significant effect on the success rate and image quality for the noniron 2D echo‐planar imaging (EPI), 3D EPI and 2D GRE (all P < 0.001) sequences. BMI had a significant effect on the iron 2D EPI (P = 0.0230) and iron 2D SE (P = 0.0040) sequences.ConclusionAll five techniques showed good diagnostic performance in staging liver fibrosis. The SE MRE sequences had higher success rates and better image quality than GRE MRE in 3T clinical hepatic imaging.Level of Evidence: 3Technical Efficacy: Stage 5J. Magn. Reson. Imaging 2018;47:976–987.

AbstractMotivationNext generation sequencing technologies have enabled the study of the human microbiome through direct sequencing of microbial DNA, resulting in an enormous amount of microbiome sequencing data. One unique characteristic of microbiome data is the phylogenetic tree that relates all the bacterial species. Closely related bacterial species have a tendency to exhibit a similar relationship with the environment or disease. Thus, incorporating the phylogenetic tree information can potentially improve the detection power for microbiome-wide association studies, where hundreds or thousands of tests are conducted simultaneously to identify bacterial species associated with a phenotype of interest. Despite much progress in multiple testing procedures such as false discovery rate (FDR) control, methods that take into account the phylogenetic tree are largely limited.ResultsWe propose a new FDR control procedure that incorporates the prior structure information and apply it to microbiome data. The proposed procedure is based on a hierarchical model, where a structure-based prior distribution is designed to utilize the phylogenetic tree. By borrowing information from neighboring bacterial species, we are able to improve the statistical power of detecting associated bacterial species while controlling the FDR at desired levels. When the phylogenetic tree is mis-specified or non-informative, our procedure achieves a similar power as traditional procedures that do not take into account the tree structure. We demonstrate the performance of our method through extensive simulations and real microbiome datasets. We identified far more alcohol-drinking associated bacterial species than traditional methods.Availability and implementationR package StructFDR is available from CRAN.Supplementary informationSupplementary data are available at Bioinformatics online.