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.

Genetic mapping of microbial and host traits reveals production of immunomodulatory lipids by Akkermansia muciniphila in the murine gut.

The molecular bases of how host genetic variation impacts the gut microbiome remain largely unknown. Here we used a genetically diverse mouse population and applied systems genetics strategies to identify interactions between host and microbe phenotypes including microbial functions, using faecal metagenomics, small intestinal transcripts and caecal lipids that influence microbe-host dynamics. Quantitative trait locus (QTL) mapping identified murine genomic regions associated with variations in bacterial taxa; bacterial functions including motility, sporulation and lipopolysaccharide production and levels of bacterial- and host-derived lipids. We found overlapping QTL for the abundance of Akkermansia muciniphila and caecal levels of ornithine lipids. Follow-up in vitro and in vivo studies revealed that A. muciniphila is a major source of these lipids in the gut, provided evidence that ornithine lipids have immunomodulatory effects and identified intestinal transcripts co-regulated with these traits including Atf3, which encodes for a transcription factor that plays vital roles in modulating metabolism and immunity. Collectively, these results suggest that ornithine lipids are potentially important for A. muciniphila-host interactions and support the role of host genetics as a determinant of responses to gut microbes.

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.

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.

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.

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.

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.

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.

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.

Identification of the Bile Acid Transporter Slco1a6 as a Candidate Gene That Broadly Affects Gene Expression in Mouse Pancreatic Islets.

We surveyed gene expression in six tissues in an F2 intercross between mouse strains C57BL/6J (abbreviated B6) and BTBR T(+) tf/J (abbreviated BTBR) made genetically obese with the Leptin(ob) mutation. We identified a number of expression quantitative trait loci (eQTL) affecting the expression of numerous genes distal to the locus, called trans-eQTL hotspots. Some of these trans-eQTL hotspots showed effects in multiple tissues, whereas some were specific to a single tissue. An unusually large number of transcripts (∼8% of genes) mapped in trans to a hotspot on chromosome 6, specifically in pancreatic islets. By considering the first two principal components of the expression of genes mapping to this region, we were able to convert the multivariate phenotype into a simple Mendelian trait. Fine mapping the locus by traditional methods reduced the QTL interval to a 298-kb region containing only three genes, including Slco1a6, one member of a large family of organic anion transporters. Direct genomic sequencing of all Slco1a6 exons identified a nonsynonymous coding SNP that converts a highly conserved proline residue at amino acid position 564 to serine. Molecular modeling suggests that Pro564 faces an aqueous pore within this 12-transmembrane domain-spanning protein. When transiently overexpressed in HEK293 cells, BTBR organic anion transporting polypeptide (OATP)1A6-mediated cellular uptake of the bile acid taurocholic acid (TCA) was enhanced compared to B6 OATP1A6. Our results suggest that genetic variation in Slco1a6 leads to altered transport of TCA (and potentially other bile acids) by pancreatic islets, resulting in broad gene regulation.

Quantification of mitochondrial acetylation dynamics highlights prominent sites of metabolic regulation.

Lysine acetylation is rapidly becoming established as a key post-translational modification for regulating mitochondrial metabolism. Nonetheless, distinguishing regulatory sites from among the thousands identified by mass spectrometry and elucidating how these modifications alter enzyme function remain primary challenges. Here, we performed multiplexed quantitative mass spectrometry to measure changes in the mouse liver mitochondrial acetylproteome in response to acute and chronic alterations in nutritional status, and integrated these data sets with our compendium of predicted Sirt3 targets. These analyses highlight a subset of mitochondrial proteins with dynamic acetylation sites, including acetyl-CoA acetyltransferase 1 (Acat1), an enzyme central to multiple metabolic pathways. We performed in vitro biochemistry and molecular modeling to demonstrate that acetylation of Acat1 decreases its activity by disrupting the binding of coenzyme A. Collectively, our data reveal an important new target of regulatory acetylation and provide a foundation for investigating the role of select mitochondrial protein acetylation sites in mediating acute and chronic metabolic transitions.

A quantitative map of the liver mitochondrial phosphoproteome reveals posttranslational control of ketogenesis.

Mitochondria are dynamic organelles that play a central role in a diverse array of metabolic processes. Elucidating mitochondrial adaptations to changing metabolic demands and the pathogenic alterations that underlie metabolic disorders represent principal challenges in cell biology. Here, we performed multiplexed quantitative mass spectrometry-based proteomics to chart the remodeling of the mouse liver mitochondrial proteome and phosphoproteome during both acute and chronic physiological transformations in more than 50 mice. Our analyses reveal that reversible phosphorylation is widespread in mitochondria, and is a key mechanism for regulating ketogenesis during the onset of obesity and type 2 diabetes. Specifically, we have demonstrated that phosphorylation of a conserved serine on Hmgcs2 (S456) significantly enhances its catalytic activity in response to increased ketogenic demand. Collectively, our work describes the plasticity of this organelle at high resolution and provides a framework for investigating the roles of proteome restructuring and reversible phosphorylation in mitochondrial adaptation.

Tsc2, a positional candidate gene underlying a quantitative trait locus for hepatic steatosis.

Nonalchoholic fatty liver disease (NAFLD) is the most common cause of liver dysfunction and is associated with metabolic diseases, including obesity, insulin resistance, and type 2 diabetes. We mapped a quantitative trait locus (QTL) for NAFLD to chromosome 17 in a cross between C57BL/6 (B6) and BTBR mouse strains made genetically obese with the Lep(ob/ob) mutation. We identified Tsc2 as a gene underlying the chromosome 17 NAFLD QTL. Tsc2 functions as an inhibitor of mammalian target of rapamycin, which is involved in many physiological processes, including cell growth, proliferation, and metabolism. We found that Tsc2(+/-) mice have increased lipogenic gene expression in the liver in an insulin-dependent manner. The coding single nucleotide polymorphism between the B6 and BTBR strains leads to a change in the ability to inhibit the expression of lipogenic genes and de novo lipogenesis in AML12 cells and to promote the proliferation of Ins1 cells. This difference is due to a different affinity of binding to Tsc1, which affects the stability of Tsc2.

Positional cloning of a type 2 diabetes quantitative trait locus; tomosyn-2, a negative regulator of insulin secretion.

We previously mapped a type 2 diabetes (T2D) locus on chromosome 16 (Chr 16) in an F2 intercross from the BTBR T (+) tf (BTBR) Lep(ob/ob) and C57BL/6 (B6) Lep(ob/ob) mouse strains. Introgression of BTBR Chr 16 into B6 mice resulted in a consomic mouse with reduced fasting plasma insulin and elevated glucose levels. We derived a panel of sub-congenic mice and narrowed the diabetes susceptibility locus to a 1.6 Mb region. Introgression of this 1.6 Mb fragment of the BTBR Chr 16 into lean B6 mice (B6.16(BT36-38)) replicated the phenotypes of the consomic mice. Pancreatic islets from the B6.16(BT36-38) mice were defective in the second phase of the insulin secretion, suggesting that the 1.6 Mb region encodes a regulator of insulin secretion. Within this region, syntaxin-binding protein 5-like (Stxbp5l) or tomosyn-2 was the only gene with an expression difference and a non-synonymous coding single nucleotide polymorphism (SNP) between the B6 and BTBR alleles. Overexpression of the b-tomosyn-2 isoform in the pancreatic ß-cell line, INS1 (832/13), resulted in an inhibition of insulin secretion in response to 3 mM 8-bromo cAMP at 7 mM glucose. In vitro binding experiments showed that tomosyn-2 binds recombinant syntaxin-1A and syntaxin-4, key proteins that are involved in insulin secretion via formation of the SNARE complex. The B6 form of tomosyn-2 is more susceptible to proteasomal degradation than the BTBR form, establishing a functional role for the coding SNP in tomosyn-2. We conclude that tomosyn-2 is the major gene responsible for the T2D Chr 16 quantitative trait locus (QTL) we mapped in our mouse cross. Our findings suggest that tomosyn-2 is a key negative regulator of insulin secretion.

Loss of PDGF-B activity increases hepatic vascular permeability and enhances insulin sensitivity.

Hepatic vasculature is not thought to pose a permeability barrier for diffusion of macromolecules from the bloodstream to hepatocytes. In contrast, in extrahepatic tissues, the microvasculature is critically important for insulin action, because transport of insulin across the endothelial cell layer is rate limiting for insulin-stimulated glucose disposal. However, very little is known concerning the role in this process of pericytes, the mural cells lining the basolateral membrane of endothelial cells. PDGF-B is a growth factor involved in the recruitment and function of pericytes. We studied insulin action in mice expressing PDGF-B lacking the proteoglycan binding domain, producing a protein with a partial loss of function (PDGF-B(ret/ret)). Insulin action was assessed through measurements of insulin signaling and insulin and glucose tolerance tests. PDGF-B deficiency enhanced hepatic vascular transendothelial transport. One outcome of this change was an increase in hepatic insulin signaling. This correlated with enhanced whole body glucose homeostasis and increased insulin clearance from the circulation during an insulin tolerance test. In obese mice, PDGF-B deficiency was associated with an 80% reduction in fasting insulin and drastically reduced insulin secretion. These mice did not have significantly higher glucose levels, reflecting a dramatic increase in insulin action. Our findings show that, despite already having a high permeability, hepatic transendothelial transport can be further enhanced. To the best of our knowledge, this is the first study to connect PDGF-B-induced changes in hepatic sinusoidal transport to changes in insulin action, demonstrating a link between PDGF-B signaling and insulin sensitivity.

Cholecystokinin is up-regulated in obese mouse islets and expands beta-cell mass by increasing beta-cell survival.

An absolute or functional deficit in beta-cell mass is a key factor in the pathogenesis of diabetes. We model obesity-driven beta-cell mass expansion by studying the diabetes-resistant C57BL/6-Leptin(ob/ob) mouse. We previously reported that cholecystokinin (Cck) was the most up-regulated gene in obese pancreatic islets. We now show that islet cholecystokinin (CCK) is up-regulated 500-fold by obesity and expressed in both alpha- and beta-cells. We bred a null Cck allele into the C57BL/6-Leptin(ob/ob) background and investigated beta-cell mass and metabolic parameters of Cck-deficient obese mice. Loss of CCK resulted in decreased islet size and reduced beta-cell mass through increased beta-cell death. CCK deficiency and decreased beta-cell mass exacerbated fasting hyperglycemia and reduced hyperinsulinemia. We further investigated whether CCK can directly affect beta-cell death in cell culture and isolated islets. CCK was able to directly reduce cytokine- and endoplasmic reticulum stress-induced cell death. In summary, CCK is up-regulated by islet cells during obesity and functions as a paracrine or autocrine factor to increase beta-cell survival and expand beta-cell mass to compensate for obesity-induced insulin resistance.

A gene expression network model of type 2 diabetes links cell cycle regulation in islets with diabetes susceptibility.

Insulin resistance is necessary but not sufficient for the development of type 2 diabetes. Diabetes results when pancreatic beta-cells fail to compensate for insulin resistance by increasing insulin production through an expansion of beta-cell mass or increased insulin secretion. Communication between insulin target tissues and beta-cells may initiate this compensatory response. Correlated changes in gene expression between tissues can provide evidence for such intercellular communication. We profiled gene expression in six tissues of mice from an obesity-induced diabetes-resistant and a diabetes-susceptible strain before and after the onset of diabetes. We studied the correlation structure of mRNA abundance and identified 105 co-expression gene modules. We provide an interactive gene network model showing the correlation structure between the expression modules within and among the six tissues. This resource also provides a searchable database of gene expression profiles for all genes in six tissues in lean and obese diabetes-resistant and diabetes-susceptible mice, at 4 and 10 wk of age. A cell cycle regulatory module in islets predicts diabetes susceptibility. The module predicts islet replication; we found a strong correlation between (2)H(2)O incorporation into islet DNA in vivo and the expression pattern of the cell cycle module. This pattern is highly correlated with that of several individual genes in insulin target tissues, including Igf2, which has been shown to promote beta-cell proliferation, suggesting that these genes may provide a link between insulin resistance and beta-cell proliferation.

Positional cloning of Sorcs1, a type 2 diabetes quantitative trait locus.

We previously mapped the type 2 diabetes mellitus-2 locus (T2dm2), which affects fasting insulin levels, to distal chromosome 19 in a leptin-deficient obese F2 intercross derived from C57BL/6 (B6) and BTBR T+ tf/J (BTBR) mice. Introgression of a 7-Mb segment of the B6 chromosome 19 into the BTBR background (strain 1339A) replicated the reduced insulin linked to T2dm2. The 1339A mice have markedly impaired insulin secretion in vivo and disrupted islet morphology. We used subcongenic strains derived from 1339A to localize the T2dm2 quantitative trait locus (QTL) to a 242-kb segment comprising the promoter, first exon and most of the first intron of the Sorcs1 gene. This was the only gene in the 1339A strain for which we detected amino acid substitutions and expression level differences between mice carrying B6 and BTBR alleles of this insert, thereby identifying variation within the Sorcs1 gene as underlying the phenotype associated with the T2dm2 locus. SorCS1 binds platelet-derived growth factor, a growth factor crucial for pericyte recruitment to the microvasculature, and may thus have a role in expanding or maintaining the islet vasculature. Our identification of the Sorcs1 gene provides insight into the pathway underlying the pathophysiology of obesity-induced type 2 diabetes mellitus.