Using Heterogeneous Stocks for Fine-Mapping Genetically Complex Traits

Pierre CL, Missfeldt Sanches T, Bimschleger HV, Gao J, Cheng R, Nguyen K, Holl KL, Polesskaya O, Ishiwari K, Chen H, Robinson TE, Flagel SB, Solberg Woods LC, Palmer AA, Meyer PJ. Genetic Loci Influencing Cue-Reactivity in Heterogeneous Stock Rats

Addiction vulnerability is associated with the tendency to attribute incentive salience to reward predictive cues. Both addiction and the attribution of incentive salience are influenced by environmental and genetic factors. To characterize the genetic contributions to incentive salience attribution, we performed a genome-wide association study (GWAS) in a cohort of 1596 heterogeneous stock (HS) rats. Rats underwent a Pavlovian conditioned approach task that characterized the responses to food-associated stimuli ("cues"). Responses ranged from cue-directed "sign-tracking" behavior to food-cup directed "goal-tracking" behavior (12 measures, SNP heritability: 0.051-0.215). Next, rats performed novel operant responses for unrewarded presentations of the cue using the conditioned reinforcement procedure. GWAS identified 14 quantitative trait loci (QTLs) for 11 of the 12 traits across both tasks. Interval sizes of these QTLs varied widely. Seven traits shared a QTL on chromosome 1 that contained a few genes (e.g., Tenm4, Mir708) that have been associated with substance use disorders and other psychiatric disorders in humans. Other candidate genes (e.g., Wnt11, Pak1) in this region had coding variants and expression-QTLs in mesocorticolimbic regions of the brain. We also conducted a Phenome-Wide Association Study (PheWAS) on addiction-related behaviors in HS rats and found that the QTL on chromosome 1 was also associated with nicotine self-administration in a separate cohort of HS rats. These results provide a starting point for the molecular genetic dissection of incentive motivational processes and provide further support for a relationship between the attribution of incentive salience and drug abuse-related traits.

Novel insights into the genetic architecture and mechanisms of host/microbiome interactions from a multi-cohort analysis of outbred laboratory rats

The intestinal microbiome influences health and disease. Its composition is affected by host genetics and environmental exposures. Understanding host genetic effects is critical but challenging in humans, due to the difficulty of detecting, mapping and interpreting them. To address this, we analysed host genetic effects in four cohorts of outbred laboratory rats exposed to distinct but controlled environments. We found that polygenic host genetic effects were consistent across environments. We identified three replicated microbiome-associated loci. One involved a sialyltransferase gene and Paraprevotella and we found a similar association, between ST6GAL1 and Paraprevotella, in a human cohort. Given Paraprevotella's known immunity-potentiating functions, this suggests ST6GAL1's effects on IgA nephropathy and COVID-19 breakthrough infections may be mediated by Paraprevotella. Moreover, we found evidence of indirect genetic effects on microbiome phenotypes, which substantially increased their total genetic variance. Finally, we identified a novel mechanism whereby indirect genetic effects can contribute to "missing heritability".

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

Pierre CL, Sumner S, Tripi J, Wang T, Chen H, Flagel S, Ishiwari K, Meyer P, Polesskaya O, Saba L, Solberg Woods LC, Palmer AA, Im HK. RatXcan: A framework for cross-species integration of genome-wide association and gene expression data

Genome-wide association studies (GWAS) have implicated specific alleles and genes as risk factors for numerous complex traits. However, translating GWAS results into biologically and therapeutically meaningful discoveries remains extremely challenging. Most GWAS results identify noncoding regions of the genome, suggesting that differences in gene regulation are the major driver of trait variability. To better integrate GWAS results with gene regulatory polymorphisms, we previously developed PrediXcan (also known as "transcriptome-wide association studies" or TWAS), which maps SNPs to predicted gene expression using GWAS data. In this study, we developed RatXcan, a framework that extends this methodology to outbred heterogeneous stock (HS) rats. RatXcan accounts for the close familial relationships among HS rats by modeling the relatedness with a random effect that encodes the genetic relatedness. RatXcan also corrects for polygenic-driven inflation because of the equivalence between a relatedness random effect and the infinitesimal polygenic model. To develop RatXcan, we trained transcript predictors for 8,934 genes using reference genotype and expression data from five rat brain regions. We found that the cis genetic architecture of gene expression in both rats and humans was sparse and similar across brain tissues. We tested the association between predicted expression in rats and two example traits (body length and BMI) using phenotype and genotype data from 5,401 densely genotyped HS rats and identified a significant enrichment between the genes associated with rat and human body length and BMI. Thus, RatXcan represents a valuable tool for identifying the relationship between gene expression and phenotypes across species and paves the way to explore shared biological mechanisms of complex traits.

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.

Construction and evaluation of a new rat reference genome assembly, GRCr8, from long reads and long-range scaffolding

We report the construction and analysis of a new reference genome assembly for Rattus norvegicus, the laboratory rat, a widely used experimental animal model organism. The assembly has been adopted as the rat reference assembly by the Genome Reference Consortium and is named GRCr8. The assembly has employed 40× Pacific Biosciences (PacBio) HiFi sequencing coverage and scaffolding using optical mapping and Hi-C. We used genomic DNA from a male BN/NHsdMcwi (BN) rat of the same strain and from the same colony as the prior reference assembly, mRatBN7.2. The assembly is at chromosome level with 98.7% of the sequence assigned to chromosomes. All chromosomes have increased in size compared with the prior assembly and k-mer analysis indicates that the subject animal is fully inbred and that the genome is represented as a single haploid assembly. Notable increases are observed in Chromosomes 3, 11, and 12 in the prospective rDNA regions. In addition, Chr Y has increased threefold in size and is more consistent with the rat karyotype than previous assemblies. Several other chromosomes have grown by the incorporation of sizable discrete new blocks. These contain highly repetitive sequences and encode numerous previously unannotated genes. In addition, centromeric sequences are incorporated in most chromosomes. Genome annotation has been performed by NCBI RefSeq, which confirms improvement in assembly quality and adds more than 1100 new protein coding genes. PacBio Iso-Seq data have been acquired from multiple tissues of the subject animal and are released concurrently with the new assembly to aid further analyses.

Large-scale characterization of cocaine addiction-like behaviors reveals that escalation of intake, aversion-resistant responding, and breaking-points are highly correlated measures of the same construct

Addiction is commonly characterized by escalation of drug intake, compulsive drug seeking, and continued use despite harmful consequences. However, the factors contributing to the transition from moderate drug use to these problematic patterns remain unclear, particularly regarding the role of sex. Many preclinical studies have been limited by small sample sizes, low genetic diversity, and restricted drug access, making it challenging to model significant levels of intoxication or dependence and translate findings to humans. To address these limitations, we characterized addiction-like behaviors in a large sample of >500 outbred heterogeneous stock (HS) rats using an extended cocaine self-administration paradigm (6 hr/daily). We analyzed individual differences in escalation of intake, progressive ratio (PR) responding, continued use despite adverse consequences (contingent foot shocks), and irritability-like behavior during withdrawal. Principal component analysis showed that escalation of intake, progressive ratio responding, and continued use despite adverse consequences loaded onto a single factor that was distinct from irritability-like behaviors. Categorizing rats into resilient, mild, moderate, and severe addiction-like phenotypes showed that females exhibited higher addiction-like behaviors, with a lower proportion of resilient individuals compared to males. These findings suggest that, in genetically diverse rats with extended drug access, escalation of intake, continued use despite adverse consequences, and PR responding are highly correlated measures of a shared underlying construct. Furthermore, our results highlight sex differences in resilience to addiction-like behaviors.

How life events may confer vulnerability to addiction: the role of epigenetics

Substance use disorder (SUD) represents a large and growing global health problem. Despite the strong addictive potency of drugs of abuse, only a minority of those exposed develop SUDs. While certain life experiences (e.g., childhood trauma) may increase subsequent vulnerability to SUDs, mechanisms underlying these effects are not yet well understood. Given the chronic and relapsing nature of SUDs, and the length of time that can elapse between prior life events and subsequent drug exposure, changes in SUD vulnerability almost certainly involve long-term epigenetic dysregulation. To validate this idea, functional effects of specific epigenetic modifications in brain regions mediating reinforcement learning (e.g., nucleus accumbens, prefrontal cortex) have been investigated in a variety of animal models of SUDs. In addition, the effects of epigenetic modifications produced by prior life experiences on subsequent SUD vulnerability have been studied, but mostly in a correlational manner. Here, we review how epigenetic mechanisms impact SUD-related behavior in animal models and summarize our understanding of the relationships among life experiences, epigenetic regulation, and future vulnerability to SUDs. Despite variations in study design, epigenetic modifications that most consistently affect SUD-related behavior are those that produce predominantly unidirectional effects on gene regulation, such as DNA methylation and histone phosphorylation. Evidence explicitly linking environmentally induced epigenetic modifications to subsequent SUD-related behavior is surprisingly sparse. We conclude by offering several directions for future research to begin to address this critical research gap.

Does tolerance to ethanol-induced ataxia explain the sensitized response to ethanol?

Under conditions of repeated exposure to ethanol, a sensitized locomotor stimulant response develops in some strains of mice. It has been hypothesized that the sensitized response is a consequence of tolerance development to the sedative/incoordinating effects of ethanol. Conversely, ethanol-induced sensitization and tolerance may be independent effects of repeated ethanol exposure. A published study in C57BL/6J by DBA/2J recombinant inbred strains concluded that the two phenomena are not genetically related and thus perhaps mechanistically distinct. To extend evaluation beyond the genetic variance found in C57BL/6J and DBA/2J mice and examine phenotypic associations, we simultaneously measured ethanol-induced sensitization and tolerance in a genetically diverse panel of 15 standard inbred mouse strains and a genetically heterogeneous stock that was produced by the intercrossing of eight inbred mouse strains. Changes in activity counts and ataxia ratio across repeated ethanol treatments indexed sensitization and tolerance, respectively. Photocell beam breaks provided the measure of activity, and foot slip errors corrected for activity in a grid test provided a measure of coordination. The results were strain and individual dependent. The genetic correlation between magnitude of sensitization and tolerance was not significant in the panel of inbred strains, but when individual data were correlated, without regard to strain, there was a significant correlation. This relationship was also significant in the genetically heterogeneous population of mice. However, magnitude of tolerance explained only 10% of the variance in sensitization among individuals of the inbred strain population, whereas it explained 44% of the variance among individuals of the eight-strain cross. When repeated exposures to ethanol were disassociated from the test apparatus, this relationship in the eight-strain cross disappeared. Furthermore, days to peak sensitization and tolerance across days did not perfectly mirror each other. Overall, our data do not support shared genetic mechanisms in sensitization and tolerance development but suggest a partial relationship among individuals that could be related to drug-environment associations.

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.

Long-access heroin self-administration induces region specific reduction of grey matter volume and microglia reactivity in the rat

In opioid use disorder (OUD) patients, a decrease in brain grey matter volume (GMV) has been reported. It is unclear whether this is the consequence of prolonged exposure to opioids or is a predisposing causal factor in OUD development. To investigate this, we conducted a structural MRI longitudinal study in NIH Heterogeneous Stock rats exposed to heroin self-administration and age-matched naïve controls housed in the same controlled environment. Structural MRI scans were acquired before (MRI1) and after (MRI2) a prolonged period of long access heroin self-administration resulting in escalation of drug intake. Heroin intake resulted in reduced GMV in various cortical and sub-cortical brain regions. In drug-naïve controls no difference was found between MRI1 and MRI2. Notably, the degree of GMV reduction in the medial prefrontal cortex (mPFC) and the insula positively correlated with the amount of heroin consumed and the escalation of heroin use. In a preliminary gene expression analysis, we identified a number of transcripts linked to immune response and neuroinflammation. This prompted us to hypothesize a link between changes in microglia homeostasis and loss of GMV. For this reason, we analyzed the number and morphology of microglial cells in the mPFC and insula. The number of neurons and their morphology was also evaluated. The primary motor cortex, where no GMV change was observed, was used as negative control. We found no differences in the number of neurons and microglia cells following heroin. However, in the same regions where reduced GMV was detected, we observed a shift towards a rounder shape and size reduction in microglia, suggestive of their homeostatic change towards a reactive state. Altogether these findings suggest that escalation of heroin intake correlates with loss of GMV in specific brain regions and that this phenomenon is linked to changes in microglial morphology.

Genomic Loci Influencing Cue-Reactivity in Heterogeneous Stock Rats

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

Long-access heroin self-administration induces region specific reduction of grey matter volume and microglia reactivity in the rat

In opioid use disorder (OUD) patients, a decrease in brain grey matter volume (GMV) has been reported. It is unclear whether this is the consequence of prolonged exposure to opioids or is a predisposing causal factor in OUD development. To investigate this, we conducted a structural MRI longitudinal study in NIH Heterogeneous Stock rats exposed to heroin self-administration and age-matched naïve controls housed in the same controlled environment. Structural MRI scans were acquired before (MRI 1 ) and after (MRI 2 ) a prolonged period of long access heroin self-administration resulting in escalation of drug intake. Heroin intake resulted in reduced GMV in various cortical and sub-cortical brain regions. In drug-naïve controls no difference was found between MRI 1 and MRI 2 . Notably, the degree of GMV reduction in the medial prefrontal cortex (mPFC) and the insula positively correlated with the amount of heroin consumed and the escalation of heroin use. In a preliminary gene expression analysis, we identified a number of transcripts linked to immune response and neuroinflammation. This prompted us to hypothesize a link between changes in microglia homeostasis and loss of GMV. For this reason, we analyzed the number and morphology of microglial cells in the mPFC and insula. The number of neurons and their morphology was also evaluated. The primary motor cortex, where no GMV change was observed, was used as negative control. We found no differences in the number of neurons and microglia cells following heroin. However, in the same regions where reduced GMV was detected, we observed a shift towards a rounder shape and size reduction in microglia, suggestive of their homeostatic change towards a reactive state. Altogether these findings suggest that escalation of heroin intake correlates with loss of GMV in specific brain regions and that this phenomenon is linked to changes in microglial morphology.

Long access heroin self-administration significantly alters gut microbiome composition and structure

Introduction

It is well known that chronic opioid use disorder is associated with alterations in gastrointestinal (GI) function that include constipation, reduced motility, and increased bacterial translocation due to compromised gut barrier function. These signs of disrupted GI function can be associated with alterations in the gut microbiome. However, it is not known if long-access opioid self-administration has effects on the gut microbiome.

Methods

We used 16S rRNA gene sequencing to investigate the gut microbiome in three independent cohorts (N=40 for each) of NIH heterogeneous stock rats before onset of long-access heroin self-administration (i.e., naïve status), at the end of a 15-day period of self-administration, and after post-extinction reinstatement. Measures of microbial α- and β-diversity were evaluated for all phases. High-dimensional class comparisons were carried out with MaAsLin2. PICRUSt2 was used for predicting functional pathways impacted by heroin based on marker gene sequences.

Results

Community α-diversity was not altered by heroin at any of the three phases by comparison to saline-yoked controls. Analyses of β-diversity showed that the heroin and saline-yoked groups clustered significantly apart from each other using the Bray-Curtis (community structure) index. Heroin caused significant alterations at the ASV level at the self-administration and extinction phases. At the phylum level, the relative abundance of Firmicutes was increased at the self-administration phase. Deferribacteres was decreased in heroin whereas Patescibacteria was increased in heroin at the extinction phase. Potential biomarkers for heroin emerged from the MaAsLin2 analysis. Bacterial metabolomic pathways relating to degradation of carboxylic acids, nucleotides, nucleosides, carbohydrates, and glycogen were increased by heroin while pathways relating to biosynthesis of vitamins, propionic acid, fatty acids, and lipids were decreased.

Discussion

These findings support the view that long access heroin self-administration significantly alters the structure of the gut microbiome by comparison to saline-yoked controls. Inferred metabolic pathway alterations suggest the development of a microbial imbalance favoring gut inflammation and energy expenditure. Potential microbial biomarkers and related functional pathways likely invoked by heroin self-administration could be targets for therapeutic intervention.

Opioid trail: Tracking contributions to opioid use disorder from host genetics to the gut microbiome

Opioid use disorder (OUD) is a worldwide public health crisis with few effective treatment options. Traditional genetics and neuroscience approaches have provided knowledge about biological mechanisms that contribute to OUD-related phenotypes, but the complexity and magnitude of effects in the brain and body remain poorly understood. The gut-brain axis has emerged as a promising target for future therapeutics for several psychiatric conditions, so characterizing the relationship between host genetics and the gut microbiome in the context of OUD will be essential for development of novel treatments. In this review, we describe evidence that interactions between host genetics, the gut microbiome, and immune signaling likely play a key role in mediating opioid-related phenotypes. Studies in humans and model organisms consistently demonstrated that genetic background is a major determinant of gut microbiome composition. Furthermore, the gut microbiome is susceptible to environmental influences such as opioid exposure. Additional work focused on gene by microbiome interactions will be necessary to gain improved understanding of their effects on OUD-related behaviors.

Genome-wide association study of delay discounting in Heterogenous Stock rats

Delay discounting refers to the behavioral tendency to devalue rewards as a function of their delay in receipt. Heightened delay discounting has been associated with substance use disorders, as well as multiple co-occurring psychopathologies. Genetic studies in humans and animal models have established that delay discounting is a heritable trait, but only a few specific genes have been associated with delay discounting. Here, we aimed to identify novel genetic loci associated with delay discounting through a genome-wide association study (GWAS) using Heterogenous Stock rats, a genetically diverse outbred population derived from eight inbred founder strains. We assessed delay discounting in 650 male and female rats using an adjusting amount procedure in which rats chose between smaller immediate sucrose rewards or a larger reward at variable delays. Preference switch points were calculated for each rat and both exponential and hyperbolic functions were fitted to these indifference points. Area under the curve (AUC) and the discounting parameter k of both functions were used as delay discounting measures. GWAS for AUC, exponential k, and indifference points for a short delay identified significant loci on chromosomes 20 and 14. The gene Slc35f1, which encodes a member of the solute carrier family of nucleoside sugar transporters, was the only gene within the chromosome 20 locus. That locus also contained an eQTL for Slc35f1, suggesting that heritable differences in the expression of that gene might be responsible for the association with behavior. The gene Adgrl3, which encodes a member of the latrophilin family of G-protein coupled receptors, was the only gene within the chromosome 14 locus. These findings implicate novel genes in delay discounting and highlight the need for further exploration.

Systems genetics approaches for understanding complex traits with relevance for human disease

Quantitative traits are often complex because of the contribution of many loci, with further complexity added by environmental factors. In medical research, systems genetics is a powerful approach for the study of complex traits, as it integrates intermediate phenotypes, such as RNA, protein, and metabolite levels, to understand molecular and physiological phenotypes linking discrete DNA sequence variation to complex clinical and physiological traits. The primary purpose of this review is to describe some of the resources and tools of systems genetics in humans and rodent models, so that researchers in many areas of biology and medicine can make use of the data.

Single-nucleus genomics in outbred rats with divergent cocaine addiction-like behaviors reveals changes in amygdala GABAergic inhibition

The amygdala processes positive and negative valence and contributes to addiction, but the cell-type-specific gene regulatory programs involved are unknown. We generated an atlas of single-nucleus gene expression and chromatin accessibility in the amygdala of outbred rats with high and low cocaine addiction-like behaviors following prolonged abstinence. Differentially expressed genes between the high and low groups were enriched for energy metabolism across cell types. Rats with high addiction index (AI) showed increased relapse-like behaviors and GABAergic transmission in the amygdala. Both phenotypes were reversed by pharmacological inhibition of the glyoxalase 1 enzyme, which metabolizes methylglyoxal-a GABAA receptor agonist produced by glycolysis. Differences in chromatin accessibility between high and low AI rats implicated pioneer transcription factors in the basic helix-loop-helix, FOX, SOX and activator protein 1 families. We observed opposite regulation of chromatin accessibility across many cell types. Most notably, excitatory neurons had greater accessibility in high AI rats and inhibitory neurons had greater accessibility in low AI rats.

Low prepulse inhibition predicts lower social interaction, impaired spatial working memory, reference memory and cognitive flexibility in genetically heterogeneous rats

The "Genetically Heterogeneous National Institutes of Health (NIHHS)" stock rat (hereafter HS) shows a wide phenotypic variation, as a result of having been derived from eight inbred rat strains. Thus, these rats may be a conceivable parallel model of a healthy human sample. In order to evaluate whether HS rats have face validity as an animal model of schizophrenia-relevant features, it should be demonstrated that they present behavioural traits that may model negative and cognitive symptoms of the disorder. Previous studies on HS rats have shown that prepulse inhibition (PPI, a measure of sensorimotor gating processes), which is impaired in schizophrenic patients, is correlated with their working memory performance. In this study, we evaluated whether low PPI in the HS stock rat predicts impairments of spatial working memory (SWM), spatial reference memory and cognitive flexibility in the Morris water maze (MWM) test, and we evaluated HS rats for social interaction (SI) in a social investigation task. HS rats were stratified into 2 different groups according to their PPI scores, i.e. low- and high-PPI. In the SI task, low-PPI rats showed decreased social behaviour compared to high-PPI rats. In addition, relative to high-PPI HS rats, the low-PPI group displayed poorer SWM performance, impaired cognitive flexibility (in a reversal task) and worsened long-term spatial memory. Such differential behaviours in social and cognitive paradigms provide evidence on the face validity of low-PPI HS rats as a model of negative-like and cognitive schizophrenia-relevant traits.

Pierre CL, Tripi JA, Richards JB, Wang T, Chen H, Flagel SB, Meyer P, Robinson TE, Solberg Woods LC, Palmer AA. An exponential increase in QTL detection with an increased sample size

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

Reward maximization assessed using a sequential patch depletion task in a large sample of heterogeneous stock rats

Choice behavior requires animals to evaluate both short- and long-term advantages and disadvantages of all potential alternatives. Impulsive choice is traditionally measured in laboratory tasks by utilizing delay discounting (DD), a paradigm that offers a choice between a smaller immediate reward, or a larger more delayed reward. This study tested a large sample of Heterogeneous Stock (HS) male (n = 896) and female (n = 898) rats, part of a larger genetic study, to investigate whether measures of reward maximization overlapped with traditional models of delay discounting via the patch depletion model using a Sequential Patch Depletion procedure. In this task, rats were offered a concurrent choice between two water "patches" and could elect to "stay" in the current patch or "leave" for an alternative patch. Staying in the current patch resulted in decreasing subsequent reward magnitudes, whereas the choice to leave a patch was followed by a delay and a resetting to the maximum reward magnitude. Based on the delay in a given session, different visit durations were necessary to obtain the maximum number of rewards. Visit duration may be analogous to an indifference point in traditional DD tasks. Males and females did not significantly differ on traditional measures of DD (e.g. delay gradient; AUC). When examining measures of patch utilization, females made fewer patch changes at all delays and spent more time in the patch before leaving for the alternative patch compared to males. Consistent with this, there was some evidence that females deviated from reward maximization more than males. However, when controlling for body weight, females had a higher normalized rate of reinforcement than males. Measures of reward maximization were only weakly associated with traditional DD measures and may represent distinctive underlying processes. Taken together, females performance differed from males with regard to reward maximization that were not observed utilizing traditional measures of DD, suggesting that the patch depletion model was more sensitive to modest sex differences when compared to traditional DD measures in a large sample of HS rats.

Reward Maximization Assessed Using a Sequential Patch Depletion Task in a Large Sample of Heterogeneous Stock Rats

Choice behavior requires animals to evaluate both short- and long-term advantages and disadvantages of all potential alternatives. Impulsive choice is traditionally measured in laboratory tasks by utilizing delay discounting (DD), a paradigm that offers a choice between a smaller immediate reward, or a larger more delayed reward. This study tested a large sample of Heterogeneous Stock (HS) male (n = 896) and female (n = 898) rats, part of a larger genetic study, to investigate whether measures of reward maximization overlapped with traditional models of delay discounting via the patch depletion model using a Sequential Patch Depletion procedure. In this task, rats were offered a concurrent choice between two water "patches" and could elect to "stay" in the current patch or "leave" for an alternative patch. Staying in the current patch resulted in decreasing subsequent reward magnitudes, whereas the choice to leave a patch was followed by a delay and a resetting to the maximum reward magnitude. Based on the delay in a given session, different visit durations were necessary to obtain the maximum number of rewards. Visit duration may be analogous to an indifference point in traditional DD tasks. While differences in traditional DD measures (e.g., delay gradient) have been detected between males and females, these effects were small and inconsistent. However, when examining measures of reward maximization, females made fewer patch changes at all delays and spent more time in the patch before leaving for the alternative patch compared to males. This pattern of choice resulted in males having a higher rate of reinforcement than females. Consistent with this, there was some evidence that females deviated from the optimal more, leading to less reward. Measures of reward maximization were only weakly associated with traditional DD measures and may represent distinctive underlying processes. Taken together, females performance differed from males with regard to reward maximization that were not observed utilizing traditional measures of DD, suggesting that the patch depletion model was more sensitive to modest sex differences when compared to traditional DD measures in a large sample of HS rats.

Pierre CL, Bimschleger H, Gao J, Cheng R, Mohammadi P, Chen H, Palmer AA, Polesskaya O, Jablonski MM. Genome-wide association study finds multiple loci associated with intraocular pressure in HS rats

Elevated intraocular pressure (IOP) is influenced by environmental and genetic factors. Increased IOP is a major risk factor for most types of glaucoma, including primary open angle glaucoma (POAG). Investigating the genetic basis of IOP may lead to a better understanding of the molecular mechanisms of POAG. The goal of this study was to identify genetic loci involved in regulating IOP using outbred heterogeneous stock (HS) rats. HS rats are a multigenerational outbred population derived from eight inbred strains that have been fully sequenced. This population is ideal for a genome-wide association study (GWAS) owing to the accumulated recombinations among well-defined haplotypes, the relatively high allele frequencies, the accessibility to a large collection of tissue samples, and the large allelic effect size compared to human studies. Both male and female HS rats (N = 1,812) were used in the study. Genotyping-by-sequencing was used to obtain ∼3.5 million single nucleotide polymorphisms (SNP) from each individual. SNP heritability for IOP in HS rats was 0.32, which agrees with other studies. We performed a GWAS for the IOP phenotype using a linear mixed model and used permutation to determine a genome-wide significance threshold. We identified three genome-wide significant loci for IOP on chromosomes 1, 5, and 16. Next, we sequenced the mRNA of 51 whole eye samples to find cis-eQTLs to aid in identification of candidate genes. We report 5 candidate genes within those loci: Tyr, Ctsc, Plekhf2, Ndufaf6 and Angpt2. Tyr, Ndufaf6 and Angpt2 genes have been previously implicated by human GWAS of IOP-related conditions. Ctsc and Plekhf2 genes represent novel findings that may provide new insight into the molecular basis of IOP. This study highlights the efficacy of HS rats for investigating the genetics of elevated IOP and identifying potential candidate genes for future functional testing.

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.

Novelty-induced locomotor behavior predicts heroin addiction vulnerability in male, but not female, rats

RATIONALE

The ongoing rise in opioid use disorder (OUD) has made it imperative to better model the individual variation within the human population that contributes to OUD vulnerability. Using animal models that capture such variation can be a useful tool. Individual variation in novelty-induced locomotion is predictive of substance use disorder (SUD) propensity. In this model, rats are characterized as high-responders (HR) or low-responders (LR) using a median split based on distance travelled during a locomotor test, and HR rats are generally found to exhibit a more SUD vulnerable behavioral phenotype.

OBJECTIVES

The HR/LR model has commonly been used to assess behaviors in male rats using psychostimulants, with limited knowledge of the predictive efficacy of this model in females or the use of an opioid as the reward. In the current study, we assessed several behaviors across the different phases of drug addiction (heroin taking, refraining, and seeking) in over 500 male and female heterogeneous stock rats run at two geographically separate locations. Rats were characterized as HRs or LRs within each sex for analysis.

RESULTS

Overall, females exhibit a more OUD vulnerable phenotype relative to males. Additionally, the HR/LR model was predictive of OUD-like behaviors in male, but not female rats. Furthermore, phenotypes did not differ in anxiety-related behaviors, reacquisition of heroin-taking, or punished heroin-taking behavior in either sex.

CONCLUSIONS

These results emphasize the importance of assessing females in models of individual variation in SUD and highlight limitations in using the HR/LR model to assess OUD propensity.

Keratinocyte-associated protein 3 plays a role in body weight and adiposity with differential effects in males and females

Despite the obesity crisis in the United States, the underlying genetics are poorly understood. Our lab previously identified Keratinocyte-associated protein 3, Krtcap3, as a candidate gene for adiposity through a genome-wide association study in outbred rats, where increased liver expression of Krtcap3 correlated with decreased fat mass. Here we seek to confirm that Krtcap3 expression affects adiposity traits. To do so, we developed an in vivo whole-body Krtcap3 knock-out (KO) rat model. Wild-type (WT) and KO rats were placed onto a high-fat (HFD) or low-fat diet (LFD) at 6 weeks of age and were maintained on diet for 13 weeks, followed by assessments of metabolic health. We hypothesized that Krtcap3-KO rats will have increased adiposity and a worsened metabolic phenotype relative to WT. We found that KO male and female rats have significantly increased body weight versus WT, with the largest effect in females on a HFD. KO females also ate more and had greater adiposity, but were more insulin sensitive than WT regardless of diet condition. Although KO males weighed more than WT under both diet conditions, there were no differences in eating behavior or fat mass. Interestingly, KO males on a HFD were more insulin resistant than WT. This study confirms that Krtcap3 plays a role in body weight regulation and demonstrates genotype- and sex-specific effects on food intake, adiposity, and insulin sensitivity. Future studies will seek to better understand these sex differences, the role of diet, and establish a mechanism for Krtcap3 in obesity.

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.

Glucocorticoid Receptor-Regulated Enhancers Play a Central Role in the Gene Regulatory Networks Underlying Drug Addiction

Substance abuse and addiction represent a significant public health problem that impacts multiple dimensions of society, including healthcare, the economy, and the workforce. In 2021, over 100,000 drug overdose deaths were reported in the US, with an alarming increase in fatalities related to opioids and psychostimulants. Understanding the fundamental gene regulatory mechanisms underlying addiction and related behaviors could facilitate more effective treatments. To explore how repeated drug exposure alters gene regulatory networks in the brain, we combined capped small (cs)RNA-seq, which accurately captures nascent-like initiating transcripts from total RNA, with Hi-C and single nuclei (sn)ATAC-seq. We profiled initiating transcripts in two addiction-related brain regions, the prefrontal cortex (PFC) and the nucleus accumbens (NAc), from rats that were never exposed to drugs or were subjected to prolonged abstinence after oxycodone or cocaine intravenous self-administration (IVSA). Interrogating over 100,000 active transcription start regions (TSRs) revealed that most TSRs had hallmarks of bonafide enhancers and highlighted the KLF/SP1, RFX, and AP1 transcription factors families as central to establishing brain-specific gene regulatory programs. Analysis of rats with addiction-like behaviors versus controls identified addiction-associated repression of transcription at regulatory enhancers recognized by nuclear receptor subfamily 3 group C (NR3C) factors, including glucocorticoid receptors. Cell-type deconvolution analysis using snATAC-seq uncovered a potential role of glial cells in driving the gene regulatory programs associated with addiction-related phenotypes. These findings highlight the power of advanced transcriptomics methods to provide insight into how addiction perturbs gene regulatory programs in the brain.

Ten Points to Improve Reproducibility and Translation of Animal Research

Findings from animal experiments are often difficult to transfer to humans. In this perspective article I discuss two questions. First, why are the results of animal experiments often so difficult to transfer to humans? And second, what can be done to improve translation from animal experiments to humans? Translation failures are often the result of poor methodology. It is not merely the fact that low statistical power of basic and preclinical studies undermine a "real effect," but the accuracy with which data from animal studies are collected and described, and the resulting robustness of the data is generally very low and often does not allow translation to a much more heterogeneous human condition. Equally important is the fact that the vast majority of publications in the biomedical field in the last few decades have reported positive findings and have thus generated a knowledge bias. Further contributions to reproducibility and translation failures are discussed in this paper, and 10 points of recommendation to improve reproducibility and translation are outlined. These recommendations are: (i) prior to planning an actual study, a systematic review or potential preclinical meta-analysis should be considered. (ii) An a priori power calculation should be carried out. (iii) The experimental study protocol should be pre-registered. (iv) The execution of the study should be in accordance with the most recent ARRIVE guidelines. (v) When planning the study, the generalizability of the data to be collected should also be considered (e.g., sex or age differences). (vi) "Method-hopping" should be avoided, meaning that it is not necessary to use the most advanced technology but rather to have the applied methodology under control. (vii) National or international networks should be considered to carry out multicenter preclinical studies or to obtain convergent evidence. (viii) Animal models that capture DSM-5 or ICD-11 criteria should be considered in the context of research on psychiatric disorders. (ix) Raw data of publication should be made publicly available and should be in accordance with the FAIR Guiding Principles for scientific data management. (x) Finally, negative findings should be published to counteract publication bias. The application of these 10 points of recommendation, especially for preclinical confirmatory studies but also to some degree for exploratory studies, will ultimately improve the reproducibility and translation of animal research.

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.

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.

Individual differences in addiction-like behaviors and choice between cocaine versus food in Heterogeneous Stock rats

RATIONALE AND OBJECTIVES

Recent studies reported that when given a mutually exclusive choice between cocaine and palatable food, most rats prefer the non-drug reward over cocaine. However, these studies used rat strains with limited genetic and behavioral diversity. Here, we used a unique outbred strain of rats (Heterogeneous Stock, HS) that mimic the genetic variability of humans.

METHODS

We first identified individual differences in addiction-like behaviors (low and high). Next, we tested choice between cocaine and palatable food using a discrete choice procedure. We characterized the individual differences using an addiction score that incorporates key features of addiction: escalated intake, highly motivated responding (progressive ratio), and responding despite adverse consequences (footshock punishment). We assessed food versus cocaine choice at different drug-free days (without pre-choice cocaine self-administration) during acquisition of cocaine self-administration or after escalation of cocaine self-administration. We also assessed drug versus food choice immediately after 1-, 2-, or 6-h cocaine self-administration.

RESULTS

Independent of the addiction score, without pre-choice cocaine (1 or more abstinence days), HS rats strongly preferred the palatable food over cocaine, even if the food reward was delayed or its size was reduced. However, rats with high but not low addiction score modestly increased cocaine choice immediately after 1-, 2-, or 6-h cocaine self-administration.

CONCLUSIONS

Like other strains, HS rats strongly prefer palatable food over cocaine. Individual differences in addiction score were associated with increased drug choice in the presence but not absence (abstinence) of cocaine. The HS strain may be useful in studies on mechanisms of addiction vulnerability.

The Cocaine and Oxycodone Biobanks, Two Repositories from Genetically Diverse and Behaviorally Characterized Rats for the Study of Addiction

The rat oxycodone and cocaine biobanks contain samples that vary by genotypes (by using genetically diverse genotyped HS rats), phenotypes (by measuring addiction-like behaviors in an advanced SA model), timepoints (samples are collected longitudinally before, during, and after SA, and terminally at three different timepoints in the addiction cycle: intoxication, withdrawal, and abstinence or without exposure to drugs through age-matched naive rats), samples collected (organs, cells, biofluids, feces), preservation (paraformaldehyde-fixed, snap-frozen, or cryopreserved) and application (proteomics, transcriptomics, microbiomics, metabolomics, epigenetics, anatomy, circuitry analysis, biomarker discovery, etc.Substance use disorders (SUDs) are pervasive in our society and have substantial personal and socioeconomical costs. A critical hurdle in identifying biomarkers and novel targets for medication development is the lack of resources for obtaining biological samples with a detailed behavioral characterization of SUD. Moreover, it is nearly impossible to find longitudinal samples. As part of two ongoing large-scale behavioral genetic studies in heterogeneous stock (HS) rats, we have created two preclinical biobanks using well-validated long access (LgA) models of intravenous cocaine and oxycodone self-administration (SA) and comprehensive characterization of addiction-related behaviors. The genetic diversity in HS rats mimics diversity in the human population and includes individuals that are vulnerable or resilient to compulsive-like responding for cocaine or oxycodone. Longitudinal samples are collected throughout the experiment, before exposure to the drug, during intoxication, acute withdrawal, and protracted abstinence, and include naive, age-matched controls. Samples include, but are not limited to, blood plasma, feces and urine, whole brains, brain slices and punches, kidney, liver, spleen, ovary, testis, and adrenal glands. Three preservation methods (fixed in formaldehyde, snap-frozen, or cryopreserved) are used to facilitate diverse downstream applications such as proteomics, metabolomics, transcriptomics, epigenomics, microbiomics, neuroanatomy, biomarker discovery, and other cellular and molecular approaches. To date, >20,000 samples have been collected from over 1000 unique animals and made available free of charge to non-profit institutions through https://www.cocainebiobank.org/ and https://www.oxycodonebiobank.org/.

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.

Integration of evidence across human and model organism studies: A meeting report

The National Institute on Drug Abuse and Joint Institute for Biological Sciences at the Oak Ridge National Laboratory hosted a meeting attended by a diverse group of scientists with expertise in substance use disorders (SUDs), computational biology, and FAIR (Findability, Accessibility, Interoperability, and Reusability) data sharing. The meeting's objective was to discuss and evaluate better strategies to integrate genetic, epigenetic, and 'omics data across human and model organisms to achieve deeper mechanistic insight into SUDs. Specific topics were to (a) evaluate the current state of substance use genetics and genomics research and fundamental gaps, (b) identify opportunities and challenges of integration and sharing across species and data types, (c) identify current tools and resources for integration of genetic, epigenetic, and phenotypic data, (d) discuss steps and impediment related to data integration, and (e) outline future steps to support more effective collaboration-particularly between animal model research communities and human genetics and clinical research teams. This review summarizes key facets of this catalytic discussion with a focus on new opportunities and gaps in resources and knowledge on SUDs.

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.

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.

Mitochondrial health is enhanced in rats with higher vs. lower intrinsic exercise capacity and extended lifespan

The intrinsic aerobic capacity of an organism is thought to play a role in aging and longevity. Maximal respiratory rate capacity, a metabolic performance measure, is one of the best predictors of cardiovascular- and all-cause mortality. Rats selectively bred for high-(HCR) vs. low-(LCR) intrinsic running-endurance capacity have up to 31% longer lifespan. We found that positive changes in indices of mitochondrial health in cardiomyocytes (respiratory reserve, maximal respiratory capacity, resistance to mitochondrial permeability transition, autophagy/mitophagy, and higher lipids-over-glucose utilization) are uniformly associated with the extended longevity in HCR vs. LCR female rats. Cross-sectional heart metabolomics revealed pathways from lipid metabolism in the heart, which were significantly enriched by a select group of strain-dependent metabolites, consistent with enhanced lipids utilization by HCR cardiomyocytes. Heart–liver–serum metabolomics further revealed shunting of lipidic substrates between the liver and heart via serum during aging. Thus, mitochondrial health in cardiomyocytes is associated with extended longevity in rats with higher intrinsic exercise capacity and, probably, these findings can be translated to other populations as predictors of outcomes of health and survival.

Network-Based Discovery of Opioid Use Vulnerability in Rats Using the Bayesian Stochastic Block Model

Opioid use disorder is a psychological condition that affects over 200,000 people per year in the U.S., causing the Centers for Disease Control and Prevention to label the crisis as a rapidly spreading public health epidemic. The behavioral relationship between opioid exposure and development of opioid use disorder (OUD) varies greatly between individuals, implying existence of sup-populations with varying degrees of opioid vulnerability. However, effective pre-clinical identification of these sub-populations remains challenging due to the complex multivariate measurements employed in animal models of OUD. In this study, we propose a novel non-linear network-based data analysis workflow that employs seven behavioral traits to identify opioid use sub-populations and assesses contributions of behavioral variables to opioid vulnerability and resiliency. Through this analysis workflow we determined how behavioral variables across heroin taking, refraining and seeking interact with one another to identify potentially heroin resilient and vulnerable behavioral sub-populations. Data were collected from over 400 heterogeneous stock rats in two geographically distinct locations. Rats underwent heroin self-administration training, followed by a progressive ratio and heroin-primed reinstatement test. Next, rats underwent extinction training and a cue-induced reinstatement test. To enter the analysis workflow, we integrated data from different cohorts of rats and removed possible batch effects. We then constructed a rat-rat similarity network based on their behavioral patterns and implemented community detection on this similarity network using a Bayesian degree-corrected stochastic block model to uncover sub-populations of rats with differing levels of opioid vulnerability. We identified three statistically distinct clusters corresponding to distinct behavioral sub-populations, vulnerable, resilient and intermediate for heroin use, refraining and seeking. We implement this analysis workflow as an open source R package, named mlsbm.

Modeling epistasis in mice and yeast using the proportion of two or more distinct genetic backgrounds: Evidence for "polygenic epistasis"

Background

The majority of quantitative genetic models used to map complex traits assume that alleles have similar effects across all individuals. Significant evidence suggests, however, that epistatic interactions modulate the impact of many alleles. Nevertheless, identifying epistatic interactions remains computationally and statistically challenging. In this work, we address some of these challenges by developing a statistical test for polygenic epistasis that determines whether the effect of an allele is altered by the global genetic ancestry proportion from distinct progenitors.

Results

We applied our method to data from mice and yeast. For the mice, we observed 49 significant genotype-by-ancestry interaction associations across 14 phenotypes as well as over 1,400 Bonferroni-corrected genotype-by-ancestry interaction associations for mouse gene expression data. For the yeast, we observed 92 significant genotype-by-ancestry interactions across 38 phenotypes. Given this evidence of epistasis, we test for and observe evidence of rapid selection pressure on ancestry specific polymorphisms within one of the cohorts, consistent with epistatic selection.

Conclusions

Unlike our prior work in human populations, we observe widespread evidence of ancestry-modified SNP effects, perhaps reflecting the greater divergence present in crosses using mice and yeast.

Many statistical tests which link genetic markers in the genome to differences in traits rely on the assumption that the same polymorphism will have identical effects in different individuals. However, there is substantial evidence indicating that this is not the case. Epistasis is the phenomenon in which multiple polymorphisms interact with one another to amplify or negate each other’s effects on a trait. We hypothesized that individual SNP effects could be changed in a polygenic manner, such that the proportion of as genetic ancestry, rather than specific markers, might be used to capture epistatic interactions. Motivated by this possibility, we develop a new statistical test that allowed us to examine the genome to identify polymorphisms which have different effects depending on the ancestral makeup of each individual. We use our test in two different populations of inbred mice and a yeast panel and demonstrate that these sorts of variable effect polymorphisms exist in 14 different physical traits in mice and 38 phenotypes in yeast as well as in murine gene expression. We use the term “polygenic epistasis” to distinguish these interactions from the more conventional two- or multi-locus interactions.

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

OBJECTIVE

Obesity is influenced by genetic and environmental factors. Despite the success of human genome-wide association studies, the specific genes that confer obesity remain largely unknown. The objective of this study was to use outbred rats to identify the genetic loci underlying obesity and related morphometric and metabolic traits.

METHODS

This study measured obesity-relevant traits, including body weight, body length, BMI, fasting glucose, and retroperitoneal, epididymal, and parametrial fat pad weight in 3,173 male and female adult N/NIH heterogeneous stock (HS) rats across three institutions, providing data for the largest rat genome-wide association study to date. Genetic loci were identified using a linear mixed model to account for the complex family relationships of the HS and using covariates to account for differences among the three phenotyping centers.

RESULTS

This study identified 32 independent loci, several of which contained only a single gene (e.g., Epha5, Nrg1, Klhl14) or obvious candidate genes (e.g., Adcy3, Prlhr). There were strong phenotypic and genetic correlations among obesity-related traits, and there was extensive pleiotropy at individual loci.

CONCLUSIONS

This study demonstrates the utility of HS rats for investigating the genetics of obesity-related traits across institutions and identify several candidate genes for future functional testing.

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.

Facilitating Complex Trait Analysis via Reduced Complexity Crosses

Genetically diverse inbred strains are frequently used in quantitative trait mapping to identify sequence variants underlying trait variation. Poor locus resolution and high genetic complexity impede variant discovery. As a solution, we explore reduced complexity crosses (RCCs) between phenotypically divergent, yet genetically similar, rodent substrains. RCCs accelerate functional variant discovery via decreasing the number of segregating variants by orders of magnitude. The simplified genetic architecture of RCCs often permit immediate identification of causal variants or rapid fine-mapping of broad loci to smaller intervals. Whole-genome sequences of substrains make RCCs possible by supporting the development of array- and targeted sequencing-based genotyping platforms, coupled with rapid genome editing for variant validation. In summary, RCCs enhance discovery-based genetics of complex traits.

Sex-dependent associations between addiction-related behaviors and the microbiome in outbred rats

BACKGROUND

Multiple factors contribute to the etiology of addiction, including genetics, sex, and a number of addiction-related behavioral traits. One behavioral trait where individuals assign incentive salience to food stimuli ("sign-trackers", ST) are more impulsive compared to those that do not ("goal-trackers", GT), as well as more sensitive to drugs and drug stimuli. Furthermore, this GT/ST phenotype predicts differences in other behavioral measures. Recent studies have implicated the gut microbiota as a key regulator of brain and behavior, and have shown that many microbiota-associated changes occur in a sex-dependent manner. However, few studies have examined how the microbiome might influence addiction-related behaviors. To this end, we sought to determine if gut microbiome composition was correlated with addiction-related behaviors determined by the GT/ST phenotype.

METHODS

Outbred male (N=101) and female (N=101) heterogeneous stock rats underwent a series of behavioral tests measuring impulsivity, attention, reward-learning, incentive salience, and locomotor response. Cecal microbiome composition was estimated using 16S rRNA gene amplicon sequencing. Behavior and microbiome were characterized and correlated with behavioral phenotypes. Robust sex differences were observed in both behavior and microbiome; further analyses were conducted within sex using the pre-established goal/sign-tracking (GT/ST) phenotype and partial least squares differential analysis (PLS-DA) clustered behavioral phenotype.

RESULTS

Overall microbiome composition was not associated to the GT/ST phenotype. However, microbial alpha diversity was significantly decreased in female STs. On the other hand, a measure of impulsivity had many significant correlations to microbiome in both males and females. Several measures of impulsivity were correlated with the genus Barnesiella in females. Female STs had notable correlations between microbiome and attentional deficient. In both males and females, many measures were correlated with the bacterial families Ruminocococcaceae and Lachnospiraceae.

CONCLUSIONS

These data demonstrate correlations between several addiction-related behaviors and the microbiome specific to sex.

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.