Explore the vast range of titles available.

We present Bisque, a tool for estimating cell type proportions in bulk expression. Bisque implements a regression-based approach that utilizes single-cell RNA-seq (scRNA-seq) or single-nucleus RNA-seq (snRNA-seq) data to generate a reference expression profile and learn gene-specific bulk expression transformations to robustly decompose RNA-seq data. These transformations significantly improve decomposition performance compared to existing methods when there is significant technical variation in the generation of the reference profile and observed bulk expression. Importantly, compared to existing methods, our approach is extremely efficient, making it suitable for the analysis of large genomic datasets that are becoming ubiquitous. When applied to subcutaneous adipose and dorsolateral prefrontal cortex expression datasets with both bulk RNA-seq and snRNA-seq data, Bisque replicates previously reported associations between cell type proportions and measured phenotypes across abundant and rare cell types. We further propose an additional mode of operation that merely requires a set of known marker genes. Traditional methods for determining cell type composition lack scalability, while single-cell technologies remain costly and noisy compared to bulk RNA-seq. Here, the authors present a highly efficient tool to measure cellular heterogeneity in bulk expression through robust integration of single-cell information.

Genetic studies of metabolites have identified thousands of variants, many of which are associated with downstream metabolic and obesogenic disorders. However, these studies have relied on univariate analyses, reducing power and limiting context-specific understanding. Here we aim to provide an integrated perspective of the genetic basis of metabolites by leveraging the Finnish Metabolic Syndrome In Men (METSIM) cohort, a unique genetic resource which contains metabolic measurements, mostly lipids, across distinct time points as well as information on statin usage. We increase effective sample size by an average of two-fold by applying the Covariates for Multi-phenotype Studies (CMS) approach, identifying 588 significant SNP-metabolite associations, including 228 new associations. Our analysis pinpoints a small number of master metabolic regulator genes, balancing the relative proportion of dozens of metabolite levels. We further identify associations to changes in metabolic levels across time as well as genetic interactions with statin at both the master metabolic regulator and genome-wide level. Genome-wide association studies of metabolites have revealed hundreds of genetic associations using univariate analyses. Here, the authors use a multivariate approach to perform association analyses for 158 serum metabolites, followed by fine mapping and GxE interaction tests with statin use and age.

Single-nucleus RNA sequencing (snRNA-seq) measures gene expression in individual nuclei instead of cells, allowing for unbiased cell type characterization in solid tissues. We observe that snRNA-seq is commonly subject to contamination by high amounts of ambient RNA, which can lead to biased downstream analyses, such as identification of spurious cell types if overlooked. We present a novel approach to quantify contamination and filter droplets in snRNA-seq experiments, called Debris Identification using Expectation Maximization (DIEM). Our likelihood-based approach models the gene expression distribution of debris and cell types, which are estimated using EM. We evaluated DIEM using three snRNA-seq data sets: (1) human differentiating preadipocytes in vitro, (2) fresh mouse brain tissue, and (3) human frozen adipose tissue (AT) from six individuals. All three data sets showed evidence of extranuclear RNA contamination, and we observed that existing methods fail to account for contaminated droplets and led to spurious cell types. When compared to filtering using these state of the art methods, DIEM better removed droplets containing high levels of extranuclear RNA and led to higher quality clusters. Although DIEM was designed for snRNA-seq, our clustering strategy also successfully filtered single-cell RNA-seq data. To conclude, our novel method DIEM removes debris-contaminated droplets from single-cell-based data fast and effectively, leading to cleaner downstream analysis. Our code is freely available for use at https://github.com/marcalva/diem .

Recent and classical work has revealed biologically and medically significant subtypes in complex diseases and traits. However, relevant subtypes are often unknown, unmeasured, or actively debated, making automated statistical approaches to subtype definition valuable. We propose reverse GWAS (RGWAS) to identify and validate subtypes using genetics and multiple traits: while GWAS seeks the genetic basis of a given trait, RGWAS seeks to define trait subtypes with distinct genetic bases. Unlike existing approaches relying on off-the-shelf clustering methods, RGWAS uses a novel decomposition, MFMR, to model covariates, binary traits, and population structure. We use extensive simulations to show that modelling these features can be crucial for power and calibration. We validate RGWAS in practice by recovering a recently discovered stress subtype in major depression. We then show the utility of RGWAS by identifying three novel subtypes of metabolic traits. We biologically validate these metabolic subtypes with SNP-level tests and a novel polygenic test: the former recover known metabolic GxE SNPs; the latter suggests subtypes may explain substantial missing heritability. Crucially, statins, which are widely prescribed and theorized to increase diabetes risk, have opposing effects on blood glucose across metabolic subtypes, suggesting the subtypes have potential translational value.

Noncommunicable diseases (NCDs) are on the rise worldwide. Obesity, cardiovascular disease, and type 2 diabetes are among a long list of “lifestyle” diseases that were rare throughout human history but are now common. The evolutionary mismatch hypothesis posits that humans evolved in environments that radically differ from those we currently experience; consequently, traits that were once advantageous may now be “mismatched” and disease causing. At the genetic level, this hypothesis predicts that loci with a history of selection will exhibit “genotype by environment” (GxE) interactions, with different health effects in “ancestral” versus “modern” environments. To identify such loci, we advocate for combining genomic tools in partnership with subsistence-level groups experiencing rapid lifestyle change. In these populations, comparisons of individuals falling on opposite extremes of the “matched” to “mismatched” spectrum are uniquely possible. More broadly, the work we propose will inform our understanding of environmental and genetic risk factors for NCDs across diverse ancestries and cultures.

Insulin is an essential hormone that regulates glucose homeostasis and metabolism. Insulin resistance (IR) arises when tissues fail to respond to insulin, and it leads to serious health problems including Type 2 Diabetes (T2D). Obesity is a major contributor to the development of IR and T2D. We previously showed that gene expression of alcohol dehydrogenase 1B (ADH1B) was inversely correlated with obesity and IR in subcutaneous adipose tissue of Mexican Americans. In the current study, a meta-analysis of the relationship between ADH1B expression and BMI in Mexican Americans, African Americans, Europeans, and Pima Indians verified that BMI was increased with decreased ADH1B expression. Using established human subcutaneous pre-adipocyte cell lines derived from lean (BMI < 30 kg m −2 ) or obese (BMI ≥ 30 kg m −2 ) donors, we found that ADH1B protein expression increased substantially during differentiation, and overexpression of ADH1B inhibited fatty acid binding protein expression. Mature adipocytes from lean donors expressed ADH1B at higher levels than obese donors. Insulin further induced ADH1B protein expression as well as enzyme activity. Knockdown of ADH1B expression decreased insulin-stimulated glucose uptake. Our findings suggest that ADH1B is involved in the proper development and metabolic activity of adipose tissues and this function is suppressed by obesity.

Background & aims We have previously reported that the serum levels of gut-derived tryptophan metabolite indole-3-propionic acid (IPA) are lower in individuals with liver fibrosis. Now, we explored the transcriptome and DNA methylome associated with serum IPA levels in human liver from obese individuals together with IPA effects on shifting the hepatic stellate cell (HSC) phenotype to inactivation in vitro. Methods A total of 116 obese individuals without type 2 diabetes (T2D) (age 46.8 ± 9.3 years; BMI: 42.7 ± 5.0 kg/m 2 ) from the Kuopio OBesity Surgery (KOBS) study undergoing bariatric surgery were included. Circulating IPA levels were measured using LC–MS, liver transcriptomics with total RNA-sequencing and DNA methylation with Infinium HumanMethylation450 BeadChip. Human hepatic stellate cells (LX-2) where used for in vitro experiments. Results Serum IPA levels were associated with the expression of liver genes enriched for apoptosis, mitophagy and longevity pathways in the liver. AKT serine/threonine kinase 1 (AKT1) was the shared and topmost interactive gene from the liver transcript and DNA methylation profile. IPA treatment induced apoptosis, reduced mitochondrial respiration as well as modified cell morphology, and mitochondrial dynamics by modulating the expression of genes known to regulate fibrosis, apoptosis, and survival in LX-2 cells. Conclusion In conclusion, these data support that IPA has a plausible therapeutic effect and may induce apoptosis and the HSC phenotype towards the inactivation state, extending the possibilities to suppress hepatic fibrogenesis by interfering with HSC activation and mitochondrial metabolism. Graphical Abstract

We hypothesized that a common SNP in the 3' untranslated region of the upstream transcription factor 1 (USF1), rs3737787, may affect lipid traits by influencing gene expression levels, and we investigated this possibility utilizing the Mexican population, which has a high predisposition to dyslipidemia. We first associated rs3737787 genotypes in Mexican Familial Combined Hyperlipidemia (FCHL) case/control fat biopsies, with global expression patterns. To identify sets of co-expressed genes co-regulated by similar factors such as transcription factors, genetic variants, or environmental effects, we utilized weighted gene co-expression network analysis (WGCNA). Through WGCNA in the Mexican FCHL fat biopsies we identified two significant Triglyceride (TG)-associated co-expression modules. One of these modules was also associated with FCHL, the other FCHL component traits, and rs3737787 genotypes. This USF1-regulated FCHL-associated (URFA) module was enriched for genes involved in lipid metabolic processes. Using systems genetics procedures we identified 18 causal candidate genes in the URFA module. The FCHL causal candidate gene fatty acid desaturase 3 (FADS3) was associated with TGs in a recent Caucasian genome-wide significant association study and we replicated this association in Mexican FCHL families. Based on a USF1-regulated FCHL-associated co-expression module and SNP rs3737787, we identify a set of causal candidate genes for FCHL-related traits. We then provide evidence from two independent datasets supporting FADS3 as a causal gene for FCHL and elevated TGs in Mexicans.

Alterations in mitochondrial metabolism in obesity may indicate disrupted communication between mitochondria and nucleus, and DNA methylation may influence this interplay. Here, we leverage data from the Finnish Twin Cohort study subcohort ( n  = 173; 86 full twin pairs, 1 singleton), including comprehensive measurements of obesity-related outcomes, mitochondrial DNA quantity and nuclear DNA methylation levels in adipose and muscle tissue, to identify one CpG at SH3BP4 significantly associated with mitochondrial DNA quantity in adipose tissue (FDR < 0.05). We also show that SH3BP4 methylation correlates with its gene expression. Additionally, we find that 14 out of the 35 obesity-related traits display significant associations with both SH3BP4 methylation and mitochondrial DNA quantity in adipose tissue. We use data from TwinsUK and the Scandinavian T2D-discordant monozygotic twin cohort, to validate the observed associations. Further analysis using ICE FALCON suggests that mitochondrial DNA quantity, insulin sensitivity and certain body fat measures are causal to SH3BP4 methylation. Examining mitochondrial DNA quantity and obesity-related traits suggests causation from mitochondrial DNA quantity to obesity, but unmeasured within-individual confounding cannot be ruled out. Our findings underscore the impact of mitochondrial DNA quantity on DNA methylation and expression of the SH3BP4 gene within adipose tissue, with potential implications for obesity. In obesity, mitochondrial communication with the nucleus may be disrupted. Here, the authors use Finnish twin data to link obesity-related changes in mitochondrial DNA quantity with the methylation and expression of SH3BP4 in adipose tissue.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.An amendment to this paper has been published and can be accessed via a link at the top of the paper.