Dataset Information

Non-random sampling leads to biased estimates of transcriptome association.

ABSTRACT: Integration of independent data resources across -omics platforms offers transformative opportunity for novel clinical and biological discoveries. However, application of emerging analytic methods in the context of selection bias represents a noteworthy and pervasive challenge. We hypothesize that combining differentially selected samples for integrated transcriptome analysis will lead to bias in the estimated association between predicted expression and the trait. Our results are based on in silico investigations and a case example focused on body mass index across four well-described cohorts apparently derived from markedly different populations. Our findings suggest that integrative analysis can lead to substantial relative bias in the estimate of association between predicted expression and the trait. The average estimate of association ranged from 51.3% less than to 96.7% greater than the true value for the biased sampling scenarios considered, while the average error was - 2.7% for the unbiased scenario. The corresponding 95% confidence interval coverage rate ranged from 46.4% to 69.5% under biased sampling, and was equal to 75% for the unbiased scenario. Inverse probability weighting with observed and estimated weights is applied as one corrective measure and appears to reduce the bias and improve coverage. These results highlight a critical need to address selection bias in integrative analysis and to use caution in interpreting findings in the presence of different sampling mechanisms between groups.

SUBMITTER: Foulkes AS

PROVIDER: S-EPMC7148323 | biostudies-literature | 2020 Apr

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Non-random sampling leads to biased estimates of transcriptome association.

Foulkes A S AS Balasubramanian R R Qian J J Reilly M P MP

Scientific reports 20200410 1

Integration of independent data resources across -omics platforms offers transformative opportunity for novel clinical and biological discoveries. However, application of emerging analytic methods in the context of selection bias represents a noteworthy and pervasive challenge. We hypothesize that combining differentially selected samples for integrated transcriptome analysis will lead to bias in the estimated association between predicted expression and the trait. Our results are based on in si ...[more]

PMID: 32277087

Similar Datasets

Project description:The site-to-site variability in species composition, known as β-diversity, is crucial to understanding spatiotemporal patterns of species diversity and the mechanisms controlling community composition and structure. However, quantifying β-diversity in microbial ecology using sequencing-based technologies is a great challenge because of a high number of sequencing errors, bias, and poor reproducibility and quantification. Herein, based on general sampling theory, a mathematical framework is first developed for simulating the effects of random sampling processes on quantifying β-diversity when the community size is known or unknown. Also, using an analogous ball example under Poisson sampling with limited sampling efforts, the developed mathematical framework can exactly predict the low reproducibility among technically replicate samples from the same community of a certain species abundance distribution, which provides explicit evidences of random sampling processes as the main factor causing high percentages of technical variations. In addition, the predicted values under Poisson random sampling were highly consistent with the observed low percentages of operational taxonomic unit (OTU) overlap (<30% and <20% for two and three tags, respectively, based on both Jaccard and Bray-Curtis dissimilarity indexes), further supporting the hypothesis that the poor reproducibility among technical replicates is due to the artifacts associated with random sampling processes. Finally, a mathematical framework was developed for predicting sampling efforts to achieve a desired overlap among replicate samples. Our modeling simulations predict that several orders of magnitude more sequencing efforts are needed to achieve desired high technical reproducibility. These results suggest that great caution needs to be taken in quantifying and interpreting β-diversity for microbial community analysis using next-generation sequencing technologies. IMPORTANCE Due to the vast diversity and uncultivated status of the majority of microorganisms, microbial detection, characterization, and quantitation are of great challenge. Although large-scale metagenome sequencing technology such as PCR-based amplicon sequencing has revolutionized the studies of microbial communities, it suffers from several inherent drawbacks, such as a high number of sequencing errors, biases, poor quantitation, and very high percentages of technical variations, which could greatly overestimate microbial biodiversity. Based on general sampling theory, this study provided the first explicit evidence to demonstrate the importance of random sampling processes in estimating microbial β-diversity, which has not been adequately recognized and addressed in microbial ecology. Since most ecological studies are involved in random sampling, the conclusions learned from this study should also be applicable to other ecological studies in general. In summary, the results presented in this study should have important implications for examining microbial biodiversity to address both basic theoretical and applied management questions.

Project description:IntroductionNon-linear Mendelian randomization is an extension of conventional Mendelian randomization that performs separate instrumental variable analyses in strata of the study population with different average levels of the exposure. The approach estimates a localized average causal effect function, representing the average causal effect of the exposure on the outcome at different levels of the exposure. The commonly used residual method for dividing the population into strata works under the assumption that the effect of the genetic instrument on the exposure is linear and constant in the study population. However, this assumption may not hold in practice.MethodsWe use the recently developed doubly ranked method to re-analyse various datasets previously analysed using the residual method. In particular, we consider a genetic score for 25-hydroxyvitamin D (25[OH]D) used in a recent non-linear Mendelian randomization analysis to assess the potential effect of vitamin D supplementation on all-cause mortality.ResultsThe effect of the genetic score on 25(OH)D concentrations varies strongly, with a five-fold difference in the estimated genetic association with the exposure in the lowest and highest decile groups. Evidence for a protective causal effect of vitamin D supplementation on all-cause mortality in low vitamin D individuals is evident for the residual method but not for the doubly ranked method. We show that the constant genetic effect assumption is more reasonable for some exposures and less reasonable for others. If the doubly ranked method indicates that this assumption is violated, then estimates from both the residual and doubly ranked methods can be biased, although bias was smaller on average in the doubly ranked method.ConclusionAnalysts wanting to perform non-linear Mendelian randomization should compare results from both the residual and doubly ranked methods, as well as consider transforming the exposure for the residual method to reduce heterogeneity in the genetic effect on the exposure.

Dataset Information

Non-random sampling leads to biased estimates of transcriptome association.

Publications

Non-random sampling leads to biased estimates of transcriptome association.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets