Abstract
The present study investigated expression and protein localization of FOXL2 messenger RNA (mRNA) in endometrium of healthy women and in patients with endometriosis during endometrial cycle. In endometriotic lesions, FOXL2 mRNA and protein were evaluated and a possible correlation with activin A mRNA expression changes was also studied. Endometrium was collected from healthy women (n = 52) and from women with endometriosis (n = 31) by hysteroscopy; endometriotic tissues were collected by laparoscopy (n = 38). FOXL2 gene expression analysis in endometrium of healthy women showed a significant expression and no significant changes in mRNA levels between proliferative and secretory phases; a similar pattern was observed in endometrium of patients with endometriosis. Immunohistochemical evaluation showed that FOXL2 protein localized in stromal and glandular cells and colocalized with SUMO-1. FOXL2 mRNA expression was 3-fold higher in endometriosis than in healthy endometrium (P < .01) and a positive correlation between FOXL2 and activin A mRNA was found (P < .05) in endometriosis. In conclusion, FOXL2 mRNA expression and its protein localization do not change during endometrial cycle in eutopic endometrium from healthy individuals or patients with endometriosis; the hyperexpression of FOXL2 in endometriotic lesions suggests an involvement of this transcriptional regulator, probably associated with activin A expression and related to the pathogenesis of endometriosis.
Keywords
Introduction
Transcription factors of the forkhead box (FOX) family are classified as helix-turn-helix DNA-binding proteins and they are components of a variety of signaling pathways involved in reproduction, downstream effector molecules capable of integrating actions of hormones, growth factors, and cytokines. 1 Previous studies suggested a possible role of FOX factors in human endometrial function. Differentiation of human endometrial stromal cells into decidual cells is associated with induction of FOXO1, engaging a transcriptional cross-talk with progesterone receptors (PRs) to coordinate cell cycle regulation and differentiation of endometrial cells. 2,3 Another member of FOX family, FOXP1, is expressed by human endometrium during menstrual cycle, reaching highest levels in mid- and late-secretory phases, and in endometriotic lesions. 4
FOXL2 is a member of the FOX family highly involved in reproductive functions, associated with premature ovarian failure (POF), and infertility. 5 FOXL2 exerts pleiotropic biological effects that include changes in subcellular localization, enhanced protein stability, and transcriptional repression of target genes by promoting the interaction of transcription factors with corepressors. 6 At pituitary level, FOXL2 modulates gonadotropin-releasing hormone (GnRH) receptor gene expression 7 and is involved in regulating activin-mediated follistatin and follicular-stimulating hormone β subunit (FSHβ) gene transcription in gonadotrope cells. 8,9 At ovarian level, FOXL2 is essential for a correct differentiation and steroidogenic activity of granulosa cells, 10,11 thus participating in the development of POF when altered. 12,13 FOXL2 activity is modulated by posttranslational modifications, such as SUMOylation, influencing the endometrial response to progesterone 14,15 and the development of ovarian follicles by repressing the steroidogenic acute regulatory gene in granulosa cells. 16
Our previous studies showed that activin is expressed in human endometrium 17 and is involved in the modulation of cell proliferation, differentiation, apoptosis, and tissue remodeling. 18 -20 Endometriotic cells also express activin A mRNA 21,22 and an effect on inflammation, 23,24 on the invasiveness of endometrial cells in cultured peritoneum, 25 and on neoangiogenesis 24 has been shown.
Following the recent evidences that FOXL2 is also expressed in bovine endometrium, 26 and that a FOXL2 mutation is present in endometrium of women with endometriosis 27,28 and with endometrial carcinoma, 29 the present study investigated (1) FOXL2 mRNA expression and protein localization in endometrium of healthy women and patients with endometriosis during endometrial cycle, also in association with SUMO; (2) in ectopic endometrial tissue from endometriotic lesions; and (3) the possible correlation between FOXL2 and activin A expression in the same specimens.
Material and Methods
Patients and Specimens
The patients (age range 25-42 years) who participated in this study were classified as (1) healthy women (n = 52), with regular menstrual cycle, undergoing laparoscopy for tubal sterilization or due to a benign, nonendometriotic ovarian cyst and (2) women with ovarian endometriosis (n = 38) who underwent laparoscopic surgery. The study was approved by the local Human Investigation Committee and informed consent was obtained from all patients before inclusion. A complete medical history was obtained and physical examination was performed for each patient. Exclusion criteria for the study included prior or current infections, endocrine disorders, or the use of hormonal treatment and GnRH analogs within the past 3 months. Patients included in the study were classified at stage III or IV endometriosis according to the American Society for Reproductive Medicine classification of endometriosis. 30
Endometrial specimens were classified as proliferative or secretory according to the last menstrual period and confirmed by transvaginal ultrasound scans by a transvaginal probe at 4.5 to 7.0 MHz and by the histological criteria of Noyes et al. 31
Endometrial specimens were obtained by hysteroscopy while specimens of endometriotic tissue were collected during the laparoscopic approach. A fragment was immediately submerged in liquid nitrogen to allow subsequent RNA extraction and real-time polymerase chain reaction (RT-PCR) and another was fixed in paraformaldehyde for immunohistochemical analysis. In case of endometrioma, endometrial tissue samples were carefully stripped from the lining cyst wall avoiding taking any of the normal ovarian cortex, as confirmed by histological evaluation.
RNA Extraction and Complementary DNA Preparation
Total RNA was extracted from frozen tissue samples as follows. Briefly, frozen samples were disrupted and homogenized. Then, liquid nitrogen and ground tissue were transferred to an appropriately sized sterile tube, allowing the liquid nitrogen to evaporate and then were immediately processed for RNA extraction using the RNeasy Protect Mini kit according to the manufacturer’s instructions (Qiagen, Hilden, Germany).
The RQ1 RNase-free DNase (Promega, Milan, Italy) was used to degrade single-stranded or double-stranded genomic DNA. Absence of residual genomic DNA was verified by polymerase chain reaction (PCR) on total RNA without reverse transcription (RT). Genomic human DNA was used as a positive control for PCR reactions. RNA quantity was assessed using an ND-1000 Nanodrop Spectrometer (Thermo Fisher Scientific, Wilmington, Delaware). RNA integrity was checked prior to downstream analysis with the FlashGel System (Lonza Group, Ltd, Switzerland) to verify total RNA quality and check for RNA degradation. Complementary DNA (cDNA) was generated from 1 μg of each RNA sample using the High Capacity cDNA Reverse Transcription kit (Applied Biosystems, Foster City, California).
Polymerase Chain Reaction and Sequencing for FOXL2
Polymerase chain reaction was performed by adding 2 µL of each cDNA sample to a 20 µL reaction containing 1× Pfx amplification buffer, 2× enhancer solution, 1 mmol/L MgSO4, 200 μmol/L deoxyribonucleotides, 300 nmol/L forward (5′-CGG CTT TGT CAT GAT GGC CAG-3′) and reverse (5′-GGC CAC GAG TTG TTG AGG AAG-3P′) primers suitable for amplifying the human FOXL2 gene, and 1 U of Platinum Pfx DNA polymerase (Invitrogen, Life Technologies Italia, Monza, Italy). Polymerase chain reaction cycling conditions were as follows: 5 minutes at 98°C; 35 cycles of 30 seconds at 98°C, 30 seconds at 60°C, and 1 minute at 72°C; 10 minutes at 72°C. Nuclease-free water was used as negative control. Polymerase chain reaction products were run on a 2% agarose gel for confirmation of predicted fragment size (624 bp) and, once purified, were labeled using dye terminator cycle sequencing (DTCS) kit (Beckman-Coulter, Milan, Italy) according to the manufacturer’s instructions. DNA sequence was determined using CEQ 8000 automated DNA analyzer (Beckman-Coulter, Milan, Italy).
Quantitative Real Time-PCR
Quantification of FOXL2 and activin A (INHBA) mRNA expression was measured in quadruplicate on 1/10 dilutions of each cDNA in 20 µL reactions using the 1× MESA GREEN qPCR MasterMix Plus for SYBR Assay No ROX (Eurogentec, Seraing, Belgium) on Rotor Gene 6000 (Corbett) with appropriate 300 nmol/L forward and reverse primers.
Nucleotide sequences of FOXL2 (NM_023067.3), INHBA (NM_002192.2), and the housekeeping hypoxanthine phosphoribosyltransferase 1 (HPRT1, NM_000194.2) genes were downloaded from GenBank library files. Polymerase chain reaction primers for real-time assays were designed using Beacon Designer 2.06 (Premier Biosoft International) or Primer3 software, giving special attention to primer length, annealing temperature, base composition, and 3′-end stability. FOXL2 primers were forward 5′-TCC AAT AAA AGA TGG TTT CGT CT-3′ and reverse 5′-TCT TCA GAT AGG GAG AGG GTG A-3′; INHBA primers were forward 5′-GAG ATC ATC ACG TTT GCC GAG-3′ and reverse 5′-GAA GAG CCA GAC TTC TGC ACG-3′; HPRT1 primers were forward 5′-CGT GAT TAG TGA TGA TGA ACC AG-3′ and reverse 5′-CGA GCA AGA CGT TGA GTC CT-3′. All the RT-PCR primer pairs used in this study span exon–exon junctions or are located on different exons, with the exceptions of those for the FOXL2 gene, which is encoded by a single exon.
During preliminary quantitative RT-PCR (qRT-PCR) assays, the optimal primer concentration was determined for each primer pair (300 nmol/L final). For each pair of primers, the efficiency of qRT-PCR (e%) and the correlation coefficients (R 2 ) were determined using serial 1:4 dilutions of template cDNA.
Products were subsequently run on 2% agarose gel to check for size specificity and eventually sequenced. Rotor gene runs consisted of a hot start at 95°C for 5 minutes, followed by 45 cycles of a 2-step protocol: 10 seconds at 95°C and 30 seconds at 59°C, with fluorescence detection at the end of each extension step. Melt curve analysis was also performed to confirm specific products and standard curves were performed to confirm efficient amplification of each gene before final analysis of all samples. HPRT1 should be considered as reliable internal control gene for qRT-PCR analysis in this study based on geNorm analysis results, hence FOXL2 relative mRNA expression levels were determined by comparison to the HPRT1 and plotted as ratio to HPRT1 expression values. Final results were expressed as fold differences in gene expression relative to the normalized calibrator, calculated by the ΔΔCt method as follows: n-fold = 2−(ΔCtsample−ΔCtcalibrator) where ΔCt values of the sample and calibrator were determined by subtracting the average threshold cycle (Ct) value of the transcript under investigation from the average Ct value of the HPRT1 gene for each sample. Raw Ct values were imported into GenEx Pro (version 4.3.5) and analyzed. All the primer pairs used in this study showed good overall efficiency (>96% and <106%) and excellent reproducibility of the amplification reactions.
Immunohistochemistry And Immunofluorescence
Immunohistochemistry was performed on 7-μm sections of paraffin-embedded samples including healthy endometrium, endometrium from patients with endometriosis, and endometrioma. Samples were washed in phosphate-buffered saline (PBS) 3 times, followed by endogenous peroxidase blockage by the use of 3% H2O2 in methanol solution. Again, samples were washed in PBS and incubated with 5% bovine serum albumin (BSA) in PBS for 1 hour to block nonspecific binding sites. The sections were incubated overnight at 4°C with the primary goat immunoglobulin (Ig) G anti-FOXL2 polyclonal antibody, diluted 1:40 in 0.1% BSA/PBS. Subsequently, sections were washed in PBS and then treated with HPR-conjugated antigoat IgG secondary antibody for 1 hour at room temperature (RT), followed by peroxidase streptavidin incubation in avidin–biotin complex reagent (Vectastain Elite Universal kit—Vector Laboratories, Burlingame, California). Diaminobenzidine (Sigma Chemicals CO, St Louis, Missouri) was subsequently used in order to visualize the immunolocalization of the primary antibodies. Tissue slices were counterstained with hematoxylin (Sigma Chemicals CO). Negative control slices were incubated with normal serum instead of primary antibodies.
For immunofluorescence, the samples were fixed in Histochoice (Amresco, Solon, Ohio) at RT for 4 hours. Sections were treated with 3% H2O2 for 1 hour and unmasked with citrate buffer solution (Thermo Fisher Scientific, Cheshire, United Kingdom) and 0.01 mol/L EDTA, pH 8. Slides were blocked with Dako Cytomation Protein Block Serum-free (Dako, Glostrup, Denmark) for 30 minutes at RT and then incubated overnight at 4°C with goat IgG anti-FOXL2 polyclonal antibody (ABCAM, Cambridge, United Kingdom) and mouse IgM anti-GMP-1 (SUMO-1) monoclonal antibody (Zymed Laboratories, San Francisco, California) diluted 1:40 and 1:20, respectively, in 0.1% BSA/PBS. After washing in PBS, the samples were incubated with secondary antibodies: Alexa Fluor 633 antimouse IgM and Alexa Fluor 488 antigoat IgG (Molecular Probes, Lifetechnologies Italia, Monza) diluted 1:500. Immunofluorescence analysis was performed using a Leica DMRB microscope (Leica Microsystems, Milano) equipped with fluorescence apparatus (488-633 excitation wavelength).
Statistical Analysis
The data were normally distributed and the results of the studies are reported as mean ± standard error of the mean. Statistical analysis was performed using the GenEx Pro and GraphPad Prism (GraphPad 4.0 Software Inc, San Diego, California). The selected reference gene HPRT1 was used to normalize Ct values, and quantities were calculated relative to the maximum Ct value. The comparison of 2 groups was performed with a Student t test while for analysis of more than 2 groups, a 1-way analysis of variance was performed, followed by the post hoc Tukey test. Correlation analyses were performed between FOXL2 and activin A expression levels, respectively, using the Pearson test. Statistical significance was set at P ≤ .05.
Results
Sequencing of FOXL2 Gene and Protein Localization in Healthy Endometrium
FOXL2 gene expression in healthy human endometrium as well as in endometrium of patients with endometriosis was demonstrated by qRT-PCR. The PCR products were first sequenced to confirm the expression of FOXL2 mRNA by the comparison with the reference sequence (NM_023067.3) from the National Center for Biotechnology Information database (data not shown).
Quantitative RT-PCR was performed to quantify mRNA levels. FOXL2 mRNA was expressed in both the groups, throughout menstrual cycle without any significant difference (Figure 1). When FOXL2 protein expression was evaluated in healthy endometrial tissue by immunohistochemistry, a FOXL2 staining was observed in the cytoplasm of glandular cells with a scattered distribution in stroma (Figure 2, healthy panel). Double-fluorescence staining revealed a colocalization of FOXL2 with SUMO-1, consistent with FOXL2 protein being SUMOylated and stabilized in the cytoplasm (Figure 2, healthy panel). Negative controls did not stain at all.

FOXL2 mRNA expression in human endometrium quantified by qRT-PCR during different phases of the endometrial cycle of healthy individuals or patients with endometriosis (healthy proliferative [HP]; healthy secretory [HS]; endometriosis proliferative [EP]; endometriosis secretive [ES]). Fold change (y axis) represents FOXL2 expression normalized to HPRT1, relative to P phase of healthy endometrium considered to be equal to 1. Bars indicate the standard error of the mean (SEM). Statistical analysis was performed using t test. mRNA indicates messenger RNA; qRT-PCR, quantitative real-time polymerase chain reaction; FOX, Forkhead box.

FOXL2 is expressed in human endometrium and colocalizes with SUMO-1. Representative micrographs of immunohistochemical staining of FOXL2 in human endometrium from women without (healthy panel) and with endometriosis (endometriosis panel) showed cytoplasmic staining in glandular and stromal cells. Immunofluorescence localization of FOXL2 and SUMO-1. Blue: DAPI; green: FOXL2; red: SUMO-1; yellow: merge (scale bars 100 µm). DAPI indicates 4′,6-diamidino-2-phenylindole; FOX, Forkhead box.
FOXL2 mRNA Expression And Protein Localization in Endometrium of Patients With Endometriosis and in Endometriotic Tissue
FOXL2 mRNA endometrial expression was not different in patients with endometriosis during endometrial cycle (Figure 1). Immunostaining of eutopic endometrium showed a cytoplasmic localization more concentrated in the glandular epithelium while diffuse in the stromal cells in endometriotic samples. SUMO-1–FOXL2 double staining showed a specific colocalization with the same pattern (Figure 2).
As shown in Figure 3, FOXL2 mRNA level was significantly higher in ovarian endometriosis than in healthy endometrium (2.63-fold induction, P < .05) or in eutopic endometrium with endometriosis (3.27-fold induction, P < .01; Figure 3, panel A).

FOXL2 mRNA expression in healthy endometrium, eutopic endometrium from women with endometriosis, ectopic endometrium from women with endometrioma. A, Fold change (y axis) represents FOXL2 expression normalized to HPRT1, relative to healthy endometrium considered to be equal to 1. A statistical significance was found for endometrioma when compared with healthy endometrium (*P < .05) and endometrium from endometriosis (°P < .01). Bars indicate the standard error of the mean (SEM). Statistical analysis was performed using t test. FOXL2 protein localization in endometriotic lesions. B, Representative micrographs of immunohistochemical staining of FOXL2 in endometrioma, cytoplasmic FOXL2 staining was observed in glandular and stromal cells (A, scale bar 100 µm; B, scale bar 25 µm). mRNA indicates messenger RNA; HPRT1, hypoxanthine phosphoribosyltransferase 1; FOX, Forkhead box.
The presence of FOXL2 protein in endometriotic lesions was confirmed by immunohistochemistry (Figure 3, panel B). FOXL2 mRNA expression levels correlated with the activin A expression in endometriosis (r = .635).
Discussion
The present study showed that (1) in healthy women, FOXL2 mRNA and protein are expressed and localized in endometrium, without significant changes during endometrial cycle; (2) in patients with endometriosis, the endometrial pattern is not different; and (3) FOXL2 is hyperexpressed in endometriotic tissue and correlated with activin A in endometriosis.
In accordance with data showing an expression and a possible role of FOX factors in human endometrial function, 4,26,27,32 our present data confirm an endometrial expression and reinforce the involvement of FOXL2 in reproductive function and disease. FOXL2 is mainly recognized as a regulatory gene involved in ovarian function 5,12,13 and our results showed FOXL2 mRNA expression and its protein localization in endometrium. The lack of significant cycle-related variations for endometrial FOXL2 mRNA expression suggests an independent sex steroid hormone expression of this gene. The immunohistochemical and immunofluorescence studies showed a stromal and glandular localization and the colocalization with SUMO-1, confirming data obtained in ovary, 32 and suggest that SUMO-1 pathway may modulate a local role of FOXL2 transcription factor inside endometrial cells. Indeed, it is known that members of FOX family can interact directly with PR, through a transcription factor “crosstalk,” depending on protein–protein interactions between steroid receptors and other sequence-specific transcription factorsand is induced and activated by cyclic adenosine monophosphate (cAMP) during decidualization of endometrial cells. 14,29,33 Cyclic adenosine monophosphate signaling profoundly alters the SUMO cycle in endometrial stromal cells in a time-dependent manner, leading to hypoSUMOylation of the activated PR-A and transcriptional activation of decidua-specific genes. 29
The role of FOXL2 in reproductive axis function is suggested by the evidence that FOXL2 has been shown to modulate ovarian 10,11,34 and pituitary gonadotrope activity. 7 -9 In particular, FOXL2 overexpression potentiates activin induction in gonadotropes and confers FSHβ responsiveness in heterologous cells where this promoter is normally refractory to activin induction. 9 A FOXL2 mutation has also been described in ovarian 35 and in endometrial 29 carcinoma supporting a clinical relevance. Activin induces the transcription of follistatin by the action of SMAD3 and SMAD-binding element 1. FOXL2 strongly binds a forkhead-binding element (FKHB), thus modulating the function of this partner of SMAD3. 36 Recently, it has been demonstrated that in rodents, FOXL2 is necessary for the full synergistic response between activin and progestins, interacting with PR directly and/or indirectly via SMAD3. 37 Therefore, FOXL2 may act both as a regulator of transcription and as a coordinator of the action of SMAD3 by modulating action of activin A.
The possible involvement of forkhead box family members in endometriosis is also suggested by other studies. 38 In particular, FOXL2 modulates molecules involved in inflammatory response and immunomodulation, 39 like interferon β1, interleukin (IL) 12A and IL29, and intercellular adhesion molecule 1. An effect on apoptosis is suggested by the increased expression of B-cell lymphoma 2-related protein A1, a direct transcription target of nuclear factor κB, which efficiently suppresses apoptosis, or immediate early response-3, able to protect from tumor necrosis factor-induced apoptosis. 39
The hyperexpression of FOXL2 in endometriotic tissue supports that this transcription factor is disregulated in endometriotic cells. Previous studies described a FOXL2 mutation in endometrium of women with endometriosis. 27,28 The upregulation of FOXL2 in endometriosis supports that these cells are highly proliferative and reactive in response to an inflammatory environment, 40 confirming data from a transcriptomic study showing mutations of FOXL2 in endometriotic tissue 28 and corroborates the hypothesis that the deregulation of this transcription factor may be associated with pathogenesis of endometriosis.
In conclusion, the present study showed that endometrial glandular and stromal cells express FOXL2, throughout the menstrual cycle and high FOXL2 mRNA expression was found in endometriotic tissue, suggesting a role in endometrial function and in pathogenic mechanisms of endometriosis.
Footnotes
Acknowledgments
We thank Dr Giuseppe Belmonte for technical support in protein immunolocalization and images selection.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
