Transcripts with high distal heritability mediate genetic effects on complex metabolic traits.

Although many genes are subject to local regulation, recent evidence suggests that complex distal regulation may be more important in mediating phenotypic variability. To assess the role of distal gene regulation in complex traits, we combined multi-tissue transcriptomes with physiological outcomes to model diet-induced obesity and metabolic disease in a population of Diversity Outbred mice. Using a novel high-dimensional mediation analysis, we identified a composite transcriptome signature that summarized genetic effects on gene expression and explained 30% of the variation across all metabolic traits. The signature was heritable, interpretable in biological terms, and predicted obesity status from gene expression in an independently derived mouse cohort and multiple human studies. Transcripts contributing most strongly to this composite mediator frequently had complex, distal regulation distributed throughout the genome. These results suggest that trait-relevant variation in transcription is largely distally regulated, but is nonetheless identifiable, interpretable, and translatable across species.

Multiple promoter and enhancer differences likely contribute to augmented G6PC2 expression in human versus mouse pancreatic islet alpha cells.

G6PC2 encodes a glucose-6-phosphatase catalytic subunit that opposes the action of glucokinase in pancreatic islets, thereby modulating the sensitivity of insulin and glucagon secretion to glucose. In mice, G6pc2 is expressed at ~20-fold higher levels in ß-cells than in α-cells, whereas in humans G6PC2 is expressed at only ~5-fold higher levels in ß-cells. We therefore hypothesize that G6PC2 likely influences glucagon secretion to a greater degree in humans. With a view to generating a humanized mouse that recapitulates augmented G6PC2 expression levels in α-cells, we sought to identify the genomic regions that confer differential mouse G6pc2 expression in α-cells versus ß-cells as well as the evolutionary changes that have altered this ratio in humans. Studies in islet-derived cell lines suggest that the elevated G6pc2 expression in mouse ß-cells versus α-cells is mainly due to a difference in the relative activity of the proximal G6pc2 promoter in these cell types. Similarly, the smaller difference in G6PC2 expression between α-cells and ß-cells in humans is potentially explained by a change in relative proximal G6PC2 promoter activity. However, we show that both glucocorticoid levels and multiple differences in the relative activity of eight transcriptional enhancers between mice and humans likely contribute to differential G6PC2 expression. Finally, we show that a mouse-specific non-coding RNA, Gm13613, whose expression is controlled by G6pc2 enhancer I, does not regulate G6pc2 expression, indicating that altered expression of Gm13613 in a humanized mouse that contains both the human promoter and enhancers should not affect G6PC2 function.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.