Genetic Regulation of Adipose Gene Expression and Cardio-Metabolic Traits

Systems genetics analysis of human body fat distribution genes identifies adipocyte processes

Excess abdominal fat is a sexually dimorphic risk factor for cardio-metabolic disease and is approximated by the waist-to-hip ratio adjusted for body mass index (WHRadjBMI). Whereas this trait is highly heritable, few causal genes are known. We aimed to identify novel drivers of WHRadjBMI using systems genetics. We used two independent cohorts of adipose tissue gene expression and constructed sex- and depot-specific Bayesian networks to model gene-gene interactions from 8,492 genes. Using key driver analysis, we identified genes that, in silico and putatively in vitro, regulate many others. 51-119 key drivers in each network were replicated in both cohorts. In other cell types, 23 of these genes are found in crucial adipocyte pathways: Wnt signaling or mitochondrial function. We overexpressed or down-regulated seven key driver genes in human subcutaneous pre-adipocytes. Key driver genes ANAPC2 and RSPO1 inhibited adipogenesis, whereas PSME3 increased adipogenesis. RSPO1 increased Wnt signaling activity. In differentiated adipocytes, MIGA1 and UBR1 down-regulation led to mitochondrial dysfunction. These five genes regulate adipocyte function, and we hypothesize that they regulate fat distribution.

Obesity causes mitochondrial fragmentation and dysfunction in white adipocytes due to RalA activation

Mitochondrial dysfunction is a characteristic trait of human and rodent obesity, insulin resistance and fatty liver disease. Here we show that high-fat diet (HFD) feeding causes mitochondrial fragmentation in inguinal white adipocytes from male mice, leading to reduced oxidative capacity by a process dependent on the small GTPase RalA. RalA expression and activity are increased in white adipocytes after HFD. Targeted deletion of RalA in white adipocytes prevents fragmentation of mitochondria and diminishes HFD-induced weight gain by increasing fatty acid oxidation. Mechanistically, RalA increases fission in adipocytes by reversing the inhibitory Ser637 phosphorylation of the fission protein Drp1, leading to more mitochondrial fragmentation. Adipose tissue expression of the human homolog of Drp1, DNM1L, is positively correlated with obesity and insulin resistance. Thus, chronic activation of RalA plays a key role in repressing energy expenditure in obese adipose tissue by shifting the balance of mitochondrial dynamics toward excessive fission, contributing to weight gain and metabolic dysfunction.

Kruppel-like factor 14 ameliorated obesity and related metabolic disorders by promoting adipose tissue browning

Obesity has been identified as a serious and debilitating disease that threatens human health, but the current treatment strategies still have some shortcomings. Exercise and dieting are difficult for many people to adhere to, and a series of surgical risks and pain brought about by volume reduction have made it difficult for the current weight loss effect to meet human expectations. In this study, we first found that mice with overexpression of the transcription factor Kruppel-like factor 14 (KLF14) in subcutaneous adipose tissue gained weight more slowly while consuming a high-fat diet than did control mice, and these mice also showed reduced insulin resistance and liver lipid deposition abnormalities. Mechanistically, the browning of white adipose tissue was promoted in adipose tissue with KLF14 overexpression; therefore, we preliminarily concluded that KLF14 can improve obesity by promoting the browning of white adipose tissue and energy consumption, thus ameliorating obesity and related metabolic disturbances. In summary, our results revealed that KLF14 may promote white adipose tissue browning, thus ameliorating high-fat diet-induced obesity and hepatic steatosis, as well as serum lipid levels and insulin resistance, thereby achieving a positive effect on metabolism.NEW & NOTEWORTHY Our study first explored the role of KLF14 in the development and progression of HFD-induced obesity in male mice. Its beneficial effect on adipose browning and metabolic disorders suggests that KLF14 may provide us a new therapeutic strategy for obesity and related metabolic complications. This health problem is of global concern and needs to be addressed.

CTRP13 ablation improves systemic glucose and lipid metabolism

OBJECTIVE

Tissue crosstalk mediated by secreted hormones underlies the integrative control of metabolism. We previously showed that CTRP13/C1QL3, a secreted protein of the C1q family, can improve glucose metabolism and insulin action in vitro and reduce food intake and body weight in mice when centrally delivered. A role for CTRP13 in regulating insulin secretion in isolated islets has also been demonstrated. It remains unclear, however, whether the effects of CTRP13 on cultured cells and in mice reflect the physiological function of the protein. Here, we use a loss-of-function mouse model to address whether CTRP13 is required for metabolic homeostasis.

METHODS

WT and Ctrp13 knockout (KO) mice fed a standard chow or a high-fat diet were subjected to comprehensive metabolic phenotyping. Transcriptomic analyses were carried out on visceral and subcutaneous fat, liver, and skeletal muscle to identify pathways altered by CTRP13 deficiency. RNA-seq data was further integrated with the Metabolic Syndrome in Man (METSIM) cohort data. Adjusted regression analysis was used to demonstrate that genetic variation of CTRP13 expression accounts for a significant proportion of variance between differentially expressed genes (DEGs) in adipose tissue and metabolic traits in humans.

RESULTS

Contrary to expectation, chow-fed Ctrp13-KO male mice had elevated physical activity, lower body weight, and improved lipid handling. On a high-fat diet (HFD), Ctrp13-KO mice of either sex were consistently more active and leaner. Loss of CTRP13 reduced hepatic glucose output and improved glucose tolerance, insulin sensitivity, and triglyceride clearance, though with notable sex differences. Consistent with the lean phenotype, transcriptomic analyses revealed a lower inflammatory profile in visceral fat and liver. Reduced hepatic steatosis was correlated with the suppression of lipid synthesis and enhanced lipid catabolism gene expression. Visceral fat had the largest number of DEGs and mediation analyses on the human orthologs of the DEGs suggested the potential causal contribution of CTRP13 to human metabolic syndrome.

CONCLUSIONS

Our results suggest that CTRP13 is a negative metabolic regulator, and its deficiency improves systemic metabolic profiles. Our data also suggest the reduction in circulating human CTRP13 levels seen in obesity and diabetes may reflect a compensatory physiologic response to counteract insulin resistance.

HSP47 levels determine the degree of body adiposity

Adiposity varies among individuals with the influence of diverse physiological, pathological, environmental, hormonal, and genetic factors, but a unified molecular basis remains elusive. Here, we identify HSP47, a collagen-specific chaperone, as a key determinant of body adiposity. HSP47 expression is abundant in adipose tissue; increased with feeding, overeating, and obesity; decreased with fasting, exercise, calorie restriction, bariatric surgery, and cachexia; and correlated with fat mass, BMI, waist, and hip circumferences. Insulin and glucocorticoids, respectively, up- and down-regulate HSP47 expression. In humans, the increase of HSP47 gene expression by its intron or synonymous variants is associated with higher body adiposity traits. In mice, the adipose-specific knockout or pharmacological inhibition of HSP47 leads to lower body adiposity compared to the control. Mechanistically, HSP47 promotes collagen dynamics in the folding, secretion, and interaction with integrin, which activates FAK signaling and preserves PPARγ protein from proteasomal degradation, partly related to MDM2. The study highlights the significance of HSP47 in determining the amount of body fat individually and under various circumstances.

Circular RNA ZNF800 (hsa_circ_0082096) regulates cancer stem cell properties and tumor growth in colorectal cancer

BACKGROUND

Cancer stem cells form a rare cell population in tumors that contributes to metastasis, recurrence and chemoresistance in cancer patients. Circular RNAs (circRNAs) are post-transcriptional regulators of gene expression that sponge targeted microRNA (miRNAs) to affect a multitude of downstream cellular processes. We previously showed in an expression profiling study that circZNF800 (hsa_circ_0082096) was up-regulated in cancer stem cell-enriched spheroids derived from colorectal cancer (CRC) cell lines.

METHODS

Spheroids were generated in suspension spheroidal culture. The ZNF800 mRNA, pluripotency stem cell markers and circZNF800 levels were determined by quantitative RT-PCR. CircZNF800-miRNA interactions were shown in RNA pulldown assays and the miRNA levels determined by stem-loop qRT-PCR. The effects of circZNF800 on cell proliferation were tested by EdU staining followed by flowcytometry. Expression of stem cell markers CD44/CD133, Lgr5 and SOX9 was demonstrated in immunofluorescence microscopy. To manipulate the cellular levels of circZNF800, circZNF800 over-expression was achieved via transfection of in vitro synthesized and circularized circZNF800, and knockdown attained using a CRISPR-Cas13d-circZNF800 vector system. Xenografted nude mice were used to demonstrate effects of circZNF800 over-expression and knockdown on tumor growth in vivo.

RESULTS

CircZNF800 was shown to be over-expressed in late-stage tumor tissues of CRC patients. Data showed that circZNF800 impeded expression of miR-140-3p, miR-382-5p and miR-579-3p while promoted the mRNA levels of ALK/ACVR1C, FZD3 and WNT5A targeted by the miRNAs, as supported by alignments of seed sequences between the circZNF800-miRNA, and miRNA-mRNA paired interactions. Analysis in CRC cells and biopsied tissues showed that circZNF800 positively regulated the expression of intestinal stem cell, pluripotency and cancer stem cell markers, and promoted CRC cell proliferation, spheroid and colony formation in vitro, all of which are cancer stem cell properties. In xenografted mice, circZNF800 over-expression promoted tumor growth, while circZNF800 knockdown via administration of CRISPR Cas13d-circZNF800 viral particles at the CRC tumor sites impeded tumor growth.

CONCLUSIONS

CircZNF800 is an oncogenic factor that regulate cancer stem cell properties to lead colorectal tumorigenesis, and may be used as a predictive marker for tumor progression and the CRISPR Cas13d-circZNF800 knockdown strategy for therapeutic intervention of colorectal cancer.

Cell-Type Composition Affects Adipose Gene Expression Associations With Cardiometabolic Traits

Understanding differences in adipose gene expression between individuals with different levels of clinical traits may reveal the genes and mechanisms leading to cardiometabolic diseases. However, adipose is a heterogeneous tissue. To account for cell-type heterogeneity, we estimated cell-type proportions in 859 subcutaneous adipose tissue samples with bulk RNA sequencing (RNA-seq) using a reference single-nuclear RNA-seq data set. Cell-type proportions were associated with cardiometabolic traits; for example, higher macrophage and adipocyte proportions were associated with higher and lower BMI, respectively. We evaluated cell-type proportions and BMI as covariates in tests of association between >25,000 gene expression levels and 22 cardiometabolic traits. For >95% of genes, the optimal, or best-fit, models included BMI as a covariate, and for 79% of associations, the optimal models also included cell type. After adjusting for the optimal covariates, we identified 2,664 significant associations (P ≤ 2e-6) for 1,252 genes and 14 traits. Among genes proposed to affect cardiometabolic traits based on colocalized genome-wide association study and adipose expression quantitative trait locus signals, 25 showed a corresponding association between trait and gene expression levels. Overall, these results suggest the importance of modeling cell-type proportion when identifying gene expression associations with cardiometabolic traits.

ARTICLE HIGHLIGHTS

Adipose tissue eQTL meta-analysis reveals the contribution of allelic heterogeneity to gene expression regulation and cardiometabolic traits

Complete characterization of the genetic effects on gene expression is needed to elucidate tissue biology and the etiology of complex traits. Here, we analyzed 2,344 subcutaneous adipose tissue samples and identified 34K conditionally distinct expression quantitative trait locus (eQTL) signals in 18K genes. Over half of eQTL genes exhibited at least two eQTL signals. Compared to primary signals, non-primary signals had lower effect sizes, lower minor allele frequencies, and less promoter enrichment; they corresponded to genes with higher heritability and higher tolerance for loss of function. Colocalization of eQTL with conditionally distinct genome-wide association study signals for 28 cardiometabolic traits identified 3,605 eQTL signals for 1,861 genes. Inclusion of non-primary eQTL signals increased colocalized signals by 46%. Among 30 genes with ≥2 pairs of colocalized signals, 21 showed a mediating gene dosage effect on the trait. Thus, expanded eQTL identification reveals more mechanisms underlying complex traits and improves understanding of the complexity of gene expression regulation.

A brown fat-enriched adipokine, ASRA, is a leptin receptor antagonist that stimulates appetite

The endocrine control of food intake remains incompletely understood, and whether the leptin receptor-mediated anorexigenic pathway in the hypothalamus is negatively regulated by a humoral factor is unknown. Here we identify an appetite-stimulating factor - ASRA - that acts as a leptin receptor antagonist. ASRA encodes an 8 kD protein that is abundantly and selectively expressed in adipose tissue and to a lesser extent, in liver, and is upregulated during fasting and cold. ASRA protein associates with autophagosomes and its secretion is induced by energy deficiency. Overexpression of ASRA in mice attenuates leptin receptor signaling leading to elevated blood glucose and development of severe hyperphagic obesity, whereas either adipose- or liver-specific ASRA knockout mice display increased leptin sensitivity, improved glucose homeostasis, reduced food intake, and resistance to high fat diet-induced obesity. Furthermore, ASRA is indispensable for cold-evoked feeding response. Recombinant ASRA (rASRA) protein binds to leptin receptor and suppresses leptin receptor signaling in cultured cells. In vivo, rASRA promotes food intake and increases blood glucose in a leptin receptor signaling-dependent manner. Our studies collectively show that ASRA, acting as a peripheral signal of energy deficit, stimulates appetite and regulates glucose metabolism by antagonizing leptin receptor signaling, thus revealing a previously unknown endocrine mechanism that has important implications for our understanding of leptin resistance.

Systems genetics analysis of human body fat distribution genes identifies Wnt signaling and mitochondrial activity in adipocytes

BACKGROUND

Excess fat in the abdomen is a sexually dimorphic risk factor for cardio-metabolic disease. The relative storage between abdominal and lower-body subcutaneous adipose tissue depots is approximated by the waist-to-hip ratio adjusted for body mass index (WHRadjBMI). Genome-wide association studies (GWAS) identified 346 loci near 495 genes associated with WHRadjBMI. Most of these genes have unknown roles in fat distribution, but many are expressed and putatively act in adipose tissue. We aimed to identify novel sex- and depot-specific drivers of WHRadjBMI using a systems genetics approach.

METHODS

We used two independent cohorts of adipose tissue gene expression with 362 - 444 males and 147 - 219 females, primarily of European ancestry. We constructed sex- and depot- specific Bayesian networks to model the gene-gene interactions from 8,492 adipose tissue genes. Key driver analysis identified genes that, in silico and putatively in vitro, regulate many others, including the 495 WHRadjBMI GWAS genes. Key driver gene function was determined by perturbing their expression in human subcutaneous pre-adipocytes using lenti-virus or siRNA.

RESULTS

51 - 119 key drivers in each network were replicated in both cohorts. We used single-cell expression data to select replicated key drivers expressed in adipocyte precursors and mature adipocytes, prioritized genes which have not been previously studied in adipose tissue, and used public human and mouse data to nominate 53 novel key driver genes (10 - 21 from each network) that may regulate fat distribution by altering adipocyte function. In other cell types, 23 of these genes are found in crucial adipocyte pathways: Wnt signaling or mitochondrial function. We selected seven genes whose expression is highly correlated with WHRadjBMI to further study their effects on adipogenesis/Wnt signaling (ANAPC2, PSME3, RSPO1, TYRO3) or mitochondrial function (C1QTNF3, MIGA1, PSME3, UBR1).Adipogenesis was inhibited in cells overexpressing ANAPC2 and RSPO1 compared to controls. RSPO1 results are consistent with a positive correlation between gene expression in the subcutaneous depot and WHRadjBMI, therefore lower relative storage in the subcutaneous depot. RSPO1 inhibited adipogenesis by increasing β-catenin activation and Wnt-related transcription, thus repressing PPARG and CEBPA. PSME3 overexpression led to more adipogenesis than controls. In differentiated adipocytes, MIGA1 and UBR1 downregulation led to mitochondrial dysfunction, with lower oxygen consumption than controls; MIGA1 knockdown also lowered UCP1 expression.

SUMMARY

ANAPC2, MIGA1, PSME3, RSPO1, and UBR1 affect adipocyte function and may drive body fat distribution.

Genetics of cell-type-specific post-transcriptional gene regulation during human neurogenesis

The function of some genetic variants associated with brain-relevant traits has been explained through colocalization with expression quantitative trait loci (eQTL) conducted in bulk post-mortem adult brain tissue. However, many brain-trait associated loci have unknown cellular or molecular function. These genetic variants may exert context-specific function on different molecular phenotypes including post-transcriptional changes. Here, we identified genetic regulation of RNA-editing and alternative polyadenylation (APA), within a cell-type-specific population of human neural progenitors and neurons. More RNA-editing and isoforms utilizing longer polyadenylation sequences were observed in neurons, likely due to higher expression of genes encoding the proteins mediating these post-transcriptional events. We also detected hundreds of cell-type-specific editing quantitative trait loci (edQTLs) and alternative polyadenylation QTLs (apaQTLs). We found colocalizations of a neuron edQTL in CCDC88A with educational attainment and a progenitor apaQTL in EP300 with schizophrenia, suggesting genetically mediated post-transcriptional regulation during brain development lead to differences in brain function.

Identifying novel regulatory effects for clinically relevant genes through the study of the Greek population

BACKGROUND

Expression quantitative trait loci (eQTL) studies provide insights into regulatory mechanisms underlying disease risk. Expanding studies of gene regulation to underexplored populations and to medically relevant tissues offers potential to reveal yet unknown regulatory variants and to better understand disease mechanisms. Here, we performed eQTL mapping in subcutaneous (S) and visceral (V) adipose tissue from 106 Greek individuals (Greek Metabolic study, GM) and compared our findings to those from the Genotype-Tissue Expression (GTEx) resource.

RESULTS

We identified 1,930 and 1,515 eGenes in S and V respectively, over 13% of which are not observed in GTEx adipose tissue, and that do not arise due to different ancestry. We report additional context-specific regulatory effects in genes of clinical interest (e.g. oncogene ST7) and in genes regulating responses to environmental stimuli (e.g. MIR21, SNX33). We suggest that a fraction of the reported differences across populations is due to environmental effects on gene expression, driving context-specific eQTLs, and suggest that environmental effects can determine the penetrance of disease variants thus shaping disease risk. We report that over half of GM eQTLs colocalize with GWAS SNPs and of these colocalizations 41% are not detected in GTEx. We also highlight the clinical relevance of S adipose tissue by revealing that inflammatory processes are upregulated in individuals with obesity, not only in V, but also in S tissue.

CONCLUSIONS

By focusing on an understudied population, our results provide further candidate genes for investigation regarding their role in adipose tissue biology and their contribution to disease risk and pathogenesis.

High-throughput functional dissection of noncoding SNPs with biased allelic enhancer activity for insulin resistance-relevant phenotypes

Most of the single-nucleotide polymorphisms (SNPs) associated with insulin resistance (IR)-relevant phenotypes by genome-wide association studies (GWASs) are located in noncoding regions, complicating their functional interpretation. Here, we utilized an adapted STARR-seq to evaluate the regulatory activities of 5,987 noncoding SNPs associated with IR-relevant phenotypes. We identified 876 SNPs with biased allelic enhancer activity effects (baaSNPs) across 133 loci in three IR-relevant cell lines (HepG2, preadipocyte, and A673), which showed pervasive cell specificity and significant enrichment for cell-specific open chromatin regions or enhancer-indicative markers (H3K4me1, H3K27ac). Further functional characterization suggested several transcription factors (TFs) with preferential allelic binding to baaSNPs. We also incorporated multi-omics data to prioritize 102 candidate regulatory target genes for baaSNPs and revealed prevalent long-range regulatory effects and cell-specific IR-relevant biological functional enrichment on them. Specifically, we experimentally verified the distal regulatory mechanism at IRS1 locus, in which rs952227-A reinforces IRS1 expression by long-range chromatin interaction and preferential binding to the transcription factor HOXC6 to augment the enhancer activity. Finally, based on our STARR-seq screening data, we predicted the enhancer activity of 227,343 noncoding SNPs associated with IR-relevant phenotypes (fasting insulin adjusted for BMI, HDL cholesterol, and triglycerides) from the largest available GWAS summary statistics. We further provided an open resource (http://www.bigc.online/fnSNP-IR) for better understanding genetic regulatory mechanisms of IR-relevant phenotypes.

Breaking enhancers to gain insights into developmental defects

Despite ground-breaking genetic studies that have identified thousands of risk variants for developmental diseases, how these variants lead to molecular and cellular phenotypes remains a gap in knowledge. Many of these variants are non-coding and occur at enhancers, which orchestrate key regulatory programs during development. The prevailing paradigm is that non-coding variants alter the activity of enhancers, impacting gene expression programs, and ultimately contributing to disease risk. A key obstacle to progress is the systematic functional characterization of non-coding variants at scale, especially since enhancer activity is highly specific to cell type and developmental stage. Here, we review the foundational studies of enhancers in developmental disease and current genomic approaches to functionally characterize developmental enhancers and their variants at scale. In the coming decade, we anticipate systematic enhancer perturbation studies to link non-coding variants to molecular mechanisms, changes in cell state, and disease phenotypes.

Predicting mechanisms of action at genetic loci associated with discordant effects on type 2 diabetes and abdominal fat accumulation

Obesity is a major risk factor for cardiovascular disease, stroke, and type 2 diabetes (T2D). Excessive accumulation of fat in the abdomen further increases T2D risk. Abdominal obesity is measured by calculating the ratio of waist-to-hip circumference adjusted for the body-mass index (WHRadjBMI), a trait with a significant genetic inheritance. Genetic loci associated with WHRadjBMI identified in genome-wide association studies are predicted to act through adipose tissues, but many of the exact molecular mechanisms underlying fat distribution and its consequences for T2D risk are poorly understood. Further, mechanisms that uncouple the genetic inheritance of abdominal obesity from T2D risk have not yet been described. Here we utilize multi-omic data to predict mechanisms of action at loci associated with discordant effects on abdominal obesity and T2D risk. We find six genetic signals in five loci associated with protection from T2D but also with increased abdominal obesity. We predict the tissues of action at these discordant loci and the likely effector Genes (eGenes) at three discordant loci, from which we predict significant involvement of adipose biology. We then evaluate the relationship between adipose gene expression of eGenes with adipogenesis, obesity, and diabetic physiological phenotypes. By integrating these analyses with prior literature, we propose models that resolve the discordant associations at two of the five loci. While experimental validation is required to validate predictions, these hypotheses provide potential mechanisms underlying T2D risk stratification within abdominal obesity.

Obesity-dependent increase in RalA activity disrupts mitochondrial dynamics in white adipocytes

Mitochondrial dysfunction is a characteristic trait of human and rodent obesity, insulin resistance, and fatty liver disease. Here we report that mitochondria undergo fragmentation and reduced oxidative capacity specifically in inguinal white adipose tissue after feeding mice high fat diet (HFD) by a process dependent on the small GTPase RalA. RalA expression and activity are increased in white adipocytes from mice fed HFD. Targeted deletion of Rala in white adipocytes prevents the obesity-induced fragmentation of mitochondria and produces mice resistant to HFD-induced weight gain via increased fatty acid oxidation. As a result, these mice also exhibit improved glucose tolerance and liver function. In vitro mechanistic studies revealed that RalA suppresses mitochondrial oxidative function in adipocytes by increasing fission through reversing the protein kinase A-catalyzed inhibitory Ser637phosphorylation of the mitochondrial fission protein Drp1. Active RalA recruits protein phosphatase 2A (PP2Aa) to specifically dephosphorylate this inhibitory site on Drp1, activating the protein, thus increasing mitochondrial fission. Adipose tissue expression of the human homolog of Drp1, DNML1, is positively correlated with obesity and insulin resistance in patients. Thus, chronic activation of RalA plays a key role in repressing energy expenditure in obese adipose tissue by shifting the balance of mitochondrial dynamics towards excessive fission, contributing to weight gain and related metabolic dysfunction.

Inferring cell-type-specific causal gene regulatory networks during human neurogenesis

BACKGROUND

Genetic variation influences both chromatin accessibility, assessed in chromatin accessibility quantitative trait loci (caQTL) studies, and gene expression, assessed in expression QTL (eQTL) studies. Genetic variants can impact either nearby genes (cis-eQTLs) or distal genes (trans-eQTLs). Colocalization between caQTL and eQTL, or cis- and trans-eQTLs suggests that they share causal variants. However, pairwise colocalization between these molecular QTLs does not guarantee a causal relationship. Mediation analysis can be applied to assess the evidence supporting causality versus independence between molecular QTLs. Given that the function of QTLs can be cell-type-specific, we performed mediation analyses to find epigenetic and distal regulatory causal pathways for genes within two major cell types of the developing human cortex, progenitors and neurons.

RESULTS

We find that the expression of 168 and 38 genes is mediated by chromatin accessibility in progenitors and neurons, respectively. We also find that the expression of 11 and 12 downstream genes is mediated by upstream genes in progenitors and neurons. Moreover, we discover that a genetic locus associated with inter-individual differences in brain structure shows evidence for mediation of SLC26A7 through chromatin accessibility, identifying molecular mechanisms of a common variant association to a brain trait.

CONCLUSIONS

In this study, we identify cell-type-specific causal gene regulatory networks whereby the impacts of variants on gene expression were mediated by chromatin accessibility or distal gene expression. Identification of these causal paths will enable identifying and prioritizing actionable regulatory targets perturbing these key processes during neurodevelopment.

The potential of integrating human and mouse discovery platforms to advance our understanding of cardiometabolic diseases

Cardiometabolic diseases encompass a range of interrelated conditions that arise from underlying metabolic perturbations precipitated by genetic, environmental, and lifestyle factors. While obesity, dyslipidaemia, smoking, and insulin resistance are major risk factors for cardiometabolic diseases, individuals still present in the absence of such traditional risk factors, making it difficult to determine those at greatest risk of disease. Thus, it is crucial to elucidate the genetic, environmental, and molecular underpinnings to better understand, diagnose, and treat cardiometabolic diseases. Much of this information can be garnered using systems genetics, which takes population-based approaches to investigate how genetic variance contributes to complex traits. Despite the important advances made by human genome-wide association studies (GWAS) in this space, corroboration of these findings has been hampered by limitations including the inability to control environmental influence, limited access to pertinent metabolic tissues, and often, poor classification of diseases or phenotypes. A complementary approach to human GWAS is the utilisation of model systems such as genetically diverse mouse panels to study natural genetic and phenotypic variation in a controlled environment. Here, we review mouse genetic reference panels and the opportunities they provide for the study of cardiometabolic diseases and related traits. We discuss how the post-GWAS era has prompted a shift in focus from discovery of novel genetic variants to understanding gene function. Finally, we highlight key advantages and challenges of integrating complementary genetic and multi-omics data from human and mouse populations to advance biological discovery.

Long-read proteogenomics to connect disease-associated sQTLs to the protein isoform effectors of disease

A major fraction of loci identified by genome-wide association studies (GWASs) lead to alterations in alternative splicing, but interpretation of how such alterations impact proteins is hindered by the technical limitations of short-read RNA-seq, which cannot directly link splicing events to full-length transcript or protein isoforms. Long-read RNA-seq represents a powerful tool to define and quantify transcript isoforms, and recently, infer protein isoform existence. Here we present a novel approach that integrates information from GWAS, splicing QTL (sQTL), and PacBio long-read RNA-seq in a disease-relevant model to infer the effects of sQTLs on the ultimate protein isoform products they encode. We demonstrate the utility of our approach using bone mineral density (BMD) GWAS data. We identified 1,863 sQTLs from the Genotype-Tissue Expression (GTEx) project in 732 protein-coding genes which colocalized with BMD associations (H 4 PP ≥ 0.75). We generated deep coverage PacBio long-read RNA-seq data (N=∼22 million full-length reads) on human osteoblasts, identifying 68,326 protein-coding isoforms, of which 17,375 (25%) were novel. By casting the colocalized sQTLs directly onto protein isoforms, we connected 809 sQTLs to 2,029 protein isoforms from 441 genes expressed in osteoblasts. Using these data, we created one of the first proteome-scale resources defining full-length isoforms impacted by colocalized sQTLs. Overall, we found that 74 sQTLs influenced isoforms likely impacted by nonsense mediated decay (NMD) and 190 that potentially resulted in the expression of new protein isoforms. Finally, we identified colocalizing sQTLs in TPM2 for splice junctions between two mutually exclusive exons, and two different transcript termination sites, making it impossible to interpret without long-read RNA-seq data. siRNA mediated knockdown in osteoblasts showed two TPM2 isoforms with opposing effects on mineralization. We expect our approach to be widely generalizable across diverse clinical traits and accelerate system-scale analyses of protein isoform activities modulated by GWAS loci.

Antimicrobial protein REG3A regulates glucose homeostasis and insulin resistance in obese diabetic mice

Innate immune mediators of pathogen clearance, including the secreted C-type lectins REG3 of the antimicrobial peptide (AMP) family, are known to be involved in the regulation of tissue repair and homeostasis. Their role in metabolic homeostasis remains unknown. Here we show that an increase in human REG3A improves glucose and lipid homeostasis in nutritional and genetic mouse models of obesity and type 2 diabetes. Mice overexpressing REG3A in the liver show improved glucose homeostasis, which is reflected in better insulin sensitivity in normal weight and obese states. Delivery of recombinant REG3A protein to leptin-deficient ob/ob mice or wild-type mice on a high-fat diet also improves glucose homeostasis. This is accompanied by reduced oxidative protein damage, increased AMPK phosphorylation and insulin-stimulated glucose uptake in skeletal muscle tissue. Oxidative damage in differentiated C2C12 myotubes is greatly attenuated by REG3A, as is the increase in gp130-mediated AMPK activation. In contrast, Akt-mediated insulin action, which is impaired by oxidative stress, is not restored by REG3A. These data highlight the importance of REG3A in controlling oxidative protein damage involved in energy and metabolic pathways during obesity and diabetes, and provide additional insight into the dual function of host-immune defense and metabolic regulation for AMP.

Genetic Regulation of SMC Gene Expression and Splicing Predict Causal CAD Genes

BACKGROUND

Coronary artery disease (CAD) is the leading cause of death worldwide. Recent meta-analyses of genome-wide association studies have identified over 175 loci associated with CAD. The majority of these loci are in noncoding regions and are predicted to regulate gene expression. Given that vascular smooth muscle cells (SMCs) play critical roles in the development and progression of CAD, we aimed to identify the subset of the CAD loci associated with the regulation of transcription in distinct SMC phenotypes.

METHODS

We measured gene expression in SMCs isolated from the ascending aortas of 151 heart transplant donors of various genetic ancestries in quiescent or proliferative conditions and calculated the association of their expression and splicing with ~6.3 million imputed single-nucleotide polymorphism markers across the genome.

RESULTS

We identified 4910 expression and 4412 splicing quantitative trait loci (sQTLs) representing regions of the genome associated with transcript abundance and splicing. A total of 3660 expression quantitative trait loci (eQTLs) had not been observed in the publicly available Genotype-Tissue Expression dataset. Further, 29 and 880 eQTLs were SMC-specific and sex-biased, respectively. We made these results available for public query on a user-friendly website. To identify the effector transcript(s) regulated by CAD loci, we used 4 distinct colocalization approaches. We identified 84 eQTL and 164 sQTL that colocalized with CAD loci, highlighting the importance of genetic regulation of mRNA splicing as a molecular mechanism for CAD genetic risk. Notably, 20% and 35% of the eQTLs were unique to quiescent or proliferative SMCs, respectively. One CAD locus colocalized with a sex-specific eQTL (TERF2IP), and another locus colocalized with SMC-specific eQTL (ALKBH8). The most significantly associated CAD locus, 9p21, was an sQTL for the long noncoding RNA CDKN2B-AS1, also known as ANRIL, in proliferative SMCs.

CONCLUSIONS

Collectively, our results provide evidence for the molecular mechanisms of genetic susceptibility to CAD in distinct SMC phenotypes.

Biochemical Activation and Regulatory Functions of Trans-Regulatory KLF14 and Its Association with Genetic Polymorphisms

Krüpple-Like family of transcription factor-14 (KLF14) is a master trans-regulatory gene that has multiple biological regulatory functions and is involved in many pathological mechanisms. It controls the expressions of several other genes which are involved in multiple regulatory functions. KLF14 plays a significant role in lipid metabolism, glucose regulation and insulin sensitivity. Cell apoptosis, proliferation, and differentiation are regulated by the KLF14 gene, and up-regulation of KLF14 prevents cancer progression. KLF14 has been used as an epigenetic biomarker for the estimation of chronological age due to the presence of different age-related CpG sites on genes that become methylated with age. Different genome-wide association studies have identified several KLF14 variants in adipose tissues. These single nucleotide polymorphisms in KLF14 have been associated with dyslipidemia, insulin resistance, and glucose intolerance. Moreover, the prevalence of genetic polymorphism is different in different populations due to ethnic differences and epigenetic modifications. In addition, environmental and physiological factors such as diet, age, gender, and obesity are also responsible for genetic mutations in KLF14.

Multi-omic integration reveals cell-type-specific regulatory networks of insulin resistance in distinct ancestry populations

Our knowledge of the cell-type-specific mechanisms of insulin resistance remains limited. To dissect the cell-type-specific molecular signatures of insulin resistance, we performed a multiscale gene network analysis of adipose and muscle tissues in African and European ancestry populations. In adipose tissues, a comparative analysis revealed ethnically conserved cell-type signatures and two adipocyte subtype-enriched modules with opposite insulin sensitivity responses. The modules enriched for adipose stem and progenitor cells as well as immune cells showed negative correlations with insulin sensitivity. In muscle tissues, the modules enriched for stem cells and fibro-adipogenic progenitors responded to insulin sensitivity oppositely. The adipocyte and muscle fiber-enriched modules shared cellular-respiration-related genes but had tissue-specific rearrangements of gene regulations in response to insulin sensitivity. Integration of the gene co-expression and causal networks further pinpointed key drivers of insulin resistance. Together, this study revealed the cell-type-specific transcriptomic networks and signaling maps underlying insulin resistance in major glucose-responsive tissues. A record of this paper's transparent peer review process is included in the supplemental information.

A brown fat-enriched adipokine Adissp controls adipose thermogenesis and glucose homeostasis

The signaling mechanisms underlying adipose thermogenesis have not been fully elucidated. Particularly, the involvement of adipokines that are selectively expressed in brown adipose tissue (BAT) and beige adipocytes remains to be investigated. Here we show that a previously uncharacterized adipokine (UPF0687 protein / human C20orf27 homolog) we named as Adissp (Adipose-secreted signaling protein) is a key regulator for white adipose tissue (WAT) thermogenesis and glucose homeostasis. Adissp expression is adipose-specific and highly BAT-enriched, and its secretion is stimulated by β3-adrenergic activation. Gain-of-functional studies collectively showed that secreted Adissp promotes WAT thermogenesis, improves glucose homeostasis, and protects against obesity. Adipose-specific Adissp knockout mice are defective in WAT browning, and are susceptible to high fat diet-induced obesity and hyperglycemia. Mechanistically, Adissp binds to a putative receptor on adipocyte surface and activates protein kinase A independently of β-adrenergic signaling. These results establish BAT-enriched Adissp as a major upstream signaling component in thermogenesis and offer a potential avenue for the treatment of obesity and diabetes.

Obesity, Fat Distribution and Risk of Cancer in Women and Men: A Mendelian Randomisation Study

Obesity and upper-body fat distribution are independent, cardiometabolic risk factors but whether they also display comparable associations with cancer risk is unknown. We investigated the causal relationships between body mass index (BMI) and BMI-adjusted waist-to-hip ratio (WHRadjBMI) and cancer risk and searched for potential drivers linking these traits to carcinogenesis using two-sample and multivariable Mendelian randomisation. In women, genetically instrumented higher BMI was associated with lower breast (OR = 0.87, 95% CI 0.81-0.93) and higher endometrial (OR = 1.75, 95% CI 1.55-1.96) cancer risk whilst WHRadjBMI was associated with higher colon cancer risk (OR = 1.22, 95% CI 1.07-1.42). In men, elevated BMI was associated with lower prostate cancer risk (OR = 0.91, 95% CI 0.85-0.98). Mechanistically, testosterone and insulin mediated 21% and 35%, respectively of the total, genetically determined association of BMI with endometrial cancer risk whilst HDL cholesterol and IGF-1 mediated 40% and 22%, respectively of the association between BMI and breast cancer risk. In men, testosterone mediated 21% of the association between BMI and prostate cancer risk. Colon cancer aside, the total amount of body fat might be more important than its location in modulating cancer susceptibility due to differential effects of obesity and fat distribution on adiposity-associated cancer drivers.

Lineage-selective super enhancers mediate core regulatory circuitry during adipogenic and osteogenic differentiation of human mesenchymal stem cells

Human mesenchymal stem cells (hMSCs) can be differentiated into osteoblasts and adipocytes. During these processes, super enhancers (SEs) play important roles. Here, we performed comprehensive characterization of the SEs changes associated with adipogenic and osteogenic differentiation of hMSCs, and revealed that SEs changed more dramatically compared with typical enhancers. We identified a set of lineage-selective SEs, whose target genes were enriched with cell type-specific functions. Functional experiments in lineage-selective SEs demonstrated their specific roles in directed differentiation of hMSCs. We also found that some key transcription factors regulated by lineage-selective SEs could form core regulatory circuitry (CRC) to regulate each other's expression and control the hMSCs fate determination. In addition, we found that GWAS SNPs of osteoporosis and obesity were significantly enriched in osteoblasts-selective SEs or adipocytes-selective SEs, respectively. Taken together, our studies unveiled important roles of lineage-selective SEs in hMSCs differentiation into osteoblasts and adipocytes.

ACE2 expression in adipose tissue is associated with cardio-metabolic risk factors and cell type composition-implications for COVID-19

BACKGROUND

COVID-19 severity varies widely. Although some demographic and cardio-metabolic factors, including age and obesity, are associated with increasing risk of severe illness, the underlying mechanism(s) are uncertain.

SUBJECTS/METHODS

In a meta-analysis of three independent studies of 1471 participants in total, we investigated phenotypic and genetic factors associated with subcutaneous adipose tissue expression of Angiotensin I Converting Enzyme 2 (ACE2), measured by RNA-Seq, which acts as a receptor for SARS-CoV-2 cellular entry.

RESULTS

Lower adipose tissue ACE2 expression was associated with multiple adverse cardio-metabolic health indices, including type 2 diabetes (T2D) (P = 9.14 × 10-6), obesity status (P = 4.81 × 10-5), higher serum fasting insulin (P = 5.32 × 10-4), BMI (P = 3.94 × 10-4), and lower serum HDL levels (P = 1.92 × 10-7). ACE2 expression was also associated with estimated proportions of cell types in adipose tissue: lower expression was associated with a lower proportion of microvascular endothelial cells (P = 4.25 × 10-4) and higher proportion of macrophages (P = 2.74 × 10-5). Despite an estimated heritability of 32%, we did not identify any proximal or distal expression quantitative trait loci (eQTLs) associated with adipose tissue ACE2 expression.

CONCLUSIONS

Our results demonstrate that individuals with cardio-metabolic features known to increase risk of severe COVID-19 have lower background ACE2 levels in this highly relevant tissue. Reduced adipose tissue ACE2 expression may contribute to the pathophysiology of cardio-metabolic diseases, as well as the associated increased risk of severe COVID-19.

Large-scale analysis of circulating glutamate and adipose gene expression in relation to abdominal obesity

Circulating levels of the amino acid glutamate are associated with central fat accumulation, yet the pathophysiology of this relationship remains unknown. We aimed to (i) refine and validate the association between circulating glutamate and abdominal obesity in a large twin cohort, and (ii) investigate whether transcriptomic profiles in adipose tissue could provide insight into the biological mechanisms underlying the association. First, in a cohort of 4665 individuals from the TwinsUK resource, we identified individuals with abdominal obesity and compared prevalence of the latter across circulating glutamate quintiles. Second, we used transcriptomic signatures generated from adipose tissue, both subcutaneous and visceral, to investigate associations with circulating glutamate levels. Individuals in the top circulating glutamate quintile had a sevenfold higher prevalence of abdominal obesity compared to those in the bottom quintile. The adipose tissue transcriptomic analyses identified GLUL, encoding Glutamate-Ammonia Ligase, as being associated with circulating glutamate and abdominal obesity, with pronounced signatures in the visceral depot. In conclusion, circulating glutamate is positively associated with the prevalence of abdominal obesity which relates to dysregulated GLUL expression specifically in visceral adipose tissue.

Lifestyle and metabolic factors for nonalcoholic fatty liver disease: Mendelian randomization study

The risk factors for nonalcoholic fatty liver disease (NAFLD) have not been clearly identified. We conducted a Mendelian randomization (MR) study to explore this. Independent genetic variants strongly associated with 5 lifestyle and 9 metabolic factors were selected as instrumental variables from corresponding genome-wide association studies (GWASs). Summary-level data for NAFLD were obtained from a GWAS meta-analysis of 8434 cases and 770,180 non-cases (discovery dataset) and another GWAS meta-analysis of 1483 cases and 17,781 non-cases (replication dataset). Univariable and multivariable MR analyses were performed. There were associations with NAFLD for lifetime smoking index (odds ratio (OR) 1.59, 95% confidence interval (CI) 1.31–1.93 per SD-increase), body mass index (BMI, OR 1.33, 95% CI 1.23–1.43 per SD-increase), waist circumference (OR 1.82; 95% CI 1.48–2.24 per SD-increase), type 2 diabetes (OR 1.21, 95% CI 1.15–1.27 per unit increase in log-transformed odds), systolic blood pressure (OR 1.17; 95% CI 1.07–1.26 per 10 mmHg increase), high-density lipoprotein cholesterol (OR 0.84, 95% CI 0.77–0.90 per SD-increase), and triglycerides (OR 1.23, 95% CI 1.15–1.33 per SD-increase). The associations for type 2 diabetes, systolic blood pressure, triglycerides, but not for high-density lipoprotein cholesterol remained strong after adjusting for genetically-predicted BMI. Genetic liability to type 2 diabetes mediated 51.4% (95% CI 13.4–89.3%) of the BMI-effects on NAFLD risk. There were suggestive inverse associations of genetically-predicted alcohol, coffee, and caffeine consumption, and vigorous physical activity with NAFLD risk. This study identified several lifestyle and metabolic factors that may be causally implicated in NAFLD.

The immune checkpoint B7-H3 (CD276) regulates adipocyte progenitor metabolism and obesity development

The immune checkpoint B7-H3 (CD276) is a member of the B7 family that has been studied in the tumor microenvironment and immunotherapy, but its potential role in metabolism remains largely unknown. Here, we show that B7-H3 is highly expressed in mouse and human adipose tissue at steady state, with the highest levels in adipocyte progenitor cells. B7-H3 is rapidly down-regulated upon the initiation of adipocyte differentiation. Combined RNA sequencing and metabolic studies reveal that B7-H3 stimulates glycolytic and mitochondrial activity of adipocyte progenitors. Loss of B7-H3 in progenitors results in impaired oxidative metabolism program and increased lipid accumulation in derived adipocytes. Consistent with these observations, mice knocked out for B7-H3 develop spontaneous obesity, metabolic dysfunction, and adipose tissue inflammation. Our results reveal an unexpected metabolic role for B7-H3 in adipose tissue and open potential new avenues for the treatment of metabolic diseases by targeting the B7-H3 pathway.

Transcriptome-wide analyses of adipose tissue in outbred rats reveal genetic regulatory mechanisms relevant for human obesity

Transcriptomic analysis in metabolically active tissues allows a systems genetics approach to identify causal genes and networks involved in metabolic disease. Outbred heterogeneous stock (HS) rats are used for genetic mapping of complex traits, but to-date, a systems genetics analysis of metabolic tissues has not been done. We investigated whether adiposity-associated genes and gene co-expression networks in outbred heterogeneous stock (HS) rats overlap those found in humans. We analyzed RNAseq data from adipose tissue of 415 male HS rats, correlated these transcripts with body weight (BW) and compared transcriptome signatures to two human cohorts: the "African American Genetics of Metabolism and Expression" and "Metabolic Syndrome in Men". We used weighted gene co-expression network analysis to identify adiposity-associated gene networks and mediation analysis to identify genes under genetic control whose expression drives adiposity. We identified 554 orthologous "consensus genes" whose expression correlates with BW in the rat and with body mass index (BMI) in both human cohorts. Consensus genes fell within eight co-expressed networks and were enriched for genes involved in immune system function, cell growth, extracellular matrix organization and lipid metabolic processes. We identified 19 consensus genes for which genetic variation may influence BW via their expression, including those involved in lipolysis (e.g., Hcar1), inflammation (e.g., Rgs1), adipogenesis (e.g., Tmem120b) or no previously known role in obesity (e.g., St14, Msa4a6). Strong concordance between HS rat and human BW/BMI associated transcripts demonstrates translational utility of the rat model, while identification of novel genes expands our knowledge of the genetics underlying obesity.

Differential Genetic and Epigenetic Effects of the KLF14 Gene on Body Shape Indices and Metabolic Traits

The KLF14 gene is a key metabolic transcriptional transregulator with monoallelic maternal expression. KLF14 variants are only associated with adipose tissue gene expression, and KLF14 promoter methylation is strongly associated with age. This study investigated whether age, sex, and obesity mediate the effects of KLF14 variants and DNA methylation status on body shape indices and metabolic traits. In total, the data of 78,742 and 1636 participants from the Taiwan Biobank were included in the regional plot association analysis for KLF14 variants and KLF14 methylation, respectively. Regional plot association studies revealed that the KLF14 rs4731702 variant and the nearby strong linkage disequilibrium polymorphisms were the lead variants for lipid profiles, blood pressure status, insulin resistance surrogate markers, and metabolic syndrome mainly in female participants and for body shape indices mainly in obese women. Significant age-dependent associations between KLF14 promoter methylation levels and body shape indices, and metabolic traits were also noted predominantly in female participants. KLF14 variants and KLF14 hypermethylation status were associated with metabolically healthy and unhealthy phenotypes, respectively, in obese individuals, and only the KLF14 variants demonstrated a significant association with both higher adiposity and lower cardiometabolic risk in the same allele, revealing uncoupled excessive adiposity from its cardiometabolic comorbidities, especially in obese women. Variations of KLF14 are associated with body shape indices, metabolic traits, insulin resistance, and metabolically healthy status. Differential genetic and epigenetic effects of KLF14 are age-, sex- and obesity-dependent. These results provided a personalized reference for the management of cardiometabolic diseases in precision medicine.

Adipocyte-Specific Modulation of KLF14 Expression in Mice Leads to Sex-Dependent Impacts on Adiposity and Lipid Metabolism

Genome-wide association studies identified single nucleotide polymorphisms on chromosome 7 upstream of KLF14 to be associated with metabolic syndrome traits and increased risk for type 2 diabetes (T2D). The associations were more significant in women than in men. The risk allele carriers expressed lower levels of the transcription factor KLF14 in adipose tissues than nonrisk allele carriers. To investigate how adipocyte KLF14 regulates metabolic traits in a sex-dependent manner, we characterized high-fat diet-fed male and female mice with adipocyte-specific Klf14 deletion or overexpression. Klf14 deletion resulted in increased fat mass in female mice and decreased fat mass in male mice. Female Klf14-deficient mice had overall smaller adipocytes in subcutaneous fat depots but larger adipocytes in parametrial depots, indicating a shift in lipid storage from subcutaneous to visceral fat depots. They had reduced metabolic rates and increased respiratory exchange ratios consistent with increased use of carbohydrates as an energy source. Fasting- and isoproterenol-induced adipocyte lipolysis was defective in female Klf14-deficient mice, and concomitantly, adipocyte triglycerides lipase mRNA levels were downregulated. Female Klf14-deficient mice cleared blood triglyceride and nonesterified fatty acid less efficiently than wild-type. Finally, adipocyte-specific overexpression of Klf14 resulted in lower total body fat in female but not male mice. Taken together, consistent with human studies, adipocyte KLF14 deficiency in female but not in male mice causes increased adiposity and redistribution of lipid storage from subcutaneous to visceral adipose tissues. Increasing KLF14 abundance in adipocytes of females with obesity and T2D may provide a novel treatment option to alleviate metabolic abnormalities.

Genetically regulated multi-omics study for symptom clusters of posttraumatic stress disorder highlights pleiotropy with hematologic and cardio-metabolic traits

Posttraumatic stress disorder (PTSD) is a psychiatric disorder that may arise in response to severe traumatic event and is diagnosed based on three main symptom clusters (reexperiencing, avoidance, and hyperarousal) per the Diagnostic Manual of Mental Disorders (version DSM-IV-TR). In this study, we characterized the biological heterogeneity of PTSD symptom clusters by performing a multi-omics investigation integrating genetically regulated gene, splicing, and protein expression in dorsolateral prefrontal cortex tissue within a sample of US veterans enrolled in the Million Veteran Program (N total = 186,689). We identified 30 genes in 19 regions across the three PTSD symptom clusters. We found nine genes to have cell-type specific expression, and over-representation of miRNA-families - miR-148, 30, and 8. Gene-drug target prioritization approach highlighted cyclooxygenase and acetylcholine compounds. Next, we tested molecular-profile based phenome-wide impact of identified genes with respect to 1678 phenotypes derived from the Electronic Health Records of the Vanderbilt University biorepository (N = 70,439). Lastly, we tested for local genetic correlation across PTSD symptom clusters which highlighted metabolic (e.g., obesity, diabetes, vascular health) and laboratory traits (e.g., neutrophil, eosinophil, tau protein, creatinine kinase). Overall, this study finds comprehensive genomic evidence including clinical and regulatory profiles between PTSD, hematologic and cardiometabolic traits, that support comorbidities observed in epidemiologic studies of PTSD.

Sex hormones, adiposity, and metabolic traits in men and women: a Mendelian randomisation study

OBJECTIVE

Epidemiological and clinical studies have highlighted important roles for sex hormones in the regulation of fat distribution and systemic metabolism. We investigated the bidirectional associations between bioavailable serum testosterone (BioT) in both sexes and oestradiol (E2) in men and adiposity and metabolic traits using Mendelian randomisation (MR).

DESIGN AND METHODS

As genetic instruments for sex hormones, we selected all the genome-wide significant, independent signals from a genome-wide association studies (GWAS) in up to 425 097 European ancestry UK Biobank participants. European population-specific, summary-level data for adiposity, metabolic, and blood pressure traits were obtained from the largest publicly available GWAS. Sex-specific, two-sample MR analyses were used to estimate the associations of sex hormones with these traits and vice versa.

RESULTS

In women, higher BioT was associated with obesity, upper-body fat distribution, and low HDL-cholesterol although, based on analyses modelling the sex hormone-binding globulin-independent effects of BioT, the last two associations might be indirect. Conversely, obesity and android fat distribution were associated with elevated serum BioT. In men, higher BioT was associated with lower hip circumference and lower fasting glucose. Reciprocally, obesity was associated with lower BioT and higher E2, while upper-body fat distribution and raised triglycerides were associated with lower E2.

CONCLUSIONS

Adipose tissue and metabolic dysfunction are associated with deranged sex hormone levels in both sexes. In women, elevated BioT might be a cause of obesity. Conversely, in men, higher BioT appears to have beneficial effects on adiposity and glucose metabolism.

On Variance of the Treatment Effect in the Treated When Estimated by Inverse Probability Weighting

In the analysis of observational studies, inverse probability weighting (IPW) is commonly used to consistently estimate the average treatment effect (ATE) or the average treatment effect in the treated (ATT). The variance of the IPW ATE estimator is often estimated by assuming that the weights are known and then using the so-called "robust" (Huber-White) sandwich estimator, which results in conservative standard errors (SEs). Here we show that using such an approach when estimating the variance of the IPW ATT estimator does not necessarily result in conservative SE estimates. That is, assuming the weights are known, the robust sandwich estimator may be either conservative or anticonservative. Thus, confidence intervals for the ATT using the robust SE estimate will not be valid, in general. Instead, stacked estimating equations which account for the weight estimation can be used to compute a consistent, closed-form variance estimator for the IPW ATT estimator. The 2 variance estimators are compared via simulation studies and in a data analysis of the association between smoking and gene expression.

β3-Adrenergic receptor downregulation leads to adipocyte catecholamine resistance in obesity

The dysregulation of energy homeostasis in obesity involves multihormone resistance. Although leptin and insulin resistance have been well characterized, catecholamine resistance remains largely unexplored. Murine β3-adrenergic receptor expression in adipocytes is orders of magnitude higher compared with that of other isoforms. While resistant to classical desensitization pathways, its mRNA (Adrb3) and protein expression are dramatically downregulated after ligand exposure (homologous desensitization). β3-Adrenergic receptor downregulation also occurs after high-fat diet feeding, concurrent with catecholamine resistance and elevated inflammation. This downregulation is recapitulated in vitro by TNF-α treatment (heterologous desensitization). Both homologous and heterologous desensitization of Adrb3 were triggered by induction of the pseudokinase TRIB1 downstream of the EPAC/RAP2A/PI-PLC pathway. TRIB1 in turn degraded the primary transcriptional activator of Adrb3, CEBPα. EPAC/RAP inhibition enhanced catecholamine-stimulated lipolysis and energy expenditure in obese mice. Moreover, adipose tissue expression of genes in this pathway correlated with body weight extremes in a cohort of genetically diverse mice and with BMI in 2 independent cohorts of humans. These data implicate a signaling axis that may explain reduced hormone-stimulated lipolysis in obesity and resistance to therapeutic interventions with β3-adrenergic receptor agonists.

A human-based multi-gene signature enables quantitative drug repurposing for metabolic disease

Insulin resistance (IR) contributes to the pathophysiology of diabetes, dementia, viral infection, and cardiovascular disease. Drug repurposing (DR) may identify treatments for IR; however, barriers include uncertainty whether in vitro transcriptomic assays yield quantitative pharmacological data, or how to optimise assay design to best reflect in vivo human disease. We developed a clinical-based human tissue IR signature by combining lifestyle-mediated treatment responses (>500 human adipose and muscle biopsies) with biomarkers of disease status (fasting IR from >1200 biopsies). The assay identified a chemically diverse set of >130 positively acting compounds, highly enriched in true positives, that targeted 73 proteins regulating IR pathways. Our multi-gene RNA assay score reflected the quantitative pharmacological properties of a set of epidermal growth factor receptor-related tyrosine kinase inhibitors, providing insight into drug target specificity; an observation supported by deep learning-based genome-wide predicted pharmacology. Several drugs identified are suitable for evaluation in patients, particularly those with either acute or severe chronic IR.

Subcutaneous adipose tissue splice quantitative trait loci reveal differences in isoform usage associated with cardiometabolic traits

Alternate splicing events can create isoforms that alter gene function, and genetic variants associated with alternate gene isoforms may reveal molecular mechanisms of disease. We used subcutaneous adipose tissue of 426 Finnish men from the METSIM study and identified splice junction quantitative trait loci (sQTLs) for 6,077 splice junctions (FDR < 1%). In the same individuals, we detected expression QTLs (eQTLs) for 59,443 exons and 15,397 genes (FDR < 1%). We identified 595 genes with an sQTL and exon eQTL but no gene eQTL, which could indicate potential isoform differences. Of the significant sQTL signals, 2,114 (39.8%) included at least one proxy variant (linkage disequilibrium r2 > 0.8) located within an intron spanned by the splice junction. We identified 203 sQTLs that colocalized with 141 genome-wide association study (GWAS) signals for cardiometabolic traits, including 25 signals for lipid traits, 24 signals for body mass index (BMI), and 12 signals for waist-hip ratio adjusted for BMI. Among all 141 GWAS signals colocalized with an sQTL, we detected 26 that also colocalized with an exon eQTL for an exon skipped by the sQTL splice junction. At a GWAS signal for high-density lipoprotein cholesterol colocalized with an NR1H3 sQTL splice junction, we show that the alternative splice product encodes an NR1H3 transcription factor that lacks a DNA binding domain and fails to activate transcription. Together, these results detect splicing events and candidate mechanisms that may contribute to gene function at GWAS loci.

A Wars2 mutant mouse shows a sex and diet specific change in fat distribution, reduced food intake and depot-specific upregulation of WAT browning

Background: Increased waist-to-hip ratio (WHR) is associated with increased mortality and risk of type 2 diabetes and cardiovascular disease. The TBX15-WARS2 locus has consistently been associated with increased WHR. Previous study of the hypomorphic Wars2 ^V117L/V117L mouse model found phenotypes including severely reduced fat mass, and white adipose tissue (WAT) browning, suggesting Wars2 could be a potential modulator of fat distribution and WAT browning. Methods: To test for differences in browning induction across different adipose depots of Wars2 ^V117L/V117L mice, we measured multiple browning markers of a 4-month old chow-fed cohort in subcutaneous and visceral WAT and brown adipose tissue (BAT). To explain previously observed fat mass loss, we also tested for the upregulation of plasma mitokines FGF21 and GDF15 and for differences in food intake in the same cohort. Finally, to test for diet-associated differences in fat distribution, we placed Wars2 ^V117L/V117L mice on low-fat or high-fat diet (LFD, HFD) and assessed their body composition by Echo-MRI and compared terminal adipose depot weights at 6 months of age. Results: The chow-fed Wars2 ^V117L/V117L mice showed more changes in WAT browning marker gene expression in the subcutaneous inguinal WAT depot (iWAT) than in the visceral gonadal WAT depot (gWAT). These mice also demonstrated reduced food intake and elevated plasma FGF21 and GDF15, and mRNA from heart and BAT. When exposed to HFD, the Wars2 ^V117L/V117L mice showed resistance to diet-induced obesity and a male and HFD-specific reduction of gWAT: iWAT ratio. Conclusion: Severe reduction of Wars2 gene function causes a systemic phenotype which leads to upregulation of FGF21 and GDF15, resulting in reduced food intake and depot-specific changes in browning and fat mass.

Multi-Omic Approaches to Identify Genetic Factors in Metabolic Syndrome

Metabolic syndrome (MetS) is a highly heritable disease and a major public health burden worldwide. MetS diagnosis criteria are met by the simultaneous presence of any three of the following: high triglycerides, low HDL/high LDL cholesterol, insulin resistance, hypertension, and central obesity. These diseases act synergistically in people suffering from MetS and dramatically increase risk of morbidity and mortality due to stroke and cardiovascular disease, as well as certain cancers. Each of these component features is itself a complex disease, as is MetS. As a genetically complex disease, genetic risk factors for MetS are numerous, but not very powerful individually, often requiring specific environmental stressors for the disease to manifest. When taken together, all sequence variants that contribute to MetS disease risk explain only a fraction of the heritable variance, suggesting additional, novel loci have yet to be discovered. In this article, we will give a brief overview on the genetic concepts needed to interpret genome-wide association studies (GWAS) and quantitative trait locus (QTL) data, summarize the state of the field of MetS physiological genomics, and to introduce tools and resources that can be used by the physiologist to integrate genomics into their own research on MetS and any of its component features. There is a wealth of phenotypic and molecular data in animal models and humans that can be leveraged as outlined in this article. Integrating these multi-omic QTL data for complex diseases such as MetS provides a means to unravel the pathways and mechanisms leading to complex disease and promise for novel treatments. © 2022 American Physiological Society. Compr Physiol 12:1-40, 2022.

Inferred expression regulator activities suggest genes mediating cardiometabolic genetic signals

Expression QTL (eQTL) analyses have suggested many genes mediating genome-wide association study (GWAS) signals but most GWAS signals still lack compelling explanatory genes. We have leveraged an adipose-specific gene regulatory network to infer expression regulator activities and phenotypic master regulators (MRs), which were used to detect activity QTLs (aQTLs) at cardiometabolic trait GWAS loci. Regulator activities were inferred with the VIPER algorithm that integrates enrichment of expected expression changes among a regulator's target genes with confidence in their regulator-target network interactions and target overlap between different regulators (i.e., pleiotropy). Phenotypic MRs were identified as those regulators whose activities were most important in predicting their respective phenotypes using random forest modeling. While eQTLs were typically more significant than aQTLs in cis, the opposite was true among candidate MRs in trans. Several GWAS loci colocalized with MR trans-eQTLs/aQTLs in the absence of colocalized cis-QTLs. Intriguingly, at the 1p36.1 BMI GWAS locus the EPHB2 cis-aQTL was stronger than its cis-eQTL and colocalized with the GWAS signal and 35 BMI MR trans-aQTLs, suggesting the GWAS signal may be mediated by effects on EPHB2 activity and its downstream effects on a network of BMI MRs. These MR and aQTL analyses represent systems genetic methods that may be broadly applied to supplement standard eQTL analyses for suggesting molecular effects mediating GWAS signals.

Therapeutic radiation exposure of the abdomen during childhood induces chronic adipose tissue dysfunction

BACKGROUNDChildhood cancer survivors who received abdominal radiotherapy (RT) or total body irradiation (TBI) are at increased risk for cardiometabolic disease, but the underlying mechanisms are unknown. We hypothesize that RT-induced adipose tissue dysfunction contributes to the development of cardiometabolic disease in the expanding population of childhood cancer survivors.METHODSWe performed clinical metabolic profiling of adult childhood cancer survivors previously exposed to TBI, abdominal RT, or chemotherapy alone, alongside a group of healthy controls. Study participants underwent abdominal s.c. adipose biopsies to obtain tissue for bulk RNA sequencing. Transcriptional signatures were analyzed using pathway and network analyses and cellular deconvolution.RESULTSIrradiated adipose tissue is characterized by a gene expression signature indicative of a complex macrophage expansion. This signature includes activation of the TREM2-TYROBP network, a pathway described in diseases of chronic tissue injury. Radiation exposure of adipose is further associated with dysregulated adipokine secretion, specifically a decrease in insulin-sensitizing adiponectin and an increase in insulin resistance-promoting plasminogen activator inhibitor-1. Accordingly, survivors exhibiting these changes have early signs of clinical metabolic derangement, such as increased fasting glucose and hemoglobin A1c.CONCLUSIONChildhood cancer survivors exposed to abdominal RT or TBI during treatment exhibit signs of chronic s.c. adipose tissue dysfunction, manifested as dysregulated adipokine secretion that may negatively impact their systemic metabolic health.FUNDINGThis study was supported by Rockefeller University Hospital; National Institute of General Medical Sciences (T32GM007739); National Center for Advancing Translational Sciences (UL1 TR001866); National Cancer Institute (P30CA008748); American Cancer Society (133831-CSDG-19-117-01-CPHPS); American Diabetes Association (1-17-ACE-17); and an anonymous donor (MSKCC).

Glycogen metabolism links glucose homeostasis to thermogenesis in adipocytes

Adipocytes increase energy expenditure in response to prolonged sympathetic activation via persistent expression of uncoupling protein 1 (UCP1)^1,2. Here we report that the regulation of glycogen metabolism by catecholamines is critical for UCP1 expression. Chronic β-adrenergic activation leads to increased glycogen accumulation in adipocytes expressing UCP1. Adipocyte-specific deletion of a scaffolding protein, protein targeting to glycogen (PTG), reduces glycogen levels in beige adipocytes, attenuating UCP1 expression and responsiveness to cold or β-adrenergic receptor-stimulated weight loss in obese mice. Unexpectedly, we observed that glycogen synthesis and degradation are increased in response to catecholamines, and that glycogen turnover is required to produce reactive oxygen species leading to the activation of p38 MAPK, which drives UCP1 expression. Thus, glycogen has a key regulatory role in adipocytes, linking glucose metabolism to thermogenesis.

Chromatin accessibility and gene expression during adipocyte differentiation identify context-dependent effects at cardiometabolic GWAS loci

Chromatin accessibility and gene expression in relevant cell contexts can guide identification of regulatory elements and mechanisms at genome-wide association study (GWAS) loci. To identify regulatory elements that display differential activity across adipocyte differentiation, we performed ATAC-seq and RNA-seq in a human cell model of preadipocytes and adipocytes at days 4 and 14 of differentiation. For comparison, we created a consensus map of ATAC-seq peaks in 11 human subcutaneous adipose tissue samples. We identified 58,387 context-dependent chromatin accessibility peaks and 3,090 context-dependent genes between all timepoint comparisons (log2 fold change>1, FDR<5%) with 15,919 adipocyte- and 18,244 preadipocyte-dependent peaks. Adipocyte-dependent peaks showed increased overlap (60.1%) with Roadmap Epigenomics adipocyte nuclei enhancers compared to preadipocyte-dependent peaks (11.5%). We linked context-dependent peaks to genes based on adipocyte promoter capture Hi-C data, overlap with adipose eQTL variants, and context-dependent gene expression. Of 16,167 context-dependent peaks linked to a gene, 5,145 were linked by two or more strategies to 1,670 genes. Among GWAS loci for cardiometabolic traits, adipocyte-dependent peaks, but not preadipocyte-dependent peaks, showed significant enrichment (LD score regression P<0.005) for waist-to-hip ratio and modest enrichment (P < 0.05) for HDL-cholesterol. We identified 659 peaks linked to 503 genes by two or more approaches and overlapping a GWAS signal, suggesting a regulatory mechanism at these loci. To identify variants that may alter chromatin accessibility between timepoints, we identified 582 variants in 454 context-dependent peaks that demonstrated allelic imbalance in accessibility (FDR<5%), of which 55 peaks also overlapped GWAS variants. At one GWAS locus for palmitoleic acid, rs603424 was located in an adipocyte-dependent peak linked to SCD and exhibited allelic differences in transcriptional activity in adipocytes (P = 0.003) but not preadipocytes (P = 0.09). These results demonstrate that context-dependent peaks and genes can guide discovery of regulatory variants at GWAS loci and aid identification of regulatory mechanisms.

Cardiovascular and metabolic diseases are widespread, and an increased understanding of genetic mechanisms behind these diseases could improve treatment. Chromatin accessibility and gene expression in relevant cell contexts can guide identification of regulatory elements and genetic mechanisms for disease traits. A relevant context for cardiovascular and metabolic disease traits is adipocyte differentiation. To identify regulatory elements and genes that display differences in activity during adipocyte differentiation, we profiled chromatin accessibility and gene expression in a human cell model of preadipocytes and adipocytes. We identified chromatin regions that change accessibility during differentiation and predicted genes they may affect. We also linked these chromatin regions to genetic variants associated with risk of disease. At one genomic region linked to fatty acids, a chromatin region more accessible in adipocytes linked to a fatty acid synthesis gene and exhibited allelic differences in transcriptional activity in adipocytes but not preadipocytes. These results demonstrate that chromatin regions and genes that change during cell context can guide discovery of regulatory variants and aid identification of disease mechanisms.

Genes in human obesity loci are causal obesity genes in C. elegans

Obesity and its associated metabolic syndrome are a leading cause of morbidity and mortality. Given the disease's heavy burden on patients and the healthcare system, there has been increased interest in identifying pharmacological targets for the treatment and prevention of obesity. Towards this end, genome-wide association studies (GWAS) have identified hundreds of human genetic variants associated with obesity. The next challenge is to experimentally define which of these variants are causally linked to obesity, and could therefore become targets for the treatment or prevention of obesity. Here we employ high-throughput in vivo RNAi screening to test for causality 293 C. elegans orthologs of human obesity-candidate genes reported in GWAS. We RNAi screened these 293 genes in C. elegans subject to two different feeding regimens: (1) regular diet, and (2) high-fructose diet, which we developed and present here as an invertebrate model of diet-induced obesity (DIO). We report 14 genes that promote obesity and 3 genes that prevent DIO when silenced in C. elegans. Further, we show that knock-down of the 3 DIO genes not only prevents excessive fat accumulation in primary and ectopic fat depots but also improves the health and extends the lifespan of C. elegans overconsuming fructose. Importantly, the direction of the association between expression variants in these loci and obesity in mice and humans matches the phenotypic outcome of the loss-of-function of the C. elegans ortholog genes, supporting the notion that some of these genes would be causally linked to obesity across phylogeny. Therefore, in addition to defining causality for several genes so far merely correlated with obesity, this study demonstrates the value of model systems compatible with in vivo high-throughput genetic screening to causally link GWAS gene candidates to human diseases.

Brain-trait-associated variants impact cell-type-specific gene regulation during neurogenesis

Interpretation of the function of non-coding risk loci for neuropsychiatric disorders and brain-relevant traits via gene expression and alternative splicing quantitative trait locus (e/sQTL) analyses is generally performed in bulk post-mortem adult tissue. However, genetic risk loci are enriched in regulatory elements active during neocortical differentiation, and regulatory effects of risk variants may be masked by heterogeneity in bulk tissue. Here, we map e/sQTLs, and allele-specific expression in cultured cells representing two major developmental stages, primary human neural progenitors (n = 85) and their sorted neuronal progeny (n = 74), identifying numerous loci not detected in either bulk developing cortical wall or adult cortex. Using colocalization and genetic imputation via transcriptome-wide association, we uncover cell-type-specific regulatory mechanisms underlying risk for brain-relevant traits that are active during neocortical differentiation. Specifically, we identified a progenitor-specific eQTL for CENPW co-localized with common variant associations for cortical surface area and educational attainment.

The Link between Type 2 Diabetes Mellitus and the Polymorphisms of Glutathione-Metabolizing Genes Suggests a New Hypothesis Explaining Disease Initiation and Progression

The present study investigated whether type 2 diabetes (T2D) is associated with polymorphisms of genes encoding glutathione-metabolizing enzymes such as glutathione synthetase (GSS) and gamma-glutamyl transferase 7 (GGT7). A total of 3198 unrelated Russian subjects including 1572 T2D patients and 1626 healthy subjects were enrolled. Single nucleotide polymorphisms (SNPs) of the GSS and GGT7 genes were genotyped using the MassArray-4 system. We found that the GSS and GGT7 gene polymorphisms alone and in combinations are associated with T2D risk regardless of sex, age, and body mass index, as well as correlated with plasma glutathione, hydrogen peroxide, and fasting blood glucose levels. Polymorphisms of GSS (rs13041792) and GGT7 (rs6119534 and rs11546155) genes were associated with the tissue-specific expression of genes involved in unfolded protein response and the regulation of proteostasis. Transcriptome-wide association analysis has shown that the pancreatic expression of some of these genes such as EDEM2, MYH7B, MAP1LC3A, and CPNE1 is linked to the genetic risk of T2D. A comprehensive analysis of the data allowed proposing a new hypothesis for the etiology of type 2 diabetes that endogenous glutathione deficiency might be a key condition responsible for the impaired folding of proinsulin which triggered an unfolded protein response, ultimately leading to beta-cell apoptosis and disease development.

Four-and-a-half LIM domain protein 2 (FHL2) deficiency protects mice from diet-induced obesity and high FHL2 expression marks human obesity

OBJECTIVE

Four-and-a-Half-LIM-domain-protein 2 (FHL2) modulates multiple signal transduction pathways but has not been implicated in obesity or energy metabolism. In humans, methylation and expression of the FHL2 gene increases with age, and high FHL2 expression is associated with increased body weight in humans and mice. This led us to hypothesize that FHL2 is a determinant of diet-induced obesity.

METHODS

FHL2-deficient (FHL2-/-) and wild type male mice were fed a high-fat diet. Metabolic phenotyping of these mice, as well as transcriptional analysis of key metabolic tissues was performed. Correlation of the expression of FHL2 and relevant genes was assessed in datasets from white adipose tissue of individuals with and without obesity.

RESULTS

FHL2 Deficiency protects mice from high-fat diet-induced weight gain, whereas glucose handling is normal. We observed enhanced energy expenditure, which may be explained by a combination of changes in multiple tissues; mild activation of brown adipose tissue with increased fatty acid uptake, increased cardiac glucose uptake and browning of white adipose tissue. Corroborating our findings in mice, expression of FHL2 in human white adipose tissue positively correlates with obesity and negatively with expression of browning-associated genes.

CONCLUSION

Our results position FHL2 as a novel regulator of obesity and energy expenditure in mice and human. Given that FHL2 expression increases during aging, we now show that low FHL2 expression associates with a healthy metabolic state.

Genetic architecture modulates diet-induced hepatic mRNA and miRNA expression profiles in Diversity Outbred mice

Genetic approaches in model organisms have consistently demonstrated that molecular traits such as gene expression are under genetic regulation, similar to clinical traits. The resulting expression quantitative trait loci (eQTL) have revolutionized our understanding of genetic regulation and identified numerous candidate genes for clinically relevant traits. More recently, these analyses have been extended to other molecular traits such as protein abundance, metabolite levels, and miRNA expression. Here, we performed global hepatic eQTL and microRNA expression quantitative trait loci (mirQTL) analysis in a population of Diversity Outbred mice fed two different diets. We identified several key features of eQTL and mirQTL, namely differences in the mode of genetic regulation (cis or trans) between mRNA and miRNA. Approximately 50% of mirQTL are regulated by a trans-acting factor, compared to ∼25% of eQTL. We note differences in the heritability of mRNA and miRNA expression and variance explained by each eQTL or mirQTL. In general, cis-acting variants affecting mRNA or miRNA expression explain more phenotypic variance than trans-acting variants. Finally, we investigated the effect of diet on the genetic architecture of eQTL and mirQTL, highlighting the critical effects of environment on both eQTL and mirQTL. Overall, these data underscore the complex genetic regulation of two well-characterized RNA classes (mRNA and miRNA) that have critical roles in the regulation of clinical traits and disease susceptibility