Impact of the Phytoestrogen Content of Laboratory Animal Feed on the Gene Expression Profile of the Reproductive System in the Immature Female Rat
2004; National Institute of Environmental Health Sciences; Volume: 112; Issue: 15 Linguagem: Inglês
10.1289/ehp.6848
ISSN1552-9924
AutoresJorge M. Naciff, Gary J. Overmann, Suzanne M. Torontali, Gregory J. Carr, Jay P. Tiesman, George P. Daston,
Tópico(s)Nutrition, Genetics, and Disease
ResumoVol. 112, No. 15 ResearchOpen AccessImpact of the Phytoestrogen Content of Laboratory Animal Feed on the Gene Expression Profile of the Reproductive System in the Immature Female Rat Jorge M. Naciff, Gary J. Overmann, Suzanne M. Torontali, Gregory J. Carr, Jay P. Tiesman, and George P. Daston Jorge M. Naciff , Gary J. Overmann , Suzanne M. Torontali , Gregory J. Carr , Jay P. Tiesman , and George P. Daston Published:1 November 2004https://doi.org/10.1289/ehp.6848Cited by:24AboutSectionsPDF ToolsDownload CitationsTrack Citations ShareShare onFacebookTwitterLinked InReddit AbstractThe effect of the dietary background of phytoestrogens on the outcome of rodent bioassays used to identify and assess the reproductive hazard of endocrine-disrupting chemicals is controversial. Phytoestrogens, including genistein, daidzein, and coumestrol, are fairly abundant in soybeans and alfalfa, common ingredients of laboratory animal diets. These compounds are weak agonists for the estrogen receptor (ER) and, when administered at sufficient doses, elicit an estrogenic response in vivo. In this study, we assessed the potential estrogenic effects of dietary phytoestrogens at the gene expression level, together with traditional biologic end points, using estrogen-responsive tissues of the immature female rat. We compared the gene expression profile of the uterus and ovaries, as a pool, obtained using a uterotrophic assay protocol, from intact prepubertal rats fed a casein-based diet (free from soy and alfalfa) or a regular rodent diet (Purina 5001) containing soy and alfalfa. Estrogenic potency of the phytoestrogen-containing diet was determined by analyzing uterine wet weight gain, luminal epithelial cell height, and gene expression profile in the uterus and ovaries. These were compared with the same parameters evaluated in animals exposed to a low dose of a potent ER agonist [0.1 μg/kg/day 17α-ethynyl estradiol (EE) for 4 days]. Exposure to dietary phytoestrogens or to a low dose of EE did not advance vaginal opening, increase uterine wet weight, or increase luminal epithelial cell height in animals fed either diet. Although there are genes whose expression differs in animals fed the soy/alfalfa-based diet versus the casein diet, those genes are not associated with estrogenic stimulation. The expression of genes well known to be estrogen regulated, such as progesterone receptor, intestinal calcium-binding protein, and complement component 3, is not affected by consumption of the soy/alfalfa-based diet when assessed by microarray or quantitative reverse transcriptase–polymerase chain reaction analysis. Our results indicate that although diet composition has an impact on gene expression in uterus and ovaries, it does not contribute to the effects of an ER agonist.The effect of the dietary background of phytoestrogens on the outcome of rodent bioassays used to identify and assess the reproductive hazard of endocrine-disrupting chemicals is controversial. Phytoestrogens, including genistein, daidzein, and coumestrol, are fairly abundant in soybeans and alfalfa, common ingredients of laboratory animal diets. In fact, soy and alfalfa are commonly used as protein sources in the manufacture of most rodent diets. Some of these ingredients are known to contain endocrine modulators, such as the phytoestrogens genistein and daidzein (abundant in soybeans and its products) and their respective glycosides (genistin and daidzin), and coumestrol (found in alfalfa). These phytoestrogens are able to bind to both estrogen receptor (ER) isoforms, ER-αand ER-β, in vitro (Beck et al. 2003; Casanova et al. 1999). They have a higher affinity for ER-β (Boue et al. 2003), but they activate both ER isoforms, although with less potency than estradiol. Both genistein and daidzein have much weaker affinities than does 17β-estradiol for the rat ERs: genistein binds 3- and 100-fold weaker, and daidzein binds 60- and 1,000-fold weaker to rat ER-β and ER-α, respectively (Boue et al. 2003; Casanova et al. 1999). These two phytoestrogens are able to elicit estrogenic responses in vivo (Boettger-Tong et al. 1998; Brown and Setchell 2001; Degen et al. 2002; Jefferson et al. 2002; Levy et al. 1995; Odum et al. 2001; Thigpen et al. 1997). The selective interaction of phytoestrogens with human ER-αand ER-β is similar in vitro to that described for the rat (Casanova et al. 1999; Kuiper et al. 1998).Genistein is also known to have other activities, such as inhibition of different enzymes, among them tyrosine kinases (Akiyama et al. 1987), nitric oxide synthase (Duarte et al. 1997), and topoisomerase II (Okura et al. 1988), and decreasing calcium-channel activity in neurons (Potier and Rovira 1999). It also decreases lipid peroxidation (Arora et al. 1998) and diacylglycerol synthesis (Dean et al. 1989). Therefore, the multiple biologic activities of phytoestrogens raise the question of whether they have the potential to influence the outcome and/or interpretation of bioassays used to identify chemicals with estrogenic potential. In particular, questions have been raised about the presence of phytoestrogens in diets fed to animals used in bioassays designed to screen chemicals that may act as weak regulators of ERs and to screen low doses of potent regulators of ERs (Thigpen et al. 1997, 2002). One such bioassay is the uterotrophic assay, designed to evaluate both ER agonists and antagonists.By using a version of the uterotrophic assay in the immature rat, one of the tier I screening assays recommended for detecting the estrogenic properties of endocrine-disrupting chemicals [Organisation for Economic Co-operation and Development (OECD) 2001; U.S. Environmental Protection Agency (U.S. EPA) 1998], we have identified a set of genes from the uterus and ovaries of prepubertal rats for which expression is regulated by estrogen exposure in a dose-dependent manner and which have the potential to be used as biomarkers for estrogen activity (Naciff et al. 2003). Gene expression changes induced by estrogen stimulation are more sensitive than the classical end points (i.e., uterine weight increase) for evaluating estrogenicity (Naciff et al. 2003). Given that components of the rodent diet commonly used in reproductive toxicology studies include chemicals with known estrogenic activity, understanding the influence of diet and dietary components on estrogen response is an important issue. In this study, we used gene expression profiling to evaluate the effect of two diets with different phytoestrogen content on the transcript profile of two organs that are responsive to estrogen stimulation: the uterus and the ovaries of prepubertal rats.Materials and MethodsChemicals.17α-Ethynyl estradiol (EE) and peanut oil were obtained from Sigma Chemical Company (St. Louis, MO).Animals and treatments.Fifteen-day-old female Sprague-Dawley rats were obtained (Charles River VAF/Plus; Charles River Laboratories, Raleigh, NC) in groups of 10 pups per surrogate mother. We chose this rat strain because it is commonly used in reproductive and developmental toxicity studies. The rats were acclimated to the local vivarium conditions (24°C; 12-hr light/12-hr dark cycle) for 5 days and were fed a casein-based diet (soy- and alfalfa-free diet; Purina 5K96, Purina Mills, St. Louis, MO). Starting on post-natal day (PND)20 and during the experimental phase of the protocol, all rats were singly housed in 20 × 32 × 20 cm plastic cages. To test the diet effect, there were two animal groups (n = 20): one group was fed a standard laboratory rodent diet (Purina 5001, Purina Mills), and the other group was maintained on the casein-based diet. The Purina 5001 diet contains phytoestrogens, mostly genistein and daidzein derived from soy and alfalfa, at levels that may have an impact on the gene expression profile (total daidzein + genistein = 0.49 mg/g; Thigpen et al. 1999), particularly in tissues regulated by estrogens such as reproductive tissues. However, those levels are not uterotrophic when evaluated by the traditional end points, uterine weight gain and increase in luminal epithelial cell height. The casein-based diet is essentially phytoestrogen free, consistently containing < 1 ppm aglycone equivalents of genistein, daidzein, and glycitein, and was fed to the four groups of animals from PND16 onward in order to remove any possible effects of the regular rodent diet (Purina 5001) previously fed to the rats by the animal supplier. All the animals were allowed free access to water and specific pelleted commercial diet (Purina 5001 or casein-based 5K96). The experimental protocol was carried out according to Procter and Gamble's animal care approved protocols, and animals were maintained in accordance with the NIH Guide for the Care and Use of Laboratory Animals (Institute of Laboratory Animal Resources 1996).Starting on PND20, each diet group was divided into two subgroups of 10 animals. One subgroup from each diet subgroup was dosed by subcutaneous injection with 0.1 μg/kg/day EE in peanut oil. This dose is not sufficient to induce a uterotrophic response in juvenile rats (Kanno et al. 2002; Naciff et al. 2003). Animals received 5 mL/kg body weight of dose solution once a day for 4 days. A 4-day dosing regime was selected to optimize detection of any effect of EE exposure at this low dose, both at the histologic level and at the gene expression level. The dose was administered between 0800 and 0900 hr each day. Controls, fed with the appropriate diet, received 5 mL/kg of peanut oil once a day for 4 days. Doses were administered on a microgram per kilogram body weight basis and adjusted daily for weight changes. Body weight (nearest 1.0 g) and the volume of the dose administered (nearest 0.1 mL) were recorded daily. The exact time of the last dose was recorded, to establish a 24-hr waiting period before tissue collection. The animals were sacrificed by CO2 asphyxiation 24 hr after the last dosing, on PND24. The body of the uterus, cut just above its junction with the cervix, with the ovaries attached, was carefully dissected free of adhering fat and mesentery and was weighed as a whole. Then, the ovaries were dissected free, and the uterine and ovarian wet weight was recorded. Both the uterus and ovaries were placed into RNAlater (50–100 mg/mL of solution; Ambion, Austin, TX) at room temperature.Histology.Reproductive tissues from two animals in each dose group were fixed in 10% neutral buffered formalin immediately after weighing and then dehydrated and embedded in paraffin. Serial 4–5 μm cross sections were made through the ovaries, oviducts, and uterine horns, which were stained with hematoxylin and eosin. The evaluation of the morphologic changes induced by the two different diets with or without EE exposure in the uterus was performed as described previously (Naciff et al. 2003).Expression profiling.We used 10 μg total RNA, extracted from uterus and ovaries from individual animals (combining only the tissues from the same animal), to prepare biotin-labeled cRNA, as previously described (Naciff et al. 2002, 2003). Labeled cRNA samples were hybridized to the Affymetrix GeneChip Test 3 Array (Affymetrix Inc., Santa Clara, CA) to assess the overall quality of each sample. After determining the target cRNA quality, we selected individual samples of pooled uteri/ovaries from five or six individual females (replicates) from each diet group, from controls, and from EE-treated subgroups (with high quality cRNA) and hybridized them to Affymetrix Rat Genome U34A high-density oligonucleotide microarrays for 16 hr. The microarrays were washed and stained by streptavidin-phycoerythrin to detect bound cRNA. The signal intensity was amplified by second staining with biotin-labeled anti-streptavidin antibody and followed by streptavidin-phycoerythrin staining. Fluorescent images were read using the Hewlett-Packard G2500A gene array scanner (Affymetrix Inc.). Affymetrix image files for the 20 chip hybridizations, and the absolute analysis results of each diet group are available from the authors upon request.Real-time reverse transcriptase-polymerase chain reaction.In order to corroborate the changes in gene expression identified by the oligonucleotide microarrays, we used a real-time (kinetic) quantitative reverse transcriptase-polymerase chain reaction (QRT-PCR) approach, as previously described (Naciff et al. 2002). This approach allowed us to evaluate the "basal level" of expression of individual genes in samples derived from animals exposed to the two different diets used in our study, as well as changes induced by low-dose EE exposure (0.1 μg/kg/day). We compared the transcript level of selected genes in samples derived from animals in all experimental groups. To confirm the amplification specificity from each primer pair, the amplified PCR products were size-fractioned by electrophoresis in a 4% agarose gel in Tris borate ethylene diamine tetraacetic acid buffer and photographed after staining with ethidium bromide. Table 1 shows the nucleotide sequences for the primers used to test the indicated gene products. Preliminary experiments were done with each primer pair to determine the overall quality and specificity of the primer design. After QRT-PCR, we observed only the expected products at the correct molecular weight.Table 1 Primers used to verify the array-based gene expression changes induced by the two different diets, by QRT-PCR.Gene nameGenBank accession no.aForward primerReverse primerAmplicon size (bp)Complement component 3 (CC3)M298665′-CGTGAGCAGCACAGAAGAGA-3′5′-CCAGGTGGTGATGGAATCTT-3′204Progesterone receptor (PgR)L169225′-CATGTCAGTGGACAGATGCT-3′5′-ACTTCAGACATCATTTCCGG-3′428Intestinal calcium-binding protein (icabp)K009945′-ATCCAAACCAGCTGTCCAAG-3′5′-TGTCGGAGCTCCTTCTTCTG-3′19611-β -Hydroxylsteroid dehydrogenase type 2 (11β HSD)U224245′-ATGGCATTGCCTGACCTTAG-3′5′-CTCAGTGCTCGGGGTAGAAG-3′194Vascular α-actin (VaACTIN)X068015′-GACACCAGGGAGTGATGGTT-3′5′-GTTAGCAAGGTCGGATGCTC-3′202Cyclophilin BAF0712255′-CAAGCCACTGAAGGATGTCA-3′5′-AAAATCAGGCCTGTGGAATG-3′239Cytochrome P450 subfamily XVII (Cyp17)M212085′-AAGTGGATCCTGGCTTTCCT-3′5′-CAATGCTGGAGTCGACGTTA-3′211AA924771 EST Rattus norvegicusAA9247725′-TTTGCTGTGCATGGGATTTA-3′5′-CCCTGCAGGATGTGAGAAGT-3′202aFrom GenBank (2004).Data analysis.We addressed potential interindividual variability by using independent samples of each experimental group (n = 5 for each set) for analysis. For the uterine/ovarian weight determination, the luminal epithelial cell height, and the gene expression analysis, we compared the data from the animals fed with the casein-based diet with the data from the animals fed the normal rodent diet (Purina 5001). For gene expression analysis, scanned output files of Affymetrix micro-arrays were visually inspected for hybridization artifacts and then analyzed using Affymetrix Microarray Suite (version 5.0) and Data Mining Tool (version 3.0) software, as described by the manufacturer (Affymetrix 2002; Lockhart et al. 1996). Arrays were scaled to an average intensity of 1,500 units and analyzed independently. The Affymetrix Rat Genome U34A microarrays used in this study have 8,740 probe sets corresponding to approximately 7,000 annotated rat genes and 1,740 expressed sequence tags (ESTs).For each transcript in the diet and dose groups, we conducted pairwise comparisons with vehicle controls fed the casein-based diet, using two-sample t-tests: first, we compared the two diet groups, and then we compared each treatment group with its respective diet control. We then conducted analysis of variance (ANOVA) for general diet and treatment effects on the signal value (which serves as a relative indicator of the level of expression of a transcript) and the log of the signal value. General diet effects were evaluated by ANOVA and a nonparametric test for dose–response trend, the Jonkheere-Terpstra test. Genes for which any of the tests had p ≤ 0.001 was taken as evidence that the expression of those genes was modified by the diet or by EE exposure. For the combined analysis of the two sets (casein-based or Purina 5001 diet), stratified nonparametric tests were conducted that were focused in detecting genes showing a diet response, or where there was a consistent treatment effect versus vehicle for the EE-treated group (0.1 μg kg/day). Here, we used linear models, with terms for both study and treatment effects, on average differences (signal values) and their log transformation, as well as stratified forms of the Wilcoxon-Mann-Whitney nonparametric statistic and a stratified form of the Jonkheere-Terpstra nonparametric statistic for diet response. Fold-change summary values for genes were calculated as a signed ratio of mean signal values (for each diet and EE-treated group compared with the appropriate control). Because fold-change values can become artificially large or undefined when mean signal values approach zero, all the values < 100 were made equal to 100 before calculating the mean signal values that are used in the fold-change calculation. All statistical analyses use the measured signal values, even if they were smaller than 100 units.ResultsEffect of diet on uterine/ovarian and uterine wet weight and uterine luminal epithelial cell height.Both diets, Purina 5001 and casein-based 5K96, were well tolerated by all the animals. We observed no evidence of overt toxicity and no clinical signs of toxicity. No difference was determined in body weights between animals fed either diet (Table 2). We did not detect premature vaginal opening in any of the animals in either diet group or in animals exposed to EE. There were no differences in wet uterine weight or in absolute and relative uterine weight (Table 2) between the two diet groups, even when the animals were exposed to low doses of EE.Table 2 Diet effect on body, uterine, and ovarian weight and luminal epithelial cell height of the juvenile (PND24) rat.Casein-based diet (5K96) Purina 5001 diet Body weight (g)Ovarian weight (mg)Uterine weight (mg)Epithelial cell height (μm)Body weight (g)Ovarian weight (mg)Uterine weight (mg)Epithelial cell height (μm)Peanut oil Mean ± SD (absolute)68.1 ± 4.832.0 ± 2.656.1 ± 8.213.3 ± 1.370.1 ± 4.934.8 ± 3.359.6 ± 10.814.0 ± 2.2 Mean ± SD (relative)a0.50 ± 0.040.58 ± 0.040.49 ± 0.070.56 ± 0.08 p-Valueb0.050.180.410.260.1 EE (μg/kg/day) Mean ± SD (absolute)68.5 ± 5.436.2 ± 2.161.9 ± 11.213.1 ± 1.671.1 ± 5.837.1 ± 1.866.4 ± 13.214.5 ± 1.3 Mean ± SD (relative)a0.52 ± 0.10.93 ± 0.10.52 ± 0.10.93 ± 0.2 p-Valueb0.050.140.360.18During the experimental phase, PND20 female rats were fed with a standard laboratory rodent diet (Purina 5001) or with a soy- and alfalfa-free diet (casein-based diet, 5K96) for 5 days (from PND20 to PND24). Epithelial cell height values were obtained from tissue sections from the midregion of each uterine horn, at equivalent areas, and with clear representation of the epithelium lining the lumen along the uterus (as shown in Figure 1). Epithelial cell height was determined by obtaining five measurements from five areas from two animals for each group. These values were used to determine the mean cell height SD for each treatment group, and the corresponding p-value.aRelative weight (mg/g body weight).bTwo-tailed t-test comparing 5K96 with Purina 5001, in control or treated animals; n = 15 for each diet group (controls) and n = 10 for EE-treated groups.The gross anatomy of the uterus and ovaries of animals fed either diet was identical, and no signs of accumulation of fluid in the uterine lumen were noted in any of the animals. We observed no differences in uterine weight gain (wet weight) or uterine epithelial cell height (Figure 1), and we found no change in the number of uterine glands. The classical morphologic changes induced by estrogen stimulation (hypertrophy of luminal epithelial, stromal, and myometrial cells; thickening of stromal layer; and some stromal inflammatory reaction) were not observed in any of the animals exposed to the two different diets, even when exposed to 0.1 μg kg/day EE (Figure 1).Figure 1 Representative uterine transversal sections from equivalent regions of vehicle-treated control immature rats (PND24; A, B, E, F) or animals treated with 0.1 μg/kg/day EE (C, D, G, H) fed with a casein-based diet (5k96; A, C, E, G) or a standard rodent diet (Purina 5001; B, D, F, H). Abbreviations: g, gland; Le, uterine lumen. See "Materials and Methods" for details. The rodent diet (B, D, F, H) containing quantifiable amounts of phytoestrogens did not have an impact on the histologic characteristics of the uterus, compared with tissues obtained from animals fed a relatively phytoestrogen-free casein-based diet (A, C, E, G). Bar = 0.08 mm for A–D; bar = 0.01 mm for E–H.Effect of diet on gene expression profile of the uterus/ovaries.In order to compare the gene expression profiles induced by the different diets (different phytoestrogen content) and the EE dose tested, we compared the average value of the signal values, a relative indicator of the level of expression of a transcript, between the two groups of independent controls. We then compared the appropriate diet-control group with the respective EE group (0.1 μg/kg/day), for all the 8,740 transcripts represented on the array.In comparing the expression profile identified in the uterus/ovaries of animals fed a casein-based diet versus the ones fed a soy/alfalfa-containing diet, we identified the expression of 29 genes that were significantly different (p ≤0.001). A list of those genes, along with their accession numbers, gene symbols, and the average fold changes, is shown in Table 3. The number of genes whose expression is modified by the diet's composition is relatively small, and the average fold change on the expression of these genes affected by the rodent standard diet, compared with the casein-based diet, is relatively low in the uterus and ovaries. Although robust expression differences for specific genes can be attributed to the composition of the diet, this list does not include genes well known to be estrogen regulated, such as progesterone receptor (PgR), intestinal calcium-binding protein (icabp), and complement component 3 (CC3).Table 3 Genes whose expression is modified by exposure to diet in the uterus/ovaries of the immature rat.GenBank accession no.aGene nameGene symbolAverage fold changebp-ValuecX67948Aquaporin 1 (aquaporin channel forming integral protein)AQP11.60.000159U56839Purinergic receptor P2Y, G-protein coupled 2P2ry21.40.000448AF017756GSK-3beta interacting protein rAxinAxin1.40.000130AA859529Diacylglycerol acyltransferaseDgat1.30.000470L06096Inositol 1,4,5-triphosphate receptor 3Itpr31.30.000420U90887Arginase type IIArg21.30.000728U78977ATPase, class II, type 9AAtp9a1.30.000022AA892562EST196365, high homology to nucleolar protein NAP57 and dyskeratosis congenita 1, dyskerinDkc11.30.000747AI639534ESTs, similar to properdin (factor P)1.30.000446AI231213ESTs, high homology to kangai 1 (suppression of tumorigenicity 6), prostateKai11.20.000561D10874Vacuolar H(+)-transporting ATPase,1.20.000865X56133Mitochondrial H+-ATP synthase alpha subunitAtp5a1−1.10.000854D13417Transcription factor HES-1 homolog of hairy and enhancer of split 1, (Drosophila)Hes1−1.20.000045Z71925Polymerase (RNA) II (DNA directed) polypeptide GPolr2g−1.20.000379AA818487ESTs, high homology to cyclophilin BPpib−1.20.000253AI112237ESTs, moderately similar to JE0384 NADH dehydrogenase−1.20.000192AA818858Peptidylprolyl isomerase A (cyclophilin A)Ppia−1.30.000943AA686579ESTs, similar to ubiquitin-like protein SMT3C precursor−1.30.000954U64705Protein synthesis initiation factor 4AII gene and E3 small nucleolar RNA gene−1.30.000405S69316GRP94/endoplasmin (5 and 3 regions)−1.30.000120M15481Insulin-like growth factor 1IGF-1−1.30.000068S69315GRP94/endoplasmin (5 and 3 regions)−1.40.000174D173103-alpha-Hydroxysteroid dehydrogenase (3-alpha-HSD)−1.40.000397X67859Autoantigen or Sjogren syndrome antigen BSsb−1.40.000103AA685903ESTs, similar to glucose regulated protein, 94 kDaGRP94−1.50.000878S68578Gonadotropin-releasing hormone receptorGrhr−1.50.000322AI009141EST203592, Rattus norvegicus−1.80.000468AF076619Growth factor receptor bound protein 14 or molecular adapter rGrb14 (Grb14), an inhibitor of insulin actionsGrb14−2.10.000158aFrom GenBank (2004).bThe average fold change was determined by comparing the average signal value of the indicated transcripts obtained from the uterus/ovaries from five females fed with the casein-based diet (5K96) versus the average signal value obtained from the same tissues from five females fed the standard rodent diet (Purina 5001).cTranscripts listed are those showing a robust response to the different diet (p < 0.001) using the stratified form of the Jonkheere-Terpstra nonparametric statistic to identify the diet response.One hypothesis is that if the soy/alfalfa-based diet was not estrogenic on its own, perhaps it would have sufficient potency to measurably enhance the effect of a sub-uterotrophic dose of EE. Although the expression of most genes from the prepubertal uterus/ovaries that respond to estrogen exposure is not altered by the diet composition, there are some that show a variable, nonstatistically significant response. For comparison, we calculated the relative fold change induced by diet for genes that showed a clear dose response to 1–10 μg kg/day EE (Naciff et al. 2003). Presumably, those genes have the potential to represent the response to weak or low levels of estrogen stimulation (expected from the dietary phytoestrogens) and are shown in Table 4. The fold change represents the ratio of the relative expression level of each gene in tissues from animals fed Purina 5001 versus those fed the 5K96 diet (as indicated in Table 3). For comparison, in Table 4 the relative expression level of the same transcripts under EE exposure is also shown. Analyzing the effect of EE exposure on the expression of the same set of genes, comparing the relative expression level of each gene in the tissues from animals exposed to EE versus their respective controls that were fed the same diet, and taking into account that lack of statistical significance for those genes listed, the average (n = 5) response to EE exposure is very similar, if not equal, even for those showing a relative large fold change, regardless of the diet fed to the animals. These results suggest that the response to the diet's composition is independent from the EE effect at this dose level of exposure. At higher doses of EE, the contribution of the dietary phytoestrogens is considered negligible because EE is a potent ER agonist able to interact with both isoforms of this receptor and with higher affinity than any of the dietary phytoestrogens (Kuiper et al. 1997). The changes induced by exposure to higher EE doses have been reported (Naciff et al. 2003). Corroboration of the microarray results by QRT-PCR for a selected group of genes is shown in Table 5. With the exception of CC3, which is undetectable by microarray analysis of samples derived from the animals exposed to two different diets (Table 5), the expression levels of the other genes are very similar, determined by either QRT-PCR or microarray analysis.Table 4 Diet effect on genes whose expression is modified by exposure to of 0.1 μg/kg EE in the uterus/ovaries of the immature rat.Average fold changebGenBank accession no.aGene nameGene symbolPurina 5001/5K96EE vs. control, 5K96EE vs. control, Purina 5001M29866Complement component 3CC3A14.77.0Y08358Eotaxin or small inducible cytokine A11Scya11A2.72.9AI013389ESTs, similar to calcium-binding protein, intestinal, vitamin D-dependentCalb32.02.41.9K00994Intestinal calcium-binding proteinicabp1.15.22.1U49062CD24 antigenCd241.11.31.2L14004Polymeric immunoglobulin receptorpigr1.51.41.1AA859661ESTs, similar to glutaminyl-peptide cyclotransferase precursorA2.11.7M57718Cytochrome P450 IV A1CYP4A1AAAU22424Hydroxysteroid dehydrogenase, 11-βtype 2Hsd11b21.02.21.6L07114Apolipoprotein B editing proteinApobec1AAAS79730Opioid receptor-like ORL1 receptorOprl11.21.31.4M88469f-SpondinSponf1.71.71.1X66845Dynein, cytoplasmic, intermediate chain 1Dncic1AAAL46593Small proline-rich protein geneSprr2.42.5−1.0L00191Fibronectin, encoding three mRNAs, exons 1, 2, 3fn−1.21.61.2M22323Gamma-enteric smooth muscle actinActg21.41.51.3D15069AdrenomedullinAdm1.92.41.2AA893870EST197673 Rattus norvegicusA1.72.0AI232078Transforming growth factor-β (TGF-β) masking proteinLtbp1−1.11.21.1U82612Fibronectin (fn-1) genefn-11.51.61.1X05834Fibronectin (fn-3) genefn-31.01.41.2L00382Skeletal muscle β-tropomyosin and fibroblast tropomyosin 1tpm11.21.71.3AA800908EST190405 Rattus norvegicus1.21.61.4M25758Phosphatidylinositol transfer proteinPitpn1.11.31.3AA799773ESTs, Rattus norvegicus1.31.31.2AA892829EST, similar to mouse bifunctional 3'-phosphoadenosine (PPS1)PPS11.21.31.1AB010963Potassium large conductance calcium-activated channelKcnmb11.21.41.2AF083269Actin-related protein complex 1bArpc1b1.01.3−1.0AA891760EST195563 Rattus norvegicusAAAAJ005394Collagen α1 type VCol5a1−1.11.61.2L11930Cyclase-associated protein homologueCap11.11.31.2X07467Glucose-6-phosphate dehydrogenaseG6pd1.21.2−1.1U26310TensinTns1.21.31.1AA891542EST195345 Rattus norvegicus, similar to mouse heat shock protein hsp40-3Dnajb51.31.31.1U44948Cysteine-rich protein 2 or smooth muscle cell LIM protein (SmLIM)Csrp2−1.1−1.6−1.4S61868Ryudocan or heparan sulfate proteoglycan core protein or syndecan-4SDC4−1.2−1.5−1.1L41254Corticosteroid-induced protein or FXYD domain-containing ion transport regulator 4Fxyd4−1.1−1.4−1.1AF023087Nerve growth factor induced factor A, or early growth response 1Egr1−1.5−1.71.1U07181Lactate dehydrogenase BLdhb−1.3−1.2−1.1X89225Solute carrier family 3, member 2Slc3a2−1.3−1.5−1.3AF054826Vesicle-associated membrane protein 5Vamp5−1.2−1.8−1.2X75253Phosp
Referência(s)