Transcriptional profiling with a blood pressure QTL interval-specific oligonucleotide array

QTL mapping of rat blood pressure loci on RNO1 within a homologous region linked to human hypertension on HSA15

Fine-mapping of regions linked to the inheritance of hypertension is accomplished by genetic dissection of blood pressure quantitative trait loci (BP QTLs) in rats. The goal of the current study was to further fine-map two genomic regions on rat chromosome 1 with opposing blood pressure effects (BP QTL1b1 and BP QTL1b1a), the homologous region of which on human chromosome 15 harbors BP QTLs. Two new substrains were constructed and studied from the previously reported BP QTL1b1, one having significantly lower systolic BP by 17 mmHg than that of the salt-sensitive (S) rat (P = 0.007). The new limits of BP QTL1b1 were between 134.09 Mb and 135.40 Mb with a 43% improvement from the previous 2.31 Mb to the current 1.31 Mb interval containing 4 protein-coding genes (Rgma, Chd2, Fam174b, and St8sia2), 2 predicted miRNAs, and 4 lncRNAs. One new substrain was constructed and studied from the previously reported BPQTL1b1a having a significantly higher systolic BP by 22 mmHg (P = 0.006) than that of the S rat. The new limits of BPQTL1b1a were between 133.53 Mb and 134.52 Mb with a 32% improvement from the previous1.45 Mb to the current 990.21 Kb interval containing 1 protein-coding gene, Mctp2, and a lncRNA. The congenic segments of these two BP QTLs overlapped between 134.09 Mb and 134.52 Mb. No exonic variants were detected in any of the genes. These findings reiterate complexity of genetic regulation of BP within QTL regions, where elements beyond protein-coding sequences could be factors in controlling BP.

COUP-TFII revisited: Its role in metabolic gene regulation

Chicken Ovalbumin Upstream Promoter Transcription Factor II (COUP-TFII) is an orphan member of the nuclear receptor family of transcriptional regulators. Although hormonal activation of COUP-TFII has not yet been identified, rodent genetic models have uncovered vital and diverse roles for COUP-TFII in biological processes. These include control of cardiac function and angiogenesis, reproduction, neuronal development, cell fate and organogenesis. Recently, an emerging body of evidence has demonstrated COUP-TFII involvement in various metabolic systems such as adipogenesis, lipid metabolism, hepatic gluconeogenesis, insulin secretion, and regulation of blood pressure. The potential relevance of these observations to human pathology has been corroborated by the identification of single nucleotide polymorphism in the human COUP-TFII promoter controlling insulin sensitivity. Of particular interest to metabolism is the ability of COUP-TFII to interact with the Glucocorticoid Receptor (GR). This interaction is known to control gluconeogenesis, principally through direct binding of COUP-TFII/GR complexes to the promoters of gluconeogenic enzyme genes. However, it is likely that this interaction is critical to other metabolic processes, since GR, like COUP-TFII, is an essential regulator of adipogenesis, insulin sensitivity, and blood pressure. This review will highlight these unique roles of COUP-TFII in metabolic gene regulation.

Sodium bicarbonate cotransporter NBCe2 gene variants increase sodium and bicarbonate transport in human renal proximal tubule cells

RATIONALE

Salt sensitivity of blood pressure affects >30% of the hypertensive and >15% of the normotensive population. Variants of the electrogenic sodium bicarbonate cotransporter NBCe2 gene, SLC4A5, are associated with increased blood pressure in several ethnic groups. SLC4A5 variants are also highly associated with salt sensitivity, independent of hypertension. However, little is known about how NBCe2 contributes to salt sensitivity, although NBCe2 regulates renal tubular sodium bicarbonate transport. We hypothesized that SLC4A5 rs10177833 and rs7571842 increase NBCe2 expression and human renal proximal tubule cell (hRPTC) sodium transport and may be a cause of salt sensitivity of blood pressure.

OBJECTIVE

To characterize the hRPTC ion transport of wild-type (WT) and homozygous variants (HV) of SLC4A5.

METHODS AND RESULTS

The expressions of NBCe2 mRNA and protein were not different between hRPTCs carrying WT or HV SLC4A5 before or after dopaminergic or angiotensin (II and III) stimulation. However, luminal to basolateral sodium transport, NHE3 protein, and Cl-/HCO3- exchanger activity in hRPTCs were higher in HV than WT SLC4A5. Increasing intracellular sodium enhanced the apical location of NBCe2 in HV hRPTCs (4.24±0.35% to 11.06±1.72% (P<0.05, N = 3, 2-way ANOVA, Holm-Sidak test)) as determined by Total Internal Reflection Fluorescence Microscopy (TIRFM). In hRPTCs isolated from kidney tissue, increasing intracellular sodium enhanced bicarbonate-dependent pH recovery rate and increased NBCe2 mRNA and protein expressions to a greater extent in HV than WT SLC4A5 (+38.00±6.23% vs HV normal salt (P<0.01, N = 4, 2-way ANOVA, Holm-Sidak test)). In hRPTCs isolated from freshly voided urine, bicarbonate-dependent pH recovery was also faster in those from salt-sensitive and carriers of HV SLC4A5 than from salt-resistant and carriers of WT SLC4A5. The faster NBCe2-specific bicarbonate-dependent pH recovery rate in HV SCL4A5 was normalized by SLC4A5- but not SLC4A4-shRNA. The binding of purified hepatocyte nuclear factor type 4A (HNF4A) to DNA was increased in hRPTCs carrying HV SLC4A5 rs7571842 but not rs10177833. The faster NBCe2-specific bicarbonate-dependent pH recovery rate in HV SCL4A5 was abolished by HNF4A antagonists.

CONCLUSION

NBCe2 activity is stimulated by an increase in intracellular sodium and is hyper-responsive in hRPTCs carrying HV SLC4A5 rs7571842 through an aberrant HNF4A-mediated mechanism.

Towards Precision Medicine for Hypertension: A Review of Genomic, Epigenomic, and Microbiomic Effects on Blood Pressure in Experimental Rat Models and Humans

Compelling evidence for the inherited nature of essential hypertension has led to extensive research in rats and humans. Rats have served as the primary model for research on the genetics of hypertension resulting in identification of genomic regions that are causally associated with hypertension. In more recent times, genome-wide studies in humans have also begun to improve our understanding of the inheritance of polygenic forms of hypertension. Based on the chronological progression of research into the genetics of hypertension as the "structural backbone," this review catalogs and discusses the rat and human genetic elements mapped and implicated in blood pressure regulation. Furthermore, the knowledge gained from these genetic studies that provide evidence to suggest that much of the genetic influence on hypertension residing within noncoding elements of our DNA and operating through pervasive epistasis or gene-gene interactions is highlighted. Lastly, perspectives on current thinking that the more complex "triad" of the genome, epigenome, and the microbiome operating to influence the inheritance of hypertension, is documented. Overall, the collective knowledge gained from rats and humans is disappointing in the sense that major hypertension-causing genes as targets for clinical management of essential hypertension may not be a clinical reality. On the other hand, the realization that the polygenic nature of hypertension prevents any single locus from being a relevant clinical target for all humans directs future studies on the genetics of hypertension towards an individualized genomic approach.

Microbiotal-Host Interactions and Hypertension

Hypertension, or elevated blood pressure (BP), has been extensively researched over decades and clearly demonstrated to be caused due to a combination of host genetic and environmental factors. Although much research remains to be conducted to pin-point the precise genetic elements on the host genome that control BP, new lines of evidence are emerging to indicate that, besides the host genome, the genomes of all indigenous commensal micro-organisms, collectively referred to as the microbial metagenome or microbiome, are important, but largely understudied, determinants of BP. Unlike the rigid host genome, the microbiome or the "second genome" can be altered by diet or microbiotal transplantation in the host. This possibility is attractive from the perspective of exploiting the microbiotal composition for clinical management of inherited hypertension. Thus, focusing on the limited current literature supporting a role for the microbiome in BP regulation, this review highlights the need to further explore the role of the co-existence of host and the microbiota as an organized biological unit called the "holobiont" in the context of BP regulation.

High-resolution mapping of a novel rat blood pressure locus on chromosome 9 to a region containing the Spp2 gene and colocalization of a QTL for bone mass

Through linkage analysis of the Dahl salt-sensitive (S) rat and the spontaneously hypertensive rat (SHR), a blood pressure (BP) quantitative trait locus (QTL) was previously located on rat chromosome 9. Subsequent substitution mapping studies of this QTL revealed multiple BP QTLs within the originally identified logarithm of odds plot by linkage analysis. The focus of this study was on a 14.39 Mb region, the distal portion of which remained unmapped in our previous studies. High-resolution substitution mapping for a BP QTL in the setting of a high-salt diet indicated that an SHR-derived congenic segment of 787.9 kb containing the gene secreted phosphoprotein-2 (Spp2) lowered BP and urinary protein excretion. A nonsynonymous G/T polymorphism in the Spp2 gene was detected between the S and S.SHR congenic rats. A survey of 45 strains showed that the T allele was rare, being detected only in some substrains of SHR and WKY. Protein modeling prediction through SWISSPROT indicated that the predicted protein product of this variant was significantly altered. Importantly, in addition to improved cardiovascular and renal function, high salt-fed congenic animals carrying the SHR T variant of Spp2 had significantly lower bone mass and altered bone microarchitecture. Total bone volume and volume of trabecular bone, cortical thickness, and degree of mineralization of cortical bone were all significantly reduced in congenic rats. Our study points to opposing effects of a congenic segment containing the prioritized candidate gene Spp2 on BP and bone mass.

Treg Cell Differentiation: From Thymus to Peripheral Tissue

Regulatory T cells (Tregs) are crucial mediators of self-tolerance in the periphery. They differentiate in the thymus, where interactions with thymus-resident antigen-presenting cells, an instructive cytokine milieu, and stimulation of the T cell receptor lead to the selection into the Treg lineage and the induction of Foxp3 gene expression. Once mature, Treg cells leave the thymus and migrate into either the secondary lymphoid tissues, e.g., lymph nodes and spleen, or peripheral nonlymphoid tissues. There is growing evidence that Treg cells go beyond the classical modulation of immune responses and also play important functional roles in nonlymphoid peripheral tissues. In this review, we summarize recent findings about the thymic Treg lineage differentiation as well as the further specialization of Treg cells in the secondary lymphoid and in the peripheral nonlymphoid organs.

Screening for genes involved in antibody response to sheep red blood cells in the chicken, Gallus gallus

Antibody response, an important trait in both agriculture and biomedicine, plays a part in protecting animals from infection. Dissecting molecular basis of antibody response may improve artificial selection for natural disease resistance in livestock and poultry. A number of genetic markers associated with antibody response have been identified in the chicken and mouse by linkage-based association studies, which only define genomic regions by genetic markers but do not pinpoint genes for antibody response. In contrast, global expression profiling has been applied to define the molecular bases of a variety of biological traits through identification of differentially expressed genes (DEGs). Here, we employed Affimetrix GeneChip Chicken Genome Arrays to identify differentially expressed genes for antibody response to sheep red blood cells (SRBC) using chickens challenged with and without SRBC or chickens with high and low anti-SRBC titers. The DEGs include those with known (i.e., MHC class I and IgH genes) or unknown function in antibody response. Classification test of these genes suggested that the response of the chicken to intravenous injection of SRBC involved multiple biological processes, including response to stress or other different stimuli, sugar, carbohydrate or protein binding, and cell or soluble fraction, in addition to antibody response. This preliminary study thus provides an insight into molecular basis of antibody response to SRBC in the chicken.

Dr Lewis Kitchener Dahl, the Dahl rats, and the "inconvenient truth" about the genetics of hypertension

Lewis K. Dahl is regarded as an iconic figure in the field of hypertension research. During the 1960s and 1970s he published several seminal articles in the field that shed light on the relationship between salt and hypertension. Further, the Dahl rat models of hypertension that he developed by a selective breeding strategy are among the most widely used models for hypertension research. To this day, genetic studies using this model are ongoing in our laboratory. While Dr. Dahl is known for his contributions to the field of hypertension, very little, if any, of his personal history is documented. This article details a short biography of Dr. Lewis Dahl, the history behind the development of the Dahl rats and presents an overview of the results obtained through the genetic analysis of the Dahl rat as an experimental model to study the inheritance of hypertension.

Mutation within the hinge region of the transcription factor Nr2f2 attenuates salt-sensitive hypertension

Genome-wide association studies (GWAS) have prioritized a transcription factor, nuclear receptor 2 family 2 (NR2F2), as being associated with essential hypertension in humans. Here we provide evidence that validates this association and indicates that Nr2f2 is a genetic determinant of blood pressure (BP). Using the zinc-finger nuclease technology, the generation of a targeted Nr2f2-edited rat model is reported. The resulting gene-edited rats have a 15 bp deletion in exon 2 leading to a five-amino-acid deletion in the hinge region of the mutant Nr2f2 protein. Both systolic and diastolic blood pressures of the Nr2f2(mutant) rats are significantly lower than controls. Because the hinge region of Nr2f2 is required for interaction with Friend of Gata2 (Fog2), protein-protein interaction is examined. Interaction of Nr2f2(mutant) protein with Fog2 is greater than that with the wild-type Nr2f2, indicating that the extent of interaction between these two transcription factors critically influences BP.

Multiple blood pressure loci with opposing blood pressure effects on rat chromosome 1 in a homologous region linked to hypertension on human chromosome 15

Genetic dissection of blood pressure (BP) quantitative trait loci (QTLs) in rats has facilitated the fine-mapping of regions linked to the inheritance of hypertension. The goal of the current study was to further fine-map one such genomic region on rat chromosome 1 (BPQTL1b1), the homologous region of which on human chromosome 15 harbors BP QTLs, as reported by four independent studies. Of the six substrains constructed and studied, the systolic BP of two of the congenic strains were significantly lower by 36 and 27 mm Hg than that of the salt-sensitive (S) rat (P < 0.0001, P = 0.0003, respectively). The congenic segments of these two strains overlapped between 135.12 and 138.78 Mb and contained eight genes and two predicted miRNAs. None of the annotations had variants within expressed sequences. These data taken together with the previous localization resolved QTL1b1 with a 70% improvement from the original 7.39 Mb to the current 2.247 Mb interval. Furthermore, the systolic BP of one of the congenic substrains was significantly higher by 20 mm Hg (P < 0.0001) than the BP of the S rat. The limits of this newly identified QTL with a BP increasing effect (QTL1b1a) were between 134.12 and 135.76 Mb, spanning 1.64 Mb, containing two protein-coding genes, Mctp2 and Rgma, and a predicted miRNA. There were four synonymous variants within Mctp2. These data provide evidence for two independent BP QTLs with opposing BP effects within the previously identified BP QTL1b1 region. Additionally, these findings illustrate the complexity underlying the genetic mechanisms of BP regulation, wherein inherited elements beyond protein-coding sequences or known regulatory regions could be operational.

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.

Regulation of renin expression by the orphan nuclear receptors Nr2f2 and Nr2f6

Understanding the transcriptional mechanisms of renin expression is key to understanding the regulation of the renin-angiotensin system. We previously identified the nuclear receptors RAR/RXR and Nr2f6 (EAR2) as positive and negative transcriptional regulators of renin expression, respectively (Liu X, Huang X, Sigmund CD. Circ Res 92: 1033-1040, 2003). Both mediate their effects through a hormone response element (HRE) within the renin enhancer. Here, we determined whether another nuclear receptor, Nr2f2 (Coup-TFII, Arp-1), identified in a screen of proteins that bind the HRE, also regulates renin expression. Luciferase assays indicate that Nr2f2 negatively regulates the renin promoter more potently than Nr2f6. Gel-shift and chromatin immunoprecipitation (ChIP) indicate that Nr2f2 and Nr2f6 can bind directly to the renin enhancer through the HRE. Surprisingly, baseline expression of endogenous renin was not effected when Nr2f2 was knocked down in As4.1 cells, whereas knockdown of Nr2f6 increased renin expression twofold. Interestingly, however, knockdown of Nr2f2 augmented the induction of renin expression caused by retinoic acid. These data indicate that both Nr2f6 and Nr2f2 can negatively regulate the renin promoter, under baseline conditions and in response to physiological queues, respectively. Therefore, Nr2f2 may require an initiating signal that results in a change at the chromatin level or activation of another transcription factor to exert its effects. We conclude that both Nr2f2 and Nr2f6 negatively regulate renin promoter activity, but may do so by divergent mechanisms.

QTL and gene expression analyses identify genes affecting carcass weight and marbling on BTA14 in Hanwoo (Korean Cattle)

Causal mutations affecting quantitative trait variation can be good targets for marker-assisted selection for carcass traits in beef cattle. In this study, linkage and linkage disequilibrium analysis (LDLA) for four carcass traits was undertaken using 19 markers on bovine chromosome 14. The LDLA analysis detected quantitative trait loci (QTL) for carcass weight (CWT) and eye muscle area (EMA) at the same position at around 50 cM and surrounded by the markers FABP4SNP2774C>G and FABP4_μsat3237. The QTL for marbling (MAR) was identified at the midpoint of markers BMS4513 and RM137 in a 3.5-cM marker interval. The most likely position for a second QTL for CWT was found at the midpoint of tenth marker bracket (FABP4SNP2774C>G and FABP4_μsat3237). For this marker bracket, the total number of haplotypes was 34 with a most common frequency of 0.118. Effects of haplotypes on CWT varied from a -5-kg deviation for haplotype 6 to +8 kg for haplotype 23. To determine which genes contribute to the QTL effect, gene expression analysis was performed in muscle for a wide range of phenotypes. The results demonstrate that two genes, LOC781182 (p = 0.002) and TRPS1 (p = 0.006) were upregulated with increasing CWT and EMA, whereas only LOC614744 (p = 0.04) has a significant effect on intramuscular fat (IMF) content. Two genetic markers detected in FABP4 were the most likely QTL position in this QTL study, but FABP4 did not show a significant effect on both traits (CWT and EMA) in gene expression analysis. We conclude that three genes could be potential causal genes affecting carcass traits CWT, EMA, and IMF in Hanwoo.

Vascular responses in aortic rings of a consomic rat panel derived from the Fawn Hooded Hypertensive strain

The present experiments, utilizing the high-throughput vascular protocol of PhysGen (Program for Genomic Applications) characterized the responses of aortic rings to vasoconstrictor (phenylephrine) and vasodilator (acetylcholine, sodium nitroprusside, and reduced tissue bath Po(2)) stimuli in consomic rat strains derived from a cross between the Fawn Hooded Hypertensive rat (FHH/EurMcwi) and the Brown Norway normotensive (BN/NHsdMcwi) rat. The effects of substituting individual BN chromosomes into the FHH genetic background were determined in animals that were maintained on a low-salt (0.4% NaCl) diet or switched to a high-salt (4% NaCl) diet for 3 wk. Sex-specific differences were evaluated in male and female consomic rats on similar dietary salt intake. Multiple chromosomes affected various vascular reactivity phenotypes in the FHH × BN consomic panel, and substantial salt-dependent changes in vascular reactivity and sex-specific differences in aortic reactivity were observed in individual consomic strains. However, compared with earlier studies of consomic rats derived from a cross between the BN rat and the Dahl salt-sensitive (SS) rat, only 3-7% of the vascular phenotypes were affected in a similar manner by substituting specific BN chromosomeschromosomes into the FHH genetic background versus the SS genetic background. The findings of the present study stress the potential value of consomic rat panels in gaining insight into genetic factors influencing vascular reactivity and suggest that the chromosomes that appear to be involved in the determination of aortic ring reactivity in different rodent models of hypertension are highly strain- and sex specific.

Genome-wide searching of rare genetic variants in WTCCC data

Although they have demonstrated success in searching for common variants for complex diseases, genome-wide association (GWA) studies are less successful in detecting rare genetic variants because of the poor statistical power of most of current methods. We developed a two-stage method that can apply to GWA studies for detecting rare variants. Here we report the results of applying this two-stage method to the Wellcome Trust Case Control Consortium (WTCCC) dataset that include seven complex diseases: bipolar disorder, cardiovascular disease, hypertension (HT), rheumatoid arthritis, Crohn's disease, type 1 diabetes and type 2 diabetes (T2D). We identified 24 genes or regions that reach genome wide significance. Eight of them are novel and were not reported in the WTCCC study. The cumulative risk (or protective) haplotype frequency for each of the 8 genes or regions is small, being at most 11%. For each of the novel genes, the risk (or protective) haplotype set cannot be tagged by the common SNPs available in chips (r (2) < 0.32). The gene identified in HT was further replicated in the Framingham Heart Study, and is also significantly associated with T2D. Our analysis suggests that searching for rare genetic variants is feasible in current GWA studies and candidate gene studies, and the results can severe as guides to future resequencing studies to identify the underlying rare functional variants.

Persistent gene expression changes in ventral tegmental area of adolescent but not adult rats in response to chronic nicotine

Because adolescent brains are undergoing extensive developmental changes, they may be uniquely sensitive to effects of addictive drugs like nicotine. We exposed adolescent and adult rats to nicotine infusion for two weeks, and then used whole genome microarray analysis to determine effects on gene expression in the ventral tegmental area. We examined brains immediately after two weeks of nicotine or saline, and also four weeks after termination of nicotine exposure. After identifying genes with a significant agextreatment interaction, we employed template matching to find specific patterns of expression across age and treatment. Of those genes that were transiently regulated (up- or down-regulated immediately following the end of nicotine treatment, but back to saline baseline 30 days later), two-thirds were specific to adult animals, while only 30% were specific to adolescents and 4% were shared across the two ages. In contrast, significant genes that were persistently regulated (altered following nicotine treatment and still altered 30 days later) were more likely (59%) to be adolescent, with only 32% in adults and 8% shared. The greatest number of significant genes was late-regulated (no change immediately after nicotine, but regulated 30 days later). Again, most were in adolescents (54%), compared to adults (10%) or shared (36%). Pathway analysis revealed that adolescent-specific genes were over-represented in several biological functions and canonical pathways, including nervous system development and function and long-term potentiation. Furthermore, adolescent-specific genes formed extensive interaction networks, unlike those specific for adults or shared. This age-specific expression pattern may relate to the heightened vulnerability of adolescents to the effects of addictive drugs. In particular, the propensity of adolescents to show persistent alterations in gene expression corresponds to the persistence of drug dependence among smokers who began their habit as adolescents. These findings support a model whereby adolescent brains are uniquely vulnerable to long-term changes in gene expression in the brain's reward pathway caused by early exposure to nicotine.

Serious limitations of the QTL/microarray approach for QTL gene discovery

BACKGROUND

It has been proposed that the use of gene expression microarrays in nonrecombinant parental or congenic strains can accelerate the process of isolating individual genes underlying quantitative trait loci (QTL). However, the effectiveness of this approach has not been assessed.

RESULTS

Thirty-seven studies that have implemented the QTL/microarray approach in rodents were reviewed. About 30% of studies showed enrichment for QTL candidates, mostly in comparisons between congenic and background strains. Three studies led to the identification of an underlying QTL gene. To complement the literature results, a microarray experiment was performed using three mouse congenic strains isolating the effects of at least 25 biometric QTL. Results show that genes in the congenic donor regions were preferentially selected. However, within donor regions, the distribution of differentially expressed genes was homogeneous once gene density was accounted for. Genes within identical-by-descent (IBD) regions were less likely to be differentially expressed in chromosome 2, but not in chromosomes 11 and 17. Furthermore, expression of QTL regulated in cis (cis eQTL) showed higher expression in the background genotype, which was partially explained by the presence of single nucleotide polymorphisms (SNP).

CONCLUSIONS

The literature shows limited successes from the QTL/microarray approach to identify QTL genes. Our own results from microarray profiling of three congenic strains revealed a strong tendency to select cis-eQTL over trans-eQTL. IBD regions had little effect on rate of differential expression, and we provide several reasons why IBD should not be used to discard eQTL candidates. In addition, mismatch probes produced false cis-eQTL that could not be completely removed with the current strains genotypes and low probe density microarrays. The reviewed studies did not account for lack of coverage from the platforms used and therefore removed genes that were not tested. Together, our results explain the tendency to report QTL candidates as differentially expressed and indicate that the utility of the QTL/microarray as currently implemented is limited. Alternatives are proposed that make use of microarray data from multiple experiments to overcome the outlined limitations.

Dynamic convergence and divergence of renal genomic and biological pathways in protection from Dahl salt-sensitive hypertension

Chromosome 13 consomic and congenic rat strains were analyzed to investigate the pattern of genomic pathway utilization involved in protection against salt-sensitive hypertension and renal injury. Introgression of the entire Brown-Norway chromosome 13 (consomic SS-13(BN)) or nonoverlapping segments of this chromosome (congenic strains, 16 Mbp in D13Rat151-D13Rat197 or 14 Mbp in D13Rat111-D13Got22) into the genome of the Dahl salt-sensitive rat attenuated salt-induced hypertension and proteinuria. mRNA abundance profiles in the renal cortex and the renal medulla from rats receiving 0.4% or 8% NaCl diets revealed two important features of pathway recruitment in these rat strains. First, the two congenic strains shared alterations in several pathways compared with Dahl salt-sensitive rats, despite the fact that the genomic segments introgressed in the two congenic strains did not overlap. Second, even though the genomic segment introgressed in each congenic strain was a part of the chromosome introgressed in the consomic strain, pathways altered in each congenic strain were not simply a subset of those altered in the consomic. Supporting the relevance of the mRNA data, differential expression of oxidative stress-related genes among the four strains of rats was associated with differences in urinary excretion of lipid peroxidation products. The findings suggest that different genetic alterations might converge to influence shared pathways in protection from hypertension, and that, depending on the genomic context, the same genetic alteration might diverge to affect different pathways.

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.

Positional identification of variants of Adamts16 linked to inherited hypertension

A previously reported blood pressure (BP) quantitative trait locus on rat Chromosome 1 was isolated in a short congenic segment spanning 804.6 kb. The 804.6 kb region contained only two genes, LOC306664 and LOC306665. LOC306664 is predicted to translate into A Disintegrin-like and Metalloproteinase with Thrombospondin Motifs-16 (Adamts16). LOC306665 is a novel gene. All predicted exons of both LOC306664 and LOC306665 were sequenced. Non-synonymous variants were identified in only one of these genes, LOC306664. These variants were naturally existing polymorphisms among inbred, outbred and wild rats. The full-length rat transcript of Adamts16 was detected in multiple tissues. Similar to ADAMTS16 in humans, expression of Adamts16 was prominent in the kidney. Renal transcriptome analysis suggested that a network of genes related to BP was differential between congenic and S rats. These genes were also differentially expressed between kidney cell lines with or without knock-down of Adamts16. Adamts16 is conserved between rats and humans. It is a candidate gene within the homologous region on human Chromosome 5, which is linked to systolic and diastolic BP in the Quebec Family Study. Multiple variants, including an Ala to Pro variant in codon 90 (rs2086310) of human ADAMTS16, were associated with human resting systolic BP (SBP). Replication study in GenNet confirmed the association of two variants of ADAMTS16 with SBP, including rs2086310. Overall, our report represents a high resolution positional cloning and translational study for Adamts16 as a candidate gene controlling BP.

Genome-wide identification of allelic expression in hypertensive rats

BACKGROUND

Identification of genes involved in complex cardiovascular disease traits has proven challenging. Inbred animal models can facilitate genetic studies of disease traits. The spontaneously hypertensive rat (SHR) is an inbred model of hypertension that exists in several closely related but genetically distinct lines.

METHODS AND RESULTS

We used renal gene-expression profiling across 3 distinct SHR lines to identify genes that show different expression in SHR than in the genetically related normotensive control strain, Wistar-Kyoto. To ensure robust discovery of genes showing SHR-specific expression differences, we considered only those genes in which differential expression is replicated in multiple animals of each of multiple hypertensive rat lines at multiple time points during the ontogeny of hypertension. Mutation analysis was performed on the identified genes to uncover allelic variation. We identified those genes in which all SHR lines share a single allele of the gene when normotensive controls (Wistar-Kyoto) have fixed the alternative allele. We then identified which of the differentially expressed genes show expression that is controlled by the alleleic variation present in and around the gene. Allelic expression was demonstrated by observing the effect on gene expression of alleles inherited in the freely segregating F(2) progeny of a cross between SHR and Wistar-Kyoto animals.

CONCLUSIONS

The result of these studies is the identification of several genes (Ptprj, Ela1, Dapk-2, and Gstt2) in which each of 4 SHR lines examined have fixed the same allele and in which each of 2 Wistar-Kyoto lines have a contrasting allele for which the inherited allele influences the level of gene expression. We further show that alleles of these genes lie in extensive haplotype blocks that have been inherited identical by descent in the hypertensive lines.

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.

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.

Haplotypic analysis of Wellcome Trust Case Control Consortium data

We applied a recently developed multilocus association testing method (localized haplotype clustering) to Wellcome Trust Case Control Consortium data (14,000 cases of seven common diseases and 3,000 shared controls genotyped on the Affymetrix 500 K array). After rigorous data quality filtering, we identified three disease-associated loci with strong statistical support from localized haplotype cluster tests but with only marginal significance in single marker tests. These loci are chromosomes 10p15.1 with type 1 diabetes (p = 5.1 x 10(-9)), 12q15 with type 2 diabetes (p = 1.9 x 10(-7)) and 15q26.2 with hypertension (p = 2.8 x 10(-8)). We also detected the association of chromosome 9p21.3 with type 2 diabetes (p = 2.8 x 10(-8)), although this locus did not pass our stringent genotype quality filters. The association of 10p15.1 with type 1 diabetes and 9p21.3 with type 2 diabetes have both been replicated in other studies using independent data sets. Overall, localized haplotype cluster analysis had better success detecting disease associated variants than a previous single-marker analysis of imputed HapMap SNPs. We found that stringent application of quality score thresholds to genotype data substantially reduced false-positive results arising from genotype error. In addition, we demonstrate that it is possible to simultaneously phase 16,000 individuals genotyped on genome-wide data (450 K markers) using the Beagle software package.

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.

Fine-mapping and comprehensive transcript analysis reveals nonsynonymous variants within a novel 1.17 Mb blood pressure QTL region on rat chromosome 10

The presence of blood pressure (BP) quantitative trait loci (QTL) on rat chromosome 10 has been clearly demonstrated by linkage analysis and substitution mapping. 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 to fine-map a chromosome 10 region (QTL1) linked to blood pressure. Comparison of seven congenic substrains refined QTL1 to a 1.17 Mb segment flanked by D10Mco88 and D10Mco89, which are located at 71,513,116 and 72,684,774 bp, respectively. The newly defined QTL1, containing 18 genes, is captured in its entirety within a single congenic substrain. A thorough transcript analysis revealed that 3 of these 18 genes, Ccl5, Ddx52, and RGD1559577, had nonsynonymous allelic variations between the S rat and the LEW rat. None of the detected transcripts within the newly defined QTL1 are implicated directly in BP control in humans or model organisms. Therefore, the present work defines a novel blood pressure QTL with three potential quantitative trait nucleotides.

Combined application of behavior genetics and microarray analysis to identify regional expression themes and gene-behavior associations

In this report we link candidate genes to complex behavioral phenotypes by using a behavior genetics approach. Gene expression signatures were generated for the prefrontal cortex, ventral striatum, temporal lobe, periaqueductal gray, and cerebellum in eight inbred strains from priority group A of the Mouse Phenome Project. Bioinformatic analysis of regionally enriched genes that were conserved across all strains revealed both functional and structural specialization of particular brain regions. For example, genes encoding proteins with demonstrated anti-apoptotic function were over-represented in the cerebellum, whereas genes coding for proteins associated with learning and memory were enriched in the ventral striatum, as defined by the Expression Analysis Systematic Explorer (EASE) application. Association of regional gene expression with behavioral phenotypes was exploited to identify candidate behavioral genes. Phenotypes that were investigated included anxiety, drug-naive and ethanol-induced distance traveled across a grid floor, and seizure susceptibility. Several genes within the glutamatergic signaling pathway (i.e., NMDA/glutamate receptor subunit 2C, calmodulin, solute carrier family 1 member 2, and glutamine synthetase) were identified in a phenotype-dependent and region-specific manner. In addition to supporting evidence in the literature, many of the genes that were identified could be mapped in silico to surrogate behavior-related quantitative trait loci. The approaches and data set described herein serve as a valuable resource to investigate the genetic underpinning of complex behaviors.

Cd14,Gbp1, andPla2g2a: three major candidate genes for experimental IBD identified by combining QTL and microarray analyses

Induction of inflammatory bowel (IBD)-like disease in mice by a targeted mutation in the Il10 gene ( Il10−/−) is inbred strain dependent. C3H/HeJBir (C3) mice are colitis susceptible, whereas C57BL/6J (B6) mice are resistant. Genetic dissection of this susceptibility revealed 10 colitogenic quantitative trait loci (QTL). The aim of this study was to identify valuable candidate genes by a combination of QTL mapping and microarray analyses. Sixteen genes were differentially expressed between B6- and C3- Il10−/−mice and were located within the QTL intervals. Three major candidate genes ( Cd14, Gbp1, Pla2g2a) showed prominent expression differences between B6- and C3- Il10−/−as well as between B6 and C3 wild-type mice, which was confirmed by semiquantitative or real-time RT-PCR. Because strain differences are known for Gbp1 and Pla2g2a, further analyses focused on Cd14. Western blot analysis revealed strain differences also on the protein level. Cd14 expression in animals with defective and intact Toll-like receptor (TLR)4 signaling (C3, C3H/HeN, B6, B6- Tlr4tm1Aki) make the TLR4 defect of C3 mice unlikely to be the reason for higher Cd14 expression. Less Cd14 expression in germ-free mice indicates a contribution of the microflora on Cd14 expression. Stimulation of naive peritoneal macrophages with bacterial antigens showed lower CD14 surface expression in B6 than in C3 mice. In conclusion, the large number of candidate genes was reduced to three major candidates that play an important role in inflammatory processes and immune response. Strain differences for them are already known or are shown in this study.