Search

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

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., Competing Interests: Gd, LC, MK, MB, BB, LM, DC, SS, LT, MF, AM, NV, DO, BS, JR, JL, JK, AK, SS, LS, KS, SB, ES, AA, SP, JM, CC, LC, AV, KM, PS, BL, BP, AC, OP, LS, AP, OG No competing interests declared, (© 2023, de Guglielmo, Carrette et al.)

Objective: The behavioral and diagnostic heterogeneity within human opioid use disorder (OUD) diagnosis is not readily captured in current animal models, limiting translational relevance of the mechanistic research that is conducted in experimental animals. We hypothesize that a non-linear clustering of OUD-like behavioral traits will capture population heterogeneity and yield subpopulations of OUD vulnerable rats with distinct behavioral and neurocircuit profiles., Methods: Over 900 male and female heterogeneous stock rats, a line capturing genetic and behavioral heterogeneity present in humans, were assessed for several measures of heroin use and rewarded and non-rewarded seeking behaviors. Using a non-linear stochastic block model clustering analysis, rats were assigned to OUD vulnerable, intermediate and resilient clusters. Additional behavioral tests and circuit analyses using c-fos protein activation were conducted on the vulnerable and resilient subpopulations., Results: OUD vulnerable rats exhibited greater heroin taking and seeking behaviors relative to those in the intermediate and resilient clusters. Akin to human OUD diagnosis, further vulnerable rat sub-clustering revealed subpopulations with different combinations of behavioral traits, including sex differences. Lastly, heroin cue-induced neuronal patterns of circuit activation differed between resilient and vulnerable phenotypes. Behavioral sex differences were recapitulated in patterns of circuitry activation, including males preferentially engaging extended amygdala stress circuitry, and females cortico-striatal drug cue-seeking circuitry., Conclusion: Using a non-linear clustering approach in rats, we captured behavioral diagnostic heterogeneity reflective of human OUD diagnosis. OUD vulnerability and resiliency were associated with distinct neuronal activation patterns, posing this approach as a translational tool in assessing neurobiological mechanisms underpinning OUD.

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 and multiple co-occurring psychopathologies. Human and animal genetic studies have established that delay discounting is heritable, but only a few associated genes have been identified. We aimed to identify novel genetic loci associated with delay discounting through a genome-wide association study (GWAS) using Heterogeneous Stock (HS) rats, a genetically diverse outbred population derived from eight inbred founder strains. We assessed delay discounting in 650 male and female HS rats using an adjusting amount procedure in which rats chose between smaller immediate sucrose rewards or a larger reward at various delays. Preference switch points were calculated 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 one indifference point identified significant loci on chromosomes 20 and 14. The gene Slc35f1, which encodes a member of the solute carrier family, was the sole gene within the chromosome 20 locus. That locus also contained an eQTL for Slc35f1, suggesting that heritable differences in the expression might be responsible for the association with behavior. Adgrl3, which encodes a latrophilin subfamily G-protein coupled receptor, was the sole gene within the chromosome 14 locus. These findings implicate novel genes in delay discounting and highlight the need for further exploration., (© 2024 The Author(s). Genes, Brain and Behavior published by International Behavioural and Neural Genetics Society and John Wiley & Sons Ltd.)

The increased prevalence of opioid use disorder (OUD) makes it imperative to disentangle the biological mechanisms contributing to individual differences in OUD vulnerability. OUD shows strong heritability, however genetic variants contributing toward vulnerability remain poorly defined. We performed a genome-wide association study using over 850 male and female heterogeneous stock (HS) rats to identify genes underlying behaviors associated with OUD such as nociception, as well as heroin-taking, extinction and seeking behaviors. By using an animal model of OUD, we were able to identify genetic variants associated with distinct OUD behaviors while maintaining a uniform environment, an experimental design not easily achieved in humans. Furthermore, we used a novel non-linear network-based clustering approach to characterize rats based on OUD vulnerability to assess genetic variants associated with OUD susceptibility. Our findings confirm the heritability of several OUD-like behaviors, including OUD susceptibility. Additionally, several genetic variants associated with nociceptive threshold prior to heroin experience, heroin consumption, escalation of intake, and motivation to obtain heroin were identified. Tom1 , a microglial component, was implicated for nociception. Several genes involved in dopaminergic signaling, neuroplasticity and substance use disorders, including Brwd1 , Pcp4, Phb1l2 and Mmp15 were implicated for the heroin traits. Additionally, an OUD vulnerable phenotype was associated with genetic variants for consumption and break point, suggesting a specific genetic contribution for OUD-like traits contributing to vulnerability. Together, these findings identify novel genetic markers related to the susceptibility to OUD-relevant behaviors in HS rats.

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., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Inc. All rights reserved.)

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

Organisms must regulate their behavior flexibly in the face of environmental challenges. Failure can lead to a host of maladaptive behavioral traits associated with a range of neuropsychiatric disorders, including attention deficit hyperactivity disorder, autism, and substance use disorders. This maladaptive dysregulation of behavior is influenced by genetic and environmental factors. For example, environmental enrichment produces beneficial neurobehavioral effects in animal models of such disorders. The present study determined the effects of environmental enrichment on a range of measures related to behavioral regulation using a large cohort of male, outbred heterogeneous stock (HS) rats as subjects. Subjects were reared from late adolescence onwards either in pairs in standard housing with minimal enrichment (n = 200) or in groups of 16 in a highly enriched environment consisting of a large multi-level cage filled with toys, running wheels, and shelters (n = 64). Rats were subjected to a battery of tests, including: (i) locomotor response to novelty, (ii) light reinforcement, (iii) social reinforcement, (iv) reaction time, (v) a patch-depletion foraging test, (vi) Pavlovian conditioned approach, (vii) conditioned reinforcement, and (viii) cocaine conditioned cue preference. Results indicated that rats housed in the enriched environment were able to filter out irrelevant stimuli more effectively and thereby regulate their behavior more efficiently than standard-housing rats. The dramatic impact of environmental enrichment suggests that behavioral studies using standard housing conditions may not generalize to more complex environments that may be more ethologically relevant., (© 2024. The Author(s).)

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., (© 2023. The Author(s).)

We previously identified Keratinocyte-associated protein 3 , Krtcap3 , as a novel adiposity gene, but subsequently found that its impact on adiposity may depend on environmental stress. To more thoroughly understand the connection between Krtcap3 , adiposity, and stress, we exposed wild-type (WT) and Krtcap3 knock-out (KO) rats to chronic stress then measured adiposity and behavioral outcomes. We found that KO rats displayed lower basal stress than WT rats under control conditions and exhibited metabolic and behavioral responses to chronic stress exposure. Specifically, stress-exposed KO rats gained more weight, consumed more food when socially isolated, and displayed more anxiety-like behaviors relative to control KO rats. Meanwhile, there were minimal differences between control and stressed WT rats. At study conclusion stress-exposed KO rats had increased corticosterone (CORT) relative to control KO rats with no differences between WT rats. In addition, KO rats, independent of prior stress exposure, had an increased CORT response to removal of their cage-mate (psychosocial stress), which was only seen in WT rats when exposed to chronic stress. Finally, we found differences in expression of the glucocorticoid receptor, Nr3c1 , in the pituitary and colon between control and stress-exposed KO rats that were not present in WT rats. These data support that Krtcap3 expression affects stress response, potentially via interactions with Nr3c1 , with downstream effects on adiposity and behavior. Future work is necessary to more thoroughly understand the role of Krtcap3 in the stress response., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Szalanczy, Fitzpatrick, Beeson, Bui, Dyson, Eller, Landry, Scott, Grzybowski, Klotz, Geurts, Weiner, Redei and Solberg Woods.)

Obesity and major depressive disorder (MDD) are both significant health issues that have been increasing in prevalence and are associated with multiple comorbidities. Obesity and MDD have been shown to be bidirectionally associated, and they are both influenced by genetics and environmental factors. However, the molecular mechanisms that link these two diseases are not yet fully understood. It is possible that these diseases are connected through the actions of the cAMP/protein kinase A (PKA) pathway. Within this pathway, adenylate cyclase 3 ( Adcy3 ) has emerged as a key player in both obesity and MDD. Numerous genetic variants in Adcy3 have been identified in humans in association with obesity. Rodent knockout studies have also validated the importance of this gene for energy homeostasis. Furthermore, Adcy3 has been identified as a top candidate gene and even a potential blood biomarker for MDD. Adcy3 and the cAMP/PKA pathway may therefore serve as an important genetic and functional link between these two diseases. In this mini-review, we discuss the role of both Adcy3 and the cAMP/PKA pathway, including specific genetic mutations, in both diseases. Understanding the role that Adcy3 mutations play in obesity and MDD could open the door for precision medicine approaches and treatments for both diseases that target this gene.

We previously identified keratinocyte-associated protein 3, Krtcap3 , as an obesity-related gene in female rats where a whole body Krtcap3 knockout (KO) led to increased adiposity compared to wild-type (WT) controls when fed a high-fat diet (HFD). We sought to replicate this work to better understand the function of Krtcap3 but were unable to reproduce the adiposity phenotype. In the current work, WT female rats ate more compared to WT in the prior study, with corresponding increases in body weight and fat mass, while there were no changes in these measures in KO females between the studies. The prior study was conducted before the COVID-19 pandemic, while the current study started after initial lockdown orders and was completed during the pandemic in a generally less stressful environment. We hypothesize that the environmental changes impacted stress levels and may explain the failure to replicate our results. Analysis of corticosterone (CORT) at euthanasia showed a significant study-by-genotype interaction where WT had significantly higher CORT relative to KO in study 1 , with no differences in study 2 . These data suggest that decreasing Krtcap3 expression may alter the environmental stress response to influence adiposity. We also found that KO rats in both studies, but not WT, experienced a dramatic increase in CORT after their cage mate was removed, suggesting a separate connection to social behavioral stress. Future work is necessary to confirm and elucidate the finer mechanisms of these relationships, but these data indicate the possibility of Krtcap3 as a novel stress gene. NEW & NOTEWORTHY Obesity is linked to both genetics and environmental factors such as stress. Krtcap3 has previously been identified as a gene associated with adiposity, and our work here demonstrates that environmental stress may influence the role of Krtcap3 on both food intake and adiposity. Obesity is strongly influenced by stress in humans, so the identification of novel genes that link stress and obesity will greatly advance our understanding of the disease.

A vexing observation in genome-wide association studies (GWASs) is that parallel analyses in different species may not identify orthologous genes. Here, we demonstrate that cross-species translation of GWASs can be greatly improved by an analysis of co-localization within molecular networks. Using body mass index (BMI) as an example, we show that the genes associated with BMI in humans lack significant agreement with those identified in rats. However, the networks interconnecting these genes show substantial overlap, highlighting common mechanisms including synaptic signaling, epigenetic modification, and hormonal regulation. Genetic perturbations within these networks cause abnormal BMI phenotypes in mice, too, supporting their broad conservation across mammals. Other mechanisms appear species specific, including carbohydrate biosynthesis (humans) and glycerolipid metabolism (rodents). Finally, network co-localization also identifies cross-species convergence for height/body length. This study advances a general paradigm for determining whether and how phenotypes measured in model species recapitulate human biology., Competing Interests: Declaration of interests T.I. is a co-founder, member of the advisory board, and has an equity interest in Data4Cure and Serinus Biosciences. T.I. is a consultant for and has an equity interest in Ideaya BioSciences and Light Horse Therapeutics. A.A.P. is on the Scientific Advisory Board of Vivid Genomics and has stock options. The terms of these arrangements have been reviewed and approved by the University of California San Diego in accordance with its conflict of interest policies., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

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., Competing Interests: Conflicts of interest The author(s) declare no conflict of interest., (© The Author(s) 2023. Published by Oxford University Press on behalf of The Genetics Society of America. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

There is an ongoing debate about the value of animal research in psychiatry with valid lines of reasoning stating the limits of individual animal models compared to human psychiatric illnesses. Human depression is not a homogenous disorder; therefore, one cannot expect a single animal model to reflect depression heterogeneity. This limited review presents arguments that the Wistar Kyoto (WKY) rats show intrinsic depression traits. The phenotypes of WKY do not completely mirror those of human depression but clearly indicate characteristics that are common with it. WKYs present despair- like behavior, passive coping with stress, comorbid anxiety, and enhanced drug use compared to other routinely used inbred or outbred strains of rats. The commonly used tests identifying these phenotypes reflect exploratory, escape-oriented, and withdrawal-like behaviors. The WKYs consistently choose withdrawal or avoidance in novel environments and freezing behaviors in response to a challenge in these tests. The physiological response to a stressful environment is exaggerated in WKYs. Selective breeding generated two WKY substrains that are nearly isogenic but show clear behavioral differences, including that of depression-like behavior. WKY and its substrains may share characteristics of subgroups of depressed individuals with social withdrawal, low energy, weight loss, sleep disturbances, and specific cognitive dysfunction. The genomes of the WKY and WKY substrains contain variations that impact the function of many genes identified in recent human genetic studies of depression. Thus, these strains of rats share characteristics of human depression at both phenotypic and genetic levels, making them a model of depression traits., (Copyright© Bentham Science Publishers; For any queries, please email at epub@benthamscience.net.)

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., (© 2022 by the American Diabetes Association.)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Szalanczy, Goff, Seshie, Deal, Grzybowski, Klotz, Chuang Key, Geurts and Solberg Woods.)

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

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Carrette, Corral, Boomhower, Brennan, Crook, Ortez, Shankar, Simpson, Maturin, Solberg Woods, Palmer, de Guglielmo and George.)

Although obesity has been a longstanding health crisis, the genetic architecture of the disease remains poorly understood. Genome-wide association studies have identified many genomic loci associated with obesity, with genes being enriched in the brain, particularly in the hypothalamus. This points to the role of the central nervous system (CNS) in predisposition to obesity, and we emphasize here several key genes along the satiety signaling pathway involved in genetic susceptibility. Interest has also risen regarding the chronic, low-grade obesity-associated inflammation, with a growing concern toward inflammation in the hypothalamus as a precursor to obesity. Recent studies have found that genetic variation in inflammatory genes play a role in obesity susceptibility, and we highlight here several key genes. Despite the interest in the genetic variants of these pathways individually, there is a lack of research that investigates the relationship between the two. Understanding the interplay between genetic variation in obesity genes enriched in the CNS and inflammation genes will advance our understanding of obesity etiology and heterogeneity, improve genetic risk prediction analyses, and highlight new drug targets for the treatment of obesity. Additionally, this increased knowledge will assist in physician's ability to develop personalized nutrition and medication strategies for combating the obesity epidemic. Though it often seems to present universally, obesity is a highly individual disease, and there remains a need in the field to develop methods to treat at the individual level., (Copyright © 2021. Published by Elsevier Inc.)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Gunturkun, Wang, Chitre, Garcia Martinez, Holl, St. Pierre, Bimschleger, Gao, Cheng, Polesskaya, Solberg Woods, Palmer and Chen.)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Allen, Kuhn, Cannella, Crow, Roberts, Lunerti, Ubaldi, Hardiman, Solberg Woods, Ciccocioppo, Kalivas and Chung.)

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/., (Copyright © 2021 Carrette et al.)

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., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society of Toxicology. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

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

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

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.

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., (© 2020 The Obesity Society.)

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.

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., Competing Interests: The authors declare no competing interest.

Cues in the environment can elicit complex emotional states, and thereby maladaptive behavior, as a function of their ascribed value. Here we capture individual variation in the propensity to attribute motivational value to reward-cues using the sign-tracker/goal-tracker animal model. Goal-trackers attribute predictive value to reward-cues, and sign-trackers attribute both predictive and incentive value. Using chemogenetics and microdialysis, we show that, in sign-trackers, stimulation of the neuronal pathway from the prelimbic cortex (PrL) to the paraventricular nucleus of the thalamus (PVT) decreases the incentive value of a reward-cue. In contrast, in goal-trackers, inhibition of the PrL-PVT pathway increases both the incentive value and dopamine levels in the nucleus accumbens shell. The PrL-PVT pathway, therefore, exerts top-down control over the dopamine-dependent process of incentive salience attribution. These results highlight PrL-PVT pathway as a potential target for treating psychopathologies associated with the attribution of excessive incentive value to reward-cues, including addiction., Competing Interests: PC, IC, YK, AP, BK, SL, JN, SF, LS, MS, SF No competing interests declared, (© 2019, Campus et al.)

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

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

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.

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.

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

Background Histologic examination of fixed renal tissue is widely used to assess morphology and the progression of disease. Commonly reported metrics include glomerular number and injury. However, characterization of renal histology is a time-consuming and user-dependent process. To accelerate and improve the process, we have developed a glomerular localization pipeline for trichrome-stained kidney sections using a machine learning image classification algorithm. Methods We prepared 4- μ m slices of kidneys from rats of various genetic backgrounds that were subjected to different experimental protocols and mounted the slices on glass slides. All sections used in this analysis were trichrome stained and imaged in bright field at a minimum resolution of 0.92 μ m per pixel. The training and test datasets for the algorithm comprised 74 and 13 whole renal sections, respectively, totaling over 28,000 glomeruli manually localized. Additionally, because this localizer will be ultimately used for automated assessment of glomerular injury, we assessed bias of the localizer for preferentially identifying healthy or damaged glomeruli. Results Localizer performance achieved an average precision and recall of 96.94% and 96.79%, respectively, on whole kidney sections without evidence of bias for or against glomerular injury or the need for manual preprocessing. Conclusions This study presents a novel and robust application of convolutional neural nets for the localization of glomeruli in healthy and damaged trichrome-stained whole-renal section mounts and lays the groundwork for automated glomerular injury scoring., (Copyright © 2018 by the American Society of Nephrology.)

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

There is little clinical data of how hypertension may influence individuals with nephron deficiency in the context of being born with a single kidney. We recently developed a new rat model (the heterogeneous stock-derived model of unilateral renal agenesis rat) that is born with a single kidney and exhibits progressive kidney injury and decline in kidney function with age. We hypothesized that DOCA-salt would induce a greater increase in blood pressure and therefore accelerate the progression of kidney injury in rats born with a solitary kidney compared with rats that have undergone unilateral nephrectomy. Time course evaluation of blood pressure, kidney injury, and renal hemodynamics was performed in the following six groups of animals from weeks 13 to 18: 1) DOCA-treated rats with a solitary kidney (DOCA+S group), 2) placebo-treated rats with a solitary kidney, 3) DOCA-treated control rats with two kidneys (DOCA+C group), 4) placebo-treated control rats with two kidneys, 5) DOCA-treated rats with two kidneys that underwent uninephrectomy (DOCA+UNX8 group), and 6) placebo-treated rats with two kidneys that underwent uninephrectomy. DOCA+S rats demonstrated a significant rise (P < 0.05) in blood pressure (192 ± 4 mmHg), proteinuria (205 ± 31 mg/24 h), and a decline in glomerular filtration rate (600 ± 42 μl·min(-1)·g kidney weight(-1)) relative to the DOCA+UNX8 (173 ± 3 mmHg, 76 ± 26 mg/24 h, and 963 ± 36 μl·min(-1)·g kidney weight(-1)) and DOCA+C (154 ± 2 mmHg, 7 ± 1 mg/24 h, and 1,484 ± 121 μl·min(-1)·g kidney weight(-1)) groups. Placebo-treated groups showed no significant change among the three groups. An assessment of renal injury markers via real-time PCR/Western blot analysis and histological analysis was concordant with the measured physiological parameters. In summary, congenital solitary kidney rats are highly susceptible to the induction of hypertension compared with uninephrectomized rats, suggesting that low nephron endowment is an important driver of elevated blood pressure, hastening nephron injury through the transmission of elevated systemic blood pressure and thereby accelerating decline in kidney function., (Copyright © 2016 the American Physiological Society.)

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

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

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

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

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

Chronic kidney disease is a growing medical concern, with an estimated 25.6 million people in the United States exhibiting some degree of kidney injury and/or decline in kidney function. Animal models provide great insight into the study of the genetics of complex diseases. In particular, heterogeneous stock (HS) rats represent a unique genetic resource enabling rapid fine-mapping of complex traits. However, they have not been explored as a model to study renal phenotypes. To evaluate the usefulness of HS rats in the genetics of renal traits, a time course evaluation (weeks 8-40) was performed for several renal phenotypes. As expected, a large degree of variation was seen for most renal traits. By week 24, three (of 40) rats exhibited marked proteinuria that increased gradually until week 40 and ranged from 33.7 to 80.2 mg/24 h. Detailed histological analysis confirmed renal damage in these rats. In addition, several rats consistently exhibited significant hematuria (5/41). Interestingly, these rats were not the same rats that exhibited proteinuria, indicating that susceptibility to different types of kidney injury is likely segregating within the HS population. One HS rat exhibited unilateral renal agenesis (URA), which was accompanied by a significant degree of proteinuria and glomerular and tubulointerstitial injury. The parents of this HS rat were identified and bred further. Additional offspring of this pair were observed to exhibit URA at frequency between 40% and 60%. In summary, these novel data demonstrate that HS rats exhibit variation in proteinuria and other kidney-related traits, confirming that the model harbors susceptibility alleles for kidney injury and providing the basis for further genetic studies.