Genomic Loci Influencing Cue-Reactivity in Heterogeneous Stock Rats

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

Evaluation of an 8-week high school science communication course designed to read, write, and present scientific research

The development of science writing and presentation skills is necessary for a successful science career. Too often these skills are not included in pre- or postsecondary science, technology, engineering, and mathematics (STEM) education, leading to a disconnect between high schoolers' expectations for college preparedness and the skills needed to succeed in college. The Young Scientist Program Summer Focus recruits high school students from historically marginalized backgrounds to participate in 8-week summer internships at Washington University in St. Louis. Students conduct hands-on biomedical research projects under the mentorship of Washington University scientists (graduate students, postdoctorates, lab staff). Here, we present the curriculum for a science communication course that accompanies this early research experience. The course is designed to strengthen students' communication skills (critical reading, writing, presenting, and peer review) through a combination of weekly lectures and active learning methods. It prepares students for the capstone of their summer internship: writing a scientific paper and presenting their results at a closing symposium. We administered pre- and postprogram surveys to four Summer Focus cohorts to determine whether the course met its learning objectives. We found significant improvements in students' self-confidence in reading, interpreting, and communicating scientific data. Thus, this course provides a successful model for introducing science literacy and communication skills that are necessary for any career in STEM. We provide a detailed outline of the course structure and content so that this training can be incorporated into any undergraduate and graduate research programs.NEW & NOTEWORTHY 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.

Genome-wide association studies of human and rat BMI converge on synapse, epigenome, and hormone signaling networks

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.

Automated quantitative trait locus analysis (AutoQTL)

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 complicated decisions related to analysis of complex traits and generate solutions to describe relationships that exist in genetic data.

RESULTS

Using a publicly available 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. AutoQTL also detects evidence of non-additive effects including deviations from additivity and 2-way epistatic interactions in 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 complement standard approaches and have 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 and feature engineering strategies.

Genome-wide association study finds multiple loci associated with intraocular pressure in HS rats

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

Genetic, epigenetic, and environmental mechanisms govern allele-specific gene expression

Allele-specific expression (ASE) is a phenomenon in which one allele is preferentially expressed over the other. Genetic and epigenetic factors cause ASE by altering the final composition of a gene's product, leading to expression imbalances that can have functional consequences on phenotypes. Environmental signals also impact allele-specific expression, but how they contribute to this cross talk remains understudied. Here, we explored how genotype, parent-of-origin, tissue, sex, and dietary fat simultaneously influence ASE biases. Male and female mice from a F₁ reciprocal cross of the LG/J and SM/J strains were fed a high or low fat diet. We harnessed strain-specific variants to distinguish between two ASE classes: parent-of-origin-dependent (unequal expression based on parental origin) and sequence-dependent (unequal expression based on nucleotide identity). We present a comprehensive map of ASE patterns in 2853 genes across three tissues and nine environmental contexts. We found that both ASE classes are highly dependent on tissue and environmental context. They vary across metabolically relevant tissues, between males and females, and in response to dietary fat. We also found 45 genes with inconsistent ASE biases that switched direction across tissues and/or environments. Finally, we integrated ASE and QTL data from published intercrosses of the LG/J and SM/J strains. Our ASE genes are often enriched in QTLs for metabolic and musculoskeletal traits, highlighting how this orthogonal approach can prioritize candidate genes. Together, our results provide novel insights into how genetic, epigenetic, and environmental mechanisms govern allele-specific expression, which is an essential step toward deciphering the genotype-to-phenotype map.

Parent-of-origin effects propagate through networks to shape metabolic traits

Parent-of-origin effects are unexpectedly common in complex traits, including metabolic and neurological traits. Parent-of-origin effects can be modified by the environment, but the architecture of these gene-by-environmental effects on phenotypes remains to be unraveled. Previously, quantitative trait loci (QTL) showing context-specific parent-of-origin effects on metabolic traits were mapped in the F16 generation of an advanced intercross between LG/J and SM/J inbred mice. However, these QTL were not enriched for known imprinted genes, suggesting another mechanism is needed to explain these parent-of-origin effects phenomena. We propose that non-imprinted genes can generate complex parent-of-origin effects on metabolic traits through interactions with imprinted genes. Here, we employ data from mouse populations at different levels of intercrossing (F0, F1, F2, F16) of the LG/J and SM/J inbred mouse lines to test this hypothesis. Using multiple populations and incorporating genetic, genomic, and physiological data, we leverage orthogonal evidence to identify networks of genes through which parent-of-origin effects propagate. We identify a network comprised of three imprinted and six non-imprinted genes that show parent-of-origin effects. This epistatic network forms a nutritional responsive pathway and the genes comprising it jointly serve cellular functions associated with growth. We focus on two genes, Nnat and F2r, whose interaction associates with serum glucose levels across generations in high-fat-fed females. Single-cell RNAseq reveals that Nnat expression increases and F2r expression decreases in pre-adipocytes along an adipogenic trajectory, a result that is consistent with our observations in bulk white adipose tissue.

SNPs, short tandem repeats, and structural variants are responsible for differential gene expression across C57BL/6 and C57BL/10 substrains

Mouse substrains are an invaluable model for understanding disease. We compared C57BL/6J, which is the most commonly used inbred mouse strain, with eight C57BL/6 and five C57BL/10 closely related inbred substrains. Whole-genome sequencing and RNA-sequencing analysis yielded 352,631 SNPs, 109,096 indels, 150,344 short tandem repeats (STRs), 3,425 structural variants (SVs), and 2,826 differentially expressed genes (DE genes) among these 14 strains; 312,981 SNPs (89%) distinguished the B6 and B10 lineages. These SNPs were clustered into 28 short segments that are likely due to introgressed haplotypes rather than new mutations. Outside of these introgressed regions, we identified 53 SVs, protein-truncating SNPs, and frameshifting indels that were associated with DE genes. Our results can be used for both forward and reverse genetic approaches and illustrate how introgression and mutational processes give rise to differences among these widely used inbred substrains.

Content and Performance of the MiniMUGA Genotyping Array: A New Tool To Improve Rigor and Reproducibility in Mouse Research

The laboratory mouse is the most widely used animal model for biomedical research, due in part to its well-annotated genome, wealth of genetic resources, and the ability to precisely manipulate its genome. Despite the importance of genetics for mouse research, genetic quality control (QC) is not standardized, in part due to the lack of cost-effective, informative, and robust platforms. Genotyping arrays are standard tools for mouse research and remain an attractive alternative even in the era of high-throughput whole-genome sequencing. Here, we describe the content and performance of a new iteration of the Mouse Universal Genotyping Array (MUGA), MiniMUGA, an array-based genetic QC platform with over 11,000 probes. In addition to robust discrimination between most classical and wild-derived laboratory strains, MiniMUGA was designed to contain features not available in other platforms: (1) chromosomal sex determination, (2) discrimination between substrains from multiple commercial vendors, (3) diagnostic SNPs for popular laboratory strains, (4) detection of constructs used in genetically engineered mice, and (5) an easy-to-interpret report summarizing these results. In-depth annotation of all probes should facilitate custom analyses by individual researchers. To determine the performance of MiniMUGA, we genotyped 6899 samples from a wide variety of genetic backgrounds. The performance of MiniMUGA compares favorably with three previous iterations of the MUGA family of arrays, both in discrimination capabilities and robustness. We have generated publicly available consensus genotypes for 241 inbred strains including classical, wild-derived, and recombinant inbred lines. Here, we also report the detection of a substantial number of XO and XXY individuals across a variety of sample types, new markers that expand the utility of reduced complexity crosses to genetic backgrounds other than C57BL/6, and the robust detection of 17 genetic constructs. We provide preliminary evidence that the array can be used to identify both partial sex chromosome duplication and mosaicism, and that diagnostic SNPs can be used to determine how long inbred mice have been bred independently from the relevant main stock. We conclude that MiniMUGA is a valuable platform for genetic QC, and an important new tool to increase the rigor and reproducibility of mouse research.

Brown Adipose Expansion and Remission of Glycemic Dysfunction in Obese SM/J Mice

We leverage the SM/J mouse to understand glycemic control in obesity. High-fat-fed SM/J mice initially develop poor glucose homeostasis relative to controls. Strikingly, their glycemic dysfunction resolves by 30 weeks of age despite persistent obesity. The mice dramatically expand their brown adipose depots as they resolve glycemic dysfunction. This occurs naturally and spontaneously on a high-fat diet, with no temperature or genetic manipulation. Removal of the brown adipose depot impairs insulin sensitivity, indicating that the expanded tissue is functioning as an insulin-stimulated glucose sink. We describe morphological, physiological, and transcriptomic changes that occur during the brown adipose expansion and remission of glycemic dysfunction, and focus on Sfrp1 (secreted frizzled-related protein 1) as a compelling candidate that may underlie this phenomenon. Understanding how the expanded brown adipose contributes to glycemic control in SM/J mice will open the door for innovative therapies aimed at improving metabolic complications in obesity.

Spontaneous restoration of functional β-cell mass in obese SM/J mice

Maintenance of functional β-cell mass is critical to preventing diabetes, but the physiological mechanisms that cause β-cell populations to thrive or fail in the context of obesity are unknown. High fat-fed SM/J mice spontaneously transition from hyperglycemic-obese to normoglycemic-obese with age, providing a unique opportunity to study β-cell adaptation. Here, we characterize insulin homeostasis, islet morphology, and β-cell function during SM/J's diabetic remission. As they resolve hyperglycemia, obese SM/J mice dramatically increase circulating and pancreatic insulin levels while improving insulin sensitivity. Immunostaining of pancreatic sections reveals that obese SM/J mice selectively increase β-cell mass but not α-cell mass. Obese SM/J mice do not show elevated β-cell mitotic index, but rather elevated α-cell mitotic index. Functional assessment of isolated islets reveals that obese SM/J mice increase glucose-stimulated insulin secretion, decrease basal insulin secretion, and increase islet insulin content. These results establish that β-cell mass expansion and improved β-cell function underlie the resolution of hyperglycemia, indicating that obese SM/J mice are a valuable tool for exploring how functional β-cell mass can be recovered in the context of obesity.

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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

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

Genome-Wide Association Study in Two Cohorts from a Multi-generational Mouse Advanced Intercross Line Highlights the Difficulty of Replication Due to Study-Specific Heterogeneity

There has been extensive discussion of the "Replication Crisis" in many fields, including genome-wide association studies (GWAS). We explored replication in a mouse model using an advanced intercross line (AIL), which is a multigenerational intercross between two inbred strains. We re-genotyped a previously published cohort of LG/J x SM/J AIL mice (F34; n = 428) using a denser marker set and genotyped a new cohort of AIL mice (F39-43; n = 600) for the first time. We identified 36 novel genome-wide significant loci in the F34 and 25 novel loci in the F39-43 cohort. The subset of traits that were measured in both cohorts (locomotor activity, body weight, and coat color) showed high genetic correlations, although the SNP heritabilities were slightly lower in the F39-43 cohort. For this subset of traits, we attempted to replicate loci identified in either F34 or F39-43 in the other cohort. Coat color was robustly replicated; locomotor activity and body weight were only partially replicated, which was inconsistent with our power simulations. We used a random effects model to show that the partial replications could not be explained by Winner's Curse but could be explained by study-specific heterogeneity. Despite this heterogeneity, we performed a mega-analysis by combining F34 and F39-43 cohorts (n = 1,028), which identified four novel loci associated with locomotor activity and body weight. These results illustrate that even with the high degree of genetic and environmental control possible in our experimental system, replication was hindered by study-specific heterogeneity, which has broad implications for ongoing concerns about reproducibility.

The Genetic Relationship Between Alcohol Consumption and Aspects of Problem Drinking in an Ascertained Sample

BACKGROUND

Genomewide association studies (GWAS) have begun to identify loci related to alcohol consumption, but little is known about whether this genetic propensity overlaps with specific indices of problem drinking in ascertained samples.

METHODS

In 6,731 European Americans who had been exposed to alcohol, we examined whether polygenic risk scores (PRS) from a GWAS of weekly alcohol consumption in the UK Biobank predicted variance in 6 alcohol-related phenotypes: alcohol use, maximum drinks within 24 hours (MAXD), total score on the Self-Rating of the Effects of Ethanol Questionnaire (SRE-T), DSM-IV alcohol dependence (DSM4AD), DSM-5 alcohol use disorder symptom counts (DSM5AUDSX), and reduction/cessation of problematic drinking. We also examined the extent to which an single nucleotide polymorphism (rs1229984) in ADH1B, which is strongly associated with both alcohol consumption and dependence, contributed to the polygenic association with these phenotypes and whether PRS interacted with sex, age, or family history of alcoholism to predict alcohol-related outcomes. We performed mixed-effect regression analyses, with family membership and recruitment site included as random effects, as well as survival modeling of age of onset of DSM4AD.

RESULTS

PRS for alcohol consumption significantly predicted variance in 5 of the 6 outcomes: alcohol use (Δmarginal R2  = 1.39%, Δ area under the curve [AUC] = 0.011), DSM4AD (Δmarginal R2  = 0.56%; ΔAUC = 0.003), DSM5AUDSX (Δmarginal R2  = 0.49%), MAXD (Δmarginal R2  = 0.31%), and SRE-T (Δmarginal R2  = 0.22%). PRS were also associated with onset of DSM4AD (hazard ratio = 1.11, p = 2.08e-5). The inclusion of rs1229984 attenuated the effects of the alcohol consumption PRS, particularly for DSM4AD and DSM5AUDSX, but the PRS continued to exert an independent effect for all 5 alcohol measures (Δmarginal R2 after controlling for ADH1B = 0.14 to 1.22%). Interactions between PRS and sex, age, or family history were nonsignificant.

CONCLUSIONS

Genetic propensity for typical alcohol consumption was associated with alcohol use and was also associated with 4 of the additional 5 outcomes, though the variance explained in this sample was modest. Future GWAS that focus on the multifaceted nature of AUD, which goes beyond consumption, might reveal additional information regarding the polygenic underpinnings of problem drinking.

Genome wide association analysis in a mouse advanced intercross line

The LG/J x SM/J advanced intercross line of mice (LG x SM AIL) is a multigenerational outbred population. High minor allele frequencies, a simple genetic background, and the fully sequenced LG and SM genomes make it a powerful population for genome-wide association studies. Here we use 1,063 AIL mice to identify 126 significant associations for 50 traits relevant to human health and disease. We also identify thousands of cis- and trans-eQTLs in the hippocampus, striatum, and prefrontal cortex of ~200 mice. We replicate an association between locomotor activity and Csmd1, which we identified in an earlier generation of this AIL, and show that Csmd1 mutant mice recapitulate the locomotor phenotype. Our results demonstrate the utility of the LG x SM AIL as a mapping population, identify numerous novel associations, and shed light on the genetic architecture of mammalian behavior.