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

The influence of the NOD Nss1/Idd5 loci on sialadenitis and gene expression in salivary glands of congenic mice

The nonobese diabetic (NOD) Nss1 and Idd5 loci have been associated with sialadenitis development in mice. In this study the NOD Nss1 and Idd5 loci were backcrossed onto the healthy control strain B10.Q by using the speed congenic breeding strategy, resulting in three congenic strains: B10.Q.Nss1, B10.Q.Nss1/Idd5 heterozygous and B10.Q.Nss1/Idd5 homozygous. We investigated the effects of the Nss1 and Idd5 loci on sialadenitis and gene expression in NOD congenic mice. One submandibular salivary gland from each mouse was used for histological analysis of sialadenitis, whereas the contralateral salivary gland was used for gene expression profiling with the Applied Biosystems Mouse Genome Survey chip v.1.0. The results were validated using quantitative reverse transcriptase PCR. The NOD Nss1 and Idd5 loci had clear influence on the onset and progression of sialadenitis in congenic mice. Double congenic mice exhibited the most severe phenotype. We successfully identified several genes that are located in the NOD congenic regions to be differentially expressed between the congenic strains and the control strain. Several of these were found to be co-regulated, such as Stat1, complement component C1q genes and Tlr12. Also, a vast contingency of interferon-regulated genes (such as Ltb, Irf7 and Irf8) and cytokine and chemokine genes (such as Ccr7 and Ccl19) were differentially expressed between the congenic strains and the control strain. Over-representation of inflammatory signalling pathways was observed among the differentially expressed genes. We have found that the introgression of the NOD loci Nss1 and Idd5 on a healthy background caused sialadenitis in NOD congenic mouse strains, and we propose that genes within these loci are important factors in the pathogenesis. Furthermore, gene expression profiling has revealed several differentially expressed genes within and outside the NOD loci that are similar to genes found to be differentially expressed in patients with Sjögren's syndrome, and as such are interesting candidates for investigation to enhance our understanding of disease mechanisms and to develop future therapies

Effects of environmental and physiological covariates on sex differences in unconditioned and conditioned anxiety and fear in a large sample of genetically heterogeneous (N/Nih-HS) rats

Physiological and environmental variables, or covariates, can account for an important portion of the variability observed in behavioural/physiological results from different laboratories even when using the same type of animals and phenotyping procedures. We present the results of a behavioural study with a sample of 1456 genetically heterogeneous N/Nih-HS rats, including males and females, which are part of a larger genome-wide fine-mapping QTL (Quantitative Trait Loci) study. N/Nih-HS rats have been derived from 8 inbred strains and provide very small distance between genetic recombinations, which makes them a unique tool for fine-mapping QTL studies. The behavioural test battery comprised the elevated zero-maze test for anxiety, novel-cage (open-field like) activity, two-way active avoidance acquisition (related to conditioned anxiety) and context-conditioned freezing (i.e. classically conditioned fear). Using factorial analyses of variance (ANOVAs) we aimed to analyse sex differences in anxiety and fear in this N/Nih-HS rat sample, as well as to assess the effects of (and interactions with) other independent factors, such as batch, season, coat colour and experimenter. Body weight was taken as a quantitative covariate and analysed by covariance analysis (ANCOVA). Obliquely-rotated factor analyses were also performed separately for each sex, in order to evaluate associations among the most relevant variables from each behavioural test and the common dimensions (i.e. factors) underlying the different behavioural responses. ANOVA analyses showed a consistent pattern of sex effects, with females showing less signs of anxiety and fear than males across all tests. There were also significant main effects of batch, season, colour and experimenter on almost all behavioural variables, as well as "sex × batch", "sex × season" and "sex × experimenter" interactions. Body weight showed significant effects in the ANCOVAs of most behavioural measures, but sex effects were still present in spite of (and after controlling for) these "body weight" effects. Factor analyses of relevant variables from each test showed a two-fold factor structure in both sexes, with the first factor mainly representing anxiety and conditioned fear in males, while in females the first factor was dominated by loadings of activity measures. Thus, besides showing consistent sex differences in anxiety-, fear- and activity-related responses in N/Nih-HS rats, the present study shows that females' behaviour is predominantly influenced by activity while males are more influenced by anxiety. Moreover, the results point out that, besides "sex" effects, physiological variables such as colour and body weight, and environmental factors as batch/season or "experimenter", have to be taken into account in both behavioural and quantitative genetic studies because of their demonstrated influences on phenotypic outcomes

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

Ansiedad en ratas genéticamente heterogéneas : hacia la identificación de genes para caracteres conductuales cuantitativos

El uso de roedores genéticamente heterogéneos constituye una estrategia única para la identificación y el 'mapeo fino' a alta resolución de locus genéticos (QTL) con influencia sobre fenotipos cuantitativos biológicos y conductuales, permitiendo la identificación de genes individuales (genes cuantitativos) con acción sobre aquéllos. Presentamos el primer estudio de este tipo con ratas genéticamente heterogéneas (N/Nih-HS; derivadas de ocho cepas de ratas consanguíneas), que se evalúan en varias pruebas conductuales que miden ansiedad/miedo no aprendidos ('Caja blanca/negra', 'Laberinto en 'cero' elevado') o aprendidos (conducta de petrificación en contexto condicionado, adquisición de la evitación activa en dos sentidos en la caja de vaivén). Las ratas N/Nih-HS presentan una conducta más parecida a la de la cepa consanguínea RLA-I (ansiosas) que a la de la RHA-I (poco ansiosas). Se hallan correlaciones significativas entre variables incondicionadas de ansiedad y de la adquisición de la evitación activa en dos sentidos, confirmadas parcialmente por análisis de regresión múltiple. Las ratas N/Nih-HS relativamente 'evitadoras' exhiben niveles más bajos de ansiedad no aprendida que las poco 'evitadoras'. Los resultados de esta evaluación del comportamiento de las ratas N/Nih-HS son discutidos en términos de su potencial utilidad para la investigación neurogenética de la ansiedad y el miedoThe use of genetically heterogeneous (outbred) rodents is a unique resource for the identification and fine mapping of genetic loci (QTL) influencing biological and behavioural quantitative phenotypes, allowing the identification of quantitative genes. We present the first study of this kind carried out with genetically heterogeneous rats (N/Nih-HS; derivated from an eight-way cross of inbred strains), whose behaviour is assessed in tests evoking unlearned (Black/white box, Elevated 'zero' maze) or learned (context conditioned freezing, two-way active avoidance acquisition in the shuttlebox) anxious/fearful responses. The behavioural profile of N/Nih-HS rats is more similar to that of RLA-I (anxious) rats rather than to RHA-I (low anxious) rats. Significant correlations are found among unconditioned anxiety variables and two-way active avoidance acquisition in the shuttlebox; these are partially confirmed by multiple regression analysis. 'High avoider' N/Nih-HS rats show lower unlearned anxiety levels than 'low avoiders'. Results of this behavioural assessment of the N/Nih-HS rats are discussed in terms of their potential usefulness for present and future neurobehavioural and genetic studies of fearfulness and anxiet

NANOG Reporter Cell Lines Generated by Gene Targeting in Human Embryonic Stem Cells

Background: Pluripotency and self-renewal of human embryonic stem cells (hESCs) is mediated by a complex interplay between extra- and intracellular signaling pathways, which regulate the expression of pluripotency-specific transcription factors. The homeodomain transcription factor NANOG plays a central role in maintaining hESC pluripotency, but the precise role and regulation of NANOG are not well defined. Methodology/Principal Findings: To facilitate the study of NANOG expression and regulation in viable hESC cultures, we generated fluorescent NANOG reporter cell lines by gene targeting in hESCs. In these reporter lines, the fluorescent reporter gene was co-expressed with endogenous NANOG and responded to experimental induction or repression of the NANOG promoter with appropriate changes in expression levels. Furthermore, NANOG reporter lines facilitated the separation of hESC populations based on NANOG expression levels and their subsequent characterization. Gene expression arrays on isolated hESC subpopulations revealed genes with differential expression in NANOG high and NANOG low hESCs, providing candidates for NANOG downstream targets hESCs. Conclusion/Significance: The newly derived NANOG reporter hESC lines present novel tools to visualize NANOG expression in viable hESCs. In future applications, these reporter lines can be used to elucidate the function and regulation of NANO

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

FGF4 and Retinoic Acid Direct Differentiation of hESCs into PDX1-Expressing Foregut Endoderm in a Time- and Concentration-Dependent Manner

BACKGROUND: Retinoic acid (RA) and fibroblast growth factor 4 (FGF4) signaling control endoderm patterning and pancreas induction/expansion. Based on these findings, RA and FGFs, excluding FGF4, have frequently been used in differentiation protocols to direct differentiation of hESCs into endodermal and pancreatic cell types. In vivo, these signaling pathways act in a temporal and concentration-dependent manner. However, in vitro, the underlying basis for the time of addition of growth and differentiation factors (GDFs), including RA and FGFs, as well as the concentration is lacking. Thus, in order to develop robust and reliable differentiation protocols of ESCs into mature pancreatic cell types, including insulin-producing beta cells, it will be important to mechanistically understand each specification step. This includes differentiation of mesendoderm/definitive endoderm into foregut endoderm--the origin of pancreatic endoderm. METHODOLOGY/PRINCIPAL FINDINGS: Here, we provide data on the individual and combinatorial role of RA and FGF4 in directing differentiation of ActivinA (AA)-induced hESCs into PDX1-expressing cells. FGF4's ability to affect endoderm patterning and specification in vitro has so far not been tested. By testing out the optimal concentration and timing of addition of FGF4 and RA, we present a robust differentiation protocol that on average generates 32% PDX1(+) cells. Furthermore, we show that RA is required for converting AA-induced hESCs into PDX1(+) cells, and that part of the underlying mechanism involves FGF receptor signaling. Finally, further characterization of the PDX1(+) cells suggests that they represent foregut endoderm not yet committed to pancreatic, posterior stomach, or duodenal endoderm. CONCLUSION/SIGNIFICANCE: In conclusion, we show that RA and FGF4 jointly direct differentiation of PDX1(+) foregut endoderm in a robust and efficient manner. RA signaling mediated by the early induction of RARbeta through AA/Wnt3a is required for PDX1 expression. Part of RA's activity is mediated by FGF signaling

Genome-wide analysis identifies 12 loci influencing human reproductive behavior.

The genetic architecture of human reproductive behavior-age at first birth (AFB) and number of children ever born (NEB)-has a strong relationship with fitness, human development, infertility and risk of neuropsychiatric disorders. However, very few genetic loci have been identified, and the underlying mechanisms of AFB and NEB are poorly understood. We report a large genome-wide association study of both sexes including 251,151 individuals for AFB and 343,072 individuals for NEB. We identified 12 independent loci that are significantly associated with AFB and/or NEB in a SNP-based genome-wide association study and 4 additional loci associated in a gene-based effort. These loci harbor genes that are likely to have a role, either directly or by affecting non-local gene expression, in human reproduction and infertility, thereby increasing understanding of these complex traits