Systems genetics approach uncovers associations between host amylase locus, gut microbiome and metabolic traits in hyperlipidemic mice

The molecular basis for how host genetic variation impacts gut microbial community and bacterial metabolic niches remain largely unknown. We leveraged 90 inbred hyperlipidemic mouse strains from the Hybrid Mouse Diversity Panel (HMDP), previously studied for a variety of cardio-metabolic traits. Metagenomic analysis of cecal DNA followed by genome-wide association analysis identified genomic loci that were associated with microbial enterotypes in the gut. Among these we detected a genetic locus surrounding multiple amylase genes that was associated with abundances of Firmicutes (Lachnospiraceae family) and Bacteroidetes (Muribaculaceae family) taxa encoding distinct starch and sugar metabolism functions. We also found that lower amylase gene number in the mouse genome was associated with higher gut Muribaculaceae levels. Previous work suggests that modulation of host amylase activity impacts the availability of carbohydrates to the host and potentially to gut bacteria. The genetic variants described above were associated with distinct gut microbial communities (enterotypes) with different predicted metabolic capacities for carbohydrate degradation. Mendelian randomization analysis revealed host phenotypes, including liver fibrosis and plasma HDL-cholesterol levels, that were associated with gut microbiome enterotypes. This work reveals novel relationships between host genetic variation, gut microbial enterotypes and host physiology/disease phenotypes in mice.

Genetic Analysis of Obesity-Induced Diabetic Nephropathy in BTBR Mice

Diabetic nephropathy (DN) is the leading cause of end-stage renal disease in the U.S. and has a significant impact on human suffering. Leptin-deficient BTBR (BTBRob/ob) mice develop hallmark features of obesity-induced DN, whereas leptin-deficient C57BL/6J (B6ob/ob) mice do not. To identify genetic loci that underlie this strain difference, we constructed an F2 intercross between BTBRob/ob and B6ob/ob mice. We isolated kidneys from 460 F2 mice and histologically scored them for percent mesangial matrix and glomerular volume (∼50 glomeruli per mouse), yielding ∼45,000 distinct measures in total. The same histological measurements were made in kidneys from B6 and BTBR mice, either lean or obese (Lepob/ob), at 4 and 10 weeks of age, allowing us to assess the contribution of strain, age, and obesity to glomerular pathology. All F2 mice were genotyped for ∼5,000 single nucleotide polymorphisms (SNPs), ∼2,000 of which were polymorphic between B6 and BTBR, enabling us to identify a quantitative trait locus (QTL) on chromosome 7, with a peak at ∼30 Mbp, for percent mesangial matrix, glomerular volume, and mesangial volume. The podocyte-specific gene nephrin (Nphs1) is physically located at the QTL and contains high-impact SNPs in BTBR, including several missense variants within the extracellular immunoglobulin-like domains.

ARTICLE HIGHLIGHTS

Erratum: Proteomic pathways to metabolic disease and type 2 diabetes in the pancreatic islet

[This corrects the article DOI: 10.1016/j.isci.2021.103099.].

D2H2: diabetes data and hypothesis hub

Motivation

There is a rapid growth in the production of omics datasets collected by the diabetes research community. However, such published data are underutilized for knowledge discovery. To make bioinformatics tools and published omics datasets from the diabetes field more accessible to biomedical researchers, we developed the Diabetes Data and Hypothesis Hub (D2H2).

Results

D2H2 contains hundreds of high-quality curated transcriptomics datasets relevant to diabetes, accessible via a user-friendly web-based portal. The collected and processed datasets are curated from the Gene Expression Omnibus (GEO). Each curated study has a dedicated page that provides data visualization, differential gene expression analysis, and single-gene queries. To enable the investigation of these curated datasets and to provide easy access to bioinformatics tools that serve gene and gene set-related knowledge, we developed the D2H2 chatbot. Utilizing GPT, we prompt users to enter free text about their data analysis needs. Parsing the user prompt, together with specifying information about all D2H2 available tools and workflows, we answer user queries by invoking the most relevant tools via the tools' API. D2H2 also has a hypotheses generation module where gene sets are randomly selected from the bulk RNA-seq precomputed signatures. We then find highly overlapping gene sets extracted from publications listed in PubMed Central with abstract dissimilarity. With the help of GPT, we speculate about a possible explanation of the high overlap between the gene sets. Overall, D2H2 is a platform that provides a suite of bioinformatics tools and curated transcriptomics datasets for hypothesis generation.

Availability and implementation

D2H2 is available at: https://d2h2.maayanlab.cloud/ and the source code is available from GitHub at https://github.com/MaayanLab/D2H2-site under the CC BY-NC 4.0 license.

Identification of genetic drivers of plasma lipoprotein size in the Diversity Outbred mouse population

Despite great progress in understanding lipoprotein physiology, there is still much to be learned about the genetic drivers of lipoprotein abundance, composition, and function. We used ion mobility spectrometry to survey 16 plasma lipoprotein subfractions in 500 Diversity Outbred mice maintained on a Western-style diet. We identified 21 quantitative trait loci (QTL) affecting lipoprotein abundance. To refine the QTL and link them to disease risk in humans, we asked if the human homologs of genes located at each QTL were associated with lipid traits in human genome-wide association studies. Integration of mouse QTL with human genome-wide association studies yielded candidate gene drivers for 18 of the 21 QTL. This approach enabled us to nominate the gene encoding the neutral ceramidase, Asah2, as a novel candidate driver at a QTL on chromosome 19 for large HDL particles (HDL-2b). To experimentally validate Asah2, we surveyed lipoproteins in Asah2-/- mice. Compared to wild-type mice, female Asah2-/- mice showed an increase in several lipoproteins, including HDL. Our results provide insights into the genetic regulation of circulating lipoproteins, as well as mechanisms by which lipoprotein subfractions may affect cardiovascular disease risk in humans.

Missense variants in SORT1 are associated with LDL-C in an Amish population

Common noncoding variants at the human 1p13.3 locus associated with SORT1 expression are among those most strongly associated with low-density lipoprotein cholesterol (LDL-C) in human genome-wide association studies. However, validation studies in mice and cell lines have produced variable results regarding the directionality of the effect of SORT1 on LDL-C. This, together with the fact that the 1p13.3 variants are associated with expression of several genes, has raised the question of whether SORT1 is the causal gene at this locus. Using whole exome sequencing in members of an Amish population, we identified coding variants in SORT1 that are associated with increased (rs141749679, K302E) and decreased (rs149456022, Q225H) LDL-C. Further, analysis of plasma lipoprotein particle subclasses by ion mobility in a subset of rs141749679 (K302E) carriers revealed higher levels of large LDL particles compared to noncarriers. In contrast to the effect of these variants in the Amish, the sortilin K302E mutation introduced into a C57BL/6J mouse via CRISPR/Cas9 resulted in decreased non-high-density lipoprotein cholesterol, and the sortilin Q225H mutation did not alter cholesterol levels in mice. This is indicative of different effects of these mutations on cholesterol metabolism in the two species. To our knowledge, this is the first evidence that naturally occurring coding variants in SORT1 are associated with LDL-C, thus supporting SORT1 as the gene responsible for the association of the 1p13.3 locus with LDL-C.

Heme biosynthesis regulates BCAA catabolism and thermogenesis in brown adipose tissue

With age, people tend to accumulate body fat and reduce energy expenditure 1 . Brown (BAT) and beige adipose tissue dissipate heat and increase energy expenditure via the activity of the uncoupling protein UCP1 and other thermogenic futile cycles 2,3 . The activity of brown and beige depots inversely correlates with BMI and age 4-11 , suggesting that promoting thermogenesis may be an effective approach for combating age-related metabolic disease 12-15 . Heme is an enzyme cofactor and signaling molecule that we recently showed to regulate BAT function 16 . Here, we show that heme biosynthesis is the primary contributor to intracellular heme levels in brown adipocytes. Inhibition of heme biosynthesis leads to mitochondrial dysfunction and reduction in UCP1. Although supplementing heme can restore mitochondrial function in heme-synthesis-deficient cells, the downregulation of UCP1 persists due to the accumulation of the heme precursors, particularly propionyl-CoA, which is a product of branched-chain amino acids (BCAA) catabolism. Cold exposure promotes BCAA uptake in BAT, and defects in BCAA catabolism in this tissue hinder thermogenesis 17 . However, BCAAs' contribution to the TCA cycle in BAT and WAT never exceeds 2% of total TCA flux 18 . Our work offers a way to integrate current literature by describing heme biosynthesis as an important metabolic sink for BCAAs.

An Enhancer Within Abcb11 Regulates G6pc2 in C57BL/6 Mouse Pancreatic Islets

G6PC2 is predominantly expressed in pancreatic islet β-cells where it encodes a glucose-6-phosphatase catalytic subunit that modulates the sensitivity of insulin secretion to glucose by opposing the action of glucokinase, thereby regulating fasting blood glucose (FBG). Prior studies have shown that the G6pc2 promoter alone is unable to confer sustained islet-specific gene expression in mice, suggesting the existence of distal enhancers that regulate G6pc2 expression. Using information from both mice and humans and knowledge that single nucleotide polymorphisms (SNPs) both within and near G6PC2 are associated with variations in FBG in humans, we identified several putative enhancers 3' of G6pc2. One region, herein referred to as enhancer I, resides in the 25th intron of Abcb11 and binds multiple islet-enriched transcription factors. CRISPR-mediated deletion of enhancer I in C57BL/6 mice had selective effects on the expression of genes near the G6pc2 locus. In isolated islets, G6pc2 and Spc25 expression were reduced ∼50%, and Gm13613 expression was abolished, whereas Cers6 and nostrin expression were unaffected. This partial reduction in G6pc2 expression enhanced islet insulin secretion at basal glucose concentrations but did not affect FBG or glucose tolerance in vivo, consistent with the absence of a phenotype in G6pc2 heterozygous C57BL/6 mice.

ARTICLE HIGHLIGHTS

PPTC7 maintains mitochondrial protein content by suppressing receptor-mediated mitophagy

PPTC7 is a resident mitochondrial phosphatase essential for maintaining proper mitochondrial content and function. Newborn mice lacking Pptc7 exhibit aberrant mitochondrial protein phosphorylation, suffer from a range of metabolic defects, and fail to survive beyond one day after birth. Using an inducible knockout model, we reveal that loss of Pptc7 in adult mice causes marked reduction in mitochondrial mass and metabolic capacity with elevated hepatic triglyceride accumulation. Pptc7 knockout animals exhibit increased expression of the mitophagy receptors BNIP3 and NIX, and Pptc7-/- mouse embryonic fibroblasts (MEFs) display a major increase in mitophagy that is reversed upon deletion of these receptors. Our phosphoproteomics analyses reveal a common set of elevated phosphosites between perinatal tissues, adult liver, and MEFs, including multiple sites on BNIP3 and NIX, and our molecular studies demonstrate that PPTC7 can directly interact with and dephosphorylate these proteins. These data suggest that Pptc7 deletion causes mitochondrial dysfunction via dysregulation of several metabolic pathways and that PPTC7 may directly regulate mitophagy receptor function or stability. Overall, our work reveals a significant role for PPTC7 in the mitophagic response and furthers the growing notion that management of mitochondrial protein phosphorylation is essential for ensuring proper organelle content and function.

Novel regulators of islet function identified from genetic variation in mouse islet Ca2+ oscillations

Insufficient insulin secretion to meet metabolic demand results in diabetes. The intracellular flux of Ca2+ into β-cells triggers insulin release. Since genetics strongly influences variation in islet secretory responses, we surveyed islet Ca2+ dynamics in eight genetically diverse mouse strains. We found high strain variation in response to four conditions: (1) 8 mM glucose; (2) 8 mM glucose plus amino acids; (3) 8 mM glucose, amino acids, plus 10 nM glucose-dependent insulinotropic polypeptide (GIP); and (4) 2 mM glucose. These stimuli interrogate β-cell function, α- to β-cell signaling, and incretin responses. We then correlated components of the Ca2+ waveforms to islet protein abundances in the same strains used for the Ca2+ measurements. To focus on proteins relevant to human islet function, we identified human orthologues of correlated mouse proteins that are proximal to glycemic-associated single-nucleotide polymorphisms in human genome-wide association studies. Several orthologues have previously been shown to regulate insulin secretion (e.g. ABCC8, PCSK1, and GCK), supporting our mouse-to-human integration as a discovery platform. By integrating these data, we nominate novel regulators of islet Ca2+ oscillations and insulin secretion with potential relevance for human islet function. We also provide a resource for identifying appropriate mouse strains in which to study these regulators.

Pancreatic islet α cell function and proliferation requires the arginine transporter SLC7A2

Interrupting glucagon signaling decreases gluconeogenesis and the fractional extraction of amino acids by liver from blood resulting in lower glycemia. The resulting hyperaminoacidemia stimulates α cell proliferation and glucagon secretion via a liver-α cell axis. We hypothesized that α cells detect and respond to circulating amino acids levels via a unique amino acid transporter repertoire. We found that Slc7a2ISLC7A2 is the most highly expressed cationic amino acid transporter in α cells with its expression being three-fold greater in α than β cells in both mouse and human. Employing cell culture, zebrafish, and knockout mouse models, we found that the cationic amino acid arginine and SLC7A2 are required for α cell proliferation in response to interrupted glucagon signaling. Ex vivo and in vivo assessment of islet function in Slc7a2-/- mice showed decreased arginine-stimulated glucagon and insulin secretion. We found that arginine activation of mTOR signaling and induction of the glutamine transporter SLC38A5 was dependent on SLC7A2, showing that both's role in α cell proliferation is dependent on arginine transport and SLC7A2. Finally, we identified single nucleotide polymorphisms in SLC7A2 associated with HbA1c. Together, these data indicate a central role for SLC7A2 in amino acid-stimulated α cell proliferation and islet hormone secretion.

What the BTBR/J mouse has taught us about diabetes and diabetic complications

Human and mouse genetics have delivered numerous diabetogenic loci, but it is mainly through the use of animal models that the pathophysiological basis for their contribution to diabetes has been investigated. More than 20 years ago, we serendipidously identified a mouse strain that could serve as a model of obesity-prone type 2 diabetes, the BTBR (Black and Tan Brachyury) mouse (BTBR T+ Itpr3tf/J, 2018) carrying the Lepob mutation. We went on to discover that the BTBR-Lepob mouse is an excellent model of diabetic nephropathy and is now widely used by nephrologists in academia and the pharmaceutical industry. In this review, we describe the motivation for developing this animal model, the many genes identified and the insights about diabetes and diabetes complications derived from >100 studies conducted in this remarkable animal model.

Lipidomic QTL in Diversity Outbred mice identifies a novel function for α/β hydrolase domain 2 (Abhd2) as an enzyme that metabolizes phosphatidylcholine and cardiolipin

We and others have previously shown that genetic association can be used to make causal connections between gene loci and small molecules measured by mass spectrometry in the bloodstream and in tissues. We identified a locus on mouse chromosome 7 where several phospholipids in liver showed strong genetic association to distinct gene loci. In this study, we integrated gene expression data with genetic association data to identify a single gene at the chromosome 7 locus as the driver of the phospholipid phenotypes. The gene encodes α/β-hydrolase domain 2 (Abhd2), one of 23 members of the ABHD gene family. We validated this observation by measuring lipids in a mouse with a whole-body deletion of Abhd2. The Abhd2KO mice had a significant increase in liver levels of phosphatidylcholine and phosphatidylethanolamine. Unexpectedly, we also found a decrease in two key mitochondrial lipids, cardiolipin and phosphatidylglycerol, in male Abhd2KO mice. These data suggest that Abhd2 plays a role in the synthesis, turnover, or remodeling of liver phospholipids.

Loss of ZNF148 enhances insulin secretion in human pancreatic β cells

Insulin secretion from pancreatic β cells is essential to the maintenance of glucose homeostasis. Defects in this process result in diabetes. Identifying genetic regulators that impair insulin secretion is crucial for the identification of novel therapeutic targets. Here, we show that reduction of ZNF148 in human islets, and its deletion in stem cell-derived β cells (SC-β cells), enhances insulin secretion. Transcriptomics of ZNF148-deficient SC-β cells identifies increased expression of annexin and S100 genes whose proteins form tetrameric complexes involved in regulation of insulin vesicle trafficking and exocytosis. ZNF148 in SC-β cells prevents translocation of annexin A2 from the nucleus to its functional place at the cell membrane via direct repression of S100A16 expression. These findings point to ZNF148 as a regulator of annexin-S100 complexes in human β cells and suggest that suppression of ZNF148 may provide a novel therapeutic strategy to enhance insulin secretion.

Genetic architecture of heart mitochondrial proteome influencing cardiac hypertrophy

Mitochondria play an important role in both normal heart function and disease etiology. We report analysis of common genetic variations contributing to mitochondrial and heart functions using an integrative proteomics approach in a panel of inbred mouse strains called the Hybrid Mouse Diversity Panel (HMDP). We performed a whole heart proteome study in the HMDP (72 strains, n=2-3 mice) and retrieved 848 mitochondrial proteins (quantified in ≥50 strains). High-resolution association mapping on their relative abundance levels revealed three trans-acting genetic loci on chromosomes (chr) 7, 13 and 17 that regulate distinct classes of mitochondrial proteins as well as cardiac hypertrophy. DAVID enrichment analyses of genes regulated by each of the loci revealed that the chr13 locus was highly enriched for complex-I proteins (24 proteins, P=2.2E-61), the chr17 locus for mitochondrial ribonucleoprotein complex (17 proteins, P=3.1E-25) and the chr7 locus for ubiquinone biosynthesis (3 proteins, P=6.9E-05). Follow-up high resolution regional mapping identified NDUFS4, LRPPRC and COQ7 as the candidate genes for chr13, chr17 and chr7 loci, respectively, and both experimental and statistical analyses supported their causal roles. Furthermore, a large cohort of Diversity Outbred mice was used to corroborate Lrpprc gene as a driver of mitochondrial DNA (mtDNA)-encoded gene regulation, and to show that the chr17 locus is specific to heart. Variations in all three loci were associated with heart mass in at least one of two independent heart stress models, namely, isoproterenol-induced heart failure and diet-induced obesity. These findings suggest that common variations in certain mitochondrial proteins can act in trans to influence tissue-specific mitochondrial functions and contribute to heart hypertrophy, elucidating mechanisms that may underlie genetic susceptibility to heart failure in human populations.

Lipidomic QTL in Diversity Outbred mice identifies a novel function for α/β hydrolase domain 2 ( Abhd2 ) as an enzyme that metabolizes phosphatidylcholine and cardiolipin

We and others have previously shown that genetic association can be used to make causal connections between gene loci and small molecules measured by mass spectrometry in the bloodstream and in tissues. We identified a locus on mouse chromosome 7 where several phospholipids in liver showed strong genetic association to distinct gene loci. In this study, we integrated gene expression data with genetic association data to identify a single gene at the chromosome 7 locus as the driver of the phospholipid phenotypes. The gene encodes α/β-hydrolase domain 2 ( Abhd2 ), one of 23 members of the ABHD gene family. We validated this observation by measuring lipids in a mouse with a whole-body deletion of Abhd2 . The Abhd2 KO mice had a significant increase in liver levels of phosphatidylcholine and phosphatidylethanolamine. Unexpectedly, we also found a decrease in two key mitochondrial lipids, cardiolipin and phosphatidylglycerol, in male Abhd2 KO mice. These data suggest that Abhd2 plays a role in the synthesis, turnover, or remodeling of liver phospholipids.

Pptc7 maintains mitochondrial protein content by suppressing receptor-mediated mitophagy

Pptc7 is a resident mitochondrial phosphatase essential for maintaining proper mitochondrial content and function. Newborn mice lacking Pptc7 exhibit aberrant mitochondrial protein phosphorylation, suffer from a range of metabolic defects, and fail to survive beyond one day after birth. Using an inducible knockout model, we reveal that loss of Pptc7 in adult mice causes marked reduction in mitochondrial mass concomitant with elevation of the mitophagy receptors Bnip3 and Nix. Consistently, Pptc7-/- mouse embryonic fibroblasts (MEFs) exhibit a major increase in mitophagy that is reversed upon deletion of these receptors. Our phosphoproteomics analyses reveal a common set of elevated phosphosites between perinatal tissues, adult liver, and MEFs-including multiple sites on Bnip3 and Nix. These data suggest that Pptc7 deletion causes mitochondrial dysfunction via dysregulation of several metabolic pathways and that Pptc7 may directly regulate mitophagy receptor function or stability. Overall, our work reveals a significant role for Pptc7 in the mitophagic response and furthers the growing notion that management of mitochondrial protein phosphorylation is essential for ensuring proper organelle content and function.

Sample multiplexing-based targeted pathway proteomics with real-time analytics reveals the impact of genetic variation on protein expression

Targeted proteomics enables hypothesis-driven research by measuring the cellular expression of protein cohorts related by function, disease, or class after perturbation. Here, we present a pathway-centric approach and an assay builder resource for targeting entire pathways of up to 200 proteins selected from >10,000 expressed proteins to directly measure their abundances, exploiting sample multiplexing to increase throughput by 16-fold. The strategy, termed GoDig, requires only a single-shot LC-MS analysis, ~1 µg combined peptide material, a list of up to 200 proteins, and real-time analytics to trigger simultaneous quantification of up to 16 samples for hundreds of analytes. We apply GoDig to quantify the impact of genetic variation on protein expression in mice fed a high-fat diet. We create several GoDig assays to quantify the expression of multiple protein families (kinases, lipid metabolism- and lipid droplet-associated proteins) across 480 fully-genotyped Diversity Outbred mice, revealing protein quantitative trait loci and establishing potential linkages between specific proteins and lipid homeostasis.

Sorting through the extensive and confusing roles of sortilin in metabolic disease

Sortilin is a post-Golgi trafficking receptor homologous to the yeast vacuolar protein sorting receptor 10 (VPS10). The VPS10 motif on sortilin is a 10-bladed β-propeller structure capable of binding more than 50 proteins, covering a wide range of biological functions including lipid and lipoprotein metabolism, neuronal growth and death, inflammation, and lysosomal degradation. Sortilin has a complex cellular trafficking itinerary, where it functions as a receptor in the trans-Golgi network, endosomes, secretory vesicles, multivesicular bodies, and at the cell surface. In addition, sortilin is associated with hypercholesterolemia, Alzheimer's disease, prion diseases, Parkinson's disease, and inflammation syndromes. The 1p13.3 locus containing SORT1, the gene encoding sortilin, carries the strongest association with LDL-C of all loci in human genome-wide association studies. However, the mechanism by which sortilin influences LDL-C is unclear. Here, we review the role sortilin plays in cardiovascular and metabolic diseases and describe in detail the large and often contradictory literature on the role of sortilin in the regulation of LDL-C levels.

QTLViewer: an interactive webtool for genetic analysis in the Collaborative Cross and Diversity Outbred mouse populations

The Collaborative Cross and the Diversity Outbred mouse populations are related multiparental populations, derived from the same 8 isogenic founder strains. They carry >50 M known genetic variants, which makes them ideal tools for mapping genetic loci that regulate phenotypes, including physiological and molecular traits. Mapping quantitative trait loci requires statistical and computational training, which can present a barrier to access for some researchers. The QTLViewer software allows users to graphically explore Collaborative Cross and Diversity Outbred quantitative trait locus mapping and related analyses performed through the R/qtl2 package. Additionally, the QTLViewer website serves as a repository for published Collaborative Cross and Diversity Outbred studies, increasing the accessibility of these genetic resources to the broader scientific community.

β Cell-specific deletion of Zfp148 improves nutrient-stimulated β cell Ca2+ responses

Insulin secretion from pancreatic β cells is essential for glucose homeostasis. An insufficient response to the demand for insulin results in diabetes. We previously showed that β cell-specific deletion of Zfp148 (β-Zfp148KO) improves glucose tolerance and insulin secretion in mice. Here, we performed Ca2+ imaging of islets from β‑Zfp148KO and control mice fed both a chow and a Western-style diet. β-Zfp148KO islets demonstrated improved sensitivity and sustained Ca2+ oscillations in response to elevated glucose levels. β-Zfp148KO islets also exhibited elevated sensitivity to amino acid-induced Ca2+ influx under low glucose conditions, suggesting enhanced mitochondrial phosphoenolpyruvate-dependent (PEP-dependent), ATP-sensitive K+ channel closure, independent of glycolysis. RNA-Seq and proteomics of β-Zfp148KO islets revealed altered levels of enzymes involved in amino acid metabolism (specifically, SLC3A2, SLC7A8, GLS, GLS2, PSPH, PHGDH, and PSAT1) and intermediary metabolism (namely, GOT1 and PCK2), consistent with altered PEP cycling. In agreement with this, β-Zfp148KO islets displayed enhanced insulin secretion in response to l-glutamine and activation of glutamate dehydrogenase. Understanding pathways controlled by ZFP148 may provide promising strategies for improving β cell function that are robust to the metabolic challenge imposed by a Western diet.

MetaNetwork Enhances Biological Insights from Quantitative Proteomics Differences by Combining Clustering and Enrichment Analyses

Interpreting proteomics data remains challenging due to the large number of proteins that are quantified by modern mass spectrometry methods. Weighted gene correlation network analysis (WGCNA) can identify groups of biologically related proteins using only protein intensity values by constructing protein correlation networks. However, WGCNA is not widespread in proteomic analyses due to challenges in implementing workflows. To facilitate the adoption of WGCNA by the proteomics field, we created MetaNetwork, an open-source, R-based application to perform sophisticated WGCNA workflows with no coding skill requirements for the end user. We demonstrate MetaNetwork's utility by employing it to identify groups of proteins associated with prostate cancer from a proteomic analysis of tumor and adjacent normal tissue samples. We found a decrease in cytoskeleton-related protein expression, a known hallmark of prostate tumors. We further identified changes in module eigenproteins indicative of dysregulation in protein translation and trafficking pathways. These results demonstrate the value of using MetaNetwork to improve the biological interpretation of quantitative proteomics experiments with 15 or more samples.

Coding variants identified in diabetic patients alter PICK1 BAR domain function in insulin granule biogenesis

Bin/amphiphysin/Rvs (BAR) domains are positively charged crescent-shaped modules that mediate curvature of negatively charged lipid membranes during remodeling processes. The BAR domain proteins PICK1, ICA69, and the arfaptins have recently been demonstrated to coordinate the budding and formation of immature secretory granules (ISGs) at the trans-Golgi network. Here, we identify 4 coding variants in the PICK1 gene from a whole-exome screening of Danish patients with diabetes that each involve a change in positively charged residues in the PICK1 BAR domain. All 4 coding variants failed to rescue insulin content in INS-1E cells upon knock down of endogenous PICK1. Moreover, 2 variants showed dominant-negative properties. In vitro assays addressing BAR domain function suggested that the coding variants compromised BAR domain function but increased the capacity to cause fission of liposomes. Live confocal microscopy and super-resolution microscopy further revealed that PICK1 resides transiently on ISGs before egress via vesicular budding events. Interestingly, this egress of PICK1 was accelerated in the coding variants. We propose that PICK1 assists in or complements the removal of excess membrane and generic membrane trafficking proteins, and possibly also insulin, from ISGs during the maturation process; and that the coding variants may cause premature budding, possibly explaining their dominant-negative function.

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

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

Identification of direct transcriptional targets of NFATC2 that promote β cell proliferation

The transcription factor NFATC2 induces β cell proliferation in mouse and human islets. However, the genomic targets that mediate these effects have not been identified. We expressed active forms of Nfatc2 and Nfatc1 in human islets. By integrating changes in gene expression with genomic binding sites for NFATC2, we identified approximately 2200 transcriptional targets of NFATC2. Genes induced by NFATC2 were enriched for transcripts that regulate the cell cycle and for DNA motifs associated with the transcription factor FOXP. Islets from an endocrine-specific Foxp1, Foxp2, and Foxp4 triple-knockout mouse were less responsive to NFATC2-induced β cell proliferation, suggesting the FOXP family works to regulate β cell proliferation in concert with NFATC2. NFATC2 induced β cell proliferation in both mouse and human islets, whereas NFATC1 did so only in human islets. Exploiting this species difference, we identified approximately 250 direct transcriptional targets of NFAT in human islets. This gene set enriches for cell cycle-associated transcripts and includes Nr4a1. Deletion of Nr4a1 reduced the capacity of NFATC2 to induce β cell proliferation, suggesting that much of the effect of NFATC2 occurs through its induction of Nr4a1. Integration of noncoding RNA expression, chromatin accessibility, and NFATC2 binding sites enabled us to identify NFATC2-dependent enhancer loci that mediate β cell proliferation.

Proteomic pathways to metabolic disease and type 2 diabetes in the pancreatic islet

Pancreatic islets are essential for maintaining physiological blood glucose levels, and declining islet function is a hallmark of type 2 diabetes. We employ mass spectrometry-based proteomics to systematically analyze islets from 9 genetic or diet-induced mouse models representing a broad cross-section of metabolic health. Quantifying the islet proteome to a depth of >11,500 proteins, this study represents the most detailed analysis of mouse islet proteins to date. Our data highlight that the majority of islet proteins are expressed in all strains and diets, but more than half of the proteins vary in expression levels, principally due to genetics. Associating these varied protein expression levels on an individual animal basis with individual phenotypic measures reveals islet mitochondrial function as a major positive indicator of metabolic health regardless of strain. This compendium of strain-specific and dietary changes to mouse islet proteomes represents a comprehensive resource for basic and translational islet cell biology.

Identification of sample mix-ups and mixtures in microbiome data in Diversity Outbred mice

In a Diversity Outbred mouse project with genotype data on 500 mice, including 297 with microbiome data, we identified three sets of sample mix-ups (two pairs and one trio) as well as at least 15 microbiome samples that appear to be mixtures of pairs of mice. The microbiome data consisted of shotgun sequencing reads from fecal DNA, used to characterize the gut microbial communities present in these mice. These sequence reads included sufficient reads derived from the host mouse to identify the individual. A number of microbiome samples appeared to contain a mixture of DNA from two mice. We describe a method for identifying sample mix-ups in such microbiome data, as well as a method for evaluating sample mixtures in this context.

Application of 2D IR Bioimaging: Hyperspectral Images of Formalin-Fixed Pancreatic Tissues and Observation of Slow Protein Degradation

We used two-dimensional IR bioimaging to study the structural heterogeneity of formalin-fixed mouse pancreas. Images were generated from the hyperspectral data sets by plotting quantities associated with the amide I vibrational mode, which is created by the backbone carbonyl stretch. Images that measure the fundamental vibrational frequencies, cross peaks, and anharmonic shifts are presented. Histograms are generated for each quantity, providing averaged values and distributions around the mean that serve as metrics for protein structures. Images were generated from tissue that had been stored in a formalin fixation for 3, 8, and 48 weeks. Over this period, all three metrics show that that the β-sheet content of the samples increased, consistent with protein aggregation. Our results indicate that formalin fixation does not entirely arrest the degradation of a protein structure in pancreas tissue.

INFIMA leverages multi-omics model organism data to identify effector genes of human GWAS variants

Genome-wide association studies reveal many non-coding variants associated with complex traits. However, model organism studies largely remain as an untapped resource for unveiling the effector genes of non-coding variants. We develop INFIMA, Integrative Fine-Mapping, to pinpoint causal SNPs for diversity outbred (DO) mice eQTL by integrating founder mice multi-omics data including ATAC-seq, RNA-seq, footprinting, and in silico mutation analysis. We demonstrate INFIMA's superior performance compared to alternatives with human and mouse chromatin conformation capture datasets. We apply INFIMA to identify novel effector genes for GWAS variants associated with diabetes. The results of the application are available at http://www.statlab.wisc.edu/shiny/INFIMA/ .

Reversal of hypertriglyceridemia in diabetic BTBR ob/ob mice does not prevent nephropathy

The etiology of diabetic nephropathy in type 2 diabetes is multifactorial. Sustained hyperglycemia is a major contributor, but additional contributions come from the hypertension, obesity, and hyperlipidemia that are also commonly present in patients with type 2 diabetes and nephropathy. The leptin deficient BTBR ob/ob mouse is a model of type 2 diabetic nephropathy in which hyperglycemia, obesity, and hyperlipidemia, but not hypertension, are present. We have shown that reversal of the constellation of these metabolic abnormalities with leptin replacement can reverse the morphologic and functional manifestations of diabetic nephropathy. Here we tested the hypothesis that reversal specifically of the hypertriglyceridemia, using an antisense oligonucleotide directed against ApoC-III, an apolipoprotein that regulates the interactions of VLDL (very low density lipoproteins) with the LDL receptor, is sufficient to ameliorate the nephropathy of Type 2 diabetes. Antisense treatment resulted in reduction of circulating ApoC-III protein levels and resulted in substantial lowering of triglycerides to near-normal levels in diabetic mice versus controls. Antisense treatment did not ameliorate proteinuria or pathologic manifestations of diabetic nephropathy, including podocyte loss. These studies indicate that pathologic manifestations of diabetic nephropathy are unlikely to be reduced by lipid-lowering therapeutics alone, but does not preclude a role for such interventions to be used in conjunction with other therapeutics commonly employed in the treatment of diabetes and its complications.

A large-scale genome-lipid association map guides lipid identification

Despite the crucial roles of lipids in metabolism, we are still at the early stages of comprehensively annotating lipid species and their genetic basis. Mass spectrometry-based discovery lipidomics offers the potential to globally survey lipids and their relative abundances in various biological samples. To discover the genetics of lipid features obtained through high-resolution liquid chromatography-tandem mass spectrometry, we analysed liver and plasma from 384 diversity outbred mice, and quantified 3,283 molecular features. These features were mapped to 5,622 lipid quantitative trait loci and compiled into a public web resource termed LipidGenie. The data are cross-referenced to the human genome and offer a bridge between genetic associations in humans and mice. Harnessing this resource, we used genome-lipid association data as an additional aid to identify a number of lipids, for example gangliosides through their association with B4galnt1, and found evidence for a group of sex-specific phosphatidylcholines through their shared locus. Finally, LipidGenie's ability to query either mass or gene-centric terms suggests acyl-chain-specific functions for proteins of the ABHD family.

Secretion of Recombinant Interleukin-22 by Engineered Lactobacillus reuteri Reduces Fatty Liver Disease in a Mouse Model of Diet-Induced Obesity

The incidence of metabolic syndrome continues to rise globally. In mice, intravenous administration of interleukin-22 (IL-22) ameliorates various disease phenotypes associated with diet-induced metabolic syndrome. In patients, oral treatment is favored over intravenous treatment, but methodologies to deliver IL-22 via the oral route are nonexistent. The goal of this study was to assess to what extent engineered Lactobacillus reuteri secreting IL-22 could ameliorate nonalcoholic fatty liver disease. We used a mouse model of diet-induced obesity and assessed various markers of metabolic syndrome following treatment with L. reuteri and a recombinant derivative. Mice that received an 8-week treatment of wild-type probiotic gained less weight and had a smaller fat pad than the control group, but these phenotypes were not further enhanced by recombinant L. reuteri However, L. reuteri secreting IL-22 significantly reduced liver weight and triglycerides at levels that exceeded those of the probiotic wild-type treatment group. Our findings are interesting in light of the observed phenotypes associated with reduced nonalcoholic liver disease, in humans the most prevalent chronic liver disease, following treatment of a next-generation probiotic that is administered orally. Once biological and environmental containment strategies are in place, therapeutic applications of recombinant Lactobacillus reuteri are on the horizon.IMPORTANCE In humans, nonalcoholic fatty liver disease (NAFLD) is the most prevalent liver disease due to the increased prevalence of obesity. While treatment of NAFLD is often geared toward lifestyle changes, such as diet and exercise, the use of dietary supplements such as probiotics is underinvestigated. Here, we report that probiotic Lactobacillus reuteri reduces fatty liver in a mouse model of diet-induced obesity. This phenotype was further enhanced upon delivery of recombinant interleukin-22 by engineered Lactobacillus reuteri These observations pave the road to a better understanding of probiotic mechanisms driving the reduction of diet-induced steatosis and to development of next-generation probiotics for use in the clinic. Ultimately, these studies may lead to rational selection of (engineered) probiotics to ameliorate fatty liver disease.

Genetic determinants of gut microbiota composition and bile acid profiles in mice

The microbial communities that inhabit the distal gut of humans and other mammals exhibit large inter-individual variation. While host genetics is a known factor that influences gut microbiota composition, the mechanisms underlying this variation remain largely unknown. Bile acids (BAs) are hormones that are produced by the host and chemically modified by gut bacteria. BAs serve as environmental cues and nutrients to microbes, but they can also have antibacterial effects. We hypothesized that host genetic variation in BA metabolism and homeostasis influence gut microbiota composition. To address this, we used the Diversity Outbred (DO) stock, a population of genetically distinct mice derived from eight founder strains. We characterized the fecal microbiota composition and plasma and cecal BA profiles from 400 DO mice maintained on a high-fat high-sucrose diet for ~22 weeks. Using quantitative trait locus (QTL) analysis, we identified several genomic regions associated with variations in both bacterial and BA profiles. Notably, we found overlapping QTL for Turicibacter sp. and plasma cholic acid, which mapped to a locus containing the gene for the ileal bile acid transporter, Slc10a2. Mediation analysis and subsequent follow-up validation experiments suggest that differences in Slc10a2 gene expression associated with the different strains influences levels of both traits and revealed novel interactions between Turicibacter and BAs. This work illustrates how systems genetics can be utilized to generate testable hypotheses and provide insight into host-microbe interactions.

Gene loci associated with insulin secretion in islets from non-diabetic mice

Genetic susceptibility to type 2 diabetes is primarily due to β-cell dysfunction. However, a genetic study to directly interrogate β-cell function ex vivo has never been previously performed. We isolated 233,447 islets from 483 Diversity Outbred (DO) mice maintained on a Western-style diet, and measured insulin secretion in response to a variety of secretagogues. Insulin secretion from DO islets ranged >1,000-fold even though none of the mice were diabetic. The insulin secretory response to each secretagogue had a unique genetic architecture; some of the loci were specific for one condition, whereas others overlapped. Human loci that are syntenic to many of the insulin secretion QTL from mouse are associated with diabetes-related SNPs in human genome-wide association studies. We report on three genes, Ptpn18, Hunk and Zfp148, where the phenotype predictions from the genetic screen were fulfilled in our studies of transgenic mouse models. These three genes encode a non-receptor type protein tyrosine phosphatase, a serine/threonine protein kinase, and a Krϋppel-type zinc-finger transcription factor, respectively. Our results demonstrate that genetic variation in insulin secretion that can lead to type 2 diabetes is discoverable in non-diabetic individuals.

Exploiting Prophage-Mediated Lysis for Biotherapeutic Release by Lactobacillus reuteri

Lactobacillus reuteri has the potential to be developed as a microbial therapeutic delivery platform because of an established safety profile, health-promoting properties, and available genome editing tools. Here, we show that L. reuteri VPL1014 exhibits a low mutation rate compared to other Gram-positive bacteria, which we expect will contribute to the stability of genetically modified strains. VPL1014 encodes two biologically active prophages, which are induced during gastrointestinal transit. We hypothesized that intracellularly accumulated recombinant protein can be released following bacteriophage-mediated lysis. To test this, we engineered VPL1014 to accumulate leptin, our model protein, inside the cell. In vitro prophage induction of recombinant VPL1014 released leptin into the extracellular milieu, which corresponded to bacteriophage production. We also employed a plasmid system that does not require antibiotic in the growth medium for plasmid maintenance. Collectively, these data provide new avenues to exploit native prophages to deliver therapeutic molecules.IMPORTANCE Lactic acid bacteria (LAB) have been explored as potential biotherapeutic vehicles for the past 20 years. To secrete a therapeutic in the extracellular milieu, one typically relies on the bacterial secretion pathway, i.e., the Sec pathway. Overexpression of a secreted protein can overload the secretory pathway and impact the organism's fitness, and optimization of the signal peptide is also required to maximize the efficiency of the release of mature protein. Here, we describe a previously unexplored approach to release therapeutics from the probiotic Lactobacillus reuteri We demonstrate that an intracellularly accumulated recombinant protein is released following prophage activation. Since we recently demonstrated that prophages are activated during gastrointestinal transit, we propose that this method will provide a straightforward and efficient approach to deliver therapeutics in vivo.

Dietary Fructose and Microbiota-Derived Short-Chain Fatty Acids Promote Bacteriophage Production in the Gut Symbiont Lactobacillus reuteri

The mammalian intestinal tract contains a complex microbial ecosystem with many lysogens, which are bacteria containing dormant phages (prophages) inserted within their genomes. Approximately half of intestinal viruses are derived from lysogens, suggesting that these bacteria encounter triggers that promote phage production. We show that prophages of the gut symbiont Lactobacillus reuteri are activated during gastrointestinal transit and that phage production is further increased in response to a fructose-enriched diet. Fructose and exposure to short-chain fatty acids activate the Ack pathway, involved in generating acetic acid, which in turn triggers the bacterial stress response that promotes phage production. L. reuteri mutants of the Ack pathway or RecA, a stress response component, exhibit decreased phage production. Thus, prophages in a gut symbiont can be induced by diet and metabolites affected by diet, which provides a potential mechanistic explanation for the effects of diet on the intestinal phage community.

Increased transport of acetyl-CoA into the endoplasmic reticulum causes a progeria-like phenotype

The membrane transporter AT-1/SLC33A1 translocates cytosolic acetyl-CoA into the lumen of the endoplasmic reticulum (ER), participating in quality control mechanisms within the secretory pathway. Mutations and duplication events in AT-1/SLC33A1 are highly pleiotropic and have been linked to diseases such as spastic paraplegia, developmental delay, autism spectrum disorder, intellectual disability, propensity to seizures, and dysmorphism. Despite these known associations, the biology of this key transporter is only beginning to be uncovered. Here, we show that systemic overexpression of AT-1 in the mouse leads to a segmental form of progeria with dysmorphism and metabolic alterations. The phenotype includes delayed growth, short lifespan, alopecia, skin lesions, rectal prolapse, osteoporosis, cardiomegaly, muscle atrophy, reduced fertility, and anemia. In terms of homeostasis, the AT-1 overexpressing mouse displays hypocholesterolemia, altered glycemia, and increased indices of systemic inflammation. Mechanistically, the phenotype is caused by a block in Atg9a-Fam134b-LC3β and Atg9a-Sec62-LC3β interactions, and defective reticulophagy, the autophagic recycling of the ER. Inhibition of ATase1/ATase2 acetyltransferase enzymes downstream of AT-1 restores reticulophagy and rescues the phenotype of the animals. These data suggest that inappropriately elevated acetyl-CoA flux into the ER directly induces defects in autophagy and recycling of subcellular structures and that this diversion of acetyl-CoA from cytosol to ER is causal in the progeria phenotype. Collectively, these data establish the cytosol-to-ER flux of acetyl-CoA as a novel event that dictates the pace of aging phenotypes and identify intracellular acetyl-CoA-dependent homeostatic mechanisms linked to metabolism and inflammation.

Islet proteomics reveals genetic variation in dopamine production resulting in altered insulin secretion

The mouse is a critical model in diabetes research, but most research in mice has been limited to a small number of mouse strains and limited genetic variation. Using the eight founder strains and both sexes of the Collaborative Cross (C57BL/6J (B6), A/J, 129S1/SvImJ (129), NOD/ShiLtJ (NOD), NZO/HILtJ (NZO), PWK/PhJ (PWK), WSB/EiJ (WSB), and CAST/EiJ (CAST)), we investigated the genetic dependence of diabetes-related metabolic phenotypes and insulin secretion. We found that strain background is associated with an extraordinary range in body weight, plasma glucose, insulin, triglycerides, and insulin secretion. Our whole-islet proteomic analysis of the eight mouse strains demonstrates that genetic background exerts a strong influence on the islet proteome that can be linked to the differences in diabetes-related metabolic phenotypes and insulin secretion. We computed protein modules consisting of highly correlated proteins that enrich for biological pathways and provide a searchable database of the islet protein expression profiles. To validate the data resource, we identified tyrosine hydroxylase (Th), a key enzyme in catecholamine synthesis, as a protein that is highly expressed in β-cells of PWK and CAST islets. We show that CAST islets synthesize elevated levels of dopamine, which suppresses insulin secretion. Prior studies, using only the B6 strain, concluded that adult mouse islets do not synthesize l-3,4-dihydroxyphenylalanine (l-DOPA), the product of Th and precursor of dopamine. Thus, the choice of the CAST strain, guided by our islet proteomic survey, was crucial for these discoveries. In summary, we provide a valuable data resource to the research community, and show that proteomic analysis identified a strain-specific pathway by which dopamine synthesized in β-cells inhibits insulin secretion.

Intracellular lipid metabolism impairs β cell compensation during diet-induced obesity

The compensatory proliferation of insulin-producing β cells is critical to maintaining glucose homeostasis at the early stage of type 2 diabetes. Failure of β cells to proliferate results in hyperglycemia and insulin dependence in patients. To understand the effect of the interplay between β cell compensation and lipid metabolism upon obesity and peripheral insulin resistance, we eliminated LDL receptor-related protein 1 (LRP1), a pleiotropic mediator of cholesterol, insulin, energy metabolism, and other cellular processes, in β cells. Upon high-fat diet exposure, LRP1 ablation significantly impaired insulin secretion and proliferation of β cells. The diminished insulin signaling was partly contributed to by the hypersensitivity to glucose-induced, Ca2+-dependent activation of Erk and the mTORC1 effector p85 S6K1. Surprisingly, in LRP1-deficient islets, lipotoxic sphingolipids were mitigated by improved lipid metabolism, mediated at least in part by the master transcriptional regulator PPARγ2. Acute overexpression of PPARγ2 in β cells impaired insulin signaling and insulin secretion. Elimination of Apbb2, a functional regulator of LRP1 cytoplasmic domain, also impaired β cell function in a similar fashion. In summary, our results uncover the double-edged effects of intracellular lipid metabolism on β cell function and viability in obesity and type 2 diabetes and highlight LRP1 as an essential regulator of these processes.

Perilipin 5 and liver fatty acid binding protein function to restore quiescence in mouse hepatic stellate cells

Hepatic stellate cell (HSC) activation occurs along with decreased Perilipin5 (Plin5) and liver fatty acid-binding protein (L-Fabp) expression and coincident lipid droplet (LD) depletion. Conversely, the activated phenotype is reversible in WT HSCs upon forced expression of Plin5. Here, we asked if L-Fabp expression is required for Plin5-mediated rescue of the quiescent phenotype. Lentiviral Plin5 transduction of passaged L-Fabp^-/- HSCs failed to reverse activation markers or restore lipogenic gene expression and LD formation. However, adenoviral L-Fabp infection of lentiviral Plin5 transduced L-Fabp^-/- HSCs restored both the quiescent phenotype and LD formation, an effect also mediated by adenoviral intestine-Fabp or adipocyte-Fabp. Expression of exogenous Plin5 in activated WT HSCs induced a transcriptional program of lipogenic gene expression including endogenous L-Fabp, but none of the other FABPs. We further demonstrated that selective, small molecule inhibition of endogenous L-Fabp also eliminated the ability of exogenous Plin5 to rescue LD formation and reverse activation of WT HSCs. This functional coordination of L-Fabp with Plin5 was 5'-AMP-activated protein kinase (AMPK)-dependent and was eliminated by AMPK inhibition. Taken together, our results indicate that L-Fabp is required for Plin5 to activate a transcriptional program that restores LD formation and reverses HSC activation.

Nat1 Deficiency Is Associated with Mitochondrial Dysfunction and Exercise Intolerance in Mice

We recently identified human N-acetyltransferase 2 (NAT2) as an insulin resistance (IR) gene. Here, we examine the cellular mechanism linking NAT2 to IR and find that Nat1 (mouse ortholog of NAT2) is co-regulated with key mitochondrial genes. RNAi-mediated silencing of Nat1 led to mitochondrial dysfunction characterized by increased intracellular reactive oxygen species and mitochondrial fragmentation as well as decreased mitochondrial membrane potential, biogenesis, mass, cellular respiration, and ATP generation. These effects were consistent in 3T3-L1 adipocytes, C2C12 myoblasts, and in tissues from Nat1-deficient mice, including white adipose tissue, heart, and skeletal muscle. Nat1-deficient mice had changes in plasma metabolites and lipids consistent with a decreased ability to utilize fats for energy and a decrease in basal metabolic rate and exercise capacity without altered thermogenesis. Collectively, our results suggest that Nat1 deficiency results in mitochondrial dysfunction, which may constitute a mechanistic link between this gene and IR.

Host Genotype and Gut Microbiome Modulate Insulin Secretion and Diet-Induced Metabolic Phenotypes

Genetic variation drives phenotypic diversity and influences the predisposition to metabolic disease. Here, we characterize the metabolic phenotypes of eight genetically distinct inbred mouse strains in response to a high-fat/high-sucrose diet. We found significant variation in diabetes-related phenotypes and gut microbiota composition among the different mouse strains in response to the dietary challenge and identified taxa associated with these traits. Follow-up microbiota transplant experiments showed that altering the composition of the gut microbiota modifies strain-specific susceptibility to diet-induced metabolic disease. Animals harboring microbial communities with enhanced capacity for processing dietary sugars and for generating hydrophobic bile acids showed increased susceptibility to metabolic disease. Notably, differences in glucose-stimulated insulin secretion between different mouse strains were partially recapitulated via gut microbiota transfer. Our results suggest that the gut microbiome contributes to the genetic and phenotypic diversity observed among mouse strains and provide a link between the gut microbiome and insulin secretion.

Targeted Mass Spectrometry Approach Enabled Discovery of O-Glycosylated Insulin and Related Signaling Peptides in Mouse and Human Pancreatic Islets

O-Linked glycosylation often involves the covalent attachment of sugar moieties to the hydroxyl group of serine or threonine on proteins/peptides. Despite growing interest in glycoproteins, little attention has been directed to glycosylated signaling peptides, largely due to lack of enabling analytical tools. Here we explore the occurrence of naturally O-linked glycosylation on the signaling peptides extracted from mouse and human pancreatic islets using mass spectrometry (MS). A novel targeted MS-based method is developed to increase the likelihood of capturing these modified signaling peptides and to provide improved sequence coverage and accurate glycosite localization, enabling the first large-scale discovery of O-glycosylation on signaling peptides. Several glycosylated signaling peptides with multiple glycoforms are identified, including the first report of glycosylated insulin-B chain and insulin-C peptide and BigLEN. This discovery may reveal potential novel functions as glycosylation could influence their conformation and biostability. Given the importance of insulin and its related peptide hormones and previous studies of glycosylated insulin analogues, this natural glycosylation may provide important insights into diabetes research and therapeutic treatments.

The Transcription Factor Nfatc2 Regulates β-Cell Proliferation and Genes Associated with Type 2 Diabetes in Mouse and Human Islets

Human genome-wide association studies (GWAS) have shown that genetic variation at >130 gene loci is associated with type 2 diabetes (T2D). We asked if the expression of the candidate T2D-associated genes within these loci is regulated by a common locus in pancreatic islets. Using an obese F2 mouse intercross segregating for T2D, we show that the expression of ~40% of the T2D-associated genes is linked to a broad region on mouse chromosome (Chr) 2. As all but 9 of these genes are not physically located on Chr 2, linkage to Chr 2 suggests a genomic factor(s) located on Chr 2 regulates their expression in trans. The transcription factor Nfatc2 is physically located on Chr 2 and its expression demonstrates cis linkage; i.e., its expression maps to itself. When conditioned on the expression of Nfatc2, linkage for the T2D-associated genes was greatly diminished, supporting Nfatc2 as a driver of their expression. Plasma insulin also showed linkage to the same broad region on Chr 2. Overexpression of a constitutively active (ca) form of Nfatc2 induced β-cell proliferation in mouse and human islets, and transcriptionally regulated more than half of the T2D-associated genes. Overexpression of either ca-Nfatc2 or ca-Nfatc1 in mouse islets enhanced insulin secretion, whereas only ca-Nfatc2 was able to promote β-cell proliferation, suggesting distinct molecular pathways mediating insulin secretion vs. β-cell proliferation are regulated by NFAT. Our results suggest that many of the T2D-associated genes are downstream transcriptional targets of NFAT, and may act coordinately in a pathway through which NFAT regulates β-cell proliferation in both mouse and human islets.

Genome-wide association studies (GWAS) and linkage studies provide a powerful way to establish a causal connection between a gene locus and a physiological or pathophysiological phenotype. We wondered if candidate genes associated with type 2 diabetes in human populations, in addition to being causal for the disease, could also be intermediate traits in a pathway leading to disease. In addition, we wished to know if there were any regulatory loci that could coordinately drive the expression of these genes in pancreatic islets and thus complete a pathway; i.e. Driver → GWAS candidate expression → type 2 diabetes. Using data from a mouse intercross between a diabetes-susceptible and a diabetes-resistant mouse strain, we found that the expression of ~40% of >130 candidate GWAS genes genetically mapped to a hot spot on mouse chromosome 2. Using a variety of statistical methods, we identified the transcription factor Nfatc2 as the candidate driver. Follow-up experiments showed that overexpression of Nfatc2 does indeed affect the expression of the GWAS genes and regulates β-cell proliferation and insulin secretion. The work shows that in addition to being causal, GWAS candidate genes can be intermediate traits in a pathway leading to disease. Model organisms can be used to explore these novel causal pathways.

Diet-Microbiota Interactions Mediate Global Epigenetic Programming in Multiple Host Tissues

Histone-modifying enzymes regulate transcription and are sensitive to availability of endogenous small-molecule metabolites, allowing chromatin to respond to changes in environment. The gut microbiota produces a myriad of metabolites that affect host physiology and susceptibility to disease; however, the underlying molecular events remain largely unknown. Here we demonstrate that microbial colonization regulates global histone acetylation and methylation in multiple host tissues in a diet-dependent manner: consumption of a "Western-type" diet prevents many of the microbiota-dependent chromatin changes that occur in a polysaccharide-rich diet. Finally, we demonstrate that supplementation of germ-free mice with short-chain fatty acids, major products of gut bacterial fermentation, is sufficient to recapitulate chromatin modification states and transcriptional responses associated with colonization. These findings have profound implications for understanding the complex functional interactions between diet, gut microbiota, and host health.

Histone chaperone ASF1B promotes human β-cell proliferation via recruitment of histone H3.3

Anti-silencing function 1 (ASF1) is a histone H3-H4 chaperone involved in DNA replication and repair, and transcriptional regulation. Here, we identify ASF1B, the mammalian paralog to ASF1, as a proliferation-inducing histone chaperone in human β-cells. Overexpression of ASF1B led to distinct transcriptional signatures consistent with increased cellular proliferation and reduced cellular death. Using multiple methods of monitoring proliferation and mitotic progression, we show that overexpression of ASF1B is sufficient to induce human β-cell proliferation. Co-expression of histone H3.3 further augmented β-cell proliferation, whereas suppression of endogenous H3.3 attenuated the stimulatory effect of ASF1B. Using the histone binding-deficient mutant of ASF1B (V94R), we show that histone binding to ASF1B is required for the induction of β-cell proliferation. In contrast to H3.3, overexpression of histone H3 variants H3.1 and H3.2 did not have an impact on ASF1B-mediated induction of proliferation. Our findings reveal a novel role of ASF1B in human β-cell replication and show that ASF1B and histone H3.3A synergistically stimulate human β-cell proliferation.

NeuCode Proteomics Reveals Bap1 Regulation of Metabolism

We introduce neutron-encoded (NeuCode) amino acid labeling of mice as a strategy for multiplexed proteomic analysis in vivo. Using NeuCode, we characterize an inducible knockout mouse model of Bap1, a tumor suppressor and deubiquitinase whose in vivo roles outside of cancer are not well established. NeuCode proteomics revealed altered metabolic pathways following Bap1 deletion, including profound elevation of cholesterol biosynthetic machinery coincident with reduced expression of gluconeogenic and lipid homeostasis proteins in liver. Bap1 loss increased pancreatitis biomarkers and reduced expression of mitochondrial proteins. These alterations accompany a metabolic remodeling with hypoglycemia, hypercholesterolemia, hepatic lipid loss, and acinar cell degeneration. Liver-specific Bap1 null mice present with fully penetrant perinatal lethality, severe hypoglycemia, and hepatic lipid deficiency. This work reveals Bap1 as a metabolic regulator in liver and pancreas, and it establishes NeuCode as a reliable proteomic method for deciphering in vivo biology.

Global Identification of Protein Post-translational Modifications in a Single-Pass Database Search

Bottom-up proteomics database search algorithms used for peptide identification cannot comprehensively identify post-translational modifications (PTMs) in a single-pass because of high false discovery rates (FDRs). A new approach to database searching enables global PTM (G-PTM) identification by exclusively looking for curated PTMs, thereby avoiding the FDR penalty experienced during conventional variable modification searches. We identified over 2200 unique, high-confidence modified peptides comprising 26 different PTM types in a single-pass database search.

Induction of miR-132 and miR-212 Expression by Glucagon-Like Peptide 1 (GLP-1) in Rodent and Human Pancreatic β-Cells

Better understanding how glucagon-like peptide 1 (GLP-1) promotes pancreatic β-cell function and/or mass may uncover new treatment for type 2 diabetes. In this study, we investigated the potential involvement of microRNAs (miRNAs) in the effect of GLP-1 on glucose-stimulated insulin secretion. miRNA levels in INS-1 cells and isolated rodent and human islets treated with GLP-1 in vitro and in vivo (with osmotic pumps) were measured by real-time quantitative PCR. The role of miRNAs on insulin secretion was studied by transfecting INS-1 cells with either precursors or antisense inhibitors of miRNAs. Among the 250 miRNAs surveyed, miR-132 and miR-212 were significantly up-regulated by GLP-1 by greater than 2-fold in INS-1 832/3 cells, which were subsequently reproduced in freshly isolated rat, mouse, and human islets, as well as the islets from GLP-1 infusion in vivo in mice. The inductions of miR-132 and miR-212 by GLP-1 were correlated with cAMP production and were blocked by the protein kinase A inhibitor H-89 but not affected by the exchange protein activated by cAMP activator 8-pCPT-2'-O-Me-cAMP-AM. GLP-1 failed to increase miR-132 or miR-212 expression levels in the 832/13 line of INS-1 cells, which lacks robust cAMP and insulin responses to GLP-1 treatment. Overexpression of miR-132 or miR-212 significantly enhanced glucose-stimulated insulin secretion in both 832/3 and 832/13 cells, and restored insulin responses to GLP-1 in INS-1 832/13 cells. GLP-1 increases the expression of miRNAs 132 and 212 via a cAMP/protein kinase A-dependent pathway in pancreatic β-cells. Overexpression of miR-132 or miR-212 enhances glucose and GLP-1-stimulated insulin secretion.