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

The Rat Genome Database: Genetic, Genomic, and Phenotypic Data Across Multiple Species

The laboratory rat, Rattus norvegicus, is an important model of human health and disease, and experimental findings in the rat have relevance to human physiology and disease. The Rat Genome Database (RGD, https://rgd.mcw.edu) is a model organism database that provides access to a wide variety of curated rat data including disease associations, phenotypes, pathways, molecular functions, biological processes, cellular components, and chemical interactions for genes, quantitative trait loci, and strains. We present an overview of the database followed by specific examples that can be used to gain experience in employing RGD to explore the wealth of functional data available for the rat and other species. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Navigating the Rat Genome Database (RGD) home page Basic Protocol 2: Using the RGD search functions Basic Protocol 3: Searching for quantitative trait loci Basic Protocol 4: Using the RGD genome browser (JBrowse) to find phenotypic annotations Basic Protocol 5: Using OntoMate to find gene-disease data Basic Protocol 6: Using MOET to find gene-ontology enrichment Basic Protocol 7: Using OLGA to generate gene lists for analysis Basic Protocol 8: Using the GA tool to analyze ontology annotations for genes Basic Protocol 9: Using the RGD InterViewer tool to find protein interaction data Basic Protocol 10: Using the RGD Variant Visualizer tool to find genetic variant data Basic Protocol 11: Using the RGD Disease Portals to find disease, phenotype, and other information Basic Protocol 12: Using the RGD Phenotypes & Models Portal to find qualitative and quantitative phenotype data and other rat strain-related information Basic Protocol 13: Using the RGD Pathway Portal to find disease and phenotype data via molecular pathways.

Chromosomal Substitution Strategies to Localize Genomic Regions Related to Complex Traits

Chromosomal substitution strategies provide a powerful tool to anonymously reveal the relationship between DNA sequence variants and a normal or disease phenotype of interest. Even in this age of CRISPR-Cas9 genome engineering, the knockdown or overexpression of a gene provides relevant information to our understanding of complex disease only when a close association of an allelic variant with the phenotype has first been established. Limitations of genetic linkage approaches led to the development of more efficient breeding strategies to substitute chromosomal segments from one animal strain into the genetic background of a different strain, enabling a direct comparison of the phenotypes of the strains with variant(s) that differ only at a defined locus. This substitution can be a whole chromosome (consomic), a part of a chromosome (congenic), or as small as only a single or several alleles (subcongenics). In contrast to complete knockout of a specific candidate gene of interest, which simply studies the effects of complete elimination of the gene, the substitution of naturally occurring variants can provide special insights into the functional actions of wild-type alleles. Strategies for production of these inbred strains are reviewed, and a number of examples are used to illustrate the utility of these model systems. Consomic/congenic strains provide a number of experimental advantages in the study of functions of genes and their variants, which are emphasized in this article, such as replication of experimental studies; determination of temporal relationships throughout a life; rigorously controlled experiments in which relations between genotype and phenotype can be tested with the confounding effects of heterogeneous genetic backgrounds, both targeted and multilayered; and "omic" studies performed at many levels of functionality, from molecules to organelles, cells to organs, and organs to organismal behavior across the life span. The application of chromosomal substitution strategies and development of consomic/congenic rat and mouse strains have greatly expanded our knowledge of genomic variants and their phenotypic relationship to physiological functions and to complex diseases such as hypertension and cancer. © 2020 American Physiological Society. Compr Physiol 10:365-388, 2020.

Wars2 is a determinant of angiogenesis

Coronary flow (CF) measured ex vivo is largely determined by capillary density that reflects angiogenic vessel formation in the heart in vivo. Here we exploit this relationship and show that CF in the rat is influenced by a locus on rat chromosome 2 that is also associated with cardiac capillary density. Mitochondrial tryptophanyl-tRNA synthetase (Wars2), encoding an L53F protein variant within the ATP-binding motif, is prioritized as the candidate at the locus by integrating genomic data sets. WARS2(L53F) has low enzyme activity and inhibition of WARS2 in endothelial cells reduces angiogenesis. In the zebrafish, inhibition of wars2 results in trunk vessel deficiencies, disordered endocardial-myocardial contact and impaired heart function. Inhibition of Wars2 in the rat causes cardiac angiogenesis defects and diminished cardiac capillary density. Our data demonstrate a pro-angiogenic function for Wars2 both within and outside the heart that may have translational relevance given the association of WARS2 with common human diseases.

Blood supply to the heart is crucial for cardiac function. Here, the authors show that the mitochondrial tryptophanyl-tRNA synthetase, WARS2, drives blood vessel generation in zebrafish and rats and that inhibition of Wars2 diminishes blood vessel growth both within and outside in the heart, suggesting a new target for manipulating angiogenesis.

Epigenomics of hypertension

Multiple genes and pathways are involved in the pathogenesis of hypertension. Epigenomic studies of hypertension are beginning to emerge and hold great promise of providing novel insights into the mechanisms underlying hypertension. Epigenetic marks or mediators including DNA methylation, histone modifications, and noncoding RNA can be studied at a genome or near-genome scale using epigenomic approaches. At the single gene level, several studies have identified changes in epigenetic modifications in genes expressed in the kidney that correlate with the development of hypertension. Systematic analysis and integration of epigenetic marks at the genome-wide scale, demonstration of cellular and physiological roles of specific epigenetic modifications, and investigation of inheritance are among the major challenges and opportunities for future epigenomic and epigenetic studies of hypertension.

Modeling inflammatory bowel disease: the zebrafish as a way forward

The zebrafish has proved to be an informative model of vertebrate development and, more recently, an emerging model of human disease. The realization of the full potential of the zebrafish as a disease model lies in two interdependent areas. The first is an appreciation that the often overlooked strength of this species lies in allowing the design of experiments that address the interplay of genetics and the environment in a manipulable manner. The second is in the application and further development of gene targeting approaches. These twin features will be addressed in this review in the context of modeling inflammatory bowel disease.

Enoyl coenzyme a hydratase domain-containing 2, a potential novel regulator of myocardial ischemia injury

Background

We reported previously that Brown Norway (BN) rats are more resistant to myocardial ischemia/reperfusion (I/R) injury than are Dahl S (SS) rats. To identify the unique genes differentially expressed in the hearts of these rats, we used DNA microarray analysis and observed that enoyl coenzyme A hydratase–containing domain 2 (ECHDC2) is highly expressed (≈18‐fold) in the SS hearts compared with the BN hearts.

Methods and Results

RT‐PCR, Western blot, and immunohistochemistry analyses verified that ECHDC2 was highly expressed in SS hearts compared with the BN hearts. ECHDC2 gene locates at chromosome 5 of rat and is expressed in mitochondria of the heart, mainly in cardiomyocytes but not in cardiofibroblasts. Overexpression of ECHDC2 in cells increased susceptibility to I/R injury while knockdown of ECHDC2 enhanced resistance to I/R injury. Furthermore, we observed that left anterior descending coronary artery ligation–induced myocardial infarction was more severe in the SS hearts than in the BN hearts or SSBN5 hearts, which was built on SS rats but had the substitution of chromosome 5 from BN rats. We also demonstrated that ECHDC2 did not alter mitochondrial O₂ consumption, metabolic intermediates and ATP production. By gas chromatography–mass spectrometry, we found that ECHDC2 overexpression increased the levels of the cellular branched chain amino acids leucine and valine.

Conclusion

ECHDC2, a mitochondrial protein, may be involved in regulating cell death and myocardial injury. Its deficiency in BN rats contributes to their increased resistance to myocardial I/R compared with SS rats. ECHDC2 increases branched chain amino acid metabolism and appears to be a novel regulator linking cell metabolism with cardiovascular disease.

Physiological genomics analysis for Alzheimer's disease

Alzheimer's disease is a common kind of dementia. This disorder can be detected in all countries around the world. This neurological disorder affects millions of population and becomes an important concern in modern neurology. There are many researches on the pathogenesis of Alzheimer's disease. Although it has been determined for a long time, there is no clear-cut that this is a case with genetic disorder or not. A physiological genomics is a new application that is useful for track function to genes within the human genome and can be applied for answering the problem of underlying pathobiology of complex diseases. The physiogenomics can be helpful for study of systemic approach on the pathophysiology, and genomics might provide useful information to better understand the pathogenesis of Alzheimer's disease. The present advent in genomics technique makes it possible to trace for the underlying genomics of disease. In this work, physiological genomics analysis for Alzheimer's disease was performed. The standard published technique is used for assessment. According to this work, there are 20 identified physiogenomics relationship on several chromosomes. Considering the results, the HADH2 gene on chromosome X, APBA1 gene on chromosome 9, AGER gene on chromosome 6, GSK3B gene on chromosome 3, CDKHR1 gene on chromosome 17, APPBP1 gene on chromosome 16, APBA2 gene on chromosome 15, GAL gene on chromosome 11, and APLP2 gene on chromosome 11 have the highest physiogenomics score (9.26) while the CASP3 gene on chromosome 4 and the SNCA gene on chromosome 4 have the lowest physiogenomics score (7.44). The results from this study confirm that Alzheimer's disease has a polygenomic origin.

Neurogenetics of aggressive behavior: studies in rodents

Aggressive behavior is observed in many animal species, such as insects, fish, lizards, frogs, and most mammals including humans. This wide range of conservation underscores the importance of aggressive behavior in the animals' survival and fitness, and the likely heritability of this behavior. Although typical patterns of aggressive behavior differ between species, there are several concordances in the neurobiology of aggression among rodents, primates, and humans. Studies with rodent models may eventually help us to understand the neurogenetic architecture of aggression in humans. However, it is important to recognize the difference between the ecological and ethological significance of aggressive behavior (species-typical aggression) and maladaptive violence (escalated aggression) when applying the findings of aggression research using animal models to human or veterinary medicine. Well-studied rodent models for aggressive behavior in the laboratory setting include the mouse (Mus musculus), rat (Rattus norvegicus), hamster (Mesocricetus auratus), and prairie vole (Microtus ochrogaster). The neural circuits of rodent aggression have been gradually elucidated by several techniques, e.g., immunohistochemistry of immediate-early gene (c-Fos) expression, intracranial drug microinjection, in vivo microdialysis, and optogenetics techniques. Also, evidence accumulated from the analysis of gene-knockout mice shows the involvement of several genes in aggression. Here, we review the brain circuits that have been implicated in aggression, such as the hypothalamus, prefrontal cortex (PFC), dorsal raphe nucleus (DRN), nucleus accumbens (NAc), and olfactory system. We then discuss the roles of glutamate and γ-aminobutyric acid (GABA), excitatory and inhibitory amino acids in the brain, as well as their receptors, in controlling aggressive behavior, focusing mainly on recent findings. At the end of this chapter, we discuss how genes can be identified that underlie individual differences in aggression, using the so-called forward genetics approach.

Increased proliferative cells in the medullary thick ascending limb of the loop of Henle in the Dahl salt-sensitive rat

Studies of transcriptome profiles have provided new insights into mechanisms underlying the development of hypertension. Cell type heterogeneity in tissue samples, however, has been a significant hindrance in these studies. We performed a transcriptome analysis in medullary thick ascending limbs of the loop of Henle isolated from Dahl salt-sensitive rats. Genes differentially expressed between Dahl salt-sensitive rats and salt-insensitive consomic SS.13(BN) rats on either 0.4% or 7 days of 8.0% NaCl diet (n=4) were highly enriched for genes located on chromosome 13, the chromosome substituted in the SS.13(BN) rat. A pathway involving cell proliferation and cell cycle regulation was identified as one of the most highly ranked pathways based on differentially expressed genes and by a Bayesian model analysis. Immunofluorescent analysis indicated that just 1 week of a high-salt diet resulted in a severalfold increase in proliferative medullary thick ascending limb cells in both rat strains, and that Dahl salt-sensitive rats exhibited a significantly greater proportion of medullary thick ascending limb cells in a proliferative state than in SS.13(BN) rats (15.0±1.4% versus 10.1±0.6%; n=7-9; P<0.05). The total number of cells per medullary thick ascending limb section analyzed was not different between the 2 strains. The study revealed alterations in regulatory pathways in Dahl salt-sensitive rats in tissues highly enriched for a single cell type, leading to the unexpected finding of a greater increase in the number of proliferative medullary thick ascending limb cells in Dahl salt-sensitive rats on a high-salt diet.

Chromosome substitution strains: gene discovery, functional analysis, and systems studies

Laboratory mice are valuable in biomedical research in part because of the extraordinary diversity of genetic resources that are available for studies of complex genetic traits and as models for human biology and disease. Chromosome substitution strains (CSSs) are important in this resource portfolio because of their demonstrated use for gene discovery, genetic and epigenetic studies, functional characterizations, and systems analysis. CSSs are made by replacing a single chromosome in a host strain with the corresponding chromosome from a donor strain. A complete CSS panel involves a total of 22 engineered inbred strains, one for each of the 19 autosomes, one each for the X and Y chromosomes, and one for mitochondria. A genome survey simply involves comparing each phenotype for each of the CSSs with the phenotypes of the host strain. The CSS panels that are available for laboratory mice have been used to dissect a remarkable variety of phenotypes and to characterize an impressive array of disease models. These surveys have revealed considerable phenotypic diversity even among closely related progenitor strains, evidence for strong epistasis and for heritable epigenetic changes. Perhaps most importantly, and presumably because of their unique genetic constitution, CSSs, and congenic strains derived from them, the genetic variants underlying quantitative trait loci (QTLs) are readily identified and functionally characterized. Together these studies show that CSSs are important resource for laboratory mice.

Inflammatory priming predisposes mice to age-related retinal degeneration

Disruption of cellular processes affected by multiple genes and accumulation of numerous insults throughout life dictate the progression of age-related disorders, but their complex etiology is poorly understood. Postmitotic neurons, such as photoreceptor cells in the retina and epithelial cells in the adjacent retinal pigmented epithelium, are especially susceptible to cellular senescence, which contributes to age-related retinal degeneration (ARD). The multigenic and complex etiology of ARD in humans is reflected by the relative paucity of effective compounds for its early prevention and treatment. To understand the genetic differences that drive ARD pathogenesis, we studied A/J mice, which develop ARD more pronounced than that in other inbred mouse models. Although our investigation of consomic strains failed to identify a chromosome associated with the observed retinal deterioration, pathway analysis of RNA-Seq data from young mice prior to retinal pathological changes revealed that increased vulnerability to ARD in A/J mice was due to initially high levels of inflammatory factors and low levels of homeostatic neuroprotective factors. The genetic signatures of an uncompensated preinflammatory state and ARD progression identified here aid in understanding the susceptible genetic loci that underlie pathogenic mechanisms of age-associated disorders, including several human blinding diseases.

Genetic divergence and the genetic architecture of complex traits in chromosome substitution strains of mice

Background

The genetic architecture of complex traits strongly influences the consequences of inherited mutations, genetic engineering, environmental and genetic perturbations, and natural and artificial selection. But because most studies are under-powered, the picture of complex traits is often incomplete. Chromosome substitution strains (CSSs) are a unique paradigm for these genome surveys because they enable statistically independent, powerful tests for the phenotypic effects of each chromosome on a uniform inbred genetic background. A previous CSS survey in mice and rats revealed many complex trait genes (QTLs), large phenotypic effects, extensive epistasis, as well as systems properties such as strongly directional phenotypic changes and genetically-determined limits on the range of phenotypic variation. However, the unusually close genetic relation between the CSS progenitor strains in that study raised questions about the impact of genetic divergence: would greater divergence between progenitor strains, with the corresponding changes in gene regulation and protein function, lead to significantly more distinctive phenotypic features, or alternatively would epistasis and systems constraints, which are pervasive in CSSs, limit the range of phenotypic variation regardless of the extent of DNA sequence variation?

Results

We analyzed results for an extensive survey of traits in two new panels of CSSs where the donor strains were derived from inbred strains with more distant origins and discovered a strong similarity in genetic and systems properties among the three CSS panels, regardless of divergence time.

Conclusion

Our results argue that DNA sequence differences between host and donor strains did not substantially affect the architecture of complex traits, and suggest instead that strong epistasis buffered the phenotypic effects of genetic divergence, thereby constraining the range of phenotypic variation.

Physiogenomic analysis of statin-treated patients: domain-specific counter effects within the ACACB gene on low-density lipoprotein cholesterol?

AIM

Administered at maximal dosages, the most common statins--atorvastatin, simvastatin and rosuvastatin--lower low-density lipoprotein cholesterol (LDLC) by an average of 37-57% in patients with primary hypercholesterolemia. We hypothesized novel genetic underpinnings for variation in LDLC levels in the context of statin therapy.

MATERIALS & METHODS

Genotyping of 384 SNPs in 202 volunteers from a lipid outpatient clinic was accomplished and LDLC levels obtained from chart records. The SNPs were distributed across 222 genes representing physiological pathways such as general metabolism, cholesterol biochemistry, cardiovascular function, inflammation, neurobiology and cell proliferation. We discovered significant associations with LDLC levels for the rs34274 SNP (p < 0.0002) and for rs2241220 (p < 0.008) in the acetyl-coenzyme A carboxylase beta (ACACB) gene. When corrected for multiple testing, the false-discovery rate associated with rs34274 was 0.076 (significance threshold: 0.10) and for rs2241220 the false-discovery rate was 0.93 (not significant). The acetyl coenzyme A carboxylase beta enzyme synthesizes malonyl coenzyme A, an essential substrate for hepatic fatty acid synthesis and an inhibitor of fatty acid oxidation.

RESULTS

The SNPs were in weak linkage disequilibrium (D = 0.302). Minor alleles at these sites demonstrate opposing influences on LDLC; the C>T substitution at rs34724 is a risk marker and the C>T substitution at rs2241220 a protective marker for LDLC levels. These SNPs hypothetically influence enzymatic activity through different mechanisms, rs34274 through the PII promoter and rs2241220 via alteration of the protein's responsiveness to allosteric influence.

CONCLUSION

Physiogenomic evidence suggests a novel link between LDLC levels and the regulation of fatty acid metabolism. The findings complement previously discovered novel SNP relationships to myalgia (pain) and myositis (serum creatine kinase activity). By genotyping for myositis, myalgia and LDLC levels, a physiogenomic model may be developed to help clinicians maximize effectiveness and minimize side effects in prescribing statins.

Genetic resistance to diet-induced obesity in chromosome substitution strains of mice

Discovery of genes that confer resistance to diseases such as diet-induced obesity could have tremendous therapeutic impact. We previously demonstrated that the C57BL/6J-Chr(A/J)/NaJ panel of chromosome substitution strains (CSSs) is a unique model for studying resistance to diet-induced obesity. In the present study, three replicate CSS surveys showed remarkable consistency, with 13 A/J-derived chromosomes reproducibly conferring resistance to high-fat-diet-induced obesity. Twenty CSS intercrosses, one derived from each of the 19 autosomes and chromosome X, were used to determine the number and location of quantitative trait loci (QTLs) on individual chromosomes and localized six QTLs. However, analyses of mean body weight in intercross progeny versus C57BL/6J provided strong evidence that many QTLs discovered in the CSS surveys eluded detection in these CSS intercrosses. Studies of the temporal effects of these QTLs suggest that obesity resistance was dynamic, with QTLs acting at different ages or after different durations of diet exposure. Thus, these studies provide insight into the genetic architecture of complex traits such as resistance to diet-induced obesity in the C57BL/6J-Chr(A/J)/NaJ CSSs. Because some of the QTLs detected in the CSS intercrosses were not detected using a traditional C57BL/6J x A/J intercross, our results demonstrate that surveys of CSSs and congenic strains derived from them are useful complementary tools for analyzing complex traits.

Receptor tyrosine kinase inhibitors in rodent pulmonary hypertension

Pulmonary hypertension (PH) is a disorder characterized by vascular remodeling and proliferation, a phenotype dependent upon unimpeded growth factor and kinase pathway activation with strong similarities to malignant tumors. This chapter details our novel application of the multikinase inhibitor, sorafenib, in rodent models of PH to improved hemodynamic parameters and attenuates PH structural changes1. Sorafenib is a Raf kinase inhibitor and our biochemical and genomic evidence supported the potential involvement of the MAPK cascade system and TGFB3 in PH development and the response to therapy. Integration of expression genomic analyses coupled with intense bioinformatics identified gene expression and ontology signatures in the development of PH and implicated the role of cytoskeletal protein such as caldesmon or nmMLCK as potentially key participants in PH-induced vascular remodeling and proliferation. Our studies suggest the PKI sorafenib as a potentially novel treatment for severe PH with the MAPK cascade a potential canonical target profoundly effecting vascular cytoskeletal -rearrangements and remodeling1.

Novel strategies and targets for the management of hypertension

Hypertension, as the sole or comorbid component of a constellation of disorders of the cardiovascular (CV) system, is present in over 90% of all patients with CV disease and affects nearly 74 million individuals in the United States. The number of medications available to treat hypertension has dramatically increased during the past 3 decades to some 50 medications as new targets involved in the normal regulation of blood pressure have been identified, resulting in the development of new agents in those classes with improved therapeutic profiles (e.g., renin-angiotensin-aldosterone system; RAAS). Despite these new agents, hypertension is not adequately managed in approximately 30% of patients, who are compliant with prescriptive therapeutics, suggesting that new agents and/or strategies to manage hypertension are still needed. Some of the newest classes of agents have targeted other components of the RAS, for example, the selective renin inhibitors, but recent advances in vascular biology have provided novel potential targets that may provide avenues for new agent development. These newer targets include downstream signaling participants in pathways involved in contraction, growth, hypertrophy, and relaxation. However, perhaps the most unique approach to the management of hypertension is a shift in strategy of using existing agents with respect to the time of day at which the agent is taken. This new strategy, termed "chronotherapy," has shown considerable promise in effectively managing hypertensive patients. Therefore, there remains great potential for future development of safe and effective agents and strategies to manage a disorder of the CV system of epidemic proportion.

Central and baroreflex control of heart period during the wake-sleep cycle in consomic rats with different genetic susceptibility to hypertension

1. In spontaneously hypertensive rats (SHR), the contributions of the baroreflex and central autonomic commands to the control of heart period (HP) vary among wake-sleep states and are impaired during quiet wakefulness and rapid eye movement sleep (REMS), respectively. 2. Dahl salt-sensitive (SS) rats are genetically susceptible to salt-sensitive hypertension, the development of which depends on diet. Substitution of chromosome 13 of SS rats with that of Brown Norway rats confers salt-resistance to consomic SS-13BN rats. 3. In the present study, we tested whether differences in the central and baroreflex contributions to HP control occur among wake-sleep states in SS and SS-13BN rats and reflect genetic susceptibility to hypertension. Rats (n = 5 per group) were fed a prohypertensive diet late during development to minimize hypertension in SS rats and were instrumented with an arterial catheter and electrodes for discriminating wake-sleep states. 4. The cross-correlation function between HP and blood pressure indicated that, in SS and SS-13BN rats, the contributions of the baroreflex and central commands to the control of HP differed significantly among wake-sleep states, with central commands outweighing the baroreflex in REMS. However, these contributions did not differ significantly between SS and SS-13BN rats in any wake-sleep state. 5. The data suggest that differences in the central and baroreflex contributions to HP control among wake-sleep states, which have been demonstrated in SHR, can be generalized to other rat models used in hypertension research. Impairments in the baroreflex and central autonomic control of HP during quiet wakefulness and REMS, respectively, cannot be generalized as an index of genetic susceptibility to hypertension.

Identification of candidate genes and gene networks specifically associated with analgesic tolerance to morphine

Chronic morphine administration may alter the expression of hundreds to thousands of genes. However, only a subset of these genes is likely involved in analgesic tolerance. In this report, we used a behavior genetics strategy to identify candidate genes specifically linked to the development of morphine tolerance. Two inbred genotypes [C57BL/6J (B6), DBA2/J (D2)] and two reciprocal congenic genotypes (B6D2, D2B6) with the proximal region of chromosome 10 (Chr10) introgressed into opposing backgrounds served as the behavior genetic filter. Tolerance after therapeutically relevant doses of morphine developed most rapidly in the B6 followed by the B6D2 genotype and did not develop in the D2 mice and only slightly in the D2B6 animals indicating a strong influence of the proximal region of Chr10 in the development of tolerance. Gene expression profiling and pattern matching identified 64, 53, 86, and 123 predisposition genes and 81, 96, 106, and 82 tolerance genes in the periaqueductal gray (PAG), prefrontal cortex, temporal lobe, and ventral striatum, respectively. A potential gene network was identified in the PAG in which 19 of the 34 genes were strongly associated with tolerance. Eleven of the network genes were found to reside in quantitative trait loci previously associated with morphine-related behaviors, whereas seven were predictive of tolerance (morphine-naive condition). Overall, the genes modified by chronic morphine administration show a strong presence in canonical pathways representative of neuroadaptation. A potentially significant role for the micro-RNA and epigenetic mechanisms in response to chronic administration of pharmacologically relevant doses of morphine was highlighted by candidate genes Dicer and H19.

Genetic architecture of complex traits: large phenotypic effects and pervasive epistasis

The genetic architecture of complex traits underlying physiology and disease in most organisms remains elusive. We still know little about the number of genes that underlie these traits, the magnitude of their effects, or the extent to which they interact. Chromosome substitution strains (CSSs) enable statistically powerful studies based on testing engineered inbred strains that have single, unique, and nonoverlapping genetic differences, thereby providing measures of phenotypic effects that are attributable to individual chromosomes. Here, we report a study of phenotypic effects and gene interactions for 90 blood, bone, and metabolic traits in a mouse CSS panel and 54 traits in a rat CSS panel. Two key observations emerge about the genetic architecture of these traits. First, the traits tend to be highly polygenic: across the genome, many individual chromosome substitutions each had significant phenotypic effects and, within each of the chromosomes studied, multiple distinct loci were found. Second, strong epistasis was found among the individual chromosomes. Specifically, individual chromosome substitutions often conferred surprisingly large effects (often a substantial fraction of the entire phenotypic difference between the parental strains), with the result that the sum of these individual effects often dramatically exceeded the difference between the parental strains. We suggest that strong, pervasive epistasis may reflect the presence of several phenotypically-buffered physiological states. These results have implications for identification of complex trait genes, developmental and physiological studies of phenotypic variation, and opportunities to engineer phenotypic outcomes in complex biological systems.

Genetic mapping of mammary tumor traits to rat chromosome 10 using a novel panel of consomic rats

Identification of novel breast cancer susceptibility and resistance genes in genetically diverse human populations is challenging, and so inbred rats have been used to identify novel mammary cancer susceptibility quantitative trait loci (QTLs) with conventional mapping approaches. An alternative approach for QTL mapping is to use chromosome substitution (consomic) rat strains, which has the advantage of rapid generation of congenic from consomic animals. Using a novel rat strain pair, SS and BN, we identified rat mammary cancer QTLs in one of two consomic rat strains tested. Female rats of inbred parental (SS and BN) and two consomic (SS-10 BN and SS-12 BN) strains were treated with 7,12-dimethylbenz[a]anthracene orally. The phenotypes of tumor incidence, latency, and multiplicity were evaluated. SS rats were highly susceptible to mammary adenocarcinoma development, whereas BN rats were completely resistant. Statistical comparison of the phenotypes between the susceptible parental and the two consomic strains identified QTLs residing within chromosome 10 controlling mammary tumor latency and multiplicity. The study shows that SS-BN consomic rat strains can be used to map mammary tumor QTLs. This novel approach should accelerate positional cloning of mammary cancer susceptibility and resistant genes in the rat and the identification of homologous genes in humans.

Distribution and functional impact of DNA copy number variation in the rat

The abundance and dynamics of copy number variants (CNVs) in mammalian genomes poses new challenges in the identification of their impact on natural and disease phenotypes. We used computational and experimental methods to catalog CNVs in rat and found that they share important functional characteristics with those in human. In addition, 113 one-to-one orthologous genes overlap CNVs in both human and rat, 80 of which are implicated in human disease. CNVs are nonrandomly distributed throughout the genome. Chromosome 18 is a cold spot for CNVs as well as evolutionary rearrangements and segmental duplications, suggesting stringent selective mechanisms underlying CNV genesis or maintenance. By exploiting gene expression data available for rat recombinant inbred lines, we established the functional relationship of CNVs underlying 22 expression quantitative trait loci. These characteristics make the rat an excellent model for studying phenotypic effects of structural variation in relation to human complex traits and disease.

Molecular networks in Dahl salt-sensitive hypertension based on transcriptome analysis of a panel of consomic rats

The Dahl salt-sensitive (SS) rat is a widely used model of human salt-sensitive hypertension and renal injury. We studied the molecular networks that underlie the complex disease phenotypes in the SS model, using a design that involved two consomic rat strains that were protected from salt-induced hypertension and one that was not protected. Substitution of Brown Norway (BN) chromosome 13 or 18, but not 20, into the SS genome was found to significantly attenuate salt-induced hypertension and albuminuria. Gene expression profiles were examined in the kidneys of SS and consomic SS-13(BN), SS-18(BN), and SS-20(BN) rats with a total of 240 cDNA microarrays. The substituted chromosome was overrepresented in genes differentially expressed between a consomic strain and SS rats on a 0.4% salt diet. F5, Serpinc1, Slc19a2, and genes represented by three other expressed sequence tags (ESTs), which are located on chromosome 13, were found to be differentially expressed between SS-13(BN) and all other strains examined. Likewise, Acaa2, B4galt6, Colec12, Hsd17b4, and five other ESTs located on chromosome 18 exhibited expression patterns unique to SS-18(BN). On exposure to a 4% salt diet, there were 184 ESTs in the renal cortex and 346 in the renal medulla for which SS-13(BN) and SS-18(BN) shared one expression pattern, while SS and SS-20(BN) shared another, mirroring the phenotypic segregation among the four strains. Molecular networks that might contribute to the development of Dahl salt-sensitive hypertension and albuminuria were constructed with an approach that merged biological knowledge-driven analysis and data-driven Bayesian probabilistic analysis.

Genomic assessment of a multikinase inhibitor, sorafenib, in a rodent model of pulmonary hypertension

Pulmonary hypertension (PH) and cancer pathology share growth factor- and MAPK stress-mediated signaling pathways resulting in endothelial and smooth muscle cell dysfunction and angioproliferative vasculopathy. In this study, we assessed sorafenib, an antineoplastic agent and inhibitor of multiple kinases important in angiogenesis [VEGF receptor (VEGFR)-1-3, PDGF receptor (PDGFR)-beta, Raf-1 kinase] as a potential PH therapy. Two PH rat models were used: a conventional hypoxia-induced PH model and an augmented PH model combining dual VEGFR-1 and -2 inhibition (SU-5416, single 20 mg/kg injection) with hypoxia. In addition to normoxia-exposed control animals, four groups were maintained at 10% inspired O(2) fraction for 3.5 wk (hypoxia/vehicle, hypoxia/SU-5416, hypoxia/sorafenib, and hypoxia/SU-5416/sorafenib). Compared with normoxic control animals, rats exposed to hypoxia/SU-5416 developed hemodynamic and histological evidence of severe PH while rats exposed to hypoxia alone displayed only mild elevations in hemodynamic values (pulmonary vascular and right ventricular pressures). Sorafenib treatment (daily gavage, 2.5 mg/kg) prevented hemodynamic changes and demonstrated dramatic attenuation of PH-associated vascular remodeling. Compared with normoxic control rats, expression profiling (Affymetrix platform) of lung RNA obtained from hypoxia [false discovery rate (FDR) 6.5%]- and hypoxia/SU-5416 (FDR 1.6%)-challenged rats yielded 1,019 and 465 differentially regulated genes (fold change >1.4), respectively. A novel molecular signature consisting of 38 differentially expressed genes between hypoxia/SU-5416 and hypoxia/SU-5416/sorafenib (FDR 6.7%) was validated by either real-time RT-PCR or immunoblotting. Finally, immunoblotting studies confirmed the upregulation of the MAPK cascade in both PH models, which was abolished by sorafenib. In summary, sorafenib represents a novel potential treatment for severe PH with the MAPK cascade a potential canonical target.

Integrating genomic and clinical medicine: searching for susceptibility genes in complex lung diseases

The integration of molecular, genomic, and clinical medicine in the post-genome era provides the promise of novel information on genetic variation and pathophysiologic cascades. The current challenge is to translate these discoveries rapidly into viable biomarkers that identify susceptible populations and into the development of precisely targeted therapies. In this article, we describe the application of comparative genomics, microarray platforms, genetic epidemiology, statistical genetics, and bioinformatic approaches within examples of complex pulmonary pathobiology. Our search for candidate genes, which are gene variations that drive susceptibility to and severity of enigmatic acute and chronic lung disorders, provides a logical framework to understand better the evolution of genomic medicine. The dissection of the genetic basis of complex diseases and the development of highly individualized therapies remain lofty but achievable goals.

Physiogenomic analysis of localized FMRI brain activity in schizophrenia

The search for genetic factors associated with disease is complicated by the complexity of the biological pathways linking genotype and phenotype. This analytical complexity is particularly concerning in diseases historically lacking reliable diagnostic biological markers, such as schizophrenia and other mental disorders. We investigate the use of functional magnetic resonance imaging (fMRI) as an intermediate phenotype (endophenotype) to identify physiogenomic associations to schizophrenia. We screened 99 subjects, 30 subjects diagnosed with schizophrenia, 13 unaffected relatives of schizophrenia patients, and 56 unrelated controls, for gene polymorphisms associated with fMRI activation patterns at two locations in temporal and frontal lobes previously implied in schizophrenia. A total of 22 single nucleotide polymorphisms (SNPs) in 15 genes from the dopamine and serotonin neurotransmission pathways were genotyped in all subjects. We identified three SNPs in genes that are significantly associated with fMRI activity. SNPs of the dopamine beta-hydroxylase (DBH) gene and of the dopamine receptor D4 (DRD4) were associated with activity in the temporal and frontal lobes, respectively. One SNP of serotonin-3A receptor (HTR3A) was associated with temporal lobe activity. The results of this study support the physiogenomic analysis of neuroimaging data to discover associations between genotype and disease-related phenotypes.

Closely linked non-additive blood pressure quantitative trait loci

There is enough evidence through linkage and substitution mapping to indicate that rat chromosome 1 harbors multiple blood pressure (BP) quantitative trait loci (QTLs). Of these, BP QTL1b was previously reported from our laboratory using congenic strains derived by introgressing normotensive alleles from the LEW rat onto the genetic background of the hypertensive Dahl salt-sensitive (S) rat. The region spanned by QTL1b is quite large (20.92 Mb), thus requiring further mapping with improved resolution so as to facilitate systematic identification of the underlying genetic determinant(s). Using congenic strains containing the LEW rat chromosomal segments on the Dahl salt-sensitive (S) rat background, further iterations of congenic substrains were constructed and characterized. Collective data obtained from this new iteration of congenic substrains provided evidence for further fragmentation of QTL1b with improved resolution. At least two separate genetic determinants of blood pressure underlie QTL1b. These are within 7.40 Mb and 7.31 Mb and are known as the QTL1b1 region and the QTL1b2 region, respectively. A genetic interaction was detected between the two BP QTLs. Interestingly, five of the previously reported differentially expressed genes located within the newly mapped QTL1b1 region remained differentially expressed. The congenic strain S.LEW(D1Mco36-D1Mco101), which harbors the QTL1b1 region alone but not the QTL1b2 region, serves as a genetic tool for further dissection of the QTL1b1 region and validation of Nr2f2 as a positional candidate gene. Overall, this study represents an intermediary yet obligatory progression towards the identification of genetic elements controlling BP.

Repeatability of published microarray gene expression analyses

Given the complexity of microarray-based gene expression studies, guidelines encourage transparent design and public data availability. Several journals require public data deposition and several public databases exist. However, not all data are publicly available, and even when available, it is unknown whether the published results are reproducible by independent scientists. Here we evaluated the replication of data analyses in 18 articles on microarray-based gene expression profiling published in Nature Genetics in 2005-2006. One table or figure from each article was independently evaluated by two teams of analysts. We reproduced two analyses in principle and six partially or with some discrepancies; ten could not be reproduced. The main reason for failure to reproduce was data unavailability, and discrepancies were mostly due to incomplete data annotation or specification of data processing and analysis. Repeatability of published microarray studies is apparently limited. More strict publication rules enforcing public data availability and explicit description of data processing and analysis should be considered.

Comparative genomics for detecting human disease genes

Originally, comparative genomics was geared toward defining the synteny of genes between species. As the human genome project accelerated, there was an increase in the number of tools and means to make comparisons culminating in having the genomic sequence for a large number of organisms spanning the evolutionary tree. With this level of resolution and a long history of comparative biology and comparative genetics, it is now possible to use comparative genomics to build or select better animal models and to facilitate gene discovery. Comparative genomics takes advantage of the functional genetic information from other organisms, (vertebrates and invertebrates), to apply it to the study of human physiology and disease. It allows for the identification of genes and regulatory regions, and for acquiring knowledge about gene function. In this chapter, the current state of comparative genomics and the available tools are discussed in the context of developing animal model systems that reflect the clinical picture.

Dissection of a genetic locus influencing renal function in the rat and its concordance with kidney disease loci on human chromosome 1q21

Previously, we conducted a genome scan on a population derived from the Dahl salt-sensitive hypertensive (S) and the spontaneously hypertensive rat (SHR) using urinary albumin excretion (UAE) as our primary measure of renal function. We identified 10 quantitative trait loci (QTL) linked to several renal and/or cardiovascular traits. In particular, linkage and subsequent congenic strain analysis demonstrated that the loci on chromosome 2 had a large and significant effect on UAE compared with the S rat. The present work sought to characterize the chromosome 2 congenic strain [S.SHR] by conducting a time-course analysis (week 4-20), including evaluating additional renal parameters, histology, electron microscopy, and gene expression/ pathway analysis. Throughout the time course the congenic strain consistently maintained a threefold reduction in UAE compared with S rats and was supported by the histological findings of significantly reduced glomerular, tubular and interstitial changes. Gene expression/pathway analysis performed at week 4, 12, and 20 revealed that pathways involved in cellular assembly and organization, cellular movement, and immune response were controlled differently between the S and congenic. When all the data are considered, the chromosome 2 congenic appears to attenuate renal damage primarily through an altered fibrotic response. Recombinant progeny testing was employed to reduce the QTL to approximately 1.5 cM containing several interesting candidate genes. The concordance of this rat QTL with renal disease loci on human chromosome 1q21 demonstrate that elucidating the causative gene and mechanism of the rat QTL may be of particular importance for understanding kidney disease in humans.

Physiogenomic strategies and resources to associate genes with rat models of heart, lung and blood disorders

As is the case for many human disorders, cardiovascular disease is a complex ailment exhibiting a multifactorial mode of transmission. Rat models have been developed to aid in the analysis of this complex genetic and phenotypic disorder. The purpose of this brief review is to describe current gene expression profiling strategies that have been implemented to search for candidate causative genes of disease phenotypes in animal models. Strategies include integrating gene expression information with linkage analysis, expression profiling chromosome-substituted and/or congenic rat strains, correlating gene expression with physiological data across a panel of rodent strains, and linking expression quantitative trait loci to physiological quantitative trait loci. A primary goal of these strategies is to narrow and prioritize the search for causal genes of physiological interest. Also discussed are ways to harness two recent publicly available resources that have been created to investigate the role of genes and environment on cardiovascular physiology and pathophysiology.

Peroxisome proliferator-activated receptors as molecular targets in relation to obesity and Type 2 diabetes

The three isotypes of peroxisome proliferator-activated receptors (PPARs) are currently perceived as major regulatory nodes (or hubs) of metabolic pathway networks, linking most prevalent diseases including Type 2 diabetes, obesity, dyslipidemia and atherosclerosis. The integrative functions of PPARs are also reflected in their ecogenetic profile, when the variants underlying pharmacogenetic interactions were also shown to modulate the effect of lifestyle factors. Despite their extensive clinical use, there are many outstanding issues, especially concerning their safety. Critical pharmacogenomic assessment is warranted for the new potent ligands of multiple PPAR isoforms as many have displayed serious side-effects in a limited number of treated subjects. Nevertheless, the advent of genomic, transcriptomic and system biology-level approaches, integrating knowledge from model systems and human biology, should greatly facilitate the transition to individualized PPAR-based therapies.

Allograft vasculopathy versus atherosclerosis

Over the last 4 decades, heart transplantation (HTx) has evolved as a mainstream therapy for heart failure. Approximately half of patients needing HTx have organ failure consequent to atherosclerosis. Despite advances in immunosuppressive drugs, long-term success of HTx is limited by the development of a particular type of coronary atherosclerosis, referred to as cardiac allograft vasculopathy (CAV). Although the exact pathogenesis of CAV remains to be established, there is strong evidence that CAV involves immunologic mechanisms operating in a milieu of nonimmunologic risk factors. The immunologic events constitute the principal initiating stimuli, resulting in endothelial injury and dysfunction, altered endothelial permeability, with consequent myointimal hyperplasia and extracellular matrix synthesis. Lipid accumulation in allograft arteries is prominent, with lipoprotein entrapment in the subendothelial tissue, through interactions with proteoglycans. The apparent endothelial "intactness" in human coronary arteries of the transplanted heart suggest that permeability and function of the endothelial barrier altered. Various insults to the vascular bed result in vascular smooth muscle cell (SMC) activation. Activated SMCs migrate from the media into the intima, proliferate, and elaborate cytokines and extracellular matrix proteins, resulting in luminal narrowing and impaired vascular function. Arteriosclerosis is a broad term that is used to encompass all diseases that lead to arterial hardening, including native atherosclerosis, postangioplasty restenosis, vein bypass graft occlusion, and CAV. These diseases exhibit many similarities; however, they are distinct from one another in numerous ways as well. The present review summarizes the current understanding of the risk factors and the pathophysiological similarities and differences between CAV and atherosclerosis.

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.