Adverse effects of subchronic exposure to cooking oil fumes on the gonads and the GPR30‑mediated signaling pathway in female rats

Background Cooking oil fumes (COFs) are composed of particulate matter, polycyclic aromatic hydrocarbons, volatile organic compounds, aldehydes, and ketones, and are currently a global health concern. Some agents in COFs are mutagenic and carcinogenic. However, only a few reports have addressed the hazardous effects of COF exposure on the female reproductive system. In this study, we explored the effects of subchronic exposure to COFs on female gonads in vivo and the possible involvement of the G-protein-coupled receptor 30 signaling pathway. Methods COFs were generated by heating commercially available canola oil in an iron pot. Adult female Wistar rats at 2 months of age were exposed to COFs at 32 mg/m 3 for 0, 0.5, 1, 2, or 4 h/day for 56 days. The estrous cycle in rats was studied twice at 7:00 a.m. and 7:00 p.m. on the 43rd treatment day until the current estrous cycle was complete. The rat body weight was measured before the experiment and at day 56 post-exposure. At the end of the experiment, rat blood was collected for gonadal hormone assay, and ovaries were collected for histology and mRNA isolation. The mRNA levels of GPR30, EGFR, STAT3, and ERK were determined by quantitative RT-PCR. Results At a concentration of 32.21 ± 5.11 mg/m3, COF exposure extended the estrous cycle in rats, and ovary coefficient decreased. COFs showed various effects on the sex hormone levels and follicles, depending on its exposure level. Exposure to COFs led to the changes in mRNA levels of the G-protein-coupled receptor 30 (GPR30), epidermal growth factor receptor (EGFR), signal transducer and activator of transcription 3 (STAT3), and extracellular signal-regulated kinase (ERK). Conclusion This study indicated that cooking oil fume exposure disrupted the estrous cycle, sex hormone patterns, and follicle development in female rats in a dose-dependent manner. These adverse effects of cooking oil fumes on female reproductive health were correlated with the G-protein-coupled receptor 30-mediated signaling pathway.


Introduction
Cooking oil fumes (COFs) are pollutant mixtures of particulate matter, polycyclic aromatic hydrocarbons, volatile organic compounds, aldehydes, and ketones (Huang et al. 2011;Zhu et al. 2001). Our previous studies showed that different cooking oils and cooking methods contain alkanes and antioxidants (Youqiong et al. 2012). COFs are emitted from domestic and commercial cooking activities. People worldwide are exposed to COFs almost every day, and COFs are currently a global health concern because of their multiple health effects.
COFs are mutagenic and carcinogenic agents (Ke et al. 2009;Lee and Gany 2013;Wu et al. 1999). Furthermore, the complex components of COFs may exert adverse effects on the reproductive system. Benzo(a)pyrene (BaP), a vital chemical component of COFs, has generated a substantial amount of interest, because it can disrupt gonad and endocrine function. The ovaries contain the necessary enzymes for BaP metabolism, fostering its reaction with reactive metabolites and disrupting normal ovarian functions (Archibong et al. 2012). Long-term, low-dose formaldehyde (one of COF components) exposure in rats substantially reduced the number and size of mature follicles, and aggravated vascular congestion and interstitial edema in their ovaries . Acetaldehyde (one of COF components) was reported to exert toxic effects, impairing the differentiation of granulosa cells and reducing ovulation and oocyte quality (Kawai et al. 2012).
Certain chemicals in COFs are endocrine disruptors, even at low doses, that can possibly cause vital changes in cell functions, such as non-genomic effects. Tectoridin, a major isoflavone isolated from the rhizome of Belamcanda chinensis and a known phytoestrogen, exerts its estrogenic effect mainly via the G protein-coupled receptor 30 (GPR30) and the extracellular signal-regulated kinase (ERK)-mediated rapid non-genomic estrogen signaling pathways. This endocrine disruption property of tectoridin (which is observed in other phytoestrogens, including genistein, bisphenol A, nonylphenol, and kepone) differentiates it from genistein. The phytoestrogen genistein exerts estrogenic effects via an estrogen receptor-dependent genomic pathway and a GPR30-dependent non-genomic pathway (Kang et al. 2009; Thomas and Dong 2006). G protein-coupled receptors (GPRs) bind to estrogen and participate in the rapid non-genomic signaling events that are widely observed following the stimulation of cells and tissues with estrogen. GPR30 is responsible for the receptor-mediated estrogeninduced activation of ERK, via the transactivation of the epidermal growth factor receptor (EGFR) and the release of heparan-bound epidermal growth factor (EGF) (Prossnitz and Barton 2014). The engagement of surface EGFR by cognate ligands drives the survival and proliferation of follicular dendritic cells through signaling to the nucleus, mainly via the mitogen-activated protein kinase and signal transducer and activator of transcription (STAT) pathways (Vermi et al. 2013).
Frying is a common activity, and females are more prone to COF exposure. However, only a few reports have addressed the hazardous effects of COF exposure on the female reproductive system. We studied these effects in the present study and hypothesized that the effects of COFs on the gonads are possibly related to the GPR30-mediated signaling pathway. The effects of subchronic COF exposure on female gonads in vivo were explored by analyzing changes in the estrous cycle, gonad weight, serum gonadal hormone levels, and follicle type. The possible effects of COFs on the GPR30 signaling pathway were studied by examining the mRNA levels of GPR30, EGFR, STAT3, and ERK in the ovarian tissue.

Generating COFs and monitoring
Pressed canola oil was bought from a local market, and 100 g was poured into an iron pot (diameter 45 cm, made of cast iron, without non-stick paint) and heated in a thermostatic oven placed atop a table at the center of a chamber (length × width × height: 1.5 m × 1.5 m × 3 m). The temperature was controlled at 260 ± 5 °C. A light ceiling fan was used to ensure proper circulation of the fumes in the chamber.
COFs were produced at a rate of approximately 218.4 mg/ min under the aforementioned conditions in our pilot experiments. The ventilation rate in the chamber was 6.75 m 3 / min, and the setting was designed with an exhaust fan, such that the COF level was constant at approximately 32 mg/m 3 while the rats were treated and the exhaust was on. Generally, the concentration reached a plateau 10 min after heating. The dust sampler with glass fiber filter membrane and the atmospheric sampler equipped with activated carbon tube were used to collect fume from the 30 cm of the oil smoke generator. The discharge of the dust sampler was set to 15 L/min and the collection time was set to 20 min. The flow of the volatile organic compounds in the cooking oil smoke by the activated carbon tube was set to be 0.1 L/min, and the sampling time was 30 min. After weighing, the fiber glass filters were detected by high-performance liquid chromatography (HPLC), and activated carbon tubes were used for GC-MS analysis. The COF level was estimated as the quotient of the increased weight of the filters and the sampling volume of air (mg/m 3 ).

3
With reference to the EPA Method Update TO-17 1667 of the US Environmental Protection Agency (USEPA 2009) and a related report, major COF compounds were analyzed. The chamber air was sampled once a week using a conditioned Tenax adsorption tube and pump at a flow rate of 15 mL/min for 4 h. During the GC-MS analysis, the sample was desorbed with an inert gas at a flow rate of 20 mL/ min for 1 min at 150 °C and then at 50 mL/min for 5 min at 250 °C. The desorbed gas was cooled at a cold trap temperature (− 10 °C), injected at a flow rate of 50 mL/min into the GC-MS system (Agilent 7890 A/Agilent 5975C, US), and analyzed under the following operating conditions: DB-5MS column (30 m × 320 μm × 1.00 μm); helium gas: flow rate of 1.1 mL/min; injector temperature: 250 °C; temperature program: the oven temperature was held at 40 °C for 5 min, programmed to 150 °C at 6 °C/min, held for another 5 min, programmed to 250 °C at 5 °C/min, and held for 5 min; interface temperature: 280 °C; ion source temperature: 150 °C; electron multiplier voltage: 1329 V. Constant current: 1.5 mL/min; and range of ion scanning: 35-300 U, 0-4-min solvent delay. Some COF compounds were identified by comparing the retention times of the peaks with those of the reference compounds eluted under the same chromatographic condition. The peak areas of the internal standards were employed for performing quantitative tests of the specific components of COFs under the experimental exposure scenario. The other COF components were identified by comparing their spectra with those available in the National Institute of Standards and Technology (NIST) digital library (MS Search 2.0 database) and performing map analysis according to a matching factor of more than 850%. Then, with reference of the designed mass and peak area of methylbenzene, an area-normalized semi-quantitative method was used for quantitative tests of those components of COFs under the experimental exposure scenario.

Animal treatment
Adult female Wistar rats at 2 months of age were randomly grouped into five polyethylene cages (n = 10 per treatment or control group). Rats were acclimatized to the animal care facility for 1 week prior to the initiation of experiments. They were maintained in an environmentally controlled room with natural lighting at 22 °C and 50-60% humidity. Rats had ad libitum access to commercial chow and water. Acclimated rats were weighed, and their ability to maintain regular 4-5 day estrous cycles was evaluated through daily vaginal histological examinations, which entailed using light microscopy to determine the stages of the estrous cycle (Goldman et al. 2007).
The treatment entailed subchronic COF exposure (through inhalation) at 32 mg/m 3 for 0, 0.5, 1, 2, or 4 h/day, and this exposure lasted continuously for 56 days (Archibong et al. 2012). The rats were washed every day with detergent and warm clean water after COF exposure and transferred to clean cages and cells, so that they received no additional exposure. This study was approved by the animal care and use committee of our institute, and all animals were used in compliance with the National Institutes of Health guide for the care and use of laboratory animals.

Estrous cycle evaluation
The estrous cycle in rats was studied twice at 7:00 a.m. and 7:00 p.m. on the 43rd treatment day until the current estrous cycle was complete (Goldman et al. 2007).

Body and gonad weight
The individual body weights of the rats were measured using an electronic balance at the beginning and on the 56th exposure day (the final day) of the treatment. The rats were then sacrificed at the estrus stage of the first cycle within 12 h after treatment. The ovaries and uterus were harvested, washed in chilled (4 °C) isotonic saline to remove the blood on the surface, and weighed. The organ-weight-to-bodyweight ratio was computed. The time variance between the final exposure and euthanasia or tissue collection was 38.7 ± 43.9 h among the rats.

Gonadal hormone assay
Blood was sampled via a heart puncture using heparinized Pasteur pipettes. Serum was harvested from each blood sample through centrifugation at 1500×g at 4 °C for 10 min and stored at − 20 °C until ovarian steroids [estradiol (E 2 ) and progesterone (P 4 )] and gonadotropins [luteinizing hormone (LH) and follicle-stimulating hormone (FSH)] were assayed. The hormone levels were measured using a radioimmunoassay kit (160220, Beijing North Institute of Biological Technology Co Ltd) on an automatic radioimmunity analyzer (SN-695B, Shanghai Hesuo Rihuan Photoelectric Instrument Co Ltd) according to the manufacturer's protocol. The sensitivity for E 2 , P 4 , LH, and FSH were 15 pg/mL, 0.3 ng/mL, 1 mIU/mL, and 1 mIU/mL, respectively, and the degree of accuracy was > 85% for all hormones.

Follicle categorization
Ovaries at the estrus stage were fixed in 4% paraformaldehyde at 4 °C, and the slides were stored overnight before being dehydrated in ethanol at a series of concentrations, rinsed in xylene, and embedded in paraffin. From these slides, 5-μm-thick ovarian sections were prepared, with a distance of greater than 20 μm between each section. Six representative sections from each ovary were selected and deparaffinized in xylene, dehydrated through a series of ethanol concentrations, and then stained with hematoxylin and eosin (H&E) using standard protocols. The number and percentage of follicles at the estrous stage in each section were histologically analyzed using 10× magnification and photographed using a photomicroscope (Olympus BX50). Two independent observers, who were blinded to the treatment assignment, conducted the histopathological examination. Follicles were classified (Luo et al. 2008) as follows: primordial follicle (an oocyte surrounded by one layer of flattened granulosa cells), primary follicle (an oocyte surrounded by one layer of cuboidal granulosa cells), secondary follicle (two or three layers of cuboidal granulosa cells with no antral space), or antral follicle (more than four layers of granulosa cells with one or more independent antral spaces or with a cumulus granulosa cell layer). Antral follicles were considered atretic if they contained at least 20 apoptotic granulosa cells, disorganized granulosa cells, a degenerating oocyte, or a fragmentation of the oocyte nucleus (Juliani et al. 2008).
Pieces ≤ 0.5-1 mm of fresh ovary tissue samples were prefixed in a stationary liquid of 3% glutaraldehyde-paraformaldehyde and were then post-fixed with 1% osmium tetroxide and 1.5% potassium ferrocyanide (Sigma). The samples were dehydrated with increasing concentrations of alcohol and finally with acetone. The specimens were embedded in Araldite resin 618. The follicles were examined with transmission electron microscopy (Philips EM208).

Quantitative RT-PCR assay for determining the mRNA levels of GPR30, EGFR, STAT3, and ERK
Total RNA was extracted from the ovary tissue using a TRIzol reagent (Invitrogen) according to the manufacturer's instructions. The concentration and purity of the isolated total RNA was determined spectrophotometrically at 260 nm and 280 nm using a spectrophotometer (Thermo Fisher Scientific). Using primers (Invitrogen), 1 μg of total RNA was reverse transcribed to cDNA using the GoScriptTM Reverse Transcription System (Progerm) according to the manufacturer's protocol. The target fragments were quantified through real-time PCR using the GoTaq ® Colorless Master Mix (Invitrogen) with the cDNA template. Each sample was tested in duplicate (Zhang et al. 2012). Specific primer sets (designed by Invitrogen) are described in Table 1. The transcript level of the housekeeping gene glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was used as an internal reference gene. PCR amplification efficiencies were determined by amplifying various quantities of target cDNA for each reaction (Filippi et al. 2009). The PCR products were separated using electrophoresis. Images were captured using a charge-coupled device camera, and the relative mRNA concentrations were analyzed using Quantity One 4.6.2 (Bio-Rad) and compared with the GAPDH level, which served as a Kong et al. (2011).

Statistical analysis
The normal distribution of COF concentrations was evaluated using the one-sample KS test. All differences in laboratory measurements among the dose groups were analyzed using one-way analysis of variance (ANOVA). Data are expressed as the mean and standard deviation ( x ± s).

Airborne COFs in the treatment chamber
There were 56 airborne COF samples used in the treatment for the concentration analysis. The COF concentrations were normally distributed. The arithmetic mean was 32.21 mg/m 3 with a standard deviation of 5.11 mg/m 3 . The GC-MS and HPLC analysis results (Table 2) indicated that aldehydes and alkyls were the main chemicals in the COFs. BaP was detectable at a high concentration under the experimental conditions. Moreover, antioxidant additives, such as butylated hydroxytoluene and tertiary butylhydroquinone, were found in the COFs.

Effects of COFs on the estrous cycle
The rats in the control group exhibited a normal 5-day estrous cycle (Table 3). The estrous cycles of all rats exposed to COFs were extended compared with those of the control group; however, this difference was not significant. Subchronic (43 days) exposure of rats to a high dose (2 h/day) and a very high dose (4 h/day) of COFs significantly prolonged the metaestrous stage (p < 0.05 and p < 0.01, respectively).

Effects of COFs on body and gonad weight
No remarkable differences in activities, including eating or drinking, were observed in the rats of the control or treatment groups during the treatment (Table 4). No significant difference was observed in the body weight among the groups before treatment. The rats in all groups had increased  body weights; however, the rats in the very high COFexposed groups weighed less than those in the control group. COF exposure led to a significant decrease in the absolute weight of the rat ovaries in very high group (p < 0.05). Moreover, the relative weights of the ovaries decreased in the rats in the very high-exposure groups. COF exposure (at all doses) did not exert any effects on the uterus weight.

Effects of COFs on gonadal hormones
Irrespective of the dose, subchronic COF exposure did not affect serum FSH levels in the female rats at the estrus stage (Table 5). The LH levels of the rats increased in the high-dose and very-high-dose group, but the E 2 levels of the rats decreased in the very-high-dose group (p < 0.05/p < 0.01). The P 4 levels of the rats decreased in the very-high group (p < 0.01).

Effects of COFs on follicle type
To investigate whether exposure to COFs caused direct or indirect damages to the follicle structure and type in ovaries, ovarian tissues were observed under the light mirror with H&E staining. Exposure to COFs resulted in increased follicular atresia and decreased mature follicles, while primordial/primary follicles were more observed in the control  Fig. 1 Optical mirror images of ovarian. Rats were inhaled with normal air (a), low dose (b), medium dose (c), high dose (d), and very high dose (e) of COFs for 56 days. Ovarian structure was observed under the microscope with HE staining. Primordial/primary follicle is more seen in a (blue arrows). Secondary follicles are more seen in b (green arrows). Corpus lutea are more seen in c (the picture is corpus luteum, blue *). Graafian follicles are rarely seen in d (blue #). Atretic follicle is more seen in e (red arrows). Scale bars indicate 5 μm rats ( Fig. 1), indicating that COFs can affect follicular development. The numbers of follicles in different grades were compared with those in ovaries. The COF treatment showed diverse effects on the percentage of follicle types. Compared to the control group, low COF exposure led to an increase in the Secondary follicle ratio during the estrous stage, but no such change resulted from high COF exposure (Table 6). Moreover, the atretic follicles increased in number (doubled or more) after COF exposure in high-dose groups. Even the ratio of the corpus luteum slightly increased in the rats of the very-high-dose group; however, no such change was observed in any other groups. The Graafian follicle

Effects of COFs on follicular cell ultrastructure
To investigate whether exposure to COFs caused direct or indirect damages to the follicle structure in the ovaries, ovarian tissues were observed under TEM. Exposure to COFs resulted in follicular cell damage and apoptosis in ovaries, while no pathological abnormality was observed in the control rats (Fig. 2), indicating that COFs can affect follicular structure.

The mRNA levels of GPR30, EGFR, STAT3, and ERK
To investigate whether the signal pathway changes in endocrine disruption are involved in COF-induced apoptosis and inhibition of follicular growth, the expression levels of GPR30, EGFR, STAT3, and ERK mRNA were measured in ovarian tissue from rats with or without exposure to COFs (Fig. 3). The results showed that GPR30, EGFR, STAT3, and ERK mRNA were expressed in each group. There were statistically significant differences in the mRNA expression levels of GPR30, EGFR, STAT3, and ERK in each group (F = 18.861, p < 0.01; F = 14.040, p < 0.01; F = 35.431, p < 0.01; and F = 30.599, p < 0.01). The mRNA expression levels of GPR30 and EGFR decreased significantly in the high-and very-high-dose group (p < 0.01). The expression levels of STAT3 mRNA in the ovaries of rats treated with COFs (low dose, moderate dose, high dose, and very high dose) were significantly higher than that of the control rats, and the expression level of ERK mRNA increased significantly in the ovaries of rats treated with COFs (moderate dose, high dose, and very high dose) (p < 0.01).

Discussion
The present in vivo study found that COF exposure decreased the overall weight, including that of the ovaries, interrupted follicle development, hormone-level disorder, and increased the STAT3 and ERK mRNA levels in the ovaries at the estrous stage in female rats. However, the GRP30 and EGRF mRNA levels fluctuated with the COF dose at this stage.

Exposure levels
During the subchronic exposure experiment, the COF concentration was constant at 32.21 ± 5.11 mg/m 3 . This observation was similar to that of Wu et al. in which the COF concentration was reported to be 25.1 ± 2.3 mg/m 3 , with lard heated at 250 ± 10 °C for 30 min in an iron pot (Wu et al. 1999). The COF concentration in our study was comparable to that of pan-fried beef in the western kitchen model (19.5 mg/m 3 and 42.8 mg/m 3 COFs generated using an electric hob or a gas hob, respectively) (Svedahl et al. 2009). However, the COF concentration in our study was much higher than those from 19 restaurants in Norway (the highest was 6.6 mg/m 3 in a small local restaurant, and the arithmetic mean for COF concentration for all of the kitchens was 0.62 mg/m 3 ) (Svendsen et al. 2002). The chemical compositions of COFs are complicated and differ according to the oil type, food type, cooking method, and heating temperature used. Under the experimental conditions, we observed that aldehydes and alkyls were the main components in the COFs, and this observation was similar to the observations of a previous study . Furthermore, BaP was detected at a high concentration under the experimental conditions. Because the cooking duration is not more than 4 h Fig. 3 Impact of subchronic exposure to COFs on mRNA levels of GPR30, EGFR, STAT3, and ERK (n = 10, relative values as compared with GAPDH, ×10 −2 ). Rats were inhaled with normal air (0), low dose (0.5), medium dose (1.0), high dose (2.0), and very high dose (4.0) of COFs for 56 days. The mRNA expression levels of GPR30, EGFR, STAT3, and ERK in rat ovaries were analyzed by real-time PCR. Values are expressed as mean ± SD (n = 10). a p < 0.05 and b p < 0.01 vs. control in a typical family, the highest exposure time was set at 4 h/ day. The experiment was performed for 56 days, so that the subchronic health effects of COF exposure on gonads could be observed in the female rats.

COF exposure and the estrous cycle
The estrous cycle is generally assessed according to changes in the vaginal epithelial cell structure to document the effects on reproductive cycles and provide an index of the functional status of the hypothalamic-pituitary-ovarian axis (Goldman et al. 2007). In our study, the rats in the control group had normal estrous cycles, which have been reported previously. However, we observed that the extension of the estrous cycle was associated with the COF exposure levels in a dose-dependent manner. Although only the metaestrous stage was significantly extended in the high-and very-highdose groups, it seemed that the cycle extension was mainly caused by the extension of the diestrous, metaestrous, and estrous stages, as indicated in Table 3. This effect was possibly caused by the BaP chemicals of COFs. BaP exposure extended the length of the estrous cycle by approximately 24 h among rats exposed to BaP at 100 μg/m 3 (Archibong et al. 2012). This observation indicated an imbalance in the ovarian steroids imposed by BaP. Consequently, irregular changes in the length of the estrous cycle occurred at the proestrous stage of the estrous cycle, in which estrogen is maximally secreted by the ovarian follicles (Archibong et al. 2012). A significant decrease in E 2 levels was observed at the proestrous stage among BaP-exposed rats. In our study, the E 2 levels decreased in the very-high-dose group. Consequently, it was hypothesized that the extension of the estrous cycle duration in the COF-exposed rats was caused by the reduction in the level of circulating E 2 , which is required to trigger the onset of the estrous stage in a timely manner, and the delay in attaining the serum E 2 peak (Goldman et al. 2007).

COF exposure and sex hormones
The levels of four hormones, including follicle-stimulating hormone (FSH), luteinizing hormone (LH), estrogen (E 2 ), and progestin (P 4 ), are considered as an important evaluation index of the influence of foreign factors on the estrous cycle and follicle development. The secretion of reproductive hormones (P 4 and E 2 ) from the ovaries is regulated by the release of LH and FSH from the anterior pituitary gland. Physiologically, increased P 4 and E 2 levels lead to decreased FSH and LH release via a negative feedback mechanism on the hypothalamic-anterior-pituitary axis. However, this typical feedback signaling was not observed in all time variance and exposure groups, which was possibly related to the endocrine disruption effects of the COFs. We observed that LH levels increased in the high-dose and very-high-dose group, whereas E 2 levels decreased in the high-dose and very-high-dose group. P 4 in the high-dose group decreased significantly compared to the control group. This was possibly because COF exposure interrupted the hormonal patterns of the rats, as some studies have suggested (Archibong et al. 2012;McLean et al. 1977). Serum concentrations of E 2 and LH, for example, decreased at the proestrous stage in BaPexposed rats, whereas those of P 4 decreased at the diestrus I stage (Archibong et al. 2012). The increase in the concentrations of LH observed at the proestrous stage in the control rats free of BaP exposure suggested that a surge in this gonadotropin was imminent (Archibong et al. 2012). BaP delays the surge of LH that is normally observed throughout the proestrous stage. However, the LH surge may have occurred much later in the BaP-exposed rats at the proestrous stage, because these rats ovulated, albeit producing fewer eggs compared with the unexposed control rats (Archibong et al. 2012;McLean et al. 1977).
In the present study, E 2 and P 4 levels decreased, while LH levels increased in the high-dose and very-high-dose group. The results indicated that high-dose exposure to COFs may inhibit the secretion of P 4 and E 2 in the ovaries, and stimulate secretion of LH in the pituitaries through negative feedback. In addition, LH turned the granular cells into luteal cells such that they secreted pregnendione, leading to increased P 4 levels. The results seemed that the sex hormones in COFs-exposed rats appeared in a pattern different from the normal one.

COF exposure and body weight, uterus, and ovaries
The results indicated that COF exposure increased the body weight compared with normal conditions. The curve of the increased body weights was consistent with that of the changes in the serum E 2 levels. Specifically, E 2 reached high levels in the low-dose group, and decreased slightly in the high-and very-high-dose groups. Conversely, the curve of body-weight gain varied with that of the serum LH levels. Estrogen promotes the accumulation of subcutaneous fat (Brown and Clegg 2010). Therefore, the body weight increased. By contrast, the anorexigenic effects of estradiol in adult female rats reduced their food intake during the estrus stage and led to a subsequent decrease in weight gain in the very-high-dose group (Eckel 2011).
The ovaries are a crucial target organ of many reproductive toxicants. A significant increase in the absolute weight of the fresh ovaries was observed in the rats in the lowand moderate-dose groups. However, a high-dose and very high-dose of COF exposure reduced the weight of the ovaries, although the difference was not significant; it suggests that high dose of COF is a direct cause of ovarian toxicity. The effects of COFs on the ovaries were exerted in a dose-dependent manner. Very high levels of COF exposure caused direct toxicity and reduced the ovary weight, and disturbed the development and functioning of the ovaries at the estrous stage. However, the uterus may be less sensitive to COF exposure than the ovaries, because no significant weight change in the rat uterus was observed in this study.

COF exposure and the follicle types and ultrastructure
Enumerating ovarian follicles is an efficient technique for estimating the extent of ovarian toxicity in female rats. Follicle development is regulated by sex hormones. Generally, increased FSH and E 2 levels improve follicle development and maturity in the estrus period of rats. However, COFs may disrupt this process and affect their regulation. Our present study showed that follicle development in the highdose COFs was influenced. The ratio of secondary follicles decreased in the very high-dose group, possibly because their development was inhibited. However, follicle aging or maturation was accelerated in all dose groups.
The LH and FSH levels changed differentially throughout the ovulation cycles, with FSH being predominant in the early follicular phase and LH being predominant in the late follicular phase. Sequential FSH and LH predominance is crucial for appropriate follicular maturation, sex steroid production, and subsequent ovulation (McCartney 2010;Zhou et al. 2013). E 2 modulates steroidogenesis, promotes granulosa cell proliferation, and regulates follicular development. The surge in LH concentration in mammals is a functional endocrine event, initiating the stages of follicular and oocyte maturation that culminate in the release of an ovum or multiple ova. The asynchrony between the developing follicles and preovulatory LH surge can be induced by COF exposure, which can disturb follicle development. Therefore, we observed ovarian structures using H&E staining and transmission electron microscopy to understand whether follicles in different stages had been destroyed. The morphological results showed that the numbers of follicle types were disrupted and that the structures were damaged by high-dose COF exposure, e.g., mitochondria were swollen and even vacuolated. The chromatin aggregation and more apoptotic bodies were present. Such as BaP was reported to reduce the secretory functions of the ovaries to interrupt the development of oocytes and cause atrophy of the ovaries (Archibong et al. 2012).
Ovarian follicular degeneration or atresia is a hormonally controlled apoptotic process (Kaipia and Hsueh 1997). The ultrastructural changes of ovarian cells in this experiment include swelling and edema of mitochondria, and the increase of lipid droplets and apoptotic bodies. Follicle atresia is mainly caused by the apoptosis of granulosa cells (GC). GC is an important cell for the synthesis and secretion of E 2 hormone. The apoptosis and damage of granulosa cells will lead to the decrease of E 2 hormone secretion. Steroid hormones are mainly produced in mitochondria, which are closely related to lipid droplets in cytoplasm. Lipid droplets may help to store the substrates of steroid hormones and transport them to mitochondria. A large number of lipid droplets and mitochondria damage were found in the highdose group, suggesting that COFs may damage the ovary by impairing the synthesis and transport of steroid hormones. Therefore, morphological changes in these cells can disrupt ovarian secretion, leading to abnormal follicular development and atresia, and the corresponding estrus cycle also changed, which is consistent with the results of this study.
COF exposure and the mRNA levels of GPR30, EGFR, STAT3, and ERK GPR30 is highly expressed in the granulosa cells and theca cells of ovaries, GPR30 may be involved in the rapid non-genomic effect of estrogen and genomic transcription reaction, and play an important role in cell proliferation and apoptosis. GPR30 binds to estrogen, and then rapidly activates adenylyl cyclase and mitogen-activated protein kinases, such as ERK, resulting in the mobilization of the intracellular Ca 2+ concentration in cells (Giuliano et al. 2013). GPR-mediated transactivation of EGFR by estrogen provides a cross-talk mechanism between estrogen and serum growth factors; that is, it induces the EGF-like effects of estrogen (Filardo 2002). STAT3 is a transcription factor activated by various cytokines and growth factors, including interferons and EGF, and regulates the genes involved in cell growth and division, cell proliferation, and apoptosis (Klampfer 2006).
Studies have shown that many environmental estrogens, which can be combined with estrogen receptor (ERα and ERβ), can be combined with estrogen G protein-coupled receptor (GPR30), such as bisphenol A (8-50 times), nonylphenol and open fluffy (3-4 times), and iris flavone glycosides. The affinity with GPR30 is higher than that of ERα and ERβ. The combination of GPR30 is mainly produced by the estrogen effect (O'Dowd et al. 1998;Chen and Wang 2012). It has been reported that methylbenzene, phenol, alkyl, and phthalate esters in COFs are toxic to female gonads. In our study, we observed that COF exposure significantly decreased the E 2 hormone level, GPR30, and EGFR expression levels of the rats in the high-and veryhigh-dose groups. GPR30 and EGFR were positively correlated with E 2 and corpus luteum ratios in the female rats at the estrous stage. E 2 , acting through ERα, GPR30, or both, was reported to activate the rapid EGFR/ERK/c-Jun pathway in rat pachytene spermatocytes (Chimento et al. 2010). Therefore, GPR30 might be involved in mediating the effect of COFs on E 2 hormone level. 1 3 COF exposure increased STAT3 and ERK mRNA expression and LH levels in a dose-dependent manner, while the ratio of the corpus luteum decreased significantly. Furthermore, STAT3 and ERK were positively correlated with the ratios of the atretic follicles. The results indicated that COF exposure at any level possibly accelerated the aging of the follicle via the STAT3 and ERK signaling pathways.

Conclusions
Subchronic exposure to COFs could disrupt the estrous cycle, sex hormone patterns, and follicle development in female rats, mainly by disturbing the function of the hypothalamic-pituitary-ovarian axis. The study explored the role of the GPR30-mediated signaling pathway in the adverse effects of COFs on the female gonads. We found that COFs interrupted the mRNA expression levels of the related molecular signals, and led to hormonal disorders, abnormal cell proliferation and ovarian injury in rats. The findings of this study are of value in understanding the mechanism of the health effects of COFs in humans. The results suggest that experimental studies for determining the health effects of COFs in vitro and epidemiologic studies may be necessary.