Celine L. St. Pierre

Background: Rats are a powerful model organism for genome-wide association studies (GWAS) of physiological and behavioral traits. Heterogeneous Stock (HS) rats, which were created by intercrossing 8 inbred rat strains, are heavily used for this purpose because of their genetic diversity and accumulated recombinations due to over 100 generations of outbreeding, which improves mapping precision. However, functional genomics resources for post-GWAS interpretation of genomic findings have been limited. Results: We processed 3,160 RNA-seq samples, 1,932 of which were previously unpublished, from 13 tissues as part of RatGTEx. We used Pantry to quantify six modalities of transcriptome regulation (expression, isoform ratio, intron excision ratio, alternative transcription start site (TSS), alternative poly(A), and RNA stability) in each sample, which allowed us to identify 144,876 independent quantitative trait loci (QTL) for 22,599 genes. We performed a transcriptome-wide association study (TWAS) with 19,796 unique genes across tissues and transcriptional modalities and 308 traits collected from 14 studies using FUSION, which identified 7,473 significant gene-trait associations. Conclusions: This work presents a significant expansion of the RatGTEx resource in terms of samples, modes of transcriptome regulation, and genetic analyses. Our findings provide molecular evidence for the involvement of specific genes in mediating GWAS signals through a range of transcriptome regulatory mechanisms. These gene-to-phenotype relationships can be easily applied to genetic analyses in humans and other species. Further experiments can use this evidence to confirm these molecular mediators and target them for the treatment of behavioral and physiological traits and disorders. To facilitate use of these results, we have updated results on RatGTEx.org and provided a new interactive TWAS web portal ( twas.ratgtex.org ).

Locomotor response to a novel environment is a heritable, addiction-relevant behavioral phenotype in rodents. We characterized its genetic basis in 9,341 male and female heterogeneous stock (HS) rats from eight studies. Locomotor activity was heritable across studies, while pairwise genetic correlations were positive but variable, indicating shared yet nonidentical genetic influences. Given nonuniform genetic architecture, we compared three complementary approaches to locus discovery that analyzed SNP effects 1) on a pooled locomotor phenotype (MEGA), 2) across studies under a fixed-effect model (METAL), or 3) on a shared genetic factor (Genomic SEM). To ensure comparability, all three locus-discovery approaches were applied to a common subset of seven studies. Together, these approaches identified 38 loci, with both convergence and method-specific discovery across analyses. Genomic SEM identified 31 loci, compared to 8 for MEGA and 13 for METAL. To interpret the loci, we integrated brain molecular QTLs, predicted coding and splice consequences, and open-field locomotor phenotypes from knockouts of orthologous mouse genes. Our results highlighted established neurobehavioral candidates spanning locomotor and reward-related biology ( Negr1 and Lmo4 ), serotonergic signaling ( Htr1b ), and GABAergic signaling ( Gabrg2 ), as well as relatively novel genes including Tcerg1l, Pcdh10, and Tmem132e . Our results illustrate how Genomic SEM can integrate genetically correlated yet heterogeneous behavioral measures and enable discovery of loci underlying their shared genetic architecture, highlighting its utility for complex-trait studies in rodents.

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.

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, h ² = .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.

Strong communication skills are necessary for a successful scientific career. Here, we describe the curriculum for a science communication course designed to accompany high school students participating in a summer biomedical research program. The course aims to improve their scientific literacy and communication skills. Students learn to read and understand scientific literature, write a paper about their summer research project, present their results, and provide feedback to peers. We found significant improvements in students’ self-confidence in reading, interpreting, and communicating scientific data after completing the course. This successful model serves as a guide for students participating in their first research experience and provides the skills for success in future science, technology, engineering, and mathematics education and careers. The curriculum presented here can be easily adapted for any research program, including undergraduate summer research experiences and graduate student laboratory rotations.

Power analyses are often used to determine the number of animals required for a genome wide association analysis (GWAS). These analyses are typically intended to estimate the sample size needed for at least one 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 dataset that consisted of 3,173 male and female adult N/NIH heterogeneous stock (HS) rats to explore the relationship between sample size and the number of significant loci discovered. Our simulations examined the number of loci identified in sub-samples of the full dataset. The sub-sampling analysis was conducted for four 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 sub-sampled 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.

Background: Quantitative Trait Locus (QTL) analysis and Genome-Wide Association Studies (GWAS) have the power to identify variants that capture significant levels of phenotypic variance in complex traits. However, effort and time are required to select the best methods and optimize parameters and pre-processing steps. Although machine learning approaches have been shown to greatly assist in optimization and data processing, applying them to QTL analysis and GWAS is challenging due to the complexity of large, heterogenous datasets. Here, we describe proof-of-concept for an automated machine learning approach, AutoQTL, with the ability to automate many complex decisions related to analysis of complex traits and generate diverse solutions to describe relationships that exist in genetic data. Results: Using a dataset of 18 putative QTL from a large-scale GWAS of body mass index in the laboratory rat, Rattus norvegicus , AutoQTL captures the phenotypic variance explained under a standard additive model while also providing evidence of non-additive effects including deviations from additivity and 2-way epistatic interactions from simulated data via multiple optimal solutions. Additionally, feature importance metrics provide different insights into the inheritance models and predictive power of multiple GWAS-derived putative QTL. Conclusions: This proof-of-concept illustrates that automated machine learning techniques can be applied to genetic data and has the potential to detect both additive and non-additive effects via various optimal solutions and feature importance metrics. In the future, we aim to expand AutoQTL to accommodate omics-level datasets with intelligent feature selection strategies.

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 genome-wide association studies (GWASs) 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. Contribution to the field statement: Glaucoma is the leading cause of irreversible blindness worldwide. Intraocular pressure (IOP) is the only known modifiable risk factor. This study describes results of the genome-wide association study (GWAS) performed in outbred rats that identifies known and novel genes involved in IOP regulation. To our knowledge, this is the first GWAS performed for IOP in a rat model. Identifying novel candidate genes in the rat model provides insight into the risk factors for glaucoma in humans and potential pharmacological targets for regulating IOP. The rat model is advantageous for studying natural variations in IOP, controlling environmental exposures, and providing easier access to tissue that can be used in phenotyping and gene expression in future studies.

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. Author Summary: Understanding how genetic variation affects phenotypic variation is critical to leveraging the wealth of genetic studies to make biologically and therapeutically useful discoveries. Since most of the genetic loci associated with complex diseases are regulatory in nature—meaning that they do not alter protein coding but rather subtly affect gene expression—transcriptome-wide association studies have been developed. However, these apply only to human data where large samples of unrelated individuals are available. For animal models, relatedness is much higher, causing higher false-positive rates. We propose a computationally efficient method to address this problem and find shared biology between humans and rats. Taken together, our development paves the way to further explore shared biological mechanisms of complex traits across species.

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