Abstract
The expression of type I GNRH receptor (GNRHR-I) and the direct role of GNRH-I on corpora lutea (CL) function were studied in the pseudopregnant rabbit model. Immunohistochemistry evidenced GNRHR-I and GNRH-I in luteal cells at early (day 4 pseudopregnancy)-, mid (day 9)-, and late (day 13)-luteal stages. Real-time RT-PCR and western blotting revealed GNRHR-I mRNA and protein at the three luteal stages. Buserelin in vivo treatment at days 9 and 13 decreased plasma progesterone levels for 48 and 24 h respectively. In in vitro cultured CL, buserelin reduced progesterone secretion, increased prostaglandin F2α (PGF2α) secretion and cyclo-oxygenase-2 (COX-2) and nitric oxide synthase (NOS) activities at days 9 and 13, and decreased PGE2 at day 13. Co-incubation with antagonists for GNRH-I (antide), inositol 1,4,5-trisphosphate (IP3, 2-amino-ethoxydiphenylborate), and diacylglycerol (DAG, 1-hexadecyl-2-acetyl glycerol) or inhibitors for phospholipase C (PLC, compound 48/80), and protein kinase C (PKC, staurosporine) counteracted the buserelin effects. Buserelin co-incubated with COX inhibitor (acetylsalicylic acid) increased progesterone and decreased PGF2α and NOS activity at days 9 and 13, whereas co-incubation with NOS inhibitor (N-nitro-l-arginine methyl ester) increased progesterone at the same luteal stages. These results suggest that GNRHR-I is constitutively expressed in rabbit CL independently of luteal stage, whereas GNRH-I down-regulates directly CL progesterone production via PGF2α at mid- and late-luteal stages of pseudopregnancy, utilizing its cognate type I receptor with a post-receptorial mechanism that involves PLC, IP3, DAG, PKC, COX-2, and NOS.
Introduction
GNRH is a hypothalamic neuronal secretory decapeptide that plays a central regulatory role in mammalian reproductive physiology. GNRH influences reproductive processes mainly by regulating pituitary gonadotropin synthesis and release, which in turn modulate steroidogenesis and gametogenesis (Conn & Crowley 1994). Even if the hypothalamus and pituitary are, respectively, the principal source of and target site for GNRH, several studies have reported an extrahypothalamic origin of GNRH, as well as an extrapituitary presence of GNRH receptor (GNRHR) in numerous peripheral tissues including reproductive organs, such as testis, prostate, ovary, oviduct, placenta, and mammary gland (Ramakrishnappa et al. 2005).
There are currently 14 identified naturally occurring GNRH structural variants across the vertebrate species (Somoza et al. 2002). However, it is evident that at least three forms of GNRH are present in the majority of the vertebrate species studied: hypothalamic GNRH (GNRH-I), mid-brain GNRH (GNRH-II), and telencephalic GNRH (GNRH-III) (Gorbman & Sower 2003). In extrapituitary reproductive tissues, GNRH is considered to act in an autocrine or a paracrine manner and to regulate ovarian steroidogenesis by exerting a stimulatory as well as an inhibitory effect on the production of steroid hormones and apoptosis in ovarian follicles and corpora lutea (CL; Dubois et al. 2002, Ramakrishnappa et al. 2005). Although three types of GNRHR have been identified thus far, mammals have only two types of receptors: GNRHR-I and -II (Millar 2005). All known GNRHRs are transmembrane G protein-coupled receptors that produce their effects by activating phospholipase A2 (PLA2), phospholipase C (PLC), PLD, or adenylate cyclase (AC) cell signaling pathways (Millar 2005). Phospholipase activation may generate arachidonic acid (AA) that is converted into prostaglandins (PGs) by the action of the cyclo-oxygenase-1 and -2 (COX-1 and -2) as well as by that of other PG synthase enzyme types (Naor 2009).
Various studies have shown that GNRH exerts both inhibitory and stimulatory effects in the gonads, with either down- or up-regulation of the ovarian cellular steroid production (Ramakrishnappa et al. 2005). Few studies have addressed the direct effects of GNRH in rabbit CL, with conflicting results: earlier works pointed to CL regression with serum progesterone decline (Hilliard et al. 1976, Rippel et al. 1976), while a later study indicated that there was no effect on steroid production by ovary cultured in vitro (Eisenberg et al. 1984). Therefore, the main objectives of the present study were to examine the in vivo and in vitro modulatory effects of GNRH on rabbit CL during early-, mid-, and late-luteal stages of pseudopregnancy. With this end in view, experiments were devised to clarify the presence of GNRHR-I and GNRH-I and their possible post-receptorial mechanisms, as well as the effects of GNRH-I on the production of progesterone, PGE2, and PGF2α and on the activities of COX-1 and -2 and nitric oxide synthase (NOS), the enzymes involved in the regulation of rabbit CL life span (Boiti et al. 2000, 2005, Zerani et al. 2007).
Materials and Methods
Reagents
Reagent for isolation of total RNA (TRIzol) was purchased from Invitrogen. iSCRIPT cDNA and iQ SYBR Green SuperMix were purchased from Bio-Rad. QIAquick PCR Purification Kit for sequencing PCR product was from Qiagen. Random hexamer primers, DNase I Amplification Grade, real-time PCR primers for GNRHR-I and 18S were supplied by Invitrogen.
Tritiated progesterone, PGF2α, PGE2, and AA were purchased from Amersham Biosciences Ltd, while non-radioactive hormones and specific antisera were from Sigma–Aldrich. The NOS detect TM assay kit was purchased from Alexis Corp. (Läufelfingen, Germany).
Goat polyclonal anti-GNRHR-I primary antibody (sc-8681) and biotin rabbit anti-goat IgG secondary antibody (sc-2774) were obtained from Santa Cruz Biotechnology (Santa Cruz, CA, USA). Mouse monoclonal anti GNRH-I primary antibody (MAB 5456) was purchased by Chemicon International (Temecula, CA, USA) and biotin goat anti-mouse IgG secondary antibody was obtained by Zymed Laboratories (South San Francisco, CA, USA). The avidin–biotin complex (ABC, Vector Elite Kit) and the chromogen 3,3′-diaminobenzidine tetrachloride (DAB) were purchased from Vector Laboratories (Burlingame, CA, USA). The mouse monoclonal anti-β-tubulin antibody was from Sigma–Aldrich. PageRuler Ladder for western blot (WB) analysis was obtained from Fermentas (Burlington, Ontario, Canada). The enhanced chemiluminescence detection system for WB (Immobilon Western Chemiluminescent HRP Substrate) was purchased from Millipore (Billerica, MA, USA). The HRP-conjugated rabbit anti-goat IgG antibody as well as the Restore Western Blot Stripping Buffer was obtained from Thermo Fisher Scientific (Rockford, IL, USA). The Protran nitrocellulose membranes were from Whatman (Dassel, Germany). Biomax film for chemiluminescent blot imaging was from Kodak Laboratories. The bands were quantified using Quantity One software (Bio-Rad Laboratories).
The kit for the protein assay was purchased from Bio-Rad Laboratories. Medium 199 and Earle's Balanced Salt Solution were from Gibco. HEPES, BSA, buserelin, antide, eCG, 4-bromophenacyl bromide, compound 48/80, propranolol, 2-O-methyladenosine, acetylsalicylic acid (ASA), and N-nitro-l-arginine methyl were purchased from Sigma, 2-amino-ethoxydiphenylborate was from Tocris Cookson Ltd (Northpoint Fourth Way Avonmouth, UK), and 1-hexadecyl-2-acetyl glycerol was from Biomol Research Laboratories (Plymouth Meeting, PA, USA), whereas all other pure grade chemicals and reagents were obtained locally.
Animals and hormonal regimen
Sexually mature New Zealand White female rabbits (3.5–4 kg body weight and 5 months of age) raised in premises owned by the University of Perugia were used for all experiments. The rabbits were housed individually in wire mesh cages under controlled light (14 h light:10 h darkness; lights off at 2100 h) and temperature (18–24 °C) conditions. Each animal had free access to food and water. Pseudopregnancy was induced with 20 IU eCG followed by 0.8 μg buserelin 3 days later (Stradaioli et al. 1997). The day of buserelin injection was designated as day 0. The animals (n=5 for each luteal stage) were killed by cervical dislocation in accordance with the guidelines and principles for the care and use of research animals. The protocols involving the use of the animals for these experiments were approved by the Bioethic Committee of the University of Perugia.
Tissue collection
Upon killing, reproductive tracts were promptly removed and thoroughly washed with saline. Within a few minutes, the CL were excised from ovaries and, after careful dissection of non-luteal tissue by fine forceps under stereoscopic magnification, immediately processed for in vitro experiments or frozen at −80 °C, after rinsing with RNase-free PBS, for later evaluation of gene and protein expression. For the immunohistochemical detection of GNRHR-I and GNRH-I, the ovaries of two additional animals for each luteal stage were fixed by immersion in 4% (w/v) formaldehyde in PBS (pH 7.4) for 24 h at room temperature, and subsequently processed for embedding in paraffin following routine tissue preparation procedures. The pituitary gland of a control rabbit was taken and processed as described above.
Immunohistochemistry
Tissue sections were deparaffinized in xylene, rehydrated through graded ethanol, and finally rinsed in distilled water; then they were submitted to antigen retrieval by microwaving in citrate buffer solution (10 mmol/l, pH 6) at 700 W for 10 min. The slides were cooled at room temperature. After rinsing with Tris-buffered saline (TBS), the specimens were dipped in 3% H2O2 (v/v) in methanol for 1 h to quench the endogenous peroxidase activity and rinsed in TBS. Background labeling was prevented by incubating the sections with normal rabbit serum (for GNRHR-I) and normal goat serum (for GNRH-I), diluted 1:10, for 30 min at room temperature. Subsequently, the sections were incubated for 1 h at 37 °C and then overnight at 4 °C in a moist chamber with goat polyclonal anti-GNRHR-I primary antibody, raised against 18 C-terminal residues, diluted 1:200 in TBS containing 0.1% BSA, whereas the ovary sections were treated also with mouse monoclonal anti-GNRH primary antibody diluted as described above. The next day, the slides were rinsed three times in TBS for 5 min each, treated again with normal serum and then incubated for 30 min at room temperature with biotin rabbit anti-goat (for GNRHR-I) or biotin goat anti-rabbit (for GNRH-I) secondary antibodies diluted 1:200 in TBS. After TBS washes, the slides were exposed to ABC for 30 min and rinsed again with TBS. The peroxidase activity sites were visualized using DAB kit as chromogen; then the specimens were rinsed twice with distilled water for 5 min each, counterstained with hematoxylin, washed in running tap water and, finally, dehydrated passing through graded ethanol (70, 95, and 100% v/v), cleared in xylene and mounted with medium for light microscopy. Tissue sections in which the primary antibody was omitted or substituted by goat or mouse IgG were used as negative controls of non-specific staining. Pituitary tissue sections from a control rabbit were included as positive controls for the GNRHR-I reaction. The intensity of GNRHR-I and GNRH-I immunostaining in follicles and CL was assessed and compared microdensitometrically. The image analysis system (IAAS 2000 image analyzer, Delta Sistemi, Rome, Italy) was calibrated by taking the background developed in sections incubated with non-immune serum as ‘zero’ and a conventional value of maximum intensity of staining (Tomassoni et al. 2010) as ‘100’. Data are expressed in arbitrary units.
GNRHR-I, FSHR, and LHR real-time RT-PCR
Total RNA was extracted from CL of three rabbits for each luteal stage as described previously (Boiti et al. 2005). According to the protocol provided by the manufacturer, 5 μg total RNA were reverse transcribed in 20 μl of iSCRIPT cDNA using random hexamer. Genomic DNA contamination was checked by developing the PCR without reverse transcriptase. Serial experiments were performed to optimize the quantitative reaction, efficiency, and threshold cycle (Ct) values. The optimal 25 μl PCR volume contained 12.5 μl iQ SYBR Green SuperMix, 1 μl forward and reverse primers (stock concentration 10 μmol/l), and 25 μl water. The primers used are listed in Table 1. All reagents were mixed as a master mix and distributed into a 96-well PCR plate before adding 2 μl cDNA (diluted 10-fold with water). For every PCR run were included reaction controls without template, as negative controls, and control without reverse transcriptase in RT, to assess whether the RNA was free by genomic DNA contamination. Samples’ amplification fidelity was also verified by agarose gel electrophoresis. PCR was performed on iCycler iQ (Bio-Rad) with an initial incubation at 95 °C for 1.5 min, followed by 40 cycles at 95 °C for 15 s, 53 °C for 30 s, during which fluorescence data were collected. The Ct value of each trace was automatically computed. PCR products were purified and sequenced by Qiaquick PCR Purification Kit according to the manufacturer's protocol. The 18S Ct housekeeping gene, was determined to normalize sample variations in the amount of starting cDNA. Standard curves generated by plotting the Ct value against the log cDNA standard dilution (1/10 dilution) in nuclease-free water were used to compare the relative amount of target genes. The slope of this graph was used to determine the reaction efficiency. Quantification of mRNA samples was performed by iCycler system software. The melting curve analysis, performed immediately after the PCR end cycle, was used to determine the specificity of each primer set. A melt-curve protocol was performed by repeating 80 °C heating for 10 s, from 55 °C with 0.5C increments, during which fluorescence data were collected. The melting temperatures of the products were determined graphically with fluorescence change rate (d(RFU)/dT) versus temperature.
Primers used for gene quantification by real-time PCR
Gene | Accession number | bp | Primers | |
---|---|---|---|---|
GNRHR-I | AY781779 Oryctolagus cuniculus | Forward | 191 | TGATCCACCTCACAAATGGA |
Reverse | ATGAAGGACCCGTGTCAGAG | |||
FSHR | AY429104.1 Homo sapiens | Forward | 150 | GAGGAATGCCATTGAACTGAGG |
Reverse | GGAAGGTTGGAGAACACATCTG | |||
LHR | S57793 Homo sapiens | Forward | 118 | CTGGAGAAGATGCACAATGG |
Reverse | CAATTAGCCTCTGAATGGACTC | |||
18S | X03205.1 Homo sapiens | Forward | 148 | CGATCAGATACCGTCGTAGT |
Reverse | TTCCTTTAAGTTTCAGCTTTGC |
Western blotting
Total luteal proteins were extracted from a pool of ten CL of three rabbits for each luteal stage. The CL were homogenized in 1 ml ice-cold RIPA buffer (PBS containing 1% Igepal CA-630, 0.5% sodium deoxycholate, 0.1% SDS v/v) as described previously (Zerani et al. 2004). After incubation at 4 °C for 20 min, the homogenates were centrifuged at 12 000 g for 60 min at 4 °C. The supernatants were collected, and their protein concentrations were measured using the protein assay kit with BSA as a standard. Equal amounts of proteins (20 μg) were separated by discontinuous 10% SDS-PAGE (w/v) with 4% stacking gel (w/v) for 40 min at 200 V and 500 mA, after which proteins were transferred onto a nitrocellulose membrane for 1 h at 100 V and 350 mA. After the transfer, non-specific binding of antiserum was prevented by incubation with 0.05% Tween-20 (v/v), 5% non-fat dry milk (w/v), and 5% BSA (w/v) in TBS. Blocked membranes were then probed with antibody against GNRHR (1:500) overnight at 4 °C. After washing with TBS, the membranes were incubated with HRP-labeled rabbit anti-goat IgG secondary antibody (1:20 000) at room temperature for 1 h. All antibody incubations were performed in TBS containing 5% BSA (w/v) and 0.05% Tween-20 (v/v). After washing with TBS, the immunocomplexes were visualized by using an enhanced chemiluminescence detection kit according to the manufacturer's protocol and exposed to X-ray film. Blot images were acquired, and the intensities of the bands were quantified by densitometric analysis using Quantity One software (Bio-Rad Laboratories). Membranes were then reblotted with mouse anti-β-tubulin monoclonal antibody (0.05 μg/ml) and processed as described above to verify equal loading of proteins. Data were normalized by reference to β-tubulin after membrane stripping for 15 min with stripping solution. Values were expressed as arbitrary units of relative abundance of the specific protein normalized with that of β-tubulin used as a loading control.
In vivo experiments
At days 4, 9, and 13 of pseudopregnancy, two groups of rabbits (n=6/group per day) were treated i.m. with either saline (0.2 ml) or buserelin (0.8 μg). From each rabbit, blood samples were collected by venipuncture of the marginal ear vein 0, 24, and 48 h after injection. The samples, collected in EDTA vacutainers, were centrifuged at 3000 g for 15 min, and the frozen plasma was stored until assayed for progesterone concentrations to assess the functional status of the ovarian CL. For the purposes of this work, functional luteolysis was defined as a 50% decrease in plasma progesterone from pre-treatment values, while complete luteolysis as the failure of CL to secrete progesterone so that blood levels fall below 1.0 ng/ml, which are found in estrous rabbits (Browning et al. 1980).
In vitro experiments
A method reported previously was used for the in vitro study (Boiti et al. 2001). Day 4, 9, or 13 CL were randomly distributed (one CL/well) into incubation wells (Becton Dickinson Co., Clifton, NJ, USA) in 1 ml culture medium 199 with Earle's Balanced Salt Solution containing 2.2 mg/ml sodium bicarbonate, 2.3 mg/ml HEPES, and 3% BSA (w/v), referred to here as M199. Before treatment, the CL were quartered inside each well using fine forceps. Each set of incubation wells was divided into 12 experimental groups, each were formed of five wells as follows: I) medium alone as control; II) GNRH-I agonist (buserelin, 200 nmol/l); III) buserelin+GNRH-I antagonist (antide, 100 nmol/l); IV) buserelin+PLA2 inhibitor (4-bromophenacyl bromide, 2 μmol/l); V) buserelin+PLC inhibitor (compound 48/80, 2 μmol/l); VI) buserelin+PLD inhibitor (propranolol, 10 μmol/l); VII) buserelin+AC inhibitor (2-O-methyladenosine, 2 μmol/l); VIII) buserelin+inositol 1,4,5-trisphosphate (IP3) antagonist (2-aminoethyl diphenylborinate, 100 μmol/l); IX) buserelin+diacylglycerol (DAG) antagonist (1-hexadecyl-2-acetyl glycerol, 100 μmol/l); X) buserelin+protein kinase C (PKC) inhibitor (staurosporine 20 μmol/l); XI) buserelin+COX inhibitor (ASA, 100 μmol/l); XII) buserelin+NOS inhibitor (N-nitro-l-arginine methyl ester, 100 μmol/l). The culture plates were incubated at 37 °C in air with 5% CO2. The media of each well were collected after 4 h of incubation and stored immediately at −20 °C for later determination of progesterone, PGF2α, and PGE2. Preliminary evidence led us to choose the incubation conditions and the minimum effective dose for buserelin used in the present in vitro study (Fig. 1).
Hormone determination
Progesterone, PGF2α, and PGE2 were determined following the RIA protocols reported previously (Boiti et al. 2001). Intra- and inter-assay coefficients of variation and minimum detectable dose were progesterone: 5 and 9%, 10 pg; PGF2α: 8 and 12%, 19 pg; PGE2: 7 and 13%, 18 pg.
Enzyme activity determination
Luteal COX-1 and -2 activities were determined by measuring the disappearance of the radiolabeled substrate [3H]AA using a method reported previously (Zerani et al. 2007). Selective COX-2 inhibitor (NS-398 1 μmol/l) and non-selective COX inhibitor (ASA, 1 mmol/l) were used. For each sample, COX-1 activity was determined by calculating the rate of loss of [3H]AA incubated with selective COX-2 inhibitor. Conversely, the COX-2 activity of each corresponding sample was determined by calculating the rate of loss of [3H]AA incubated without selective COX-2 inhibitor and subtracting from this value that of COX-1. The values for COX-1 and -2 were corrected by subtracting the [3H]AA disappearance values due to other enzymatic activities and to non-enzymatic reactions.
NOS activity was determined by monitoring the conversion of [3H]l-arginine to [3H]l-citrulline with a commercial NOS assay kit, according to the experimental protocol described by Boiti et al. (2000).
Statistical analysis
Data on gene and protein expressions, hormone levels, and enzyme activities were examined by ANOVA followed by the Student–Newman–Keuls t-test. Differences were considered significant at P<0.01.
Results
GNRHR-I and GNRH-I immunolocalization
GNRHR-I
Immunopositivity for GNRHR-I was detected in the cytoplasm of luteal cells and endothelial cells (Fig. 2, photo a); additionally, it was detected in the ooplasm of oocytes and cytoplasm of follicular and thecal cells present both in antral and pre-antral follicles (Fig. 2, photos b and c). Immunosignals were also evident in the ooplasm of oocytes and cytoplasm of follicular cells present in primordial follicles (Fig. 2, photo d), ovarian epithelium (Fig. 2, photo d), and endocrine interstitial cells (Fig. 2, photo a). The intervening fibroblast cells within the CL were always unreactive. Differences in the presence and intensity of immunostaining were not observed during the three luteal stages. The control staining procedure failed to disclose appreciable reactivity at any of the sites described for GNRHR (inserts of Fig. 2, photos a, b and d). The cytoplasm of pituitary cells exhibited intense immunostaining (Fig. 2, photo e). In the ovary collected at day 4 of pseudopregnancy, the intensity of GNRHR-I immunostaining was higher in CL than in follicles (39.92±5.46 and 29.53±3.99 respectively, P<0.01), whereas no differences were found in the ovaries at later stages of pseudopregnancy (data not shown).
GNRH-I
Positive immunosignals for GNRH-I were evidenced in the cytoplasm of luteal cells at all stages of pseudopregnancy (Fig. 3, photo a). The cytoplasm of follicular and thecal cells and the ooplasm within antral follicles were immunopositive (Fig. 3, photo b) as well as the follicular cells and the ooplasm within pre-antral and primordial follicles (Fig. 3, photo c). The ovarian epithelium (Fig. 3, photo b), endothelial cells (Fig. 3, photo a), and interstitial cells displayed an intense immunoreactivity. The connective tissue within the CL was not immunostained, whereas that of the ovarian cortex immunoreacted to GNRH-I (Fig. 3, photo c). The control sections were not stained at any sites examined for GNRH-I (insert of Fig. 3, photo b). The intensity of GNRH-I immunostaining did not show any difference between CL and follicles at any luteal stage considered (data not shown).
Expression of mRNA and protein for luteal GNRHR-I, FSHR, and LHR
GNRHR-I transcript was expressed in CL at different stages of pseudopregnancy with no changes in its relative abundance (Fig. 4, panels A and B). The protein abundance of GNRHR-I did not differ between days 4, 9, and 13 of pseudopregnancy (Fig. 4, panels C and D). The expressions of FSHR and LHR transcripts were higher (P<0.01) in day 13 CL than in day 4 and 9 CL (Fig. 5).
In vivo experiments
Buserelin administration at days 9 and 13 of pseudopregnancy decreased (P<0.01) the peripheral plasma progesterone concentrations for 48 and 24 h respectively (Fig. 6), whereas buserelin injection at day 4 did not affect progesterone levels (Fig. 6, panel B).
In vitro experiments
Hormone releases
Basal in vitro progesterone secretion was higher in day 9 CL (P<0.01) than in day 4 and 13 CL, secretion in day 13 CL was higher (P<0.01) than in day 4 CL (Fig. 7, top panel). Buserelin reduced (P<0.01) progesterone in vitro secretion by both day 9 and 13 CL, but had no effect on day 4 CL (Fig. 7, top panel). The co-incubation of buserelin with antide, compound 48/80, 2-aminoethyl diphenylborinate, 1-hexadecyl-2-acetyl glycerol, or staurosporine counteracted (P<0.01) the effects of buserelin alone on day 9 and 13 CL (Fig. 7, top panel). The co-treatments with ASA or N-nitro-l-arginine methyl ester inhibitor increased (P<0.01) progesterone secretion by both day 9 and 13 CL (Fig. 7, top panel). The co-treatment with 4-bromophenacyl bromide, propranolol, or 2-O-methyladenosine did not affect the buserelin-decreased progesterone secretion on day 9 and 13 CL (data not shown).
Basal in vitro PGF2α secretion was higher in day 13 CL (P<0.01) than in day 4 and 9 CL, and secretion in day 9 CL was higher (P<0.01) than in day 4 CL (Fig. 7, bottom panel). Buserelin increased (P<0.01) PGF2α in vitro secretion by both day 9 and 13 CL (Fig. 7, bottom panel). The co-incubation of buserelin with antide, compound 48/80, 2-aminoethyl diphenylborinate, 1-hexadecyl-2-acetyl glycerol, or staurosporine counteracted (P<0.01) the effects of buserelin alone on day 9 and 13 CL (Fig. 7, bottom panel). The co-treatment with ASA decreased (P<0.01) PGF2α secretion in all CL types, while the co-treatment with N-nitro-l-arginine methyl ester did not affect the buserelin-induced PGF2α secretion (Fig. 7, bottom panel).
Basal in vitro PGE2 secretion was higher in day 4 CL (P<0.01) than in day 9 and 13 CL (Fig. 7, middle panel). Buserelin reduced (P<0.01) PGE2 in vitro secretion by day 13 CL (Fig. 7, middle panel). The co-incubation of buserelin with antide, compound 48/80, 2-aminoethyl diphenylborinate, 1-hexadecyl-2-acetyl glycerol, or staurosporine counteracted (P<0.01) the effects of buserelin alone on day 13 CL (Fig. 7, middle panel). The co-treatment with ASA decreased (P<0.01) PGE2 secretion in all CL types, while the co-treatment with N-nitro-l-arginine methyl ester did not affect the buserelin-induced PGE2 secretion (Fig. 7, middle panel).
Enzyme activities
Basal COX-1 activity was higher in day 13 CL (P<0.01) than in day 4 and 9 CL (Fig. 8, top panel). All treatments did not affect COX-1 activity in all CL types, except for ASA that decreased (P<0.01) the enzymatic activity (Fig. 8, top panel).
Basal COX-2 activity of day 13 CL was higher (P<0.01) than that of day 4 and 9 CL; COX-2 activity of day 9 CL was higher (P<0.01) than that of day 4 CL (Fig. 8, middle panel). Buserelin increased (P<0.01) COX-2 activity in CL of all stages (Fig. 8, middle panel). The co-incubation of buserelin with antide, compound 48/80, 2-aminoethyl diphenylborinate, 1-hexadecyl-2-acetyl glycerol, or staurosporine counteracted (P<0.01) the effects of buserelin alone (Fig. 8, middle panel). The co-treatments with ASA decreased (P<0.01) COX-2 activity in all CL types, while the co-treatments with N-nitro-l-arginine methyl ester did not affect the buserelin-induced COX-2 activity increase (Fig. 8, middle panel).
Basal NOS activity of day 4 CL was higher (P<0.01) than that of day 9 CL (Fig. 8, bottom panel). Buserelin increased (P<0.01) NOS activity in both day 9 and 13 CL (Fig. 8, bottom panel). The co-incubation of buserelin with antide, compound 48/80, 2-aminoethyl diphenylborinate, 1-hexadecyl-2-acetyl glycerol, or staurosporine counteracted (P<0.01) the effects of buserelin alone at day 9 and 13 (Fig. 8, bottom panel). The co-treatments with N-nitro-l-arginine methyl ester decreased (P<0.01) NOS activity in CL of all stages, while that with ASA in both day 9 and 13 CL (Fig. 8, bottom panel).
All hormone release and enzyme activity data are summarized in the synoptic Table 2.
In vitro effects of GNRH-I agonist (buserelin, GNRH-I) alone or co-incubated with antagonists for GNRH-I (antide, GNRH-Ia), IP3 (2-aminoethyl diphenylborinate, IP3a), and DAG (1-hexadecyl-2-acetyl glycerol, DAGa), or inhibitors for PLC (compound 48/80, PLCi), PKC (staurosporine, PKCi), COX (acetylsalicylic acid, COXi), and NOS (N-nitro-l-arginine methyl ester, NOSi) on progesterone, PGE2, and PGF2α release and COX-1, COX-2, and NOS activities by rabbit CL at days 4, 9, and 13 of pseudopregnancy. Values are means±s.d. of five replicates
Progesterone | PGF2α | PGE2 | COX-1 | COX-2 | NOS | |
---|---|---|---|---|---|---|
Day 4 | ||||||
Control | 1997±273aα | 707±225aα | 796±200aα | 7859±1135aα | 3145±825aα | 19 202±1702aα |
GNRH-I | 1755±152a | 640±167a | 860±105a | 6728±1178a | 9705±2101b | 18 910±1613a |
GNRH-I+GNRH-Ia | 1671±336a | 702±172a | 899±358a | 6264±2667a | 3478±492a | 17 096±4393a |
GNRH-I+PLCi | 2099±248a | 597±112a | 699±158a | 7030±1912a | 3044±508a | 16 475±1788a |
GNRH-I+IP3a | 2280±332a | 715±183a | 899±140a | 8152±1881a | 2587±437a | 20 172±3336a |
GNRH-I+DAGa | 1860±389a | 627±114a | 785±187a | 7961±803a | 3219±567a | 17 438±2315a |
GNRH-I+PKCi | 2001±305a | 617±152a | 725±247a | 7526±1640a | 2746±305a | 18 762±3655a |
GNRH-I+COXi | 1696±371a | 133±47b | 137±47b | 537±323b | 372±252c | 18 472±2255a |
GNRH-I+NOSi | 1947±317a | 637±185a | 899±184a | 6509±2542a | 9497±1072b | 1246±340b |
Day 9 | ||||||
Control | 5197±383aβ | 2089±297aβ | 161±71aβ | 7220±1855aα | 8116±1830aβ | 14 212±2764aβ |
GNRH-I | 2048±240b | 3825±346b | 103±44a | 7442±1505a | 41 452±5559b | 35 421±4201b |
GNRH-I+GNRH-Ia | 5479±1693a | 2170±447a | 179±50a | 8013±962a | 6472±2044a | 15 608±3236b |
GNRH-I+PLCi | 4851±394a | 1842±160a | 174±32a | 8135±1197a | 8135±1094a | 14 197±2369b |
GNRH-I+IP3a | 4458±452a | 2051±536a | 181±19a | 7061±1073a | 7213±1241a | 19 971±5048b |
GNRH-I+DAGa | 4708±250a | 1825±431a | 191±27a | 8048±951a | 6956±1586a | 15 230±3553b |
GNRH-I+PKCi | 4580±664a | 2107±133a | 162±39a | 6053±1139a | 8237±1039a | 15 246±4083b |
GNRH-I+COXi | 7604±633c | 202±27c | 78±10b | 767±197b | 318±144c | 7065±1667c |
GNRH-I+NOSi | 7919±558c | 4343±521b | 205±52a | 6436±3051a | 43 295±5218b | 1158±304d |
Day 13 | ||||||
Control | 3789±182aγ | 3304±245aγ | 248±46aβ | 16 440±5308aβ | 26 599±4624aγ | 17 549±950aαβ |
GNRH-I | 1894±176b | 4614±266b | 113±37b | 16 636±5068a | 45 060±3325b | 36 432±3697b |
GNRH-I+GNRH-Ia | 3680±877a | 2768±485a | 255±97a | 15 640±2021a | 24 809±3494a | 17 212±4669a |
GNRH-I+PLCi | 3922±356a | 2819±492a | 272±47a | 18 757±4722a | 23 776±4026a | 17 616±1326a |
GNRH-I+IP3a | 4165±434a | 2816±296a | 247±41a | 18 126±3678a | 27 674±5558a | 20 435±2380a |
GNRH-I+DAGa | 3850±391a | 3104±262a | 268±31a | 19 989±2879a | 26 756±4166a | 18 426±1684a |
GNRH-I+PKCi | 3892±576a | 2981±369a | 282±60a | 20 608±3325a | 31 061±2442a | 18 326±5029a |
GNRH-I+COXi | 6108±613c | 185±34c | 86±31b | 796±482b | 292±71c | 5640±1982c |
GNRH-I+NOSi | 6161±693c | 4922±366b | 261±70a | 15 182±3319a | 40 792±4471b | 963±164d |
Different letters indicate significantly different values (P<0.01): Greek letters among control CL days; Latin letters among experimental groups of the same CL day.
Discussion
The present study, for the first time, indicates the presence of a GNRH/GNRHR system that modulates the steroidogenic luteal function of pseudopregnant rabbits.
The first evidence for the presence of ovarian GNRHR, through ligand-specific binding sites, was observed in rat granulosa and luteal cells (Reeves et al. 1980, Pieper et al. 1981, Olofsson et al. 1995). Subsequently, the presence of GNRHR mRNA was also identified in human granulosa luteal cells (Kang et al. 2000) and in both follicles and CL of bovine ovary (Ramakrishnappa et al. 2001). Among the different stages of CL, GNRHR transcripts were clearly detectable in stage III luteal tissues, whereas the expression levels were slightly lower in stage II, and barely or no detectable signal was evident from stages I and IV luteal tissues (Milvae et al. 1984).
Our study showed the presence of GNRHR-I at both transcript and protein levels in CL of pseudopregnant rabbits. Apparently, both mRNA and protein abundances did not change at different luteal stages, suggesting that in the rabbit this receptor is not dynamically regulated during the CL lifespan. In both early- and mid-luteal stages, luteal GNRHR-I, FSHR, and LHR transcript levels were expressed at approximately the same relative abundance. By converse, the GNRHR-I to FSHR or to LHR mRNAs ratio decreased in CL of late-luteal stage as a consequence of the two- to three-fold up-regulation of these two ovarian markers. Positive staining for GNRHR-I was clearly evidenced in the cytoplasm of luteal cells during all stages of pseudopregnancy.
The presence of GNRH-I has been revealed in the ovary of rats (Aten et al. 1986, Schirman-Hildelsheim et al. 2005, Sengupta et al. 2008), monkeys (Chakrabarti et al. 2008), and humans (Aten et al. 1987). The present investigation confirms that also the ovary of rabbits expresses GNRH-I in the cytoplasm of different cell types, including the luteal cells, independently of the luteal stage examined.
The presence of GNRHR-I and GNRH-I, similarly to what found in other species, supports a direct role of GNRH in the regulation of luteal functions in rabbits. GNRHR-I immunoreactivity was also observed in other rabbit ovary cells as oocytes, follicular and thecal cells, so stressing the potential paracrine and/or autocrine role of the GNRH/GNRHR system in modulating ovary function in mammals (Ramakrishnappa et al. 2005). Interestingly, in the early-luteal stage, the microdensitometry results revealed the presence of larger GNRHR-I amounts in luteal than in follicle cells suggesting that GNRH-I is involved also in the luteotropic processes.
GNRH has been shown to exert a mixture of both inhibitory and stimulatory effects in the gonads, causing either inhibitory or stimulatory effects on ovarian cellular steroid output (Ramakrishnappa et al. 2005). In bovine, two GNRH antagonists, buserelin and deslorelin, showed in vivo and in vitro inhibitory effects on progesterone release (Milvae et al. 1984, D'Occhio & Aspden 1999). Conversely, GNRH stimulatory in vitro effects on progesterone production were demonstrated in monkey and human granulosa cells (Casper et al. 1984, Liu et al. 1991, Bussenot et al. 1993). Finally, other authors documented that GNRH has no effect on progesterone production in human granulosa–lutein cells (Casper et al. 1984). In rabbits, GNRH agonist directly affected the ovary by inducing ovulation in both hypophysectomized animals and in vitro perfused system (Koos & LeMaire 1985). At day 11 of pseudopregnancy, GNRH injection induced luteinization of ovarian follicles and functional regression of CL as assessed by decline in serum progesterone levels within 3 days (Rippel et al. 1976). Depending on the dosage administered between days 8 and 13 of gestation, GNRH was found to either cause abortion by luteolysis or support fetal survival by stimulating luteotropic activity (Eisenberg et al. 1984). Conversely, other authors found that GNRH did not exhibit any direct inhibitory effect on steroid production in extracorporeal perfused ovary (Eisenberg et al. 1984).
Our data strengthen the idea that GNRH-I has a direct role in the down-regulation of rabbit CL progesterone production. The in vitro experiments showed that buserelin reduced progesterone production by cultured CL collected at days 9 and 13 of pseudopregnancy, whereas this steroid synthesis was not affected by GNRH-I agonist in day 4 CL. In addition, this GNRH-I role is also supported by the data of the co-incubation with antide, a GNRH-I antagonist, that counteracted the buserelin-induced progesterone decrease. The in vivo experiments mirrored the above-reported in vitro data; in fact, blood progesterone levels were down-regulated when buserelin was administered at mid- and late-pseudopregnancy. Even if we cannot exclude that this in vivo buserelin-induced circulating progesterone decrease is also caused by a pituitary gonadotropin down-regulation, the in vivo and in vitro coincident lack of buserelin effects during early-pseudopregnancy suggests that this progesterone decline is primarily or exclusively due to the direct action of GNRH-I agonist on pseudopregnant CL.
Recently, we demonstrated, through in vitro studies on CL of pseudopregnant rabbits, that PGF2α and PGE2 affect progesterone release differently, depending on the luteal stage (Boiti et al. 2001). PGE2 depressed NOS activity and increased progesterone production in day 4 CL but was totally ineffective in day 9 CL; on the contrary, PGF2α up-regulated NOS activity and induced functional luteolysis in day 9 CL but had no effect on CL collected at the early-luteal phase. This physiological mechanism protects the growing CL from functional luteolysis to occur in the early-luteal stage until day 6 of pseudopregnancy, when CL shift from refractoriness to partial responsiveness to exogenous PGF2α, and acquire luteolytic competence (Boiti et al. 1998, 2001, 2003). The data on the buserelin ineffectiveness at early-luteal stage suggest that this luteolytic-protecting mechanism includes also the GNRH-I/GNRHR-I system, at least in vitro.
Many studies indicate that GNRH, complexing with cognate receptor, mainly activates PLC via Gq/11 family G proteins. PLC catalyzes the hydrolysis of phosphatidylinositol 4,5-bisphosphate to IP3 and DAG (Millar 2005). Our data confirmed that this GNRH post-receptorial mechanism is present also in the rabbit CL. In fact, the effects of buserelin in decreasing progesterone release on day 9 and 13 CL were counteracted by the co-incubation with PLC and PKC inhibitors or IP3 and DAG antagonists, while the inhibitors of AC, PLA2, and PLD had no influences.
A recent study showed that in the gonadotrope cell line LβT2, GNRH produced a marked time-dependent induction of COX-2 expression and an increase in PGE2 and PGF2α production (Naor et al. 2007). These authors hypothesize a novel GNRH signaling pathway mediated by PGF2α, which acts through an autocrine/paracrine modality. The present work suggests that the GNRH-I effects on rabbit CL progesterone production are mediated by the activation of the COX-2/PGF2α/NOS pathway. In fact, buserelin increased in vitro COX-2 activity, PGF2α production, and NOS activity in both day 9 and 13 CL, coinciding with decreasing progesterone release. In addition, the co-treatment with antide, PLC and PKC inhibitors, or IP3 and DAG antagonists counteracted the buserelin actions on COX-2, PGF2α, NOS, and progesterone, while COX inhibitor counteracted those on NOS and progesterone, and, in turn, NOS inhibitor counteracted those on progesterone only. The decrease in PGE2 induced by buserelin in day 13 CL is probably due to the increased conversion of this PG into PGF2α through PGE2-9-ketoreductase enzyme activation and/or the inactivation of PGE2 synthase. A similar intracellular mechanism was reported for the PGF2α auto-amplifying loop utilized by the rabbit CL during the luteolytic process (Zerani et al. 2007). Our results also indicated that early CL refractoriness to luteolytic factors extends to GNRH-I as well, given that the all in vitro effects of buserelin on hormone productions and enzymatic activities occur only in day 9 and 13 CL.
In conclusion, this study provides the evidence that GNRH-I, in autocrine, paracrine, and/or endocrine manner, directly down-regulates progesterone production in rabbit CL that have acquired luteolytic competence, with a receptorial/post-receptorial mechanism as summarized in Fig. 9. GNRH-I couples its cognate type I receptor activating PLC and, in turn, releasing IP3 and DAG; these two intracellular messengers, through the PKC activation, start up COX-2 activity, with the consequent rise of PGF2α production; this PG induces (via paracrine, autocrine, and/or intracrine mechanism) an increase in NOS activity (Boiti et al. 2003); finally NO down-regulates progesterone levels (Boiti et al. 2000).
However, though the present data throw new light on the knowledge of the physiological mechanisms regulating luteal activity, further studies are needed to better understand the GNRH/GNRHR system fine-tuning that controls rabbit CL life span, including the presence of type II GNRHR or the involvement of PGF2α synthase and PGE2-9-ketoreductase, key enzymes of the PGF2α synthesis.
Declaration of interest
The authors declare that there is no conflict of interest that could be perceived as prejudicing the impartiality of the research reported.
Funding
This work was supported by a grant from the Italian Ministero dell'Università e della Ricerca (Contract n. 2005078028) and by donation of the BC Red-water Trust.
Acknowledgements
The authors gratefully acknowledge the revision of the English text by Sheila Beatty. The authors wish to thank Mrs G Mancini for her assistance in tissue collection.
References
Aten RF, Williams AT, Behrman HR & Wolin DL 1986 Ovarian gonadotropin-releasing hormone-like protein(s): demonstration and characterization. Endocrinology 118 961–967 doi:10.1210/endo-118-3-961.
Aten RF, Polan ML, Bayless R & Behrman HR 1987 A gonadotropin releasing hormone (GnRH)-like protein in human ovaries: similarity to the GnRH-like ovarian protein of the rat. Journal of Clinical Endocrinology and Metabolism 64 1288–1293 doi:10.1210/jcem-64-6-1288.
Boiti C, Canali C, Zerani M & Gobbetti A 1998 Changes in refractoriness of rabbit corpora lutea to a prostaglandin F-2 alpha analogue, alfaprostol, during pseudopregnancy. Prostaglandins and Other Lipid Mediators 56 255–264 doi:10.1016/S0090-6980(98)00053-7.
Boiti C, Zerani M, Zampini D & Gobbetti A 2000 Nitric oxide synthase activity and progesterone release by isolated corpora lutea of rabbits in early- and mid-luteal phase of pseudopregnancy are differently modulated by prostaglandin E-2 and prostaglandin F-2α via adenylate cyclase and phospholipase C. Journal of Endocrinology 164 179–186 doi:10.1677/joe.0.1640179.
Boiti C, Zampini D, Zerani M, Guelfi G & Gobbetti A 2001 Prostaglandin receptors and role of G protein-activated pathways on corpora lutea of pseudopregnant rabbit in vitro. Journal of Endocrinology 168 141–151 doi:10.1677/joe.0.1680141.
Boiti C, Guelfi G, Zampini D, Brecchia G, Gobbetti A & Zerani M 2003 Regulation of nitric-oxide synthase isoforms and role of nitric oxide during prostaglandin F2alpha-induced luteolysis in rabbits. Reproduction 125 807–816 doi:10.1530/rep.0.1250807.
Boiti C, Guelfi G, Brecchia G, Dall'Aglio C, Ceccarelli P, Maranesi M, Mariottini C, Zampini D, Gobbetti A & Zerani M 2005 Role of endothelin-1 system in the luteolytic process of pseudopregnant rabbits. Endocrinology 146 1293–1300 doi:10.1210/en.2004-1099.
Browning JY, Keyes PL & Wolf RC 1980 Comparison of serum progesterone, 20 alpha-dihydroprogesterone, and estradiol-17 beta in pregnant and pseudopregnant rabbits: evidence for postimplantation recognition of pregnancy. Biology of Reproduction 23 1014–1019 doi:10.1095/biolreprod23.5.1014.
Bussenot I, Azoulay-Barjonet C & Parinaud J 1993 Modulation of the steroidogenesis of cultured human granulosa–lutein cells by gonadotropin-releasing hormone analogs. Journal of Clinical Endocrinology and Metabolism 76 1376–1379 doi:10.1210/jc.76.5.1376.
Casper RF, Erickson GF & Yen SS 1984 Studies of the effect of GnRH and its agonist on human luteal steroidogenesis in vitro. Fertility and Sterility 42 39–43.
Chakrabarti N, Subbarao T, Sengupta A, Xu F, Stouffer RI & Sridaran R 2008 Expression of mRNA and proteins for GnRH I and II and their receptors in primate corpus luteum during menstrual cycle. Molecular Reproduction and Development 75 1567–1577 doi:10.1002/mrd.20898.
Conn PM & Crowley WF Jr 1994 Gonadotropin-releasing hormone and its analogues. Annual Review of Medicine 45 391–405 doi:10.1146/annurev.med.45.1.391.
D'Occhio MJ & Aspden WJ 1999 Endocrine and reproductive responses of male and female cattle to agonists of gonadotropin-releasing hormone. Journal of Reproduction and Fertility. Supplement 54 101–114.
Dubois EA, Zandbergen MA, Peute J & Goos HJ 2002 Evolutionary development of three gonadotropin-releasing hormone (GnRH) systems in vertebrates. Brain Research Bulletin 57 413–418 doi:10.1016/S0361-9230(01)00676-1.
Eisenberg E, Kitai H, Kobayashi Y, Santulli R & Wallach EE 1984 Gonadotropin-releasing hormone: effects on the in vitro perfused rabbit ovary. Biology of Reproduction 30 1216–1221 doi:10.1095/biolreprod30.5.1216.
Gorbman A & Sower SA 2003 Evolution of the role of GnRH in animal (Metazoan) biology. General and Comparative Endocrinology 134 207–213 doi:10.1016/j.ygcen.2003.09.018.
Hilliard J, Pang CN & Sawyer CH 1976 Effects of luteinizing hormone-releasing hormone on fetal survival in pregnant rabbits. Fertility and Sterility 27 421–425.
Kang S, Cheng K, Nathwani P, Choi K & Leung PC 2000 Autocrine role of gonadotropin-releasing hormone and its receptor in ovarian cancer cell growth. Endocrinology 13 297–304 doi:10.1385/ENDO:13:3:297.
Koos RD & LeMaire WJ 1985 The effects of a gonadotropin-releasing hormone agonist on ovulation and steroidogenesis during perfusion of rabbit and rat ovaries in vitro. Endocrinology 116 628–632 doi:10.1210/endo-116-2-628.
Liu Y, Hu ZY, Feng Q & Zou RJ 1991 Paradoxical effect of a GnRH agonist on steroidogenesis in cultured monkey granulosa cells. Science in China. Series B 34 1452–1460.
Millar RP 2005 GnRHs and GnRH receptors. Brain Research Bulletin 88 5–28 doi:10.1016/j.anireprosci.2005.05.032.
Milvae RA, Murphy BD & Hansel W 1984 Prolongation of the bovine estrous cycle with a gonadotropin-releasing hormone analog. Biology of Reproduction 31 664–670 doi:10.1095/biolreprod31.4.664.
Naor Z 2009 Signaling by G-protein-coupled receptor (GPCR): studies on the GnRH receptor. Frontiers in Neuroendocrinology 30 10–29 doi:10.1016/j.yfrne.2008.07.001.
Naor Z, Jabbour HN, Naidich M, Pawson AJ, Morgan K, Battersby S, Millar MR, Brown P & Millar RP 2007 Reciprocal cross talk between gonadotropin-releasing hormone (GnRH) and prostaglandin receptors regulates GnRH receptor expression and differential gonadotropin secretion. Molecular Endocrinology 21 524–537 doi:10.1210/me.2006-0253.
Olofsson JI, Conti CC & Leung PCK 1995 Homologous and heterologous regulation of gonadotropin-releasing hormone receptor gene expression in preovulatory rat granulosa cells. Endocrinology 136 974–980 doi:10.1210/en.136.3.974.
Pieper DR, Richards JS & Marshall JC 1981 Ovarian gonadotropin-releasing hormone (GnRH) receptors: characterization, distribution, and induction by GnRH. Endocrinology 108 1148–1155 doi:10.1210/endo-108-4-1148.
Ramakrishnappa N, Merwe GKVD & Rajamahendran R Gonadotropin-releasing hormone receptor messenger ribonucleic acid expression in bovine ovary Biology of Reproduction 64 Supplement 2001 229.
Ramakrishnappa N, Rajamahendran R, Lin YM & Leung PCK 2005 GnRH in non-hypothalamic reproductive tissues. Animal Reproduction Science 88 95–113 doi:10.1016/j.anireprosci.2005.05.009.
Reeves JJ, Seguin C, Lefebvre FA, Kelly PA & Labrie F 1980 Similar luteinizing hormone-releasing hormone binding sites in rat anterior pituitary and ovary. PNAS 77 5567–5571 doi:10.1073/pnas.77.9.5567.
Rippel RH, Johnson ES & Kimura ET 1976 Regression of corpora lutea in the rabbit after injection of a gonadotropin-releasing peptide. Proceedings of the Society for Experimental Biology and Medicine 152 29–32.
Schirman-Hildelsheim TD, Bar T, Ben-Aroya N & Koch Y 2005 Differential gonadotropin releasing hormone (GnRH) and GnRH receptor messenger ribonucleic acid expression pattern in different tissues of the female rat across the estrous cycle. Endocrinology 146 3401–3408 doi:10.1210/en.2005-0240.
Sengupta A, Chakrabarti N & Sridaran R 2008 Presence of immunoreactive gonadotropin releasing hormone (GnRH) and its receptor (GnRHR) in rat ovary during pregnancy. Molecular Reproduction and Development 75 1031–1044 doi:10.1002/mrd.20834.
Somoza GM, Miranda LA, Strobl-Mazzulla P & Guilgur LG 2002 Gonadotropin-releasing hormone (GnRH): from fish to mammalian brains. Cellular and Molecular Neurobiology 22 326–336 doi:10.1023/A:1021888420271.
Stradaioli G, Verini-Supplizi A, Monaci M, Canali C & Boiti C 1997 Effects of different doses of PMSG on ovarian response and in vitro embryo development in rabbit. World Rabbit Science 3 143–148.
Tomassoni D, Bramanti V & Amenta F 2010 Expression of aquaporins 1 and 4 in the brain of spontaneously hypertensive rats. Brain Research 1325 155–163 doi:10.1016/j.brainres.2010.02.023.
Zerani M, Boiti C, Zampini D, Brecchia G, Dall'Aglio C, Ceccarelli P & Gobbetti A 2004 Ob receptor in rabbit ovary and leptin in vitro regulation of corpora lutea. Journal of Endocrinology 183 279–288 doi:10.1677/joe.1.05507.
Zerani M, Dall'Aglio C, Maranesi M, Gobbetti A, Brecchia G, Mercati F & Boiti C 2007 Intraluteal regulation of prostaglandin F2α-induced prostaglandin biosynthesis in pseudopregnant rabbits. Reproduction 133 1005–1016 doi:10.1530/REP-06-0107.