A high-resolution consensus linkage map of the rat, integrating radiation hybrid and genetic maps

Rat Genome Assemblies, Annotation, and Variant Repository

The first and only published version of the rat reference genome sequence was RGSC3.1, accomplished by the Rat Genome Sequencing Project Consortium. Here we present the history of the community effort in the correction of sequence errors and filling missing gaps in the process of refining and providing researchers with a high-quality rat reference sequence. The genome assembly improvements, addition of different evidence resources over time, such as RNA-Seq data, and software development methodologies had a positive impact on the gene model annotations. Over the years we observed a great increase in the numbers of genes, protein coding sequences, predicted transcripts and transcript features. Before the sequencing of the rat genome was possible, first biochemical and next genomic markers like RAPD, AFLP, RFLP, and SSLP were fundamental in research studies involving cross-breeding between different rat strains, in finding the level of polymorphism, linkage mapping, and phylogeny. Linkage maps provide information on recombination rates, give insight into intra- and interspecies gene rearrangements, and help to identify Mendelian loci and Quantitative Trait Loci (QTL). In the 1990s many reports were published on the construction of rat linkage maps that incorporated increasing numbers of markers and facilitated the localization of disease loci. Current genetic monitoring and linkage mapping relies on single nucleotide polymorphisms (SNPs). The Rat Genome Database collects information on genetic variation from the worldwide community of rat researchers and provides tools for searching and retrieving these data. As of today we show details about almost 605 million variants coming from many studies in our Variant Visualizer tool.

A High-Resolution Genetic Map for the Laboratory Rat

An accurate and high-resolution genetic map is critical for mapping complex traits, yet the resolution of the current rat genetic map is far lower than human and mouse, and has not been updated since the original Jensen-Seaman map in 2004. For the first time, we have refined the rat genetic map to sub-centimorgan (cM) resolution (<0.02 cM) by using 95,769 genetic markers and 870 informative meioses from a cohort of 528 heterogeneous stock (HS) rats. Global recombination rates in the revised sex-averaged map (0.66 cM/Mb) did not differ compared to the historical map (0.65 cM/Mb); however, substantial refinement was made to the localization of highly recombinant regions within the revised map. Also for the first time, sex-specific rat genetic maps were generated, which revealed both genomewide and fine-scale variation in recombination rates between male and female rats. Reanalysis of multiple quantitative trait loci (QTL) using the historical and refined rat genetic maps demonstrated marked changes to QTL localization, shape, and effect size. As a resource to the rat research community, we have provided revised centimorgan positions for all physical positions within the rat genome and commonly used genetic markers for trait mapping, including 44,828 SSLP markers and the RATDIV genotyping array. Collectively, this study provides a substantial improvement to the rat genetic map and an unprecedented resource for analysis of complex traits and recombination in the rat.

Variation and Evolution of the Meiotic Requirement for Crossing Over in Mammals

The segregation of homologous chromosomes at the first meiotic division is dependent on the presence of at least one well-positioned crossover per chromosome. In some mammalian species, however, the genomic distribution of crossovers is consistent with a more stringent baseline requirement of one crossover per chromosome arm. Given that the meiotic requirement for crossing over defines the minimum frequency of recombination necessary for the production of viable gametes, determining the chromosomal scale of this constraint is essential for defining crossover profiles predisposed to aneuploidy and understanding the parameters that shape patterns of recombination rate evolution across species. Here, I use cytogenetic methods for in situ imaging of crossovers in karyotypically diverse house mice (Mus musculus domesticus) and voles (genus Microtus) to test how chromosome number and configuration constrain the distribution of crossovers in a genome. I show that the global distribution of crossovers in house mice is thresholded by a minimum of one crossover per chromosome arm, whereas the crossover landscape in voles is defined by a more relaxed requirement of one crossover per chromosome. I extend these findings in an evolutionary metaanalysis of published recombination and karyotype data for 112 mammalian species and demonstrate that the physical scale of the genomic crossover distribution has undergone multiple independent shifts from one crossover per chromosome arm to one per chromosome during mammalian evolution. Together, these results indicate that the chromosomal scale constraint on crossover rates is itself a trait that evolves among species, a finding that casts light on an important source of crossover rate variation in mammals.

The Role of Prolactin in Bone Metastasis and Breast Cancer Cell-Mediated Osteoclast Differentiation

BACKGROUND

Metastasis to the bone is a deleterious aspect of breast cancer and is a preferred site that results in bone loss. Hormones such as prolactin (PRL) have not yet been studied for their role in modulating the secondary tumor bone microenvironment.

METHODS

We used quantitative immunohistochemistry with 134 samples of human primary breast cancer and 17 matched primary breast cancers and bone metastases. A Cox proportional hazards regression model was fitted to evaluate the associations between high prolactin receptor (PRLR) expression and time to bone metastasis, adjusting for estrogen receptor status, lymph node status, and chemotherapy status. We assessed osteoclast differentiation, osteoclast size, and measured pit formation in dentine slices. Statistical tests were two-sided.

RESULTS

High PRLR expression in the primary breast tumor was associated with a shorter time to metastasis that includes bone (PRLRAQUA Max-per 100 unit hazard ratio = 1.04, 95% confidence interval = 1.00 to 1.07, P = .03). We observed the PRLR in rare samples of bone metastases and matched primary breast cancer. PRL treatment of breast cancer cells induced osteoclast differentiation and bone lysis via secreted factors and was abrogated by a PRLR antagonist (delta1-9-G129R-hPRL). We demonstrated that sonic hedgehog is a PRL-regulated cytokine in breast cancer cells and part of the mechanism that induces osteoclast differentiation.

CONCLUSIONS

Our evidence indicates that PRL-PRLR can escalate the impact of breast cancer on bone metastasis and that the presence of the PRLR in the tumor microenvironment of breast cancer bone metastasis has the potential to modulate the microenvironment to induce lytic osteoclast formation.

A domestic cat X chromosome linkage map and the sex-linked orange locus: mapping of orange, multiple origins and epistasis over nonagouti

A comprehensive genetic linkage map of the domestic cat X chromosome was generated with the goal of localizing the genomic position of the classic X-linked orange (O) locus. Microsatellite markers with an average spacing of 3 Mb were selected from sequence traces of the cat 1.9x whole genome sequence (WGS), including the pseudoautosomal region 1 (PAR1). Extreme variation in recombination rates (centimorgans per megabase) was observed along the X chromosome, ranging from a virtual absence of recombination events in a region estimated to be >30 Mb to recombination frequencies of 15.7 cM/Mb in a segment estimated to be <0.3 Mb. This detailed linkage map was applied to position the X-linked orange gene, placing this locus on the q arm of the X chromosome, as opposed to a previously reported location on the p arm. Fine mapping placed the locus between markers at positions 106 and 116.8 Mb in the current 1.9x-coverage sequence assembly of the cat genome. Haplotype analysis revealed potential recombination events that could reduce the size of the candidate region to 3.5 Mb and suggested multiple origins for the orange phenotype in the domestic cat. Furthermore, epistasis of orange over nonagouti was demonstrated at the genetic level.

Assessment of AFLP marker behaviour in enriching STS radiation hybrid maps

Radiation hybrid (RH) mapping provides a powerful tool to build high-resolution maps of genomes. Here, we demonstrate the use of the AFLP technique for high-throughput typing of RH cell lines. Cattle were used as the model species because an RH panel was available to investigate the behaviour of AFLP markers within the microsatellite- and STS-based maps of this species. A total of 747 AFLP markers were typed on the TM112 RH radiation panel and 651 of these were assigned by two-point analysis to the 29 bovine autosomes and sex chromosomes. AFLP markers were added to the 1222 microsatellite and STS markers that were included in earlier RH maps. Multipoint maps were constructed for seven example chromosomes, which retained 248 microsatellite and STS markers, and added 123 AFLP markers at LOD 4. The addition of the AFLP markers increased the number of markers by 42.1% and the map length by 10.4%. The AFLP markers showed lower retention frequency (RF) values than the STS markers. The comparison of RF values in AFLP markers and their corresponding AFLP-derived STSs demonstrated that the lower RF values were due to the lower detection sensitivity of the AFLP technique. Despite these differences, AFLP and AFLP-derived STS markers mapped to identical or similar positions. These results demonstrate that it is possible to merge AFLP and microsatellite markers in the same map. The application of AFLP technology could permit the rapid construction of RH maps in species for which extensive genome information and large numbers of SNP and microsatellite markers are not available.

Mapping diabetes QTL in an intercross derived from a congenic strain of the Brown Norway and Goto-Kakizaki rats

Genetic studies in experimental crosses derived from the inbred Goto-Kakizaki (GK) rat model of spontaneous diabetes mellitus have identified quantitative trait loci (QTL) for diabetes phenotypes in a large region of rat Chromosome (RNO) 1. To test the impact of GK variants on QTL statistical and biological features, we combined genetic and physiologic studies in a cohort of F(2) hybrids derived from a QTL substitution congenic strain (QTLSCS) carrying a 110-cM GK haplotype of RNO1 introgressed onto the genetic background of the Brown Norway (BN) strain. Glucose intolerance and altered insulin secretion in QTLSCS rats when compared with BN controls were consistent with original QTL features in a GK x BN F(2) cross. Segregating GK alleles in the QTLSCS F(2) cross account for most of these phenotypic differences between QTLSCS and BN rats. However, significant QTL for diabetes traits in both the QTLSCS and GK x BN F(2) cohorts account for a similar small proportion of their variance. Comparing results from these experimental systems provides indirect estimates of the contribution of genetic interactions and environmental factors to QTL architecture as well as locus and biological targets for future post-QTL mapping studies in congenic substrains.

Chromosomal Mapping of Quantitative Trait Loci Controlling Elastin Content in Rat Aorta

Extracellular matrix molecules such as elastin and collagens provide mechanical support to the vessel wall. In addition to its structural role, elastin is a regulator that maintains homeostasis through biologic signaling. Genetically determined minor modifications in elastin and collagen in the aorta could influence the onset and evolution of arterial pathology, such as hypertension and its complications. We previously demonstrated that the inbred Brown Norway (BN) rat shows an aortic elastin deficit in both abdominal and thoracic segments, partly because of a decrease in tropoelastin synthesis when compared with the LOU rat, that elastin gene polymorphisms in these strains do not significantly account for. After a genome-wide search for quantitative trait loci (QTL) influencing the aortic elastin, collagen, and cell protein contents in an F2 population derived from BN and LOU rats, we identified on chromosomes 2 and 14, 3 QTL specifically controlling elastin levels, and a further highly significant QTL on chromosome 17 linked to the level of cell proteins. We also mapped 3 highly significant QTL linked to body weight (on chromosomes 1 and 3) and heart weight (on chromosome 1) in the cross. This study demonstrates the polygenic control of the content of key components of the arterial wall. Such information represents a first step in understanding possible mechanisms involved in dysregulation of these parameters in arterial pathology.

Quantitative trait locus dissection in congenic strains of the Goto-Kakizaki rat identifies a region conserved with diabetes loci in human chromosome 1q

Genetic studies in human populations and rodent models have identified regions of human chromosome 1q21–25 and rat chromosome 2 showing evidence of significant and replicated linkage to diabetes-related phenotypes. To investigate the relationship between the human and rat diabetes loci, we fine mapped the rat locus Nidd/ gk2 linked to hyperinsulinemia in an F2 cross derived from the diabetic (type 2) Goto-Kakizaki (GK) rat and the Brown Norway (BN) control rat, and carried out its genetic and pathophysiological characterization in BN.GK congenic strains. Evidence of glucose intolerance and enhanced insulin secretion in a congenic strain allowed us to localize the underlying diabetes gene(s) in a rat chromosomal interval of ∼3–6 cM conserved with an 11-Mb region of human 1q21–23. Positional diabetes candidate genes were tested for transcriptional changes between congenics and controls and sequence variations in a panel of inbred rat strains. Congenic strains of the GK rats represent powerful novel models for accurately defining the pathophysiological impact of diabetes gene(s) at the locus Nidd/ gk2 and improving functional annotations of diabetes candidates in human 1q21–23.

Chromosomal Mapping of Genetic Loci Controlling Absence Epilepsy Phenotypes in the WAG/Rij Rat

PURPOSE

The WAG/Rij rat is among the most appropriate models for the study of spontaneous childhood absence epilepsy, without complex neurologic disorders that are associated with some mouse models for absence epilepsy. Previous studies have allowed the identification of distinct types of spike-wave discharges (SWDs) characterizing seizures in this strain. The purpose of this study was to investigate the genetic basis of electroencephalographic (EEG) properties of SWDs.

METHODS

An intercross was derived from WAG/Rij and ACI inbred strains that are known to differ substantially in the number of SWDs. Phenotypic analyses based on 23-h EEG recording in all progenies allowed the quantification of type I and type II SWD phenotypes. A genome-wide scan was performed with 145 microsatellite markers, which were used to test for evidence of genetic linkage to SWD quantitative phenotypes.

RESULTS

We were able to map quantitative trait loci independently, controlling type I and type II SWD variables to rat chromosomes 5 and 9. Strongest linkages were obtained for D5Mgh15 and total duration of type II SWD (lod, 3.64) and for D9Rat103 and the average duration of type I SWD (lod, 3.91). These loci were denoted T2swd/wag and T1swd/wag, respectively.

CONCLUSIONS

The independent genetic control of type I and type II SWDs underlines the complexity of the molecular mechanisms participating in SWDs. The identification of these genetic loci represents an important step in our fundamental knowledge of the architecture of SWDs and may provide new insights for resolving the genetic heterogeneity of absence epilepsy.

High-density rat radiation hybrid maps containing over 24,000 SSLPs, genes, and ESTs provide a direct link to the rat genome sequence

The laboratory rat is a major model organism for systems biology. To complement the cornucopia of physiological and pharmacological data generated in the rat, a large genomic toolset has been developed, culminating in the release of the rat draft genome sequence. The rat draft sequence used a variety of assembly packages, as well as data from the Radiation Hybrid (RH) map of the rat as part of their validation. As part of the Rat Genome Project, we have been building a high-density RH map to facilitate data integration from multiple maps and now to help validate the genome assembly. By incorporating vectors from our lab and several other labs, we have doubled the number of simple sequence length polymorphisms (SSLPs), genes, expressed sequence tags (ESTs), and sequence-tagged sites (STSs) compared to any other genome-wide rat map, a total of 24,437 elements. During the process, we also identified a novel approach for integrating the RH placement results from multiple maps. This new integrated RH map contains approximately 10 RH-mapped elements per Mb on the genome assembly, enabling the RH maps to serve as a scaffold for a variety of data visualization tools.

Integrative genomics: in silico coupling of rat physiology and complex traits with mouse and human data

Integration of the large variety of genome maps from several organisms provides the mechanism by which physiological knowledge obtained in model systems such as the rat can be projected onto the human genome to further the research on human disease. The release of the rat genome sequence provides new information for studies using the rat model and is a key reference against which existing and new rat physiological results can be aligned. Previously, we described comparative maps of the rat, mouse, and human based on EST sequence comparisons combined with radiation hybrid maps. Here, we use new data and introduce the Integrated Genomics Environment, an extensive database of curated and integrated maps, markers, and physiological results. These results are integrated by using VCMapview, a java-based map integration and visualization tool. This unique environment allows researchers to relate results from cytogenetic, genetic, and radiation hybrid studies to the genome sequence and compare regions of interest between human, mouse, and rat. Integrating rat physiology with mouse genetics and clinical results from human by using the respective genomes provides a novel route to capitalize on comparative genomics and the strengths of model organism biology.

Integration of the rat recombination and EST maps in the rat genomic sequence and comparative mapping analysis with the mouse genome

Inbred strains of the laboratory rat are widely used for identifying genetic regions involved in the control of complex quantitative phenotypes of biomedical importance. The draft genomic sequence of the rat now provides essential information for annotating rat quantitative trait locus (QTL) maps. Following the survey of unique rat microsatellite (11,585 including 1648 new markers) and EST (10,067) markers currently available, we have incorporated a selection of 7952 rat EST sequences in an improved version of the integrated linkage-radiation hybrid map of the rat containing 2058 microsatellite markers which provided over 10,000 potential anchor points between rat QTL and the genomic sequence of the rat. A total of 996 genetic positions were resolved (avg. spacing 1.77 cM) in a single large intercross and anchored in the rat genomic sequence (avg. spacing 1.62 Mb). Comparative genome maps between rat and mouse were constructed by successful computational alignment of 6108 mapped rat ESTs in the mouse genome. The integration of rat linkage maps in the draft genomic sequence of the rat and that of other species represents an essential step for translating rat QTL intervals into human chromosomal targets.

Polygenic Control of Idiopathic Generalized Epilepsy Phenotypes in the Genetic Absence Rats from Strasbourg (GAERS)

PURPOSE

Generalized nonconvulsive absence seizures are characterized by the occurrence of synchronous and bilateral spike-and-wave discharges (SWDs) on electroencephalographic recordings, concomitant with behavioral arrest. The GAERS (genetic absence rats from Strasbourg) strain, a well-characterized inbred model for idiopathic generalized epilepsy, spontaneously develops EEG paroxysms that resemble those of typical absence seizures. The purpose of this study was to investigate the genetic control of SWD variables by using a combination of genetic analyses and electrophysiological measurements in an experimental cross derived from GAERS and Brown Norway (BN) rats.

METHODS

SWD subphenotypes were quantified on EEG recordings performed at both 3 and 6 months in a cohort of 118 GAERS x BN F2 animals. A genome-wide scan of the F2 progenies was carried out with 146 microsatellite markers that were used to test each marker locus for evidence of genetic linkage to the SWD quantitative traits.

RESULTS

We identified three quantitative trait loci (QTLs) in chromosomes 4, 7, and 8 controlling specific SWD variables in the cross, including frequency, amplitude, and severity of SWDs. Age was a major factor influencing the detection of genetic linkage to the various components of the SWDs.

CONCLUSIONS

The identification of these QTLs demonstrates the polygenic control of SWDs in the GAERS strain. Genetic linkages to specific SWD features underline the complex mechanisms contributing to SWD development in idiopathic generalized epilepsy.

Comparative recombination rates in the rat, mouse, and human genomes

Levels of recombination vary among species, among chromosomes within species, and among regions within chromosomes in mammals. This heterogeneity may affect levels of diversity, efficiency of selection, and genome composition, as well as have practical consequences for the genetic mapping of traits. We compared the genetic maps to the genome sequence assemblies of rat, mouse, and human to estimate local recombination rates across these genomes. Humans have greater overall levels of recombination, as well as greater variance. In rat and mouse, the size of the chromosome and proximity to telomere have less effect on local recombination rate than in human. At the chromosome level, rat and mouse X chromosomes have the lowest recombination rates, whereas human chromosome X does not show the same pattern. In all species, local recombination rate is significantly correlated with several sequence variables, including GC%, CpG density, repetitive elements, and the neutral mutation rate, with some pronounced differences between species. Recombination rate in one species is not strongly correlated with the rate in another, when comparing homologous syntenic blocks of the genome. This comparative approach provides additional insight into the causes and consequences of genomic heterogeneity in recombination.

Genetic analysis for diabetes in a new rat model of nonobese type 2 diabetes, Spontaneously Diabetic Torii rat

The Spontaneously Diabetic Torii (SDT) rat has recently been established as a new rat model of nonobese type 2 diabetes. In this study, we characterized diabetic features in SDT rats, and performed quantitative trait locus (QTL) analysis for glucose intolerance using 319 male (BNxSDT)xSDT backcrosses. Male SDT rats exhibited glucose intolerance at 20 weeks, and spontaneously developed diabetes with the incidence of 100% at 38 weeks, and glucose intolerance is well associated with the development of diabetes. The QTL analysis identified three highly significant QTLs (Gisdt1, Gisdt2, and Gisdt3) for glucose intolerance on rat chromosomes 1, 2, and X, respectively. The SDT allele for these QTLs significantly exacerbated glucose intolerance. Furthermore, synergistic interactions among these QTLs were detected. These findings indicate that diabetic features in SDT rats are inherited as polygenic traits and that SDT rats would provide insights into genetics of human type 2 diabetes.

Further dissection of a genomic locus associated with behavioral activity in the Wistar-Kyoto hyperactive rat, an animal model of hyperkinesis

Molecular genetic studies of attention-deficit hyperactivity disorder (ADHD) are a major focus of current research since this syndrome has been shown to be highly heritable.(1) Our approach has been to search for quantitative trait loci (QTL) in a genetic animal model of hyperkinesis, the Wistar-Kyoto hyperactive (WKHA) rat, by a whole-genome scan analysis. In a previous article, we reported the detection of a major QTL associated with behavioral activity in an F2 cross between WKHA and Wistar-Kyoto (WKY) rat strains.(2) Here, we extend our analysis of this cross by adding new genetic markers, now defining a 10 cM interval on rat chromosome 8 associated with ambulatory and exploratory activities. Then we present a replication of this QTL detection, at least for exploratory activity, by a new genetic mapping analysis of an activity QTL in an F2 cross between the WKHA and Brown Norway (BN) rat strains. Overall, the results provide compelling evidence for the presence of gene(s) influencing activity at this locus. The QTL interval has been refined such that the human orthologous region could be defined and tested in human populations for association with ADHD. Ultimately, the improved dissection of this genomic locus should allow the identification of the causal genes.

Mapping of genetic loci predisposing to hypertriglyceridaemia in the hereditary hypertriglyceridaemic rat: analysis of genetic association with related traits of the insulin resistance syndrome

AIMS/HYPOTHESIS

Hypertriglyceridaemia is an important risk factor for coronary heart disease, especially in the context of the insulin resistance syndrome where it often occurs with hypertension. The two phenotypes are also associated in the hereditary hypertriglyceridaemic (hHTg) rat. The aim of this study was to map quantitative trait loci that affect plasma triglyceride concentration in the hHTg rat and determine whether they co-localize with loci for blood pressure.

METHODS

Second filial generation progeny (n=189) from a cross of the hHTg rat with the Brown Norway rat were phenotyped for fasting plasma triglyceride, glucose and insulin concentrations, and direct unrestrained resting blood pressure. A partial genome-scan was conducted using 153 microsatellite markers that were polymorphic between the two strains.

RESULTS

A major locus (lod score 6.5) influencing plasma triglyceride concentration in a co-dominant fashion was mapped to chromosome 4 between D1Mit 5 and D1Mit17. Chromosome 8 contained multiple peaks with a lod score greater than 4.0 influencing triglyceride concentration. Importantly, none of the triglyceride loci had an effect on blood pressure. The triglyceride locus on chromosome 4 co-localized with a locus for fasting plasma insulin (lod score 4.1), although the effect on insulin concentration was in the opposite direction to that on triglyceride.

CONCLUSION/INTERPRETATION

We have mapped the major loci that affect plasma triglyceride concentration in the hHTg rat. These loci do not influence blood pressure suggesting that these commonly associated phenotypes of the insulin resistance syndrome are not be due to pleiotropic effects of the same gene(s).

Human chromosome 17 in essential hypertension

Hypertension affects up to 30% of the adult population in Western societies and is a major risk factor for kidney disease, stroke and coronary heart disease. It is a complex trait thought to be influenced by a number of genes and environmental factors, although the precise aetiology remains unknown at this time. A number of methods have been successfully used to identify mutations that cause Mendelian traits and these are now being applied to the investigation of complex diseases. This review summarises the data gathered, using such approaches, that suggest there is a gene or genes on chromosome 17 causing human essential hypertension. Studies in rodent models are discussed first, followed by studies of human hypertension that include the investigation of pseudohypoaldosteronism type II, a monogenic trait that manifests with hypertension alongside other phenotypic variables. In addition, candidate gene studies, genome screens and linkage studies based on comparative mapping are outlined. To date no gene has been identified on human chromosome 17 that influences blood pressure and causes human essential hypertension. However, results of ongoing fine mapping and candidate gene studies in both rodents and man are eagerly awaited.

Congenic mapping of a blood pressure QTL on chromosome 16 of Dahl rats

A Chromosome (Chr) 16 segment of the Dahl salt-sensitive (S) rat was shown by linkage to contain a blood pressure (BP) quantitative trait locus (QTL). To verify and further narrow down the region harboring the QTL, we made two congenic strains by replacing two segments of the S rats with the homologous segments of the Lewis (LEW) rats. The construction of these congenic strains was facilitated by a genome-wide marker screening. The two congenic strains contained a segment in common, and BPs of both were significantly lower than that of the S strain. Consequently, a BP QTL could be localized to the overlapping region of about 49.4 centiRay (cR) including the telomere on a radiation hybrid map. Heart weights, left and right ventricular weights, kidney weights, and aortic weights over length were all significantly decreased in the congenic strains compared with the S strain. Thus, there appeared to exist an association between the effects of the QTL on BP and on cardiac, renal, and vascular hypertrophy.