Genomic structure of nucleotide diversity among Lyon rat models of metabolic syndrome

A revamped rat reference genome improves the discovery of genetic diversity in laboratory rats

The seventh iteration of the reference genome assembly for Rattus norvegicus-mRatBN7.2-corrects numerous misplaced segments and reduces base-level errors by approximately 9-fold and increases contiguity by 290-fold compared with its predecessor. Gene annotations are now more complete, improving the mapping precision of genomic, transcriptomic, and proteomics datasets. We jointly analyzed 163 short-read whole-genome sequencing datasets representing 120 laboratory rat strains and substrains using mRatBN7.2. We defined ∼20.0 million sequence variations, of which 18,700 are predicted to potentially impact the function of 6,677 genes. We also generated a new rat genetic map from 1,893 heterogeneous stock rats and annotated transcription start sites and alternative polyadenylation sites. The mRatBN7.2 assembly, along with the extensive analysis of genomic variations among rat strains, enhances our understanding of the rat genome, providing researchers with an expanded resource for studies involving rats.

A revamped rat reference genome improves the discovery of genetic diversity in laboratory rats

The seventh iteration of the reference genome assembly for Rattus norvegicus-mRatBN7.2-corrects numerous misplaced segments and reduces base-level errors by approximately 9-fold and increases contiguity by 290-fold compared to its predecessor. Gene annotations are now more complete, significantly improving the mapping precision of genomic, transcriptomic, and proteomics data sets. We jointly analyzed 163 short-read whole genome sequencing datasets representing 120 laboratory rat strains and substrains using mRatBN7.2. We defined ~20.0 million sequence variations, of which 18.7 thousand are predicted to potentially impact the function of 6,677 genes. We also generated a new rat genetic map from 1,893 heterogeneous stock rats and annotated transcription start sites and alternative polyadenylation sites. The mRatBN7.2 assembly, along with the extensive analysis of genomic variations among rat strains, enhances our understanding of the rat genome, providing researchers with an expanded resource for studies involving rats.

Body Composition and Metabolic Changes in a Lyon Hypertensive Congenic Rat and Identification of Ercc6l2 as a Positional Candidate Gene

Central obesity is genetically complex, and its exponential increase in the last decades have made it a critical public health issue. The Lyon Hypertensive (LH) rat is a well-characterized hypertensive model that also exhibits spontaneous and profound differences in body weight and adiposity, relative to its metabolically healthy control, the Lyon Normotensive (LN) rat. The mechanisms underlying the body weight differences between these strains are not well-understood, thus a congenic model (LH17LNa) was developed where a portion of the proximal arm of LN chromosome 17 is introgressed on the LH genomic background to assess the contribution of LN alleles on obesity features. Male and female LH17LNa rats were studied, but male congenics did not significantly differ from LH in this study. Female LH17LNa rats exhibited decreases in total body growth, as well as major alterations to their body composition and adiposity. The LH17LNa female rats also showed decreases in metabolic rate, and a reduction in food intake. The increased adiposity in the female LH17LNa rats was specific to abdominal white adipose tissue, and this phenomenon was further explained by significant hypertrophy in those adipocytes, with no evidence of adipocyte hyperplasia. Sequencing of the parental strains identified a novel frameshift mutation in the candidate gene Ercc6l2, which is involved in transcription-coupled DNA repair, and is implicated in premature aging. The discovery of the significance of Ercc6l2 in the context of female-specific adipocyte biology could represent a novel role of DNA repair failure syndromes in obesity pathogenesis.

Genome variation and conserved regulation identify genomic regions responsible for strain specific phenotypes in rat

Background

The genomes of laboratory rat strains are characterised by a mosaic haplotype structure caused by their unique breeding history. These mosaic haplotypes have been recently mapped by extensive sequencing of key strains. Comparison of genomic variation between two closely related rat strains with different phenotypes has been proposed as an effective strategy for the discovery of candidate strain-specific regions involved in phenotypic differences. We developed a method to prioritise strain-specific haplotypes by integrating genomic variation and genomic regulatory data predicted to be involved in specific phenotypes. Specifically, we aimed to identify genomic regions associated with Metabolic Syndrome (MetS), a disorder of energy utilization and storage affecting several organ systems.

Results

We compared two Lyon rat strains, Lyon Hypertensive (LH) which is susceptible to MetS, and Lyon Low pressure (LL), which is susceptible to obesity as an intermediate MetS phenotype, with a third strain (Lyon Normotensive, LN) that is resistant to both MetS and obesity. Applying a novel metric, we ranked the identified strain-specific haplotypes using evolutionary conservation of the occupancy three liver-specific transcription factors (HNF4A, CEBPA, and FOXA1) in five rodents including rat. Consideration of regulatory information effectively identified regions with liver-associated genes and rat orthologues of human GWAS variants related to obesity and metabolic traits. We attempted to find possible causative variants and compared them with the candidate genes proposed by previous studies. In strain-specific regions with conserved regulation, we found a significant enrichment for published evidence to obesity—one of the metabolic symptoms shown by the Lyon strains—amongst the genes assigned to promoters with strain-specific variation.

Conclusions

Our results show that the use of functional regulatory conservation is a potentially effective approach to select strain-specific genomic regions associated with phenotypic differences among Lyon rats and could be extended to other systems.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4351-9) contains supplementary material, which is available to authorized users.

Contribution of independent and pleiotropic genetic effects in the metabolic syndrome in a hypertensive rat

Hypertension is a major risk factor for cardiovascular disease, Type 2 diabetes, and end organ failure, and is often found concomitant with disorders characteristic of the Metabolic Syndrome (MetS), including obesity, dyslipidemia, and insulin resistance. While the associated features often occur together, the pathway(s) or mechanism(s) linking hypertension in MetS are not well understood. Previous work determined that genetic variation on rat chromosome 17 (RNO17) contributes to several MetS-defining traits (including hypertension, obesity, and dyslipidemia) in the Lyon Hypertensive (LH) rat, a genetically determined MetS model. We hypothesized that at least some of the traits on RNO17 are controlled by a single gene with pleiotropic effects. To address this hypothesis, consomic and congenic strains were developed, whereby a defined fragment of RNO17 from the LH rat was substituted with the control Lyon Normotensive (LN) rat, and MetS phenotypes were measured in the resultant progeny. Compared to LH rats, LH-17^LN consomic rats have significantly reduced body weight, blood pressure, and lipid profiles. A congenic strain (LH-17^LNc), with a substituted fragment at the distal end of RNO17 (17q12.3; 74–97 Mb; rn4 assembly), showed differences from the LH rat in blood pressure and serum total cholesterol and triglycerides. Interestingly, there was no difference in body weight between the LH-17^LNc and the parental LH rat. These data indicate that blood pressure and serum lipids are regulated by a gene(s) in the distal congenic interval, and could be due to pleiotropy. The data also indicate that body weight is not determined by the same gene(s) at this locus. Interestingly, only two small haplotypes spanning a total of approximately 0.5 Mb differ between the LH and LN genomes in the congenic interval. Genes in these haplotypes are strong candidate genes for causing dyslipidemia in the LH rat. Overall, MetS, even in a simplified genetic model such as the LH-17^LN rat, is likely due to both independent and pleiotropic gene effects.

Strain-Specific Single-Nucleotide Polymorphisms in Hypertensive ISIAH Rats

Single-nucleotide polymorphisms (SNPs) in the coding and regulatory regions of genes can affect transcription rate and translation efficiency, modify protein function, and, in some cases, cause the development of diseases. In the current study, the RNA-Seq approach has been used to discover strain-specific SNPs in ISIAH (inherited stress-induced arterial hypertension) rats, which are known as a model of stress-induced arterial hypertension. The comparison of the ISIAH SNPs with genome sequencing data available for another 42 rat strains and substrains, 11 of them known as hypertensive, showed a considerable genetic distance between the genotypes of ISIAH and all other rat strains and substrains. The study revealed 1849 novel SNPs specific for ISIAH rats and 158 SNPs present only in the genotypes of hypertensive rats. Amino acid substitutions with possible deleterious effect on protein function were detected. Several of them were found in the genes associated with hypertension. These SNPs may be considered as novel molecular targets for further studies aimed at assessing their potential in the therapy of stress-induced hypertension.

The Genomic Scrapheap Challenge; Extracting Relevant Data from Unmapped Whole Genome Sequencing Reads, Including Strain Specific Genomic Segments, in Rats

Unmapped next-generation sequencing reads are typically ignored while they contain biologically relevant information. We systematically analyzed unmapped reads from whole genome sequencing of 33 inbred rat strains. High quality reads were selected and enriched for biologically relevant sequences; similarity-based analysis revealed clustering similar to previously reported phylogenetic trees. Our results demonstrate that on average 20% of all unmapped reads harbor sequences that can be used to improve reference genomes and generate hypotheses on potential genotype-phenotype relationships. Analysis pipelines would benefit from incorporating the described methods and reference genomes would benefit from inclusion of the genomic segments obtained through these efforts.

Genomic landscape of rat strain and substrain variation

Background

Since the completion of the rat reference genome in 2003, whole-genome sequencing data from more than 40 rat strains have become available. These data represent the broad range of strains that are used in rat research including commonly used substrains. Currently, this wealth of information cannot be used to its full extent, because the variety of different variant calling algorithms employed by different groups impairs comparison between strains. In addition, all rat whole genome sequencing studies to date used an outdated reference genome for analysis (RGSC3.4 released in 2004).

Results

Here we present a comprehensive, multi-sample and uniformly called set of genetic variants in 40 rat strains, including 19 substrains. We reanalyzed all primary data using a recent version of the rat reference assembly (RGSC5.0 released in 2012) and identified over 12 million genomic variants (SNVs, indels and structural variants) among the 40 strains. 28,318 SNVs are specific to individual substrains, which may be explained by introgression from other unsequenced strains and ongoing evolution by genetic drift. Substrain SNVs may have a larger predicted functional impact compared to older shared SNVs.

Conclusions

In summary we present a comprehensive catalog of uniformly analyzed genetic variants among 40 widely used rat inbred strains based on the RGSC5.0 assembly. This represents a valuable resource, which will facilitate rat functional genomic research. In line with previous observations, our genome-wide analyses do not show evidence for contribution of multiple ancestral founder rat subspecies to the currently used rat inbred strains, as is the case for mouse. In addition, we find that the degree of substrain variation is highly variable between strains, which is of importance for the correct interpretation of experimental data from different labs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1594-1) contains supplementary material, which is available to authorized users.

Systems biology with high-throughput sequencing reveals genetic mechanisms underlying the metabolic syndrome in the Lyon hypertensive rat

BACKGROUND

The metabolic syndrome (MetS) is a collection of co-occurring complex disorders including obesity, hypertension, dyslipidemia, and insulin resistance. The Lyon hypertensive and Lyon normotensive rats are models of MetS sensitivity and resistance, respectively. To identify genetic determinants and mechanisms underlying MetS, an F2 intercross between Lyon hypertensive and Lyon normotensive was comprehensively studied.

METHODS AND RESULTS

Multidimensional data were obtained including genotypes of 1536 single-nucleotide polymorphisms, 23 physiological traits, and >150 billion nucleotides of RNA-seq reads from the livers of F2 intercross offspring and parental rats. Phenotypic and expression quantitative trait loci (eQTL) were mapped. Application of systems biology methods identified 17 candidate MetS genes. Several putative causal cis-eQTL were identified corresponding with phenotypic QTL loci. We found an eQTL hotspot on rat chromosome 17 that is causally associated with multiple MetS-related traits and found RGD1562963, a gene regulated in cis by this eQTL hotspot, as the most likely eQTL driver gene directly affected by genetic variation between Lyon hypertensive and Lyon normotensive rats.

CONCLUSIONS

Our study sheds light on the intricate pathogenesis of MetS and demonstrates that systems biology with high-throughput sequencing is a powerful method to study the pathogenesis of complex genetic diseases.