FoxK1 associated gene regulatory network in hepatic insulin action and its relationship to FoxO1 and insulin receptor mediated transcriptional regulation

OBJECTIVE

Insulin acts on the liver via changes in gene expression to maintain glucose and lipid homeostasis. This study aimed to the Forkhead box protein K1 (FOXK1) associated gene regulatory network as a transcriptional regulator of hepatic insulin action and to determine its role versus FoxO1 and possible actions of the insulin receptor at the DNA level.

METHODS

Genome-wide analysis of FoxK1 binding were studied by chromatin immunoprecipitation sequencing and compared to those for IR and FoxO1. These were validated by knockdown experiments and gene expression analysis.

RESULTS

Chromatin immunoprecipitation (ChIP) sequencing shows that FoxK1 binds to the proximal promoters and enhancers of over 4000 genes, and insulin enhances this interaction for about 75% of them. These include genes involved in cell cycle, senescence, steroid biosynthesis, autophagy, and metabolic regulation, including glucose metabolism and mitochondrial function and are enriched in a TGTTTAC consensus motif. Some of these genes are also bound by FoxO1. Comparing this FoxK1 ChIP-seq data to that of the insulin receptor (IR) reveals that FoxK1 may act as the transcription factor partner for some of the previously reported roles of IR in gene regulation, including for LARS1 and TIMM22, which are involved in rRNA processing and cell cycle.

CONCLUSION

These data demonstrate that FoxK1 is an important regulator of gene expression in response to insulin in liver and may act in concert with FoxO1 and IR in regulation of genes in metabolism and other important biological pathways.

Sustained hyperglycemia specifically targets translation of mRNAs for insulin secretion

Pancreatic β cells are specialized for coupling glucose metabolism to insulin peptide production and secretion. Acute glucose exposure robustly and coordinately increases translation of proinsulin and proteins required for secretion of mature insulin peptide. By contrast, chronically elevated glucose levels that occur during diabetes impair β cell insulin secretion and have been shown experimentally to suppress insulin translation. Whether translation of other genes critical for insulin secretion is similarly downregulated by chronic high glucose is unknown. Here, we used high-throughput ribosome profiling and nascent proteomics in MIN6 insulinoma cells to elucidate the genome-wide impact of sustained high glucose on β cell mRNA translation. Before induction of ER stress or suppression of global translation, sustained high glucose suppressed glucose-stimulated insulin secretion and downregulated translation of not only insulin, but also mRNAs related to insulin secretory granule formation, exocytosis, and metabolism-coupled insulin secretion. Translation of these mRNAs was also downregulated in primary rat and human islets following ex vivo incubation with sustained high glucose and in an in vivo model of chronic mild hyperglycemia. Furthermore, translational downregulation decreased cellular abundance of these proteins. Our study uncovered a translational regulatory circuit during β cell glucose toxicity that impairs expression of proteins with critical roles in β cell function.

Sustained hyperglycemia specifically targets translation of mRNAs for insulin secretion

Pancreatic β-cells are specialized for coupling glucose metabolism to insulin peptide production and secretion. Acute glucose exposure robustly and coordinately increases translation of proinsulin and proteins required for secretion of mature insulin peptide. By contrast, chronically elevated glucose levels that occur during diabetes impair β-cell insulin secretion and have been shown experimentally to suppress insulin translation. Whether translation of other genes critical for insulin secretion are similarly downregulated by chronic high glucose is unknown. Here, we used high-throughput ribosome profiling and nascent proteomics in MIN6 insulinoma cells to elucidate the genome-wide impact of sustained high glucose on β-cell mRNA translation. Prior to induction of ER stress or suppression of global translation, sustained high glucose suppressed glucose-stimulated insulin secretion and downregulated translation of not only insulin, but also of mRNAs related to insulin secretory granule formation, exocytosis, and metabolism-coupled insulin secretion. Translation of these mRNAs was also downregulated in primary rat and human islets following ex-vivo incubation with sustained high glucose and in an in vivo model of chronic mild hyperglycemia. Furthermore, translational downregulation decreased cellular abundance of these proteins. Our findings uncover a translational regulatory circuit during β-cell glucose toxicity that impairs expression of proteins with critical roles in β-cell function.

Fructose induced KHK-C can increase ER stress independent of its effect on lipogenesis to drive liver disease in diet-induced and genetic models of NAFLD

Non-alcoholic fatty liver disease (NAFLD) is a liver manifestation of metabolic syndrome, and is estimated to affect one billion individuals worldwide. An increased intake of a high-fat diet (HFD) and sugar-sweetened beverages are risk-factors for NAFLD development, but how their combined intake promotes progression to a more severe form of liver injury is unknown. Here we show that fructose metabolism via ketohexokinase (KHK) C isoform leads to unresolved endoplasmic reticulum (ER) stress when coupled with a HFD intake. Conversely, a liver-specific knockdown of KHK in mice consuming fructose on a HFD is adequate to improve the NAFLD activity score and exert a profound effect on the hepatic transcriptome. Overexpression of KHK-C in cultured hepatocytes is sufficient to induce ER stress in fructose free media. Upregulation of KHK-C is also observed in mice with genetically induced obesity or metabolic dysfunction, whereas KHK knockdown in these mice improves metabolic function. Additionally, in over 100 inbred strains of male or female mice hepatic KHK expression correlates positively with adiposity, insulin resistance, and liver triglycerides. Similarly, in 241 human subjects and their controls, hepatic Khk expression is upregulated in early, but not late stages of NAFLD. In summary, we describe a novel role of KHK-C in triggering ER stress, which offers a mechanistic understanding of how the combined intake of fructose and a HFD propagates the development of metabolic complications.

Continuous glucose monitoring in patients with post-bariatric hypoglycaemia reduces hypoglycaemia and glycaemic variability

AIM

To determine whether continuous glucose monitoring (CGM) can reduce hypoglycaemia in patients with post-bariatric hypoglycaemia (PBH).

MATERIALS AND METHODS

In an open-label, nonrandomized, pre-post design with sequential assignment, CGM data were collected in 22 individuals with PBH in two sequential phases: (i) masked (no access to sensor glucose or alarms); and (ii) unmasked (access to sensor glucose and alarms for low or rapidly declining sensor glucose). Twelve participants wore the Dexcom G4 device for a total of 28 days, while 10 wore the Dexcom G6 device for a total of 20 days.

RESULTS

Participants with PBH spent a lower percentage of time in hypoglycaemia over 24 hours with unmasked versus masked CGM (<3.3 mM/L, or <60 mg/dL: median [median absolute deviation {MAD}] 0.7 [0.8]% vs. 1.4 [1.7]%, P = 0.03; <3.9 mM/L, or <70 mg/dL: median [MAD] 2.9 [2.5]% vs. 4.7 [4.8]%; P = 0.04), with similar trends overnight. Sensor glucose data from the unmasked phase showed a greater percentage of time spent between 3.9 and 10 mM/L (70-180 mg/dL) (median [MAD] 94.8 [3.9]% vs. 90.8 [5.2]%; P = 0.004) and lower glycaemic variability over 24 hours (median [MAD] mean amplitude of glycaemic excursion 4.1 [0.98] vs. 4.4 [0.99] mM/L; P = 0.04). During the day, participants also spent a greater percentage of time in normoglycaemia with unmasked CGM (median [MAD] 94.2 [4.8]% vs. 90.9 [6.2]%; P = 0.005), largely due to a reduction in hyperglycaemia (>10 mM/L, or 180 mg/dL: median [MAD] 1.9 [2.2]% vs. 3.9 [3.6]%; P = 0.02).

CONCLUSIONS

Real-time CGM data and alarms are associated with reductions in low sensor glucose, elevated sensor glucose, and glycaemic variability. This suggests CGM allows patients to detect hyperglycaemic peaks and imminent hypoglycaemia, allowing dietary modification and self-treatment to reduce hypoglycaemia. The use of CGM devices may improve safety in PBH, particularly for patients with hypoglycaemia unawareness.

ScreenDMT reveals linoleic acid diols replicably associate with BMI and stimulate adipocyte calcium fluxes

Activating brown adipose tissue (BAT) improves systemic metabolism, making it a promising target for metabolic syndrome. BAT is activated by 12, 13-dihydroxy-9Z-octadecenoic acid (12, 13-diHOME), which we previously identified to be inversely associated with BMI and which directly improves metabolism in multiple tissues. Here we profile plasma lipidomics from a cohort of 83 people and test which lipids' association with BMI replicates in a concordant direction using our novel tool ScreenDMT, whose power and validity we demonstrate via mathematical proofs and simulations. We find that the linoleic acid diols 12, 13-diHOME and 9, 10-diHOME both replicably inversely associate with BMI and mechanistically activate calcium fluxes in mouse brown and white adipocytes in vitro, which implicates this pathway and 9, 10-diHOME as candidate therapeutic targets. ScreenDMT can be applied to test directional mediation, directional replication, and qualitative interactions, such as identifying biomarkers whose association is shared (replication) or opposite (qualitative interaction) across diverse populations.

Leucine-973 is a crucial residue differentiating insulin and IGF-1 receptor signaling

Insulin and IGF-1 receptors (IR and IGF1R) are highly homologous and share similar signaling systems, but each has a unique physiological role, with IR primarily regulating metabolic homeostasis and IGF1R regulating mitogenic control and growth. Here, we show that replacement of a single amino acid at position 973, just distal to the NPEY motif in the intracellular juxtamembrane region, from leucine, which is highly conserved in IRs, to phenylalanine, the highly conserved homologous residue in IGF1Rs, resulted in decreased IRS-1/PI3K/Akt/mTORC1 signaling and increased Shc/Gab1/MAPK cell cycle signaling. As a result, cells expressing L973F-IR exhibited decreased insulin-induced glucose uptake, increased cell growth, and impaired receptor internalization. Mice with knockin of the L973F-IR showed similar alterations in signaling in vivo, and this led to decreased insulin sensitivity, a modest increase in growth, and decreased weight gain when mice were challenged with a high-fat diet. Thus, leucine-973 in the juxtamembrane region of the IR acts as a crucial residue differentiating IR signaling from IGF1R signaling.

Fructose Induced KHK-C Increases ER Stress and Modulates Hepatic Transcriptome to Drive Liver Disease in Diet-Induced and Genetic Models of NAFLD

Non-alcoholic fatty liver disease (NAFLD) is a liver manifestation of metabolic syndrome, and is estimated to affect one billion individuals worldwide. An increased intake of a high-fat diet (HFD) and sugar-sweetened beverages are risk-factors for NAFLD development, but how their combined intake promotes progression to a more severe form of liver injury is unknown. Here we show that fructose metabolism via ketohexokinase (KHK) C isoform increases endoplasmic reticulum (ER) stress in a dose dependent fashion, so when fructose is coupled with a HFD intake it leads to unresolved ER stress. Conversely, a liver-specific knockdown of KHK in C57BL/6J male mice consuming fructose on a HFD is adequate to improve the NAFLD activity score and exert a profound effect on the hepatic transcriptome. Overexpression of KHK-C in cultured hepatocytes is sufficient to induce ER stress in fructose free media. Upregulation of KHK-C is also observed in genetically obesity ob/ob, db/db and lipodystrophic FIRKO male mice, whereas KHK knockdown in these mice improves metabolic function. Additionally, in over 100 inbred strains of male or female mice hepatic KHK expression correlates positively with adiposity, insulin resistance, and liver triglycerides. Similarly, in 241 human subjects and their controls, hepatic Khk expression is upregulated in early, but not late stages of NAFLD. In summary, we describe a novel role of KHK-C in triggering ER stress, which offers a mechanistic understanding of how the combined intake of fructose and a HFD propagates the development of metabolic complications.

Unique ligand and kinase-independent roles of the insulin receptor in regulation of cell cycle, senescence and apoptosis

Insulin acts through the insulin receptor (IR) tyrosine kinase to exert its classical metabolic and mitogenic actions. Here, using receptors with either short or long deletion of the β-subunit or mutation of the kinase active site (K1030R), we have uncovered a second, previously unrecognized IR signaling pathway that is intracellular domain-dependent, but ligand and tyrosine kinase-independent (LYK-I). These LYK-I actions of the IR are linked to changes in phosphorylation of a network of proteins involved in the regulation of extracellular matrix organization, cell cycle, ATM signaling and cellular senescence; and result in upregulation of expression of multiple extracellular matrix-related genes and proteins, down-regulation of immune/interferon-related genes and proteins, and increased sensitivity to apoptosis. Thus, in addition to classical ligand and tyrosine kinase-dependent (LYK-D) signaling, the IR regulates a second, ligand and tyrosine kinase-independent (LYK-I) pathway, which regulates the cellular machinery involved in senescence, matrix interaction and response to extrinsic challenges.

E2F1 transcription factor mediates a link between fat and islets to promote β cell proliferation in response to acute insulin resistance

Prevention or amelioration of declining β cell mass is a potential strategy to cure diabetes. Here, we report the pathways utilized by β cells to robustly replicate in response to acute insulin resistance induced by S961, a pharmacological insulin receptor antagonist. Interestingly, pathways that include CENP-A and the transcription factor E2F1 that are independent of insulin signaling and its substrates appeared to mediate S961-induced β cell multiplication. Consistently, pharmacological inhibition of E2F1 blocks β-cell proliferation in S961-injected mice. Serum from S961-treated mice recapitulates replication of β cells in mouse and human islets in an E2F1-dependent manner. Co-culture of islets with adipocytes isolated from S961-treated mice enables β cells to duplicate, while E2F1 inhibition limits their growth even in the presence of adipocytes. These data suggest insulin resistance-induced proliferative signals from adipocytes activate E2F1, a potential therapeutic target, to promote β cell compensation.

Transcriptomic Regulation of Muscle Mitochondria and Calcium Signaling by Insulin/IGF-1 Receptors Depends on FoxO Transcription Factors

Insulin and IGF-1, acting through the insulin receptor (IR) and IGF-1 receptor (IGF1R), maintain muscle mass and mitochondrial function, at least part of which occurs via their action to regulate gene expression. Here, we show that while muscle-specific deletion of IR or IGF1R individually results in only modest changes in the muscle transcriptome, combined deletion of IR/IGF1R (MIGIRKO) altered > 3000 genes, including genes involved in mitochondrial dysfunction, fibrosis, cardiac hypertrophy, and pathways related to estrogen receptor, protein kinase A (PKA), and calcium signaling. Functionally, this was associated with decreased mitochondrial respiration and increased ROS production in MIGIRKO muscle. To determine the role of FoxOs in these changes, we performed RNA-Seq on mice with muscle-specific deletion of FoxO1/3/4 (M-FoxO TKO) or combined deletion of IR, IGF1R, and FoxO1/3/4 in a muscle quintuple knockout (M-QKO). This revealed that among IR/IGF1R regulated genes, >97% were FoxO-dependent, and their expression was normalized in M-FoxO TKO and M-QKO muscle. FoxO-dependent genes were related to oxidative phosphorylation, inflammatory signaling, and TCA cycle. Metabolomic analysis showed accumulation of TCA cycle metabolites in MIGIRKO, which was reversed in M-QKO muscle. Likewise, calcium signaling genes involved in PKA signaling and sarcoplasmic reticulum calcium homeostasis were markedly altered in MIGIRKO muscle but normalized in M-QKO. Thus, combined loss of insulin and IGF-1 action in muscle transcriptionally alters mitochondrial function and multiple regulatory and signaling pathways, and these changes are mediated by FoxO transcription factors.

Signaling defects associated with insulin resistance in nondiabetic and diabetic individuals and modification by sex

Insulin resistance is present in one-quarter of the general population, predisposing these people to a wide range of diseases. Our aim was to identify cell-intrinsic determinants of insulin resistance in this population using induced pluripotent stem cell-derived (iPSC-derived) myoblasts (iMyos). We found that these cells exhibited a large network of altered protein phosphorylation in vitro. Integrating these data with data from type 2 diabetic iMyos revealed critical sites of conserved altered phosphorylation in IRS-1, AKT, mTOR, and TBC1D1 in addition to changes in protein phosphorylation involved in Rho/Rac signaling, chromatin organization, and RNA processing. There were also striking differences in the phosphoproteome in cells from men versus women. These sex-specific and insulin-resistance defects were linked to functional differences in downstream actions. Thus, there are cell-autonomous signaling alterations associated with insulin resistance within the general population and important differences between men and women, many of which also occur in diabetes, that contribute to differences in physiology and disease.

Using single-nucleus RNA-sequencing to interrogate transcriptomic profiles of archived human pancreatic islets

Background

Human pancreatic islets are a central focus of research in metabolic studies. Transcriptomics is frequently used to interrogate alterations in cultured human islet cells using single-cell RNA-sequencing (scRNA-seq). We introduce single-nucleus RNA-sequencing (snRNA-seq) as an alternative approach for investigating transplanted human islets.

Methods

The Nuclei EZ protocol was used to obtain nuclear preparations from fresh and frozen human islet cells. Such preparations were first used to generate snRNA-seq datasets and compared to scRNA-seq output obtained from cells from the same donor. Finally, we employed snRNA-seq to obtain the transcriptomic profile of archived human islets engrafted in immunodeficient animals.

Results

We observed virtually complete concordance in identifying cell types and gene proportions as well as a strong association of global and islet cell type gene signatures between scRNA-seq and snRNA-seq applied to fresh and frozen cultured or transplanted human islet samples.

Conclusions

We propose snRNA-seq as a reliable strategy to probe transcriptomic profiles of freshly harvested or frozen sources of transplanted human islet cells especially when scRNA-seq is not ideal.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-021-00941-8.

Glycemic Patterns Are Distinct in Post-Bariatric Hypoglycemia After Gastric Bypass (PBH-RYGB)

CONTEXT

Severe hypoglycemia with neuroglycopenia, termed post-bariatric hypoglycemia (PBH). typically occurs postprandially, but it is also reported after activity or mid-nocturnally.

OBJECTIVE

To quantify glycemia, glycemic variability, and magnitude/duration of low sensor glucose (SG) values in patients with PBH after Roux-en-Y gastric bypass (PBH-RYGB).

METHODS

This retrospective analysis of data from an academic medical center included individuals with PBH-RYGB (n = 40), reactive hypoglycemia without gastrointestinal surgery (Non-Surg Hypo, n = 20), prediabetes (Pre-DM, n = 14), newly diagnosed T2D (n = 5), and healthy controls (HC, n = 38). Masked continuous glucose monitoring (Dexcom G4) was used to assess patterns over 24 hours, daytime (6 am-midnight) and nighttime (midnight-6 am). Prespecified measures included mean and median SG, variability, and percent time at thresholds of sensor glucose.

RESULTS

Mean and median SG were similar for PBH-RYGB and HC (mean: 99.8 ± 18.6 vs 96.9 ± 10.2 mg/dL; median: 93.0 ± 14.8 vs 94.5 ± 7.4 mg/dL). PBH-RYGB had a higher coefficient of variation (27.3 ± 6.8 vs 17.9 ± 2.4%, P < 0.0001) and range (154.5 ± 50.4 vs 112.0 ± 26.7 mg/dL, P < 0.0001). Nadir was lowest in PBH-RYGB (42.5 ± 3.7 vs HC 49.0 ± 11.9 mg/dL, P = 0.0046), with >2-fold greater time with SG < 70 mg/dL vs HC (7.7 ± 8.4 vs 3.2 ± 4.1%, P = 0.0013); these differences were greater at night (12.6 ± 16.9 vs 1.0 ± 1.5%, P < 0.0001). Non-Surg Hypo also had 4-fold greater time with SG < 70 at night vs HC (SG < 70: 4.0 ± 5.9% vs 1.0 ± 1.5%), but glycemic variability was not increased.

CONCLUSION

Patients with PBH-RYGB experience higher glycemic variability and frequency of SG < 70 compared to HC, especially at night. These data suggest that additional pathophysiologic mechanisms beyond prandial changes contribute to PBH.

Insulin Resistance and Gender Define a Cell Autonomous Supernetworkof Protein Phosphorylation

Many hormones and growth factors, including insulin, act through networks of protein phosphorylation. Insulin resistance is an important factor in the pathophysiology of many metabolic disorders. The aim of this study was to uncover the cell autonomous determinants of insulin action and protein phosphorylation using induced pluripotent stem cell (iPSC)-derived myoblasts (iMyos) in vitro. Here, we show that iMyos from non-diabetic individuals in the highest quintile of insulin resistance show impaired insulin signaling, defective insulin-stimulated glucose uptake and decreased glycogen synthase activity compared to iMyos from the insulin sensitive individuals, indicating these cells mirror in vitro the alterations seen in vivo. Global phosphoproteomic analysis uncovered a large network of proteins whose phosphorylation was altered in association with insulin resistance, most outside the canonical insulin-signaling cascade. More surprisingly, we also observed striking differences in the phosphoproteomic signature of iMyos derived from male versus female subjects, involving multiple pathways regulating diverse cellular functions, including DNA and RNA processing, GTPase signaling, and SUMOylation/ubiquitination. These findings provide new insights into the cell autonomous mechanisms underlying insulin resistance in the non-diabetic population and provide evidence of a major, previously unrecognized, supernetwork of cell signaling differences in males and females that must be considered in understanding the molecular basis of sex-based differences in normal physiology and disease.

A Cell-Autonomous Signature of Dysregulated Protein Phosphorylation Underlies Muscle Insulin Resistance in Type 2 Diabetes

Skeletal muscle insulin resistance is the earliest defect in type 2 diabetes (T2D), preceding and predicting disease development. To what extent this reflects a primary defect or is secondary to tissue cross talk due to changes in hormones or circulating metabolites is unknown. To address this question, we have developed an in vitro disease-in-a-dish model using iPS cells from T2D patients differentiated into myoblasts (iMyos). We find that T2D iMyos in culture exhibit multiple defects mirroring human disease, including an altered insulin signaling, decreased insulin-stimulated glucose uptake, and reduced mitochondrial oxidation. More strikingly, global phosphoproteomic analysis reveals a multidimensional network of signaling defects in T2D iMyos going beyond the canonical insulin-signaling cascade, including proteins involved in regulation of Rho GTPases, mRNA splicing and/or processing, vesicular trafficking, gene transcription, and chromatin remodeling. These cell-autonomous defects and the dysregulated network of protein phosphorylation reveal a new dimension in the cellular mechanisms underlying the fundamental defects in T2D.

Skeletal muscle expression of adipose-specific phospholipase in peripheral artery disease

Flow-limiting atherosclerotic lesions of arteries supplying the limbs are a cause of symptoms in patients with peripheral artery disease (PAD). Musculoskeletal metabolic factors also contribute to the pathophysiology of claudication, which is manifest as leg discomfort that impairs walking capacity. Accordingly, we conducted a case-control study to determine whether skeletal muscle metabolic gene expression is altered in PAD. Calf skeletal muscle gene expression of patients with PAD and healthy subjects was analyzed using microarrays. The top-ranking gene differentially expressed between PAD and controls (FDR < 0.001) was PLA2G16, which encodes adipose-specific phospholipase A2 (AdPLA) and is implicated in the maintenance of insulin sensitivity and regulation of lipid metabolism. Differential expression was confirmed by qRT-PCR; PLA2G16 was downregulated by 68% in patients with PAD (p < 0.001). Expression of Pla2g16 was then measured in control (db/+) and diabetic (db/db) mice that underwent unilateral femoral artery ligation. There was significantly reduced expression of Pla2g16 in the ischemic leg of both control and diabetic mice (by 51%), with significantly greater magnitude of reduction in the diabetic mice (by 79%). We conclude that AdPLA is downregulated in humans with PAD and in mice with hindlimb ischemia. Reduced AdPLA may contribute to impaired walking capacity in patients with PAD via its effects on skeletal muscle metabolism. Further studies are needed to fully characterize the role of AdPLA in PAD and to investigate its potential as a therapeutic target for alleviating symptoms of claudication.

Integrated metabolomics reveals altered lipid metabolism in adipose tissue in a model of extreme longevity

Adipose tissue plays an essential role in metabolic health. Ames dwarf mice are exceptionally long-lived and display metabolically beneficial phenotypes in their adipose tissue, providing an ideal model for studying the intersection between adipose tissue and longevity. To this end, we assessed the metabolome and lipidome of adipose tissue in Ames dwarf mice. We observed distinct lipid profiles in brown versus white adipose tissue of Ames dwarf mice that are consistent with increased thermogenesis and insulin sensitivity, such as increased cardiolipin and decreased ceramide concentrations. Moreover, we identified 5-hydroxyeicosapentaenoic acid (5-HEPE), an ω-3 fatty acid metabolite, to be increased in Ames dwarf brown adipose tissue (BAT), as well as in circulation. Importantly, 5-HEPE is increased in other models of BAT activation and is negatively correlated with body weight, insulin resistance, and circulating triglyceride concentrations in humans. Together, these data represent a novel lipid signature of adipose tissue in a mouse model of extreme longevity.

Association of Periconception Paternal Body Mass Index With Persistent Changes in DNA Methylation of Offspring in Childhood

IMPORTANCE

While prenatal nutrition and maternal obesity are recognized as important contributors to epigenetic changes and childhood obesity, the role of paternal obesity in the epigenome of offspring has not been well studied.

OBJECTIVES

To test whether periconception paternal body mass index (BMI) is associated with DNA methylation patterns in newborns, to examine associations between maternal and paternal BMI and the epigenome of offspring, and to examine persistence of epigenetic marks at ages 3 and 7 years.

DESIGN, SETTING, AND PARTICIPANTS

Project Viva is a prebirth cohort study of mothers and children including 2128 live births that enrolled mothers from April 1999 to July 2002 and followed offspring to adolescence. This study analyzed the subset of participants with available data on paternal BMI and DNA methylation in offspring blood in the newborn period, at age 3 years, and at age 7 years. Data were analyzed from July 2017 to October 2019.

EXPOSURES

The primary exposure was paternal periconception BMI; associations were adjusted for maternal prepregnancy BMI and stratified according to maternal BMI above or below 25.

MAIN OUTCOMES AND MEASURES

The primary outcome was genome-wide DNA methylation patterns in offspring blood collected at birth, age 3 years, and age 7 years.

RESULTS

A total of 429 father-mother-infant triads were included. The mean (SD) periconception paternal BMI was 26.4 (4.0) and mean maternal prepregnancy BMI was 24.5 (5.2); 268 fathers had BMI greater than or equal to 25 (mean [SD], 28.5 [3.3]) and 161 had BMI less than 25 (mean [SD], 22.8 [1.8]). Paternal BMI greater than or equal to 25 was associated with increased offspring birth weight compared with paternal BMI less than 25 (mean [SD] z score, 0.38 [0.91] vs 0.11 [0.96]; P = .004). Cord blood DNA methylation at 9 CpG sites was associated with paternal BMI independent of maternal BMI (q < .05). Methylation at cg04763273, between TFAP2C and BMP7, decreased by 5% in cord blood with every 1-unit increase in paternal BMI (P = 3.13 × 10-8); hypomethylation at this site persisted at ages 3 years and 7 years. Paternal BMI was associated with methylation at cg01029450 in the promoter region of the ARFGAP3 gene; methylation at this site was also associated with lower infant birth weight (β = -0.0003; SD = 0.0001; P = .03) and with higher BMI z score at age 3 years.

CONCLUSIONS AND RELEVANCE

In this study, paternal BMI was associated with DNA methylation, birth weight, and childhood BMI z score in offspring.

Characterization of Glycolytic Enzymes and Pyruvate Kinase M2 in Type 1 and 2 Diabetic Nephropathy

OBJECTIVE

Elevated glycolytic enzymes in renal glomeruli correlated with preservation of renal function in the Medalist Study, individuals with ≥50 years of type 1 diabetes. Specifically, pyruvate kinase M2 (PKM2) activation protected insulin-deficient diabetic mice from hyperglycemia-induced glomerular pathology. This study aims to extend these findings in a separate cohort of individuals with type 1 and type 2 diabetes and discover new circulatory biomarkers for renal protection through proteomics and metabolomics of Medalists' plasma. We hypothesize that increased glycolytic flux and improved mitochondrial biogenesis will halt the progression of diabetic nephropathy.

RESEARCH DESIGN AND METHODS

Immunoblots analyzed selected glycolytic and mitochondrial enzymes in postmortem glomeruli of non-Medalists with type 1 diabetes (n = 15), type 2 diabetes (n = 19), and no diabetes (n = 5). Plasma proteomic (SOMAscan) (n = 180) and metabolomic screens (n = 214) of Medalists with and without stage 3b chronic kidney disease (CKD) were conducted and significant markers validated by ELISA.

RESULTS

Glycolytic (PKM1, PKM2, and ENO1) and mitochondrial (MTCO2) enzymes were significantly elevated in glomeruli of CKD- versus CKD+ individuals with type 2 diabetes. Medalists' plasma PKM2 correlated with estimated glomerular filtration rate (r ² = 0.077; P = 0.0002). Several glucose and mitochondrial enzymes in circulation were upregulated with corresponding downregulation of toxic metabolites in CKD-protected Medalists. Amyloid precursor protein was also significantly upregulated, tumor necrosis factor receptors downregulated, and both confirmed by ELISA.

CONCLUSIONS

Elevation of enzymes involved in the metabolism of intracellular free glucose and its metabolites in renal glomeruli is connected to preserving kidney function in both type 1 and type 2 diabetes. The renal profile of elevated glycolytic enzymes and reduced toxic glucose metabolites is reflected in the circulation, supporting their use as biomarkers for endogenous renal protective factors in people with diabetes.

Plasma FGF-19 Levels are Increased in Patients with Post-Bariatric Hypoglycemia

BACKGROUND

Hypoglycemia is an increasingly recognized complication of bariatric surgery. Mechanisms contributing to glucose lowering remain incompletely understood. We aimed to identify differentially abundant plasma proteins in patients with post-bariatric hypoglycemia (PBH) after Roux-en-Y gastric bypass (RYGB), compared to asymptomatic post-RYGB.

METHODS

Proteomic analysis of blood samples collected after overnight fast and mixed meal challenge in individuals with PBH, asymptomatic RYGB, severe obesity, or overweight recruited from outpatient hypoglycemia or bariatric clinics.

RESULTS

The top-ranking differentially abundant protein at 120 min after mixed meal was fibroblast growth factor 19 (FGF-19), an intestinally derived hormone regulated by bile acid-FXR signaling; levels were 2.4-fold higher in PBH vs. asymptomatic post-RYGB (mean + SEM, 1094 ± 141 vs. 428 ± 45, P < 0.001, FDR < 0.01). FGF-19 ELISA confirmed 3.5-fold higher concentrations in PBH versus asymptomatic (360 ± 70 vs. 103 ± 18, P = 0.025). To explore potential links between increased FGF-19 and GLP-1, residual samples from other human studies in which GLP-1 was modulated were assayed. FGF-19 levels did not change in response to infusion of GLP-1 and PYY in overweight/obese individuals. Infusion of the GLP-1 receptor antagonist exendin 9-39 in recently operated asymptomatic post-RYGB did not alter FGF-19 levels after mixed meal. By contrast, GLP-1 receptor antagonist infusion yielded a significant increase in FGF-19 levels after oral glucose in individuals with PBH. While plasma bile acids did not differ between PBH and asymptomatic post-RYGB, these data suggest unique interrelationships between GLP-1 and FGF-19 in PBH.

CONCLUSIONS

Taken together, these data support FGF-19 as a potential contributor to insulin-independent pathways driving postprandial hypoglycemia in PBH.

SGLT2 inhibition reprograms systemic metabolism via FGF21-dependent and -independent mechanisms

Pharmacologic inhibition of the renal sodium/glucose cotransporter-2 induces glycosuria and reduces glycemia. Given that SGLT2 inhibitors (SGLT2i) reduce mortality and cardiovascular risk in type 2 diabetes, improved understanding of molecular mechanisms mediating these metabolic effects is required. Treatment of obese but nondiabetic mice with the SGLT2i canagliflozin (CANA) reduces adiposity, improves glucose tolerance despite reduced plasma insulin, increases plasma ketones, and improves plasma lipid profiles. Utilizing an integrated transcriptomic-metabolomics approach, we demonstrate that CANA modulates key nutrient-sensing pathways, with activation of 5' AMP-activated protein kinase (AMPK) and inhibition of mechanistic target of rapamycin (mTOR), independent of insulin or glucagon sensitivity or signaling. Moreover, CANA induces transcriptional reprogramming to activate catabolic pathways, increase fatty acid oxidation, reduce hepatic steatosis and diacylglycerol content, and increase hepatic and plasma levels of FGF21. Given that these phenotypes mirror the effects of FGF21 to promote lipid oxidation, ketogenesis, and reduction in adiposity, we hypothesized that FGF21 is required for CANA action. Using FGF21-null mice, we demonstrate that FGF21 is not required for SGLT2i-mediated induction of lipid oxidation and ketogenesis but is required for reduction in fat mass and activation of lipolysis. Taken together, these data demonstrate that SGLT2 inhibition triggers a fasting-like transcriptional and metabolic paradigm but requires FGF21 for reduction in adiposity.

FoxO Transcription Factors Are Critical Regulators of Diabetes-Related Muscle Atrophy

Insulin deficiency and uncontrolled diabetes lead to a catabolic state with decreased muscle strength, contributing to disease-related morbidity. FoxO transcription factors are suppressed by insulin and thus are key mediators of insulin action. To study their role in diabetic muscle wasting, we created mice with muscle-specific triple knockout of FoxO1/3/4 and induced diabetes in these M-FoxO-TKO mice with streptozotocin (STZ). Muscle mass and myofiber area were decreased 20-30% in STZ-Diabetes mice due to increased ubiquitin-proteasome degradation and autophagy alterations, characterized by increased LC3-containing vesicles, and elevated levels of phosphorylated ULK1 and LC3-II. Both the muscle loss and markers of increased degradation/autophagy were completely prevented in STZ FoxO-TKO mice. Transcriptomic analyses revealed FoxO-dependent increases in ubiquitin-mediated proteolysis pathways in STZ-Diabetes, including regulation of Fbxo32 (Atrogin1), Trim63 (MuRF1), Bnip3L, and Gabarapl. These same genes were increased 1.4- to 3.3-fold in muscle from humans with type 1 diabetes after short-term insulin deprivation. Thus, FoxO-regulated genes play a rate-limiting role in increased protein degradation and muscle atrophy in insulin-deficient diabetes.

TGF-β2 is an exercise-induced adipokine that regulates glucose and fatty acid metabolism

Exercise improves health and well-being across diverse organ systems, and elucidating mechanisms underlying the beneficial effects of exercise can lead to new therapies. Here, we show that transforming growth factor-β2 (TGF-β2) is secreted from adipose tissue in response to exercise and improves glucose tolerance in mice. We identify TGF-β2 as an exercise-induced adipokine in a gene expression analysis of human subcutaneous adipose tissue biopsies after exercise training. In mice, exercise training increases TGF-β2 in scWAT, serum, and its secretion from fat explants. Transplanting scWAT from exercise-trained wild type mice, but not from adipose tissue-specific Tgfb2−/− mice, into sedentary mice improves glucose tolerance. TGF-β2 treatment reverses the detrimental metabolic effects of high fat feeding in mice. Lactate, a metabolite released from muscle during exercise, stimulates TGF-β2 expression in human adipocytes. Administration of the lactate-lowering agent dichloroacetate during exercise training in mice decreases circulating TGF-β2 levels and reduces exercise-stimulated improvements in glucose tolerance. Thus, exercise training improves systemic metabolism through inter-organ communication with fat via a lactate-TGF-β2-signaling cycle.

The Crohn's disease polymorphism, ATG16L1 T300A, alters the gut microbiota and enhances the local Th1/Th17 response

Inflammatory bowel disease (IBD) is driven by dysfunction between host genetics, the microbiota, and immune system. Knowledge gaps remain regarding how IBD genetic risk loci drive gut microbiota changes. The Crohn's disease risk allele ATG16L1 T300A results in abnormal Paneth cells due to decreased selective autophagy, increased cytokine release, and decreased intracellular bacterial clearance. To unravel the effects of ATG16L1 T300A on the microbiota and immune system, we employed a gnotobiotic model using human fecal transfers into ATG16L1 T300A knock-in mice. We observed increases in Bacteroides ovatus and Th1 and Th17 cells in ATG16L1 T300A mice. Association of altered Schaedler flora mice with B. ovatus specifically increased Th17 cells selectively in ATG16L1 T300A knock-in mice. Changes occur before disease onset, suggesting that ATG16L1 T300A contributes to dysbiosis and immune infiltration prior to disease symptoms. Our work provides insight for future studies on IBD subtypes, IBD patient treatment and diagnostics.

Developmental and functional heterogeneity of white adipocytes within a single fat depot

Recent studies suggest that, even within a single adipose depot, there may be distinct subpopulations of adipocytes. To investigate this cellular heterogeneity, we have developed multiple conditionally immortalized clonal preadipocyte lines from white adipose tissue of mice. Analysis of these clones reveals at least three white adipocyte subpopulations. These subpopulations have differences in metabolism and differentially respond to inflammatory cytokines, insulin, and growth hormones. These also have distinct gene expression profiles and can be tracked by differential expression of three marker genes: Wilms' tumor 1, transgelin, and myxovirus 1. Lineage tracing analysis with dual-fluorescent reporter mice indicates that these adipocyte subpopulations have differences in gene expression and metabolism that mirror those observed in the clonal cell lines. Furthermore, preadipocytes and adipocytes from these subpopulations differ in their abundance in different fat depots. Thus, white adipose tissue, even in a single depot, is comprised of distinct subpopulations of white adipocytes with different physiological phenotypes. These differences in adipocyte composition may contribute to the differences in metabolic behavior and physiology of different fat depots.

Metabolic Effects of Betaine: A Randomized Clinical Trial of Betaine Supplementation in Prediabetes

CONTEXT

Plasma betaine correlates with insulin sensitivity in humans. Betaine supplementation improves metabolic effects in mice fed a high-fat diet.

OBJECTIVE

To assess metabolic effects of oral betaine in obese participants with prediabetes.

DESIGN

A 12-week, parallel arm, randomized, double-masked, placebo-controlled trial.

SETTING

University-affiliated hospital.

PARTICIPANTS AND INTERVENTIONS

Persons with obesity and prediabetes (N = 27) were randomly assigned to receive betaine 3300 mg orally twice daily for 10 days, then 4950 mg twice daily for 12 weeks, or placebo.

MAIN OUTCOME MEASURES

Changes from baseline in insulin sensitivity, glycemia, hepatic fat, and endothelial function.

RESULTS

There was a 16.5-fold increase in plasma dimethylglycine [dimethylglycine (DMG); P < 0.0001] levels, but modest 1.3- and 1.5-fold increases in downstream serine and methionine levels, respectively, in the betaine vs placebo arm. Betaine tended to reduce fasting glucose levels (P = 0.08 vs placebo) but had no other effect on glycemia. Insulin area under curve after oral glucose was reduced for betaine treatment compared with placebo (P = 0.038). Insulin sensitivity, assessed by euglycemic hyperinsulinemic clamp, was not improved. Serum total cholesterol levels increased after betaine treatment compared with placebo (P = 0.032). There were no differences in change in intrahepatic triglyceride or endothelial function between groups.

CONCLUSION

DMG accumulation supports DMG dehydrogenase as rate limiting for betaine metabolism in persons with prediabetes. Betaine had little metabolic effect. Additional studies may elucidate mechanisms contributing to differences between preclinical and human responses to betaine, and whether supplementation of metabolites downstream of DMG improves metabolism.

Integrating Extracellular Flux Measurements and Genome-Scale Modeling Reveals Differences between Brown and White Adipocytes

White adipocytes are specialized for energy storage, whereas brown adipocytes are specialized for energy expenditure. Explicating this difference can help identify therapeutic targets for obesity. A common tool to assess metabolic differences between such cells is the Seahorse Extracellular Flux (XF) Analyzer, which measures oxygen consumption and media acidification in the presence of different substrates and perturbagens. Here, we integrate the Analyzer's metabolic profile from human white and brown adipocytes with a genome-scale metabolic model to predict flux differences across the metabolic map. Predictions matched experimental data for the metabolite 4-aminobutyrate, the protein ABAT, and the fluxes for glucose, glutamine, and palmitate. We also uncovered a difference in how adipocytes dispose of nitrogenous waste, with brown adipocytes secreting less ammonia and more urea than white adipocytes. Thus, the method and software we developed allow for broader metabolic phenotyping and provide a distinct approach to uncovering metabolic differences.

YAP suppresses gluconeogenic gene expression through PGC1α

Cell growth and proliferation are tightly coupled to metabolism, and dissecting the signaling molecules which link these processes is an important step toward understanding development, regeneration, and cancer. The transcriptional regulator Yes-associated protein 1 (YAP) is a key regulator of liver size, development, and function. We now show that YAP can also suppress gluconeogenic gene expression. Yap deletion in primary hepatocytes potentiates the gluconeogenic gene response to glucagon and dexamethasone, whereas constitutively active YAP suppresses it. The effects of YAP are mediated by the transcriptional coactivator peroxisome proliferator-activated receptor-gamma coactivator 1. YAP inhibits its ability to bind to and activate transcription from the promoters of its gluconeogenic targets, and the effects of YAP are blunted upon its knockdown. In vivo, constitutively active YAP lowers plasma glucose levels and increases liver size.

CONCLUSION

YAP appears to reprogram cellular metabolism, diverting substrates away from the energy-consuming process of gluconeogenesis and toward the anabolic process of growth. (Hepatology 2017;66:2029-2041).

Age-associated NF-κB signaling in myofibers alters the satellite cell niche and re-strains muscle stem cell function

Skeletal muscle is a highly regenerative tissue, but muscle repair potential is increasingly compromised with advancing age. In this study, we demonstrate that increased NF-κB activity in aged muscle fibers contributes to diminished myogenic potential of their associated satellite cells. We further examine the impact of genetic modulation of NF-κB signaling in muscle satellite cells or myofibers on recovery after damage. These studies reveal that NF-κB activity in differentiated myofibers is sufficient to drive dysfunction of muscle regenerative cells via cell-non-autonomous mechanisms. Inhibition of NF-κB, or its downstream target Phospholipase A2, in myofibers rescued muscle regenerative potential in aged muscle. Moreover, systemic administration of sodium salicylate, an FDA-approved NF-κB inhibitor, decreased inflammatory gene expression and improved repair in aged muscle. Together, these studies identify a unique NF-κB regulated, non-cell autonomous mechanism by which stem cell function is linked to lipid signaling and homeostasis, and provide important new targets to stimulate muscle repair in aged individuals.

Evidence of stress in β cells obtained with laser capture microdissection from pancreases of brain dead donors

Isolated islets used for transplantation are known to be stressed, which can result from the circumstances of death, in particular brain death, the preservation of the pancreas with its warm and cold ischemia, from the trauma of the isolation process, and the complex events that occur during tissue culture. The current study focused upon the events that occur before the islet isolation procedure. Pancreases were obtained from brain dead donors (n = 7) with mean age 50 (11) and normal pancreatic tissue obtained at surgery done for pancreatic neoplasms (n = 7), mean age 69 (9). Frozen sections were subjected to laser capture microdissection (LCM) to obtain β-cell rich islet tissue, from which extracted RNA was analyzed with microarrays. Gene expression of the 2 groups was evaluated with differential expression analysis for genes and pathways. Marked changes were found in pathways concerned with endoplasmic reticulum stress with its unfolded protein response (UPR), apoptotic pathways and components of inflammation. In addition, there were changes in genes important for islet cell identity. These findings advance our understanding of why islets are stressed before transplantation, which may lead to strategies to reduce this stress and lead to better clinical outcomes.

Adipose-derived circulating miRNAs regulate gene expression in other tissues

Adipose tissue is a major site of energy storage and has a role in the regulation of metabolism through the release of adipokines. Here we show that mice with an adipose-tissue-specific knockout of the microRNA (miRNA)-processing enzyme Dicer (ADicerKO), as well as humans with lipodystrophy, exhibit a substantial decrease in levels of circulating exosomal miRNAs. Transplantation of both white and brown adipose tissue-brown especially-into ADicerKO mice restores the level of numerous circulating miRNAs that are associated with an improvement in glucose tolerance and a reduction in hepatic Fgf21 mRNA and circulating FGF21. This gene regulation can be mimicked by the administration of normal, but not ADicerKO, serum exosomes. Expression of a human-specific miRNA in the brown adipose tissue of one mouse in vivo can also regulate its 3' UTR reporter in the liver of another mouse through serum exosomal transfer. Thus, adipose tissue constitutes an important source of circulating exosomal miRNAs, which can regulate gene expression in distant tissues and thereby serve as a previously undescribed form of adipokine.

Antibiotic effects on gut microbiota and metabolism are host dependent

Interactions of diet, gut microbiota, and host genetics play important roles in the development of obesity and insulin resistance. Here, we have investigated the molecular links between gut microbiota, insulin resistance, and glucose metabolism in 3 inbred mouse strains with differing susceptibilities to metabolic syndrome using diet and antibiotic treatment. Antibiotic treatment altered intestinal microbiota, decreased tissue inflammation, improved insulin signaling in basal and stimulated states, and improved glucose metabolism in obesity- and diabetes-prone C57BL/6J mice on a high-fat diet (HFD). Many of these changes were reproduced by the transfer of gut microbiota from antibiotic-treated donors to germ-free or germ-depleted mice. These physiological changes closely correlated with changes in serum bile acids and levels of the antiinflammatory bile acid receptor Takeda G protein-coupled receptor 5 (TGR5) and were partially recapitulated by treatment with a TGR5 agonist. In contrast, antibiotic treatment of HFD-fed, obesity-resistant 129S1 and obesity-prone 129S6 mice did not improve metabolism, despite changes in microbiota and bile acids. These mice also failed to show a reduction in inflammatory gene expression in response to the TGR5 agonist. Thus, changes in bile acid and inflammatory signaling, insulin resistance, and glucose metabolism driven by an HFD can be modified by antibiotic-induced changes in gut microbiota; however, these effects depend on important interactions with the host's genetic background and inflammatory potential.

Defects in muscle branched-chain amino acid oxidation contribute to impaired lipid metabolism

OBJECTIVE

Plasma levels of branched-chain amino acids (BCAA) are consistently elevated in obesity and type 2 diabetes (T2D) and can also prospectively predict T2D. However, the role of BCAA in the pathogenesis of insulin resistance and T2D remains unclear.

METHODS

To identify pathways related to insulin resistance, we performed comprehensive gene expression and metabolomics analyses in skeletal muscle from 41 humans with normal glucose tolerance and 11 with T2D across a range of insulin sensitivity (SI, 0.49 to 14.28). We studied both cultured cells and mice heterozygous for the BCAA enzyme methylmalonyl-CoA mutase (Mut) and assessed the effects of altered BCAA flux on lipid and glucose homeostasis.

RESULTS

Our data demonstrate perturbed BCAA metabolism and fatty acid oxidation in muscle from insulin resistant humans. Experimental alterations in BCAA flux in cultured cells similarly modulate fatty acid oxidation. Mut heterozygosity in mice alters muscle lipid metabolism in vivo, resulting in increased muscle triglyceride accumulation, increased plasma glucose, hyperinsulinemia, and increased body weight after high-fat feeding.

CONCLUSIONS

Our data indicate that impaired muscle BCAA catabolism may contribute to the development of insulin resistance by perturbing both amino acid and fatty acid metabolism and suggest that targeting BCAA metabolism may hold promise for prevention or treatment of T2D.

Maternal obesity programs mitochondrial and lipid metabolism gene expression in infant umbilical vein endothelial cells

BACKGROUND/OBJECTIVES

Maternal obesity increases risk for childhood obesity, but molecular mechanisms are not well understood. We hypothesized that primary umbilical vein endothelial cells (HUVEC) from infants of overweight and obese mothers would harbor transcriptional patterns reflecting offspring obesity risk.

SUBJECTS/METHODS

In this observational cohort study, we recruited 13 lean (pre-pregnancy body mass index (BMI) <25.0 kg m^-2) and 24 overweight-obese ('ov-ob', BMI⩾25.0 kg m^-2) women. We isolated primary HUVEC, and analyzed both gene expression (Primeview, Affymetrix) and cord blood levels of hormones and adipokines.

RESULTS

A total of 142 transcripts were differentially expressed in HUVEC from infants of overweight-obese mothers (false discovery rate, FDR<0.05). Pathway analysis revealed that genes involved in mitochondrial and lipid metabolism were negatively correlated with maternal BMI (FDR<0.05). To test whether these transcriptomic patterns were associated with distinct nutrient exposures in the setting of maternal obesity, we analyzed the cord blood lipidome and noted significant increases in the levels of total free fatty acids (lean: 95.5±37.1 μg ml^-1, ov-ob: 124.1±46.0 μg ml^-1, P=0.049), palmitate (lean: 34.5±12.7 μg ml^-1, ov-ob: 46.3±18.4 μg ml^-1, P=0.03) and stearate (lean: 20.8±8.2 μg ml^-1, ov-ob: 29.7±17.2 μg ml^-1, P=0.04), in infants of overweight-obese mothers.

CONCLUSIONS

Prenatal exposure to maternal obesity alters HUVEC expression of genes involved in mitochondrial and lipid metabolism, potentially reflecting developmentally programmed differences in oxidative and lipid metabolism.

Dietary Betaine Supplementation Increases Fgf21 Levels to Improve Glucose Homeostasis and Reduce Hepatic Lipid Accumulation in Mice

Identifying markers of human insulin resistance may permit development of new approaches for treatment and prevention of type 2 diabetes. To this end, we analyzed the fasting plasma metabolome in metabolically characterized human volunteers across a spectrum of insulin resistance. We demonstrate that plasma betaine levels are reduced in insulin-resistant humans and correlate closely with insulin sensitivity. Moreover, betaine administration to mice with diet-induced obesity prevents the development of impaired glucose homeostasis, reduces hepatic lipid accumulation, increases white adipose oxidative capacity, and enhances whole-body energy expenditure. In parallel with these beneficial metabolic effects, betaine supplementation robustly increased hepatic and circulating fibroblast growth factor (Fgf)21 levels. Betaine administration failed to improve glucose homeostasis and liver fat content in Fgf21(-/-) mice, demonstrating that Fgf21 is necessary for betaine's beneficial effects. Together, these data indicate that dietary betaine increases Fgf21 levels to improve metabolic health in mice and suggest that betaine supplementation merits further investigation as a supplement for treatment or prevention of type 2 diabetes in humans.

Metabolic modeling predicts metabolite changes in Mycobacterium tuberculosis

BACKGROUND

Mycobacterium tuberculosis (MTB) is the causal agent of the disease tuberculosis (TB). Metabolic adaptations are thought to be critical to the survival of MTB during pathogenesis. Computational tools that can be used to study MTB metabolism in silico and prioritize resource-intensive experimental work could significantly accelerate research.

RESULTS

We have developed E-Flux-MFC, an enhancement of our original E-Flux method that enables the prediction of changes in the production of external and internal metabolites corresponding to gene expression measurements. We have used this method to simulate the changes in the metabolic state of Mycobacterium tuberculosis (MTB). We have validated the accuracy of E-Flux-MFC for predicting changes in lipids and metabolites during a hypoxia time course using previously published metabolomics and transcriptomics data. We have further validated the accuracy of the method for predicting changes in MTB lipids following the deletion and induction of two well-studied transcription factors (TFs). We have applied the method to predict the metabolic impact of the induction of each of the approximately 180 MTB TFs using a previously generated and publically available expression data set.

CONCLUSIONS

E-flux-MFC can be used to study global changes in MTB metabolites from gene expression data associated with environmental and genetic perturbations. The application of this method to a data set of MTB TF perturbations provides a resource for studying the large number of TFs whose functions remain unknown. Most TFs impact metabolites indirectly through the propagation of gene expression changes through the regulatory network rather than through their direct regulons. E-Flux-MFC is also applicable to any organism for which accurate metabolic models are available.

Lipid-mediated regulation of SKN-1/Nrf in response to germ cell absence

In Caenorhabditis elegans, ablation of germline stem cells (GSCs) extends lifespan, but also increases fat accumulation and alters lipid metabolism, raising the intriguing question of how these effects might be related. Here, we show that a lack of GSCs results in a broad transcriptional reprogramming in which the conserved detoxification regulator SKN-1/Nrf increases stress resistance, proteasome activity, and longevity. SKN-1 also activates diverse lipid metabolism genes and reduces fat storage, thereby alleviating the increased fat accumulation caused by GSC absence. Surprisingly, SKN-1 is activated by signals from this fat, which appears to derive from unconsumed yolk that was produced for reproduction. We conclude that SKN-1 plays a direct role in maintaining lipid homeostasis in which it is activated by lipids. This SKN-1 function may explain the importance of mammalian Nrf proteins in fatty liver disease and suggest that particular endogenous or dietary lipids might promote health through SKN-1/Nrf.

Metabolic modeling of muscle metabolism identifies key reactions linked to insulin resistance phenotypes

Objective

Dysregulated muscle metabolism is a cardinal feature of human insulin resistance (IR) and associated diseases, including type 2 diabetes (T2D). However, specific reactions contributing to abnormal energetics and metabolic inflexibility in IR are unknown.

Methods

We utilize flux balance computational modeling to develop the first systems-level analysis of IR metabolism in fasted and fed states, and varying nutrient conditions. We systematically perturb the metabolic network to identify reactions that reproduce key features of IR-linked metabolism.

Results

While reduced glucose uptake is a major hallmark of IR, model-based reductions in either extracellular glucose availability or uptake do not alter metabolic flexibility, and thus are not sufficient to fully recapitulate IR-linked metabolism. Moreover, experimentally-reduced flux through single reactions does not reproduce key features of IR-linked metabolism. However, dual knockdowns of pyruvate dehydrogenase (PDH), in combination with reduced lipid uptake or lipid/amino acid oxidation (ETFDH), does reduce ATP synthesis, TCA cycle flux, and metabolic flexibility. Experimental validation demonstrates robust impact of dual knockdowns in PDH/ETFDH on cellular energetics and TCA cycle flux in cultured myocytes. Parallel analysis of transcriptomic and metabolomics data in humans with IR and T2D demonstrates downregulation of PDH subunits and upregulation of its inhibitory kinase PDK4, both of which would be predicted to decrease PDH flux, concordant with the model.

Conclusions

Our results indicate that complex interactions between multiple biochemical reactions contribute to metabolic perturbations observed in human IR, and that the PDH complex plays a key role in these metabolic phenotypes.

Reconstruction and validation of a genome-scale metabolic model for the filamentous fungus Neurospora crassa using FARM

The filamentous fungus Neurospora crassa played a central role in the development of twentieth-century genetics, biochemistry and molecular biology, and continues to serve as a model organism for eukaryotic biology. Here, we have reconstructed a genome-scale model of its metabolism. This model consists of 836 metabolic genes, 257 pathways, 6 cellular compartments, and is supported by extensive manual curation of 491 literature citations. To aid our reconstruction, we developed three optimization-based algorithms, which together comprise Fast Automated Reconstruction of Metabolism (FARM). These algorithms are: LInear MEtabolite Dilution Flux Balance Analysis (limed-FBA), which predicts flux while linearly accounting for metabolite dilution; One-step functional Pruning (OnePrune), which removes blocked reactions with a single compact linear program; and Consistent Reproduction Of growth/no-growth Phenotype (CROP), which reconciles differences between in silico and experimental gene essentiality faster than previous approaches. Against an independent test set of more than 300 essential/non-essential genes that were not used to train the model, the model displays 93% sensitivity and specificity. We also used the model to simulate the biochemical genetics experiments originally performed on Neurospora by comprehensively predicting nutrient rescue of essential genes and synthetic lethal interactions, and we provide detailed pathway-based mechanistic explanations of our predictions. Our model provides a reliable computational framework for the integration and interpretation of ongoing experimental efforts in Neurospora, and we anticipate that our methods will substantially reduce the manual effort required to develop high-quality genome-scale metabolic models for other organisms.

Few organisms have been as foundational to the development of modern genetics and cellular metabolism as Neurospora crassa. Given the wealth of knowledge available for this filamentous fungus, the effort required to manually curate a high-quality genome-scale metabolic reconstruction would be daunting. To aid the reconstruction process, we developed three optimization-based algorithms. The first algorithm predicts flux while linearly accounting for metabolite dilution; the second algorithm removes blocked reactions with one compact linear program; and the third algorithm reconciles differences between in silico predictions and experimental observations of mutant viability. We have used these algorithms to develop the first genome-scale metabolic model for Neurospora. We have validated the accuracy of our model against an independent test set of more than 300 growth/no-growth phenotypes, and our model displays 93% sensitivity and specificity. Simulating the biochemical genetics experiments originally performed on Neurospora, we comprehensively predicted essential genes, nutrient rescues of auxotroph mutants and synthetic lethal interactions. With these predictions, we provide potential mechanistic insight into known mutant phenotypes, and testable hypotheses for novel mutant phenotypes. The model, the algorithms and the testable hypotheses provide a computational foundation for the study of Neurospora crassa metabolism.

How accurate can genetic predictions be?

Background

Pre-symptomatic prediction of disease and drug response based on genetic testing is a critical component of personalized medicine. Previous work has demonstrated that the predictive capacity of genetic testing is constrained by the heritability and prevalence of the tested trait, although these constraints have only been approximated under the assumption of a normally distributed genetic risk distribution.

Results

Here, we mathematically derive the absolute limits that these factors impose on test accuracy in the absence of any distributional assumptions on risk. We present these limits in terms of the best-case receiver-operating characteristic (ROC) curve, consisting of the best-case test sensitivities and specificities, and the AUC (area under the curve) measure of accuracy. We apply our method to genetic prediction of type 2 diabetes and breast cancer, and we additionally show the best possible accuracy that can be obtained from integrated predictors, which can incorporate non-genetic features.

Conclusion

Knowledge of such limits is valuable in understanding the implications of genetic testing even before additional associations are identified.

Expression profiling of uterine leiomyomata cytogenetic subgroups reveals distinct signatures in matched myometrium: transcriptional profilingof the t(12;14) and evidence in support of predisposing genetic heterogeneity

Uterine leiomyomata (UL), the most common neoplasm in reproductive-age women, are classified into distinct genetic subgroups based on recurrent chromosome abnormalities. To develop a molecular signature of UL with t(12;14)(q14-q15;q23-q24), we took advantage of the multiple UL arising as independent clonal lesions within a single uterus. We compared genome-wide expression levels of t(12;14) UL to non-t(12;14) UL from each of nine women in a paired analysis, with each sample weighted for the percentage of t(12;14) cells to adjust for mosaicism with normal cells. This resulted in a transcriptional profile that confirmed HMGA2, known to be overexpressed in t(12;14) UL, as the most significantly altered gene. Pathway analysis of the differentially expressed genes showed significant association with cell proliferation, particularly G1/S checkpoint regulation. This is consistent with the known larger size of t(12;14) UL relative to karyotypically normal UL or to UL in the deletion 7q22 subgroup. Unsupervised hierarchical clustering demonstrated that patient variability is relatively dominant to the distinction of t(12;14) UL compared with non-t(12;14) UL or of t(12;14) UL compared with del(7q) UL. The paired design we employed is therefore important to produce an accurate t(12;14) UL-specific gene list by removing the confounding effects of genotype and environment. Interestingly, myometrium not only clustered away from the tumors, but generally separated based on associated t(12;14) versus del(7q) status. Nine genes were identified whose expression can distinguish the myometrium origin. This suggests an underlying constitutional genetic predisposition to these somatic changes which could potentially lead to improved personalized management and treatment.

Upregulation of proapoptotic microRNA mir-125a after massive small bowel resection in rats

OBJECTIVE

Short bowel syndrome remains a condition of high morbidity and mortality, and current therapeutic options carry significant side effects. To identify new treatments we focused on postresection changes in microRNAs--short noncoding RNAs, which suppress target genes--and suggest a previously undiscovered role for microRNA-125a (mir-125a) in intestinal adaptation.

METHODS

Rats underwent either 80% massive small bowel resection or transection and were harvested after 48 hours. Jejunum was harvested for microRNA microarrays, laser capture microdissection, and RNA and protein analysis. Mir-125a was overexpressed in intestinal epithelium-6 (crypt-derived) cells (IEC-6) and effects on proliferation and apoptosis determined using MTS and flow cytometry. Expression of potential targets of mir-125a in rat jejunum and IEC-6 cells was determined using quantitative real-time polymerase chain reaction (RNA) and Western blotting (protein).

RESULTS

Resection upregulated mir-125a and mir-214 by 2.4-folds and 3.2-folds, respectively. Highest levels of expression were noted in the crypt fraction. Mir-125a overexpression induced apoptosis and resultant growth arrest in IEC-6 cells. The expression of the prosurvival Bcl-2 family member Mcl-1 was downregulated in both mir-125a-overexpressing IEC-6 cells and in jejunum of resected rats, confirming Mcl-1 as a previously undiscovered target of mir-125a.

CONCLUSIONS

Upregulation of mir-125a suppresses the prosurvival protein Mcl1, producing the increase in apoptosis known to accompany the proliferative changes characteristic of intestinal adaptation. Our data highlight a potential role for microRNAs as mediators of the adaptive process and may facilitate the development of new therapeutic options for short bowel syndrome.

Glycosaminoglycans and glucose prevent apoptosis in 4-methylumbelliferone-treated human aortic smooth muscle cells

Smooth muscle cells (SMCs) have a pivotal role in cardiovascular diseases and are responsible for hyaluronan (HA) deposition in thickening vessel walls. HA regulates SMC proliferation, migration, and inflammation, which accelerates neointima formation. We used the HA synthesis inhibitor 4-methylumbelliferone (4-MU) to reduce HA production in human aortic SMCs and found a significant increase of apoptotic cells. Interestingly, the exogenous addition of HA together with 4-MU reduced apoptosis. A similar anti-apoptotic effect was observed also by adding other glycosaminoglycans and glucose to 4-MU-treated cells. Furthermore, the anti-apoptotic effect of HA was mediated by Toll-like receptor 4, CD44, and PI3K but not by ERK1/2.

A transcriptional network signature characterizes lung cancer subtypes

BACKGROUND

Transcriptional networks play a central role in cancer development. The authors described a systems biology approach to cancer classification based on the reverse engineering of the transcriptional network surrounding the 2 most common types of lung cancer: adenocarcinoma (AC) and squamous cell carcinoma (SCC).

METHODS

A transcriptional network classifier was inferred from the molecular profiles of 111 human lung carcinomas. The authors tested its classification accuracy in 7 independent cohorts, for a total of 422 subjects of Caucasian, African, and Asian descent.

RESULTS

The model for distinguishing AC from SCC was a 25-gene network signature. Its performance on the 7 independent cohorts achieved 95.2% classification accuracy. Even more surprisingly, 95% of this accuracy was explained by the interplay of 3 genes (KRT6A, KRT6B, KRT6C) on a narrow cytoband of chromosome 12. The role of this chromosomal region in distinguishing AC and SCC was further confirmed by the analysis of another group of 28 independent subjects assayed by DNA copy number changes. The copy number variations of bands 12q12, 12q13, and 12q12-13 discriminated these samples with 84% accuracy.

CONCLUSIONS

These results suggest the existence of a robust signature localized in a relatively small area of the genome, and show the clinical potential of reverse engineering transcriptional networks from molecular profiles.

MicroRNA mir-16 is anti-proliferative in enterocytes and exhibits diurnal rhythmicity in intestinal crypts

BACKGROUND AND AIMS

The intestine exhibits profound diurnal rhythms in function and morphology, in part due to changes in enterocyte proliferation. The regulatory mechanisms behind these rhythms remain largely unknown. We hypothesized that microRNAs are involved in mediating these rhythms, and studied the role of microRNAs specifically in modulating intestinal proliferation.

METHODS

Diurnal rhythmicity of microRNAs in rat jejunum was analyzed by microarrays and validated by qPCR. Temporal expression of diurnally rhythmic mir-16 was further quantified in intestinal crypts, villi, and smooth muscle using laser capture microdissection and qPCR. Morphological changes in rat jejunum were assessed by histology and proliferation by immunostaining for bromodeoxyuridine. In IEC-6 cells stably overexpressing mir-16, proliferation was assessed by cell counting and MTS assay, cell cycle progression and apoptosis by flow cytometry, and cell cycle gene expression by qPCR and immunoblotting.

RESULTS

mir-16 peaked 6 hours after light onset (HALO 6) with diurnal changes restricted to crypts. Crypt depth and villus height peaked at HALO 13-14 in antiphase to mir-16. Overexpression of mir-16 in IEC-6 cells suppressed specific G1/S regulators (cyclins D1-3, cyclin E1 and cyclin-dependent kinase 6) and produced G1 arrest. Protein expression of these genes exhibited diurnal rhythmicity in rat jejunum, peaking between HALO 11 and 17 in antiphase to mir-16.

CONCLUSIONS

This is the first report of circadian rhythmicity of specific microRNAs in rat jejunum. Our data provide a link between anti-proliferative mir-16 and the intestinal proliferation rhythm and point to mir-16 as an important regulator of proliferation in jejunal crypts. This function may be essential to match proliferation and absorptive capacity with nutrient availability.

Integrative predictive model of coronary artery calcification in atherosclerosis

BACKGROUND

Many different genetic and clinical factors have been identified as causes or contributors to atherosclerosis. We present a model of preclinical atherosclerosis based on genetic and clinical data that predicts the presence of coronary artery calcification in healthy Americans of European descent 45 to 84 years of age in the Multi-Ethnic Study of Atherosclerosis (MESA).

METHODS AND RESULTS

We assessed 712 individuals for the presence or absence of coronary artery calcification and assessed their genotypes for 2882 single-nucleotide polymorphisms. With the use of these single-nucleotide polymorphisms and relevant clinical data, a Bayesian network that predicts the presence of coronary calcification was constructed. The model contained 13 single-nucleotide polymorphisms (from genes AGTR1, ALOX15, INSR, PRKAB1, IL1R2, ESR2, KCNK1, FBLN5, PPARA, VEGFA, PON1, TDRD6, PLA2G7, and 1 ancestry informative marker) and 5 clinical variables (sex, age, weight, smoking, and diabetes mellitus) and achieved 85% predictive accuracy, as measured by area under the receiver operating characteristic curve. This is a significant (P<0.001) improvement on models that use just the single-nucleotide polymorphism data or just the clinical variables.

CONCLUSIONS

We present an investigation of joint genetic and clinical factors associated with atherosclerosis that shows predictive results for both cases, as well as enhanced performance for their combination.

Identifying the molecular signature of the interstitial deletion 7q subgroup of uterine leiomyomata using a paired analysis

Uterine leiomyomata (UL), the most common neoplasm in reproductive-age women, have recurrent cytogenetic abnormalities including interstitial deletion of 7q. To develop a molecular signature, matched del(7q) and non-del(7q) tumors identified by FISH or karyotyping from 11 women were profiled with expression arrays. Our analysis using paired t tests demonstrates this matched design is critical to eliminate the confounding effects of genotype and environment that underlie patient variation. A gene list ordered by genome-wide significance showed enrichment for the 7q22 target region. Modification of the gene list by weighting each sample for percent of del(7q) cells to account for the mosaic nature of these tumors further enhanced the frequency of 7q22 genes. Pathway analysis revealed two of the 19 significant functional networks were associated with development and the most represented pathway was protein ubiquitination, which can influence tumor development by stabilizing oncoproteins and destabilizing tumor suppressor proteins. Array CGH (aCGH) studies determined the only consistent genomic imbalance was deletion of 9.5 megabases from 7q22-7q31.1. Combining the aCGH data with the del(7q) UL mosaicism-weighted expression analysis resulted in a list of genes that are commonly deleted and whose copy number is correlated with significantly decreased expression. These genes include the proliferation inhibitor HPB1, the loss of expression of which has been associated with invasive breast cancer, as well as the mitosis integrity-maintenance tumor suppressor RINT1. This study provides a molecular signature of the del(7q) UL subgroup and will serve as a platform for future studies of tumor pathogenesis.