Oestrogen receptor beta isoform expression in sporadic colorectal cancer, familial adenomatous polyposis and progressive stages of colorectal cancer
BMC Cancer volume 17, Article number: 754 (2017)
Among the sex hormones, oestrogen may play a role in colorectal cancer, particularly in conjunction with oestrogen receptor-β (ERβ). The expression of ERβ isoform variants and their correlations with familial adenomatous polyposis (FAP) syndrome and sporadic colorectal carcinomas are poorly described.
This study aimed to investigate the expression levels of the ERβ1, ERβ2, ERβ4 and ERβ5 isoform variants using quantitative RT-PCR (921 analyses) in FAP, normal mucosa, adenomatous polyps and sporadic colorectal carcinomas.
Decreased expression of ERβ isoforms was identified in sporadic polyps and in sporadic colorectal cancer as well as in polyps from FAP syndrome patients compared with normal tissues (p < 0.001). In FAP patients, ERβ1 and ERβ5 isoforms showed significant down-expression in polyps (p < 0.001) compared with matched normal tissues. However, no differences were observed when sporadic colorectal carcinomas were compared to normal mucosa tissues. These findings suggest an association of the ERβ isoform variants in individuals affected by germline mutations of the APC gene. Progressively decreased expression of ERβ was found in polyps at early stages of low-grade dysplasia, followed by T1-T2 and T3-T4 tumours (p < 0.05). In sporadic colorectal cancer, the loss of expression was an independent predictor of recurrence, and ERβ1 and ERβ5 expression levels were associated with better disease-free survival (p = 0.002).
These findings may provide a better understanding of oestrogens and their potential preventive and therapeutic effects on sporadic colorectal cancer and cancers associated with FAP syndrome.
Colorectal cancer (CRC), the second and third most commonly diagnosed cancer in women and men, respectively, presents a variable geographical incidence that has been attributed to differences in diet and environmental exposure . These variations occur in conjunction with a genetically determined background of susceptibility. The incidence of CRC is substantially higher in men and women greater than 50 years of age, although the basis for this difference is unknown. Sex steroids have been considered a contributing factor because high parity, early age at first pregnancy, oral contraceptive use and oestrogen replacement therapy are associated with decreased risk of CRC [2,3,4,5]. In addition, younger women (<45 years of age) with CRC have better overall survival than men of the same age .
Familial adenomatous polyposis syndrome (FAP), a precancerous colorectal condition, is an inherited disease caused by a germline mutation in the adenomatous polyposis coli (APC) gene . In FAP, many neoplastic lesions of different stages are typically found within the same individual, thus providing a better understanding of the adenoma-carcinoma progression. The fact that these polyps appear in the second decade of life (i.e., during puberty) and show increased size and number during adolescence (i.e., at the peak of production) also suggests that sex steroids may act as a cofactor in the development of colorectal polyps and in carcinogenesis linked to FAP .
Evidence from several studies suggests that among the sex hormones, oestrogen may play a role in colorectal cancer, particularly through oestrogen receptor-β (ERβ) [9,10,11,12,13]. However, since the discovery of this receptor in 1996 , research efforts have been focused on defining its biological function, which remains poorly understood. Several alternative splicing isoforms of the oestrogen receptor beta gene (ESR2) occur (transcript variants a, b, d, f, g, h, i, k, l; ERβ isoforms 1–5), and these isoforms have been detected in both normal and malignant cells [15,16,17,18]. The expression of ESR2 splicing variants in colon cancer cells was first reported in 2001 . To our knowledge, this is the first study in which the expression levels of ESR2 variants in FAP, adenomatous polyps and sporadic colorectal carcinomas have been evaluated and compared with normal mucosa, tumour stage and prognosis data. The goal of this study was to provide a better understanding of oestrogens and their potential effects on sporadic and hereditary colorectal tumours.
Fresh frozen tissue samples from 98 patients with sporadic CRC were retrospectively collected in the Institutional biobank between the years 2005 and 2016; 48 of the patients had stage T1/T2 tumours, and 50 had stage T3/T4 tumours. As a reference, paired histologically normal mucosa was obtained from 52 cases. Formalin-fixed, paraffin-embedded (FFPE) tissue samples were obtained from 52 sporadic polyps. In addition, 64 FFPE polyps and 41 FFPE adjacent normal tissues were obtained from 41 FAP patients, between the years 2000 and 2016. All specimens were submitted to macrodissection and histology confirmation. Tumor samples were composed of at least 70% of epithelial cells, and normal tissues presented more than 90% of epithelial cells. In total, 307 samples were evaluated, with 921 analyses of three transcripts and their associated isoform pairs.
Patients undergoing neoadjuvant treatment and those with other forms of hereditary CRC or inflammatory bowel disease were excluded from the study. The Institutional Review Board of the A.C. Camargo Cancer Center, Sao Paulo, Brazil, approved this study (CEP01453/10).
The histopathological classification of the tumours and the clinical staging followed the recommendations of the WHO International Classification of Diseases for Oncology  and the Tumour-Node-Metastasis staging system (TNM) , respectively. The medical records of the patients were examined to obtain detailed clinical and pathological data. FAP was diagnosed in individuals with more than 100 adenomatous colon polyps or with APC mutations.
The routine follow-up of patients with sporadic colorectal tumours treated at the institution included quarterly clinic visits during the first two years, with laboratory tests, assessment of tumour markers, chest X-ray and ultrasound or abdominal sectional imaging (tomography) performed at alternating visits. These procedures were performed every six months from the third to the fifth year and annually thereafter. Colonoscopy was performed at the first, third and fifth years postoperatively. In the event of recurrence, tests were requested for re-staging and treatment planning. Positron emission tomography (PET-CT) was performed in cases considered to have higher risk of metastasis.
The mean follow-up time was 70.1 months, and the median was 58.6 months (1.67 to 188.7 months). The FAP patients were carriers of the classic phenotype (>100 adenomatous polyps) and 53.6% were male with a median age at syndrome diagnosis of 31 years (19 to 56 year old). The polyps were located in the colon, and histology analysis showed low-grade dysplasia. Demographic, clinical and pathological characteristics and treatment of individuals with sporadic cancer (N = 98) are summarised in Table 1.
Total RNA was extracted from macrodissected fresh frozen tissues and FFPE tissue samples using a QIAsymphony RNA Kit (Qiagen) and a RecoverAll Total Nucleic Acid Isolation Kit (Ambion) according to the manufacturer instructions. The RNA quantity and quality from tissue samples (FFPE and fresh frozen) were evaluated using a NanoDrop ND-1000 Spectrophotometer (v.3.0.1, Labtrade). The RNA quality from fresh frozen tissue samples was also evaluated using the Bioanalyser Agilent RNA 6000 Nano LabChip kit (Agilent 2100 bioanalyser) (Additional file 1: Table S1).
Selection of the reference genes
The Cancer Genome Atlas (TCGA: level 3 colon tumour RNA sequencing database) was used to select the reference genes for RT-qPCR normalisation. The most stable genes were identified among the 30 potential reference transcripts (Applied Biosystems TLDA test) using 41 surrounding normal tissues and 285 colon adenocarcinomas. A standard deviation (SD) ranking was conducted for each gene in all samples, and the p-value (unpaired t-test for unequal variances) was calculated comparing the normal and tumour samples. The lowest SDs and highest p-values were ranked; the PUM1, POP4 and EIF2B1 genes had the highest rankings (lowest variation between samples and without differences between the normal and tumour groups). The expression levels of ESR2 were evaluated by RT-qPCR using the PUM1, POP4 and EIF2B1 genes as a reference.
Reverse transcription quantitative polymerase chain reaction (RT-qPCR)
Three pair of primers flanking the ESR1 gene were designed. Primer pair 1 amplifies the transcript variants a, b, d, f, k and l (ERβ1, ERβ2 and ERβ4 isoforms); primer pair 2 amplifies transcripts a and g (ERβ1 and ERβ5 isoforms), and primer pair 3 amplifies the transcript variants b, l and k, which are translated into the ERβ2 isoform (Table 2 and Fig. 1). Three endogenous references (PUM1, EIF2B1, and POP4) were used in the RT-qPCR assays (Table 1). The primers were designed using Primer-Blast (available at http://www.ncbi.nlm.nih.gov/tools/primer-blast/), and Primer 3 was selected based on the study by Yang et al., 2016 .
Total RNA samples were digested with 1 U of DNase I (amplification grade, Life Technologies) in 10X DNase I Reaction Buffer. The reactions were performed in a PTC-100 thermal cycler (Peltier-Effect Cycling - MJ Research) for 15 min at room temperature; the enzyme was then inactivated by heating at 70 °C for 10 min. The cDNA synthesis was performed in a final volume of 20 μL containing 5X first-strand buffer (250 mM Tris–HCl pH 8.3, 375 mM KCl and 15 mM MgCl2), 10 mM of each dNTP, 0.5 μg/μL oligo (dT), 0.1 M dithiothreitol and 200 U of reverse transcriptase (SuperScriptTM III, Invitrogen). Reverse transcription was performed at 42 °C for 60 min followed by inactivation at 70 °C for 15 min. The ABI Prism 7900 Sequence Detection System automatic thermocycler was used for RT-qPCR amplification. All reactions were assembled by robotic pipetting into 384-well plates with QIAgility (Qiagen, Courtaboeuf, France) in a total volume of 12.5 μL of Power SYBR Green Master Mix (Applied Biosystems, Foster City, CA, USA), 20 ng of cDNA and 200 nM of each primer. All samples were analysed in duplicate and subjected to the following cycling conditions: an initial temperature of 95 °C for 10 min and 45 cycles of 95 °C for 15 s and 60 °C for 1 min. The amplification quality was verified by analysing the dissociation curve (specificity) to discriminate primer-dimers from low levels of transcript expression. A duplicate of no template control (NTC) was included in each PCR amplification for each primer pair, all showing negative signals (Cq = 45 or small dimer amplification detection Cq > 42). The values obtained for all samples were normalised by determining the ratio of the gene of interest to the reference gene. The mean Cq (quantification cycle) was used for normalisation. Samples with a mean Cq greater than the mean + 1 SD of the geometric mean of the Cqs of the three reference genes were excluded (8 fresh frozen and 30 FFPE samples). The model proposed by Pfaffl  was used for data normalisation.
ESR2 expression in the cancer genome atlas (TCGA) database of colon cancer
ESR2 expression data in colon carcinomas versus normal tissues (log2 + 1 reads per million) was assessed from The Cancer Genome Atlas (TCGA: level 3 RNA sequencing database) (t test). An isoform specific analysis of ESR2 was implemented using Isoform Expression View (normal-tumor comparison) of the ISOexpresso tool (http://wiki.tgilab.org/ISOexpresso/) .
The findings obtained from five groups of tissues were compared according to the tissue type: fresh frozen or FFPE. The variants expression levels of fresh frozen sporadic colorectal tumours (ST: 98 samples) were compared with normal mucosa (NS: 52 samples). Similarly, the expression levels of FFPE samples including sporadic polyps (SP: 52 samples), normal mucosa from FAP patients (NFAP: 41 samples) and adenomatous polyps from FAP patients (PFAP: 64 samples) were compared. A comparison with all groups (NS, ST, SP, NFAP and PFAP) was also performed, independent of the tumour type (fresh frozen and FFPE) (Additional file).
Descriptive statistics were used to characterise the sample. A frequency distribution was used to describe categorical variables and measures of central tendency and variability for numerical or continuous variables. The chi-squared frequency test was used to compare categorical variables, and Fisher’s exact test was used to verify associations. Student’s t test (paired and unpaired) was used to compare continuous variables. Analysis of variance (ANOVA) was used to compare continuous variables of multiple groups. The Kaplan-Meier method was adopted to estimate the overall and disease-free survival probabilities. The difference between survival curves of the same variable was evaluated using the log-rank test. The cut-off for relative gene expression was determined using the log-rank test (maximally selected rank statistics in R) . Multivariate analysis was used to predict the combined effect of independent variables on survival using the Cox model with proportional risks. The data were statistically analysed using Statistical Package for Social Sciences (SPSS) version 20.0 and GraphPad Prism 5.0 (GraphPad Software Inc., La Jolla, CA). The significance level for all statistical tests was 5%.
A significant decrease of ESR2 expression levels was observed in polyps and tumours compared to normal mucosa (fresh frozen tissues) or normal tissue from FAP patients (FFPE samples), regardless of heredity (Fig. 2A). However, differences in expression levels were found according to the group tested and the isoforms.
Primer pair 1: Transcript variants a, b, d, f, k, and l: ERβ1, ERβ2 and ERβ4 isoforms
The levels of transcript variants a, b, d, f, k, and l, corresponding to the ERβ1, ERβ2 and ERβ4 isoforms, were significantly lower in sporadic tumours (ST) (p < 0.001; Mean Diff 3.285; 95% CI 2.780 to 3.791) compared with normal sporadic mucosa (NS) (fresh tissue) (Fig. 2a).
In FAP syndrome patients, lower expression levels of ERβ (a) mRNA were found in polyps compared with paired FAP mucosa samples (p < 0.05; Mean Diff 3.338; q = 3.854; 95% CI 0.3677 to 6.308). Sporadic polyps (SP) were also down-expressed when compared to FAP mucosa samples p < 0.05; Mean Diff 3.1; q = 3.883; 95% CI: 0.3621 to 5.839). No difference was observed in sporadic polyps (SP) versus PFAP (Fig. 2a).
Comparing all FFPE and fresh frozen tissue samples, no significant differences in mRNA expression levels were found in sporadic normal mucosa and normal FAP mucosa. However, invasive polyps and tumours showed lower ESR2 expression levels compared with normal mucosa from both, sporadic tumours and FAP patients (Additonal file 2: Figure S1).
Primer pair 2: Transcript variants a and g: ERβ1 and ERβ5 isoforms
No significant differences were found in the expression levels of mRNA of ERβ1 and ERβ5 isoforms compared with normal mucosa (NS) and tumour tissues (ST) from fresh frozen samples (Fig. 2a). However, in FFPE tissues, the normal mucosa from FAP syndrome patients showed high expression levels of ERβ1 and ERβ5 isoforms compared with the sporadic polyps (p < 0.001; Mean Diff 7.568; q = 7.381; 95% CI 4.054 to 11.08) and FAP polyps (p < 0.001; Mean Diff 6.4; q = 5.721; 95% CI 2.566 to 10.23) (Fig. 2a).
The analysis of all groups together (Additional file 2: Figure S1) revealed downexpression levels ERβ1 and β5 isoforms in sporadic tumours compared to normal mucosa of FAP patients (p < 0.001; Mean Diff −5.926; q = 8.04; 95% CI −8.81 to −3.04). Furthermore, the normal mucosa of FAP syndrome patients showed a higher concentration of ERβ1 and ERβ5 isoforms than the other tissues analysed.
Primer pair 3: Transcript variants b, i and k: ERβ2 isoform
ERβ2 expression levels were significantly lower in tumour tissues (p < 0.001; Mean Diff 1.718; 95% CI 1.161 to 2.275) than in sporadic normal mucosa (Fresh frozen tissue). Similarly, for FFPE tissues, lower ERβ2 mRNA expression levels were found in FAP polyps (p < 0.05; Mean Diff 2.944; q = 4.393; 95% CI 0.6458 to 5.243) and sporadic polyps (p < 0.05; Mean Diff 2.438; q = 3.945; 95% CI 0.3186 to 4.557) compared with FAP mucosa.
Considering all groups together (Additional file 2: Figure S1), sporadic normal mucosa and normal FAP mucosa showed no significant differences in ERβ2 mRNA levels. No significant differences in ERβ2 expression between sporadic polyps and tumour tissues were detected.
ERβ expression according to the cancer genome atlas (TCGA)
RNA sequencing data from 41 normal mucosa and 285 colon adenocarcinomas were extracted from the TCGA database. A significant decrease in ERβ expression level was found in sporadic tumours (p < 0.001; Mean Diff 0.984 ± 0.146; 95% CI 0.697 to 1.271) (Fig. 2b). The isoform specific analysis (http://wiki.tgilab.org/ISOexpresso/) highlighted low expression levels of the variants d (ERβ4) and b (ERβ2) in tumour samples (Additional file 3: Figure S2).
ERβ expression according to T stage (TNM)
T3/T4 sporadic tumours presented significantly decreased expression levels of ERβ when compared with T1/T2 tumours (p = 0.02; Mean Diff 0.87; 95% CI 0.12 to 1.63; Student’s t test). Decreased expression levels were observed only when Primer 1 was used, which amplifies a greater number of transcription variants (the a, b, d, f, k, and l variants, which produce the isoforms ERβ1, ERβ2 and ERβ4) (Fig. 2c). The comparison of T stage (TNM) (T1 vs T2 vs T3 vs T4) in the four classes revealed no significant differences (one-way ANOVA and Tukey’s post hoc test).
Disease-free survival and overall survival analysis
The time to occurrence of relapse ranged from 4.7 to 36.4 months (median 15.4 months and mean 15.6 months). The probability of recurrence-free survival was 77.4% for both 5- and 10-year survival. By univariate analysis, disease-free survival was greater in subjects ≥65 years of age and in those with normal serum carcinoembryonic antigen levels, T1-T2 tumours, absence of blood embolization, absence of lymph node or metastatic invasion, well differentiated tumour histological grade or Primer 2 (ERβ1 and ERβ5) expression levels ≤ −4.5. (Table 3, Additional file 4: Figure S3). The multivariate analysis revealed that T and N stage and ERβ expression levels ≤ −4.5 are independent predictors of disease-free survival (Table 4).
The cumulative probability of overall survival was 83.9% at 5 years and 77.1% at 10 years. The univariate analysis revealed greater overall survival in patients with normal serum carcinoembryonic antigen levels, T1-T2 tumours and the absence of perineural, lymph node or metastatic invasion. ERβ expression levels were not statistically significant in predicting overall survival (Table 3).
The effects of oestrogens on tissue are mediated by members of the nuclear oestrogen receptor subfamilies ERα and ERβ. In normal colorectal tissues, ERβ is predominant and appears to play an important role in the biological mechanisms of action of sex steroids [25, 26]; its expression decreases at higher Dukes’ stages, which display little or no ERα expression [9,10,11, 19]. Despite being important in the development of breast cancer, ERα is found in low levels in normal colorectal tissue . Using different strategies, several studies reported no differences between the expression levels of ERα comparing colorectal tumours and normal mucosa [12, 13, 27,28,29]. Therefore, the ERβ variant isoforms may be the key to understand the ERβ signalling in different tissues .
The expression of splicing variants in colon cancer cells was first reported by Campbell-Thomson et al. in 2001 . In normal human colon tissue, ERβ1, ERβ2 and ERβ5 are expressed, whereas ERβ3 and ERβ4 have been found only in testicles. In 2005, Wong et al. reported the expression of ERβ isoforms in colorectal carcinomas . In that study, the expression levels of ERβ1 and ERβ2 were completely lost in 22% and 49% of primary colorectal carcinomas, respectively. In contrast, a truncated isoform referred to as ERβ5 was found in all colorectal carcinomas and in tumour cell lines. Despite the limited number of cases, the authors suggested a prognostic role of ERβ1 and β2 isoforms.
A few alternative ERβ1 isoforms (ERβ2/ERβcx, ERβ3, ERβ4 and ERβ5) that result from alternative splicing of the last receptor-encoding exon have been reported. These proteins have a truncated domain or are otherwise altered in the C-terminal region . The C-terminal region of the ER contains a ligand-binding domain that is involved in transcription activation, receptor dimerization, nuclear translocation and interaction of the receptor with transcription co-regulators . Therefore, modifications in this region can affect the activity and biological function of the ER. Leung et al. (2006) reported that ERβ2–4-5 isoforms are deficient in intrinsic ligand-dependent transactivation activity, making ERβ1 (wild type) the only functional receptor isoform . However, ERβ2–4-5 isoforms form heterodimers with ERβ1 and may increase ERβ1-induced transcription activity at physiological oestrogen concentrations.
In tissues such as prostate, breast and endometrium, a relationship between oestrogen signalling and the initiation of carcinogenesis, cancer progression and therapeutic response have been suggested [16, 33,34,35,36,37,38]. Oestrogenic action and the quantitative variation in the proportion of ERβ in different tissues indicate that the ERβ subtypes have different functions [35, 39]. Although the specific roles of ERβ subtypes in colorectal cancer are still uncertain, it is plausible that the varying co-expression of different isoform groups in normal mucosa, polyps and tumours of the colon, as shown in this study, can alter carcinogenesis. These findings should be taken into account in the design of future studies and screening protocols. In addition, the frequency of screening examinations should be re-evaluated in subjects who are undergoing hormone replacement or hormone deprivation therapy, treatments that are commonly used in breast cancer. These patients have a higher risk of developing colorectal cancer, and the hormonal exposure they receive may be involved in this increased risk. The effect of the use of hormone replacement therapy in patients with CRC or FAP syndrome must therefore be carefully evaluated.
The analysis to evaluate ERβ isoforms levels was performed in fresh frozen tissues (sporadic tumours) and in FFPE samples (sporadic and FAP polyps). One limitation of our study is the use of FFPE samples, which presents RNA with lower quality compared to fresh frozen tissue samples. However, the comparison of the expression levels of different ERβ isoforms using normal mucosa (fresh frozen tissue) and normal FAP mucosa (FFPE) revealed decreased ERβ expression in tumour samples and polyps. This result was confirmed grouping all samples independently of the tissue type (fresh tissue or FFPE). To validate the data, RNA-Seq data from 285 colon adenocarcinomas and 41 non-tumour colon tissue samples extracted from the TCGA dataset were used. In agreement with our data, this analysis revealed that mRNA ERβ expression levels were significantly lower in tumour tissue than in non-tumorous tissue samples (p < 0.001; Mean Diff 0.984 ± 0.146; 95% CI 0.697 to 1.271). Nguyen-Vu et al. (2016) using RNA-seq data of 233 colon adenocarcinomas and 21 non-tumor colon tissues from The Cancer Genome Atlas (TCGA) dataset. The ERβ expression was decreased in the cancerous state compared to non-cancerous tissues .
In our study, ERβ1, ERβ2 and ERβ4 variants (primer sets 1 and 3) showed downexpression in ST compared to NS. Similarly, ERβ1, ERβ2, ERβ4 and ERβ5 variants (primer sets 1, 2 and 3) presented downexpression in FAP polyps compared with NFAP. These findings were confirmed with TCGA dataset results in ERβ2 and ERβ4 performed in sporadic tumours. However, these transcripts demonstrated very low expression by RNA sequencing and some variants was not even detected by TCGA (i.e.: ERβ1 and ERβ5). Although RNA sequencing is a suitable method to identify and discriminate mRNA variants, the sensitivity of RT-qPCR is superior, which is considered a gold standard methodology to evaluate transcripts [41,42,43].
Primer pair 2, which amplifies the transcript variants a and g (encoding ERβ1 and ERβ5) showed the highest expression levels of ESR2 in normal FAP mucosa (NFAP). FAP polyps emerge during puberty and increase in number during adolescence, which is the period of peak oestrogen production. Based on these findings, it appears that the transcript variants have the potential to act as cofactors in the development of polyps and FAP-related colorectal carcinogenesis. Interesting, comparing sporadic tumours (ST) with sporadic normal mucosa (NS), the ERβ1 and ERβ5 levels were not significant.
A primary chemoprevention trial  consisting of a double-blinded, randomised four-year study using sulindac included subjects with APC mutations as well as phenotypically unaffected subjects. The authors reported total eradication of the polyps in the placebo group, which included one patient who received an occasional administration of oral contraceptives and who had developed polyps . The patient follow-up comprised flexible rectosigmoidoscopy every four months for 48 months, and an increased prevalence and recurrence of polyps was observed after the discontinuation of oral contraceptives.
Selective oestrogen receptor modulators (SORMs or SERMs) such as raloxifene and tamoxifen have both oestrogenic and anti-oestrogenic actions depending on the tissue type [44,45,46]. For example, raloxifene is not only effective in preventing osteoporosis  but has also been shown to be as effective as the SORM archetype tamoxifen in preventing breast cancer ; however, it increases the risk of uterine adenocarcinoma . Despite the widespread clinical application of these modulators, very little is known about how they may affect colon cancer. It is important to identify oestrogen receptor variant subtypes and to determine their relative expression levels at different stages of tumorigenesis, especially in different tissues, because specific oestrogen receptor subtypes can be involved in the specific agonistic or antagonistic actions of oestrogens. The next step is therefore to determine the correlation between the predominance of specific isoforms and their interactions with various forms of oestrogen (isoflavones, SORM and hormone replacement therapy). The incidence of colorectal cancer is much lower in Asian countries than in the western world, a fact that may be explained by the high intake of phyto-oestrogens, such as soy, in Asian cultures . This intake offers a protective effect that may be associated with the mediation of ER binding to genistein, which is the main phyto-oestrogen in soy and is thought to be a tyrosine kinase inhibitor, a characteristic that may explain some of its action on the colon [51, 52].
When invasive tumours were analysed individually, no differences in expression associated with the degree of tumour differentiation were found. This result likely occurred because the majority (87.8%) of the tumours studied herein displayed a moderate degree of differentiation. Immunohistochemistry analysis using an ERβ monoclonal antibody revealed decreased expression in higher grade or more dedifferentiated tumours [11, 13].
Although the results of a comparison of individual T stages (T1 vs T2 vs T3 vs T4) were not significant, most likely due to the limited number of cases at stages T1 and T4, the grouped analysis revealed that larger (T3/T4) tumours presented decreased ERβ expression levels compared to T1/T2 tumours. These alterations were more evident and significant for the primer pair 1, which amplified the largest number of transcription variants (a, b, d, f, k and l), representing ERβ1, ERβ2 and ERβ4 isoforms. The amplification with the primers 2 and 3 resulted in a significant reduction of the expression levels in normal mucosa to very low values during the colorectal carcinogenesis, which may explain the absence of differences between T1/T2 and T3/T4 tumours. One plausible explanation is that the primer 1 is more sensitive for the detection of these differences, because contains a group with more homologous variants. The loss of ERβ expression was significantly different in adenomatous polyps with low-grade dysplasia, illustrating that differences in the expression of these receptors occur at early stages of carcinogenesis.
Using an expression level threshold of −4.5, we found that low expression levels of the ERβ1 and ERβ5 isoforms (transcript variants a and g: primer pair 2) were independent factors associated with sporadic colorectal cancer recurrence (ERβ > −4.5 presented a 4.08-fold higher risk of relapse (CI 1.69–9.84); p = 0.002). Paradoxically, this sample set of cases displayed no significant difference in loss of normal mucosa ERβ expression. With respect to polyps and smaller and larger sporadic tumours, this subgroup was different in tissues from patients with FAP syndrome; therefore, it is plausible that this isoform group may be correlated with the APC suppressor gene, which is mutated in FAP syndrome. In contrast, sporadic cases are often associated with other variables associated with somatic carcinogenesis.
Cases of FAP cancer were not included in this study, which was designed to evaluate the prognosis of FAP tumours, because FAP is a very rare syndrome and therefore involves a heterogeneous group of patients; as a result, cases of FAP have been historically treated in several ways, which may lead to uncertainties in the prognostic results. Unfortunately, the gold standard for preventive treatment of colorectal cancer in FAP individuals is still total prophylactic colectomy at a young age, preferably prior to the malignant transformation of polyps . Thus, we investigated mRNA expression of ERβ isoforms in normal mucosa and polyps. Significant loss was demonstrated in all subgroups of ERβ isoforms in adenomatous polyps of PAF compared to normal mucosa (p < 0.001).
Recently, an alternative mechanism for oestrogenic action was described where the G protein-coupled oestrogen receptor (GPER) mediates the rapid non-genomic effects of oestrogen, phyto-oestrogen and xeno-oestrogen . GPER activation may inhibit the growth of CRC cells, both in vitro and in vivo, through multiple intracellular signalling pathways. Liu et al. (2017) showed that GPER expression in tumour tissue was markedly lower than in corresponding adjacent normal mucosal tissue. Tumours expressing lower levels of GPER exhibited a significantly lower survival rate than those with higher GPER expression levels . The physiological significance of these rapid effects and their integration with nuclear responses to oestrogen are important issues that still need to be further investigated.
Overall, we demonstrated that the mRNA expression levels of ERβ isoforms are downregulated in sporadic colorectal cancer and in FAP individuals. T3/T4 tumours also presented decreased expression of ERβ. Additionally, the expression levels of ERβ1 and ERβ5 were associated with the probability of disease-free survival. These differences may inform new clinical studies aimed at preventative strategies, especially in groups of patients who are receiving hormone therapy or are under conditions of hormonal deprivation. Although the limitation of our study was the analysis of the gene expression levels in a set of FFEP specimens, our data pointed out that the deregulation of ESR2 isoform variants may be associated with colorectal cancer progression.
Oestrogen Receptor Beta
Oestrogen Receptor Beta Gene
Familial Adenomatous Polyposis
Normal Mucosa from FAP Patients
Sporadic Normal Mucosa
Adenomatous Polyps from FAP Patients
The Cancer Genome Atlas
Tumour-Node-Metastasis Staging System
Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global cancer statistics, 2012. CA Cancer J Clin. 2015;65:87–108.
Ferlay J, Soerjomataram I, Dikshit R, Eser S, Mathers C, Rebelo M, et al. Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012. Int J Cancer. 2015;136:E359–86.
Chlebowski RT, Wactawski-Wende J, Ritenbaugh C, Hubbell FA, Ascensao J, Rodabough RJ, et al. Estrogen plus progestin and colorectal cancer in postmenopausal women. N Engl J Med. 2004;350:991–1004.
Simon MS, Chlebowski RT, Wactawski-Wende J, Johnson KC, Muskovitz A, Kato I, et al. Estrogen plus progestin and colorectal cancer incidence and mortality. J Clin Oncol. 2012;30:3983–90.
Barzi A, Lenz AM, Labonte MJ, Lenz HJ. Molecular pathways: estrogen pathway in colorectal cancer. Clin Cancer Res. 2013;19:5842–8.
Hendifar A, Yang D, Lenz F, Lurje G, Pohl A, Lenz C, et al. Gender disparities in metastatic colorectal cancer survival. Clin Cancer Res. 2009;15:6391–7.
Half E, Bercovich D, Rozen P. Familial adenomatous polyposis. Orphanet J Rare Dis. 2009;4:22.
Vasen HF, Moslein G, Alonso A, Aretz S, Bernstein I, Bertario L, et al. Guidelines for the clinical management of familial adenomatous polyposis (FAP). Gut. 2008;57:704–13.
Foley EF, Jazaeri AA, Shupnik MA, Jazaeri O, Rice LW. Selective loss of estrogen receptor beta in malignant human colon. Cancer Res. 2000;60:245–8.
Jassam N, Bell SM, Speirs V, Quirke P. Loss of expression of oestrogen receptor beta in colon cancer and its association with Dukes' staging. Oncol Rep. 2005;14:17–21.
Konstantinopoulos PA, Kominea A, Vandoros G, Sykiotis GP, Andricopoulos P, Varakis I, et al. Oestrogen receptor beta (ERbeta) is abundantly expressed in normal colonic mucosa, but declines in colon adenocarcinoma paralleling the tumour's dedifferentiation. Eur J Cancer. 2003;39:1251–8.
Nussler NC, Reinbacher K, Shanny N, Schirmeier A, Glanemann M, Neuhaus P, et al. Sex-specific differences in the expression levels of estrogen receptor subtypes in colorectal cancer. Gend Med. 2008;5:209–17.
Rudolph A, Toth C, Hoffmeister M, Roth W, Herpel E, Jansen L, et al. Expression of oestrogen receptor beta and prognosis of colorectal cancer. Br J Cancer. 2012;107:831–9.
Kuiper GG, Enmark E, Pelto-Huikko M, Nilsson S, Gustafsson JA. Cloning of a novel receptor expressed in rat prostate and ovary. Proc Natl Acad Sci U S A. 1996;93:5925–30.
Poola I, Abraham J, Baldwin K. Identification of ten exon deleted ERbeta mRNAs in human ovary, breast, uterus and bone tissues: alternate splicing pattern of estrogen receptor beta mRNA is distinct from that of estrogen receptor alpha. FEBS Lett. 2002;516:133–8.
Leung YK, Lam HM, Wu S, Song D, Levin L, Cheng L, et al. Estrogen receptor beta2 and beta5 are associated with poor prognosis in prostate cancer, and promote cancer cell migration and invasion. Endocr Relat Cancer. 2010;17:675–89.
Leung YK, Mak P, Hassan S, Ho SM. Estrogen receptor (ER)-beta isoforms: a key to understanding ER-beta signaling. Proc Natl Acad Sci U S A. 2006;103:13162–7.
Shaaban AM, Green AR, Karthik S, Alizadeh Y, Hughes TA, Harkins L, et al. Nuclear and cytoplasmic expression of ERbeta1, ERbeta2, and ERbeta5 identifies distinct prognostic outcome for breast cancer patients. Clin Cancer Res. 2008;14:5228–35.
Campbell-Thompson M, Lynch IJ, Bhardwaj B. Expression of estrogen receptor (ER) subtypes and ERbeta isoforms in colon cancer. Cancer Res. 2001;61:632–40.
Fritz A, Percy C, Jack A, Shanmugaratnam K, Sobin L, Parkin DM, et al. International classification of diseases for oncology. 3rd ed. Geneva: world health. Organization. 2000;
Edge S, Byrd DR, Compton CC, Fritz AG, Greene FL, Trotti A. American joint committee on cancer: AJCC cancer staging manual. 7th ed. New York, NY: Springer; 2016.
Yang IS, Son H, Kim S, Kim S. ISOexpresso: a web-based platform for isoform-level expression analysis in human cancer. BMC Genomics. 2016;17(1):631.
Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29:e45.
Hothorn T, Lausen B. Maximally selected rank statistics in R. R News. 2002;2:3–5.
Kennelly R, Kavanagh DO, Hogan AM, Winter DC. Oestrogen and the colon: potential mechanisms for cancer prevention. Lancet Oncol. 2008;9:385–91.
Hartman J, Gustafsson JA. Estrogen receptors in colorectal cancer: goalkeepers, strikers, or bystanders? Cancer Prev Res. 2010;3:897–9.
Waliszewski P, Blaszczyk M, Wolinska-Witort E, Drews M, Snochowski M, Hurst RE. Molecular study of sex steroid receptor gene expression in human colon and in colorectal carcinomas. J Surg Oncol. 1997;64:3–11.
Kaklamanos IG, Bathe OF, Franceschi D, Lazaris AC, Davaris P, Glinatsis M, et al. Expression of receptors for estrogen and progesterone in malignant colonic mucosa as a prognostic factor for patient survival. J Surg Oncol. 1999;72:225–9.
Oshima CT, Wonraht DR, Catarino RM, Mattos D, Forones NM. Estrogen and progesterone receptors in gastric and colorectal cancer. Hepato-Gastroenterology. 1999;46:3155–8.
Wong NA, Malcomson RD, Jodrell DI, Groome NP, Harrison DJ, Saunders PT. ERbeta isoform expression in colorectal carcinoma: an in vivo and in vitro study of clinicopathological and molecular correlates. J Pathol. 2005;207:53–60.
Thomas C, Gustafsson JA. The different roles of ER subtypes in cancer biology and therapy. Nat Rev Cancer. 2011;11:597–608.
Ascenzi P, Bocedi A, Marino M. Structure-function relationship of estrogen receptor alpha and beta: impact on human health. Mol Asp Med. 2006;27:299–402.
Colditz GA. Relationship between estrogen levels, use of hormone replacement therapy, and breast cancer. J Natl Cancer Inst. 1998;90:814–23.
Hankinson SE, Colditz GA, Willett WC. Towards an integrated model for breast cancer etiology: the lifelong interplay of genes, lifestyle, and hormones. Breast Cancer Res. 2004;6:213–8.
Ellem SJ, Risbridger GP. Treating prostate cancer: a rationale for targeting local oestrogens. Nat Rev Cancer. 2007;7:621–7.
Musgrove EA, Sutherland RL. Biological determinants of endocrine resistance in breast cancer. Nat Rev Cancer. 2009;9:631–43.
Speirs V, Green AR, Hughes TA, Ellis IO, Saunders PT, Shaaban AM. Clinical importance of estrogen receptor beta isoforms in breast cancer. J Clin Oncol. 2008;26:5825–6.
Stabile LP, Dacic S, Land SR, Lenzner DE, Dhir R, Acquafondata M, et al. Combined analysis of estrogen receptor beta-1 and progesterone receptor expression identifies lung cancer patients with poor outcome. Clin Cancer Res. 2011;17:154–64.
Shang Y. Molecular mechanisms of oestrogen and SERMs in endometrial carcinogenesis. Nat Rev Cancer. 2006;6:360–8.
Nguyen-Vu T, Wang J, Mesmar F, Mukhopadhyay S, Saxena A, McCollum CW, Gustafsson JÅ, Bondesson M, Williams C. Estrogen receptor beta reduces colon cancer metastasis through a novel miR-205 - PROX1 mechanism. Oncotarget. 2016 Jul 5;7(27):42159–71.
Byron SA, Van Keuren-Jensen KR, Engelthaler DM, Carpten JD, Craig DW. Translating RNA sequencing into clinical diagnostics: opportunities and challenges. Nat Rev Genet 2016 May;17(5):257–271. doi: 10.1038/nrg.2016.10. Epub 2016 Mar 21. Review. PubMed PMID: 26996076.
Costa C, Giménez-Capitán A, Karachaliou N, Rosell R. Comprehensive molecular screening: from the RT-PCR to the RNA-seq. Transl Lung Cancer Res. 2013 Apr;2(2):87–91.
Rajkumar AP, Qvist P, Lazarus R, Lescai F, Ju J, Nyegaard M, Mors O, Børglum AD, Li Q, Christensen JH. Experimental validation of methods for differential gene expression analysis and sample pooling in RNA-seq. BMC Genomics. 2015 Jul 25;16:548.
Giardiello FM, Hylind LM, Trimbath JD, Hamilton SR, Romans KE, Cruz-Correa M, et al. Oral contraceptives and polyp regression in familial adenomatous polyposis. Gastroenterology. 2005;128:1077–80.
Katzenellenbogen BS, Katzenellenbogen JA. Defining the “S” in SERMs. Science. 2002;295:2380–1.
Kuiper GG, Lemmen JG, Carlsson B, Corton JC, Safe SH, van der Saag PT, et al. Interaction of estrogenic chemicals and phytoestrogens with estrogen receptor beta. Endocrinology. 1998;139:4252–63.
Ettinger B, Black DM, Mitlak BH, Knickerbocker RK, Nickelsen T, Genant HK, et al. Reduction of vertebral fracture risk in postmenopausal women with osteoporosis treated with raloxifene: results from a 3-year randomized clinical trial. Multiple outcomes of Raloxifene evaluation (MORE) investigators. JAMA. 1999;282:637–45.
Cummings SR, Eckert S, Krueger KA, Grady D, Powles TJ, Cauley JA, et al. The effect of raloxifene on risk of breast cancer in postmenopausal women: results from the MORE randomized trial. Multiple outcomes of Raloxifene evaluation. JAMA. 1999;281:2189–97.
Vogel VG, Costantino JP, Wickerham DL, Cronin WM, Cecchini RS, Atkins JN, et al. Effects of tamoxifen vs raloxifene on the risk of developing invasive breast cancer and other disease outcomes: the NSABP study of tamoxifen and Raloxifene (STAR) P-2 trial. JAMA. 2006;295:2727–41.
Lechner D, Kallay E, Cross HS. Phytoestrogens and colorectal cancer prevention. Vitam Horm. 2005;70:169–98.
Kuiper GG, Carlsson B, Grandien K, Enmark E, Haggblad J, Nilsson S, et al. Comparison of the ligand binding specificity and transcript tissue distribution of estrogen receptors alpha and beta. Endocrinology. 1997;138:863–70.
Wang H, Li Q, Chen H. Genistein affects histone modifications on Dickkopf-related protein 1 (DKK1) gene in SW480 human colon cancer cell line. PLoS One. 2012;7:e40955.
Prossnitz ER, Barton M. The G-protein-coupled estrogen receptor GPER in health and disease. Nat Rev Endocrinol. 2011;7:715–26.
Liu Q, Chen Z, Jiang G, et al. Epigenetic down regulation of G protein-coupled estrogen receptor (GPER) functions as a tumor suppressor in colorectal cancer. Mol Cancer. 2017;16:87.
The authors would like to thank the AC Camargo Cancer Center Biobank for providing and processing the samples. We are grateful to Vinicius Calsavara for assistance with the statistical analysis.
This study was supported by grants from the National Council for Scientific and Technological Development (CNPq 301,603/2012–0; 573,589/08–9) and the National Institute of Science and Technology in Oncogenomics (INCITO - Fundação de Amparo à Pesquisa do Estado de São Paulo - FAPESP 2008/57887–9). The funders had no role in the study design, data collection, data analysis, decision to publish or preparation of the manuscript. SRR received an investigator fellowship award from CNPq.
Availability of data and materials
Additional data and materials may be obtained from the corresponding author on reasonable request.
Ethics approval and consent to participate
The study was conducted with the approval of the Institutional Review Board of the A.C. Camargo Cancer Center, Sao Paulo, Brazil (01453/10). All participants provided written informed consent prior to specimen collection.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Geometric mean (GM) of the Cqs values of the reference transcripts and the quantification/quality information of the RNA samples. Samples with a mean Cq greater than the mean + 1 standard deviation of the GM of the Cqs of the three reference genes were excluded (XLSX 40 kb)
ESR2 expression levels in normal mucosa and polyps of FAP patients and in sporadic colon carcinomas. All groups of samples were compared without consider the sample type (Fresh frozen tissue or FFPE). NS = Sporadic normal mucosa; SP = Sporadic polyps; ST = Sporadic tumour; NFAP = Normal mucosa FAP; PFAP = Polyp FAP; *: p < 0.05; **: p < 0.01; ***: p < 0.001; (Tukey’s Multiple Comparison Test and Student’s t-test) (PNG 86 kb)
ESR2 isoforms expression levels in normal and in adenocarcinoma colon tissues deposited in the TCGA database. The isoforms most prevalent identified as uc001xgy.2 (green), uc001xgu.3 (blue), uc001xgx.3 (light blue), uc001xgw.3 (pink), and uc001xgz.2 (purple) coding (pink), and No-coding (purple) represents the transcripts variants and isoforms d (ERβ4), b (ERβ2), j (Non-coding RNA – (NC)), h (NC) and e (NC), respectively (PNG 106 kb)
Kaplan − Meier estimate and cumulative incidence curve of the disease-free survival of colorectal cancer patients as a function of ERβ1 and ERβ5 isoform expression levels (primer pair 2: transcript variants a and g). The cutoff for relative gene expression of −4.5 (log2) was determined using the log-rank test (maximally selected rank statistics in R)  (PNG 86 kb)
About this article
Cite this article
Stevanato Filho, P.R., Aguiar Júnior, S., Begnami, M.D. et al. Oestrogen receptor beta isoform expression in sporadic colorectal cancer, familial adenomatous polyposis and progressive stages of colorectal cancer. BMC Cancer 17, 754 (2017). https://doi.org/10.1186/s12885-017-3688-4