Differential regulation of gene expression by ovariectomy in mouse aorta

Amparo C. Villablanca1, Kristine A. Lewis1, Doris Tham1,2 and John C. Rutledge2

1 Division of Cardiovascular Medicine, University of California, Davis, California 95616-8636
2 Division of Endocrinology, Clinical Nutrition, and Vascular Medicine, University of California, Davis, California 95616-8636


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 References
 
The purpose of this study was to investigate the effects of ovarian hormones on gene expression in the vascular wall. Our approach employed an RT-PCR-based cloning strategy of DNA differential display analysis and verification/confirmation of differential expression by semi-quantitative PCR and real-time PCR. mRNA analysis of normal aortas from intact and ovariectomized female C57BL/6J mice, showed altered expression of 20 genes with significant (>70%) sequence homology to known genes. Eight were selected for further study based on the genes’ known function and potential relevance to vascular physiology. Differential expression of mRNA for three genes was confirmed by both semi-quantitative and real-time RT-PCR using gene-specific primers. Ovariectomy downregulated expression of elongation factor-1{alpha} (3.5-fold), ganglioside-induced differentiation associated protein (8.2-fold), and NADH:ubiquinone oxidoreductase (3.8-fold). Thus, in normal mouse aortas, ovariectomy resulted in significant differential downregulation of a number of vascular genes important to vascular cell growth and angiogenesis, cellular differentiation, and mitochondrial energy metabolism, respectively. These studies have implications for our understanding of hormonal regulation of vascular gene expression and the therapeutic targeting of specific vascular genetic sequences by female sex steroid hormones.

vascular; differential display; estrogen; hormones


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 References
 
FEMALE SEX STEROID HORMONES have complex and important vasoactive functions that are just now becoming better elucidated. Although the ovaries make a number of androgens (including androstenedione), ovarian steroids, and protein hormones, the actions of the primary ovarian hormones estradiol and estrone have particular cardiovascular significance and have been most extensively characterized. Estrogen is an important female sex steroid hormone with key actions on both physiological and pathophysiological conditions including reproduction, neoplasia, and atherosclerosis. Specifically, increased cardiovascular disease risk, linked to deficiency of female sex steroid hormones, suggests that they may have important vascular actions.

Although the mechanisms by which ovarian steroids affect vascular function are not clear, ovarian hormones are becoming increasingly recognized as complex and powerful cardiovascular mediators (4, 9, 35, 11). Possible mechanisms by which ovarian hormones, principally estrogen, act on vascular tissue include enhanced vasodilatation via increased production of nitric oxide, a potent vasodilator (7, 36), changes in plasma concentrations of lipids and reductions in vascular cholesterol accumulation (7, 12), and prevention of apoptotic cell death (2, 30).

In addition, as recently reported (13), changes in vascular gene expression also occur in response to estrogen. Vasoactive genes known to be estrogen-regulated include inflammatory mediators such as interleukins IL-4 and IL-6 (3, 8) and adhesion molecules VCAM-1 and ICAM-1 (5, 2628). In addition, we have demonstrated that the substance P receptor is upregulated by estrogen (32). More recently, using differential display mRNA gene expression profiling of human aortic endothelial cells, our laboratory has reported previously uncharacterized estradiol-sensitive dose- and time-dependent upregulation of three additional genes with potential importance to vascular function: aldose reductase, a caspase homolog, and plasminogen activator inhibitor-1 (33).

In this report, we extend our recent observations in cultured aortic endothelial cells to examine the effect of loss of ovarian hormones on the genetic profile of the entire vascular wall. Data from a variety of experimental systems strongly indicate that the loss of ovarian estrogens, as occurs in postmenopausal women and following ovariectomy, is a potent signal for regulation of gene expression (4, 9, 12, 35). It is therefore highly probable that ovariectomy is an essential trigger for significant genetic reprogramming of the vessel wall. However, this remains a relatively uninvestigated area. The current studies were undertaken to establish an initial characterization of changes in gene expression patterns in vessels in response to loss of ovarian sex hormones and to identify specific vascular genetic targets for ovarian hormones.

Based on our previous work, we reasoned that the loss of ovarian hormones would be coupled to a significant genetic reprogramming in the vascular wall of aortas. Experiments were therefore designed using the technique of mRNA gene display (14) to identify differentially expressed genes in response to ovariectomy in normal mouse aortas. This strategy enabled us to identify genes that were selectively induced and/or downregulated by ovariectomy and to establish relative changes in gene expression levels compared with intact mice. This knowledge is important to advance our current understanding of the profile of genes under hormonal regulation in the vasculature.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 References
 
Materials.
For differential display, materials were obtained from the following sources: TRI Reagent (Molecular Research, Cincinnati, OH); MessageClean DNase I (GenHunter, Nashville, TN); phenol/CHCl3 (3:1) (Sigma, St. Louis, MO); RNAimage mRNA differential display system fingerprinting (GenHunter); Sequagel-6 (National Diagnostics, Atlanta, GA); [{alpha}-33P]dATP Easytide(2,000 Ci/mmol; New England Nuclear, Boston, MA); anchor and arbitrary oligonucleotide primers (Molecular Probes, Eugene, OR), oligo-dT primer (Life Technologies-Invitrogen, Carlsbad, CA); and 1:5,000 SYBR Gold nucleic acid prestain (Molecular Probes).

For reamplification, cloning and sequencing, materials, and reagents were obtained from the following sources: AmpliTaq DNA polymerase (Perkin-Elmer, Foster City, CA); T7 promoter and M13 reverse primers (Midland Certified Reagent, Midland, TX); Invitrogen TOPO TA cloning kit, containing 10x ligation buffer, T4 DNA ligase, SOC medium, INV{alpha}F' E. coli competent cells, and pCR2.1-TOPO vector (Invitrogen, Carlsbad, CA); Luria-Bertani (LB) medium and agar (Media Room, UC Davis); ampicillin (Apothecon, Princeton, NJ); 5-bromo-4-chloro-3-indolyl-ß-D-galactopyranoside (X-Gal; Boehringer Mannheim, Indianapolis, IN); KH2PO4 and K2HPO4 (Fisher Scientific, Fairlawn, NJ); and Terrific Broth containing tryptone, yeast extract, and glycerol (Difco, Detroit, MI).

Strategy, overview, and advantages of the experimental approach used to identify and verify differential expression.
To identify genes whose expression was differentially regulated by ovariectomy in mouse aortas, we bilaterally ovariectomized normal postpubescent C57BL/6J mice, removed their aortas, and compared gene expression to aortas from C57BL/6J mice with intact ovaries. Total RNA was extracted from the aortas. Three 1-base anchored oligo-dT primers with 5' HindIII sites, used in combination with a series of 16 additional arbitrary 13-mer oligonucleotide upstream primers (also containing 5' HindIII sites), were used to reverse transcribe and amplify the mRNAs from the control and ovariectomized aorta samples. Following amplification, the mRNA subpopulations were resolved on DNA sequencing gels. Hundreds of bands representing expressed gene fragments were screened for differential expression in aortas from intact and ovariectomized mice. DNA from excised bands was eluted, directly sequenced, and queried against National Center for Biotechnology Information (NCBI) databases using the Basic Local Alignment Search Tool (BLAST) algorithm. A match of interest was defined as having >70% homology to a known gene. Following database verification, genes with potential physiological relevance to vascular function were selected for further verification. To confirm differential expression of apparently differentially displayed cDNAs, gene-specific primers were designed and used to amplify the reverse transcription products, which were subsequently analyzed by both semi-quantitative and real-time RT-PCR.

The technique of delta RT-PCR used in these studies is based on reverse transcription and amplification of 3'-terminal RNA sequences (3, 14). The method relies on a combination of thermophilic DNA polymerases with 3'-5' exonuclease (proofreading) activity to produce larger, easier to detect, higher fidelity, and not overcycled PCR products than is possible with conventional PCR (20). These studies were also designed to minimize false-positive results and optimize sensitivity. This technology ensures that differentially expressed genes are more readily identified and cloned than with prior generations of differential display. Specifically, the use of three 1-base anchored oligo-dT primers to subdivide the mRNA population reduces the redundancy and under-representation of the subpopulation of mRNAs (14). In addition, with built-in restriction sites at the 5' ends of both anchored and arbitrary 13-mers, the longer primer pairs produce highly selective and reproducible cDNA patterns (34).

Ovariectomy and isolation of mouse aorta.
Eight female C57BL/6J mice were obtained from Charles River (San Diego, CA). Four mice were bilaterally ovariectomized at 6–7 wk of age and the other four mice left intact. Both sets of mice were kept in as near identical housing conditions as possible controlling for the same environment, food, light, and temperature conditions. The ovariectomized mice were allowed to recover for a week, and both groups were killed 5–6 wk later using pentobarbital anesthesia (50 mg/kg). Under RNase-free conditions, the aortas were immediately removed by micro-excision of the entire aorta from the heart to the iliac bifurcation. The aortas were then trimmed of adherent fat, rinsed clean of blood in diethyl pyrocarbonate (DEPC)-treated PBS, and stored at -80°C for subsequent extraction of total RNA.

RNA isolation.
The aortas of ovariectomized and intact mice were homogenized in a tissue homogenizer, and total RNA was isolated using TRI-Reagent (Molecular Research, Cincinnati, OH), per the manufacturer’s instructions. Each aorta yielded ~14–16 µg of total RNA, the predicted yield per the manufacturer’s estimates. Total RNA was then treated with DNase I for removal of contaminating genomic DNA. All lab ware and solutions used to isolate RNA were rendered RNase-free prior to use. RNA samples were stored in DEPC-treated water at -80°C until ready for use.

Serum estradiol assay.
To confirm that following ovariectomy levels of circulating estradiol were low as expected, serum estradiol levels were determined in ovariectomized mice and compared with those of intact control mice using an estradiol enzyme immunoassay (EIA) kit (Cayman Chemical, Ann Arbor, MI) following the manufacturer’s instructions. Using this method, we found plasma estradiol levels were 168 pg/ml and 23 pg/ml, in intact and ovariectomized mice, respectively.

Differential mRNA display analysis of mouse aortas.
Differential display was performed using an RNAimage mRNA differential display fingerprinting system from GenHunter. Eight samples of total RNA (0.2 µg each), prepared individually from each of the aortas of the four intact and the four ovariectomized mice, were used as template for reverse transcription to cDNA prior to performing differential display. RNA samples were not pooled for any of the analyses. Fibroblast RNA provided with the GenHunter kit, and H2O samples minus the reverse transcriptase, were used for positive and negative RNA controls, respectively. To subdivide the mRNA population for each RNA sample into three subpopulations, total RNA was first converted into first-strand cDNA with MMLV reverse transcriptase, and each of three different 1-base anchored H-T11M oligo-dT primers [where H = the 5' end restriction site (HindIII), T11= an 11-bp long poly-T tail, and M = G, C, or A].

For fingerprinting, gels were hybridized with 2.0 µCi {alpha}-33P-labeled cDNA (3.5 ng/sample) prepared by reverse transcription of total tissue RNA. PCR amplification was performed with DNA polymerase and the same H-T11M oligo-dT anchor primers and "H-AP" arbitrary primers used for the reverse transcription reactions; see Table 1 for primer sequences. A total of 32 primer pairs were used to amplify the 8 aortic tissue samples. Each primer pair would be expected to yield a different fingerprint and identify a different set of differentially expressed RNAs, typically 1–2 differential bands per pair. Based on previously published axioms (15), the 16 arbitrary primers used, combined with each of the three 1-base anchored primers, would be expected to detected between 56% and 64% of the total aortic mRNAs.


View this table:
[in this window]
[in a new window]
 
Table 1. Anchor oligo-dT and arbitrary oligonucleotide primers used for RT-PCR and differential display amplification reactions

 
Multiple sequencing gels were prepared and electrophoresed with several replicates of samples from intact and ovariectomized mice. Each replicate consisted of two different DNA samples from intact aortas (control 1 and control 2) loaded adjacent to two different cDNA samples from ovariectomized aortas (OVX1 and OVX2). This greatly facilitated subsequent screening of the gels and optimized the likelihood of correctly selecting bands that appeared to be differentially displayed by ovariectomy. cDNA sequences were reamplified in a thermal cycler (Perkin-Elmer 2400) based on chance homology to "H-AP" primers following the GenHunter RNAimage protocol using standard PCR techniques (21). The amplified cDNA subpopulation of mRNAs distributed on DNA sequencing gels was exposed overnight to Kodak Biomax-MR film.

Identification, recovery, and reamplification of apparently differentially displayed cDNAs.
Bands demonstrating apparent differential expression by analysis of gel radiographs were identified by two independent observers. Differential display was suggested by the presence of bands in both samples from control aortas and absence in both samples from ovariectomized aortas, or vice versa. An example is provided in Fig. 1. Thus only transcripts that reproducibly differentially displayed were selected for further analysis. Bands that appeared to be differentially expressed were excised from the gels, then cDNA was eluted and reamplified by PCR using the fingerprinting H-T11M and H-AP primers. Bands were directly sequenced and cloned, and the clones were resequenced.



View larger version (39K):
[in this window]
[in a new window]
 
Fig. 1. Representative example of autoradiograph from mRNA differential display analysis. This is a portion of a larger polyacrylamide sequencing gel showing an example of a differentially expressed band. The gel was loaded with side-by-side cDNA samples that had been reverse transcribed from total RNA from each of two different aortas from intact mice (control, C1 and C2) compared with side-by-side samples from each of two different aortas from ovariectomized mice (OVX1 and OVX2). Arrow: band (34c) is differentially downregulated by OVX (present in both control intact aorta samples but absent in both OVX aorta samples).

 
Cloning and sequencing of cDNAs.
For cloning, the reamplified PCR products were ligated (pCR2.1-TOPO vector; Invitrogen) following transformation of One-Shot TOP10F E. coli competent cells according to the manufacturer’s protocol using standard molecular cloning techniques (25). The PCR product insert was reamplified using T7 promoter and M13 reverse primers designed to match the flanking pCR2.1-TOPO vector sequences to the PCR insert. The thermal cycling parameters were: thermal activation for 8.5 min at 94°C, and 25 cycles of PCR (melting for 30 s at 94°C, and annealing/extension for 30 s at 50°C, 1 min 45 s at 72°C, followed by an additional 10 min at 72°C). The correct insert was confirmed by electrophoresis through a 1% agarose gel stained with ethidium bromide, and the product was then purified (Centricon-100 concentrator). Differentially displayed bands that were successfully amplified and cloned were subsequently sequenced on a 4.25% acryl/bisacrylamide (5% LongRanger) gel at the University of California, Davis Sequencing Facility (ABI Prism 377 DNA Sequencer) using an ABI Prism Dye Terminator Cycle Sequencing Ready Reaction kit with AmpliTaq DNA polymerase FS, per the manufacturer’s instructions. Data was analyzed using ABI Prism Sequencing 2.1.1 software.

Sequence data was matched against human and mouse genome databases, including GenBank (http://www.ncbi.nlm.nih.gov/Education/index.html) and the Jackson Laboratory site (http://www.jax.org), for nucleotide sequences, and against protein and expressed sequence tag (EST) databases using FASTA (FastA) and NetBLAST based on heuristic algorithms (1, 22). In general, bands that produced direct sequence information and had high homology (>70%) to known human mRNAs were chosen and pursued for subsequent verification of differential expression.

Semi-quantitative RT-PCR-based verification analysis of differential gene expression.
For the initial verification analysis we used semi-quantitative RT-PCR, as total RNA from mouse aorta (~14.5 µg total RNA/aorta) was insufficient for poly(A)+ selection and detection by conventional Northern analysis. Verification was performed with gene-specific primers (29, 33), as low levels of contaminating cDNAs in probes prepared directly from differentially displayed bands can obscure accurate verification of bands (18). The online program Primer3 (24) http://www.genome.wi.mit.edu/cgi-bin/primer/primer3.cgi was used to design specific primers for each of the genes of interests and D-glyceraldehyde-3-phosphate dehydrogenase (GAPDH), a "housekeeping" gene. Primer pairs yielded PCR products with the exact length predicted by the known sequences of the five genes (3b, 5a, 18f, 22a, and 34c) selected for verification analysis (Table 2).


View this table:
[in this window]
[in a new window]
 
Table 2. Specific oligonucleotide primers used in this study for semi-quantitative PCR verification and real-time PCR confirmation of differential gene expression

 
For verification, reverse transcription of total RNA (2 µg each) from two of the intact and ovariectomized mouse aortas used in the differential display analysis was performed to yield cDNAs. Oligo-dT was used as primer in these reactions. Not subdividing the mRNA population with H-T11M oligo-dT permitted use of the same housekeeping gene primers with each of the cDNAs of interest. This was important to equalize the initial total cDNA template concentrations by PCR product signal, thus permitting direct comparison of signals for semi-quantitative analysis. In addition, several preliminary template equalization experiments were performed utilizing GAPDH primers to amplify the cDNA templates. Once the initial starting concentrations of intact and ovariectomized mouse aorta cDNAs were equalized, all subsequent cDNA template dilutions resulted in the same densitometry values for each corresponding dilution for GAPDH.

A titration analysis curve was then constructed for the five genes of interest and GAPDH by making a dilution series of cDNA template (1:10, 1:20, 1:40, and 1:80, corresponding to approximately 3.5, 1.8, 0.9, and 0.4 ng of cDNA, respectively). cDNAs were then PCR amplified, and the band intensity signals quantified by densitometry to verify that relative differences in gene expression remained constant over a wide template concentration range. The optimal conditions for amplification (the relative proportion of upstream to downstream primers, temperatures, and cycle number) were experimentally determined. To reduce spurious priming during gene amplification, PCR was performed for 10 "touchdown" cycles (6) as follows: The first "touchdown" cycle was for 30 s at 94°C, 30 s at 61°C for annealing, and 2 min at 72°C. The next "touchdown" cycle was for 30 s at 94°C, 30 s at 60.9°C for annealing, and 2 min at 72°C. For the remaining 8 cycles the annealing temperature dropped by 0.1°C after each cycle. These touchdown cycles were followed by 10–15 cycles of 30 s at 94°C, 30 s at 60°C, and 2 min at 72°C. Samples were held for a final 7-min extension cycle at 72°C, then cooled to 4°C. Correct PCR product was verified by electrophoresis through a 1% agarose gel stained with 1:5,000 SYBR Gold nucleic acid prestain.

DNA bands were visualized with a STORM Fluorimager (Molecular Dynamics, Sunnyvale, CA) and quantified by spot densitometry using STORM Fluorimager ImageQuant software. The mean densitometry signal (calculated as an average of the signal at each template dilution) was expressed as a ratio (intact: ovariectomized) for each of the five genes of interest. The 1:80 template dilution data was not used for this final analysis given difficulty in obtaining accurate readings at this very low dilution. Because there was equalization of target cDNA by GAPDH, it was not necessary to further normalize signals to GAPDH prior to directly comparing intact and ovariectomized sample signals.

Real-time RT-PCR-based confirmation analysis of differential gene expression.
Verification of differential gene expression by semi-quantitative PCR with gene-specific primers is an established approach for analysis of differential expression (18). However, the semi-quantitative RT-PCR verification experiments used the same template RNA as was used for differential display analysis. Therefore, to further confirm differential mRNA expression, a more quantitative approach, real-time PCR, was employed. For these studies, two additional mouse aortas (one ovariectomized and one intact) were analyzed. These aortas were not used in the previous differential expression and semi-quantitative PCR verification experiments. However, to mitigate against the possible influence of environmental factors on gene expression, care was taken to ensure that the additional aortas were from the same experimental animal cohort as that described earlier in MATERIALS AND METHODS.

Briefly, total RNA from the ovariectomized and intact aortas was treated with DNase (MessageClean DNase I; GenHunter) and cDNA generated by RT-PCR. A mixture of oligo-dT oligodeoxynucleotide primers (T12–18) and random hexamer primers were used for these studies following the protocol for SuperScript RNase H- reverse transcriptase (GIBCO-BRL, Rockville, MD). Real-time detection of PCR was performed using the GeneAmp 5700 sequence detection system (Applied Biosystems, Foster City, CA) according to the manufacturer’s instructions. The same gene-specific primers used for semi-quantitative PCR analysis were used for real-time PCR detection of the five genes of interest (3b, 5a, 18f, 22a, and 34c). Equal amounts of sample cDNA, used in duplicate or triplicate, were amplified with the SYBR Green I Master Mix (Applied Biosystems). The thermal cycling parameters were: thermal activation for 10 min at 95°C followed by 40 cycles of PCR (melting for 15 s at 95°C and annealing extension for 1 min at 60°C). A standard curve was constructed using a template dilution series (1:10, 1:20, 1:40, 1:80, and 1:160) of total RNA from mouse aorta. A "no template control" was included with each PCR. Amplification efficiency for each of the five genes of interest was then validated and normalized against GAPDH. To compare the relative expression of the various genes, data was expressed as the ratio of intact:ovariectomized aortas.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 References
 
Differential display of mRNAs from ovariectomized and intact mouse aortas.
The differential display autoradiographs yielded over 2,000 gene transcripts that were carefully screened for differential display by two independent observers (see MATERIALS AND METHODS). A total of 38 cDNA fragments appeared to be differentially displayed by ovariectomy and were subsequently isolated from the gels. Of these, 32 fragments were successfully cloned and sequenced. Fifteen transcripts were apparently upregulated by ovariectomy (present in aortas of ovariectomized mice but absent in aortas of control intact mice), and 17 transcripts were apparently downregulated by ovariectomy (present in the aortas of intact mice but absent in aortas of ovariectomized mice). Because differential display fingerprinting is not a quantitative tool, but is best utilized as an absent or present screening, it is not possible to categorically state that gene expression for the downregulated genes was zero following ovariectomy. However, when autoradiographs were viewed by two independent observers, there was no detectable gene expression as demonstrated by complete absence of a visual band.

Following database verification in GenBank, 20 of the apparently differentially displayed bands were found to have significant (greater than or equal to 70%) homology to known genes. These were considered for further analysis (Table 3). The remaining 12 genes had lesser degrees of homology to known genes and were therefore not considered for further study. As anticipated from GenBank searches, none of the apparently differentially expressed clones represented novel genes. However, several had short and poor matches to known genes or matches to ESTs of unknown function, potentially representing novel genes.


View this table:
[in this window]
[in a new window]
 
Table 3. Summary of OVX-regulated clones in mouse aorta with significant homology to know genes and of genes selected for verification of differential expression

 
Selection of genes for verification of differential expression.
Two main criteria were used to select genes for further study and verification: 1) whether the gene was of interest based on its known function, and 2) whether the gene’s function had potential physiological relevance to cardiovascular function (e.g., vascular cell growth, vascular permeability, vascular inflammation, etc.). This determination was made after a careful evaluation of the literature for the reported function of each of the 20 significantly homologous genes. Based on these criteria, we opted to not pursue further 12 of the 20 apparently differentially expressed genes. Of the remaining eight genes selected for final verification (3b, 5a, 18f, 18d, 21da, 22a, 32ab, and 34c), two were upregulated by ovariectomy (clones 21da and 32ab), and six were downregulated by ovariectomy (clones 3b, 5a, 18f, 18d, 22a, and 34c).

Verification of the effect of ovariectomy on aortic differential gene expression by semi-quantitative RT-PCR analysis.
Semi-quantitative RT-PCR identified significant (>2-fold greater than control) differential expression of mRNAs for five of the eight genes of interest. Coincidentally, all of the five genes were downregulated by ovariectomy (Table 3). They had cDNA sequences encoding for: Tat binding protein 1 (TBP-1; clone 3b), a heparin-binding angiogenic growth factor; skeletal muscle {alpha}-tropomyosin (clone 5a), an actin-binding protein central to the control of calcium-regulated striated muscle contraction; elongation factor 1{alpha} (EF-1{alpha}; clone 18f), a polypeptide involved in vascular endothelial cell sprouting and elongation; and ganglioside-induced differentiation associated protein (clone 22a), a gene involved in signal transduction to induce cellular differentiation. Moreover, the gene for mitochondrial nicotinamide adenine dinucleotide (NADH)-ubiquinone oxidoreductase chain 6 (clone 34c), an important mitochondrial electron transport enzyme, was also downregulated by ovariectomy. One of the downregulated genes of interest (clone 18d) and two of the upregulated genes of interest (clones 21da and 32ab) demonstrated less than twofold greater expression than GAPDH and therefore were not significantly differential following ovariectomy.

To assess relative changes in differential gene expression among the five genes verified to be differentially expressed by semi-quantitative RT-PCR, direct comparisons were made between densitometry signals for intact and ovariectomized aortas. This was possible, as the cDNA PCR template between aortic samples of intact and ovariectomized mice, at four serial template dilutions, had been equalized between GAPDH and each of the five genes of interest; see Fig. 2A for an example. The semi-quantitative RT-PCR analysis of the mRNA densitometry ratios (intact:ovariectomized aorta) for GAPDH and the five genes of interest is shown in Fig. 2, B and C. As can be noted, the five differentially expressed genes demonstrated significantly greater mRNA expression in intact compared with ovariectomized mouse aortas, consistent with downregulation of gene expression by ovariectomy. The relative differences in differential expression levels between intact and ovariectomized aortas were as follows: 3.3-fold for TBP-1 (clone 3b), 3.2-fold for skeletal muscle {alpha}-tropomyosin (clone 5a), 3.5-fold for EF-1{alpha} (clone 18f), 8.2-fold for ganglioside-induced differentiation associated protein (clone 22a), and 3.8-fold for NADH:ubiquinone oxidoreductase (clone 34c).



View larger version (40K):
[in this window]
[in a new window]
 
Fig. 2. Semi-quantitative RT-PCR verification of mRNA differential gene expression. A: template equalization. Reverse-transcribed cDNA from intact and ovariectomized (OVX) mouse aortas for each of the specific genes of interest (clone 3b, clone 5a, clone 18f, clone 22a, and clone 34c; see text) were co-amplified with gene-specific primers and primers for a housekeeping gene, GAPDH. Serial sequential 1:2 dilutions of cDNA template for the specific genes and for GAPDH were used to verify that relative differences in gene expression remained constant over a wide template concentration range. The graph demonstrates equalized cDNA PCR template from intact and OVX mouse aortas for GAPDH at four cDNA template dilutions (1:10, 1:20, 1:40, and 1:80). Similar experiments were done for each of the five specific genes of interest with similar linear curves obtained (data not shown). B: verification radiography. cDNA products from the intact and OVX aortas were PCR amplified with gene-specific primers for the genes of interest (clones 3b, 5a, 18f, 22a, and 34c). The PCR products were gel electrophoresed and band intensity, at each of four cDNA template dilutions (1:10, 1:20, 1:30, and 1:40), quantified by densitometry. The gels demonstrate greater signal intensity for the genes of interest in aortas from intact vs. OVX mice at the four cDNA template dilutions (compare four bands from intact aortas on left of gels to four bands from OVX aortas on right of gels), consistent with downregulation of gene expression by OVX. As expected, gene expression for GAPDH, the housekeeping gene control, was unaffected by OVX treatment at each of the four cDNA template dilutions. C: densitometry ratios. Densitometry ratios of differentially expressed bands for intact and OVX samples were calculated from the individual band densitometry signals at each of four cDNA template dilutions (1:10, 1:20, 1:40 and 1:80). The mean densitometry ratio for the 1:10, 1:20, and 1:40 cDNA template dilutions was then calculated for each gene of interest and GAPDH. The 1:80 template dilution data was not used for this final analysis, given significant variability and difficulty in obtaining adequate readings at this very low dilution. For each of the five specific genes of interest (clones 3a, 5a, 18f, 22a, and 34c), the mean densitometry signal ratio of intact:OVX was significantly greater (>2.0) in intact compared with OVX aortas, confirming downregulation of expression by OVX.

 
Confirmation of the effect of ovariectomy on differential aortic gene expression by real-time RT-PCR analysis.
Table 4 summarizes the mRNA expression ratios for intact:ovariectomized aortas using real-time PCR for confirmation of differential expression. They were as follows: 316-fold for EF-1{alpha} (clone 18f), 2.2-fold for ganglioside-induced differentiation associated protein (clone 22a), and 18,000-fold for NADH:ubiquinone oxidoreductase (clone 34c). The fold change for NADH:ubiquinone oxidoreductase (clone 34c) was unusually large in that this gene had no detectable message in aortas following ovariectomy compared with robust expression in intact aorta. In contrast to the verification studies using semi-quantitative PCR, real-time PCR confirmed differential expression for bands 18f, 22a, and 34c. Clones 3b and 5a were determined to not be differential in this analysis, suggesting that real-time PCR may have identified these exceptions as false positives. As expected, the gene tag for GAPDH control showed nondifferential expression. Thus differential expression for the mRNA for three genes, EF-1{alpha}, ganglioside-induced differentiation associated protein, and NADH:ubiquinone oxidoreductase, was confirmed to be significantly downregulated by both semi-quantitative and real-time PCR-based verification.


View this table:
[in this window]
[in a new window]
 
Table 4. Confirmation of differential gene expression by real-time PCR in OVX mouse aortas

 

    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 References
 
To gain insight into the actions of ovarian hormones in genetic regulation of the vascular wall, we performed mRNA differential display analysis to identify genes differentially regulated by ovarian sex steroids in murine aortas. In so doing, we have demonstrated the utility of differential display as a molecular screen to profile genes that are differentially expressed by loss of ovarian steroids in vascular tissue, and we identified vascular genes whose expression was not previously known to be significantly modulated by ovariectomy. In the systematic screening performed in this study, 20 genes with significant sequence homology to known genes demonstrated apparent differential display following ovariectomy. Eight genes were selected for verification of differential expression based on the gene’s physiological relevance to cardiovascular function. Three genes were subsequently verified and confirmed to have significant differential expression in aortas; all were downregulated by ovariectomy compared with the intact animals: elongation factor (EF)-1{alpha}, ganglioside-induced differentiation associated protein, and NADH:ubiquinone oxidoreductase. These genes encode a host of proteins with seemingly divergent functions, suggesting a link between loss of ovarian function and angiogenesis, cell growth and differentiation, and mitochondrial energy metabolism, respectively.

The EF-1{alpha} family of polypeptides, one of the differentially displayed genes identified to be downregulated in our study, promote polypeptide chain elongation during mRNA translation. EF-1{alpha} is a ubiquitously expressed important component of the protein synthesis apparatus in eukaryotic tissues (10). EF-1{alpha} is also involved in vascular endothelial cell sprouting and elongation in vitro, as well as the formation of networks of elongated endothelial cells (17). Therefore, downregulation of EF-1{alpha} by ovariectomy could be important in retarding angiogenesis.

The second gene identified to be downregulated by ovariectomy in this study was the gene for ganglioside-induced differentiation associated protein. This gene was originally identified in a search of genes involved in ganglioside-induced differentiation of a mouse neuroblastoma cell line, where it causes cell differentiation with neurite sprouting (16). The gene is expressed in most tissues and might be involved in signal transduction to induce cellular differentiation. However, it is unknown whether it could potentially have a role in the cellular differentiation that accompanies endothelial cell sprouting in angiogenesis.

The third gene downregulated by ovariectomy was mitochondrial NADH:ubiquinone oxidoreductase chain 6, an important mitochondrial electron transport enzyme. This gene regulates the redox state of the carriers of the electron-transport chain of mitochondria that correlate with changes in a number of physiological parameters. It is proposed that NADH:ubiquinone oxidoreductase inhibitors block multiple and possibly reactive oxygen species-modulated pathways (23). In addition, cytokines, such as tumor necrosis factor-{alpha} (TNF{alpha}) and IL-1{alpha}, inhibit both pyruvate dehydrogenase activity and mitochondrial function in tumor and vascular smooth muscle cells (37). If the same occurs with ovariectomy, then ovariectomy may have the potential for primary impairment of mitochondrial function and the cellular redox state.

In nonvascular animal systems ovariectomy has been previously associated with changes in the expression of genes with primary relevance to hormonal actions. For example, in rats ovariectomy is associated with a progressive increase in mRNA levels for {alpha}-2,3-sialyltransferas ({alpha}-2,3-STase), one of the enzymes that incorporates sialic acid residues into follicle stimulating hormone (FSH) and an important regulator of FSH isoforms and estrogen output (4). However, ovariectomy has also been demonstrated to be associated with changes in expression of genes pertinent to vascular physiology. For example, mRNA expression for steroid hormone receptors was demonstrated to decrease after ovariectomy in aortic and venous rat vessels (11). Furthermore, in ovariectomized female pigs, expression of mRNA for both the endothelial isoform of nitric oxide synthase (eNOS) and endothelin-1, two important vasoactive peptides, is increased compared with their expression in intact female pigs (31). Interestingly, in humans, expression of the P450c5 gene, the key enzyme involved in the conversion of progestins to androgens, demonstrates significant variation in levels of expression in postmenopausal ovaries compared with normal cycling ovaries, without significant differences in the overall gene expression levels between the two states (9). Therefore, the assumption that loss of ovarian hormones inevitably results in changes in gene expression in all systems needs to be challenged and specifically tested.

The most plausible mechanism by which ovariectomy could regulate gene expression in the blood vessel wall is by ovarian steroid-stimulated regulation of mRNA transcription. However, in other systems, changes in gene expression may be a consequence of mRNA breakdown (19). Therefore, it is possible that ovariectomy could affect both the transcriptional rate of the genes identified and/or the stability of their mRNAs. Although our study was not designed to tease out these processes, the data presented provide significant new information about vascular targets for hormone action and lay the foundation for additional future mechanistic work in this area.

In this study, only genes whose sequences are found in GenBank were investigated, yielding information on regulation of known genes by loss of ovarian steroids. It would be interesting as a next step to embark on the discovery of unknown genes regulated by ovariectomy. Novel genes may be represented by EST fragments that are in GenBank, yet have no known function or characterization. In addition, poor matches to GenBank nucleotide sequences also could be novel genes. In fact, some of the ovariectomy-regulated sequences found by this study corresponded to ESTs in the human genome database with unknown function or poor matches to known genes; perhaps some of these ESTs are to novel genes. Although those ESTs were not further characterized, cloning them would be of future interest.

The bipartite differential display technique utilized for these studies, consisting of 1) RT-PCR and electrophoresis followed by 2) gene identification and verification, takes advantage of new innovations designed to improve reliability and reduce false positives in differential display. For example, the methodology improves on previous approaches to differential display by subdividing the mRNA population in three, built in restriction sites at the 5' ends of both anchored and arbitrary primers, longer PCR primers, loading of samples in duplicate, and the use of two approaches to verify and confirm differentially expressed cDNAs. However, our study has certain limitations. The techniques, although reliable in identifying differentially expressed gene products, do not provide a functional assessment of the differentially expressed genes identified. It is also possible that some of the lower sequence homology genes (<80%) may represent mismatches. Furthermore, ovariectomy could have varying effects on vascular gene expression at differing times following ovariectomy and at different phases of the reproductive cycle. This would not be revealed by the present studies as we studied mice at one point in time. An additional limitation is that loss of ovarian steroids, hormones, and peptides following ovariectomy does not permit inferences to be made about the role of specific ovarian hormones (e.g., estrogen) or ovarian vasoactive peptides (e.g., vascular endothelial growth factor and relaxin) on regulation of gene expression. Rather, our studies are an important first step toward this analysis. Lastly, our studies do not localize gene expression changes in the vessel wall. Given that the aorta is a complex tissue consisting of several cell types (endothelial cells, smooth muscle cells, and adventitia), gene expression could clearly vary by cell type relative to cell function. Nonetheless, our studies provide novel data in this evolving area of science, lay new groundwork for further research in this area, and add to our understanding of the molecular basis for female sex steroid hormone action in the vasculature.

In summary, the results of our studies serve to extend the findings of prior studies in three ways. First, we have demonstrated that ovariectomy is associated with significant changes in vascular gene expression in the vessel wall. Second, we have identified three vascular genes previously unreported to be regulated by ovarian hormones that may act as genetic targets for hormone action. And third, in our analysis we have broadened the spectrum of gene functions potentially regulated by ovarian steroids in the vasculature.

In conclusion, the genes differentially regulated by ovariectomy identified in the present study have not previously been identified to be under the control of ovarian hormones. These differentially regulated genes encode for a wide array of cellular regulatory functions. Based on the function of the differentially expressed genes identified, loss of ovarian hormones, principally but not exclusively estrogens, may have relevance to a wide spectrum of regulatory factors in the vasculature related to mitochondrial energy metabolism, cellular differentiation, and the cellular growth and vascular sprouting that characterizes new vessel formation. Thus, understanding the genetic vascular changes accompanying loss of ovarian steroids may enhance our understanding of hormonal regulation of vascular function in health and disease.


    ACKNOWLEDGMENTS
 
This work was supported by National Institutes of Health Career Development Award KO1-HL-04142 (to A. C. Villablanca) and Grant R01-HL-55667 (J. C. Rutledge). We also are indebted to the Richard A. and Nora Eccles Harrison Endowed Chair in Diabetes Research for support of this project.


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

Address for reprint requests and other correspondence: A. C. Villablanca, Division of Cardiovascular Medicine, Univ. of California, Davis, One Shields Ave., TB 172, Davis, CA 95616-8636 (E-mail: avillablanca{at}ucdavis.edu).

10.1152/physiolgenomics.00040.2002.


    References
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 References
 

  1. Altschul SF, Gish W, Miller W, Myers EW, and Lipman DJ. Basic local alignment search tool. J Mol Biol 215: 403–410, 1990.[ISI][Medline]
  2. Alvarez RJ, Gips SJ, Moldovan N, Wilhide CC, Milliken EE, Hoang AT, Hruban RH, Silverman HS, Dang CV, and Goldschmidt-Clermont PJ. 17ß-Estradiol inhibits apoptosis of endothelial cells. Biochem Biophys Res Commun 197: 43–49, 1997.
  3. Bauer D, Warthoe P, Rohde M, and Struss M. PCR Methods Applications Manual Supplement. Cold Spring Harbor, NY: Cold Spring Harbor Laboratory, 1994.
  4. Damián-Matsumura P, Zaga V, Maldonado A, Sánchez-Hernández C, Timossi C, and Ulloa-Aguirre A. Oestrogens regulate pituitary {alpha}2,3-sialyltransferase messenger ribonucleic acid levels in the female rat. J Mol Endocrinol 23: 153–165, 1999.[Abstract/Free Full Text]
  5. Dickens GR, Matheny CJ, Morris PE, Clifton GD, and Ensom MH. A pilot study of estrogen’s effects on bronchial myocyte adhesion molecule expression. Pharmacotherapy 19: 1426–1431, 1999.[ISI][Medline]
  6. Don RH, Cox PT, Wainwright BJ, Baker K, and Mattick JS. "Touchdown" PCR to circumvent spurious priming during gene amplification (Abstract). Nucleic Acids Res 19: 4008, 1991.[ISI][Medline]
  7. Farhat MY, Lavigne MC, and Ramwell PW. The vascular protective effects of estrogen. FASEB J 10: 615–624, 1996.[Abstract/Free Full Text]
  8. Harris MT, Feldberg RS, Lau KM, Lazarus NH, and Cochrane DE. Expression of proinflammatory genes during estrogen-induced inflammation of the rat prostate. Prostate 44: 19–25, 2000.[ISI][Medline]
  9. José M, Puche C, Cabero A, Cabero L, and Meseguer A. Expression of P450c17 messenger ribonucleic acid in postmenopausal human ovary tissues. Fertil Steril 71: 528–535, 1999.[ISI][Medline]
  10. Kahns S, Lund A, Kristensen P, Knudsen CR, Clark BF, Cavallius J, and Merrick WC. The elongation factor 1 A-2 isoform from rabbit: cloning of the cDNA and characterization of the protein. Nucleic Acids Res 26: 1884–1890, 1998.[Abstract/Free Full Text]
  11. Knauthe R, Diel P, Hegele-Hartung C, Engelhaupt A, and Fritzemeier KH. Sexual dimorphism of steroid hormone receptor messenger ribonucleic acid expression and hormonal regulation in rat vascular tissue. Endocrinology 137: 3220–3227, 1996.[Abstract]
  12. Kushwaha RS and Hazzard WR. Exogenous estrogens attenuate dietary hypercholesterolemia and atherosclerosis. Metabolism 30: 359–366, 1981.[ISI][Medline]
  13. Landers JP and Spelsberg TC. New concepts in steroid hormone action: transcription factors, protooncogenes, and the cascade model for steroid regulation of gene expression. Crit Rev Eukaryot Gene Expr 2: 19–63, 1992.[Medline]
  14. Liang P and Pardee AB. Differential display of eukaryotic messenger RNA by means of the polymerase chain reaction. Science 257: 967–971, 1992.[ISI][Medline]
  15. Liang P, Zhu W, Zhang X, Guo Z, O’Connell RP, Averboukh L, Wang F, and Pardee AB. Differential display using one-base anchored oligo-dT primers. Nucleic Acids Res 22: 5763–5764, 1994.[ISI][Medline]
  16. Liu H, Nakagawa T, Kanematsu T, Uchida T, and Tsuji S. Isolation of 10 differentially expressed cDNAs in differentiated Neuro2a cells induced through controlled expression of the GD3 synthase gene. J Neurochem 72: 1781–1790, 1999.[ISI][Medline]
  17. Mantell DJ, Owens PE, Bundred NJ, Mawer EB, and Canfield AE. 1{alpha},25-dihydroxyvitamin D(3) inhibits angiogenesis in vitro and in vivo. Circ Res 87: 214–220, 2000.[Abstract/Free Full Text]
  18. Martin KJKC, O’Hare MJ, Pardee AB, Sager R. Identification and verification of differential display cDNAs using gene-specific primers and hybridization arrays. Biotechniques 24: 1018–1026, 1998.[ISI][Medline]
  19. Nakamura N, Obayashi H, Fujii M, Fukui M, Yoshimori K, Ogata M, Hasegawa G, Shigeta H, Kitagawa Y, Yoshikawa T, Kondo M, Ohta M, Nishimura M, Nishinaka T, and Nishimura CY. Induction of aldose reductase in cultured human microvascular endothelial cells by advanced glycation end products. Free Radic Biol Med 29: 17–25, 2000.[ISI][Medline]
  20. Nelson K, Brannan J, and Kretz K. The fidelity of TaqPlus DNA polymerase in PCR. Strategies Mol Biol 8: 24–25, 1995.
  21. Neulen J, Yan Z, Raczek S, Weindel K, Keck C, Weich HA, Marmé D, and Breckwoldt M. Human chorionic gonadotropin-dependent expression of vascular endothelial growth factor/vascular permeability factor in human granulosa cells: importance in ovarian hyperstimulation syndrome. J Clin Endocrinol Metab 80: 1967–1971, 1995.[Abstract]
  22. Pearson WR and Lipman DJ. Improved tools for biological sequence comparison. Proc Natl Acad Sci USA 85: 2444–2448, 1988.[Abstract]
  23. Rowlands JC and Casida JE. NADH:ubiquinone oxidoreductase inhibitors block induction of ornithine decarboxylase activity in MCF-7 human breast cancer cells. Pharmacol Toxicol 83: 214–219, 1998.[ISI][Medline]
  24. Rozen S and Skaletsky H. Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol 132: 365–386, 2000.[Medline]
  25. Sambrook J, Fritsch EF, and Maniatis T. Molecular Cloning. A Laboratory Manual. Plainview, NY: Cold Spring Harbor Laboratory Press, 1989.
  26. Scarabin PY, Alhenc-Gelas M, Oger E, and Plu-Bureau G. Hormone replacement therapy and circulating ICAM-1 in postmenopausal women: a randomised controlled trial. Thromb Haemost 81: 673–675, 1999.[ISI][Medline]
  27. Simoncini T, De Caterina R, and Genazzani AR. Selective estrogen receptor modulators: different actions on vascular cell adhesion molecule-1 (VCAM-1) expression in human endothelial cells. J Clin Endocrinol Metab 84: 815–818, 1999.[Abstract/Free Full Text]
  28. Simoncini T, Maffei S, Basta G, Barsacchi G, Genazzani AR, Liao JK, and De Caterina R. Estrogens and glucocorticoids inhibit endothelial vascular cell adhesion molecule-1 expression by different transcriptional mechanisms. Circ Res 87: 19–25, 2000.[Abstract/Free Full Text]
  29. Singer-Sam J, Robinson MO, Bellve AR, Simon MI, and Riggs AD. Measurement by quantitative PCR of changes in HPRT, PGK-1, PGK-2, APRT, MTase, and Zfy gene transcripts during mouse spermatogenesis. Nucleic Acids Res 18: 1255–1259, 1990.[Abstract]
  30. Spyridopoulos I, Sullivan AB, Kearney M, Isner JM, and Losordo DW. Estrogen-receptor-mediated inhibition of human endothelial cell apoptosis. Estradiol as a survival factor. Circulation 95: 1505–1514, 1997.[Abstract/Free Full Text]
  31. Ton C, Hwang DM, Dempsey AA, and Liew CC. Identification and primary structure of five human NADH-ubiquinone oxidoreductase subunits. Biochem Biophys Res Commun 241: 589–594, 1997.[ISI][Medline]
  32. Villablanca AC and Hanley MR. 17 B-estradiol stimulates substance P receptor gene expression. Mol Cell Endocrinol 135: 109–117, 1997.[ISI][Medline]
  33. Villablanca AC, Lewis K, and Rutledge JC. Time and dose-dependent differential upregulation of three genes by 17ß-estradiol in human aortic endothelial cells. J Appl Physiol 92: 1064–1073, 2002.[Abstract/Free Full Text]
  34. Vogeli-Lange R, Burckert N, Boller T, and Wiemken A. Screening for positive clones generated by differential display. Methods Mol Biol 85: 95–103, 1997.[Medline]
  35. Wang X, Barber DA, Lewis DA, McGregor CG, Sieck GC, Fitzpatrick LA, and Miller VM. Gender and transcriptional regulation of NO synthase and ET-1 in porcine aortic endothelial cells. Am J Physiol Heart Circ Physiol 273: H1962–H1967, 1997.[Abstract/Free Full Text]
  36. Williams JK, Adams MR, and Klopfenstein HS. Estrogen modulates responses of atherosclerotic coronary arteries. Circulation 81: 1680–1687, 1990.[Abstract]
  37. Zell R, Geck P, Werdan K, and Boekstegers P. TNF-{alpha} and IL-1{alpha} inhibit both pyruvate dehydrogenase activity and mitochondrial function in cardiomyocytes: evidence for primary impairment of mitochondrial function. Mol Cell Biochem 177: 61–67, 1997.[ISI][Medline]




This Article
Abstract
Full Text (PDF)
All Versions of this Article:
12/3/175    most recent
00040.2002v1
Citation Map
Services
Email this article to a friend
Similar articles in this journal
Similar articles in ISI Web of Science
Similar articles in PubMed
Alert me to new issues of the journal
Download to citation manager
Search for citing articles in:
ISI Web of Science (1)
Google Scholar
Articles by Villablanca, A. C.
Articles by Rutledge, J. C.
Articles citing this Article
PubMed
PubMed Citation
Articles by Villablanca, A. C.
Articles by Rutledge, J. C.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online