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.

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.

BACKGROUND AND AIMS: Gallstone disease (GSD) associates with significant morbidity and mortality. Decreased secretion of bile acids has been suggested as a driving factor for GSD. Recently, we linked the protein phosphatase 1 regulatory subunit 3 beta (PPP1R3B) rs4240624 genotype to decreased bile acid levels in bile. In this study, we investigated whether these individuals had an increased risk for GSD as well as the differences in the lipid composition of the gallbladder bile of these individuals compared to controls and patients with GSD. METHODS: Bile acids, cholesterol, and phospholipid levels in gallbladder bile samples were enzymatically measured in 46 patients (34 female, age 45.7 ± 9.8 years, BMI 41.3 ± 4.4 kg/m2) who underwent elective laparoscopic Roux-en-Y gastric bypass. The lipidome of gallbladder bile was analyzed using high-performance liquid chromatography-mass spectrometry. Gallstone status was evaluated using abdominal ultrasonography before the surgery. RESULTS: The G allele of PPP1R3B rs4240624 was significantly associated with GSD in patients with obesity. We validated this association in the UK Biobank. Bile lipidomics demonstrated that 13 of the 17 minor lipid classes measured were higher in individuals with the G allele. The concentrations of bile acids, cholesterol, and phospholipids, as well as the cholesterol saturation index, were lower in patients with GSD than in those without gallstones. GSD had an effect similar to that of PPP1R3B genotype on minor lipids. CONCLUSION: The PPP1R3B rs4240624 genotype is associated with gallstones and with changes in gallbladder bile similar to those observed in patients with gallstones, suggesting that the PPP1R3B genotype contributes to the risk of gallstones by altering the bile lipidome.

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.

Mechanisms of abdominal obesity GWAS variants have remained largely unknown. To elucidate these mechanisms, we leveraged subcutaneous adipose tissue (SAT) single nucleus RNA-sequencing and genomics data. After discovering that heritability of abdominal obesity is enriched in adipocytes, we focused on a SAT unique adipocyte marker gene, the transcription factor TBX15, and its abdominal obesity-associated deleterious missense variant, rs10494217. The allele frequency of rs10494217 revealed a north-to-south decreasing gradient, with consistent significant FST values observed for 25 different populations when compared to Finns, a population with a history of genetic isolation. Given the role of Tbx15 in mouse thermogenesis, the frequency may have increased as an adaptation to cold in Finns. Our selection analysis provided significant evidence of selection for the abdominal obesity risk allele T of rs10494217 in Finns, with a north-to-south decreasing trend in other populations, and demonstrated that latitude significantly predicts the allele frequency. We also discovered that the risk allele status significantly affects SAT adipocyte expression of multiple adipocyte marker genes in trans in two cohorts. Two of these trans genes have been connected to thermogenesis, supporting the thermogenic effect of the TBX15 missense variant as a possible cause of its selection. Adipose expression of one trans gene, a lncRNA, AC002066.1, was strongly associated with adipocyte size, implicating it in metabolically unhealthy adipocyte hypertrophy. In summary, the abdominal obesity variant rs10494217 was selected in Finns, and individuals with the risk allele have trans effects on adipocyte expression of genes relating to thermogenesis and adipocyte hypertrophy.

OBJECTIVE: Metabolic dysfunction-associated steatotic liver disease (MASLD), formerly known as non-alcoholic fatty liver disease (NAFLD), is the most prevalent liver disease globally, yet no therapies are approved. The effects of Escherichia coli Nissle 1917 expressing aldafermin, an engineered analog of the intestinal hormone FGF19, in combination with dietary change were investigated as a potential treatment for MASLD. METHODS: MASLD was induced in C57BL/6J male mice by American lifestyle-induced obesity syndrome diet and then switched to a standard chow diet for seven weeks. In addition to the dietary change, the intervention group received genetically engineered E. coli Nissle expressing aldafermin, while control groups received either E. coli Nissle vehicle or no treatment. MASLD-related plasma biomarkers were measured using an automated clinical chemistry analyzer. The liver steatosis was assessed by histology and bioimaging analysis using Fiji (ImageJ) software. The effects of the intervention in the liver were also evaluated by RNA sequencing and liquid-chromatography-based non-targeted metabolomics analysis. Pathway enrichment studies were conducted by integrating the differentially expressed genes from the transcriptomics findings with the metabolites from the metabolomics results using Ingenuity pathway analysis. RESULTS: After the intervention, E. coli Nissle expressing aldafermin along with dietary changes reduced body weight, liver steatosis, plasma aspartate aminotransferase, and plasma cholesterol levels compared to the two control groups. The integration of transcriptomics with non-targeted metabolomics analysis revealed the downregulation of amino acid metabolism and related receptor signaling pathways potentially implicated in the reduction of hepatic steatosis and insulin resistance. Moreover, the downregulation of pathways linked to lipid metabolism and changes in amino acid-related pathways suggested an overall reduction of oxidative stress in the liver. CONCLUSIONS: These data support the potential for using engineered microbial therapeutics in combination with dietary changes for managing MASLD.

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.

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.

Background

Obesity predisposes individuals to multiple cardiometabolic disorders, including type 2 diabetes (T2D). As body mass index (BMI) cannot reliably differentiate fat from lean mass, the metabolically detrimental abdominal obesity has been estimated using waist-hip ratio (WHR). Waist-hip ratio adjusted for body mass index (WHRadjBMI) in turn is a well-established sex-specific marker for abdominal fat and adiposity, and a predictor of adverse metabolic outcomes, such as T2D. However, the underlying genes and regulatory mechanisms orchestrating the sex differences in obesity and body fat distribution in humans are not well understood.

Methods

We searched for genetic master regulators of WHRadjBMI by employing integrative genomics approaches on human subcutaneous adipose RNA sequencing (RNA-seq) data (n ~ 1400) and WHRadjBMI GWAS data (n ~ 700,000) from the WHRadjBMI GWAS cohorts and the UK Biobank (UKB), using co-expression network, transcriptome-wide association study (TWAS), and polygenic risk score (PRS) approaches. Finally, we functionally verified our genomic results using gene knockdown experiments in a human primary cell type that is critical for adipose tissue function.

Results

Here, we identified an adipose gene co-expression network that contains 35 obesity GWAS genes and explains a significant amount of polygenic risk for abdominal obesity and T2D in the UKB (n = 392,551) in a sex-dependent way. We showed that this network is preserved in the adipose tissue data from the Finnish Kuopio Obesity Study and Mexican Obesity Study. The network is controlled by a novel adipose master transcription factor (TF), TBX15, a WHRadjBMI GWAS gene that regulates the network in trans. Knockdown of TBX15 in human primary preadipocytes resulted in changes in expression of 130 network genes, including the key adipose TFs, PPARG and KLF15, which were significantly impacted (FDR < 0.05), thus functionally verifying the trans regulatory effect of TBX15 on the WHRadjBMI co-expression network.

Conclusions

Our study discovers a novel key function for the TBX15 TF in trans regulating an adipose co-expression network of 347 adipose, mitochondrial, and metabolically important genes, including PPARG, KLF15, PPARA, ADIPOQ, and 35 obesity GWAS genes. Thus, based on our converging genomic, transcriptional, and functional evidence, we interpret the role of TBX15 to be a main transcriptional regulator in the adipose tissue and discover its importance in human abdominal obesity.

BACKGROUND: Obesity is a worldwide epidemic characterized by adipose tissue (AT) inflammation. AT is also a source of extracellular vesicles (EVs) that have recently been implicated in disorders related to metabolic syndrome. However, our understanding of mechanistic aspect of obesitys impact on EV secretion from human AT remains limited. METHODS: We investigated EVs from human Simpson Golabi Behmel Syndrome (SGBS) adipocytes, and from AT as well as plasma of subjects undergoing bariatric surgery. SGBS cells were treated with TNFα, palmitic acid, and eicosapentaenoic acid. Various analyses, including nanoparticle tracking analysis, electron microscopy, high-resolution confocal microscopy, and gas chromatography-mass spectrometry, were utilized to study EVs. Plasma EVs were analyzed with imaging flow cytometry. RESULTS: EVs from mature SGBS cells differed significantly in size and quantity compared to preadipocytes, disagreeing with previous findings in mouse adipocytes and indicating that adipogenesis promotes EV secretion in human adipocytes. Inflammatory stimuli also induced EV secretion, and altered EV fatty acid (FA) profiles more than those of cells, suggesting the role of EVs as rapid responders to metabolic shifts. Visceral AT (VAT) exhibited higher EV secretion compared to subcutaneous AT (SAT), with VAT EV counts positively correlating with plasma triacylglycerol (TAG) levels. Notably, the plasma EVs of subjects with obesity contained a higher number of adiponectin-positive EVs than those of lean subjects, further demonstrating higher AT EV secretion in obesity. Moreover, plasma EV counts of people with obesity positively correlated with body mass index and TNF expression in SAT, connecting increased EV secretion with AT expansion and inflammation. Finally, EVs from SGBS adipocytes and AT contained TAGs, and EV secretion increased despite signs of less active lipolytic pathways, indicating that AT EVs could be involved in the mobilization of excess lipids into circulation. CONCLUSIONS: We are the first to provide detailed FA profiles of human AT EVs. We report that AT EV secretion increases in human obesity, implicating their role in TAG transport and association with adverse metabolic parameters, thereby emphasizing their role in metabolic disorders. These findings promote our understanding of the roles that EVs play in human AT biology and metabolic disorders.