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

A revamped rat reference genome improves the discovery of genetic diversity in laboratory rats

The seventh iteration of the reference genome assembly for Rattus norvegicus-mRatBN7.2-corrects numerous misplaced segments and reduces base-level errors by approximately 9-fold and increases contiguity by 290-fold compared with its predecessor. Gene annotations are now more complete, improving the mapping precision of genomic, transcriptomic, and proteomics datasets. We jointly analyzed 163 short-read whole-genome sequencing datasets representing 120 laboratory rat strains and substrains using mRatBN7.2. We defined ∼20.0 million sequence variations, of which 18,700 are predicted to potentially impact the function of 6,677 genes. We also generated a new rat genetic map from 1,893 heterogeneous stock rats and annotated transcription start sites and alternative polyadenylation sites. The mRatBN7.2 assembly, along with the extensive analysis of genomic variations among rat strains, enhances our understanding of the rat genome, providing researchers with an expanded resource for studies involving rats.

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.

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 co-expression 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 co-expression 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 co-expressed 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, Msa4a6). 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.

Multi-Omic Approaches to Identify Genetic Factors in Metabolic Syndrome

Metabolic syndrome (MetS) is a highly heritable disease and a major public health burden worldwide. MetS diagnosis criteria are met by the simultaneous presence of any three of the following: high triglycerides, low HDL/high LDL cholesterol, insulin resistance, hypertension, and central obesity. These diseases act synergistically in people suffering from MetS and dramatically increase risk of morbidity and mortality due to stroke and cardiovascular disease, as well as certain cancers. Each of these component features is itself a complex disease, as is MetS. As a genetically complex disease, genetic risk factors for MetS are numerous, but not very powerful individually, often requiring specific environmental stressors for the disease to manifest. When taken together, all sequence variants that contribute to MetS disease risk explain only a fraction of the heritable variance, suggesting additional, novel loci have yet to be discovered. In this article, we will give a brief overview on the genetic concepts needed to interpret genome-wide association studies (GWAS) and quantitative trait locus (QTL) data, summarize the state of the field of MetS physiological genomics, and to introduce tools and resources that can be used by the physiologist to integrate genomics into their own research on MetS and any of its component features. There is a wealth of phenotypic and molecular data in animal models and humans that can be leveraged as outlined in this article. Integrating these multi-omic QTL data for complex diseases such as MetS provides a means to unravel the pathways and mechanisms leading to complex disease and promise for novel treatments. © 2022 American Physiological Society. Compr Physiol 12:1-40, 2022.

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.

Nephron-deficient HSRA rats exhibit renal injury with age but have limited renal damage from streptozotocin-induced hyperglycemia

Hypertension and diabetes are the greatest factors influencing the progression of chronic kidney disease (CKD). Investigation into the role of nephron number in CKD alone or with hypertension has revealed a strong inverse relationship between the two; however, not much is known about the connection between nephron number and diabetic kidney disease. The heterogeneous stock-derived model of unilateral renal agenesis (HSRA) rat, a novel model of nephron deficiency, provides a unique opportunity to study the association between nephron number and hypertension and diabetes on CKD. HSRA rats exhibit failure of one kidney to develop in 50-75% of offspring, whereas the remaining offspring are born with two kidneys. Rats born with one kidney (HSRA-S) develop significant renal injury with age compared with two-kidney littermates (HSRA-C). The induction of hypertension as a secondary stressor leads to significantly more renal injury in HSRA-S compared with HSRA-C rats and nephrectomized HSRA-C (HSRA-UNX) rats. The present study sought to address the hypothesis that nephron deficiency in the HSRA rat would hasten renal injury in the presence of a secondary stressor of hyperglycemia. HSRA animals did not exhibit diabetes-related traits at any age; thus, streptozotocin (STZ) was used to induce hyperglycemia in HSRA-S, HSRA-C, and HSRA-UNX rats. STZ- and vehicle-treated animals were followed for 15 wk. STZ-treated animals developed robust hyperglycemia, but in contrast to the response to hypertension, neither HSRA-S nor HSRA-UNX animals developed proteinuria compared with vehicle treatment. In total, our data indicate that hyperglycemia from STZ alone does not have a significant impact on the onset or progression of injury in young one-kidney HSRA animals.NEW & NOTEWORTHY The HSRA rat, a novel model of nephron deficiency, provides a unique opportunity to study the association between nephron number and confounding cardiovascular complications that impact kidney health. Although hypertension was previously shown to exacerbate renal injury in young HSRA animals, diabetic hyperglycemia did not lead to worse renal injury, suggesting that nephron number has limited impact on kidney injury, at least in this model.

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.

Early life stress induces hyperactivity but not increased anxiety-like behavior or ethanol drinking in outbred heterogeneous stock rats

Early life stress is known to impact vulnerability to psychopathological disorders in adulthood, including anxiety and alcohol use disorder (AUD), but the mechanisms underlying susceptibility to these outcomes are not fully understood. In the current study, we used adolescent social isolation (ASI) to determine whether Heterogeneous Stock (HS) rats, an outbred model used for genetic fine-mapping, could be used to study the genetics contributing to ASI-induced anxiety- and AUD-like behavior. We isolated (ASI) or group-housed (adolescent group-housed; AGH) 64 male HS rats at 4 weeks of age. After 5 weeks in these housing conditions, multiple anxiety and coping/despair-like behaviors were measured. All rats were then individually housed and assessed for voluntary ethanol self-administration. At euthanasia, synaptoneurosomes were isolated from a subset of brains to examine the expression of two proteins associated with alcohol drinking-related behaviors, GluA1 and SK2, in the dorsal (dHC) and ventral hippocampus (vHC). We found that ASI increased hyperactivity in the open field test relative to AGH, with no changes in other anxiety-like behaviors. Surprisingly, ASI rats demonstrated decreased immobility and increased climbing in the forced swim test relative to AGH. In contrast to prior studies by us and others, we found no difference in self-administration of 20% ethanol, with decreased ethanol self-administration in ASI relative to AGH rats at higher ethanol concentrations. Furthermore, while ASI in Long-Evans rats resulted in decreased SK2 expression in vHC synaptosomes, no differences were seen in vHC synaptosomes for SK2 or GluA1 in HS rats. These results demonstrate that HS rats are protected against many of the negative effects previously seen in response to ASI, namely anxiety-like behavior and increased ethanol self-administration. The current work suggests that a lack of change in SK2 and GluA1 expression levels in the vHC may play a role in conferring this protection.

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.

Genome-to-phenome research in rats: progress and perspectives

Because of their relatively short lifespan (<4 years), rats have become the second most used model organism to study health and diseases in humans who may live for up to 120 years. First-, second- and third-generation sequencing technologies and platforms have produced increasingly greater sequencing depth and accurate reads, leading to significant advancements in the rat genome assembly during the last 20 years. In fact, whole genome sequencing (WGS) of 47 strains have been completed. This has led to the discovery of genome variants in rats, which have been widely used to detect quantitative trait loci underlying complex phenotypes based on gene, haplotype, and sweep association analyses. DNA variants can also reveal strain, chromosome and gene functional evolutions. In parallel, phenome programs have advanced significantly in rats during the last 15 years and more than 10 databases host genome and/or phenome information. In order to discover the bridges between genome and phenome, systems genetics and integrative genomics approaches have been developed. On the other hand, multiple level information transfers from genome to phenome are executed by differential usage of alternative transcriptional start (ATS) and polyadenylation (APA) sites per gene. We used our own experiments to demonstrate how alternative transcriptome analysis can lead to enrichment of phenome-related causal pathways in rats. Development of advanced genome-to-phenome assays will certainly enhance rats as models for human biomedical research.

Identification of Candidate Risk Factor Genes for Human Idelalisib Toxicity Using a Collaborative Cross Approach

Idelalisib is a phosphatidylinositol 3-kinase inhibitor highly selective for the delta isoform that has shown good efficacy in treating chronic lymphocytic leukemia and follicular lymphoma. In clinical trials, however, idelalisib was associated with rare, but potentially serious liver and lung toxicities. In this study, we used the Collaborative Cross (CC) mouse population to identify genetic factors associated with the drug response that may inform risk management strategies for idelalisib in humans. Eight male mice (4 matched pairs) from 50 CC lines were treated once daily for 14 days by oral gavage with either vehicle or idelalisib at a dose selected to achieve clinically relevant peak plasma concentrations (150 mg/kg/day). The drug was well tolerated across all CC lines, and there were no observations of overt liver injury. Differences across CC lines were seen in drug concentration in plasma samples collected at the approximate Tmax on study Days 1, 7, and 14. There were also small but statistically significant treatment-induced alterations in plasma total bile acids and microRNA-122, and these may indicate early hepatocellular stress required for immune-mediated hepatotoxicity in humans. Idelalisib treatment further induced significant elevations in the total cell count of terminal bronchoalveolar lavage fluid, which may be analogous to pneumonitis observed in the clinic. Genetic mapping identified loci associated with interim plasma idelalisib concentration and the other 3 treatment-related endpoints. Thirteen priority candidate quantitative trait genes identified in CC mice may now guide interrogation of risk factors for adverse drug responses associated with idelalisib in humans.

High-fat diet negatively impacts both metabolic and behavioral health in outbred heterogeneous stock rats

Obesity is influenced by genetics and diet and has wide ranging comorbidities, including anxiety and depressive disorders. Outbred heterogeneous stock (HS) rats are used for fine-genetic mapping of complex traits and may be useful for understanding gene by diet interactions. In this study, HS rats were fed diets containing 60% kcal from fat (high-fat diet, HFD) or 10% kcal from fat (low-fat diet, LFD) and tested for metabolic (study 1) and behavioral (study 2) outcomes. In study 1, we measured glucose tolerance, fasting glucose and insulin, fat pad weights and despair-like behavior in the forced swim test (FST). In study 2, we assessed anxiety-like (elevated plus maze, EPM; open field test, OFT) and despair-like/coping (splash test, SpT; and FST) behaviors. Body weight and food intake were measured weekly in both studies. We found negative effects of HFD on metabolic outcomes, including increased body weight and fat pad weights, decreased glucose tolerance, and increased fasting insulin. We also found negative effects of HFD on despair-like/coping and anxiety-like behaviors. These include increased immobility in the FST, decreased open arm time in the EPM, and increased movement and rest episodes and decreased rearing in the OFT. The diet-induced changes in EPM and OFT were independent of overall locomotion. Additionally, diet-induced changes in OFT behaviors were independent of adiposity, while adiposity was a confounding factor for EPM and FST behavior. This work establishes the HS as a model to study gene by diet interactions affecting metabolic and behavioral health.

Using Genetic and Species Diversity to Tackle Kidney Disease

Progress in the identification of causal genes and understanding of the mechanism underlying kidney disease is hindered by the almost exclusive use of a few animal models with restrictive monogenic backgrounds that may be more resistant to kidney disease compared with humans and, therefore, poor models. Exploring the large genetic diversity in classical animal models, such as mice and rats, and leveraging species diversity will allow us to use the genetic advantages of zebrafish, Drosophila, and other species, to develop both new animal models that are more relevant to the study of human kidney disease and potential therapies.

Whole genome sequencing and novel candidate genes for CAKUT and altered nephrogenesis in the HSRA rat

The HSRA rat is a model of congenital abnormalities of the kidney and urogenital tract (CAKUT). Our laboratory has used this model to investigate the role of nephron number (functional unit of the kidney) in susceptibility to develop kidney disease as 50-75% offspring are born with a single kidney (HSRA-S), while 25-50% are born with two kidneys (HSRA-C). HSRA-S rats develop increased kidney injury and hypertension with age compared with nephrectomized two-kidney animals (HSRA-UNX), suggesting that even slight differences in nephron number can be an important driver in decline in kidney function. The HSRA rat was selected and inbred from a family of outbred heterogeneous stock (NIH-HS) rats that exhibited a high incidence of CAKUT. The HS model was originally developed from eight inbred strains (ACI, BN, BUF, F344, M520, MR, WKY, and WN). The genetic make-up of the HSRA is therefore a mosaic of these eight inbred strains. Interestingly, the ACI progenitor of the HS model exhibits CAKUT in 10-15% of offspring with the genetic cause being attributed to the presence of a long-term repeat (LTR) within exon 1 of the c-Kit gene. Our hypothesis is that the HSRA and ACI share this common genetic cause, but other alleles in the HSRA genome contribute to the increased penetrance of CAKUT (75% HSRA vs. 15% in ACI). To facilitate genetic studies and better characterize the model, we sequenced the whole genome of the HSRA to a depth of ~50×. A genome-wide variant analysis of high-impact variants identified a number of novel genes that could be linked to CAKUT in the HSRA model. In summary, the identification of new genes/modifiers that lead to CAKUT/loss of one kidney in the HSRA model will provide greater insight into association between kidney development and susceptibility to develop cardiovascular disease later in life.

Using Heterogeneous Stocks for Fine-Mapping Genetically Complex Traits

In this chapter we will review both the rationale and experimental design for using Heterogeneous Stock (HS) populations for fine-mapping of complex traits in mice and rats. We define an HS as an outbred population derived from an intercross between two or more inbred strains. HS have been used to perform genome-wide association studies (GWAS) for multiple behavioral, physiological, and gene expression traits. GWAS using HS require four key steps, which we review: selection of an appropriate HS population, phenotyping, genotyping, and statistical analysis. We provide advice on the selection of an HS, comment on key issues related to phenotyping, discuss genotyping methods relevant to these populations, and describe statistical genetic analyses that are applicable to genetic analyses that use HS.

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.

Social and anxiety-like behaviors contribute to nicotine self-administration in adolescent outbred rats

Both emotional and social traits interact with genetic factors to influence smoking behavior. We previously established a socially acquired nicotine intravenous self-administration model where social learning of a nicotine-associated odor cue reversed conditioned flavor aversion and promoted nicotine intake. In this study, we first phenotyped ~800 adolescent heterogeneous stock rats in open field, novel object interaction, social interaction, elevated plus maze, and marble burying behaviors. These rats were then phenotyped on socially acquired nicotine self-administration. We found 243 significant correlations between different behavioral tests. Principal component regression analysis found that ~10-20% of the variance in nicotine-related measures, such as intake during the first or the last three fixed-ratio sessions, the progressive ratio session, and reinstatement behavior, can be explained by variations in behavioral traits. Factors corresponding to social behavior and anxiety were among the strongest predictors of nicotine intake and reinstatement of nicotine-seeking behavior. We also found many sex differences in behavioral measures. These data indicated that the genetic diversity of this population, in combination with social behaviour and anxiety, are significant contributors to the divergent nicotine self-administration behavior and indicated a high probability of discovering sex-specific genetic mechanisms for nicotine intake in future genome-wide association studies.

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.

Two-pore channels and disease

Two-pore channels (TPCs) are Ca²⁺-permeable endo-lysosomal ion channels subject to multi-modal regulation. They mediate their physiological effects through releasing Ca²⁺ from acidic organelles in response to cues such as the second messenger, NAADP. Here, we review emerging evidence linking TPCs to disease. We discuss how perturbing both local and global Ca²⁺ changes mediated by TPCs through chemical and/or molecular manipulations can induce or reverse disease phenotypes. We cover evidence from models of Parkinson's disease, non-alcoholic fatty liver disease, Ebola infection, cancer, cardiac dysfunction and diabetes. A need for more drugs targeting TPCs is identified.

Joint Analysis of Strain and Parent-of-Origin Effects for Recombinant Inbred Intercrosses Generated from Multiparent Populations with the Collaborative Cross as an Example

Multiparent populations (MPP) have become popular resources for complex trait mapping because of their wider allelic diversity and larger population size compared with traditional two-way recombinant inbred (RI) strains. In mice, the collaborative cross (CC) is one of the most popular MPP and is derived from eight genetically diverse inbred founder strains. The strategy of generating RI intercrosses (RIX) from MPP in general and from the CC in particular can produce a large number of completely reproducible heterozygote genomes that better represent the (outbred) human population. Since both maternal and paternal haplotypes of each RIX are readily available, RIX is a powerful resource for studying both standing genetic and epigenetic variations of complex traits, in particular, the parent-of-origin (PoO) effects, which are important contributors to many complex traits. Furthermore, most complex traits are affected by >1 genes, where multiple quantitative trait locus mapping could be more advantageous. In this paper, for MPP-RIX data but taking CC-RIX as a working example, we propose a general Bayesian variable selection procedure to simultaneously search for multiple genes with founder allelic effects and PoO effects. The proposed model respects the complex relationship among RIX samples, and the performance of the proposed method is examined by extensive simulations.

Editor's Highlight: Candidate Risk Factors and Mechanisms for Tolvaptan-Induced Liver Injury Are Identified Using a Collaborative Cross Approach

Clinical trials of tolvaptan showed it to be a promising candidate for the treatment of Autosomal Dominant Polycystic Kidney Disease (ADPKD) but also revealed potential for idiosyncratic drug-induced liver injury (DILI) in this patient population. To identify risk factors and mechanisms underlying tolvaptan DILI, 8 mice in each of 45 strains of the genetically diverse Collaborative Cross (CC) mouse population were treated with a single oral dose of either tolvaptan or vehicle. Significant elevations in plasma alanine aminotransferase (ALT) were observed in tolvaptan-treated animals in 3 of the 45 strains. Genetic mapping coupled with transcriptomic analysis in the liver was used to identify several candidate susceptibility genes including epoxide hydrolase 2, interferon regulatory factor 3, and mitochondrial fission factor. Gene pathway analysis revealed that oxidative stress and immune response pathways were activated in response to tolvaptan treatment across all strains, but genes involved in regulation of bile acid homeostasis were most associated with tolvaptan-induced elevations in ALT. Secretory leukocyte peptidase inhibitor (Slpi) mRNA was also induced in the susceptible strains and was associated with increased plasma levels of Slpi protein, suggesting a potential serum marker for DILI susceptibility. In summary, tolvaptan induced signs of oxidative stress, mitochondrial dysfunction, and innate immune response in all strains, but variation in bile acid homeostasis was most associated with susceptibility to the liver response. This CC study has indicated potential mechanisms underlying tolvaptan DILI and biomarkers of susceptibility that may be useful in managing the risk of DILI in ADPKD patients.

Heterogeneous Stock Populations for Analysis of Complex Traits

Heterogeneous Stock (HS) populations allow for fine-resolution genetic mapping of a variety of complex traits. HS mice and rats were created from breeding together eight inbred strains, followed by maintaining the colony in a manner that minimizes inbreeding. After 50 or more generations of breeding, the resulting animals' chromosomes represent a genetic mosaic of the founders' haplotypes, with the average distance between recombination events in the centiMorgan range. This allows for genetic mapping to only a few Mb, a much smaller region than what can be identified using traditional F2 intercross or backcross mapping strategies. HS animals have been used to fine-map a variety of complex traits including anxiety and fear behaviors, diabetes, asthma, and heart disease, among others. Once a quantitative trait locus (QTL) has been identified, founder sequence and expression analysis can be used to identify underlying causal genes. In the following review, we provide an overview of how HS rats and mice have been used to identify genetic loci, and in some cases the causal genes, underlying complex traits. We discuss the creation and breeding strategies for both HS rats and mice. We then discuss the statistical analyses used to identify genetic loci, as well as strategies to identify causal genes underlying these loci. We end the chapter by discussing limitations faced when using HS populations, including several statistical challenges that have not been fully resolved.

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.

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.

Premature responding is associated with approach to a food cue in male and female heterogeneous stock rats

RATIONALE

Disorders of behavioral regulation, including attention deficit hyperactivity disorder (ADHD) and drug addiction, are in part due to poor inhibitory control, attentional deficits, and hyper-responsivity to reward-associated cues.

OBJECTIVES

To determine whether these traits are related, we tested genetically variable male and female heterogeneous stock rats in the choice reaction time (CRT) task and Pavlovian conditioned approach (PavCA). Sex differences in the response to methylphenidate during the CRT were also assessed.

METHODS

In the CRT task, rats were required to withhold responding until one of two lights indicated whether responses into a left or right port would be reinforced with water. Reaction time on correct trials and premature responses were the operational definitions of attention and response inhibition, respectively. Rats were also pretreated with oral methylphenidate (0, 2, 4 mg/kg) during the CRT task to determine whether this drug would improve performance. Subsequently, during PavCA, presentation of an illuminated lever predicted the delivery of a food pellet into a food-cup. Lever-directed approach (sign-tracking) and food-cup approach (goal-tracking) were the primary measures, and rats were categorized as "sign-trackers" and "goal-trackers" using an index based on these measures.

RESULTS

Sign-trackers made more premature responses than goal-trackers but showed no differences in reaction time. There were sex differences in both tasks, with females having higher sign-tracking, completing more CRT trials, and making more premature responses after methylphenidate administration.

CONCLUSIONS

These results indicate that response inhibition is related to reward-cue responsivity, suggesting that these traits are influenced by common genetic factors.

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.

Bayesian modeling of haplotype effects in multiparent populations

A general Bayesian model, Diploffect, is described for estimating the effects of founder haplotypes at quantitative trait loci (QTL) detected in multiparental genetic populations; such populations include the Collaborative Cross (CC), Heterogeneous Socks (HS), and many others for which local genetic variation is well described by an underlying, usually probabilistically inferred, haplotype mosaic. Our aim is to provide a framework for coherent estimation of haplotype and diplotype (haplotype pair) effects that takes into account the following: uncertainty in haplotype composition for each individual; uncertainty arising from small sample sizes and infrequently observed haplotype combinations; possible effects of dominance (for noninbred subjects); genetic background; and that provides a means to incorporate data that may be incomplete or has a hierarchical structure. Using the results of a probabilistic haplotype reconstruction as prior information, we obtain posterior distributions at the QTL for both haplotype effects and haplotype composition. Two alternative computational approaches are supplied: a Markov chain Monte Carlo sampler and a procedure based on importance sampling of integrated nested Laplace approximations. Using simulations of QTL in the incipient CC (pre-CC) and Northport HS populations, we compare the accuracy of Diploffect, approximations to it, and more commonly used approaches based on Haley–Knott regression, describing trade-offs between these methods. We also estimate effects for three QTL previously identified in those populations, obtaining posterior intervals that describe how the phenotype might be affected by diplotype substitutions at the modeled locus.

Genomes and phenomes of a population of outbred rats and its progenitors

Finding genetic variants that contribute to phenotypic variation is one of the main challenges of modern genetics. We used an outbred population of rats (Heterogeneous Stock, HS) in a combined sequence-based and genetic mapping analysis to identify sequence variants and genes contributing to complex traits of biomedical relevance. Here we describe the sequences of the eight inbred progenitors of the HS and the variants that segregate between them. We report the genotyping of 1,407 HS rats, and the collection from 2,006 rats of 195 phenotypic measures that are relevant to models of anxiety, type 2 diabetes, hypertension and osteoporosis. We make available haplotype dosages for the 1,407 genotyped rats, since genetic mapping in the HS is best carried out by reconstructing each HS chromosome as a mosaic of the progenitor genomes. Finally, we have deposited an R object that makes it easy to incorporate our sequence data into any genetic study of HS rats. Our genetic data are available for both Rnor3.4 and Rnor5.0 rat assemblies.

Propensity for social interaction predicts nicotine-reinforced behaviors in outbred rats

Social and genetic factors can influence smoking behavior. Using olfactogustatory stimuli as the sensory cue for intravenous nicotine self-administration (SA), we previously showed that social learning of nicotine contingent odor cue prevented rats from developing conditioned taste aversion and allowed them to instead establish stable nicotine SA. We hypothesized that genetic factors influenced socially acquired nicotine SA. A heterogeneous stock (HS; N/NIH) of outbred rats was trained to self-administer nicotine using the social learning protocol. Both male and female HS rats acquired nicotine SA, but females self-administered more nicotine than males. After extinction, the context previously paired with nicotine SA, in conjunction with socially transmitted drug cues, was sufficient to cause reinstatement of drug-seeking behavior. Wide variation in both nicotine intake and reinstatement was observed. Using multiple regression analysis, we found that measures of social interaction were significant predictors of nicotine intake and reinstatement of drug seeking in both males and females. Furthermore, measures of depression were predictors of nicotine intake in both males and females, anxiety was a predictor only in males and response to novelty was a predictor only in females. In males, measures of both depression and anxiety predicted nicotine reinstatement. Together, these data supported the ideas that genetically determined propensities for emotional and social phenotypes are significant determinants for nicotine-reinforced behavior, and that the HS rat is a suitable tool for dissecting genetic mechanisms that may underlie the interaction between social behavior, anxiety, depression and smoking.

QTL mapping in outbred populations: successes and challenges

Quantitative trait locus (QTL) mapping in animal populations has been a successful strategy for identifying genomic regions that play a role in complex diseases and traits. When conducted in an F2 intercross or backcross population, the resulting QTL is frequently large, often encompassing 30 Mb or more and containing hundreds of genes. To narrow the locus and identify candidate genes, additional strategies are needed. Congenic strains have proven useful but work less well when there are multiple tightly linked loci, frequently resulting in loss of phenotype. As an alternative, we discuss the use of highly recombinant outbred models for directly fine-mapping QTL to only a few megabases. We discuss the use of several currently available models such as the advanced intercross (AI), heterogeneous stocks (HS), the diversity outbred (DO), and commercially available outbred stocks (CO). Once a QTL has been fine-mapped, founder sequence and expression QTL mapping can be used to identify candidate genes. In this regard, the large number of alleles found in outbred stocks can be leveraged to identify causative genes and variants. We end this review by discussing some important statistical considerations when analyzing outbred populations. Fine-resolution mapping in outbred models, coupled with full genome sequence, has already led to the identification of several underlying causative genes for many complex traits and diseases. These resources will likely lead to additional successes in the coming years.

Using the emerging Collaborative Cross to probe the immune system

The Collaborative Cross (CC) is an emerging panel of recombinant inbred (RI) mouse strains. Each strain is genetically distinct but all descended from the same eight inbred founders. In 66 strains from incipient lines of the CC (pre-CC), as well as the 8 CC founders and some of their F1 offspring, we examined subsets of lymphocytes and antigen-presenting cells. We found significant variation among the founders, with even greater diversity in the pre-CC. Genome-wide association using inferred haplotypes detected highly significant loci controlling B-to-T cell ratio, CD8 T-cell numbers, CD11c and CD23 expression. Comparison of overall strain effects in the CC founders with strain effects at QTL in the pre-CC revealed sharp contrasts in the genetic architecture of two traits with significant loci: variation in CD23 can be explained largely by additive genetics at one locus, whereas variation in B-to-T ratio has a more complex etiology. For CD23, we found a strong QTL whose confidence interval contained the CD23 structural gene Fcer2a. Our data on the pre-CC demonstrate the utility of the CC for studying immunophenotypes and the value of integrating founder, CC and F1 data. The extreme immunophenotypes observed could have pleiotropic effects in other CC experiments.

Strong genetic influences on measures of behavioral-regulation among inbred rat strains

A fundamental challenge for any complex nervous system is to regulate behavior in response to environmental challenges. Three measures of behavioral-regulation were tested in a panel of eight inbred rat strains. These measures were: (1) sensation seeking as assessed by locomotor response to novelty and the sensory reinforcing effects of light onset, (2) attention and impulsivity, as measured by a choice reaction time task and (3) impulsivity as measured by a delay discounting task. Deficient behavioral-regulation has been linked to a number of psychopathologies, including ADHD, Schizophrenia, Autism, drug abuse and eating disorders. Eight inbred rat strains (August Copenhagen Irish, Brown Norway, Buffalo, Fischer 344, Wistar Kyoto, Spontaneous Hypertensive Rat, Lewis, Dahl Salt Sensitive) were tested. With n = 9 for each strain, we observed robust strain differences for all tasks; heritability was estimated between 0.43 and 0.66. Performance of the eight inbred rat strains on the choice reaction time task was compared to the performance of outbred Sprague Dawley (n = 28) and Heterogeneous strain rats (n = 48). The results indicate a strong genetic influence on complex tasks related to behavioral-regulation and indicate that some of the measures tap common genetically driven processes. Furthermore, our results establish the potential for future studies aimed at identifying specific alleles that influence variability for these traits. Identification of such alleles could contribute to our understanding of the molecular genetic basis of behavioral-regulation, which is of fundamental importance and likely contributes to multiple psychiatric disorders.

Rats are the smart choice: Rationale for a renewed focus on rats in behavioral genetics

Due in part to their rich behavioral repertoire rats have been widely used in behavioral studies of drug abuse-related traits for decades. However, the mouse became the model of choice for researchers exploring the genetic underpinnings of addiction after the first mouse study was published demonstrating the capability of engineering the mouse genome through embryonic stem cell technology. The sequencing of the mouse genome and more recent re-sequencing of numerous inbred mouse strains have further cemented the status of mice as the premier mammalian organism for genetic studies. As a result, many of the behavioral paradigms initially developed and optimized for rats have been adapted to mice. However, numerous complex and interesting drug abuse-related behaviors that can be studied in rats are very difficult or impossible to adapt for use in mice, impeding the genetic dissection of those traits. Now, technological advances have removed many of the historical limitations of genetic studies in rats. For instance, the rat genome has been sequenced and many inbred rat strains are now being re-sequenced and outbred rat stocks are being used to fine-map QTLs. In addition, it is now possible to create "knockout" rats using zinc finger nucleases (ZFN), transcription activator-like effector nucleases (TALENs) and related techniques. Thus, rats can now be used to perform quantitative genetic studies of sophisticated behaviors that have been difficult or impossible to study in mice. This article is part of a Special Issue entitled 'NIDA 40th Anniversary Issue'.