Katie Holl

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.

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.

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.

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

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.

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, and 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 heterozygosity. We show that these regions are consistent across the eight inbred strains, including the BN strain, which was the basis of the rat reference genome. These regions show significantly higher read depth than other regions in the genome. The evidence suggests that these regions may contain segmental duplications that are incorrectly overlaid in the reference genome. We provide masks for these suspected regions of mis-assembly as a resource for the community to flag potentially false interpretations of mapping results or functional data.

Genome wide association study of body weight, body mass index, adiposity, and fasting glucose in 3,173 outbred rats

Obesity is a global health crisis that is influenced by both genetic and environmental factors. Rodent model organisms can be used to understand the biological and genetic basis of obesity and related morphological traits. A major advantage of model organisms is that they can be studied under uniform environmental conditions, thus reducing the complex role of environment and gene by environment interactions. Furthermore, fat pads and other tissues can be dissected and weighed, so that their role in determining body weight can be precisely defined. Highly recombinant populations allow for genetic fine-mapping of complex traits, greatly reducing the number of plausible candidate genes. We performed the largest rat GWAS ever undertaken, using 3,173 male and female adult N/NIH heterogeneous stock (HS) rats, which were created by mixing 8 inbred strains. We identified 31 independent loci for body weight, body length, body mass index, fat pad weight (retroperitoneal, epididymal, and parametrial), and fasting glucose. We observed strong evidence of pleotropic effects across multiple phenotypes. Three loci contained only a single gene (Epha5, Nrg1 and Klhl14), whereas others were larger and contained many genes. We replicated a locus containing Prlhr, and a second locus containing Adcy3, which we had previously identified in a smaller HS rat study. Finally, by subsampling our dataset, we showed an exponential growth of significant loci as sample size increased towards 3,173. Our results demonstrate the potential for rodent studies to add to our understanding of the molecular genetic factors that contribute to obesity-relevant traits and emphasize the importance of sample size.Objective Obesity is influenced by genetic and environmental factors. Despite success of human genome wide association studies (GWAS), the specific genes that confer obesity remain largely unknown. The objective of this study was to use outbred rats to identify genetic loci underlying obesity and related morphometric and metabolic traits. Methods We measured obesity-relevant traits including body weight, body length, body mass index, 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 GWAS to date. Genetic loci were identified using a linear mixed model that accounted for the complex family relationships of the HS and covariate to account for differences among the three phenotyping centers. Results We identified 32 independent loci, several of which contained only a single gene (e.g. Epha5, Nrg1 and Klhl14) or obvious candidate genes (Adcy3, Prlhr). There were strong phenotypic and genetic correlations among obesity-related traits, and extensive pleiotropy at individual loci. Conclusions These studies demonstrate utility of HS rats for investigating the genetics of obesity related traits across institutions and identify several candidate genes for future functional testing.

Dissociating addiction-related endophenotypes: Incentive salience attribution, sensation-seeking and novelty-seeking are independent traits in male and female heterogeneous stock rats

Abstract 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.Many factors contribute to addiction vulnerability, and isolating these could aid in improving treatment outcomes. Recent studies suggest that the way individuals respond to reward cues can predict the propensity for reinstatement or relapse, even after relatively brief periods of drug exposure. Other phases of addiction have been associated with different behavioral traits. For example, “sensation-seeking” is predictive of the initial propensity to take drugs; whereas “novelty-seeking” predicts compulsive drug-seeking behavior. Here we used a large sample of ∼1,600 outbred heterogeneous stock rats to determine the relationship between the propensity to attribute incentive salience to reward cues and these two other addiction-related traits. Male and female adult 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.

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.

Heterogeneous stock rats: a model to study the genetics of despair-like behavior in adolescence

Major depressive disorder (MDD) is a complex illness caused by both genetic and environmental factors. Antidepressant resistance also has a genetic component. To date, however, very few genes have been identified for major depression or antidepressant resistance. In this study, we investigated whether outbred heterogeneous stock (HS) rats would be a suitable model to uncover the genetics of depression and its connection to antidepressant resistance. The Wistar Kyoto (WKY) rat, one of the eight founders of the HS, is a recognized animal model of juvenile depression and is resistant to fluoxetine antidepressant treatment. We therefore hypothesized that adolescent HS rats would exhibit variation in both despair‐like behavior and response to fluoxetine treatment. We assessed heritability of despair‐like behavior and response to sub‐acute fluoxetine using a modified forced swim test (FST) in 4‐week‐old HS rats. We also tested whether blood transcript levels previously identified as depression biomarkers in adolescent human subjects are differentially expressed in HS rats with high vs. low FST immobility. We demonstrate heritability of despair‐like behavior in 4‐week‐old HS rats and show that many HS rats are resistant to fluoxetine treatment. In addition, blood transcript levels of Amfr, Cdr2 and Kiaa1539, genes previously identified in human adolescents with MDD, are differentially expressed between HS rats with high vs. low immobility. These data demonstrate that FST despair‐like behavior will be amenable to genetic fine‐mapping in adolescent HS rats. The overlap between human and HS blood biomarkers suggest that these studies may translate to depression in humans.