Insights into Dahl salt-sensitive hypertension revealed by temporal patterns of renal medullary gene expression

Mingyu Liang, Baozhi Yuan, Elizabeth Rute, Andrew S. Greene, Michael Olivier and Allen W. Cowley, Jr.

Department of Physiology, Medical College of Wisconsin, Milwaukee, Wisconsin 53226


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS AND MATERIALS
 RESULTS
 DISCUSSION
 References
 
Dahl salt-sensitive SS and consomic, salt-resistant SS-13BN/Mcw rats possess a highly similar genetic background but exhibit substantial differences in blood pressure salt sensitivity. We used cDNA microarrays to examine sequential changes of mRNA expression of ~2,000 currently known rat genes in the renal medulla (a tissue critical for long-term blood pressure regulation) in SS and SS-13BN/Mcw rats in response to a high-salt diet (16 h, 3 days, or 2 wk). Differentially expressed genes in each between-group comparison were identified based on a threshold determined experimentally using a reference distribution that was constructed by comparing rats within the same group. A difference analysis of 54 microarrays identified 50 genes that exhibited the most distinct temporal patterns of expression between SS and SS-13BN/Mcw rats over the entire time course. Thirty of these genes could be linked to the regulation of arterial blood pressure or renal injury based on their known involvement in functional pathways such as renal tubular transport, metabolism of vasoactive substances, extracellular matrix formation, and apoptosis. Importantly, the majority of the 30 genes exhibited temporal expression patterns that would be expected to lower arterial pressure and reduce renal injury in SS-13BN/Mcw compared with SS rats. The phenotypic impact of the other 20 genes was less clear. These 50 genes are widely distributed on chromosome 13 and several other chromosomes. This suggested that primary genetic defects, although important, are unlikely to be solely responsible for the full manifestation of this type of hypertension and associated injury phenotypes. In summary, the results of this study identified a number of pathways potentially important for the amelioration of hypertension and renal injury in SS-13BN/Mcw rats, and these results generated a series of testable hypotheses related to the role of the renal medulla in the complex mechanism of salt-sensitive hypertension.

microarray; 11ß-hydroxysteroid dehydrogenase; glucagon receptor; extracellular matrix; apoptosis


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS AND MATERIALS
 RESULTS
 DISCUSSION
 References
 
ENHANCED SENSITIVITY OF BLOOD pressure to dietary salt intake is observed in a significant portion of hypertensive patients, especially in the African-American population (22). To explore the role of gene-environment interactions in salt-sensitive hypertension, we have developed an inbred consomic rat strain (SS-13BN/Mcw) in which chromosome 13 of the salt-resistant Brown Norway rat (BN) was substituted into the genetic background of the Dahl salt-sensitive rat (SS). A substantial reduction of blood pressure salt sensitivity and renal injury was found in SS-13BN/Mcw rats compared with SS rats (4).

This consomic strain can be used to develop congenic strains to identify small regions and eventually specific gene(s) on SS chromosome 13 that may contain primary genetic defects associated with blood pressure salt sensitivity. However, primary genetic defects, although very important, are unlikely to be solely responsible for the full manifestation of hypertension and associated injury phenotypes. These primary genetic defects could directly or indirectly affect the expression of genes across the genome. Functional pathways involving those genes could then importantly contribute to the development of salt-sensitive hypertension and associated injury. Due to the high degree of genetic identity and the substantial difference of blood pressure phenotypes, the comparison of SS-13BN/Mcw and SS provides an excellent tool with which to study these mechanistic pathways. To explore these pathways, a custom-made cDNA microarray representing ~2,000 rat genes has been developed (13). These genes have been selected because most of them have some known biological functions, facilitating the effort focused on linking gene expression patterns with functional pathways. It also reduces the cost of these experiments, allowing studies with more replicates to be performed. It is important to note that only about 40 genes have been mapped to rat chromosome 13 (see Rat Genome Database, http://rgd.mcw.edu), and only 13 of them are included in this microarray. Therefore, this microarray is used primarily as a tool to study the mechanistic pathways of Dahl salt-sensitive hypertension including those involving genome-wide gene interactions, rather than identify the primary genetic defects on SS chromosome 13 that may initiate Dahl salt-sensitive hypertension.

In a previous study, we used this ~2,000-gene microarray to examine gene expression profiles in the renal medulla (a tissue playing a critical role in the long-term regulation of arterial blood pressure; Refs. 5 and 9) in SS-13BN/Mcw and SS rats on either a low-salt diet or switched to a high-salt diet for 2 wk (13). It is known that SS rats exhibit a rapid reduction of medullary blood flow and develop hypertension rapidly and progressively when exposed to a high-salt diet (15). The present study was, therefore, designed to examine sequential changes of renal medullary gene expression in SS-13BN/Mcw and SS rats upon exposure to the high-salt diet for 16 h (overnight), 3 days, or 2 wk. These time points represented initial high-salt exposure without significant hypertension (16 h), the beginning of significant rise of arterial blood pressure (3 days), and established hypertension (2 wk) in SS rats (4, 15). Analysis of 54 microarrays identified 50 genes that exhibited the most distinct temporal patterns of expression between the two strains of rats over the time course studied. These expression patterns agreed well with the phenotypic differences between SS-13BN/Mcw and SS and suggested several functional pathways that might be important for the development of salt-sensitive hypertension and associated renal injury.


    METHODS AND MATERIALS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS AND MATERIALS
 RESULTS
 DISCUSSION
 References
 
Animals used.
Chromosome 13 from the salt-resistant Brown Norway (BN) rat was substituted for chromosome 13 of the Dahl salt-sensitive (SS) rat of the Medical College of Wisconsin (MCW) colony through selective breeding as described previously (4). The resulting consomic SS-13BN/Mcw rat, therefore, contained homozygous SS alleles on every chromosome, as confirmed by total genome scans, except chromosome 13, which was homozygous BN at all alleles as confirmed by a scan at higher density. Analysis of 1,541 microsatellite markers indicates that SS-13BN/Mcw has only 1.95% allelic differences compared with SS, whereas other salt-resistant strains are at least 30% different from SS. The consomic line was maintained by brother-sister mating.

Dietary salt intake protocol and collection of tissue.
Nine SS rats and nine SS-13BN/Mcw rats were maintained on a normal light cycle (light from 6:00 AM to 6:00 PM) and on a low-salt (0.4% NaCl) diet until 11–12 wk of age. Three rats of each strain were then switched to a high-salt (4% NaCl) diet for 16 h (from 5:00 PM to 9:00 AM the next day), while another three rats switched to the high-salt diet for 3 days. On the day of tissue harvest, rats were anesthetized with Inactin (25 mg/kg ip) and ketamine (20 mg/kg im). Kidneys were removed quickly, and the renal medulla (both papilla and outer medulla) was selectively dissected using a sterile surgical blade and scissors, snap frozen in liquid nitrogen, and stored at -80°C. All tissue samples were collected between 9:00 AM and noon. Tissues were collected within 15 min following the administration of anesthesia. Therefore, anesthesia should have minimal confounding effects on gene expression, although they could not be completely ruled out.

Construction and hybridization of rat known-gene cDNA microarrays.
The description of microarray experiments in the present study conforms to the "minimum information about a microarray experiment" (MIAME) standard (2) to the extent possible. This includes detailed descriptions of custom-made microarrays, experimental design, hybridization protocols, and data analysis methods.

Custom-made cDNA microarrays containing cDNA probes that represented ~80% of all currently known rat genes were used in the present study. The majority of these genes have some known biological functions in rats. These microarrays were constructed as we described previously (13) with the following modifications. cDNA probes for 1,871 genes, instead of 1,751 genes, were printed on this array, which included 1,687 clones purchased from Research Genetics (Huntsville, AL), and 184 genes cloned in our department or purchased from ATCC (Manassas, VA).

Microarrays were processed and hybridized using a two-color (Cy3 and Cy5) method with dye switching as we described previously (13).

Experimental design and overall analytical strategy.
Rats were paired for microarray hybridization, and the results were analyzed as shown in Fig. 1. Each of the four between-group comparisons comprised three pairs of rats examined by six microarrays with dye switching for each pair. Dye switching was not necessary for within-group comparisons, because the two rats of each pair were equivalent in terms of their treatment status. A total of 30 microarrays were hybridized in the present study. The expression data from the present study were combined with those from 24 microarrays hybridized in our previous study (13) and used to identify genes exhibiting the most distinct temporal patterns of expression between SS and SS-13BN/Mcw over the time course studied.



View larger version (38K):
[in this window]
[in a new window]
 
Fig. 1. Experimental design and overall analytical strategy. Renal medullary mRNA expression was compared by cDNA microarrays in 8 groups, each containing 3 individual rats. Solid lines represent between-group hybridizations performed in the present study, with each comparison comprising 6 microarray hybridizations. Dash lines represented between-group hybridizations performed in the previous study (13). Dotted lines represented within-group hybridizations performed in the present study. Within-group comparisons were performed to construct a "reference distribution" of gene expression ratios that was used to determine the threshold of differential expression as detailed in RESULTS. SS, SS rats; 13, SS-13BN/Mcw rats; LS, maintained on the low-salt (0.4% NaCl) diet; HS, switched to the high-salt (4% NaCl) diet for the time period indicated, i.e., 16 h, 3 d (3 days), and 2 wk.

 
Microarray data analysis.
Raw values of fluorescent intensity in each spot were obtained from microarray images using ImaGene 4.01 software (BioDiscovery, Los Angeles, CA). These raw data were categorized, selected, and adjusted to yield natural log-transformed, normalized ratios following the systematic method that we described previously (13). This data analysis method quantitatively identifies and systematically excludes spots with low intensities or high local background to avoid generations of disproportionate and misleading ratios. Approximately 1,000 genes passed the stringent data selection process in at least 5 of the 6 microarrays in each between-group comparison. Differentially expressed genes were identified based on a threshold that was determined experimentally using the "reference distribution" as detailed in RESULTS.

A difference calculation was used to identify genes with the most distinct temporal patterns of expression between SS and SS-13BN/Mcw rats over the time course studied. This calculation measured the sum of the squared differences between the expression levels of each gene in the two rat strains at four time points (low-salt, high-salt 16 h, 3 days, and 2 wk), which is essentially a variant of the Euclidean distance calculation. Specifically, the expression level of a gene in SS rats on the low-salt diet (SSLS) was arbitrarily set at 0. The Ln(ratio) for this gene between each high-salt time point (16 h, 3 days, or 2 wk) and SSLS was used to indicate its relative expression level in SS rats at each high-salt time point. Relative expression levels of this gene in SS-13BN/Mcw rats at those four time points were similarly obtained with the expression level on the low-salt diet (13LS) set at 0. The difference of this gene’s expression pattern between SS and SS-13BN/Mcw rats over the time course was then calculated as D = {sum}4i = 1(Xi - Yi)2, in which Xi and Yi were the expression levels of this gene (obtained as described above) in SS and SS-13BN/Mcw rats, respectively, at the four time points. To take into consideration the fact that genes were not always expressed equally between SSLS and 13LS, relative expression levels in SS-13BN/Mcw were corrected by adding the Ln(ratio) between 13LS and SSLS. A second set of D values was calculated using these corrected expression levels. This second set of D values theoretically would better reflect the true expression differences. However, it had the disadvantage of being disproportionally affected by the single Ln(ratio) between 13LS and SSLS. Therefore, we used both sets of D values to rank the genes rather than relying on one of them. A gene was considered to exhibit highly distinct temporal patterns of expression between SS and SS-13BN/Mcw if its D values appeared in the largest 10% of both sets of D values and if it did not have any missing X or Y values over the entire time course. A cutoff at the largest 10% was chosen to yield a manageable set of genes. The validity of this cutoff level was supported by the fact that the majority of the genes identified were differentially expressed at least at one particular time point.

Northern blot.
mRNA expressions of 11 genes were further examined by Northern blotting as we described previously (13, 26).

Retrieval of gene mapping information.
Rat chromosomal locations of genes of interest were retrieved from the Rat Genome Database (RGD) maintained at our institution (http://rgd.mcw.edu) and the RGD Virtual Comparative Mapping Tool. This tool enables prediction of gene locations on the rat genome based on human or mouse genomes for those genes that were not directly mapped to the rat genome by genetic linkage or radiation hybrid methods.

Verification of clone identity.
Some of the genes identified as differentially expressed in the present study had been sequence-verified in our previous study (13). For the rest of the differentially expressed genes, PCR products were generated from corresponding clones and sequenced using the BDT chemistry (Applied Biosystems, Foster City, CA). Fifty-four partial sequences were successfully obtained and, based on BLAST searches, three (5.6%) did not match the identification provided by Research Genetics. These three clones were removed from the list of differentially expressed genes.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS AND MATERIALS
 RESULTS
 DISCUSSION
 References
 
Identification of differentially expressed genes based on a threshold determined experimentally using a reference distribution.
A "reference distribution" approach was used in the present study to determine experimentally the threshold of differential expression. To estimate the magnitude of random variation, individual samples from the same group were compared against each other using microarray hybridization as shown in Fig. 1. When the Ln(ratio) values obtained from the six microarrays for within-group comparisons (calculated by alternating Cy3/Cy5 and Cy5/Cy3) were averaged for each gene, the averaged Ln(ratio) values followed a normal distribution with a standard deviation of 0.100. When the three microarrays from SS rats and the three from SS-13BN/Mcw rats were analyzed separately, both followed a normal distribution with respective standard deviations of 0.136 and 0.129. To be conservative, the distribution from the SS rats was utilized as the "reference distribution." Application of this same analysis to the results from microarrays for the between-group comparisons (Fig. 1) yielded standard deviations of 0.200, 0.187, 0.305, and 0.238 for comparisons of SSHS 16 h vs. SSLS, SSHS 3 days vs. SSLS, 13HS 16 h vs. 13LS, and 13HS 3 days vs. 13LS. An example contrasting the distribution from SSHS 16 h vs. SSLS with the reference distribution is shown in Fig. 2.



View larger version (20K):
[in this window]
[in a new window]
 
Fig. 2. The Ln(ratio) distribution from a between-group comparison is compared with the reference distribution constructed from within-group comparisons. Total numbers of genes (frequency) were scaled to be the same for the two distributions. The averaged Ln(ratio) values in both distributions were normalized to zero. Standard deviations of the reference distribution and the distribution of SSLS vs. SSHS 16-hr were 0.136 and 0.200, respectively. SSLS, SS rats maintained on the low-salt (0.4% NaCl) diet; SSHS 16-hr, SS rats switched to the high-salt (4% NaCl) diet for 16 h (overnight).

 
On the basis of the reference distribution, the threshold of differential expression was determined to be Ln(ratio) values of ±0.447 (or ratios above 1.56 or below 0.64). This threshold corresponded to the 99.9% interval of the reference distribution. When this threshold was applied to each of the four between-group comparisons, 46 genes were differentially expressed between SSHS 16 h and SSLS, 31 genes between SSHS 3 days and SSLS, 109 genes between 13HS 16 h and 13LS, and 43 genes between 13HS 3 days and 13LS. These genes, together with their clone ID and Ln(ratio) values, are presented in the Supplemental Table, which is published at the Physiological Genomics web site.1 When the data from our previous study (13) were reexamined, it was found that the threshold value used previously (two times the standard deviation) was rather close to ±0.447, supporting the validity of the standard deviation approach in the specific context of our previous study.

Ln(ratio) values of the expression of 11 genes obtained from both microarrays and Northern blots are shown in Table 1. Ln(ratio) values from microarrays significantly correlated with those from Northern blots with a correlation coefficient of 0.376 (P = 0.01). Moreover, 9 of the 11 genes, with dipeptidyl peptidase and rat prolyl oligopeptidase (rPOP) as the exceptions, were consistently shown to be expressed either higher or lower in SS-13BN/Mcw compared with SS by both microarrays and Northern blots. In this and our previous (13, 27) studies, a total of 62 genes have been examined by both microarrays and Northern blots, and significant correlation between microarray and Northern blot results has been consistently observed. This represents one of the largest validation data sets in the microarray literature, and supports the overall reliability of the microarray technique we used.


View this table:
[in this window]
[in a new window]
 
Table 1. Ln(ratio)s obtained from microarrays and Northern blots

 
Genes exhibiting the most distinct temporal patterns of expression between SS and SS-13BN/Mcw rats.
Data obtained from the present study were combined with those from the previous study (13) for identifying genes exhibiting the most distinct temporal patterns of expression between SS and SS-13BN/Mcw over the time course studied. A total of 54 microarrays were analyzed using the difference calculation described in METHODS AND MATERIALS. These analyses resulted in the identification of 50 genes that were considered to have the most distinct temporal patterns of renal medullary mRNA expression in SS compared with SS-13BN/Mcw rats over the entire time course studied. Many of these genes were considered differentially expressed at one or more particular time points according to the threshold of differential expression determined in this or the previous (13) studies.

Extensive searches and evaluations of literature and publicly available databases were performed to retrieve information regarding the biological functions of proteins encoded by these 50 genes. Fifteen of these 50 genes, shown in Fig. 3, were found to be involved in the action or metabolism of substances that regulate renal tubular reabsorption and/or vascular tone. Examples of these substances included corticosterone, aldosterone, glucagon, natriuretic peptides, and arachidonic acid metabolites.



View larger version (41K):
[in this window]
[in a new window]
 
Fig. 3. Genes with known functions related to blood pressure regulation and exhibiting distinct temporal patterns of mRNA expression between SS (•) and SS-13BN/Mcw rats ({circ}). LS, maintained on the low-salt diet; HS, switched to the high-salt diet for the time period indicated. Ln(ratio) values shown are between each condition and SS on the low-salt diet. Temporal expression patterns of the majority of these genes are consistent with the lower blood pressure seen in SS-13BN/Mcw compared with SS when both strains of rats are exposed to the high-salt diet. rPOP, rat prolyl oligopeptidase.

 
Another 15 genes with distinct temporal patterns of expression could be linked to the regulation of renal injury and perhaps arterial blood pressure. As shown in Fig. 4, these genes were involved in apoptosis, extracellular matrix formation and fibrosis, oxidative stress and toxicity, and cellular growth and proliferation.



View larger version (37K):
[in this window]
[in a new window]
 
Fig. 4. Genes with known functions related to tissue injury or blood pressure regulation and exhibiting distinct temporal patterns of mRNA expression between SS (•) and SS-13BN/Mcw rats ({circ}). Ln(ratio) values shown are between each condition and SS on the low-salt diet. Temporal expression patterns of the majority of these genes indicate lower apoptotic activity in SS-13BN/Mcw and are consistent with the lower blood pressure and less renal interstitial fibrosis seen in SS-13BN/Mcw compared with SS when both strains of rats are exposed to the high-salt diet.

 
Another 14 genes with distinct temporal patterns of expression could not be directly linked to the regulation of arterial blood pressure or renal injury based on their known biological functions. Six genes that did not have known functions also exhibited distinct temporal patterns of expression between SS and SS-13BN/Mcw rats. These 20 genes and their temporal expression patterns are shown in Fig. 5.



View larger version (38K):
[in this window]
[in a new window]
 
Fig. 5. Genes exhibiting distinct temporal patterns of mRNA expression between SS (•) and SS-13BN/Mcw rats ({circ}), but without apparent links to blood pressure regulation or tissue injury, or with unknown functions. Ln(ratio) values shown are between each condition and SS on the low-salt diet.

 
Thirty-eight of the 50 genes with distinct temporal patterns of expression (Figs. 3, 4, and 5) have been directly or comparatively mapped on rat chromosomes. Three of them, 11ß-hydroxysteroid dehydrogenase (11ß-HSD) type 1, NGF-inducible antiproliferative putative secreted protein, and thiol-specific antioxidant protein, are mapped to rat chromosome 13. The rest are distributed widely on many chromosomes.


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS AND MATERIALS
 RESULTS
 DISCUSSION
 References
 
The results of the present study identified a number of genes and pathways potentially important for the development of Dahl salt-sensitive hypertension and associated renal injury. The temporal expression patterns of the majority of the 15 genes shown in Fig. 3 would be expected to decrease renal tubular reabsorption of fluid and sodium and reduce vascular tone in SS-13BN/Mcw compared with SS. Those shown in Fig. 4 would in general suggest that the renal medulla of SS-13BN/Mcw rats has lower apoptotic activity and reduced extracellular matrix formation compared with SS. This was consistent with the substantial reduction of blood pressure salt sensitivity and renal injury in SS-13BN/Mcw that we previously observed (4). Conversely, the abnormal regulation of these genes in SS could contribute to the development of salt-sensitive hypertension and associated renal injury.

Figure 6 depicts a hypothetical scheme summarizing several of these pathways that seem to be convincingly related to the amelioration of salt-sensitive hypertension and renal injury in SS-13BN/Mcw rats compared with SS rats, according to the known functions of those genes (1, 3, 68, 1012, 1618, 20, 21, 2325). An interesting example is the 11ß-HSD pathway of glucocorticoid metabolism. In rats, 11ß-HSD type 1 has the ability to convert inert 11-dehydrocorticosterone (functionally similar to cortisone in humans) to active corticosterone (functionally similar to cortisol in humans), whereas type 2 inactivates corticosterone back to 11-dehydrocorticosterone. Cortisol is a potent stimulator of mineralocorticoid receptors located in the distal nephron (7), leading to sodium retention. Homozygous loss-of-function mutations in the 11ß-HSD type 2 gene were found to cause the human syndrome of apparent mineralocorticoid excess characterized with hypertension (16). The downregulation of 11ß-HSD type 1 and the upregulation of 11ß-HSD type 2 in SS-13BN/Mcw compared with SS shown in the present study, therefore, would be expected to decrease the level of active glucocorticoids in the kidney of SS-13BN/Mcw rats. This is consistent with lower arterial blood pressure observed in SS-13BN/Mcw compared with SS (4). Specific hypotheses could be developed to test the functional importance of this and many other pathways (Figs. 36) in Dahl salt-sensitive hypertension.



View larger version (21K):
[in this window]
[in a new window]
 
Fig. 6. Potential pathways linking gene expression patterns to the regulation of blood pressure or renal injury. Names of genes are shown in bold font; {uparrow}, upregulation or increases in SS-13BN/Mcw compared with SS; {downarrow}, downregulation or decreases in SS-13BN/Mcw compared with SS.

 
Several additional points emerged from these results. It appears that multiple pathways may contribute to the development of Dahl salt-sensitive hypertension in a temporally specific manner. For example, several genes were differentially expressed between SS and SS-13BN/Mcw only at early time points of 16 h or 3 days, whereas others exhibited differential expression only after 2 wk of high-salt exposure. It is also interesting to note that several genes related to renal injury were differentially expressed at early time points (Fig. 4). This suggested that the renal injury observed in SS rats could be the result of both the high-salt diet and hypertension, rather than hypertension alone. These results indicate the existence of a complex network of pathways that interact with each other and collectively contribute to the full manifestation of this specific type of hypertensive phenotypes. It emphasizes the importance of examining sequential changes of gene expression and the integration of multiple pathways for understanding the pathophysiology of complex diseases such as salt-sensitive hypertension.

The fact that only a few of the 50 genes with distinct temporal patterns of expression were mapped to rat chromosome 13 is reminiscent of the results of our previous study (13). This could be caused by several mechanisms. Genomic differences on chromosome 13 and the products they encode could be regulatory factors that can directly influence the expression of genes on many chromosomes. These genome-wide gene interactions could also be indirect through downstream signaling pathways or as a result of feedback regulation by functional alterations. In addition, the possibility of residual BN alleles on chromosomes other than chromosome 13 in SS-13BN/Mcw rats cannot be excluded. It remains to be determined whether the differentially expressed genes that did map to chromosome 13 in the present study indeed contain fundamental genetic differences associated with blood pressure salt sensitivity. Approaches such as generating congenic lines and sublines, gene targeting or transfer, and mutation identification will be useful in this regard (19). In the meantime, multiple pathways identified in the present study, whether located on chromosome 13 or not, should improve our understanding of the complex network of pathways underlying the pathophysiology of Dahl salt-sensitive hypertension. The importance of these genes is supported not only by the theoretical presence of such a network of pathways but also by the fact that the expression patterns of the majority of these genes were consistent with the blood pressure and injury phenotypes in these two strains of rats. This is interesting, particularly because it is not unusual for microarray studies to generate expression profiles that are difficult to reconcile with known phenotypes.

A novel "reference distribution" method was used to determine the threshold of differential expression in the present study. With this method, an apparent difference between two groups must exceed the random variation between individual samples in each group to be considered a true difference. This, indeed, is a common principle behind most statistical tests in biological and other fields of investigation. The ratio distributions for between-group comparisons had consistently larger standard deviations than the reference distribution, indicating that the magnitude of between-group differences in mRNA expression exceeded that of random variations. With the 99.9% interval used in the present study, not more than 1 gene would be expected to be identified as differentially expressed by chance alone. If one would like to reduce the number of false negatives while tolerating a larger number of false positives, then the threshold could be reduced, e.g., to a value corresponding to a 95% interval of the reference distribution, knowing that they would now have ~50 (5% of 1,000) false-positive genes.

In summary, cDNA microarrays with newly developed analytical methods have been used to examine sequential changes of renal medullary gene expressions in SS and consomic SS-13BN/Mcw rats. The results identified a number of pathways potentially important for the amelioration of hypertension and renal injury in SS-13BN/Mcw rats and generated a series of testable hypotheses related to the role of the renal medulla in the complex mechanism of salt-sensitive hypertension. These results also indicated that multiple pathways could be involved in the development of Dahl salt-sensitive hypertension and suggested that primary genetic defects, although important, are unlikely to be solely responsible for the full manifestation of this type of hypertension and associated injury phenotypes.


    ACKNOWLEDGMENTS
 
We gratefully acknowledge Regina Cole for assistance with clone resequencing, Dr. Sandra do Lia Amaral and Katherine Fredrich for assistance with microarray preparations, Gina Tadisch for maintaining rat colonies, Randy Berdan and SPS Productions for assistance with computer program development for compiling gene identification files, and Meredith Skelton for critical review of the manuscript.

This study was supported by National Heart, Lung, and Blood Institute Grants HL-66579, HL-54998, and HL-29587.


    FOOTNOTES
 
Article published online before print. See web site for date of publication (http://physiolgenomics.physiology.org).

Address for reprint requests and other correspondence: M. Liang, Dept. of Physiology, Medical College of Wisconsin, 8701 Watertown Plank Road, Milwaukee, WI 53226 (E-mail: mliang{at}mcw.edu).

10.1152/physiolgenomics.00089.2002.

1 The Supplementary Material for this article is available online at http://physiolgenomics.physiology.org/cgi/content/full/12/3/229/DC1. Back


    References
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS AND MATERIALS
 RESULTS
 DISCUSSION
 References
 

  1. Ahloulay M, Dechaux M, Hassler C, Bouby N, and Bankir L. Cyclic AMP is a hepatorenal link influencing natriuresis and contributing to glucagon-induced hyperfiltration in rats. J Clin Invest 98: 2251–2258, 1996.[Abstract/Free Full Text]
  2. Brazma A, Hingamp P, Quackenbush J, Sherlock G, Spellman P, Stoeckert C, Aach J, Ansorge W, Ball CA, Causton HC, Gaasterland T, Glenisson P, Holstege FC, Kim IF, Markowitz V, Matese JC, Parkinson H, Robinson A, Sarkans U, Schulze-Kremer S, Stewart J, Taylor R, Vilo J, and Vingron M. Minimum information about a microarray experiment (MIAME): toward standards for microarray data. Nat Genet 29: 365–371, 2001.[ISI][Medline]
  3. Chambers SM and Morris BJ. Glucagon receptor gene mutation in essential hypertension. Nat Genet 12: 122, 1996.[ISI][Medline]
  4. Cowley AW Jr, Roman RJ, Kaldunski ML, Dumas P, Dickhout JG, Greene AS, and Jacob HJ. Brown Norway chromosome 13 confers protection from high salt to consomic Dahl S rat. Hypertension 37: 456–461, 2001.[Abstract/Free Full Text]
  5. Cowley AW Jr and Roman RJ. The role of the kidney in hypertension. JAMA 275: 1581–1589, 1996.[ISI][Medline]
  6. Franco-Saenz R, Shen P, Lee SJ, Cicila GT, and Henrich WL. Regulation of the genes for 11ß-hydroxysteroid dehydrogenase type 1 and type 2 in the kidney of the Dahl rat. J Hypertens 17: 1089–1093, 1999.[ISI][Medline]
  7. Funder JW, Pearce PT, Smith R, and Smith AI. Mineralocorticoid action: target tissue specificity is enzyme, not receptor, mediated. Science 242: 583–585, 1988.[ISI][Medline]
  8. Green DR. Apoptotic pathways: paper wraps stone blunts scissors. Cell 102: 1–4, 2000.[ISI][Medline]
  9. Guyton AC. Blood pressure control: special role of the kidneys and body fluids. Science 252: 1813–1816, 1991.[ISI][Medline]
  10. Hamet P, Thorin-Trescases N, Moreau P, Dumas P, Tea BS, deBlois D, Kren V, Pravenec M, Kunes J, Sun Y, and Tremblay J. Workshop: excess growth and apoptosis: is hypertension a case of accelerated aging of cardiovascular cells? Hypertension 37: 760–766, 2001.[Abstract/Free Full Text]
  11. Haraguchi M, Border WA, Huang Y, and Noble NA. t-PA promotes glomerular plasmin generation and matrix degradation in experimental glomerulonephritis. Kidney Int 59: 2146–2155, 2001.[ISI][Medline]
  12. Li P, Nijhawan D, Budihardjo I, Srinivasula SM, Ahmad M, Alnemri ES, and Wang X. Cytochrome c and dATP-dependent formation of Apaf-1/caspase-9 complex initiates an apoptotic protease cascade. Cell 91: 479–489, 1997.[ISI][Medline]
  13. Liang M, Yuan B, Rute E, Greene AS, Zou AP, Soares P, McQuestion GD, Slocum GR, Jacob HJ, and Cowley AW Jr. Renal medullary genes in salt-sensitive hypertension: a chromosomal substitution and cDNA microarray study. Physiol Genomics 8: 139–149, 2002. First published January 2, 2002; 10.1152/physiolgenomics.00083.2001.[Abstract/Free Full Text]
  14. Marshall E. Rat genome spurs an unusual partnership. Science 291: 1872, 2001.[Free Full Text]
  15. Miyata N and Cowley AW Jr. Renal intramedullary infusion of L-arginine prevents reduction of medullary blood flow and hypertension in Dahl salt-sensitive rats. Hypertension 33: 446–450, 1999.[Abstract/Free Full Text]
  16. Mune T, Rogerson FM, Nikkila H, Agarwal AK, and White PC. Human hypertension caused by mutations in the kidney isozyme of 11ß-hydroxysteroid dehydrogenase. Nat Genet 10: 394–399, 1995.[ISI][Medline]
  17. Napirei M, Karsunky H, Zevnik B, Stephan H, Mannherz HG, and Moroy T. Features of systemic lupus erythematosus in Dnase1-deficient mice. Nat Genet 25: 135–136, 2000.[ISI][Medline]
  18. Pohlova I, Miksik I, Kunes J, and Pacha J. 11ß-Hydroxysteroid dehydrogenase activity in spontaneously hypertensive and Dahl rats. Am J Hypertens 13: 927–933, 2000.[ISI][Medline]
  19. Pravenec M, Wallace C, Aitman TJ, and Kurtz TW. Gene expression profiling in hypertension research. A critical perspective. Hypertension 41: 3–8, 2003. Published online before print December 2, 2002, 10.1161/01.HYP.0000050927.96979.41.[Abstract/Free Full Text]
  20. Seckl JR and Walker BR. Minireview: 11ß-hydroxysteroid dehydrogenase type 1: a tissue-specific amplifier of glucocorticoid action. Endocrinology 142: 1371–1376, 2001.[Abstract/Free Full Text]
  21. Stewart PM, Wallace AM, Valentino R, Burt D, Shackleton CH, and Edwards CR. Mineralocorticoid activity of liquorice: 11-ß-hydroxysteroid dehydrogenase deficiency comes of age. Lancet 2: 821–824, 1987.[ISI][Medline]
  22. Sullivan JM, Prewitt RL, and Ratts TE. Sodium sensitivity in normotensive and borderline hypertensive humans. Am J Med Sci 295: 370–377, 1988.[ISI][Medline]
  23. Takeda Y, Inaba S, Furukawa K, and Miyamori I. Renal 11ß-hydroxysteroid dehydrogenase in genetically salt-sensitive hypertensive rats. Hypertension 32: 1077–1082, 1998.[Abstract/Free Full Text]
  24. Weber MA. Vasopeptidase inhibitors. Lancet 358: 1525–1532, 2001.[ISI][Medline]
  25. Yasutomo K, Horiuchi T, Kagami S, Tsukamoto H, Hashimura C, Urushihara M, and Kuroda Y. Mutation of DNASE1 in people with systemic lupus erythematosus. Nat Genet 28: 313–314, 2001.[ISI][Medline]
  26. Yuan B and Cowley AW Jr. Evidence that reduced renal medullary nitric oxide synthase activity of Dahl S rats enables small elevations of arginine vasopressin to produce sustained hypertension. Hypertension 37: 524–528, 2001.[Abstract/Free Full Text]
  27. Yuan B, Liang M, Yang Z, Rute E, Taylor N, Olivier M, and Cowley AW Jr. Gene expression reveals vulnerability to oxidative stress and interstitial fibrosis of the renal outer medulla to non-hypertensive elevations of angiotensin II. Am J Physiol Regul Integr Comp Physiol (January 23, 2003). 10.1152/ajpregu.00257.2002.