Dataset Information

Zero-Inflated Beta Regression for Differential Abundance Analysis with Metagenomics Data.

ABSTRACT: Metagenomics data have been growing rapidly due to the advances in NGS technologies. One goal of human microbial studies is to detect abundance differences across clinical conditions. Besides small sample size and high dimension, metagenomics data are usually represented as compositions (proportions) with a large number of zeros and skewed distribution. Efficient tools for handling such compositional data need to be developed. We propose a zero-inflated beta regression approach (ZIBSeq) for identifying differentially abundant features between multiple clinical conditions. The proposed method takes the sparse nature of metagenomics data into account and handle the compositional data efficiently. Compared with other available methods, the proposed approach demonstrates better performance with large AUC values for most simulation studies. When applied to a human metagenomics data, it also identifies biologically important taxa reported from previous studies. The software in R is available upon request from the first author.

SUBMITTER: Peng X

PROVIDER: S-EPMC6109378 | biostudies-literature | 2016 Feb

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Zero-Inflated Beta Regression for Differential Abundance Analysis with Metagenomics Data.

Peng Xiaoling X Li Gang G Liu Zhenqiu Z

Journal of computational biology : a journal of computational molecular cell biology 20151216 2

Metagenomics data have been growing rapidly due to the advances in NGS technologies. One goal of human microbial studies is to detect abundance differences across clinical conditions. Besides small sample size and high dimension, metagenomics data are usually represented as compositions (proportions) with a large number of zeros and skewed distribution. Efficient tools for handling such compositional data need to be developed. We propose a zero-inflated beta regression approach (ZIBSeq) for iden ...[more]

PMID: 26675626

Similar Datasets

Project description:BackgroundTsetse flies are the major vectors of human trypanosomiasis of the form Trypanosoma brucei rhodesiense and T.b.gambiense. They are widely spread across the sub-Saharan Africa and rendering a lot of challenges to both human and animal health. This stresses effective agricultural production and productivity in Africa. Delimiting the extent and magnitude of tsetse coverage has been a challenge over decades due to limited resources and unsatisfactory technology. In a bid to overcome these limitations, this study attempted to explore modelling skills that can be applied to spatially estimate tsetse abundance in the country using limited tsetse data and a set of remote-sensed environmental variables.MethodologyEntomological data for the period 2008-2018 as used in the model were obtained from various sources and systematically assembled using a structured protocol. Data harmonisation for the purposes of responsiveness and matching was carried out. The key tool for tsetse trapping was itemized as pyramidal trap in many instances and biconical trap in others. Based on the spatially explicit assembled data, we ran two regression models; standard Poisson and Zero-Inflated Poisson (ZIP), to explore the associations between tsetse abundance in Uganda and several environmental and climatic covariates. The covariate data were constituted largely by satellite sensor data in form of meteorological and vegetation surrogates in association with elevation and land cover data. We finally used the Zero-Inflated Poisson (ZIP) regression model to predict tsetse abundance due to its superiority over the standard Poisson after model fitting and testing using the Vuong Non-Nested statistic.ResultsA total of 1,187 tsetse sampling points were identified and considered as representative for the country. The model results indicated the significance and level of responsiveness of each covariate in influencing tsetse abundance across the study area. Woodland vegetation, elevation, temperature, rainfall, and dry season normalised difference vegetation index (NDVI) were important in determining tsetse abundance and spatial distribution at varied scales. The resultant prediction map shows scaled tsetse abundance with estimated fitted numbers ranging from 0 to 59 flies per trap per day (FTD). Tsetse abundance was found to be largest at low elevations, in areas of high vegetative activity, in game parks, forests and shrubs during the dry season. There was very limited responsiveness of selected predictors to tsetse abundance during the wet season, matching the known fact that tsetse disperse most significantly during wet season.ConclusionsA methodology was advanced to enable compilation of entomological data for 10 years, which supported the generation of tsetse abundance maps for Uganda through modelling. Our findings indicate the spatial distribution of the G. f. fuscipes as; low 0-5 FTD (48%), medium 5.1-35 FTD (18%) and high 35.1-60 FTD (34%) grounded on seasonality. This approach, amidst entomological data shortages due to limited resources and absence of expertise, can be adopted to enable mapping of the vector to provide better decision support towards designing and implementing targeted tsetse and tsetse-transmitted African trypanosomiasis control strategies.

Project description:BackgroundThe rapid growth of high-throughput sequencing-based microbiome profiling has yielded tremendous insights into human health and physiology. Data generated from high-throughput sequencing of 16S rRNA gene amplicons are often preprocessed into composition or relative abundance. However, reproducibility has been lacking due to the myriad of different experimental and computational approaches taken in these studies. Microbiome studies may report varying results on the same topic, therefore, meta-analyses examining different microbiome studies to provide consistent and robust results are important. So far, there is still a lack of implemented methods to properly examine differential relative abundances of microbial taxonomies and to perform meta-analysis examining the heterogeneity and overall effects across microbiome studies.ResultsWe developed an R package 'metamicrobiomeR' that applies Generalized Additive Models for Location, Scale and Shape (GAMLSS) with a zero-inflated beta (BEZI) family (GAMLSS-BEZI) for analysis of microbiome relative abundance datasets. Both simulation studies and application to real microbiome data demonstrate that GAMLSS-BEZI well performs in testing differential relative abundances of microbial taxonomies. Importantly, the estimates from GAMLSS-BEZI are log (odds ratio) of relative abundances between comparison groups and thus are analogous between microbiome studies. As such, we also apply random effects meta-analysis models to pool estimates and their standard errors across microbiome studies. We demonstrate the meta-analysis examples and highlight the utility of our package on four studies comparing gut microbiomes between male and female infants in the first six months of life.ConclusionsGAMLSS-BEZI allows proper examination of microbiome relative abundance data. Random effects meta-analysis models can be directly applied to pool comparable estimates and their standard errors to evaluate the overall effects and heterogeneity across microbiome studies. The examples and workflow using our 'metamicrobiomeR' package are reproducible and applicable for the analyses and meta-analyses of other microbiome studies.

Dataset Information

Zero-Inflated Beta Regression for Differential Abundance Analysis with Metagenomics Data.

Publications

Zero-Inflated Beta Regression for Differential Abundance Analysis with Metagenomics Data.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets