Transcriptome-wide analyses of adipose tissue in outbred rats reveal genetic regulatory mechanisms relevant for human obesity.

Transcriptomic analysis in metabolically active tissues allows a systems genetics approach to identify causal genes and networks involved in metabolic disease. Outbred heterogeneous stock (HS) rats are used for genetic mapping of complex traits, but to-date, a systems genetics analysis of metabolic tissues has not been done. We investigated whether adiposity-associated genes and gene coexpression networks in outbred heterogeneous stock (HS) rats overlap those found in humans. We analyzed RNAseq data from adipose tissue of 415 male HS rats, correlated these transcripts with body weight (BW) and compared transcriptome signatures to two human cohorts: the "African American Genetics of Metabolism and Expression" and "Metabolic Syndrome in Men." We used weighted gene coexpression network analysis to identify adiposity-associated gene networks and mediation analysis to identify genes under genetic control whose expression drives adiposity. We identified 554 orthologous "consensus genes" whose expression correlates with BW in the rat and with body mass index (BMI) in both human cohorts. Consensus genes fell within eight coexpressed networks and were enriched for genes involved in immune system function, cell growth, extracellular matrix organization, and lipid metabolic processes. We identified 19 consensus genes for which genetic variation may influence BW via their expression, including those involved in lipolysis (e.g., Hcar1), inflammation (e.g., Rgs1), adipogenesis (e.g., Tmem120b), or no previously known role in obesity (e.g., St14 and Ms4a6a). Strong concordance between HS rat and human BW/BMI associated transcripts demonstrates translational utility of the rat model, while identification of novel genes expands our knowledge of the genetics underlying obesity.

Tpcn2 knockout mice have improved insulin sensitivity and are protected against high-fat diet-induced weight gain.

Type 2 diabetes is a complex disorder affected by multiple genes and the environment. Our laboratory has shown that in response to a glucose challenge, two-pore channel 2 ( Tpcn2) knockout mice exhibit a decreased insulin response but normal glucose clearance, suggesting they have improved insulin sensitivity compared with wild-type mice. We tested the hypothesis that improved insulin sensitivity in Tpcn2 knockout mice would protect against the negative effects of a high fat diet. Male and female Tpcn2 knockout (KO), heterozygous (Het), and wild-type (WT) mice were fed a low-fat (LF) or high-fat (HF) diet for 24 wk. HF diet significantly increases body weight in WT mice relative to those on the LF diet; this HF diet-induced increase in body weight is blunted in the Het and KO mice. Despite the protection against diet-induced weight gain, however, Tpcn2 KO mice are not protected against HF-diet-induced changes in glucose or insulin area under the curve during glucose tolerance tests in female mice, while HF diet has no significant effect on glucose tolerance in the male mice, regardless of genotype. Glucose disappearance during an insulin tolerance test is augmented in male KO mice, consistent with our previous findings suggesting enhanced insulin sensitivity in these mice. Male KO mice exhibit increased fasting plasma total cholesterol and triglyceride concentrations relative to WT mice on the LF diet, but this difference disappears in HF diet-fed mice where there is increased cholesterol and triglycerides across all genotypes. These data demonstrate that knockout of Tpcn2 may increase insulin action in male, but not female, mice. In addition, both male and female KO mice are protected against diet-induced weight gain, but this protection is likely independent from glucose tolerance, insulin sensitivity, and plasma lipid levels.

Genomic Loci Influencing Cue-Reactivity in Heterogeneous Stock Rats.

Addiction vulnerability is associated with the tendency to attribute incentive salience to reward predictive cues; both addiction and the attribution of incentive salience are influenced by environmental and genetic factors. To characterize the genetic contributions to incentive salience attribution, we performed a genome-wide association study (GWAS) in a cohort of 1,645 genetically diverse heterogeneous stock (HS) rats. We tested HS rats in a Pavlovian conditioned approach task, in which we characterized the individual responses to food-associated stimuli ("cues"). Rats exhibited either cue-directed "sign-tracking" behavior or food-cup directed "goal-tracking" behavior. We then used the conditioned reinforcement procedure to determine whether rats would perform a novel operant response for unrewarded presentations of the cue. We found that these measures were moderately heritable (SNP heritability, h2 = .189-.215). GWAS identified 14 quantitative trait loci (QTLs) for 11 of the 12 traits we examined. Interval sizes of these QTLs varied widely. 7 traits shared a QTL on chromosome 1 that contained a few genes (e.g. Tenm4, Mir708) that have been associated with substance use disorders and other mental health traits in humans. Other candidate genes (e.g. Wnt11, Pak1) in this region had coding variants and expression-QTLs in mesocorticolimbic regions of the brain. We also conducted a Phenome-Wide Association Study (PheWAS) on other behavioral measures in HS rats and found that regions containing QTLs on chromosome 1 were also associated with nicotine self-administration in a separate cohort of HS rats. These results provide a starting point for the molecular genetic dissection of incentive salience and provide further support for a relationship between attribution of incentive salience and drug abuse-related traits.

Genetic background in the rat affects endocrine and metabolic outcomes of bisphenol F exposure.

Environmental bisphenol compounds like bisphenol F (BPF) are endocrine-disrupting chemicals (EDCs) affecting adipose and classical endocrine systems. Genetic factors that influence EDC exposure outcomes are poorly understood and are unaccounted variables that may contribute to the large range of reported outcomes in the human population. We previously demonstrated that BPF exposure increased body growth and adiposity in male N/NIH heterogeneous stock (HS) rats, a genetically heterogeneous outbred population. We hypothesize that the founder strains of the HS rat exhibit EDC effects that were strain- and sex-dependent. Weanling littermate pairs of male and female ACI, BN, BUF, F344, M520, and WKY rats randomly received either vehicle (0.1% EtOH) or 1.125 mg BPF/l in 0.1% EtOH for 10 weeks in drinking water. Body weight and fluid intake were measured weekly, metabolic parameters were assessed, and blood and tissues were collected. BPF increased thyroid weight in ACI males, thymus and kidney weight in BUF females, adrenal weight in WKY males, and possibly increased pituitary weight in BN males. BUF females also developed a disruption in activity and metabolic rate with BPF exposure. These sex- and strain-specific exposure outcomes illustrate that HS rat founders possess diverse bisphenol-exposure risk alleles and suggest that BPF exposure may intensify inherent organ system dysfunction existing in the HS rat founders. We propose that the HS rat will be an invaluable model for dissecting gene EDC interactions on health.

The Power of the Heterogeneous Stock Rat Founder Strains in Modeling Metabolic Disease.

Metabolic diseases are a host of complex conditions, including obesity, diabetes mellitus, and metabolic syndrome. Endocrine control systems (eg, adrenals, thyroid, gonads) are causally linked to metabolic health outcomes. N/NIH Heterogeneous Stock (HS) rats are a genetically heterogeneous outbred population developed for genetic studies of complex traits. Genetic mapping studies in adult HS rats identified loci associated with cardiometabolic risks, such as glucose intolerance, insulin resistance, and increased body mass index. This study determined underappreciated metabolic health traits and the associated endocrine glands within available substrains of the HS rat founders. We hypothesize that the genetic diversity of the HS rat founder strains causes a range of endocrine health conditions contributing to the diversity of cardiometabolic disease risks. ACI/EurMcwi, BN/NHsdMcwi, BUF/MnaMcwi, F344/StmMcwi, M520/NRrrcMcwi, and WKY/NCrl rats of both sexes were studied from birth until 13 weeks of age. Birth weight was recorded, body weight was measured weekly, metabolic characteristics were assessed, and blood and tissues were collected. Our data show wide variation in endocrine traits and metabolic health states in ACI, BN, BUF, F344, M520, and WKY rat strains. This is the first report to compare birth weight, resting metabolic rate, endocrine gland weight, hypothalamic-pituitary-thyroid axis hormones, and brown adipose tissue weight in these rat strains. Importantly, this work unveils new potential for the HS rat population to model early life adversity and adrenal and thyroid pathophysiology. The HS population likely inherited risk alleles for these strain-specific traits, making the HS rat a powerful model to investigate interventions on endocrine and metabolic health.

Genetic Mapping of Multiple Traits Identifies Novel Genes for Adiposity, Lipids, and Insulin Secretory Capacity in Outbred Rats.

Despite the successes of human genome-wide association studies, the causal genes underlying most metabolic traits remain unclear. We used outbred heterogeneous stock (HS) rats, coupled with expression data and mediation analysis, to identify quantitative trait loci (QTLs) and candidate gene mediators for adiposity, glucose tolerance, serum lipids, and other metabolic traits. Physiological traits were measured in 1,519 male HS rats, with liver and adipose transcriptomes measured in >410 rats. Genotypes were imputed from low-coverage whole-genome sequencing. Linear mixed models were used to detect physiological and expression QTLs (pQTLs and eQTLs, respectively), using both single nucleotide polymorphism (SNP)- and haplotype-based models for pQTL mapping. Genes with cis-eQTLs that overlapped pQTLs were assessed as causal candidates through mediation analysis. We identified 14 SNP-based pQTLs and 19 haplotype-based pQTLs, of which 10 were in common. Using mediation, we identified the following genes as candidate mediators of pQTLs: Grk5 for fat pad weight and serum triglyceride pQTLs on Chr1, Krtcap3 for fat pad weight and serum triglyceride pQTLs on Chr6, Ilrun for a fat pad weight pQTL on Chr20, and Rfx6 for a whole pancreatic insulin content pQTL on Chr20. Furthermore, we verified Grk5 and Ktrcap3 using gene knockdown/out models, thereby shedding light on novel regulators of obesity.

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.

An exponential increase in QTL detection with an increased sample size.

Power analyses are often used to determine the number of animals required for a genome-wide association study (GWAS). These analyses are typically intended to estimate the sample size needed for at least 1 locus to exceed a genome-wide significance threshold. A related question that is less commonly considered is the number of significant loci that will be discovered with a given sample size. We used simulations based on a real data set that consisted of 3,173 male and female adult N/NIH heterogeneous stock rats to explore the relationship between sample size and the number of significant loci discovered. Our simulations examined the number of loci identified in subsamples of the full data set. The subsampling analysis was conducted for 4 traits with low (0.15 ± 0.03), medium (0.31 ± 0.03 and 0.36 ± 0.03), and high (0.46 ± 0.03) SNP-based heritabilities. For each trait, we subsampled the data 100 times at different sample sizes (500, 1,000, 1,500, 2,000, and 2,500). We observed an exponential increase in the number of significant loci with larger sample sizes. Our results are consistent with similar observations in human GWAS and imply that future rodent GWAS should use sample sizes that are significantly larger than those needed to obtain a single significant result.

Fine-mapping diabetes-related traits, including insulin resistance, in heterogeneous stock rats.

Type 2 diabetes (T2D) is a disease of relative insulin deficiency resulting from both insulin resistance and beta cell failure. We have previously used heterogeneous stock (HS) rats to fine-map a locus for glucose tolerance. We show here that glucose intolerance in the founder strains of the HS colony is mediated by different mechanisms: insulin resistance in WKY and an insulin secretion defect in ACI, and we demonstrate a high degree of variability for measures of insulin resistance and insulin secretion in HS rats. As such, our goal was to use HS rats to fine-map several diabetes-related traits within a region on rat chromosome 1. We measured blood glucose and plasma insulin levels after a glucose tolerance test in 782 male HS rats. Using 97 SSLP markers, we genotyped a 68 Mb region on rat chromosome 1 previously implicated in glucose and insulin regulation. We used linkage disequilibrium mapping by mixed model regression with inferred descent to identify a region from 198.85 to 205.9 that contains one or more quantitative trait loci (QTL) for fasting insulin and a measure of insulin resistance, the quantitative insulin sensitivity check index. This region also encompasses loci identified for fasting glucose and Insulin_AUC (area under the curve). A separate <3 Mb QTL was identified for body weight. Using a novel penalized regression method we then estimated effects of alternative haplotype pairings under each locus. These studies highlight the utility of HS rats for fine-mapping genetic loci involved in the underlying causes of T2D.

Body Composition and Metabolic Changes in a Lyon Hypertensive Congenic Rat and Identification of Ercc6l2 as a Positional Candidate Gene.

Central obesity is genetically complex, and its exponential increase in the last decades have made it a critical public health issue. The Lyon Hypertensive (LH) rat is a well-characterized hypertensive model that also exhibits spontaneous and profound differences in body weight and adiposity, relative to its metabolically healthy control, the Lyon Normotensive (LN) rat. The mechanisms underlying the body weight differences between these strains are not well-understood, thus a congenic model (LH17LNa) was developed where a portion of the proximal arm of LN chromosome 17 is introgressed on the LH genomic background to assess the contribution of LN alleles on obesity features. Male and female LH17LNa rats were studied, but male congenics did not significantly differ from LH in this study. Female LH17LNa rats exhibited decreases in total body growth, as well as major alterations to their body composition and adiposity. The LH17LNa female rats also showed decreases in metabolic rate, and a reduction in food intake. The increased adiposity in the female LH17LNa rats was specific to abdominal white adipose tissue, and this phenomenon was further explained by significant hypertrophy in those adipocytes, with no evidence of adipocyte hyperplasia. Sequencing of the parental strains identified a novel frameshift mutation in the candidate gene Ercc6l2, which is involved in transcription-coupled DNA repair, and is implicated in premature aging. The discovery of the significance of Ercc6l2 in the context of female-specific adipocyte biology could represent a novel role of DNA repair failure syndromes in obesity pathogenesis.

Genome-Wide Association Study on Three Behaviors Tested in an Open Field in Heterogeneous Stock Rats Identifies Multiple Loci Implicated in Psychiatric Disorders.

Many personality traits are influenced by genetic factors. Rodents models provide an efficient system for analyzing genetic contribution to these traits. Using 1,246 adolescent heterogeneous stock (HS) male and female rats, we conducted a genome-wide association study (GWAS) of behaviors measured in an open field, including locomotion, novel object interaction, and social interaction. We identified 30 genome-wide significant quantitative trait loci (QTL). Using multiple criteria, including the presence of high impact genomic variants and co-localization of cis-eQTL, we identified 17 candidate genes (Adarb2, Ankrd26, Cacna1c, Cacng4, Clock, Ctu2, Cyp26b1, Dnah9, Gda, Grxcr1, Eva1a, Fam114a1, Kcnj9, Mlf2, Rab27b, Sec11a, and Ube2h) for these traits. Many of these genes have been implicated by human GWAS of various psychiatric or drug abuse related traits. In addition, there are other candidate genes that likely represent novel findings that can be the catalyst for future molecular and genetic insights into human psychiatric diseases. Together, these findings provide strong support for the use of the HS population to study psychiatric disorders.

Bisphenol F Exposure in Adolescent Heterogeneous Stock Rats Affects Growth and Adiposity.

Bisphenol F (BPF) is increasingly substituting bisphenol A in manufacturing polycarbonates and consumer products. The cardiometabolic effects of BPF in either humans or model organisms are not clear, and no studies to date have investigated the role of genetic background on susceptibility to BPF-induced cardiometabolic traits. The primary goal of this project was to determine if BPF exposure influences growth and adiposity in male N:NIH heterogeneous stock (HS) rats, a genetically heterogeneous population. Littermate pairs of male HS rats were randomly exposed to either vehicle (0.1% ethanol) or 1.125 µg/ml BPF in 0.1% ethanol for 5 weeks in drinking water starting at 3 weeks-of-age. Water consumption and body weight was measured weekly, body composition was determined using nuclear magnetic resonance, urine and feces were collected in metabolic cages, and blood and tissues were collected at the end of the study. BPF-exposed rats showed significantly increased body growth and abdominal adiposity, risk factors for cardiometabolic disease. Urine output was increased in BPF-exposed rats, driving a trend in increased creatinine clearance. We also report the first relationship between a bisphenol metabolizing enzyme and a bisphenol-induced phenotype. Preliminary heritability estimates of significant phenotypes suggest that BPF exposure may alter trait variation. These findings support BPF exposure as a cardiometabolic disease risk factor and indicate that the HS rat will be a useful model for dissecting gene by BPF interactions on metabolic health.

Sensitivity to food and cocaine cues are independent traits in a large sample of heterogeneous stock rats.

Sensitivity to cocaine and its associated stimuli ("cues") are important factors in the development and maintenance of addiction. Rodent studies suggest that this sensitivity is related, in part, to the propensity to attribute incentive salience to food cues, which, in turn, contributes to the maintenance of cocaine self-administration, and cue-induced relapse of drug-seeking. Whereas each of these traits has established links to drug use, the relatedness between the individual traits themselves has not been well characterized in preclinical models. To this end, the propensity to attribute incentive salience to a food cue was first assessed in two distinct cohorts of 2716 outbred heterogeneous stock rats (HS; formerly N:NIH). We then determined whether each cohort was associated with performance in one of two paradigms (cocaine conditioned cue preference and cocaine contextual conditioning). These measure the unconditioned locomotor effects of cocaine, as well as conditioned approach and the locomotor response to a cocaine-paired floor or context. There was large individual variability and sex differences among all traits, but they were largely independent of one another in both males and females. These findings suggest that these traits may contribute to drug-use via independent underlying neuropsychological processes.

Sept8/SEPTIN8 involvement in cellular structure and kidney damage is identified by genetic mapping and a novel human tubule hypoxic model.

Chronic kidney disease (CKD), which can ultimately progress to kidney failure, is influenced by genetics and the environment. Genes identified in human genome wide association studies (GWAS) explain only a small proportion of the heritable variation and lack functional validation, indicating the need for additional model systems. Outbred heterogeneous stock (HS) rats have been used for genetic fine-mapping of complex traits, but have not previously been used for CKD traits. We performed GWAS for urinary protein excretion (UPE) and CKD related serum biochemistries in 245 male HS rats. Quantitative trait loci (QTL) were identified using a linear mixed effect model that tested for association with imputed genotypes. Candidate genes were identified using bioinformatics tools and targeted RNAseq followed by testing in a novel in vitro model of human tubule, hypoxia-induced damage. We identified two QTL for UPE and five for serum biochemistries. Protein modeling identified a missense variant within Septin 8 (Sept8) as a candidate for UPE. Sept8/SEPTIN8 expression increased in HS rats with elevated UPE and tubulointerstitial injury and in the in vitro hypoxia model. SEPTIN8 is detected within proximal tubule cells in human kidney samples and localizes with acetyl-alpha tubulin in the culture system. After hypoxia, SEPTIN8 staining becomes diffuse and appears to relocalize with actin. These data suggest a role of SEPTIN8 in cellular organization and structure in response to environmental stress. This study demonstrates that integration of a rat genetic model with an environmentally induced tubule damage system identifies Sept8/SEPTIN8 and informs novel aspects of the complex gene by environmental interactions contributing to CKD risk.

Fine-mapping a locus for glucose tolerance using heterogeneous stock rats.

Heterogeneous stock (HS) animals provide the ability to map quantitative trait loci at high resolution [<5 Megabase (Mb)] in a relatively short time period. In the current study, we hypothesized that the HS rat colony would be useful for fine-mapping a region on rat chromosome 1 that has previously been implicated in glucose regulation. We administered a glucose tolerance test to 515 HS rats and genotyped these animals with 69 microsatellite markers, spaced an average distance of <1 Mb apart, on a 67 Mb region of rat chromosome 1. Using regression modeling of inferred haplotypes based on a hidden Markov model reconstruction and mixed model analysis in which we accounted for the complex family structure of the HS, we identified one sharp peak within this region. Using positional bootstrapping, we determined the most likely location of this locus is from 205.04 to 207.48 Mb. This work demonstrates the utility of HS rats for fine-mapping complex traits and emphasizes the importance of taking into account family structure when using highly recombinant populations.

Genome-Wide Association Study in 3,173 Outbred Rats Identifies Multiple Loci for Body Weight, Adiposity, and Fasting Glucose.

OBJECTIVE: Obesity is influenced by genetic and environmental factors. Despite the success of human genome-wide association studies, the specific genes that confer obesity remain largely unknown. The objective of this study was to use outbred rats to identify the genetic loci underlying obesity and related morphometric and metabolic traits. METHODS: This study measured obesity-relevant traits, including body weight, body length, BMI, fasting glucose, and retroperitoneal, epididymal, and parametrial fat pad weight in 3,173 male and female adult N/NIH heterogeneous stock (HS) rats across three institutions, providing data for the largest rat genome-wide association study to date. Genetic loci were identified using a linear mixed model to account for the complex family relationships of the HS and using covariates to account for differences among the three phenotyping centers. RESULTS: This study identified 32 independent loci, several of which contained only a single gene (e.g., Epha5, Nrg1, Klhl14) or obvious candidate genes (e.g., Adcy3, Prlhr). There were strong phenotypic and genetic correlations among obesity-related traits, and there was extensive pleiotropy at individual loci. CONCLUSIONS: This study demonstrates the utility of HS rats for investigating the genetics of obesity-related traits across institutions and identify several candidate genes for future functional testing.

Extended regions of suspected mis-assembly in the rat reference genome.

We performed whole-genome sequencing for eight inbred rat strains commonly used in genetic mapping studies. They are the founders of the NIH heterogeneous stock (HS) outbred colony. We provide their sequences and variant calls to the rat genomics community. When analyzing the variant calls we identified regions with unusually high levels of heterozygosity. These regions are consistent across the eight inbred strains, including Brown Norway, which is the basis of the rat reference genome. These regions show higher read depths than other regions in the genome and contain higher rates of apparent tri-allelic variant sites. The evidence suggests that these regions may correspond to duplicated segments that were incorrectly overlaid as a single segment in the reference genome. We provide masks for these regions of suspected mis-assembly as a resource for the community to flag potentially false interpretations of mapping or functional results.

Incentive salience attribution, "sensation-seeking" and "novelty-seeking" are independent traits in a large sample of male and female heterogeneous stock rats.

There are a number of traits that are thought to increase susceptibility to addiction, and some of these are modeled in preclinical studies. For example, "sensation-seeking" is predictive of the initial propensity to take drugs; whereas "novelty-seeking" predicts compulsive drug-seeking behavior. In addition, the propensity to attribute incentive salience to reward cues can predict the propensity to approach drug cues, and reinstatement or relapse, even after relatively brief periods of drug exposure. The question addressed here is the extent to which these three 'vulnerability factors' are related; that is, predictive of one another. Some relationships have been reported in small samples, but here a large sample of 1,598 outbred male and female heterogeneous stock rats were screened for Pavlovian conditioned approach behavior (to obtain an index of incentive salience attribution; 'sign-tracking'), and subsequently tested for sensation-seeking and novelty-seeking. Despite the large N there were no significant correlations between these traits, in either males or females. There were, however, novel relationships between multiple measures of incentive salience attribution and, based on these findings, we generated a new metric that captures "incentive value". Furthermore, there were sex differences on measures of incentive salience attribution and sensation-seeking behavior that were not previously apparent.

Genetic Fine-Mapping and Identification of Candidate Genes and Variants for Adiposity Traits in Outbred Rats.

OBJECTIVE: Obesity is a major risk factor for multiple diseases and is in part heritable, yet the majority of causative genetic variants that drive excessive adiposity remain unknown. Here, outbred heterogeneous stock (HS) rats were used in controlled environmental conditions to fine-map novel genetic modifiers of adiposity. METHODS: Body weight and visceral fat pad weights were measured in male HS rats that were also genotyped genome-wide. Quantitative trait loci (QTL) were identified by genome-wide association of imputed single-nucleotide polymorphism (SNP) genotypes using a linear mixed effect model that accounts for unequal relatedness between the HS rats. Candidate genes were assessed by protein modeling and mediation analysis of expression for coding and noncoding variants, respectively. RESULTS: HS rats exhibited large variation in adiposity traits, which were highly heritable and correlated with metabolic health. Fine-mapping of fat pad weight and body weight revealed three QTL and prioritized five candidate genes. Fat pad weight was associated with missense SNPs in Adcy3 and Prlhr and altered expression of Krtcap3 and Slc30a3, whereas Grid2 was identified as a candidate within the body weight locus. CONCLUSIONS: These data demonstrate the power of HS rats for identification of known and novel heritable mediators of obesity traits.

Identification of a novel gene for diabetic traits in rats, mice, and humans.

The genetic basis of type 2 diabetes remains incompletely defined despite the use of multiple genetic strategies. Multiparental populations such as heterogeneous stocks (HS) facilitate gene discovery by allowing fine mapping to only a few megabases, significantly decreasing the number of potential candidate genes compared to traditional mapping strategies. In the present work, we employed expression and sequence analysis in HS rats (Rattus norvegicus) to identify Tpcn2 as a likely causal gene underlying a 3.1-Mb locus for glucose and insulin levels. Global gene expression analysis on liver identified Tpcn2 as the only gene in the region that is differentially expressed between HS rats with glucose intolerance and those with normal glucose regulation. Tpcn2 also maps as a cis-regulating expression QTL and is negatively correlated with fasting glucose levels. We used founder sequence to identify variants within this region and assessed association between 18 variants and diabetic traits by conducting a mixed-model analysis, accounting for the complex family structure of the HS. We found that two variants were significantly associated with fasting glucose levels, including a nonsynonymous coding variant within Tpcn2. Studies in Tpcn2 knockout mice demonstrated a significant decrease in fasting glucose levels and insulin response to a glucose challenge relative to those in wild-type mice. Finally, we identified variants within Tpcn2 that are associated with fasting insulin in humans. These studies indicate that Tpcn2 is a likely causal gene that may play a role in human diabetes and demonstrate the utility of multiparental populations for positionally cloning genes within complex loci.