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

The Rat Genome Database (RGD) facilitates genomic and phenotypic data integration across multiple species for biomedical research

Model organism research is essential for discovering the mechanisms of human diseases by defining biologically meaningful gene to disease relationships. The Rat Genome Database (RGD, ( https://rgd.mcw.edu )) is a cross-species knowledgebase and the premier online resource for rat genetic and physiologic data. This rich resource is enhanced by the inclusion and integration of comparative data for human and mouse, as well as other human disease models including chinchilla, dog, bonobo, pig, 13-lined ground squirrel, green monkey, and naked mole-rat. Functional information has been added to records via the assignment of annotations based on sequence similarity to human, rat, and mouse genes. RGD has also imported well-supported cross-species data from external resources. To enable use of these data, RGD has developed a robust infrastructure of standardized ontologies, data formats, and disease- and species-centric portals, complemented with a suite of innovative tools for discovery and analysis. Using examples of single-gene and polygenic human diseases, we illustrate how data from multiple species can help to identify or confirm a gene as involved in a disease and to identify model organisms that can be studied to understand the pathophysiology of a gene or pathway. The ultimate aim of this report is to demonstrate the utility of RGD not only as the core resource for the rat research community but also as a source of bioinformatic tools to support a wider audience, empowering the search for appropriate models for human afflictions.

The porcine fat mass and obesity associated (FTO) gene is associated with fat deposition in Italian Duroc pigs

In humans, common variants in the fat mass and obesity associated (FTO) gene are associated with body mass index and obesity. Here we sequenced exon 4, parts of introns 3 and 4 and two portions of the 3'-untranslated region of the porcine FTO gene in a panel of nine pigs of different breeds and identified three SNPs. Allele frequencies of the g.276T>G (AM931150) mutation were studied in seven pig breeds. This mutation was used to linkage-map FTO to SSC6. Association analyses between the g.276T>G polymorphism and several traits [pH of semimembranosus muscle and estimated breeding values (EBV) for average daily gain, back fat thickness, lean cuts, ham weight and feed:gain ratio] were carried out in 257 sib-tested Italian Large White pigs. Only feed:gain ratio showed P<0.05. A selective genotyping approach was applied, analysing two extreme and divergent groups of Italian Large White pigs selected on the basis of back fat thickness EBV (50 with most positive and 50 with most negative values). Fisher's exact test (two-tailed) was not significant when comparing the allele frequencies of these two groups. The same approach was used in the Italian Duroc breed for which two extreme and divergent groups of animals were selected according to visible intermuscular fat EBV. Differences of allele frequencies between these two groups were highly significant (P<0.00001, P<0.001 and P<0.0001, considering all animals or only two- or three-generation unrelated animals respectively), indicating association between the analysed FTO marker and intermuscular fat deposition.

Multiple blood pressure loci on rat chromosome 13 attenuate development of hypertension in the Dahl S hypertensive rat

Previous studies have indicated that substitution of chromosome 13 of the salt-resistant Brown Norway BN/SsNHsdMcwi (BN) rat into the genomic background of the Dahl salt-sensitive SS/JrHsdMcwi (SS) rat attenuates the development of salt-sensitive hypertension and renal damage. To identify the regions within chromosome 13 that attenuate the development of hypertension during a high-salt diet in the SS rat, we phenotyped a series of overlapping congenic lines covering chromosome 13, generated from an intercross between the consomic SS-13BN rat and the SS rat. Blood pressure was determined in chronically catheterized rats after 2 wk of high-salt diet (8% NaCl) together with microalbuminuria as an index of renal damage. Four discrete regions were identified, ranging in size from 4.5 to 16 Mbp, each of which independently provided significant protection from hypertension during high-salt diet, reducing blood pressure by 20–29 mmHg. Protection was more robust in female than male rats in some of the congenic strains, suggesting a sex interaction with some of the genes determining blood pressure during high-salt diet. Among the 23 congenic strains, several regions overlapped. When three of the “protective” regions were combined onto one broad congenic strain, no summation effect was seen, obtaining the same decrease in blood pressure as with each one independently. We conclude from these studies that there are four regions within chromosome 13 containing genes that interact epistatically and influence arterial pressure.

Proteomic profiling of animal models mimicking skeletal muscle disorders

Over the last few decades of biomedical research, animal models of neuromuscular diseases have been widely used for determining pathological mechanisms and for testing new therapeutic strategies. With the emergence of high-throughput proteomics technology, the identification of novel protein factors involved in disease processes has been decisively improved. This review outlines the usefulness of the proteomic profiling of animal disease models for the discovery of new reliable biomarkers, for the optimization of diagnostic procedures and the development of new treatment options for skeletal muscle disorders. Since inbred animal strains show genetically much less interindividual differences as compared to human patients, considerably lower experimental repeats are capable of producing meaningful proteomic data. Thus, animal model proteomics can be conveniently employed for both studying basic mechanisms of molecular pathogenesis and the effects of drugs, genetic modifications or cell-based therapies on disease progression. Based on the results from comparative animal proteomics, a more informed decision on the design of clinical proteomics studies could be reached. Since no one animal model represents a perfect pathobiochemical replica of all of the symptoms seen in complex human disorders, the proteomic screening of novel animal models can also be employed for swift and enhanced protein biochemical phenotyping.

The Rat Genome Database, update 2007--easing the path from disease to data and back again

The Rat Genome Database (RGD, http://rgd.mcw.edu) is one of the core resources for rat genomics and recent developments have focused on providing support for disease-based research using the rat model. Recognizing the importance of the rat as a disease model we have employed targeted curation strategies to curate genes, QTL and strain data for neurological and cardiovascular disease areas. This work has centered on rat but also includes data for mouse and human to create 'disease portals' that provide a unified view of the genes, QTL and strain models for these diseases across the three species. The disease curation efforts combined with normal curation activities have served to greatly increase the content of the database, particularly for biological information, including gene ontology, disease, pathway and phenotype ontology annotations. In addition to improving the features and database content, community outreach has been expanded to demonstrate how investigators can leverage the resources at RGD to facilitate their research and to elicit suggestions and needs for future developments. We have published a number of papers that provide additional information on the ontology annotations and the tools at RGD for data mining and analysis to better enable researchers to fully utilize the database.

Physiogenomic resources for rat models of heart, lung and blood disorders

Cardiovascular disorders are influenced by genetic and environmental factors. The TIGR rodent expression web-based resource (TREX) contains over 2,200 microarray hybridizations, involving over 800 animals from 18 different rat strains. These strains comprise genetically diverse parental animals and a panel of chromosomal substitution strains derived by introgressing individual chromosomes from normotensive Brown Norway (BN/NHsdMcwi) rats into the background of Dahl salt sensitive (SS/JrHsdMcwi) rats. The profiles document gene-expression changes in both genders, four tissues (heart, lung, liver, kidney) and two environmental conditions (normoxia, hypoxia). This translates into almost 400 high-quality direct comparisons (not including replicates) and over 100,000 pairwise comparisons. As each individual chromosomal substitution strain represents on average less than a 5% change from the parental genome, consomic strains provide a useful mechanism to dissect complex traits and identify causative genes. We performed a variety of data-mining manipulations on the profiles and used complementary physiological data from the PhysGen resource to demonstrate how TREX can be used by the cardiovascular community for hypothesis generation.

Interaction between Rf-1 and Rf-4 quantitative trait loci increases susceptibility to renal damage in double congenic rats

BACKGROUND

Five quantitative trait loci (QTLs), Rf-1 to Rf-5, were found in Fawn-Hooded hypertensive (FHH) rats influencing susceptibility to renal damage. Previously, we found that single transfer of the Rf-1 QTL from FHH rats onto the renal-resistant August x Copenhagen Irish (ACI) strain caused a small increase in renal susceptibility. To investigate the separate role of the Rf-4 QTL and its interaction with Rf-1, we generated a single congenic strain carrying Rf-4 and a double congenic carrying both Rf-1 and Rf-4.

METHODS

Differences in renal susceptibility between ACI, Rf-1A, and Rf-4 single congenics and Rf-1A+4 double congenics were assessed using four different treatments: control (two-kidney), two-kidney with l-arginine analogue N-nitro-l-arginine methyl ester (L-NAME)-induced hypertension, unilateral nephrectomy, and unilateral nephrectomy + L-NAME. In separate experiments, renal blood flow (RBF) autoregulation was compared between two-kidney ACI and congenic rats.

RESULTS

Compared to ACI, Rf-1A rats developed more renal damage, while Rf-4 rats did not. The most severe renal damage was found in the Rf-1A+4 double congenic rats. Analysis of variance (ANOVA) demonstrated a significant interaction between the Rf-1A and Rf-4 QTLs. The magnitude of the interaction varied with the type and duration of the treatment. The RBF autoregulation was impaired in Rf-1A single and Rf-1A+4 double congenics, while in Rf-4 single congenics it was similar to that of ACI controls.

CONCLUSION

These findings indicate that the Rf-1 QTL directly influences renal susceptibility and autoregulation. In contrast, the Rf-4 QTL shows no direct effects, but significantly increases susceptibility to renal damage via an interaction with Rf-1.

Integrating Genetic and Gene Expression Data to Study the Metabolic Syndrome and Diabetes in Mice

Increasingly, the mouse is becoming the standard model for mammalian physiology and disease. It can be genetically analyzed and manipulated with relative ease. Moreover, the endogenous genetic variation that exists among inbred mouse strains can be exploited to identify genetic control of complex physiologic processes involved in diabetes and the metabolic syndrome, among other conditions relevant to human disease. Recent advances in genetics and gene expression technology have greatly increased the knowledge to be derived from this approach when applied to traditional genetic studies.

Tools and strategies for physiological genomics: the Rat Genome Database

The broad goal of physiological genomics research is to link genes to their functions using appropriate experimental and computational techniques. Modern genomics experiments enable the generation of vast quantities of data, and interpretation of this data requires the integration of information derived from many diverse sources. Computational biology and bioinformatics offer the ability to manage and channel this information torrent. The Rat Genome Database (RGD; http://rgd.mcw.edu) has developed computational tools and strategies specifically supporting the goal of linking genes to their functional roles in rat and, using comparative genomics, to human and mouse. We present an overview of the database with a focus on these unique computational tools and describe strategies for the use of these resources in the area of physiological genomics.

Integrative genomics based identification of potential human hepatocarcinogenesis-associated cell cycle regulators: RHAMM as an example

DNA microarray has been widely used to examine gene expression profile of different human tumors. The information generated from microarray analysis usually represents the overall range of cancer-associated abnormality associated with gene regulation. In order to identify key regulatory genes involved in carcinogenesis of human cancer, hypothesis driven data mining of the microarray data plus experimental validation becomes a critical approach in the post-genome era. Here, we present an integrative genomic analysis of published microarray data and homolog gene database. Over 20,000 genes were examined to reveal 16 genes specific to vertebrates, cell cycle G2/M regulated, and overexpressed in human HCC. Using Affymetrix microarray analysis, we found that all 16 genes were up-regulated in human HCC. Among these 16 genes, we experimentally validated the up-regulation of receptor for hyaluronan-mediated motility (RHAMM) in different cell model systems. We first confirmed elevation of RHAMM in the G2/M phase of synchronized HeLa cells. We also found that RHAMM had an elevated level of expression in all the HCC samples we examined and it was induced during the G2/M phase of regenerating mouse hepatocytes after partial hepatectomy. Thus, the expression of RHAMM appears to be tightly regulated during mammalian cell cycle G2/M progression. The ectopic overexpression of RHAMM in 293T cells resulted in the accumulation of cells at G2/M phase. RHAMM-induced mitotic arrest of cells was predominantly in the prophase. Taken together, using an integrated functional genomic approach, we have uncovered a set of genes that may play specific roles in cell cycle progression and in HCC development. To elucidate the function of these genes in cell cycle regulation may shed light on the control mechanism of human HCC in the future.

The Rat Genome Database (RGD): developments towards a phenome database

The Rat Genome Database (RGD) (http://rgd.mcw.edu) aims to meet the needs of its community by providing genetic and genomic infrastructure while also annotating the strengths of rat research: biochemistry, nutrition, pharmacology and physiology. Here, we report on RGD's development towards creating a phenome database. Recent developments can be categorized into three groups. (i) Improved data collection and integration to match increased volume and biological scope of research. (ii) Knowledge representation augmented by the implementation of a new ontology and annotation system. (iii) The addition of quantitative trait loci data, from rat, mouse and human to our advanced comparative genomics tools, as well as the creation of new, and enhancement of existing, tools to enable users to efficiently browse and survey research data. The emphasis is on helping researchers find genes responsible for disease through the use of rat models. These improvements, combined with the genomic sequence of the rat, have led to a successful year at RGD with over two million page accesses that represent an over 4-fold increase in a year. Future plans call for increased annotation of biological information on the rat elucidated through its use as a model for human pathobiology. The continued development of toolsets will facilitate integration of these data into the context of rat genomic sequence, as well as allow comparisons of biological and genomic data with the human genomic sequence and of an increasing number of organisms.

The Visual Language of Synteny

The study of polygenic disorders such as cardiovascular and metabolic diseases requires access to vast amounts of experimental and in silico data. Where animal models of disease are being used, visualization of syntenic genome regions is one of the most important tools supporting data analysis. We define what is required to visualize synteny in terms of the data being displayed, the screen layout, and user interaction. We then describe a prototype visualization tool, SyntenyVista, which provides integrated access to quantitative trait loci, microarray, and gene datasets. We believe that SyntenyVista is a significant step towards an improved representation of comparative genomics data.

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.