Genomes and phenomes of a population of outbred rats and its progenitors

Finding genetic variants that contribute to phenotypic variation is one of the main challenges of modern genetics. We used an outbred population of rats (Heterogeneous Stock, HS) in a combined sequence-based and genetic mapping analysis to identify sequence variants and genes contributing to complex traits of biomedical relevance. Here we describe the sequences of the eight inbred progenitors of the HS and the variants that segregate between them. We report the genotyping of 1,407 HS rats, and the collection from 2,006 rats of 195 phenotypic measures that are relevant to models of anxiety, type 2 diabetes, hypertension and osteoporosis. We make available haplotype dosages for the 1,407 genotyped rats, since genetic mapping in the HS is best carried out by reconstructing each HS chromosome as a mosaic of the progenitor genomes. Finally, we have deposited an R object that makes it easy to incorporate our sequence data into any genetic study of HS rats. Our genetic data are available for both Rnor3.4 and Rnor5.0 rat assemblies

Dominance is common in mammals and is associated with trans-acting gene expression and alternative splicing

Background: Dominance and other non-additive genetic effects arise from the interaction between alleles, and historically these phenomena play a major role in quantitative genetics. However, most genome-wide association studies (GWAS) assume alleles act additively. // Results: We systematically investigate both dominance—here representing any non-additive within-locus interaction—and additivity across 574 physiological and gene expression traits in three mammalian stocks: F2 intercross pigs, rat heterogeneous stock, and mice heterogeneous stock. Dominance accounts for about one quarter of heritable variance across all physiological traits in all species. Hematological and immunological traits exhibit the highest dominance variance, possibly reflecting balancing selection in response to pathogens. Although most quantitative trait loci (QTLs) are detectable as additive QTLs, we identify 154, 64, and 62 novel dominance QTLs in pigs, rats, and mice respectively that are undetectable as additive QTLs. Similarly, even though most cis-acting expression QTLs are additive, gene expression exhibits a large fraction of dominance variance, and trans-acting eQTLs are enriched for dominance. Genes causal for dominance physiological QTLs are less likely to be physically linked to their QTLs but instead act via trans-acting dominance eQTLs. In addition, thousands of eQTLs are associated with alternatively spliced isoforms with complex additive and dominant architectures in heterogeneous stock rats, suggesting a possible mechanism for dominance. // Conclusions: Although heritability is predominantly additive, many mammalian genetic effects are dominant and likely arise through distinct mechanisms. It is therefore advantageous to consider both additive and dominance effects in GWAS to improve power and uncover causality

High-resolution genome screen for bone mineral density in heterogeneous stock rat

We previously demonstrated that skeletal mass, structure, and biomechanical properties vary considerably in heterogeneous stock (HS) rat strains. In addition, we observed strong heritability for several of these skeletal phenotypes in the HS rat model, suggesting that it represents a unique genetic resource for dissecting the complex genetics underlying bone fragility. The purpose of this study was to identify and localize genes associated with bone mineral density in HS rats. We measured bone phenotypes from 1524 adult male and female HS rats between 17 and 20 weeks of age. Phenotypes included dual-energy X-ray absorptiometry (DXA) measurements for bone mineral content and areal bone mineral density (aBMD) for femur and lumbar spine (L3-L5), and volumetric BMD measurements by CT for the midshaft and distal femur, femur neck, and fifth lumbar vertebra (L5). A total of 70,000 polymorphic single-nucleotide polymorphisms (SNPs) distributed throughout the genome were selected from genotypes obtained from the Affymetrix rat custom SNPs array for the HS rat population. These SNPs spanned the HS rat genome with a mean linkage disequilibrium coefficient between neighboring SNPs of 0.95. Haplotypes were estimated across the entire genome for each rat using a multipoint haplotype reconstruction method, which calculates the probability of descent for each genotyped locus from each of the eight founder HS strains. The haplotypes were tested for association with each bone density phenotype via a mixed model with covariate adjustment. We identified quantitative trait loci (QTLs) for BMD phenotypes on chromosomes 2, 9, 10, and 13 meeting a conservative genomewide empiric significance threshold (false discovery rate [FDR] = 5%; p < 3 × 10(-6)). Importantly, most QTLs were localized to very small genomic regions (1-3 megabases [Mb]), allowing us to identify a narrow set of potential candidate genes including both novel genes and genes previously shown to have roles in skeletal development and homeostasis

Fine mapping of bone structure and strength QTLs in heterogeneous stock rat

We previously demonstrated that skeletal structure and strength phenotypes vary considerably in heterogeneous stock (HS) rats. These phenotypes were found to be strongly heritable, suggesting that the HS rat model represents a unique genetic resource for dissecting the complex genetic etiology underlying bone fragility. The purpose of this study was to identify and localize genes associated with bone structure and strength phenotypes using 1524 adult male and female HS rats between 17 to 20 weeks of age. Structure measures included femur length, neck width, head width; femur and lumbar spine (L3-5) areas obtained by DXA; and cross-sectional areas (CSA) at the midshaft, distal femur and femoral neck, and the 5th lumbar vertebra measured by CT. In addition, measures of strength of the whole femur and femoral neck were obtained. Approximately 70,000 polymorphic SNPs distributed throughout the rat genome were selected for genotyping, with a mean linkage disequilibrium coefficient between neighboring SNPs of 0.95. Haplotypes were estimated across the entire genome for each rat using a multipoint haplotype reconstruction method, which calculates the probability of descent at each locus from each of the 8 HS founder strains. The haplotypes were then tested for association with each structure and strength phenotype via a mixed model with covariate adjustment. We identified quantitative trait loci (QTLs) for structure phenotypes on chromosomes 3, 8, 10, 12, 17 and 20, and QTLs for strength phenotypes on chromosomes 5, 10 and 11 that met a conservative genome-wide empiric significance threshold (FDR=5%; P<3×10(-6)). Importantly, most QTLs were localized to very narrow genomic regions (as small as 0.3 Mb and up to 3 Mb), each harboring a small set of candidate genes, both novel and previously shown to have roles in skeletal development and homeostasis

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

BackgroundMultiple factors contribute to the etiology of addiction, including genetics, sex, and a number of addiction-related behavioral traits. One behavioral trait where individuals assign incentive salience to food stimuli ("sign-trackers", ST) are more impulsive compared to those that do not ("goal-trackers", GT), as well as more sensitive to drugs and drug stimuli. Furthermore, this GT/ST phenotype predicts differences in other behavioral measures. Recent studies have implicated the gut microbiota as a key regulator of brain and behavior, and have shown that many microbiota-associated changes occur in a sex-dependent manner. However, few studies have examined how the microbiome might influence addiction-related behaviors. To this end, we sought to determine if gut microbiome composition was correlated with addiction-related behaviors determined by the GT/ST phenotype.MethodsOutbred male (N=101) and female (N=101) heterogeneous stock rats underwent a series of behavioral tests measuring impulsivity, attention, reward-learning, incentive salience, and locomotor response. Cecal microbiome composition was estimated using 16S rRNA gene amplicon sequencing. Behavior and microbiome were characterized and correlated with behavioral phenotypes. Robust sex differences were observed in both behavior and microbiome; further analyses were conducted within sex using the pre-established goal/sign-tracking (GT/ST) phenotype and partial least squares differential analysis (PLS-DA) clustered behavioral phenotype.ResultsOverall microbiome composition was not associated to the GT/ST phenotype. However, microbial alpha diversity was significantly decreased in female STs. On the other hand, a measure of impulsivity had many significant correlations to microbiome in both males and females. Several measures of impulsivity were correlated with the genus Barnesiella in females. Female STs had notable correlations between microbiome and attentional deficient. In both males and females, many measures were correlated with the bacterial families Ruminocococcaceae and Lachnospiraceae.ConclusionsThese data demonstrate correlations between several addiction-related behaviors and the microbiome specific to sex

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

Background: Multiple factors contribute to the etiology of addiction, including genetics, sex, and a number of addiction-related behavioral traits. One behavioral trait where individuals assign incentive salience to food stimuli (â sign-trackersâ , ST) are more impulsive compared to those that do not (â goal-trackersâ , GT), as well as more sensitive to drugs and drug stimuli. Furthermore, this GT/ST phenotype predicts differences in other behavioral measures. Recent studies have implicated the gut microbiota as a key regulator of brain and behavior, and have shown that many microbiota-associated changes occur in a sex-dependent manner. However, few studies have examined how the microbiome might influence addiction-related behaviors. To this end, we sought to determine if gut microbiome composition was correlated with addiction-related behaviors determined by the GT/ST phenotype. Methods: Outbred male (N=101) and female (N=101) heterogeneous stock rats underwent a series of behavioral tests measuring impulsivity, attention, reward-learning, incentive salience, and locomotor response. Cecal microbiome composition was estimated using 16S rRNA gene amplicon sequencing. Behavior and microbiome were characterized and correlated with behavioral phenotypes. Robust sex differences were observed in both behavior and microbiome; further analyses were conducted within sex using the pre-established goal/sign-tracking (GT/ST) phenotype and partial least squares differential analysis (PLS-DA) clustered behavioral phenotype. Results: Overall microbiome composition was not associated to the GT/ST phenotype. However, microbial alpha diversity was significantly decreased in female STs. On the other hand, a measure of impulsivity had many significant correlations to microbiome in both males and females. Several measures of impulsivity were correlated with the genus Barnesiella in females. Female STs had notable correlations between microbiome and attentional deficient. In both males and females, many measures were correlated with the bacterial families Ruminocococcaceae and Lachnospiraceae. Conclusions: These data demonstrate correlations between several addiction-related behaviors and the microbiome specific to sex

Natural Polymorphisms in Tap2 Influence Negative Selection and CD4 : CD8 Lineage Commitment in the Rat

Contains fulltext : 136368.pdf (publisher's version ) (Open Access)Genetic variation in the major histocompatibility complex (MHC) affects CD4ratioCD8 lineage commitment and MHC expression. However, the contribution of specific genes in this gene-dense region has not yet been resolved. Nor has it been established whether the same genes regulate MHC expression and T cell selection. Here, we assessed the impact of natural genetic variation on MHC expression and CD4ratioCD8 lineage commitment using two genetic models in the rat. First, we mapped Quantitative Trait Loci (QTLs) associated with variation in MHC class I and II protein expression and the CD4ratioCD8 T cell ratio in outbred Heterogeneous Stock rats. We identified 10 QTLs across the genome and found that QTLs for the individual traits colocalized within a region spanning the MHC. To identify the genes underlying these overlapping QTLs, we generated a large panel of MHC-recombinant congenic strains, and refined the QTLs to two adjacent intervals of approximately 0.25 Mb in the MHC-I and II regions, respectively. An interaction between these intervals affected MHC class I expression as well as negative selection and lineage commitment of CD8 single-positive (SP) thymocytes. We mapped this effect to the transporter associated with antigen processing 2 (Tap2) in the MHC-II region and the classical MHC class I gene(s) (RT1-A) in the MHC-I region. This interaction was revealed by a recombination between RT1-A and Tap2, which occurred in 0.2% of the rats. Variants of Tap2 have previously been shown to influence the antigenicity of MHC class I molecules by altering the MHC class I ligandome. Our results show that a restricted peptide repertoire on MHC class I molecules leads to reduced negative selection of CD8SP cells. To our knowledge, this is the first study showing how a recombination between natural alleles of genes in the MHC influences lineage commitment of T cells

A multiple-phenotype imputation method for genetic studies

Genetic association studies have yielded a wealth of biologic discoveries. However, these have mostly analyzed one trait and one SNP at a time, thus failing to capture the underlying complexity of these datasets. Joint genotypephenotype analyses of complex, high-dimensional datasets represent an important way to move beyond simple GWAS with great potential. The move to high-dimensional phenotypes will raise many new statistical problems. In this paper we address the central issue of missing phenotypes in studies with any level of relatedness between samples. We propose a multiple phenotype mixed model and use a computationally efficient variational Bayesian algorithm to fit the model. On a variety of simulated and real datasets from a range of organisms and trait types, we show that our method outperforms existing state-of-the-art methods from the statistics and machine learning literature and can boost signals of associatio

Fine-mapping complex traits in heterogeneous stock rats

The fundamental theme my thesis explores is the relationship between genetic variation and phenotypic variation. It addresses three main questions. What is the genetic architecture of traits in the HS? How can sequence information help identifying the sequence variants and genes responsible for phenotypic variation? Are the genetic factors contributing to phenotypic variation in the rat homologous to those contributing to variation in the same phenotype in the mouse?To address these questions, I analysed data collected by the EURATRANS consortium on 1,407 Heterogeneous Stock (HS) rats descended from eight inbred strains through sixty generations of outbreeding. The HS rats were genotyped at 803,485 SNPs and 160 measures relevant to a number of models of disease (e.g. anxiety, type 2 diabetes, multiple sclerosis) were collected. The eight founders of the Stock were genotyped and sequenced. I identified loci in the genome that contribute to phenotypic variation (Quantitative Trait Loci, QTLs), and integrated sequence information with the mapping results to identify the genetic variants underlying the QTLs.I made some important observations about the nature of genetic architecture in rats, and how this compares to mice and humans. I also showed how sequence information can be used to improve mapping resolution, and in some cases to identify causal variants. However, I report an unexpected observation: at the majority of QTLs, the genetic effect cannot be accounted for by a single variant. This finding suggests that genetic variation cannot be reduced to sequence variation. This complexity will need to be taken into account by studies that aim at unravelling the genetic basis of complex traits.</p

Identifying genes for neurobehavioural traits in rodents: progress and pitfalls

Identifying genes and pathways that contribute to differences in neurobehavioural traits is a key goal in psychiatric research. Despite considerable success in identifying quantitative trait loci (QTLs) associated with behaviour in laboratory rodents, pinpointing the causal variants and genes is more challenging. For a long time, the main obstacle was the size of QTLs, which could encompass tens if not hundreds of genes. However, recent studies have exploited mouse and rat resources that allow mapping of phenotypes to narrow intervals, encompassing only a few genes. Here, we review these studies, showcase the rodent resources they have used and highlight the insights into neurobehavioural traits provided to date. We discuss what we see as the biggest challenge in the field – translating QTLs into biological knowledge by experimentally validating and functionally characterizing candidate genes – and propose that the CRISPR/Cas genome-editing system holds the key to overcoming this obstacle. Finally, we challenge traditional views on inbred versus outbred resources in the light of recent resource and technology developments