Heterogeneous Stock Populations for Analysis of Complex Traits

Establishing the hybrid rat diversity program: a resource for dissecting complex traits

Rat models have been a major model for studying complex disease mechanisms, behavioral phenotypes, environmental factors, and for drug development and discovery. Inbred rat strains control for genetic background and allow for repeated, reproducible, cellular and whole animal phenotyping. The Hybrid Rat Diversity Panel (HRDP) was designed to be a powerful panel of inbred rats with genomic, physiological, and behavioral data to serve as a resource for systems genetics. The HRDP consists of 96-98 inbred rat strains aimed to maximize power to detect specific genetic loci associated with complex traits while maximizing the genetic diversity among strains. The panel consists of 32-34 genetically diverse inbred strains and two panels of recombinant inbred panels. To establish the HRDP program, embryo resuscitation and breeding were done to establish colonies for distribution. Whole genome sequencing was performed to achieve 30X coverage. Genomic, phenotype, and strain information is available through the Hybrid Rat Diversity Panel Portal at the Rat Genome Database ( http://rgd.mcw.edu ).

A cost-effective, high-throughput, highly accurate genotyping method for outbred populations

Affordable sequencing and genotyping methods are essential for large-scale genome-wide association studies. While genotyping microarrays and reference panels for imputation are available for human subjects, nonhuman model systems often lack such options. Our lab previously demonstrated an efficient and cost-effective method to genotype heterogeneous stock rats using double-digest genotyping by sequencing. However, low-coverage whole-genome sequencing offers an alternative method that has several advantages. Here, we describe a cost-effective, high-throughput, high-accuracy genotyping method for N/NIH heterogeneous stock rats that can use a combination of sequencing data previously generated by double-digest genotyping by sequencing and more recently generated by low-coverage whole-genome sequencing data. Using double-digest genotyping-by-sequencing data from 5,745 heterogeneous stock rats (mean 0.21× coverage) and low-coverage whole-genome sequencing data from 8,760 heterogeneous stock rats (mean 0.27× coverage), we can impute 7.32 million biallelic single-nucleotide polymorphisms with a concordance rate > 99.76% compared to high-coverage (mean 33.26× coverage) whole-genome sequencing data for a subset of the same individuals. Our results demonstrate the feasibility of using sequencing data from double-digest genotyping by sequencing or low-coverage whole-genome sequencing for accurate genotyping and demonstrate techniques that may also be useful for other genetic studies in nonhuman subjects.

Y and mitochondrial chromosomes in the heterogeneous stock rat population

Genome-wide association studies typically evaluate the autosomes and sometimes the X Chromosome, but seldom consider the Y or mitochondrial (MT) Chromosomes. We genotyped the Y and MT Chromosomes in heterogeneous stock (HS) rats (Rattus norvegicus), an outbred population created from 8 inbred strains. We identified 8 distinct Y and 4 distinct MT Chromosomes among the 8 founders. However, only 2 types of each nonrecombinant chromosome were observed in our modern HS rat population (generations 81-97). Despite the relatively large sample size, there were virtually no significant associations for behavioral, physiological, metabolome, or microbiome traits after correcting for multiple comparisons. However, both Y and MT Chromosomes were strongly associated with the expression of a few genes located on those chromosomes, which provided a positive control. Our results suggest that within modern HS rats there are no Y and MT Chromosomes differences that strongly influence behavioral or physiological traits. These results do not address other ancestral Y and MT Chromosomes that do not appear in modern HS rats, nor do they address effects that may exist in other rat populations, or in other species.

The evolutionary history of wild and domestic brown rats (Rattus norvegicus)

The brown rat (Rattus norvegicus) occupies nearly every terrestrial habitat with a human presence and is one of our most important model organisms. Despite this prevalence, gaps remain in understanding the evolution of brown rat commensalism, their global dispersal, and mechanisms underlying contemporary adaptations to diverse environments. In this Review, we explore recent advances in the evolutionary history of brown rats and discuss key challenges, including finding and accurately dating historical specimens, disentangling histories of multiple domestication events, and synthesizing functional variation in wild rat populations with the development of laboratory strains. Advances in zooarchaeology and population genomics will usher in a new golden age of research on the evolutionary biology of brown rats, with positive feedbacks on their use as biomedical models.

A Cost-effective, High-throughput, Highly Accurate Genotyping Method for Outbred Populations

Affordable sequencing and genotyping methods are essential for large scale genome-wide association studies. While genotyping microarrays and reference panels for imputation are available for human subjects, non-human model systems often lack such options. Our lab previously demonstrated an efficient and cost-effective method to genotype heterogeneous stock rats using double-digest genotyping-by-sequencing. However, low-coverage whole-genome sequencing offers an alternative method that has several advantages. Here, we describe a cost-effective, high-throughput, high-accuracy genotyping method for N/NIH heterogeneous stock rats that can use a combination of sequencing data previously generated by double-digest genotyping-by-sequencing and more recently generated by low-coverage whole-genome-sequencing data. Using double-digest genotyping-by-sequencing data from 5,745 heterogeneous stock rats (mean 0.21x coverage) and low-coverage whole-genome-sequencing data from 8,760 heterogeneous stock rats (mean 0.27x coverage), we can impute 7.32 million bi-allelic single-nucleotide polymorphisms with a concordance rate >99.76% compared to high-coverage (mean 33.26x coverage) whole-genome sequencing data for a subset of the same individuals. Our results demonstrate the feasibility of using sequencing data from double-digest genotyping-by-sequencing or low-coverage whole-genome-sequencing for accurate genotyping, and demonstrate techniques that may also be useful for other genetic studies in non-human subjects.

Reversal of biological age in multiple rat organs by young porcine plasma fraction

Young blood plasma is known to confer beneficial effects on various organs in mice and rats. However, it was not known whether plasma from young adult pigs rejuvenates old rat tissues at the epigenetic level; whether it alters the epigenetic clock, which is a highly accurate molecular biomarker of aging. To address this question, we developed and validated six different epigenetic clocks for rat tissues that are based on DNA methylation values derived from n = 613 tissue samples. As indicated by their respective names, the rat pan-tissue clock can be applied to DNA methylation profiles from all rat tissues, while the rat brain, liver, and blood clocks apply to the corresponding tissue types. We also developed two epigenetic clocks that apply to both human and rat tissues by adding n = 1366 human tissue samples to the training data. We employed these six rat clocks to investigate the rejuvenation effects of a porcine plasma fraction treatment in different rat tissues. The treatment more than halved the epigenetic ages of blood, heart, and liver tissue. A less pronounced, but statistically significant, rejuvenation effect could be observed in the hypothalamus. The treatment was accompanied by progressive improvement in the function of these organs as ascertained through numerous biochemical/physiological biomarkers, behavioral responses encompassing cognitive functions. An immunoglobulin G (IgG) N-glycosylation pattern shift from pro- to anti-inflammatory also indicated reversal of glycan aging. Overall, this study demonstrates that a young porcine plasma-derived treatment markedly reverses aging in rats according to epigenetic clocks, IgG glycans, and other biomarkers of aging.

Y and Mitochondrial Chromosomes in the Heterogeneous Stock Rat Population

Genome-wide association studies typically evaluate the autosomes and sometimes the X Chromosome, but seldom consider the Y or mitochondrial Chromosomes. We genotyped the Y and mitochondrial chromosomes in heterogeneous stock rats (Rattus norvegicus), which were created in 1984 by intercrossing eight inbred strains and have subsequently been maintained as an outbred population for 100 generations. As the Y and mitochondrial Chromosomes do not recombine, we determined which founder had contributed these chromosomes for each rat, and then performed association analysis for all complex traits (n=12,055; intersection of 12,116 phenotyped and 15,042 haplotyped rats). We found the eight founders had 8 distinct Y and 4 distinct mitochondrial Chromosomes, however only two of each were observed in our modern heterogeneous stock rat population (Generations 81-97). Despite the unusually large sample size, the p-value distribution did not deviate from expectations; there were no significant associations for behavioral, physiological, metabolome, or microbiome traits after correcting for multiple comparisons. However, both Y and mitochondrial Chromosomes were strongly associated with expression of a few genes located on those chromosomes, which provided a positive control. Our results suggest that within modern heterogeneous stock rats there are no Y and mitochondrial Chromosomes differences that strongly influence behavioral or physiological traits. These results do not address other ancestral Y and mitochondrial Chromosomes that do not appear in modern heterogeneous stock rats, nor do they address effects that may exist in other rat populations, or in other species.

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.

Intermittent-access operant alcohol self-administration promotes binge-like drinking and drinking despite negative consequences in male and female heterogeneous stock rats

The study of Alcohol Use Disorders (AUD) in preclinical models is hampered by difficulty in training rodents to voluntarily consume high levels of alcohol. The intermittency of alcohol access/exposure is well known to modulate alcohol consumption (e.g., alcohol deprivation effect, intermittent-access two-bottle-choice) and recently, intermittent access operant self-administration procedures have been used to produce more intense and binge-like self-administration of intravenous psychostimulant and opioid drugs. In the present study, we sought to systematically manipulate the intermittency of operant self-administered alcohol access to determine the feasibility of promoting more intensified, binge-like alcohol consumption. To this end, 24 male and 23 female NIH Heterogeneous Stock rats were trained to self-administer 10% w/v ethanol, before being split into three different-access groups. Short Access (ShA) rats continued receiving 30-min training sessions, Long Access (LgA) rats received 16-h sessions, and Intermittent Access (IntA) rats received 16-h sessions, wherein the hourly alcohol-access periods were shortened over sessions, down to 2 min. IntA rats demonstrated an increasingly binge-like pattern of alcohol drinking in response to restriction of alcohol access, while ShA and LgA rats maintained stable intake. All groups were tested on orthogonal measures of alcohol-seeking and quinine-punished alcohol drinking. The IntA rats displayed the most punishment-resistant drinking. In a separate experiment, we replicated our main finding, that intermittent access promotes a more binge-like pattern of alcohol self-administration using 8 male and 8 female Wistar rats. In conclusion, intermittent access to self-administered alcohol promotes more intensified self-administration. This approach may be useful in developing preclinical models of binge-like alcohol consumption in AUD.

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.

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

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

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.

Genetic Background Matters: Population-Based Studies in Model Organisms for Translational Research

We are all similar but a bit different. These differences are partially due to variations in our genomes and are related to the heterogeneity of symptoms and responses to treatments that patients exhibit. Most animal studies are performed in one single strain with one manipulation. However, due to the lack of variability, therapies are not always reproducible when treatments are translated to humans. Panels of already sequenced organisms are valuable tools for mimicking human phenotypic heterogeneities and gene mapping. This review summarizes the current knowledge of mouse, fly, and yeast panels with insightful applications for translational research.

mRatBN7.2: familiar and unfamiliar features of a new rat genome reference assembly

Rat genomic tools have been slower to emerge than for those of humans and mice and have remained less thorough and comprehensive. The arrival of a new and improved rat reference genome, mRatBN7.2, in late 2020 is a welcome event. This assembly, like predecessor rat reference assemblies, is derived from an inbred Brown Norway rat. In this "user" survey we hope to provide other users of this assembly some insight into its characteristics and some assessment of its improvements as well as a few caveats that arise from the unique aspects of this assembly. mRatBN7.2 was generated by the Wellcome Sanger Institute as part of the large Vertebrate Genomes Project. This rat assembly has now joined human, mouse, chicken, and zebrafish in the National Center for Biotechnology Information (NCBI)'s Genome Reference Consortium, which provides ongoing curation of the assembly. Here we examine the technical procedures by which the assembly was created and assess how this assembly constitutes an improvement over its predecessor. We also indicate the technical limitations affecting the assembly, providing illustrations of how these limitations arise and the impact that results for this reference assembly.

The utility of a closed breeding colony of Peromyscus leucopus for dissecting complex traits

Deermice of the genus Peromyscus are well suited for addressing several questions of biologist interest, including the genetic bases of longevity, behavior, physiology, adaptation, and their ability to serve as disease vectors. Here, we explore a diversity outbred approach for dissecting complex traits in Peromyscus leucopus, a nontraditional genetic model system. We take advantage of a closed colony of deer-mice founded from 38 individuals and subsequently maintained for ∼40-60 generations. From 405 low-pass short-read sequenced deermice we accurate impute genotypes at 16 million single nucleotide polymorphisms. Conditional on observed genotypes simulations were conducted in which three different sized quantitative trait loci contribute to a complex trait under three different genetic models. Using a stringent significance threshold power was modest, largely a function of the percent variation attributable to the simulated quantitative trait loci, with the underlying genetic model having only a subtle impact. We additionally simulated 2,000 pseudo-individuals, whose genotypes were consistent with those observed in the genotyped cohort and carried out additional power simulations. In experiments employing more than 1,000 mice power is high to detect quantitative trait loci contributing greater than 2.5% to a complex trait, with a localization ability of ∼100 kb. We finally carried out a Genome-Wide Association Study on two demonstration traits, bleeding time and body weight, and uncovered one significant region. Our work suggests that complex traits can be dissected in founders-unknown P. leucopus colony mice and similar colonies in other systems using easily obtained genotypes from low-pass sequencing.

Genetic characterization of outbred Sprague Dawley rats and utility for genome-wide association studies

Sprague Dawley (SD) rats are among the most widely used outbred laboratory rat populations. Despite this, the genetic characteristics of SD rats have not been clearly described, and SD rats are rarely used for experiments aimed at exploring genotype-phenotype relationships. In order to use SD rats to perform a genome-wide association study (GWAS), we collected behavioral data from 4,625 SD rats that were predominantly obtained from two commercial vendors, Charles River Laboratories and Harlan Sprague Dawley Inc. Using double-digest genotyping-by-sequencing (ddGBS), we obtained dense, high-quality genotypes at 291,438 SNPs across 4,061 rats. This genetic data allowed us to characterize the variation present in Charles River vs. Harlan SD rats. We found that the two populations are highly diverged (F_ST > 0.4). Furthermore, even for rats obtained from the same vendor, there was strong population structure across breeding facilities and even between rooms at the same facility. We performed multiple separate GWAS by fitting a linear mixed model that accounted for population structure and using meta-analysis to jointly analyze all cohorts. Our study examined Pavlovian conditioned approach (PavCA) behavior, which assesses the propensity for rats to attribute incentive salience to reward-associated cues. We identified 46 significant associations for the various metrics used to define PavCA. The surprising degree of population structure among SD rats from different sources has important implications for their use in both genetic and non-genetic studies.

Outbred Sprague Dawley rats are among the most commonly used rats for neuroscience, physiology and pharmacological research; in the year 2020, 4,188 publications contained the keyword “Sprague Dawley”. Rats identified as “Sprague Dawley” are sold by several commercial vendors, including Charles River Laboratories and Harlan Sprague Dawley Inc. (now Envigo). Despite their widespread use, little is known about the genetic diversity of SD. We genotyped more than 4,000 SD rats, which we used for a genome-wide association study (GWAS) and to characterize genetic differences between SD rats from Charles River Laboratories and Harlan. Our analysis revealed extensive population structure both between and within vendors. The GWAS for Pavlovian conditioned approach (PavCA) identified a number of genome-wide significant loci for that complex behavioral trait. Our results demonstrate that, despite sharing an identical name, SD rats that are obtained from different vendors are very different. Future studies should carefully define the exact source of SD rats being used and may exploit their genetic diversity for genetic studies of complex traits.

Leptin Protects Against the Development and Expression of Cocaine Addiction-Like Behavior in Heterogeneous Stock Rats

Cocaine affects food intake, metabolism and bodyweight. It has been hypothesized that feeding hormones like leptin play a role in this process. Preclinical studies have shown a mutually inhibitory relationship between leptin and cocaine, with leptin also decreasing the rewarding effects of cocaine intake. But prior studies have used relatively small sample sizes and did not investigate individual differences in genetically heterogeneous populations. Here, we examined whether the role of individual differences in bodyweight and blood leptin level are associated with high or low vulnerability to addiction-like behaviors using data from 306 heterogeneous stock rats given extended access to intravenous self-administration of cocaine and 120 blood samples from 60 of these animals, that were stored in the Cocaine Biobank. Finally, we tested a separate cohort to evaluate the causal effect of exogenous leptin administration on cocaine seeking. Bodyweight was reduced due to cocaine self-administration in males during withdrawal and abstinence, but was increased in females during abstinence. However, bodyweight was not correlated with addiction-like behavior vulnerability. Blood leptin levels after ∼6 weeks of cocaine self-administration did not correlate with addiction-like behaviors, however, baseline blood leptin levels before any access to cocaine negatively predicted addiction-like behaviors 6 weeks later. Finally, leptin administration in a separate cohort of 59 animals reduced cocaine seeking in acute withdrawal and after 7 weeks of protracted abstinence. These results demonstrate that high blood leptin level before access to cocaine may be a protective factor against the development of cocaine addiction-like behavior and that exogenous leptin reduces the motivation to take and seek cocaine. On the other hand, these results also show that blood leptin level and bodyweight changes in current users are not relevant biomarkers for addiction-like behaviors.

Pierre C, Bimschleger H, Gao J, Cheng R, Polesskaya O, Solberg Woods LC, Palmer AA, Chen H. 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.

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.

On the Use of Heterogeneous Stock Mice to Map Transcriptomes Associated With Excessive Ethanol Consumption

We and many others have noted the advantages of using heterogeneous (HS) animals to map genes and gene networks associated with both behavioral and non-behavioral phenotypes. Importantly, genetically complex Mus musculus crosses provide substantially increased resolution to examine old and new relationships between gene expression and behavior. Here we report on data obtained from two HS populations: the HS/NPT derived from eight inbred laboratory mouse strains and the HS-CC derived from the eight collaborative cross inbred mouse strains that includes three wild-derived strains. Our work has focused on the genes and gene networks associated with risk for excessive ethanol consumption, individual variation in ethanol consumption and the consequences, including escalation, of long-term ethanol consumption. Background data on the development of HS mice is provided, including advantages for the detection of expression quantitative trait loci. Examples are also provided of using HS animals to probe the genes associated with ethanol preference and binge ethanol consumption.

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.

Prenatal opioid exposure and vulnerability to future substance use disorders in offspring

The heightened incidence of opioid use during pregnancy has resulted in unprecedented rates of neonates prenatally exposed to opioids. Prenatal opioid exposure (POE) results in significantly adverse medical, developmental, and behavioral outcomes in offspring. Of growing interest is whether POE contributes to future vulnerability to substance use disorders. The effects of POE on brain development is difficult to assess in humans, as the timing, dose, and route of drug exposure together with complex genetic and environmental factors affect susceptibility to addiction. Preclinical models of POE have allowed us to avoid methodological difficulties and confounding factors of POE in humans. Here, we review the effects of maternal opioid exposure on the developing brain with an emphasis on the neurobiological basis of drug addiction and on preclinical models of POE and their limitations. These studies have indicated that POE increases self-administration of drugs, reward-driven behaviors in the conditioned place paradigm, and locomotor sensitization. While addiction is multifaceted and vulnerability to drug addiction is still inconclusive in human studies of prenatally exposed infants, animal studies do provide a noteworthy corroboration of negative behavioral outcomes.

Effects of Pair Housing on Patency of Jugular Catheters in Rats (Rattus norvegicus)

Chronic vascular access devices are widely used in a variety of species for repeated blood sampling or substance administration. Jugular catheters are commonly used for studying addiction-related behaviors in rats. Rats with catheters have historically been individually housed for the duration of the study to prevent cage mates from damaging the catheter. The 2 goals of this study were to determine 1) the effects of pair housing on catheter patency and 2) the effects of pair housing on catheter patency of rats in a study of opioid self-administration and cue-induced reinstatement of opioid-seeking behavior. The latter study also represented an opportunity for experimental refinement as it evaluated the temporary use of a barrier that allowed for pair-housed rats to be physically separated. Male Heterogeneous Stock (HS; n = 24) and Sprague-Dawley (SD; n = 121) rats were allocated to either single- or pair-housed condition. To assess the effect of social housing on catheter patency, rats (HS, n = 24; SD, n = 36) were monitored in their assigned housing condition for one month, with scheduled evaluation of catheter patency and structural damage. To examine the effect of social housing on catheter patency during a study of opioid self-administration and cue-induced reinstatement of opioid-seeking behavior, rats (SD, n = 85) were monitored in their assigned housing condition with similar routine patency evaluations. Catheter patency rates between single- and pairhoused rats were not statistically different in the first experiment, and pair-housed animals were successfully maintained on an infusion study in the second experiment. The use of a barrier between pair-housed rats after surgery allowed continued social contact with no observed adverse effects. These results suggest that, pair housing is a viable option for rats with chronic vascular implants, and may improve their wellbeing by allowing them to display species-typical social behaviors.

Ratpost: a searchable database of protein-inactivating sequence variations in 40 sequenced rat-inbred strains

Rat-inbred strains are essential as scientific tools. We have analyzed the publicly available genome sequences of 40 rat-inbred strains and provide an overview of sequence variations leading to amino acid changes in protein-coding genes, premature STOP codons or loss of STOP codons, and short in-frame insertions and deletions of all protein-coding genes across all these inbred lines. We provide an overview of the predicted impact on protein function of all these affected proteins in the database, by comparing their sequence with the sequences of the rat reference strain BN/SsNHsdMcwi. We also investigate the flaws of the protein-coding sequences of this reference strain itself, by comparing them with a consensus genome. These data can be retrieved via a searchable website (Ratpost.be) and allow a global, better interpretation of genetic background effects and a source of naturally defective alleles in these 40 sequenced classical and high-priority rat-inbred strains.

Adapting Genotyping-by-Sequencing and Variant Calling for Heterogeneous Stock Rats

The heterogeneous stock (HS) is an outbred rat population derived from eight inbred rat strains. HS rats are ideally suited for genome wide association studies; however, only a few genotyping microarrays have ever been designed for rats and none of them are currently in production. To address the need for an efficient and cost effective method of genotyping HS rats, we have adapted genotype-by-sequencing (GBS) to obtain genotype information at large numbers of single nucleotide polymorphisms (SNPs). In this paper, we have outlined the laboratory and computational steps we took to optimize double digest genotype-by-sequencing (ddGBS) for use in rats. We evaluated multiple existing computational tools and explain the workflow we have used to call and impute over 3.7 million SNPs. We have also compared various rat genetic maps, which are necessary for imputation, including a recently developed map specific to the HS. Using our approach, we obtained concordance rates of 99% with data obtained using data from a genotyping array. The principles and computational pipeline that we describe could easily be adapted for use in other species for which reliable reference genome sets are available.

Nociceptin attenuates the escalation of oxycodone self-administration by normalizing CeA-GABA transmission in highly addicted rats

Approximately 25% of patients who are prescribed opioids for chronic pain misuse them, and 5 to 10% develop an opioid use disorder. Although the neurobiological target of opioids is well known, the molecular mechanisms that are responsible for the development of addiction-like behaviors in some but not all individuals are poorly known. To address this issue, we used a unique outbred rat population (heterogeneous stock) that better models the behavioral and genetic diversity that is found in humans. We characterized individual differences in addiction-like behaviors using an addiction index that incorporates the key criteria of opioid use disorder: escalated intake, highly motivated responding, and hyperalgesia. Using in vitro electrophysiological recordings in the central nucleus of the amygdala (CeA), we found that rats with high addiction-like behaviors (HA) exhibited a significant increase in γ-aminobutyric acid (GABA) transmission compared with rats with low addiction-like behaviors (LA) and naive rats. The superfusion of CeA slices with nociceptin/orphanin FQ peptide (N/OFQ; 500 nM), an endogenous opioid-like peptide, normalized GABA transmission in HA rats. Intra-CeA levels of N/OFQ were lower in HA rats than in LA rats. Intra-CeA infusions of N/OFQ (1 μg per site) reversed the escalation of oxycodone self-administration in HA rats but not in LA rats. These results demonstrate that the downregulation of N/OFQ levels in the CeA may be responsible for hyper-GABAergic tone in the CeA that is observed in individuals who develop addiction-like behaviors. Based on these results, we hypothesize that small molecules that target the N/OFQ system might be useful for the treatment of opioid use disorder.

Dopamine D3 Receptor Antagonism Reverses the Escalation of Oxycodone Self-administration and Decreases Withdrawal-Induced Hyperalgesia and Irritability-Like Behavior in Oxycodone-Dependent Heterogeneous Stock Rats

Prescription opioids, such as oxycodone, are highly effective analgesics for clinical pain management, but approximately 25% of patients who are prescribed opioids misuse them, and 5%–10% develop an opioid use disorder (OUD). Effective therapies for the prevention and treatment of opioid abuse and addiction need to be developed. The present study evaluated the effects of the highly selective dopamine D₃ receptor antagonist VK4-116 ([R]-N-[4-(4-[3-chloro-5-ethyl-2-methoxyphenyl]piperazin-1-yl)-3-hydroxybutyl]-1H-indole-2-carboxamide) on oxycodone addictive-like behaviors. We used a model of extended access to oxycodone self-administration and tested the effects of VK4-116 on the escalation of oxycodone self-administration and withdrawal-induced hyperalgesia and irritability-like behavior in male and female rats. Pretreatment with VK4-116 (5–25 mg/kg, i.p.) dose-dependently decreased the escalation of oxycodone self-administration and reduced withdrawal-induced hyperalgesia and irritability-like behavior in opioid-dependent rats. These findings demonstrate a key role for D₃ receptors in both the motivation to take opioids and negative emotional states that are associated with opioid withdrawal and suggest that D₃ receptor antagonism may be a viable therapeutic approach for the treatment of OUD.

Understanding Addiction Using Animal Models

Drug addiction is a neuropsychiatric disorder with grave personal consequences that has an extraordinary global economic impact. Despite decades of research, the options available to treat addiction are often ineffective because our rudimentary understanding of drug-induced pathology in brain circuits and synaptic physiology inhibits the rational design of successful therapies. This understanding will arise first from animal models of addiction where experimentation at the level of circuits and molecular biology is possible. We will review the most common preclinical models of addictive behavior and discuss the advantages and disadvantages of each. This includes non-contingent models in which animals are passively exposed to rewarding substances, as well as widely used contingent models such as drug self-administration and relapse. For the latter, we elaborate on the different ways of mimicking craving and relapse, which include using acute stress, drug administration or exposure to cues and contexts previously paired with drug self-administration. We further describe paradigms where drug-taking is challenged by alternative rewards, such as appetitive foods or social interaction. In an attempt to better model the individual vulnerability to drug abuse that characterizes human addiction, the field has also established preclinical paradigms in which drug-induced behaviors are ranked by various criteria of drug use in the presence of negative consequences. Separation of more vulnerable animals according to these criteria, along with other innate predispositions including goal- or sign-tracking, sensation-seeking behavior or impulsivity, has established individual genetic susceptibilities to developing drug addiction and relapse vulnerability. We further examine current models of behavioral addictions such as gambling, a disorder included in the DSM-5, and exercise, mentioned in the DSM-5 but not included yet due to insufficient peer-reviewed evidence. Finally, after reviewing the face validity of the aforementioned models, we consider the most common standardized tests used by pharmaceutical companies to assess the addictive potential of a drug during clinical trials.

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.

Rat Models of Metabolic Syndrome

Metabolic syndrome is a complex disorder that comprises several other complex disorders, including obesity, hypertension, dyslipidemia, and diabetes. There are several rat models that encompass component features of MetS. Some models are inbred strains selected for one or more traits underlying MetS; others are population models with genetic risk for MetS traits, are induced by environmental stressors such as diet, are spontaneous monogenic mutant models, or are congenic strains derived from a combination of these models. Together they can be studied to identify the genetic and physiological underpinnings of MetS to identify candidate genes or mechanisms for study in human MetS subjects.

R/qtl2: Software for Mapping Quantitative Trait Loci with High-Dimensional Data and Multiparent Populations

R/qtl2 is an interactive software environment for mapping quantitative trait loci (QTL) in experimental populations. The R/qtl2 software expands the scope of the widely-used R/qtl software package to include multiparental populations, better handles modern high-dimensional data....

R/qtl2 is an interactive software environment for mapping quantitative trait loci (QTL) in experimental populations. The R/qtl2 software expands the scope of the widely used R/qtl software package to include multiparent populations derived from more than two founder strains, such as the Collaborative Cross and Diversity Outbred mice, heterogeneous stocks, and MAGIC plant populations. R/qtl2 is designed to handle modern high-density genotyping data and high-dimensional molecular phenotypes, including gene expression and proteomics. R/qtl2 includes the ability to perform genome scans using a linear mixed model to account for population structure, and also includes features to impute SNPs based on founder strain genomes and to carry out association mapping. The R/qtl2 software provides all of the basic features needed for QTL mapping, including graphical displays and summary reports, and it can be extended through the creation of add-on packages. R/qtl2, which is free and open source software written in the R and C++ programming languages, comes with a test framework.

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.

Comparative Microbiome Signatures and Short-Chain Fatty Acids in Mouse, Rat, Non-human Primate, and Human Feces

Gut microbiome plays a fundamental role in several aspects of host health and diseases. There has been an exponential surge in the use of animal models that can mimic different phenotypes of the human intestinal ecosystem. However, data on host species-specific signatures of gut microbiome and its metabolites like short-chain fatty acids (SCFAs; i.e., acetate, propionate, and butyrate) and lactate in these models and their similarities/differences from humans remain limited, due to high variability in protocols and analyses. Here, we analyze the fecal microbiota composition and the fecal levels of SCFAs and lactate in three of the most-widely used animal models, i.e., mice, rats, and non-human primates (NHPs) and compare them with human subjects, using data generated on a single platform with same protocols. The data show several species-specific similarities and differences in the gut microbiota and fecal organic acids between these species groups. Based on β-diversity, the gut microbiota in humans seems to be closer to NHPs than to mice and rats; however, among rodents, mice microbiota appears to be closer to humans than rats. The phylum-level analyses demonstrate higher Firmicutes-Bacteroidetes ratio in humans and NHPs vs. mice and rats. Human microbiota is dominated by Bacteroides followed by Ruminococcaceae and Clostridiales. Mouse gut is predominated by members of the family S24-7 followed by those from the order Clostridiales, whereas rats and NHPs have higher abundance of Prevotella compared with mice and humans. Also, fecal levels of lactate are higher in mice and rats vs. NHPs and humans, while acetate is highest in human feces. These data of host species-specific gut microbiota signatures in some of the most widely used animal models in context to the human microbiota might reflect disparities in host factors, e.g., diets, genetic origin, gender and age, and hence call for prospective studies investigating the features of gut microbiome in such animal models by controlling for these host elements.

Genome-wide association for testis weight in the diversity outbred mouse population

Testis weight is a genetically mediated trait associated with reproductive efficiency across numerous species. We sought to evaluate the genetically diverse, highly recombinant Diversity Outbred (DO) mouse population as a tool to identify and map quantitative trait loci (QTLs) associated with testis weight. Testis weights were recorded for 502 male DO mice and the mice were genotyped on the GIGAMuga array at ~ 143,000 SNPs. We performed a genome-wide association analysis and identified one significant and two suggestive QTLs associated with testis weight. Using bioinformatic approaches, we developed a list of candidate genes and identified those with known roles in testicular size and development. Candidates of particular interest include the RNA demethylase gene Alkbh5, the cyclin-dependent kinase inhibitor gene Cdkn2c, the dynein axonemal heavy chain gene Dnah11, the phospholipase D gene Pld6, the trans-acting transcription factor gene Sp4, and the spermatogenesis-associated gene Spata6, each of which has a human ortholog. Our results demonstrate the utility of DO mice in high-resolution genetic mapping of complex traits, enabling us to identify developmentally important genes in adult mice. Understanding how genetic variation in these genes influence testis weight could aid in the understanding of mechanisms of mammalian reproductive function.