Molecular analysis of the vaginal response to estrogens in the ovariectomized rat and postmenopausal woman
© Jelinsky et al; licensee BioMed Central Ltd. 2008
Received: 25 March 2008
Accepted: 25 June 2008
Published: 25 June 2008
Vaginal atrophy (VA) is the thinning of the vaginal epithelial lining, typically the result of lowered estrogen levels during menopause. Some of the consequences of VA include increased susceptibility to bacterial infection, pain during sexual intercourse, and vaginal burning or itching. Although estrogen treatment is highly effective, alternative therapies are also desired for women who are not candidates for post-menopausal hormone therapy (HT). The ovariectomized (OVX) rat is widely accepted as an appropriate animal model for many estrogen-dependent responses in humans; however, since reproductive biology can vary significantly between mammalian systems, this study examined how well the OVX rat recapitulates human biology.
We analyzed 19 vaginal biopsies from human subjects pre and post 3-month 17β-estradiol treated by expression profiling. Data were compared to transcriptional profiling generated from vaginal samples obtained from ovariectomized rats treated with 17β-estradiol for 6 hrs, 3 days or 5 days. The level of differential expression between pre- vs. post- estrogen treatment was calculated for each of the human and OVX rat datasets. Probe sets corresponding to orthologous rat and human genes were mapped to each other using NCBI Homologene.
A positive correlation was observed between the rat and human responses to estrogen. Genes belonging to several biological pathways and GO categories were similarly differentially expressed in rat and human. A large number of the coordinately regulated biological processes are already known to be involved in human VA, such as inflammation, epithelial development, and EGF pathway activation.
At the transcriptional level, there is evidence of significant overlap of the effects of estrogen treatment between the OVX rat and human VA samples.
In many women, declining estrogen levels at menopause lead to a number of physiological changes, including vasomotor instability (hot flushes), mood swings, decreased libido, vaginal atrophy (VA), serum lipid changes and the increased risk of developing osteoporosis. While progestin and/or estrogen therapy are clinically proven to ameliorate some of these symptoms, these medications are contraindicated in some populations, and other women choose not to take them. Because of these issues, recent work has focused on developing more specific, and in many cases, nonhormonal alternatives to traditional post-menopausal hormone therapy. These alternatives have focused largely on reduction of hot flushes [1, 2] and the treatment/prevention of osteoporosis [3–5].
VA affects up to 40% of postmenopausal women, with those with an active sex life being more apt to complain of symptoms (see  for review). Its hallmarks are thinning epithelia and increased vaginal pH. Symptoms include dryness, burning, itching, decreased lubrication during sexual stimulation and painful intercourse. In extreme cases the epithelium is absent and adhesions form in the vaginal vault, preventing clinical exams and drainage of fluid . Currently, estrogens (both systemically and locally administered) are the only approved medical therapies for treating VA. Lubricants can relieve some of the symptoms of VA, but do not treat the underlying cause.
Despite the fact that the vagina is a classic estrogen-responsive organ, it has received little attention compared to the adjacent uterus. Uterine responses to estrogens on the reproductive cycle have been well-characterized and several expression profiles have been generated [8–10], but no analogous studies have been performed on the vagina, aside from a recent report that examines the vaginal effects of neonatal diethylstilbestrol exposure .
Our ultimate goal is to develop nonhormonal therapies for VA in order to provide additional treatment options for post-menopausal women. Key to evaluating new compounds is the development of preclinical animal models, preferably in the rat or mouse. Drawing on the extensive studies of the rat uterine response to estrogens for the development of selective estrogen receptor modulators (SERMs) [12–15], we constructed an analogous model of rat vaginal responses. This report details the time course of mRNA changes elicited by 17β-estradiol in the ovariectomized (OVX) rat vaginal vault and compares these responses to those seen in the human vaginal epithelium after three months of 17β-estradiol patch therapy.
Animals care/tissue collection/dosing
All animal studies were approved by the Institutional Animal Care and Use Committee of Wyeth Research, Collegeville, PA and were conducted in accordance with Association for the Assessment and Accreditation of Laboratory Animal Care (AAALAC) guidelines. Sprague Dawley rats (Taconic, Germantown NY) were OVX at 8 weeks of age and rested for about 1 week before administration of daily subcutaneous doses of 17β-estradiol (20 μg/kg) or vehicle (50% DMSO/50% 1× Dulbecco's PBS). Animals were euthanized after 6 hrs, 3 days and 5 days of treatment. For the 3- and 5-day timepoints, animals were euthanized 6 hours after the last dose. At necropsy, the vaginal vault was excised and trimmed of adherent tissue (e.g. the urethra). Typically, organs were trimmed just below the cervix and above the external orifice of the vagina before flash freezing in liquid nitrogen.
For the studies using betacellulin, 8-week-old female Sprague Dawley rats were OVX and allowed to recover from surgery for 10 days prior to treatment. Rats were dosed intravaginally for 5 days with mouse betacellulin (200 nM; R&D Systems, Minneapolis MN) once daily using a carbomer-934P-based vehicle (pH 7.0) and a 0.5 inch long, 18 gauge ball-tipped gavage needle. Following euthanasia by CO2 asphyxiation and pneumothorax, the vaginal vault was removed, flattened on an index card and fixed in 10% buffered formalin (VWR International, West Chester, PA) for 16–72 hours. Vaginas were equilibrated in 15% sucrose for 3 hours prior to embedding in OCT (Sigma, St Louis MO) and cryosectioning. Ten micron sections were stained with hematoxylin and eosin (both from Sigma) for histological analysis.
Human tissue collection
A clinical study was initiated to collect 3 mm vaginal biopsy specimens from post-menopausal women clinically diagnosed with VA. The diagnosis of VA was confirmed by a vaginal maturation index (VMI) containing ≤5% superficial cells and vaginal pH > 5.0. Patients were treated for 3 months with 17β-estradiol, delivered via a skin patch, at a dose of 0.05 mg/day. All 19 subjects had increased VMI (>5%) and/or reduced pH (< or = 5) following treatment. Vaginal biopsies were collected on day 1 prior to the initiation of treatment with 17β-estradiol and on the last day of treatment (day 85). Samples were immediately frozen on dry ice and stored at -80°C until RNA isolation. This study was conducted at the Ohio State University, Division of Clinical Pharmacology (Columbus, Ohio). Written informed consent was obtained from all subjects prior to any study-related procedure. Western Institutional Review Board (Olympia, WA) approved the protocol and informed consent. The details of the clinical analyses have been reported elsewhere .
Each rat vagina was homogenized in 3 mL of TRIzol reagent (Invitrogen, Carlsbad CA) using a Polytron PT3000 grinder (Brinkmann, Westbury NY) at a setting of 12–15 for 1 min. Subsequent RNA extraction was as described previously . For human samples, RNA extraction was performed according to a modified RNeasy mini kit method. (Qiagen, Valencia, CA) .
ComplementaryDNA was amplified from 50 ng total RNA with the Ovation RNA Amplification System (Nugen, San Carlos CA). Real-time PCR reactions (qRT-PCR) were conducted using the TaqMan® Universal PCR Master Mix (Applied Biosystems). Forward and reverse primers (10 μM) and 20 μM labeled probe was mixed with 2× PCR master mix in a total volume of 15 ul and added to 50 ng cDNA in final reaction volume of 25 μl.
PCR reactions were run on an automated fluorometer (ABI Prism 7000 Sequence Detection System, Applied Biosystems) in a 96-well format. PCR conditions for all reactions were as follows: 1 cycle of 50°C for 2 min., 95°C for 10 min. followed by 40 cycles of 95°C for 15 sec. and 60°C for 1 min. Relative expression was determined by using the CT method (Pontius et al., 2003) using Sequence Detector software, version 1.6.3. Complementary DNAs of each sample were evaluated in triplicate. Results were normalized to a control gene, SFSR10 RNA, which has been shown to have very consistent expression across all patient samples by microarray analysis .
The primer sequences for each transcript are as follows:
BTC: forward primer, 5'-GGTGTGAGAGAGTTGACTTGTTTTACC-3'; reverse primer, 5'-TGCTATCAAACAAATCACCAGAATC-3'; probe, 5'-TGTCCTCTGTCTCCTCT-3'
SFSR10: forward primer, 5'-AAGGGAACACTATACCTGTCATGGA-3'; reverse primer, 5'-GGTAAGCAAAGGACCTGAAAATAATATT-3'; probe, 5'-AAGACTTTGCCTGTTCAT-'
For rat vaginal samples, 5 μg of total RNA were used to generate biotin labeled cRNA using an oligo T7 primer in a reverse transcription reaction followed by in vitro transcription reaction with biotin labeled UTP and CTP. Ten μg of cRNA were fragmented and hybridized to RAE230A and RAE230B arrays (Affymetrix, Santa Clara, CA). For human biopsies, RNA was labeled using Ovation™ Biotin System (Nugen, San Francisco, CA). Labeled cDNA (2.5 μg) was fragmented and hybridized to HGU133 Plus 2.0 arrays (Affymetrix, Santa Clara, CA) . Hybridized arrays were stained according to manufacturer's protocols on a Fluidics Station 450 and scanned on an Affymetrix scanner 3000. All array images were visually inspected for defects and quality. Arrays with excessive background, low signal intensity, or major defects within the array were eliminated from further analysis. Signal values were determined using Gene Chip Operating System 1.0 (GCOS, Affymetrix). For each array, all probe sets were normalized to a mean signal intensity value of 100. The default GCOS statistical values were used for all analysis. Signal values and absolute detection calls were imported into Genesis 2.0 (GeneLogic, Gaithersburg, MD) for analysis. The complete data set is available at the Gene Expression Omnibus as GSE11622 .
Determination of probe sets corresponding to homologous gene pairs
Mapping of RAE230 to HGU133 plus 2.0 probe sets was based on NCBI's Homologene (clusters of orthologous genes). First, all qualifiers were annotated with unigene accession numbers through BioExpress (GeneLogic). These accession numbers then were used to map orthologous sequences through NCBI Homologene . Commonly one probe set from one species maps to multiple probe sets from a second species. To resolve multiple mappings, clusters of orthologous probe sets were generated as follows. For each human probe set, all corresponding rat probe sets were identified and grouped into a "cluster". Human probe sets that map to the same set of rat probes were clustered. For each human cluster, the corresponding locuslink ID(s) were identified. Clusters with multiple Locuslink ID annotations were discarded. Approximately 10% of the human probe sets clustered to multiple Locuslink IDs. Clustering identified 9093 groups, spanning 22460 human probe set IDs and 13219 rat probe set IDs. Each cluster was annotated (assigned GO/KEGG etc categories) based on the human Locuslink ID in that cluster.
Identification of regulated genes
Identification of estrogen-regulated genes was done in two independent ways. In the first method, probe sets were considered estrogen regulated if the following conditions were met: 1) A gene was expressed greater than 50 signal units, 2) the percentage of samples with a Present (P) call as determined by GCOS default settings was greater than or equal to 66% in any sample condition, 3) the fold change between all estrogen treated and vehicle treated samples was at least 1.7, and 4) the p-value based on an Wilcoxon t-test was ≤0.01 (rat) or based on a Wilcoxon paired test was <0.01 (human). Over 5300 rat probe sets and 1668 human probe sets met these conditions. The second method used the 9093 clusters identified above. Arrays were subjected to a Loess transformation and the CyberT-statistics were calculated. False discovery rate estimates (q-values) were estimated using a null distribution of CyberT's from permuted data.
Identification of significantly regulated gene sets
Statistical identification of significantly regulated biological pathways was implemented using a modified approach of Tian et al. . 7798 gene sets were tested in the sigpathway analysis. These were derived from the follow sources: 1) GO annotations (from Affymetrix) , 2) KEGG/GenMAPP annotations (from Affymetrix) , 3) Functional categories (from MSigDB), 4) Genes that are located in close proximity along the chromosome (from MSigDB), 5) Genes containing instances of a regulatory motif of interest (from MSigDB), and 6) Genes that are co-regulated with a given gene in a variety of datasets (from MSigDB). Gene sets from MSigDB were derived from pathway databases (Biocarta, KEGG, and BioCyc) and pathway-specific microarray annotations , as well as >5,000 gene sets from Gene Ontology. Only gene sets in which the number of genes was greater or equal to 10 and less than or equal to 400 were used for further analysis. Each gene set was assigned a score of differential expression between two conditions (essentially, the sum of t-statistics for each gene within the set), and the false discovery rate (Q) associated with this score was calculated via two separate permutation tests (permuting gene labels (Q1) and sample labels (Q2)). A gene set was considered significant if its two false discovery rates are below 0.15 (Q1 <= 0.15 and Q2 <= 0.15) for all 4 comparisons. 45 gene sets passed these criteria.
Estrogen-responsive genes in the rat vaginal vault
Number of probe sets regulated by 17-β estradiol in rat vaginal tissue
Genes regulated at 5 days
Genes regulated at 3 days
Genes regulated at 6 hours
Differentially expressed genes in the human
Human vaginal biopsy specimens from 19 postmenopausal women clinically diagnosed with VA were collected. Participants in the study were treated for 3 months with 17β-estradiol, delivered via a skin patch, with vaginal biopsies collected prior to the initiation of treatment and on the last day of treatment. Microarray analysis of these biopsies resulted in 1668 genes regulated by at least 1.7 fold (772 up regulated and 896 down regulated, Wilcoxon paired t-test p < 0.01). The details of this analysis are described elsewhere . Like the rat vaginal samples, we find the human vagina to be very responsive to estradiol although we see fewer genes regulated in the human samples. The difference in the number of significant genes could be the result of several factors, such as the greater variability (in genetics, age, etc) amongst the human patients in comparison to the inbred rats, or the fact that the human biopsies contained primarily epithelial cells, while the rat samples were the entire vaginal vaults.
Mapping of RAE230 to HGU133 Plus 2.0 probe sets
Top coregulated clusters of genes
growth arrest and DNA-damage-inducible, alpha
3-hydroxy-3-methylglutaryl-Coenzyme A synthase 1
family with sequence similarity 46, member B
cellular retinoic acid binding 2
carbonyl reductase 3
secreted LY6/PLAUR domain
ATPase, H+ transporting, lysosomal, V1 subunit D
Kruppel-like factor 4
glutamic pyruvate transaminase
chloride intracellular channel 3
T-cell differentiation protein 2
S100 calcium binding protein A10
hypothetical protein FLJ12151
3-hydroxy-3-methylglutaryl-Coenzyme A reductase
musculoaponeurotic fibrosarcoma homolog B
RAB10, member RAS oncogene family
serine proteinase inhibitor, clade B 2
arginine vasopressin-induced 1
absent in melanoma 1-like
BCL2/adenovirus E1B 19 kDa interacting protein 3
basic leucine zipper and W2 domains 1
hypothetical protein MGC13102
eukaryotic translation termination factor 1
GPI-anchored metastasis-associated protein homolog
heme oxygenase 1
K+ channel tetramerisation domain containing 11
cyclin-dependent kinase 7
jun B proto-oncogene
RIO kinase 3
malic enzyme 1
tyrosine 3-monooxygenase activation protein, zeta
Top coregulated clusters of genes
thioredoxin interacting protein
retinol binding protein 1, cellular
insulin-like growth factor binding protein 3
interferon induced transmembrane protein 3 (1–8 U)
dehydrogenase/reductase (SDR family) member 8
lactate dehydrogenase B
lysosomal associated protein transmembrane 4 beta
mucin 1, transmembrane
secreted frizzled-related protein 2
interferon induced transmembrane protein 1 (9–27)
protein phosphatase 1, regulatory (inhibitor) subunit 1B
protein kinase, cAMP-dependent, catalytic, beta
transmembrane channel-like 4
integral membrane protein 2C
dehydrogenase/reductase (SDR family) member 3
integral membrane protein 2B
Niemann-Pick disease, type C2
AXL receptor tyrosine kinase
major histocompatibility complex, class II, DR alpha
Eukaryotic translation elongation factor 1 alpha 1
forkhead box Q1
interferon gamma receptor 1
insulin-like growth factor binding protein 7
basic helix-loop-helix domain containing, class B, 3
selenium binding protein 1
chromosome 11 open reading frame 8
lymphocyte-specific protein 1
Correlation between Rat and Human Data
Enriched pathways of estrogen regulated genes in human and rat
Cholesterol_Biosynthesis // Genmapp
Cellular Lipid Metabolism // GO:0044255
Mitotic Cell Cycle // GO:0000278
Cysteine Protease Inhibitor Activity // GO:0004869
Small Protein Conjugating Enzyme Activity // GO:0008639
Epidermis Development // GO:0008544
Steroid Metabolism // GO:0008202
Ubiquitin Conjugating Enzyme Activity // GO:0004840
Chr3Q13: Cytogenetic Band
Complement_Activation_Classical // GenMAPP
TGF_Beta_Signaling_Pathway // GenMAPP
Matrix_Metalloproteinases // GenMAPP
Complement Activation, Classical Pathway // GO:0006958
Cell Differentiation // GO:0030154
chr1q24: Cytogenetic band
Locomotory Behavior // GO:0007626
Transmembrane Receptor Protein Tyrosine Kinase Activity // GO:0004714
Specific Rna Polymerase Ii Transcription Factor Activity // GO:0003704
RAR_UP: Manually Curated
Response To Wounding // GO:0009611
Cell Motility // GO:0006928
Locomotion // GO:0040011
classic complement pathway
Vascular Endothelial Growth Factor Receptor Activity // GO:0005021
Transmembrane Receptor Protein Kinase Activity // GO:0019199
Immune Complement Pathway
STTTCRNTTT: IRF:Interferon regulatory factor
Response To Virus // GO:0009615
Cell Adhesion // GO:0007155
Humoral Defense Mechanism // GO:0016064
Humoral Immune Response // GO:0006959
Antigen Processing // GO:0030333
Antigen Presentation // GO:0019882
Response To External Biotic Stimulus // GO:0043207
Response To Pest, Pathogen Or Parasite // GO:0009613
Immune Response // GO:0006955
Experimental validation of betacellulin
BTC expression levels increase significantly after estrogen treatment in both rat (Figure 5B) and human. Because expression levels of BTC on the human array were relatively low, qRT-PCR was used to confirm regulation in the human (Figure 5C). BTC is upregulated in 14 of the 17 patients examined. Therefore, it appears that BTC has similar regulation in rat and human. Administration of BTC protein intravaginally in OVX rats results in a thickening and proliferation of the vaginal epithelium as detected by histology (Figure 5A, panel BTC). BTC administration was also noted in some cases to induce abnormal vaginal epithelial proliferation (Figure 5, panel BTC-abn) consistent with its mitogenic role.
The current animal model of VA involves the use of OVX rats. Since the ovaries are the major source of estrogens, removal of the ovaries mimics an estrogen-deprived state. It is well documented that as women enter menopause, the levels of circulating estrogens decrease, resulting in a number of physiological changes, including VA. VA can be reversed through estrogen hormone replacement, suggesting that loss of estrogen is the major factor leading to VA. OVX animal systems mimic many human postmenopausal systems, in both phenotype and histology, and have been well-characterized models for studying these changes. We have taken a systematic molecular approach to determine if the animal model for VA is similar to the human disorder at the transcriptional level.
We have used this experiment to identify genes and pathways that are estrogen regulated in both rat and human vaginal tissue. It is important to note that there are differences in the experimental conditions. First, we evaluated the animal model after 6 hours, 3 or 5 days of treatment. At the 5-day time point, the vaginal epithelium has fully responded to estrogen treatment by histology. We have not, however, compared the transcriptional profiles of normal vaginal tissue across the estrus cycle to vaginal tissue from an OVX animal treated with an estrogen. It is possible that there are differences between the estrogen treated, OVX-animals vs. ovary-intact animals. Second, there have been no comparable studies done in postmenopausal women to determine the time course of histological changes following estrogen therapy. As a result, we chose to look at a time point (3 months) in which we would be sure to see a significant effect. Due to differences in collection times, it is quite possible that we are looking at slightly different molecular endpoints between the human and the OVX rat. Third, it is also worth noting that the tissues being profiled are different. In the animal model we examined the entire vaginal vault, while the human samples contained primarily the vaginal epithelia. Estrogen response likely depends on the cell type composition.
Furthermore, because the rat transcriptome is less well annotated than human, inaccurate or incomplete mappings between the two chipsets may result. Our methods involved grouping together human probe sets mapping to the same gene, calculating the average of their t-statistics, and comparing this value to the average t-statistic for the orthologous rat gene. This procedure was chosen after observing that this metric yielded greater agreement between human and rat in a kidney-vs-liver comparison (see below). However, averaging the results for probe sets querying the same gene could obscure splice-variant-specific responses to estrogen, or dilute desirable signal with poorly performing probe sets. Given the above caveats, it is still encouraging that we find a significant positive correlation in regulation between the two species, suggesting similar underlying molecular mechanisms for estrogen response. This study is also consistent with other validations of this animal model in the literature. OVX animals have long been used to model the effects of estrogens on the reproductive tract and skeleton, and have been validated through a number of different criteria, including histology and response to clinically active agents.
While there is a significant positive correlation of the 17β-estradiol response in rat and human, it should be noted that the correlation is not perfect. There are many possible technical reasons for the differences as described above. To put this in perspective, we performed an identical analysis, where instead of comparing pre- vs post- estrogen treatment, we compared kidney vs liver samples in each species. A comparison between two very distinct tissues represented an idealized dataset, where the results should be the same for human and rat, and also a large degree of differential expression is expected, due to tissue-specific genes. In this regard, we considered the degree of agreement between human and rat for this tissue comparison to be the upper limit of what can be expected. The t-statistic correlation between species was 0.48 (analogous to Figure 2). While 0.48 is higher than the 0.33 we see in our experiments, this is not very surprising since a treatment response within a single tissue is a subtle effect, relative to a comparison between two different tissues. Furthermore, McCarroll et al.  studied the transcriptional response to aging in flies vs. worms. They found a t-statistic correlation 0.15 which they note to be highly significant. We expect to have a higher correlation coefficient than 15% since rat and human are evolutionarily closer than flies and worms. Therefore, we feel the ~30% correlation coefficient is in the expected value range (between 15% and 48%).
We were initially surprised that the majority of genes regulated by 17β-estradiol were suppressed rather than stimulated. Our assumption was that 17β-estradiol would exert its pro-proliferatory activities primarily by upregulating genes. However, our data are in line with those obtained from microarray studies done in MCF-7 breast cancer cells, where estrogens also have a pro-proliferatory effect: 70% of the regulated genes were suppressed .
Estrogen is the gold standard therapy for treatment of VA and the identification of estrogen regulated genes in both human and rat likely yields genes responsible for the estrogenic effect. We have chosen to validate a gene that is estrogen-responsive in both rat and human vaginal samples. Betacellulin is a member of the epidermal growth factor family of peptide ligands and is a potent mitogen for a wide variety of cell types affecting cell proliferation, differentiation, survival and migration. Estrogens are known to stimulate expression of other members of this family, such as amphiregulin . It has functional redundancy to other members of the EGF family . The increased epithelial thickening seen in the OVX vagina is consistent with its mitogenic role and suggests that BTC may have a normal role in estrogen-dependent epithelial maintenance. Additionally, BTC is able to stimulate mammary gland differentiation as well as promote the survival of breast cancer cell lines [31, 32]. However, the uncontrolled proliferation seen in the OVX treated vagina suggests that 17β-estradiol may be regulating other factors to control BTC's mitogenic effect.
Here we have identified estrogen-regulated genes in the vagina of humans and a commonly used rat animal model. In addition to demonstrating that the vagina is a highly estrogen-responsive organ at the molecular level, we show that the transcriptional response to estrogen is highly correlated between the two species. Furthermore, we hypothesized that genes with estrogen-sensitivity in both species are more likely to be responsible for the changes observed in estrogen treated OVX rats. We have shown that BTC, a gene regulated by estrogen in both human and rat, does in fact have the ability to recapitulate the estrogen induced epithelial proliferation. It is likely that many other genes also possess the ability to induce epithelial proliferation and current investigations are underway to identify them.
Vaginal maturation index
Kyoto Encyclopedia of Genes and Genomes
mitochondrial interest group database
epidermal growth factor
estrogen receptor alpha
Affymetrix gene chip operating system 1.0.
The authors thank the Wyeth Collegeville BioResources department for veterinary and technical support.
- Adelson KB, Loprinzi CL, Hershman DL: Treatment of hot flushes in breast and prostate cancer. Exper Opin Pharmaco. 2005, 6 (7): 1095-1106. 10.1517/14656522.214.171.1245.View ArticleGoogle Scholar
- Deecher DC: Physiology of thermoregulatory dysfunction and current approaches to the treatment of vasomotor symptoms. Expert Opin Investig Drugs. 2005, 14 (4): 435-448. 10.1517/135437126.96.36.1995.View ArticlePubMedGoogle Scholar
- Wada T, Nakashima T, Hiroshi N, Penninger JM: RANKL-RANK signaling in osteoclastogenesis and bone disease. Trends Mol Med. 2006, 12 (1): 17-25. 10.1016/j.molmed.2005.11.007.View ArticlePubMedGoogle Scholar
- Rubin MR, Bilezikian JP: Parathyroid hormone as an anabolic skeletal therapy. Drugs. 2005, 65 (17): 2481-2498. 10.2165/00003495-200565170-00005.View ArticlePubMedGoogle Scholar
- Grey A, Reid IR: Emerging and potential therapies for osteoporosis. Expert Opin Investig Drugs. 2005, 14 (3): 265-278. 10.1517/135437188.8.131.525.View ArticlePubMedGoogle Scholar
- Willhite LA, O'Connell MB: Urogenital atrophy: Prevention and treatment. Pharmacotherapy. 2001, 21 (4): 464-480. 10.1592/phco.21.5.464.34486.View ArticlePubMedGoogle Scholar
- Woodruff JD, Parmley TH: The Human Vagina. 1978, New York , Elsevier North-Holland Biomedical Press, 2 : 291-308. Nonobstetrical vaginal injuries, Hafez ESE, Evans TN, Human Reproductive Medicine, Hafez ESE.Google Scholar
- Bauersachs S, Ulbrich SE, Gross K, Schmidt SEM, Meyer HHD, Einspanier R, Wenigerkind H, Vermehren M, Blum H, Sinowatz F, Wolf E: Gene expression profiling of bovine endometrium during the oestrous cycle: detection of molecular pathways involved in functional changes. J Mol Endocrinol. 2005, 34 (3): 889-908. 10.1677/jme.1.01799.View ArticlePubMedGoogle Scholar
- Borthwick JM, Charnock-Jones DS, Tom BD, Hull ML, Teirney R, Phillips SC, Smith SK: Determination of the transcript profile of human endometrium. Mol Hum Reprod. 2003, 9 (1): 19-33. 10.1093/molehr/gag004.View ArticlePubMedGoogle Scholar
- Wu XX, Pang ST, Sahlin L, Blanck A, Norstedt G, Flores-Morales A: Gene expression profiling of the effects of castration and estrogen treatment in the rat uterus. Biol Reprod. 2003, 69 (4): 1308-1317. 10.1095/biolreprod.103.015420.View ArticlePubMedGoogle Scholar
- Suzuki A, Watanabe H, Mizutani T, Sato T, Ohta Y, Iguchi T: Global gene expression in mouse vaginae exposed to diethylstilbestrol at different ages. Exp Biol Med (Maywood). 2006, 231 (5): 632-640.Google Scholar
- Geiser AG, Hummel CW, Draper MW, Henck JW, Cohen IR, Rudmann DG, Donnelly KB, Adrian MD, Shepherd TA, Wallace OB, McCann DJ, Oldham SW, Bryant HU, Sato M, Dodge JA: A new selective estrogen receptor modulator with potent uterine antagonist activity, agonist activity in bone, and minimal ovarian stimulation. Endocrinology. 2005, 146 (10): 4524-4535. 10.1210/en.2005-0024.View ArticlePubMedGoogle Scholar
- Ke HZ, Brown TA, Thompson DD: Lasofoxifene (CP-336,156), a novel selective estrogen receptor modulator, in preclinical studies. J Am Aging Assoc. 2002, 25 (2): 87-99.PubMedPubMed CentralGoogle Scholar
- Komm BS, Kharode YP, Bodine PVN, Harris HA, Miller CP, Lyttle CR: Bazedoxifene acetate: A selective estrogen receptor modulator with improved selectivity. Endocrinology. 2005, 146 (9): 3999-4008. 10.1210/en.2005-0030.View ArticlePubMedGoogle Scholar
- Sorbera LA, Castaner J, Silvestre JS: Pipendoxifene - Treatment of breast cancer - Estrogen receptor modulator. Drugs of the Future. 2002, 27 (10): 942-947. 10.1358/dof.2002.027.10.703467.View ArticleGoogle Scholar
- Cotreau M, Chennathukuzhi V, Harris HA, Han L, Dorner A, Hatstat E, Zakaria M, Apseloff G, Strahs A, Crabtree JS, Winneker R, Jelinsky SA: A Study of 17b-estradiol-Regulated Genes in the Vagina of Postmenopausal Women with Vaginal Atrophy . Maturitas. 2007, 58 (4): 366-376. 10.1016/j.maturitas.2007.09.009.View ArticlePubMedGoogle Scholar
- Crabtree JS, Zhang X, Peano BJ, Zhang Z, Winneker RC, Harris HA: Development of a mouse model of mammary gland versus uterus tissue selectivity using estrogen- and progesterone-regulated gene markers. J Steroid Biochem Mol Biol. 2006, 101 (1): 11-21. 10.1016/j.jsbmb.2006.06.017.View ArticlePubMedGoogle Scholar
- Barrett T, Troup DB, Wilhite SE, Ledoux P, Rudnev D, Evangelista C, Kim IF, Soboleva A, Tomashevsky M, Edgar R: NCBI GEO: mining tens of millions of expression profiles--database and tools update. Nucleic Acids Res. 2007, 35 (Database issue): D760-5. 10.1093/nar/gkl887.View ArticlePubMedGoogle Scholar
- Wheeler DL, Barrett T, Benson DA, Bryant SH, Canese K, Church DM, DiCuccio M, Edgar R, Federhen S, Helmberg W, Kenton DL, Khovayko O, Lipman DJ, Madden TL, Maglott DR, Ostell J, Pontius JU, Pruitt KD, Schuler GD, Schriml LM, Sequeira E, Sherry ST, Sirotkin K, Starchenko G, Suzek TO, Tatusov R, Tatusova TA, Wagner L, Yaschenko E: Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2005, 33 (Database issue): D39-45. 10.1093/nar/gki062.View ArticlePubMedGoogle Scholar
- Tian L, Greenberg SA, Kong SW, Altschuler J, Kohane IS, Park PJ: Discovering statistically significant pathways in expression profiling studies. Proc Natl Acad Sci U S A. 2005, 102 (38): 13544-13549. 10.1073/pnas.0506577102.View ArticlePubMedPubMed CentralGoogle Scholar
- Affymetrix Net Affx. [http://www.affymetrix.com]
- Superarray. [http://www.superarray.com]
- Darwish H, Krisinger J, Furlow JD, Smith C, Murdoch FE, DeLuca HF: An estrogen-responsive element mediates the transcriptional regulation of calbindin D-9K gene in rat uterus. J Biol Chem. 1991, 266 (1): 551-558.PubMedGoogle Scholar
- Klinge CM: Estrogen receptor interaction with estrogen response elements. Nucleic Acids Res. 2001, 29 (14): 2905-2919. 10.1093/nar/29.14.2905.View ArticlePubMedPubMed CentralGoogle Scholar
- O'Lone R, Frith MC, Karlsson EK, Hansen U: Genomic targets of nuclear estrogen receptors. Mol Endocrinol. 2004, 18 (8): 1859-1875. 10.1210/me.2003-0044.View ArticlePubMedGoogle Scholar
- Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005, 102 (43): 15545-15550. 10.1073/pnas.0506580102.View ArticlePubMedPubMed CentralGoogle Scholar
- Dunbar AJ, Goddard C: Structure-function and biological role of betacellulin. Int J Biochem Cell Biol. 2000, 32 (8): 805-815. 10.1016/S1357-2725(00)00028-5.View ArticlePubMedGoogle Scholar
- Levin ER: Bidirectional signaling between the estrogen receptor and the epidermal growth factor receptor. Mol Endocrinol. 2003, 17 (3): 309-317. 10.1210/me.2002-0368.View ArticlePubMedGoogle Scholar
- McCarroll SA, Murphy CT, Zou S, Pletcher SD, Chin CS, Jan YN, Kenyon C, Bargmann CI, Li H: Comparing genomic expression patterns across species identifies shared transcriptional profile in aging. Nat Genet. 2004, 36 (2): 197-204. 10.1038/ng1291.View ArticlePubMedGoogle Scholar
- Frasor J, Danes JM, Komm B, Chang KC, Lyttle CR, Katzenellenbogen BS: Profiling of estrogen up- and down-regulated gene expression in human breast cancer cells: insights into gene networks and pathways underlying estrogenic control of proliferation and cell phenotype. Endocrinology. 2003, 144 (10): 4562-4574. 10.1210/en.2003-0567.View ArticlePubMedGoogle Scholar
- Beerli RR, Hynes NE: Epidermal growth factor-related peptides activate distinct subsets of ErbB receptors and differ in their biological activities. J Biol Chem. 1996, 271 (11): 6071-6076. 10.1074/jbc.271.11.6071.View ArticlePubMedGoogle Scholar
- Daly JM, Olayioye MA, Wong AM, Neve R, Lane HA, Maurer FG, Hynes NE: NDF/heregulin-induced cell cycle changes and apoptosis in breast tumour cells: role of PI3 kinase and p38 MAP kinase pathways. Oncogene. 1999, 18 (23): 3440-3451. 10.1038/sj.onc.1202700.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/1/27/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.