Genetic Drivers of Pancreatic Islet Function

Population variability in X-chromosome inactivation across 10 mammalian species

One of the two X-chromosomes in female mammals is epigenetically silenced in embryonic stem cells by X-chromosome inactivation. This creates a mosaic of cells expressing either the maternal or the paternal X allele. The X-chromosome inactivation ratio, the proportion of inactivated parental alleles, varies widely among individuals, representing the largest instance of epigenetic variability within mammalian populations. While various contributing factors to X-chromosome inactivation variability are recognized, namely stochastic and/or genetic effects, their relative contributions are poorly understood. This is due in part to limited cross-species analysis, making it difficult to distinguish between generalizable or species-specific mechanisms for X-chromosome inactivation ratio variability. To address this gap, we measure X-chromosome inactivation ratios in ten mammalian species (9531 individual samples), ranging from rodents to primates, and compare the strength of stochastic models or genetic factors for explaining X-chromosome inactivation variability. Our results demonstrate the embryonic stochasticity of X-chromosome inactivation is a general explanatory model for population X-chromosome inactivation variability in mammals, while genetic factors play a minor role.

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.

A PLURIPOTENT STEM CELL PLATFORM FOR IN VITRO SYSTEMS GENETICS STUDIES OF MOUSE DEVELOPMENT

The directed differentiation of pluripotent stem cells (PSCs) from panels of genetically diverse individuals is emerging as a powerful experimental system for characterizing the impact of natural genetic variation on developing cell types and tissues. Here, we establish new PSC lines and experimental approaches for modeling embryonic development in a genetically diverse, outbred mouse stock (Diversity Outbred mice). We show that a range of inbred and outbred PSC lines can be stably maintained in the primed pluripotent state (epiblast stem cells -- EpiSCs) and establish the contribution of genetic variation to phenotypic differences in gene regulation and directed differentiation. Using pooled in vitro fertilization, we generate and characterize a genetic reference panel of Diversity Outbred PSCs (n = 230). Finally, we demonstrate the feasibility of pooled culture of Diversity Outbred EpiSCs as "cell villages", which can facilitate the differentiation of large numbers of EpiSC lines for forward genetic screens. These data can complement and inform similar efforts within the stem cell biology and human genetics communities to model the impact of natural genetic variation on phenotypic variation and disease-risk.

An in vitro neurogenetics platform for precision disease modeling in the mouse

The power and scope of disease modeling can be markedly enhanced through the incorporation of broad genetic diversity. The introduction of pathogenic mutations into a single inbred mouse strain sometimes fails to mimic human disease. We describe a cross-species precision disease modeling platform that exploits mouse genetic diversity to bridge cell-based modeling with whole organism analysis. We developed a universal protocol that permitted robust and reproducible neural differentiation of genetically diverse human and mouse pluripotent stem cell lines and then carried out a proof-of-concept study of the neurodevelopmental gene DYRK1A. Results in vitro reliably predicted the effects of genetic background on Dyrk1a loss-of-function phenotypes in vivo. Transcriptomic comparison of responsive and unresponsive strains identified molecular pathways conferring sensitivity or resilience to Dyrk1a1A loss and highlighted differential messenger RNA isoform usage as an important determinant of response. This cross-species strategy provides a powerful tool in the functional analysis of candidate disease variants identified through human genetic studies.

GenomeMUSter mouse genetic variation service enables multitrait, multipopulation data integration and analysis

Hundreds of inbred mouse strains and intercross populations have been used to characterize the function of genetic variants that contribute to disease. Thousands of disease-relevant traits have been characterized in mice and made publicly available. New strains and populations including consomics, the collaborative cross, expanded BXD, and inbred wild-derived strains add to existing complex disease mouse models, mapping populations, and sensitized backgrounds for engineered mutations. The genome sequences of inbred strains, along with dense genotypes from others, enable integrated analysis of trait-variant associations across populations, but these analyses are hampered by the sparsity of genotypes available. Moreover, the data are not readily interoperable with other resources. To address these limitations, we created a uniformly dense variant resource by harmonizing multiple data sets. Missing genotypes were imputed using the Viterbi algorithm with a data-driven technique that incorporates local phylogenetic information, an approach that is extendable to other model organisms. The result is a web- and programmatically accessible data service called GenomeMUSter, comprising single-nucleotide variants covering 657 strains at 106.8 million segregating sites. Interoperation with phenotype databases, analytic tools, and other resources enable a wealth of applications, including multitrait, multipopulation meta-analysis. We show this in cross-species comparisons of type 2 diabetes and substance use disorder meta-analyses, leveraging mouse data to characterize the likely role of human variant effects in disease. Other applications include refinement of mapped loci and prioritization of strain backgrounds for disease modeling to further unlock extant mouse diversity for genetic and genomic studies in health and disease.

Genome-wide association studies of human and rat BMI converge on synapse, epigenome, and hormone signaling networks

A vexing observation in genome-wide association studies (GWASs) is that parallel analyses in different species may not identify orthologous genes. Here, we demonstrate that cross-species translation of GWASs can be greatly improved by an analysis of co-localization within molecular networks. Using body mass index (BMI) as an example, we show that the genes associated with BMI in humans lack significant agreement with those identified in rats. However, the networks interconnecting these genes show substantial overlap, highlighting common mechanisms including synaptic signaling, epigenetic modification, and hormonal regulation. Genetic perturbations within these networks cause abnormal BMI phenotypes in mice, too, supporting their broad conservation across mammals. Other mechanisms appear species specific, including carbohydrate biosynthesis (humans) and glycerolipid metabolism (rodents). Finally, network co-localization also identifies cross-species convergence for height/body length. This study advances a general paradigm for determining whether and how phenotypes measured in model species recapitulate human biology.

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.

GenomeMUSter mouse genetic variation service enables multi-trait, multi-population data integration and analyses

Hundreds of inbred laboratory mouse strains and intercross populations have been used to functionalize genetic variants that contribute to disease. Thousands of disease relevant traits have been characterized in mice and made publicly available. New strains and populations including the Collaborative Cross, expanded BXD and inbred wild-derived strains add to set of complex disease mouse models, genetic mapping resources and sensitized backgrounds against which to evaluate engineered mutations. The genome sequences of many inbred strains, along with dense genotypes from others could allow integrated analysis of trait - variant associations across populations, but these analyses are not feasible due to the sparsity of genotypes available. Moreover, the data are not readily interoperable with other resources. To address these limitations, we created a uniformly dense data resource by harmonizing multiple variant datasets. Missing genotypes were imputed using the Viterbi algorithm with a data-driven technique that incorporates local phylogenetic information, an approach that is extensible to other model organism species. The result is a web- and programmatically-accessible data service called GenomeMUSter ( https://muster.jax.org ), comprising allelic data covering 657 strains at 106.8M segregating sites. Interoperation with phenotype databases, analytic tools and other resources enable a wealth of applications including multi-trait, multi-population meta-analysis. We demonstrate this in a cross-species comparison of the meta-analysis of Type 2 Diabetes and of substance use disorders, resulting in the more specific characterization of the role of human variant effects in light of mouse phenotype data. Other applications include refinement of mapped loci and prioritization of strain backgrounds for disease modeling to further unlock extant mouse diversity for genetic and genomic studies in health and disease.

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.

The trouble with triples: Examining the impact of measurement error in mediation analysis

Mediation analysis is used in genetic mapping studies to identify candidate gene mediators of quantitative trait loci (QTL). We consider genetic mediation analysis of triplets-sets of three variables consisting of a target trait, the genotype at a QTL for the target trait, and a candidate mediator that is the abundance of a transcript or protein whose coding gene co-locates with the QTL. We show that, in the presence of measurement error, mediation analysis can infer partial mediation even in the absence of a causal relationship between the candidate mediator and the target. We describe a measurement error model and a corresponding latent variable model with estimable parameters that are combinations of the causal effects and measurement errors across all three variables. The relative magnitudes of the latent variable correlations determine whether or not mediation analysis will tend to infer the correct causal relationship in large samples. We examine case studies that illustrate the common failure modes of genetic mediation analysis and demonstrate how to evaluate the effects of measurement error. While genetic mediation analysis is a powerful tool for identifying candidate genes, we recommend caution when interpreting mediation analysis findings.

Which mouse multiparental population is right for your study? The Collaborative Cross inbred strains, their F1 hybrids, or the Diversity Outbred population

Multiparental populations (MPPs) encompass greater genetic diversity than traditional experimental crosses of two inbred strains, enabling broader surveys of genetic variation underlying complex traits. Two such mouse MPPs are the Collaborative Cross (CC) inbred panel and the Diversity Outbred (DO) population, which are descended from the same eight inbred strains. Additionally, the F1 intercrosses of CC strains (CC-RIX) have been used and enable study designs with replicate outbred mice. Genetic analyses commonly used by researchers to investigate complex traits in these populations include characterizing how heritable a trait is, i.e. its heritability, and mapping its underlying genetic loci, i.e. its quantitative trait loci (QTLs). Here we evaluate the relative merits of these populations for these tasks through simulation, as well as provide recommendations for performing the quantitative genetic analyses. We find that sample populations that include replicate animals, as possible with the CC and CC-RIX, provide more efficient and precise estimates of heritability. We report QTL mapping power curves for the CC, CC-RIX, and DO across a range of QTL effect sizes and polygenic backgrounds for samples of 174 and 500 mice. The utility of replicate animals in the CC and CC-RIX for mapping QTLs rapidly decreased as traits became more polygenic. Only large sample populations of 500 DO mice were well-powered to detect smaller effect loci (7.5-10%) for highly complex traits (80% polygenic background). All results were generated with our R package musppr, which we developed to simulate data from these MPPs and evaluate genetic analyses from user-provided genotypes.

The potential of integrating human and mouse discovery platforms to advance our understanding of cardiometabolic diseases

Cardiometabolic diseases encompass a range of interrelated conditions that arise from underlying metabolic perturbations precipitated by genetic, environmental, and lifestyle factors. While obesity, dyslipidaemia, smoking, and insulin resistance are major risk factors for cardiometabolic diseases, individuals still present in the absence of such traditional risk factors, making it difficult to determine those at greatest risk of disease. Thus, it is crucial to elucidate the genetic, environmental, and molecular underpinnings to better understand, diagnose, and treat cardiometabolic diseases. Much of this information can be garnered using systems genetics, which takes population-based approaches to investigate how genetic variance contributes to complex traits. Despite the important advances made by human genome-wide association studies (GWAS) in this space, corroboration of these findings has been hampered by limitations including the inability to control environmental influence, limited access to pertinent metabolic tissues, and often, poor classification of diseases or phenotypes. A complementary approach to human GWAS is the utilisation of model systems such as genetically diverse mouse panels to study natural genetic and phenotypic variation in a controlled environment. Here, we review mouse genetic reference panels and the opportunities they provide for the study of cardiometabolic diseases and related traits. We discuss how the post-GWAS era has prompted a shift in focus from discovery of novel genetic variants to understanding gene function. Finally, we highlight key advantages and challenges of integrating complementary genetic and multi-omics data from human and mouse populations to advance biological discovery.

Integrated transcriptomics contrasts fatty acid metabolism with hypoxia response in β-cell subpopulations associated with glycemic control

BACKGROUND

Understanding how heterogeneous β-cell function impacts diabetes is imperative for therapy development. Standard single-cell RNA sequencing analysis illuminates some factors driving heterogeneity, but new strategies are required to enhance information capture.

RESULTS

We integrate pancreatic islet single-cell and bulk RNA sequencing data to identify β-cell subpopulations based on gene expression and characterize genetic networks associated with β-cell function in obese SM/J mice. We identify β-cell subpopulations associated with basal insulin secretion, hypoxia response, cell polarity, and stress response. Network analysis associates fatty acid metabolism and basal insulin secretion with hyperglycemic-obesity, while expression of Pdyn and hypoxia response is associated with normoglycemic-obesity.

CONCLUSIONS

By integrating single-cell and bulk islet transcriptomes, our study explores β-cell heterogeneity and identifies novel subpopulations and genetic pathways associated with β-cell function in obesity.

Multi-omics analysis identifies drivers of protein phosphorylation

BACKGROUND

Phosphorylation of proteins is a key step in the regulation of many cellular processes including activation of enzymes and signaling cascades. The abundance of a phosphorylated peptide (phosphopeptide) is determined by the abundance of its parent protein and the proportion of target sites that are phosphorylated.

RESULTS

We quantified phosphopeptides, proteins, and transcripts in heart, liver, and kidney tissue samples of mice from 58 strains of the Collaborative Cross strain panel. We mapped ~700 phosphorylation quantitative trait loci (phQTL) across the three tissues and applied genetic mediation analysis to identify causal drivers of phosphorylation. We identified kinases, phosphatases, cytokines, and other factors, including both known and potentially novel interactions between target proteins and genes that regulate site-specific phosphorylation. Our analysis highlights multiple targets of pyruvate dehydrogenase kinase 1 (PDK1), a regulator of mitochondrial function that shows reduced activity in the NZO/HILtJ mouse, a polygenic model of obesity and type 2 diabetes.

CONCLUSIONS

Together, this integrative multi-omics analysis in genetically diverse CC strains provides a powerful tool to identify regulators of protein phosphorylation. The data generated in this study provides a resource for further exploration.

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.

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.

Western diet augments metabolic and arterial dysfunction in a sex-specific manner in outbred, genetically diverse mice

Western diet (WD), characterized by excess saturated fat and sugar intake, is a major contributor to obesity and metabolic and arterial dysfunction in humans. However, these phenotypes are not consistently observed in traditional inbred, genetically identical mice. Therefore, we sought to determine the effects of WD on visceral adiposity and metabolic/arterial function in UM-HET3 mice, an outbred, genetically diverse strain of mice. Male and female UM-HET3 mice underwent normal chow (NC) or WD for 12 weeks. Body mass and visceral adiposity were higher in WD compared to NC (P < 0.05). Female WD mice had greater visceral adiposity than male WD mice (P < 0.05). The results of glucose and insulin tolerance tests demonstrated that metabolic function was lower in WD compared to NC mice (P < 0.05). Metabolic dysfunction in WD as was driven by male mice, as metabolic function in female WD mice was unchanged (P > 0.05). Systolic blood pressure (BP) and aortic stiffness were increased in WD after 2 weeks compared to baseline and continued to increase through week 12 (P < 0.05). Systolic BP and aortic stiffness were higher from weeks 2-12 in WD compared to NC (P < 0.05). Aortic collagen content was higher in WD compared to NC (P < 0.05). Carotid artery endothelium-dependent dilation was lower in WD compared to NC (P < 0.05). These data suggest sex-related differences in visceral adiposity and metabolic dysfunction in response to WD. Despite this, arterial dysfunction was similar in male and female WD mice, indicating this model may provide unique translational insight into similar sex-related observations in humans that consume WD.

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.

A Bayesian model selection approach to mediation analysis

Genetic studies often seek to establish a causal chain of events originating from genetic variation through to molecular and clinical phenotypes. When multiple phenotypes share a common genetic association, one phenotype may act as an intermediate for the genetic effects on the other. Alternatively, the phenotypes may be causally unrelated but share genetic loci. Mediation analysis represents a class of causal inference approaches used to determine which of these scenarios is most plausible. We have developed a general approach to mediation analysis based on Bayesian model selection and have implemented it in an R package, bmediatR. Bayesian model selection provides a flexible framework that can be tailored to different analyses. Our approach can incorporate prior information about the likelihood of models and the strength of causal effects. It can also accommodate multiple genetic variants or multi-state haplotypes. Our approach reports posterior probabilities that can be useful in interpreting uncertainty among competing models. We compared bmediatR with other popular methods, including the Sobel test, Mendelian randomization, and Bayesian network analysis using simulated data. We found that bmediatR performed as well or better than these alternatives in most scenarios. We applied bmediatR to proteome data from Diversity Outbred (DO) mice, a multi-parent population, and demonstrate the power of mediation with multi-state haplotypes. We also applied bmediatR to data from human cell lines to identify transcripts that are mediated through or are expressed independently from local chromatin accessibility. We demonstrate that Bayesian model selection provides a powerful and versatile approach to identify causal relationships in genetic studies using model organism or human data.

It's What and When You Eat: An Overview of Transcriptional and Epigenetic Responses to Dietary Perturbations in Pancreatic Islets

Our ever-changing modern environment is a significant contributor to the increased prevalence of many chronic diseases, and particularly, type 2 diabetes mellitus (T2DM). Although the modern era has ushered in numerous changes to our daily living conditions, changes in "what" and "when" we eat appear to disproportionately fuel the rise of T2DM. The pancreatic islet is a key biological controller of an organism's glucose homeostasis and thus plays an outsized role to coordinate the response to environmental factors to preserve euglycemia through a delicate balance of endocrine outputs. Both successful and failed adaptation to dynamic environmental stimuli has been postulated to occur due to changes in the transcriptional and epigenetic regulation of pathways associated with islet secretory function and survival. Therefore, in this review we examined and evaluated the current evidence elucidating the key epigenetic mechanisms and transcriptional programs underlying the islet's coordinated response to the interaction between the timing and the composition of dietary nutrients common to modern lifestyles. With the explosion of next generation sequencing, along with the development of novel informatic and -omic approaches, future work will continue to unravel the environmental-epigenetic relationship in islet biology with the goal of identifying transcriptional and epigenetic targets associated with islet perturbations in T2DM.

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.

Host-pathogen genetic interactions underlie tuberculosis susceptibility in genetically diverse mice

The outcome of an encounter with Mycobacterium tuberculosis (Mtb) depends on the pathogen’s ability to adapt to the variable immune pressures exerted by the host. Understanding this interplay has proven difficult, largely because experimentally tractable animal models do not recapitulate the heterogeneity of tuberculosis disease. We leveraged the genetically diverse Collaborative Cross (CC) mouse panel in conjunction with a library of Mtb mutants to create a resource for associating bacterial genetic requirements with host genetics and immunity. We report that CC strains vary dramatically in their susceptibility to infection and produce qualitatively distinct immune states. Global analysis of Mtb transposon mutant fitness (TnSeq) across the CC panel revealed that many virulence pathways are only required in specific host microenvironments, identifying a large fraction of the pathogen’s genome that has been maintained to ensure fitness in a diverse population. Both immunological and bacterial traits can be associated with genetic variants distributed across the mouse genome, making the CC a unique population for identifying specific host-pathogen genetic interactions that influence pathogenesis.

β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.

Regulation of protein abundance in genetically diverse mouse populations

Genetically diverse mouse populations are powerful tools for characterizing the regulation of the proteome and its relationship to whole-organism phenotypes. We used mass spectrometry to profile and quantify the abundance of 6,798 proteins in liver tissue from mice of both sexes across 58 Collaborative Cross (CC) inbred strains. We previously collected liver proteomics data from the related Diversity Outbred (DO) mice and their founder strains. We show concordance across the proteomics datasets despite being generated from separate experiments, allowing comparative analysis. We map protein abundance quantitative trait loci (pQTLs), identifying 1,087 local and 285 distal in the CC mice and 1,706 local and 414 distal in the DO mice. We find that regulatory effects on individual proteins are conserved across the mouse populations, in particular for local genetic variation and sex differences. In comparison, proteins that form complexes are often co-regulated, displaying varying genetic architectures, and overall show lower heritability and map fewer pQTLs. We have made this resource publicly available to enable quantitative analyses of the regulation of the proteome.

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.

Human-relevant mechanisms and risk factors for TAK-875-Induced liver injury identified via a gene pathway-based approach in Collaborative Cross mice

Development of TAK-875 was discontinued when a small number of serious drug-induced liver injury (DILI) cases were observed in Phase 3 clinical trials. Subsequent studies have identified hepatocellular oxidative stress, mitochondrial dysfunction, altered bile acid homeostasis, and immune response as mechanisms of TAK-875 DILI and the contribution of genetic risk factors in oxidative response and mitochondrial pathways to the toxicity susceptibility observed in patients. We tested the hypothesis that a novel preclinical approach based on gene pathway analysis in the livers of Collaborative Cross mice could be used to identify human-relevant mechanisms of toxicity and genetic risk factors at the level of the hepatocyte as reported in a human genome-wide association study. Eight (8) male mice (4 matched pairs) from each of 45 Collaborative Cross lines were treated with a single oral (gavage) dose of either vehicle or 600 mg/kg TAK-875. As expected, liver injury was not detected histologically and few changes in plasma biomarkers of hepatotoxicity were observed. However, gene expression profiling in the liver identified hundreds of transcripts responsive to TAK-875 treatment across all strains reflecting alterations in immune response and bile acid homeostasis and the interaction of treatment and strain reflecting oxidative stress and mitochondrial dysfunction. Fold-change expression values were then used to develop pathway-based phenotypes for genetic mapping which identified candidate risk factor genes for TAK-875 toxicity susceptibility at the level of the hepatocyte. Taken together, these findings support our hypothesis that a gene pathway-based approach using Collaborative Cross mice could inform sensitive strains, human-relevant mechanisms of toxicity, and genetic risk factors for TAK-875 DILI. This novel preclinical approach may be helpful in understanding, predicting, and ultimately preventing clinical DILI for other drugs.

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/ .

Genetic architecture modulates diet-induced hepatic mRNA and miRNA expression profiles in Diversity Outbred mice

Genetic approaches in model organisms have consistently demonstrated that molecular traits such as gene expression are under genetic regulation, similar to clinical traits. The resulting expression quantitative trait loci (eQTL) have revolutionized our understanding of genetic regulation and identified numerous candidate genes for clinically relevant traits. More recently, these analyses have been extended to other molecular traits such as protein abundance, metabolite levels, and miRNA expression. Here, we performed global hepatic eQTL and microRNA expression quantitative trait loci (mirQTL) analysis in a population of Diversity Outbred mice fed two different diets. We identified several key features of eQTL and mirQTL, namely differences in the mode of genetic regulation (cis or trans) between mRNA and miRNA. Approximately 50% of mirQTL are regulated by a trans-acting factor, compared to ∼25% of eQTL. We note differences in the heritability of mRNA and miRNA expression and variance explained by each eQTL or mirQTL. In general, cis-acting variants affecting mRNA or miRNA expression explain more phenotypic variance than trans-acting variants. Finally, we investigated the effect of diet on the genetic architecture of eQTL and mirQTL, highlighting the critical effects of environment on both eQTL and mirQTL. Overall, these data underscore the complex genetic regulation of two well-characterized RNA classes (mRNA and miRNA) that have critical roles in the regulation of clinical traits and disease susceptibility

Hepatic transcriptional profile reveals the role of diet and genetic backgrounds on metabolic traits in female progenitor strains of the Collaborative Cross

Mice have provided critical mechanistic understandings of clinical traits underlying metabolic syndrome (MetSyn) and susceptibility to MetSyn in mice is known to vary among inbred strains. We investigated the diet- and strain-dependent effects on metabolic traits in the eight Collaborative Cross (CC) founder strains (A/J, C57BL/6J, 129S1/SvImJ, NOD/ShiLtJ, NZO/HILtJ, CAST/EiJ, PWK/PhJ, and WSB/EiJ). Liver transcriptomics analysis showed that both atherogenic diet and host genetics have profound effects on the liver transcriptome, which may be related to differences in metabolic traits observed between strains. We found strain differences in circulating trimethylamine N-oxide (TMAO) concentration and liver triglyceride content, both of which are traits associated with metabolic diseases. Using a network approach, we identified a module of transcripts associated with TMAO and liver triglyceride content, which was enriched in functional pathways. Interrogation of the module related to metabolic traits identified NADPH oxidase 4 (Nox4), a gene for a key enzyme in the production of reactive oxygen species, which showed a strong association with plasma TMAO and liver triglyceride. Interestingly, Nox4 was identified as the highest expressed in the C57BL/6J and NZO/HILtJ strains and the lowest expressed in the CAST/EiJ strain. Based on these results, we suggest that there may be genetic variation in the contribution of Nox4 to the regulation of plasma TMAO and liver triglyceride content. In summary, we show that liver transcriptomic analysis identified diet- or strain-specific pathways for metabolic traits in the Collaborative Cross (CC) founder strains.

Synergizing Mouse and Human Studies to Understand the Heterogeneity of Obesity

Obesity is routinely considered as a single disease state, which drives a "one-size-fits-all" approach to treatment. We recently convened the first annual University of North Carolina Interdisciplinary Nutrition Sciences Symposium to discuss the heterogeneity of obesity and the need for translational science to advance understanding of this heterogeneity. The symposium aimed to advance scientific rigor in translational studies from animal to human models with the goal of identifying underlying mechanisms and treatments. In this review, we discuss fundamental gaps in knowledge of the heterogeneity of obesity ranging from cellular to population perspectives. We also advocate approaches to overcoming limitations in the field. Examples include the use of contemporary mouse genetic reference population models such as the Collaborative Cross and Diversity Outbred mice that effectively model human genetic diversity and the use of translational models that integrate -omics and computational approaches from pre-clinical to clinical models of obesity. Finally, we suggest best scientific practices to ensure strong rigor that will allow investigators to delineate the sources of heterogeneity in the population with obesity. Collectively, we propose that it is critical to think of obesity as a heterogeneous disease with complex mechanisms and etiologies, requiring unique prevention and treatment strategies tailored to the individual.

Mouse Genetic Reference Populations: Cellular Platforms for Integrative Systems Genetics

Interrogation of disease-relevant cellular and molecular traits exhibited by genetically diverse cell populations enables in vitro systems genetics approaches for uncovering the basic properties of cellular function and identity. Primary cells, stem cells, and organoids derived from genetically diverse mouse strains, such as Collaborative Cross and Diversity Outbred populations, offer the opportunity for parallel in vitro/in vivo screening. These panels provide genetic resolution for variant discovery and functional characterization, as well as disease modeling and in vivo validation capabilities. Here we review mouse cellular systems genetics approaches for characterizing the influence of genetic variation on signaling networks and phenotypic diversity, and we discuss approaches for data integration and cross-species validation.

Diet-induced obesity in genetically diverse collaborative cross mouse founder strains reveals diverse phenotype response and amelioration by quercetin treatment in 129S1/SvImJ, PWK/EiJ, CAST/PhJ, and WSB/EiJ mice

Significant evidence suggests protective effects of flavonoids against obesity in animal models, but these often do not translate to humans. One explanation for this disconnect is use of a few mouse strains (notably C57BL/6 J) in obesity studies. Obesity is a multifactorial disease. The underlying causes are not fully replicated by the high-fat C57BL/6 J model, despite phenotypic similarities. Furthermore, the impact of genetic factors on the activities of flavonoids is unknown. This study was designed to explore how diverse mouse strains respond to diet-induced obesity when fed a representative flavonoid. A subset of Collaborative Cross founder strains (males and females) were placed on dietary treatments (low-fat, high-fat, high-fat with quercetin, high-fat with quercetin and antibiotics) longitudinally. Diverse responses were observed across strains and sexes. Quercetin appeared to moderately blunt weight gain in male C57 and both sexes of 129S1/SvImJ mice, and slightly increased weight gain in female C57 mice. Surprisingly, quercetin dramatically blunted weight gain in male, but not female, PWK/PhJ mice. For female mice, quercetin blunted weight gain (relative to the high-fat phase) in CAST/PhJ, PWK/EiJ and WSB/EiJ mice compared to C57. Antibiotics did not generally result in loss of protective effects of quercetin. This highlights complex interactions between genetic factors, sex, obesity stimuli, and flavonoid intake, and the need to move away from single inbred mouse models to enhance translatability to diverse humans. These data justify use of genetically diverse Collaborative Cross and Diversity Outbred models which are emerging as invaluable tools in the field of personalized nutrition.

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 Architecture Modulates Diet-Induced Hepatic mRNA and miRNA Expression Profiles in Diversity Outbred Mice

Genetic approaches in model organisms have consistently demonstrated that molecular traits such as gene expression are under genetic regulation, similar to clinical traits. The resulting expression quantitative trait loci (eQTL) have revolutionized our understanding of genetic regulation and identified numerous candidate genes for clinically relevant traits. More recently, these analyses have been extended to other molecular traits such as protein abundance, metabolite levels, and miRNA expression. Here, we performed global hepatic eQTL and microRNA expression quantitative trait loci (mirQTL) analysis in a population of Diversity Outbred mice fed two different diets. We identified several key features of eQTL and mirQTL, namely differences in the mode of genetic regulation (cis or trans) between mRNA and miRNA. Approximately 50% of mirQTL are regulated by a trans-acting factor, compared to ∼25% of eQTL. We note differences in the heritability of mRNA and miRNA expression and variance explained by each eQTL or mirQTL. In general, cis-acting variants affecting mRNA or miRNA expression explain more phenotypic variance than trans-acting variants. Lastly, we investigated the effect of diet on the genetic architecture of eQTL and mirQTL, highlighting the critical effects of environment on both eQTL and mirQTL. Overall, these data underscore the complex genetic regulation of two well-characterized RNA classes (mRNA and miRNA) that have critical roles in the regulation of clinical traits and disease susceptibility.

(Epi)genomic heterogeneity of pancreatic islet function and failure in type 2 diabetes

Background

Pancreatic Islets of Langerhans are heterogeneous tissues consisting of multiple endocrine cell types that carry out distinct yet coordinated roles to regulate blood glucose homeostasis. Islet dysfunction and specifically failure of the beta cells to secrete adequate insulin are known precursors to type 2 diabetes (T2D) onset. However, the exact genetic, (epi)genomic, and environmental mechanisms that contribute to islet failure, and ultimately to T2D pathogenesis, require further elucidation.

Scope of review

This review summarizes efforts and advances in dissection of the complex genetic underpinnings of islet function and resilience in T2D pathogenesis. In this review, we will highlight results of the latest T2D genome-wide association study (GWAS) and discuss how these data are being combined with clinical measures in patients to uncover putative T2D subtypes and with functional (epi)genomic studies in islets to understand the genetic programming of islet cell identity, function, and adaptation. Finally, we discuss new and important opportunities to address major knowledge gaps in our understanding of islet (dys)function in T2D risk and progression.

Major conclusions

Genetic variation exerts clear effects on the islet epigenome, regulatory element usage, and gene expression. Future (epi)genomic comparative analyses between T2D and normal islets should incorporate genetics to distinguish patient-specific from disease-specific differences. Incorporating genotype information into future analyses and studies will also enable more precise insights into the molecular genetics of islet deficiency and failure in T2D risk, and should ultimately contribute to a stratified view of T2D and more precise treatment strategies. Islet cellular heterogeneity continues to remain a challenge for understanding the associations between islet failure and T2D development. Further efforts to obtain purified islet cell type populations and determine the specific genetic and environmental effects on each will help address this. Beyond observation of islets at steady state conditions, more research of islet stress and stimulation responses are needed to understand the transition of these tissues from a healthy to diseased state. Together, focusing on these objectives will provide more opportunities to prevent, treat, and manage T2D.

Mouse Systems Genetics as a Prelude to Precision Medicine

Mouse models have been instrumental in understanding human disease biology and proposing possible new treatments. The precise control of the environment and genetic composition of mice allows more rigorous observations, but limits the generalizability and translatability of the results into human applications. In the era of precision medicine, strategies using mouse models have to be revisited to effectively emulate human populations. Systems genetics is one promising paradigm that may promote the transition to novel precision medicine strategies. Here, we review the state-of-the-art resources and discuss how mouse systems genetics helps to understand human diseases and to advance the development of precision medicine, with an emphasis on the existing resources and strategies.

Integrative QTL analysis of gene expression and chromatin accessibility identifies multi-tissue patterns of genetic regulation

Gene transcription profiles across tissues are largely defined by the activity of regulatory elements, most of which correspond to regions of accessible chromatin. Regulatory element activity is in turn modulated by genetic variation, resulting in variable transcription rates across individuals. The interplay of these factors, however, is poorly understood. Here we characterize expression and chromatin state dynamics across three tissues-liver, lung, and kidney-in 47 strains of the Collaborative Cross (CC) mouse population, examining the regulation of these dynamics by expression quantitative trait loci (eQTL) and chromatin QTL (cQTL). QTL whose allelic effects were consistent across tissues were detected for 1,101 genes and 133 chromatin regions. Also detected were eQTL and cQTL whose allelic effects differed across tissues, including local-eQTL for Pik3c2g detected in all three tissues but with distinct allelic effects. Leveraging overlapping measurements of gene expression and chromatin accessibility on the same mice from multiple tissues, we used mediation analysis to identify chromatin and gene expression intermediates of eQTL effects. Based on QTL and mediation analyses over multiple tissues, we propose a causal model for the distal genetic regulation of Akr1e1, a gene involved in glycogen metabolism, through the zinc finger transcription factor Zfp985 and chromatin intermediates. This analysis demonstrates the complexity of transcriptional and chromatin dynamics and their regulation over multiple tissues, as well as the value of the CC and related genetic resource populations for identifying specific regulatory mechanisms within cells and tissues.

Systems genetics applications in metabolism research

The common forms of metabolic diseases are highly complex, involving hundreds of genes, environmental and lifestyle factors, age-related changes, sex differences and gut-microbiome interactions. Systems genetics is a population-based approach to address this complexity. In contrast to commonly used 'reductionist' approaches, such as gain or loss of function, that examine one element at a time, systems genetics uses high-throughput 'omics' technologies to quantitatively assess the many molecular differences among individuals in a population and then to relate these to physiologic functions or disease states. Unlike genome-wide association studies, systems genetics seeks to go beyond the identification of disease-causing genes to understand higher-order interactions at the molecular level. The purpose of this review is to introduce the systems genetics applications in the areas of metabolic and cardiovascular disease. Here, we explain how large clinical and omics-level data and databases from both human and animal populations are available to help researchers place genes in the context of pathways and networks and formulate hypotheses that can then be experimentally examined. We provide lists of such databases and examples of the integration of reductionist and systems genetics data. Among the important applications emerging is the development of improved nutritional and pharmacological strategies to address the rise of metabolic diseases.

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.

Testing Pleiotropy vs. Separate QTL in Multiparental Populations

The high mapping resolution of multiparental populations, combined with technology to measure tens of thousands of phenotypes, presents a need for quantitative methods to enhance understanding of the genetic architecture of complex traits. When multiple traits map to a common genomic region, knowledge of the number of distinct loci provides important insight into the underlying mechanism and can assist planning for subsequent experiments. We extend the method of Jiang and Zeng (1995), for testing pleiotropy with a pair of traits, to the case of more than two alleles. We also incorporate polygenic random effects to account for population structure. We use a parametric bootstrap to determine statistical significance. We apply our methods to a behavioral genetics data set from Diversity Outbred mice. Our methods have been incorporated into the R package qtl2pleio.

Peptidome analysis reveals the involvement of endogenous peptides in mouse pancreatic dysfunction with aging

Type 2 diabetes (T2D) is a glucose regulation disorder that has significantly enhanced mortality and the global disease burden. The prevalence of T2D has increased worldwide and is higher in the elderly. The function of pancreatic islets decreases with age, which is one important reason for the occurrence of diabetes in the elderly. Recently, peptidome analysis has attracted attention. However, the role of age-related peptides in pancreatic dysfunction has not been investigated extensively. Here, we conducted a comparison of endogenous peptides between pancreas from adult and aging mice by liquid chromatography tandem mass spectrometry (LC-MS/MS). A total of 2,089 peptides originating from 1,280 protein precursors were identified, of which 232 were upregulated and 183 were downregulated in the aging mice (fold change ≥ 2 and p < 0.05), suggesting that the expression of pancreatic peptides in mice varied with age. The molecular weight of most peptides was <3.0 kDa, and the isoelectric point distribution had a bimodal characteristic. Further analysis of cleavage site patterns indicated that proteases cleaved pancreatic proteins according to their rules. Moreover, Gene Ontology and pathway analyses showed that the differentially expressed peptides potentially had specific effects on pancreatic dysfunction. Some differential peptides were located within the domains of precursor proteins that were closely associated with the development of diabetes. We believe that our research may advance the current understanding of pancreas-derived peptides and that certain peptides may be involved in the etiology of diabetes.