Human C-terminally truncated ERα variants resulting from the use of alternative exons in the ligand-binding domain

Estrogen Receptor α Isoforms Generated by Alternative Use of Cryptic Exons

Estrogen receptor α (ERα) regulates several physiological functions. In pathophysiological conditions, ERα is involved in the development and progression of estrogen-sensitive tumors. The ERα gene contains multiple 5'-untranslated exons and eight conventional coding exons and presents multiple isoforms generated by alternative promoter usage and alternative splicing. This gene also possesses non-conventional exons in the 3'- and intronic regions, and alternative use of cryptic exons contributes to further diversity of ERα mRNAs and proteins. Recently, the genomic organization of ERα genes and the splicing profiles of their transcripts were comparatively analyzed in humans, mice, and rats, and multiple ERα isoforms with distinct structures and functions were identified. These transcripts contain cryptic sequences that encode insertion-containing or truncated ERα proteins. In particular, alternative cryptic exons with in-frame stop codons yield transcripts encoding C-terminally-truncated ERα proteins. The C-terminally-truncated ERα isoforms lack part or all of the ligand-binding domain but possess an isoform-specific sequence. Some of these isoforms exhibit constitutive transactivation and resistance to estrogen receptor antagonists. Although numerous studies have reported conflicting results regarding their functions, the critical determinant for their gain-of-function has been identified structurally. Here we review recent progress in ERα variant research concerning the genomic organization of ERα genes, splicing profiles of ERα transcripts, and transactivation properties of ERα isoforms.

Alternative promoters and splicing create multiple functionally distinct isoforms of oestrogen receptor alpha in breast cancer and healthy tissues

BACKGROUND

Oestrogen receptor alpha (ER) is involved in cell growth and proliferation and functions as a transcription factor, a transcriptional coregulator, and in cytoplasmic signalling. It affects, for example, bone, endometrium, ovaries and mammary epithelium. It is a key biomarker in clinical management of breast cancer, where it is used as a prognostic and treatment-predictive factor, and a therapeutical target. Several ER isoforms have been described, but transcript annotation in public databases is incomplete and inconsistent, and functional differences are not well understood.

METHODS

We have analysed short- and long-read RNA sequencing data from breast tumours, breast cancer cell lines, and normal tissues to create a comprehensive annotation of ER transcripts and combined it with experimental studies of full-length protein and six alternative isoforms.

RESULTS

The isoforms have varying transcription factor activity, subcellular localisation, and response to the ER-targeting drugs tamoxifen and fulvestrant. Antibodies differ in ability to detect alternative isoforms, which raises concerns for the interpretation of ER-status in routine pathology.

CONCLUSIONS

Future work should investigate the effects of alternative isoforms on patient survival and therapy response. An accurate annotation of ER isoforms will aid in interpretation of clinical data and inform functional studies to improve our understanding of the ER in health and disease.

RSK-Mediated Non-canonical Activation of EphA2 by Tamoxifen

The long-term administration of tamoxifen to estrogen receptor α (ERα)-positive breast cancer patients is an established treatment that reduces mortality and recurrence. However, resistance to tamoxifen and an increased risk of endometrial cancer may occur; therefore, the mechanisms by which tamoxifen causes these adverse effects warrant further study. Tamoxifen has been shown to activate mitogen-activated protein kinase (MAPK) in an ERα-independent manner; therefore, we investigated its effects on the MAPK-mediated non-canonical activation of EphA2, a critical event regulating cell migration. Tamoxifen at slightly higher concentrations induced the rapid phosphorylation of EphA2 at Ser-897 via the MAPK/extracellular signal-regulated kinase (ERK) kinase (MEK)-ERK-ribosomal S6 kinases (RSK) pathway in HeLa cells. In addition, tamoxifen significantly enhanced the migration ability of ERα-negative MDA-MB-231 breast cancer cells in RSK- and EphA2-dependent manners. Phosphorylated EphA2 was internalized and re-localized to the plasma membrane, including lamellipodia, in an RSK-dependent manner. Collectively, the present results provide novel insights into the tumor-promoting activity of tamoxifen.

ESR1 ChIP-Seq Identifies Distinct Ligand-Free ESR1 Genomic Binding Sites in Human Hepatocytes and Liver Tissue

The estrogen receptor alpha (ESR1) is an important gene transcriptional regulator, known to mediate the effects of estrogen. Canonically, ESR1 is activated by its ligand estrogen. However, the role of unliganded ESR1 in transcriptional regulation has been gaining attention. We have recently shown that ligand-free ESR1 is a key regulator of several cytochrome P450 (CYP) genes in the liver, however ligand-free ESR1 has not been characterized genome-wide in the human liver. To address this, ESR1 ChIP-Seq was conducted in human liver samples and in hepatocytes with or without 17beta-estradiol (E2) treatment. We identified both ligand-dependent and ligand-independent binding sites throughout the genome. These two ESR1 binding categories showed different genomic localization, pathway enrichment, and cofactor colocalization, indicating different ESR1 regulatory function depending on ligand availability. By analyzing existing ESR1 data from additional human cell lines, we uncovered a potential ligand-independent ESR1 activity, namely its co-enrichment with the zinc finger protein 143 (ZNF143). Furthermore, we identified ESR1 binding sites near many gene loci related to drug therapy, including the CYPs. Overall, this study shows distinct ligand-free and ligand-bound ESR1 chromatin binding profiles in the liver and suggests the potential broad influence of ESR1 in drug metabolism and drug therapy.

Quantitative expression data of human estrogen receptor α variants in non-functioning pituitary adenomas obtained by reverse transcription-digital polymerase chain reaction analysis

Expression profiles of gonadal steroid receptor variants have been reportedly associated with malignancy in breast and prostate cancers [1,2]. However, such associations with pituitary tumors remain unclear. Therefore, the expression levels of the wild-type ESR1 (ERα66) and the ESR1 variants (ERαi34, ERαi45c, and ERαΔ5) transcripts encoding constitutively active ERα proteins with C-terminal truncation in non-functioning pituitary adenomas (NFPAs) were evaluated using reverse transcription-digital polymerase chain reaction. The results revealed that the expression levels of the variants were approximately two orders of magnitude lower than that of ERα66 in NFPAs. These data were based on our previous article entitled “Accurate assessment of estrogen receptor profiles in non-functioning pituitary adenomas using RT-digital PCR and immunohistochemistry” [3].

Identification of a novel C-terminally truncated estrogen receptor α variant (ERαi34) with constitutive transactivation and estrogen receptor antagonist resistance

Constitutively active estrogen receptor α (ERα) variants with C-terminal truncation are candidate molecules for gain of both endocrine- and chemotherapy-resistance in estrogen-sensitive tumors. Our previous reports using artificially truncated ERα constructs demonstrated that ERα variants encoded in 1-2-3-cryptic exon- and 1-2-3-4-cryptic exon-types of transcripts have potentials to display constitutive transactivation of an estrogen response element reporter. However, naturally occurring 1-2-3-cryptic exon-type ERα variants have not been cloned yet. Therefore, the present study was designed to identify naturally occurring ERα variants encoded in 1-2-3-cryptic exon-type ERα transcripts. We cloned a novel C-terminally truncated ERα variant (ERαi34) encoded in a 1-2-3-i34 transcript from MCF-7 cells and characterized its constitutive and ER antagonist-resistant transactivation in transfected cells. Stable transfection of the variant into MCF-7 cells augmented basal cell proliferation insensitive to fulvestrant. Collectively, we validated the structure-based mechanisms underlying constitutive and ER antagonist-resistant transactivation by C-terminally truncated ERα variants.

Applicability of Anti-Human Estrogen Receptor β Antibody PPZ0506 for the Immunodetection of Rodent Estrogen Receptor β Proteins

Several lines of controversial evidence concerning estrogen receptor β (ERβ) remain to be solved because of the unavailability of specific antibodies against ERβ. The recent validation analysis identified a monoclonal antibody (PPZ0506) with sufficient specificity against human ERβ. However, the specificity and cross-reactivity of PPZ0506 antibody against ERβ proteins from laboratory animals have not been confirmed. In the present study, we aimed to validate the applicability of PPZ0506 to rodent studies. The antibody exhibited specific cross-reactivity against mouse and rat ERβ proteins in immunoblot and immunocytochemical experiments using transfected cells. In immunohistochemistry for rat tissue sections, PPZ0506 showed immunoreactive signals in the ovary, prostate, and brain. These immunohistochemical profiles of rat ERβ proteins in rat tissues accord well with its mRNA expression patterns. Although the antibody was reported to show the moderate signals in human testis, no immunoreactive signals were observed in rat testis. Subsequent RT-PCR analysis revealed that this species difference in ERβ expression resulted from different expression profiles related to the alternative promoter usage between humans and rats. In conclusion, we confirmed applicability of PPZ0506 for rodent ERβ studies, and our results provide a fundamental basis for further examination of ERβ functions.

Ligand-Free Estrogen Receptor α (ESR1) as Master Regulator for the Expression of CYP3A4 and Other Cytochrome P450 Enzymes in the Human Liver

Cytochrome P450 3A4 isoform (CYP3A4) transcription is controlled by hepatic transcription factors (TFs), but how TFs dynamically interact remains uncertain. We hypothesize that several TFs form a regulatory network with nonlinear, dynamic, and hierarchical interactions. To resolve complex interactions, we have applied a computational approach for estimating Sobol's sensitivity indices (SSI) under generalized linear models to existing liver RNA expression microarray data (GSE9588) and RNA-seq data from genotype-tissue expression (GTEx), generating robust importance ranking of TF effects and interactions. The SSI-based analysis identified TFs and interacting TF pairs, triplets, and quadruplets involved in CYP3A4 expression. In addition to known CYP3A4 TFs, estrogen receptor α (ESR1) emerges as key TF with the strongest main effect and as the most frequently included TF interacting partner. Model predictions were validated using small interfering RNA (siRNA)/short hairpin RNA (shRNA) gene knockdown and clustered regularly interspaced short palindromic repeats (CRISPR)-mediated transcriptional activation of ESR1 in biliary epithelial Huh7 cells and human hepatocytes in the absence of estrogen. Moreover, ESR1 and known CYP3A4 TFs mutually regulate each other. Detectable in both male and female hepatocytes without added estrogen, the results demonstrate a role for unliganded ESR1 in CYP3A4 expression consistent with unliganded ESR1 signaling reported in other cell types. Added estrogen further enhances ESR1 effects. We propose a hierarchical regulatory network for CYP3A4 expression directed by ESR1 through self-regulation, cross regulation, and TF-TF interactions. We also demonstrate that ESR1 regulates the expression of other P450 enzymes, suggesting broad influence of ESR1 on xenobiotics metabolism in human liver. Further studies are required to understand the mechanisms underlying role of ESR1 in P450 regulation. SIGNIFICANCE STATEMENT: This study focuses on identifying key transcription factors and regulatory networks for CYP3A4, the main drug metabolizing enzymes in liver. We applied a new computational approach (Sobol's sensitivity analysis) to existing hepatic gene expression data to determine the role of transcription factors in regulating CYP3A4 expression, and used molecular genetics methods (siRNA/shRNA gene knockdown and CRISPR-mediated transcriptional activation) to test these interactions in life cells. This approach reveals a robust network of TFs, including their putative interactions and the relative impact of each interaction. We find that ESR1 serves as a key transcription factor function in regulating CYP3A4, and it appears to be acting at least in part in a ligand-free fashion.

WWOX Phosphorylation, Signaling, and Role in Neurodegeneration

Homozygous null mutation of tumor suppressor WWOX/Wwox gene leads to severe neural diseases, metabolic disorders and early death in the newborns of humans, mice and rats. WWOX is frequently downregulated in the hippocampi of patients with Alzheimer’s disease (AD). In vitro analysis revealed that knockdown of WWOX protein in neuroblastoma cells results in aggregation of TRAPPC6AΔ, TIAF1, amyloid β, and Tau in a sequential manner. Indeed, TRAPPC6AΔ and TIAF1, but not tau and amyloid β, aggregates are present in the brains of healthy mid-aged individuals. It is reasonable to assume that very slow activation of a protein aggregation cascade starts sequentially with TRAPPC6AΔ and TIAF1 aggregation at mid-ages, then caspase activation and APP de-phosphorylation and degradation, and final accumulation of amyloid β and Tau aggregates in the brains at greater than 70 years old. WWOX binds Tau-hyperphosphorylating enzymes (e.g., GSK-3β) and blocks their functions, thereby supporting neuronal survival and differentiation. As a neuronal protective hormone, 17β-estradiol (E2) binds WWOX at an NSYK motif in the C-terminal SDR (short-chain alcohol dehydrogenase/reductase) domain. In this review, we discuss how WWOX and E2 block protein aggregation during neurodegeneration, and how a 31-amino-acid zinc finger-like Zfra peptide restores memory loss in mice.

17β-estradiol regulates the RNA-binding protein Nova1, which then regulates the alternative splicing of estrogen receptor β in the aging female rat brain

Alternative RNA splicing results in the translation of diverse protein products arising from a common nucleotide sequence. These alternative protein products are often functional and can have widely divergent actions from the canonical protein. Studies in humans and other vertebrate animals have demonstrated that alternative splicing events increase with advanced age, sometimes resulting in pathological consequences. Menopause represents a critical transition for women, where the beneficial effects of estrogens are no longer evident; therefore, factors underlying increased pathological conditions in women are confounded by the dual factors of aging and declining estrogens. Estrogen receptors (ERs) are subject to alternative splicing, the spliced variants increase following menopause, and they fail to efficiently activate estrogen-dependent signaling pathways. However, the factors that regulate the alternative splicing of ERs remain unknown. We demonstrate novel evidence supporting a potential biological feedback loop where 17β-estradiol regulates the RNA-binding protein Nova1, which, in turn, regulates the alternative splicing of ERβ. These data increase our understanding of ER alternative splicing and could have potential implications for women taking hormone replacement therapy after menopause.

Characterization of rodent constitutively active estrogen receptor α variants and their constitutive transactivation mechanisms

Estrogen receptor α (ERα) mRNAs exhibit remarkable heterogeneity owing to complicated alternative splicing. Some encode C-terminally-truncated ERα proteins, which display ligand-independent transactivation or dominant-negative activity. We previously characterized C-terminally-truncated ERα mRNA variants with cryptic sequences in humans and mice, and demonstrated that helices in the ligand-binding domains (LBDs) of ERα variants contribute to ligand-independent transcriptional activity. However, existence of non-conventional coding exons and generation of constitutively active ERα variants have remained to be examined in rats. To comparatively analyze modular organization and splicing profiles of the ERα genes, the range of C-terminally-truncated ERα variants was explored in rats and mice using rapid amplification of cDNA ends and RT-PCR cloning. Furthermore, their functions were determined in transiently transfected cells using expression constructs and luciferase reporter assays. Multiple cryptic exons and C-terminally-truncated ERα variant mRNAs were identified in rats and mice. Naturally occurring and artificially truncated variants/constructs lacking helix 5 potentially exhibited gain-of-function in transfected cells. Although cryptic exons and splicing profiles were poorly conserved among humans, mice, and rats, constitutively active variants were generated from the ERα genes. Moreover, the primary mechanism of ligand-independent activation by C-terminally-truncated ERα variants is C-terminus to helix 5 deletion in the LBD. This comparative study documented the complexity of ERα genes, mRNAs, and proteins, and further determined the underlying structural basis of ligand-independent activation by C-terminally-truncated ERα variants.

The Importance of Steroid Uptake and Intracrine Action in Endometrial and Ovarian Cancers

Endometrial and ovarian cancers predominately affect women after menopause, and are more frequently observed in developed countries. These are considered to be hormone-dependent cancers, as steroid hormones, and estrogens in particular, have roles in their onset and progression. After the production of estrogens in the ovary has ceased, estrogen synthesis occurs in peripheral tissues. This depends on the cellular uptake of estrone-sulfate and dehydroepiandrosterone-sulfate, as the most important steroid precursors in the plasma of postmenopausal women. The uptake through transporter proteins, such as those of the organic anion-transporting polypeptide (OATP) and organic anion-transporter (OAT) families, is followed by the synthesis and action of estradiol E2. Here, we provide an overview of the current understanding of this intracrine action of steroid hormones, which depends on the availability of the steroid precursors and transmembrane transporters for precursor uptake, along with the enzymes for the synthesis of E2. The data is also provided relating to the selected transmembrane transporters from the OATP, OAT, SLC51, and ABC-transporter families, and the enzymes involved in the E2-generating pathways in cancers of the endometrium and ovary. Finally, we discuss these transporters and enzymes as potential drug targets.

Intronic SNP in ESR1 encoding human estrogen receptor alpha is associated with brain ESR1 mRNA isoform expression and behavioral traits

Genetic variants of ESR1 have been implicated in multiple diseases, including behavioral disorders, but causative variants remain uncertain. We have searched for regulatory variants affecting ESR1 expression in human brain, measuring allelic ESR1 mRNA expression in human brain tissues with marker SNPs in exon4 representing ESR1-008 (or ESRα-36), and in the 3'UTR of ESR1-203, two main ESR1 isoforms in brain. In prefrontal cortex from subjects with bipolar disorder, schizophrenia, and controls (n = 35 each; Stanley Foundation brain bank), allelic ESR1 mRNA ratios deviated from unity up to tenfold at the exon4 marker SNP, with large allelic ratios observed primarily in bipolar and schizophrenic subjects. SNP scanning and targeted sequencing identified rs2144025, associated with large allelic mRNA ratios (p = 1.6E10-6). Moreover, rs2144025 was significantly associated with ESR1 mRNA levels in the Brain eQTL Almanac and in brain regions in the Genotype-Tissue Expression project. In four GWAS cohorts, rs2104425 was significantly associated with behavioral traits, including: hypomanic episodes in female bipolar disorder subjects (GAIN bipolar disorder study; p = 0.0004), comorbid psychological symptoms in both males and females with attention deficit hyperactivity disorder (GAIN ADHD, p = 0.00002), psychological diagnoses in female children (eMERGE study of childhood health, subject age ≥9, p = 0.0009), and traits in schizophrenia (e.g., grandiose delusions, GAIN schizophrenia, p = 0.0004). The first common ESR1 variant (MAF 12-33% across races) linked to regulatory functions, rs2144025 appears conditionally to affect ESR1 mRNA expression in the brain and modulate traits in behavioral disorders.