Search

We consider a Bayesian functional data analysis for observations measured as extremely long sequences. Splitting the sequence into several small windows with manageable lengths, the windows may not be independent especially when they are neighboring each other. We propose to utilize Bayesian smoothing splines to estimate individual functional patterns within each window and to establish transition models for parameters involved in each window to address the dependence structure between windows. The functional difference of groups of individuals at each window can be evaluated by the Bayes factor based on Markov Chain Monte Carlo samples in the analysis. In this paper, we examine the proposed method through simulation studies and apply it to identify differentially methylated genetic regions in TCGA lung adenocarcinoma data. [ABSTRACT FROM AUTHOR]

The performance of computational methods and software to identify differentially expressed features in single‐cell RNA‐sequencing (scRNA‐seq) has been shown to be influenced by several factors, including the choice of the normalization method used and the choice of the experimental platform (or library preparation protocol) to profile gene expression in individual cells. Currently, it is up to the practitioner to choose the most appropriate differential expression (DE) method out of over 100 DE tools available to date, each relying on their own assumptions to model scRNA‐seq expression features. To model the technological variability in cross‐platform scRNA‐seq data, here we propose to use Tweedie generalized linear models that can flexibly capture a large dynamic range of observed scRNA‐seq expression profiles across experimental platforms induced by platform‐ and gene‐specific statistical properties such as heavy tails, sparsity, and gene expression distributions. We also propose a zero‐inflated Tweedie model that allows zero probability mass to exceed a traditional Tweedie distribution to model zero‐inflated scRNA‐seq data with excessive zero counts. Using both synthetic and published plate‐ and droplet‐based scRNA‐seq datasets, we perform a systematic benchmark evaluation of more than 10 representative DE methods and demonstrate that our method (Tweedieverse) outperforms the state‐of‐the‐art DE approaches across experimental platforms in terms of statistical power and false discovery rate control. Our open‐source software (R/Bioconductor package) is available at https://github.com/himelmallick/Tweedieverse. [ABSTRACT FROM AUTHOR]

Maternal genetic variants associated with offspring birth weight and adult type 2 diabetes (T2D) risk loci show some overlap. Whether T2D genetic risk influences longitudinal fetal weight and the gestational timing when these relationships begin is unknown. We investigated the associations of T2D genetic risk scores (GRS) with longitudinal fetal weight and birth weight among 1,513 pregnant women from four ancestral groups. Women had up to five ultrasonography examinations. Ancestry-matched GRS were constructed separately using 380 European- (GRSeur), 104 African- (GRSafr), and 189 East Asian- (GRSeas) related T2D loci discovered in different population groups. Among European Americans, the highest quartile GRSeur was significantly associated with 53.8 g higher fetal weight (95% CI 19.2-88.5) over the pregnancy. The associations began at gestational week 24 and continued through week 40, with a 106.8 g (95% CI 6.5-207.1) increase in birth weight. The findings were similar in analysis further adjusted for maternal glucose challenge test results. No consistent association was found using ancestry-matched or cross-ancestry GRS in non-Europeans. In conclusion, T2D genetic susceptibility may influence fetal growth starting at midsecond trimester among Europeans. Absence of similar associations in non-Europeans urges the need for further genetic T2D studies in diverse ancestries. [ABSTRACT FROM AUTHOR]

It is challenging to associate features such as human health outcomes, diet, environmental conditions, or other metadata to microbial community measurements, due in part to their quantitative properties. Microbiome multi-omics are typically noisy, sparse (zero-inflated), high-dimensional, extremely non-normal, and often in the form of count or compositional measurements. Here we introduce an optimized combination of novel and established methodology to assess multivariable association of microbial community features with complex metadata in population-scale observational studies. Our approach, MaAsLin 2 (Microbiome Multivariable Associations with Linear Models), uses generalized linear and mixed models to accommodate a wide variety of modern epidemiological studies, including cross-sectional and longitudinal designs, as well as a variety of data types (e.g., counts and relative abundances) with or without covariates and repeated measurements. To construct this method, we conducted a large-scale evaluation of a broad range of scenarios under which straightforward identification of meta-omics associations can be challenging. These simulation studies reveal that MaAsLin 2's linear model preserves statistical power in the presence of repeated measures and multiple covariates, while accounting for the nuances of meta-omics features and controlling false discovery. We also applied MaAsLin 2 to a microbial multi-omics dataset from the Integrative Human Microbiome (HMP2) project which, in addition to reproducing established results, revealed a unique, integrated landscape of inflammatory bowel diseases (IBD) across multiple time points and omics profiles. Author summary: Recently, several statistical methods have been proposed to identify phenotypic or environmental associations with features (e.g., taxa, genes, pathways, chemicals, etc.) from molecular profiles of microbial communities. Particularly for human microbiome epidemiology, however, most of these are primarily focused on univariable associations that analyze only one or a few environmental covariates. This is a critical gap to address, given the growing commonality of population-scale microbiome research and the complexity of associated study designs, including dietary, pharmaceutical, clinical, and environmental covariates, often with samples from multiple time points or tissues. Surprisingly, there have been no systematic evaluations of statistical analysis methods appropriate for such studies, nor consensus on appropriate methods for scalable microbiome epidemiology. To this end, we developed and validated a statistical model (MaAsLin) that provides both the first unified method and the first large-scale, comprehensive benchmarking of multivariable associations in population-scale microbial community studies. We hope that the MaAsLin 2 implementation, validated through extensive simulations and an application to HMP2 IBD multi-omics, will be helpful for researchers in future analysis of both human-associated and environmental microbial communities. [ABSTRACT FROM AUTHOR]

Maternal dyslipidemia during pregnancy has been associated with suboptimal fetal growth and increased cardiometabolic diseasse risk in offspring. Altered placental function driven by placental gene expression is a hypothesized mechanism underlying these associations. We tested the relationship between maternal plasma lipid concentrations and placental gene expression. Among 64 pregnant women from the NICHD Fetal Growth Studies–Singleton cohort with maternal first trimester plasma lipids we extracted RNA-Seq on placental samples obtained at birth. Placental gene co-expression networks were validated by regulatory network analysis that integrated transcription factors and gene expression, and genome-wide transcriptome analysis. Network analysis detected 24 gene co-expression modules in placenta, of which one module was correlated with total cholesterol (r = 0.27, P-value = 0.03) and LDL-C (r = 0.31, P-value = 0.01). Genes in the module (n = 39 genes) were enriched in inflammatory response pathways. Out of the 39 genes in the module, three known lipid-related genes (MPO , PGLYRP1 and LTF) and MAGEC2 were validated by the regulatory network analysis, and one known lipid-related gene (ALX4) and two germ-cell development-related genes (MAGEC2 and LUZP4) were validated by genome-wide transcriptome analysis. Placental gene expression signatures associated with unfavorable maternal lipid concentrations may be potential pathways underlying later life offspring cardiometabolic traits. Clinical Trial Registration: ClinicalTrials.gov, identifier NCT00912132. [ABSTRACT FROM AUTHOR]

Motivation The discovery of biologically interpretable and clinically actionable communities in heterogeneous omics data is a necessary first step toward deriving mechanistic insights into complex biological phenomena. Here, we present a novel clustering approach, omeClust , for community detection in omics profiles by simultaneously incorporating similarities among measurements and the overall complex structure of the data. Results We show that omeClust outperforms published methods in inferring the true community structure as measured by both sensitivity and misclassification rate on simulated datasets. We further validated omeClust in diverse, multiple omics datasets, revealing new communities and functionally related groups in microbial strains, cell line gene expression patterns and fetal genomic variation. We also derived enrichment scores attributable to putatively meaningful biological factors in these datasets that can serve as hypothesis generators facilitating new sets of testable hypotheses. Availability and implementation omeClust is open-source software, and the implementation is available online at http://github.com/omicsEye/omeClust. Supplementary information Supplementary data are available at Bioinformatics online. [ABSTRACT FROM AUTHOR]

It is well documented that maternal type 2 diabetes (T2D) is related to fetal overgrowth. Further, maternal T2D genetic risk score (GRS) has been associated with increased offspring birthweight. However, the earliest gestational age when fetal growth starts to be associated with T2D GRS is unknown. We investigated the associations between maternal T2D GRS and longitudinal ultrasound measured fetal weight. We included 1,513 women of diverse race/ethnicity from the NICHD Fetal Growth Studies-singleton cohort. Women were enrolled at gestational weeks 8-13 and randomly assigned to four ultrasound schedules to capture weekly fetal growth. Genome-wide single nucleotide polymorphisms (SNPs) previously known to be associated with T2D in European-, African-, and Asian- ancestry populations (380, 104, and 19 SNPs, respectively) were used to calculate weighted GRS. The analysis involved European GRS for Whites and Hispanics, African GRS for Blacks and Asian GRS for Asians. Associations between quartiles of GRS and fetal weight were adjusted for pre-pregnancy BMI, population structure, fetal sex and gestational weeks. Among Whites, the highest quartile of European GRS was related to a 53.8 g increment in fetal weight (95% CI 19.2-88.5 g) over the pregnancy, and across 24 to 40 weeks (from 15.8 g increase at week 24 to 196.0 g increase at week 40). Similar findings were observed in analysis that excluded women with impaired glucose tolerance. African-, European- and Asian-GRSs were not significantly associated with fetal weight among other race/ethnic groups. Genetic susceptibility to T2D based on European-specific GRS among Whites was related to increased fetal weight starting at week 24 even among normal glycemia women, suggesting an early influence of genetic factors on fetal growth. Absence of similar associations in non-Whites is likely due to the small number of established T2D genetic loci warranting the need for genetic studies of T2D in diverse populations. Disclosure: M. Ouidir: None. X. Zeng: None. S. Chatterjee: None. C. Zhang: None. F. Tekola-ayele: None. [ABSTRACT FROM AUTHOR]