Hepatokine ITIH3 protects against hepatic steatosis by downregulating mitochondrial bioenergetics and de novo lipogenesis.

Recent studies demonstrate that liver secretory proteins, also known as hepatokines, regulate normal development, obesity, and simple steatosis to non-alcoholic steatohepatitis (NASH) progression. Using a panel of ∼100 diverse inbred strains of mice and a cohort of bariatric surgery patients, we found that one such hepatokine, inter-trypsin inhibitor heavy chain 3 (ITIH3), was progressively lower in severe non-alcoholic fatty liver disease (NAFLD) disease states highlighting an inverse relationship between Itih3/ITIH3 expression and NAFLD severity. Follow-up animal and cell culture models demonstrated that hepatic ITIH3 overexpression lowered liver triglyceride and lipid droplet accumulation, respectively. Conversely, ITIH3 knockdown in mice increased the liver triglyceride in two independent NAFLD models. Mechanistically, ITIH3 reduced mitochondrial respiration and this, in turn, reduced liver triglycerides, via downregulated de novo lipogenesis. This was accompanied by increased STAT1 signaling and Stat3 expression, both of which are known to protect against NAFLD/NASH. Our findings indicate hepatokine ITIH3 as a potential biomarker and/or treatment for NAFLD.

Partitioned polygenic risk scores identify distinct types of metabolic dysfunction-associated steatotic liver disease.

Metabolic dysfunction-associated steatotic liver disease (MASLD) encompasses an excess of triglycerides in the liver, which can lead to cirrhosis and liver cancer. While there is solid epidemiological evidence of MASLD coexisting with cardiometabolic disease, several leading genetic risk factors for MASLD do not increase the risk of cardiovascular disease, suggesting no causal relationship between MASLD and cardiometabolic derangement. In this work, we leveraged measurements of visceral adiposity and identified 27 novel genetic loci associated with MASLD. Among these loci, we replicated 6 in several independent cohorts. Next, we generated two partitioned polygenic risk scores (PRS) based on the mechanism of genetic association with MASLD encompassing intra-hepatic lipoprotein retention. The two PRS suggest the presence of at least two distinct types of MASLD, one confined to the liver resulting in a more aggressive liver disease and one that is systemic and results in a higher risk of cardiometabolic disease.

Age-dependent genes in adipose stem and precursor cells affect regulation of fat cell differentiation and link aging to obesity via cellular and genetic interactions.

BACKGROUND: Age and obesity are dominant risk factors for several common cardiometabolic disorders, and both are known to impair adipose tissue function. However, the underlying cellular and genetic factors linking aging and obesity on adipose tissue function have remained elusive. Adipose stem and precursor cells (ASPCs) are an understudied, yet crucial adipose cell type due to their deterministic adipocyte differentiation potential, which impacts the capacity to store fat in a metabolically healthy manner. METHODS: We integrated subcutaneous adipose tissue (SAT) bulk (n=435) and large single-nucleus RNA sequencing (n=105) data with the UK Biobank (UKB) (n=391,701) data to study age-obesity interactions originating from ASPCs by performing cell-type decomposition, differential expression testing, cell-cell communication analyses, and construction of polygenic risk scores for body mass index (BMI). RESULTS: We found that the SAT ASPC proportions significantly decrease with age in an obesity-dependent way consistently in two independent cohorts, both showing that the age dependency of ASPC proportions is abolished by obesity. We further identified 76 genes (72 SAT ASPC marker genes and 4 transcription factors regulating ASPC marker genes) that are differentially expressed by age in SAT and functionally enriched for developmental processes and adipocyte differentiation (i.e., adipogenesis). The 76 age-perturbed ASPC genes include multiple negative regulators of adipogenesis, such as RORA, SMAD3, TWIST2, and ZNF521, form tight clusters of longitudinally co-expressed genes during human adipogenesis, and show age-based differences in cellular interactions between ASPCs and adipose cell types. Finally, our genetic data demonstrate that cis-regional variants of these genes interact with age as predictors of BMI in an obesity-dependent way in the large UKB, while no such gene-age interaction on BMI is observed with non-age-dependent ASPC marker genes, thus independently confirming our cellular ASPC results at the biobank level. CONCLUSIONS: Overall, we discover that obesity prematurely induces a decrease in ASPC proportions and identify 76 developmentally important ASPC genes that implicate altered negative regulation of fat cell differentiation as a mechanism for aging and directly link aging to obesity via significant cellular and genetic interactions.

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.

Deep learning-based phenotype imputation on population-scale biobank data increases genetic discoveries.

Biobanks that collect deep phenotypic and genomic data across many individuals have emerged as a key resource in human genetics. However, phenotypes in biobanks are often missing across many individuals, limiting their utility. We propose AutoComplete, a deep learning-based imputation method to impute or 'fill-in' missing phenotypes in population-scale biobank datasets. When applied to collections of phenotypes measured across ~300,000 individuals from the UK Biobank, AutoComplete substantially improved imputation accuracy over existing methods. On three traits with notable amounts of missingness, we show that AutoComplete yields imputed phenotypes that are genetically similar to the originally observed phenotypes while increasing the effective sample size by about twofold on average. Further, genome-wide association analyses on the resulting imputed phenotypes led to a substantial increase in the number of associated loci. Our results demonstrate the utility of deep learning-based phenotype imputation to increase power for genetic discoveries in existing biobank datasets.

Applying an evolutionary mismatch framework to understand disease susceptibility.

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.

Role of Matrix Gla Protein in Transforming Growth Factor-ß Signaling and Nonalcoholic Steatohepatitis in Mice.

BACKGROUND & AIMS: Nonalcoholic steatohepatitis (NASH) is a complex disease involving both genetic and environmental factors in its onset and progression. We analyzed NASH phenotypes in a genetically diverse cohort of mice, the Hybrid Mouse Diversity Panel, to identify genes contributing to disease susceptibility. METHODS: A "systems genetics" approach, involving integration of genetic, transcriptomic, and phenotypic data, was used to identify candidate genes and pathways in a mouse model of NASH. The causal role of Matrix Gla Protein (MGP) was validated using heterozygous MGP knockout (Mgp+/-) mice. The mechanistic role of MGP in transforming growth factor-beta (TGF-ß) signaling was examined in the LX-2 stellate cell line by using a loss of function approach. RESULTS: Local cis-acting regulation of MGP was correlated with fibrosis, suggesting a causal role in NASH, and this was validated using loss of function experiments in 2 models of diet-induced NASH. Using single-cell RNA sequencing, Mgp was found to be primarily expressed in hepatic stellate cells and dendritic cells in mice. Knockdown of MGP expression in stellate LX-2 cells led to a blunted response to TGF-ß stimulation. This was associated with reduced regulatory SMAD phosphorylation and TGF-ß receptor ALK1 expression as well as increased expression of inhibitory SMAD6. Hepatic MGP expression was found to be significantly correlated with the severity of fibrosis in livers of patients with NASH, suggesting relevance to human disease. CONCLUSIONS: MGP regulates liver fibrosis and TGF-ß signaling in hepatic stellate cells and contributes to NASH pathogenesis.

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.

Cross-Tissue Single-Nucleus RNA Sequencing Discovers Tissue-Resident Adipocytes Involved in Propanoate Metabolism in the Human Heart.

BACKGROUND: Adipocytes are crucial regulators of cardiovascular health. However, not much is known about gene expression profiles of adipocytes residing in nonfat cardiovascular tissues, their genetic regulation, and contribution to coronary artery disease. Here, we investigated whether and how the gene expression profiles of adipocytes in the subcutaneous adipose tissue differ from adipocytes residing in the heart. METHODS: We used single-nucleus RNA-sequencing data sets of subcutaneous adipose tissue and heart and performed in-depth analysis of tissue-resident adipocytes and their cell-cell interactions. RESULTS: We first discovered tissue-specific features of tissue-resident adipocytes, identified functional pathways involved in their tissue specificity, and found genes with cell type-specific expression enrichment in tissue-resident adipocytes. By following up these results, we discovered the propanoate metabolism pathway as a novel distinct characteristic of the heart-resident adipocytes and found a significant enrichment of coronary artery disease genome-wide association study risk variants among the right atrium-specific adipocyte marker genes. Our cell-cell communication analysis identified 22 specific heart adipocyte-associated ligand-receptor pairs and signaling pathways, including THBS (thrombospondin) and EPHA (ephrin type-A), further supporting the distinct tissue-resident role of heart adipocytes. Our results also suggest chamber-level coordination of heart adipocyte expression profiles as we observed a consistently larger number of adipocyte-associated ligand-receptor interactions and functional pathways in the atriums than ventricles. CONCLUSIONS: Overall, we introduce a new function and genetic link to coronary artery disease for the previously unexplored heart-resident adipocytes.

Increased body mass index is linked to systemic inflammation through altered chromatin co-accessibility in human preadipocytes.

Obesity-induced adipose tissue dysfunction can cause low-grade inflammation and downstream obesity comorbidities. Although preadipocytes may contribute to this pro-inflammatory environment, the underlying mechanisms are unclear. We used human primary preadipocytes from body mass index (BMI) -discordant monozygotic (MZ) twin pairs to generate epigenetic (ATAC-sequence) and transcriptomic (RNA-sequence) data for testing whether increased BMI alters the subnuclear compartmentalization of open chromatin in the twins' preadipocytes, causing downstream inflammation. Here we show that the co-accessibility of open chromatin, i.e. compartmentalization of chromatin activity, is altered in the higher vs lower BMI MZ siblings for a large subset ( ~ 88.5 Mb) of the active subnuclear compartments. Using the UK Biobank we show that variants within these regions contribute to systemic inflammation through interactions with BMI on C-reactive protein. In summary, open chromatin co-accessibility in human preadipocytes is disrupted among the higher BMI siblings, suggesting a mechanism how obesity may lead to inflammation via gene-environment interactions.

Cross-tissue omics analysis discovers ten adipose genes encoding secreted proteins in obesity-related non-alcoholic fatty liver disease.

BACKGROUND: Non-alcoholic fatty liver disease (NAFLD) is a fast-growing, underdiagnosed, epidemic. We hypothesise that obesity-related inflammation compromises adipose tissue functions, preventing efficient fat storage, and thus driving ectopic fat accumulation into the liver. METHODS: To identify adipose-based mechanisms and potential serum biomarker candidates (SBCs) for NAFLD, we utilise dual-tissue RNA-sequencing (RNA-seq) data in adipose tissue and liver, paired with histology-based NAFLD diagnosis, from the same individuals in a cohort of obese individuals. We first scan for genes that are differentially expressed (DE) for NAFLD in obese individuals' subcutaneous adipose tissue but not in their liver; encode proteins secreted to serum; and show preferential adipose expression. Then the identified genes are filtered to key adipose-origin NAFLD genes by best subset analysis, knockdown experiments during human preadipocyte differentiation, recombinant protein treatment experiments in human liver HepG2 cells, and genetic analysis. FINDINGS: We discover a set of genes, including 10 SBCs, that may modulate NAFLD pathogenesis by impacting adipose tissue function. Based on best subset analysis, we further follow-up on two SBCs CCDC80 and SOD3 by knockdown in human preadipocytes and subsequent differentiation experiments, which show that they modulate crucial adipogenesis genes, LPL, SREBPF1, and LEP. We also show that treatment of the liver HepG2 cells with the CCDC80 and SOD3 recombinant proteins impacts genes related to steatosis and lipid processing, including PPARA, NFE2L2, and RNF128. Finally, utilizing the adipose NAFLD DE gene cis-regulatory variants associated with serum triglycerides (TGs) in extensive genome-wide association studies (GWASs), we demonstrate a unidirectional effect of serum TGs on NAFLD with Mendelian Randomization (MR) analysis. We also demonstrate that a single SNP regulating one of the SBC genes, rs2845885, produces a significant MR result by itself. This supports the conclusion that genetically regulated adipose expression of the NAFLD DE genes may contribute to NAFLD through changes in serum TG levels. INTERPRETATION: Our results from the dual-tissue transcriptomics screening improve the understanding of obesity-related NAFLD by providing a targeted set of 10 adipose tissue-active genes as new serum biomarker candidates for the currently grossly underdiagnosed fatty liver disease. FUNDING: The work was supported by NIH grants R01HG010505 and R01DK132775. The Genotype-Tissue Expression (GTEx) Project was supported by the Common Fund of the Office of the Director of the National Institutes of Health, and by NCI, NHGRI, NHLBI, NIDA, NIMH, and NINDS. The KOBS study (J. P.) was supported by the Finnish Diabetes Research Foundation, Kuopio University Hospital Project grant (EVO/VTR grants 2005-2019), and the Academy of Finland grant (Contract no. 138006). This study was funded by the European Research Council under the European Union's Horizon 2020 research and innovation program (Grant No. 802825 to M. U. K.). K. H. P. was funded by the Academy of Finland (grant numbers 272376, 266286, 314383, and 335443), the Finnish Medical Foundation, Gyllenberg Foundation, Novo Nordisk Foundation (grant numbers NNF10OC1013354, NNF17OC0027232, and NNF20OC0060547), Finnish Diabetes Research Foundation, Finnish Foundation for Cardiovascular Research, University of Helsinki, and Helsinki University Hospital and Government Research Funds. I. S. was funded by the Instrumentarium Science Foundation. Personal grants to U. T. A. were received from the Matti and Vappu Maukonen Foundation, Ella och Georg Ehrnrooths Stiftelse and the Finnish Foundation for Cardiovascular Research.

Isolation of Nuclei from Human Snap-frozen Liver Tissue for Single-nucleus RNA Sequencing.

Single-nucleus RNA sequencing (snRNA-seq) provides a powerful tool for studying cell type composition in heterogenous tissues. The liver is a vital organ composed of a diverse set of cell types; thus, single-cell technologies could greatly facilitate the deconvolution of liver tissue composition and various downstream omics analyses at the cell-type level. Applying single-cell technologies to fresh liver biopsies can, however, be very challenging, and snRNA-seq of snap-frozen liver biopsies requires some optimization given the high nucleic acid content of the solid liver tissue. Therefore, an optimized protocol for snRNA-seq specifically targeted for the use of frozen liver samples is needed to improve our understanding of human liver gene expression at the cell-type resolution. We present a protocol for performing nuclei isolation from snap-frozen liver tissues, as well as guidance on the application of snRNA-seq. We also provide guidance on optimizing the protocol to different tissue and sample types.

Impaired Adipocyte SLC7A10 Promotes Lipid Storage in Association With Insulin Resistance and Altered BCAA Metabolism.

CONTEXT: The neutral amino acid transporter SLC7A10/ASC-1 is an adipocyte-expressed gene with reduced expression in insulin resistance and obesity. Inhibition of SLC7A10 in adipocytes was shown to increase lipid accumulation despite decreasing insulin-stimulated uptake of glucose, a key substrate for de novo lipogenesis. These data imply that alternative lipogenic substrates to glucose fuel continued lipid accumulation during insulin resistance in obesity. OBJECTIVE: We examined whether increased lipid accumulation during insulin resistance in adipocytes may involve alter flux of lipogenic amino acids dependent on SLC7A10 expression and activity, and whether this is reflected by extracellular and circulating concentrations of marker metabolites. METHODS: In adipocyte cultures with impaired SLC7A10, we performed RNA sequencing and relevant functional assays. By targeted metabolite analyses (GC-MS/MS), flux of all amino acids and selected metabolites were measured in human and mouse adipose cultures. Additionally, SLC7A10 mRNA levels in human subcutaneous adipose tissue (SAT) were correlated to candidate metabolites and adiposity phenotypes in 2 independent cohorts. RESULTS: SLC7A10 impairment altered expression of genes related to metabolic processes, including branched-chain amino acid (BCAA) catabolism, lipogenesis, and glyceroneogenesis. In 3T3-L1 adipocytes, SLC7A10 inhibition increased fatty acid uptake and cellular content of glycerol and cholesterol. SLC7A10 impairment in SAT cultures altered uptake of aspartate and glutamate, and increased net uptake of BCAAs, while increasing the net release of the valine catabolite 3- hydroxyisobutyrate (3-HIB). In human cohorts, SLC7A10 mRNA correlated inversely with total fat mass, circulating triacylglycerols, BCAAs, and 3-HIB. CONCLUSION: Reduced SLC7A10 activity strongly affects flux of BCAAs in adipocytes, which may fuel continued lipogenesis during insulin resistance, and be reflected in increased circulating levels of the valine-derived catabolite 3-HIB.

Evolutionary mismatch and the role of GxE interactions in human disease.

Globally, we are witnessing the rise of complex, non-communicable diseases (NCDs) related to changes in our daily environments. Obesity, asthma, cardiovascular disease, and type 2 diabetes are part of a long list of "lifestyle" diseases that were rare throughout human history but are now common. A key idea from anthropology and evolutionary biology-the evolutionary mismatch hypothesis-seeks to explain this phenomenon. It posits that humans evolved in environments that radically differ from the ones experienced by most people today, and thus traits that were advantageous in past environments may now be "mismatched" and disease-causing. This hypothesis is, at its core, a genetic one: it predicts that loci with a history of selection will exhibit "genotype by environment" (GxE) interactions and have differential health effects in ancestral versus modern environments. Here, we discuss how this concept could be leveraged to uncover the genetic architecture of NCDs in a principled way. Specifically, we advocate for partnering with small-scale, subsistence-level groups that are currently transitioning from environments that are arguably more "matched" with their recent evolutionary history to those that are more "mismatched". These populations provide diverse genetic backgrounds as well as the needed levels and types of environmental variation necessary for mapping GxE interactions in an explicit mismatch framework. Such work would make important contributions to our understanding of environmental and genetic risk factors for NCDs across diverse ancestries and sociocultural contexts.

Single-cell dissection of the obesity-exercise axis in adipose-muscle tissues implies a critical role for mesenchymal stem cells.

Exercise training is critical for the prevention and treatment of obesity, but its underlying mechanisms remain incompletely understood given the challenge of profiling heterogeneous effects across multiple tissues and cell types. Here, we address this challenge and opposing effects of exercise and high-fat diet (HFD)-induced obesity at single-cell resolution in subcutaneous and visceral white adipose tissue and skeletal muscle in mice with diet and exercise training interventions. We identify a prominent role of mesenchymal stem cells (MSCs) in obesity and exercise-induced tissue adaptation. Among the pathways regulated by exercise and HFD in MSCs across the three tissues, extracellular matrix remodeling and circadian rhythm are the most prominent. Inferred cell-cell interactions implicate within- and multi-tissue crosstalk centered around MSCs. Overall, our work reveals the intricacies and diversity of multi-tissue molecular responses to exercise and obesity and uncovers a previously underappreciated role of MSCs in tissue-specific and multi-tissue beneficial effects of exercise.

Hyperinsulinemia Is Highly Associated With Markers of Hepatocytic Senescence in Two Independent Cohorts.

Cellular senescence is an essentially irreversible growth arrest that occurs in response to various cellular stressors and may contribute to development of type 2 diabetes mellitus and nonalcoholic fatty liver disease (NAFLD). In this article, we investigated whether chronically elevated insulin levels are associated with cellular senescence in the human liver. In 107 individuals undergoing bariatric surgery, hepatic senescence markers were assessed by immunohistochemistry as well as transcriptomics. A subset of 180 participants from the ongoing Finnish Kuopio OBesity Surgery (KOBS) study was used as validation cohort. We found plasma insulin to be highly associated with various markers of cellular senescence in liver tissue. The liver transcriptome of individuals with high insulin revealed significant upregulation of several genes associated with senescence: p21, TGFß, PI3K, HLA-G, IL8, p38, Ras, and E2F. Insulin associated with hepatic senescence independently of NAFLD and plasma glucose. By using transcriptomic data from the KOBS study, we could validate the association of insulin with p21 in the liver. Our results support a potential role for hyperinsulinemia in induction of cellular senescence in the liver. These findings suggest possible benefits of lowering insulin levels in obese individuals with insulin resistance.

Long-range chromosomal interactions increase and mark repressed gene expression during adipogenesis.

Obesity perturbs central functions of human adipose tissue, centred on differentiation of preadipocytes to adipocytes, i.e., adipogenesis. The large environmental component of obesity makes it important to elucidate epigenetic regulatory factors impacting adipogenesis. Promoter Capture Hi-C (pCHi-C) has been used to identify chromosomal interactions between promoters and associated regulatory elements. However, long range interactions (LRIs) greater than 1 Mb are often filtered out of pCHi-C datasets, due to technical challenges and their low prevalence. To elucidate the unknown role of LRIs in adipogenesis, we investigated preadipocyte differentiation to adipocytes using pCHi-C and bulk and single nucleus RNA-seq data. We first show that LRIs are reproducible between biological replicates, and they increase >2-fold in frequency across adipogenesis. We further demonstrate that genomic loci containing LRIs are more epigenetically repressed than regions without LRIs, corresponding to lower gene expression in the LRI regions. Accordingly, as preadipocytes differentiate into adipocytes, LRI regions are more likely to contain repressed preadipocyte marker genes; whereas these same LRI regions are depleted of actively expressed adipocyte marker genes. Finally, we show that LRIs can be used to restrict multiple testing of the long-range cis-eQTL analysis to identify variants that regulate genes via LRIs. We exemplify this by identifying a putative long range cis regulatory mechanism at the LYPLAL1/TGFB2 obesity locus. In summary, we identify LRIs that mark repressed regions of the genome, and these interactions increase across adipogenesis, pinpointing developmental regions that need to be repressed in a cell-type specific way for adipogenesis to proceed.

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.

Human liver single nucleus and single cell RNA sequencing identify a hepatocellular carcinoma-associated cell-type affecting survival.

BACKGROUND: Hepatocellular carcinoma (HCC) is a common primary liver cancer with poor overall survival. We hypothesized that there are HCC-associated cell-types that impact patient survival. METHODS: We combined liver single nucleus (snRNA-seq), single cell (scRNA-seq), and bulk RNA-sequencing (RNA-seq) data to search for cell-type differences in HCC. To first identify cell-types in HCC, adjacent non-tumor tissue, and normal liver, we integrated single-cell level data from a healthy liver cohort (n = 9 non-HCC samples) collected in the Strasbourg University Hospital; an HCC cohort (n = 1 non-HCC, n = 14 HCC-tumor, and n = 14 adjacent non-tumor samples) collected in the Singapore General Hospital and National University; and another HCC cohort (n = 3 HCC-tumor and n = 3 adjacent non-tumor samples) collected in the Dumont-UCLA Liver Cancer Center. We then leveraged these single cell level data to decompose the cell-types in liver bulk RNA-seq data from HCC patients' tumor (n = 361) and adjacent non-tumor tissue (n = 49) from the Cancer Genome Atlas (TCGA) multi-center cohort. For replication, we decomposed 221 HCC and 209 adjacent non-tumor liver microarray samples from the Liver Cancer Institute (LCI) cohort collected by the Liver Cancer Institute and Zhongshan Hospital of Fudan University. RESULTS: We discovered a tumor-associated proliferative cell-type, Prol (80.4% tumor cells), enriched for cell cycle and mitosis genes. In the liver bulk tissue from the TCGA cohort, the proportion of the Prol cell-type is significantly increased in HCC and associates with a worse overall survival. Independently from our decomposition analysis, we reciprocally show that Prol nuclei/cells significantly over-express both tumor-elevated and survival-decreasing genes obtained from the bulk tissue. Our replication analysis in the LCI cohort confirmed that an increased estimated proportion of the Prol cell-type in HCC is a significant marker for a shorter overall survival. Finally, we show that somatic mutations in the tumor suppressor genes TP53 and RB1 are linked to an increase of the Prol cell-type in HCC. CONCLUSIONS: By integrating liver single cell, single nucleus, and bulk expression data from multiple cohorts we identified a proliferating cell-type (Prol) enriched in HCC tumors, associated with a decreased overall survival, and linked to TP53 and RB1 somatic mutations.

Identification of 90 NAFLD GWAS loci and establishment of NAFLD PRS and causal role of NAFLD in coronary artery disease.

The prevalence of non-alcoholic fatty liver disease (NAFLD), now also known as metabolic dysfunction-associated fatty liver disease (MAFLD), is rapidly increasing worldwide due to the ongoing obesity epidemic. However, currently the NALFD diagnosis requires non-readily available imaging technologies or liver biopsy, which has drastically limited the sample sizes of NAFLD studies and hampered the discovery of its genetic component. Here we utilized the large UK Biobank (UKB) to accurately estimate the NAFLD status in UKB based on common serum traits and anthropometric measures. Scoring all individuals in UKB for NAFLD risk resulted in 28,396 NAFLD cases and 108,652 healthy individuals at a >90% confidence level. Using this imputed NAFLD status to perform the largest NAFLD genome-wide association study (GWAS) to date, we identified 94 independent (R2 < 0.2) NAFLD GWAS loci, of which 90 have not been identified before; built a polygenic risk score (PRS) model to predict the genetic risk of NAFLD; and used the GWAS variants of imputed NAFLD for a tissue-aware Mendelian randomization analysis that discovered a significant causal effect of NAFLD on coronary artery disease (CAD). In summary, we accurately estimated the NAFLD status in UKB using common serum traits and anthropometric measures, which empowered us to identify 90 GWAS NAFLD loci, build NAFLD PRS, and discover a significant causal effect of NAFLD on CAD.