Early growth response protein 1 regulates promoter activity of α-plasma membrane calcium ATPase 2, a major calcium pump in the brain and auditory system
BMC Molecular Biology volume 18, Article number: 14 (2017)
Along with sodium/calcium (Ca2+) exchangers, plasma membrane Ca2+ ATPases (ATP2Bs) are main regulators of intracellular Ca2+ levels. There are four ATP2B paralogs encoded by four different genes. Atp2b2 encodes the protein pump with the fastest activation, ATP2B2. In mice, the Atp2b2 transcript has several alternate transcriptional start site variants: α, β, µ and δ. These variants are expressed in developmental and tissue specific manners. The α and β Atp2b2 transcripts are equally expressed in the brain. αAtp2b2 is the only transcript found in the outer hair cells of young mice (Silverstein RS, Tempel BL. in Neuroscience 141:245–257, 2006). Mutations in the coding region of the mouse Atp2b2 gene indicate a narrow window for tolerated dysfunction of the ATP2B2 protein, specifically in the auditory system. This highlights the necessity of tight regulation of this gene for normal cell physiology.
Although ATP2Bs are important regulators of Ca2+ in many cell types, little is known about their transcriptional regulation. This study identifies the proximal promoter of the αAtp2b2 transcript. Further investigations indicate that ATOH1 and EGR1 modulate promoter activity. Additionally, we report that EGR1 increases endogenous expression of Atp2b2 transcript in two cell lines. Electrophoretic mobility shift assays (EMSA) indicate that EGR1 binds to a specific site in the CpG island of the αAtp2b2 promoter.
This study furthers our understanding of Atp2b2 regulation by: (I) elucidating transcriptional regulatory mechanisms for Atp2b2, and (II) identifying transcription factors that modulate expression of Atp2b2 in the brain and peripheral auditory system and (III) allows for future studies modulating gene expression of Atp2b2.
Importance of Ca2+ and Ca2+ regulators
Ca2+ is involved in cellular excitation and is an important second messenger. In a resting cell, typical intracellular Ca2+ concentrations are between 0.1 and 1.0 μM, while extracellular levels are much higher (~2 mM). Consequently, maintaining this large electrochemical gradient is critical for normal cell physiology and necessitates an energy dependent mechanism of Ca2+ expulsion . Plasma membrane Ca2+ ATPases are ubiquitously expressed in the plasma membrane and use ATP to pump Ca2+ out of the cell. Four paralogs are found in mammalian cells: Atp2b1, Atp2b2, Atp2b3 and Atp2b4. These proteins are highly similar but differ in tissue expression and speed of activation.
Animal and disease models associated with ATP2Bs
The importance of ATP2Bs for normal cellular function is demonstrated by widespread tissue expression, a pivotal role in disease, and the lethality of ATP2B null mice . In mammals, ATP2B1 and ATP2B4 are the housekeeping Ca2+ regulators expressed throughout all tissues . Null mutations of ATP2B1 are embryonic lethal and null mutations of ATP2B4 cause problems with sperm motility . Genome wide association studies in humans have linked SNPs in ATP2B1 to increased cardiovascular risk  and a mutation in ATP2B4 was associated with spastic paraplegia . ATP2B2 is expressed predominately in the brain and mammary glands. In mice, Atp2b2 null mutations cause profound deafness and ataxia [6, 7]. Similarly, mutations in human ATP2B2 are associated with hearing loss [8, 9]. Atp2b2 is also an emerging player in autism and breast cancer [10, 11].
ATP2B2 is highly expressed in the auditory system which is extremely sensitive to small physiological changes. This makes it an ideal model for studying fine tuning of Ca2+ regulation. The deafwaddler (dfw) mutant mice, dfw 2J and dfw i5, have null mutations in Atp2b2 resulting in negligible expression of Atp2b2 transcript and protein [6, 12]. Heterozygous mutants express about 50% normal protein levels and have significantly impaired hearing sensitivity compared to wild-type animals . This finding suggests that one normal copy of Atp2b2 is insufficient to produce a normal hearing phenotype and defines Atp2b2 as haploinsufficient.
The original dfw mutation is a hypomorph that decreases, but does not completely abolish ATP2B2 activity. Studies of pump kinetics estimate that total function of ATP2B2 in homozygous dfw mutants is ~30% of normal animals . Although homozygous mutants are profoundly deaf, heterozygous dfw mice do NOT have any significant hearing loss . This is surprising considering the predicted increase in total pump function for dfw heterozygotes over the dfw 2J and dfw i5 heterozygotes is quite small. Thus, the threshold for ATP2B2 pump dysfunction in the auditory system is narrow and exhibits the necessity for tight regulation of this gene and protein.
The importance of studying endogenous regulation of ATP2Bs
Although a wealth of information exists for post-transcriptional mechanisms regulating ATP2Bs, there is very little knowledge of transcriptional initiation of these genes. The only ATP2B promoter studied to date was reported in a brief publication outlining the promoter elements of Atp2b1 . There have been no studies looking at the promoter of the Atp2b2 gene. It is clear that understanding the intricacies of ATP2B gene transcription will be important for combating the pathophysiology of diseases associated with this family of Ca2+ regulators. With over 30 different splice variants and four different genes, the ATP2B family is vast. The current study aims to characterize the promoter elements of the neuronal and hair cell-specific α transcript of Atp2b2, and identify important transcription factors that bind to and modulate expression of this gene.
OC-1 and OC-2 cells  were grown at 33 °C and 5% CO2 in Modified Eagle Media with Glutamax (Gibco: 41090-036) supplemented with 10% fetal bovine serum and 50 U/ml of γ-Interferon. Neuro-2a (N2A) cells were grown at 37 °C and 5% CO2 in Dulbecco’s Modified Eagle Medium (Gibco) supplemented with 10% fetal bovine serum and 1% streptomycin and penicillin. HeLa cells were grown at 37 °C and 5% CO2 in Dulbecco’s Modified Eagle Medium (Gibco) supplemented with 10% fetal bovine serum and 1% penicillin/streptomycin.
Cloning promoter constructs
PCR was used to amplify the ~5.5 kbp proximal promoter in three pieces off of genomic DNA from the DBA/2J inbred mouse strain. Next, exons of the 5′ untranslated region (UTR) were amplified from complimentary DNA (cDNA) of dfw 2J brainstem. A synthetic NcoI cut site was inserted in the 5′ UTR at the translational start site. These four pieces were cloned into TOPO® vectors (ThermoFisher Scientific) and sequenced. The four pieces were digested out of the TOPO® vectors with the following restriction enzymes: I—HindIII and AvrII, II—AvrII and BssSI, III—BssSI and Hpy99I, IV—Hpy99I and NcoI (Table 1). These pieces were ligated into the pGL3 Firefly Luciferase Vector (Promega Corporation: E1751).
Cloning transcription factor constructs
Transcription factor cDNAs were subcloned from Addgene plasmids (Addgene: ATOH1: 33333, EGR1: 11729 and GATA3: 1332) or amplified from cochlear or brainstem tissue of CBA/CaJ mice (USF1 and POU4F3). All cDNAs were ligated into the pcDNA vector with an internal ribosome entry site (IRES) followed by EGFP. The EGR1 cDNA contains an N-terminal Flag® tag. EGFP expression was noted in cells transiently transfected with the transcription factor constructs. Transcription factor cDNA was quantified in the OC-1 and N2A cell lines at baseline and after transfection with the plasmids. All of the transcription factors were expressed at similar levels after transfection (Additional file 1: Table S1).
OC-1 and OC-2 cells were plated 24 h before all transfections at a density of 4 × 104 cells per well in a 24 well plate. Cells were harvested 48 h after transfection with passive lysis buffer and assayed using the dual luciferase reporter assay system (Promega: E1910) with a manual luminometer (Promega). Three biological replicates were assayed in each experiment. Data is the average of at least three experiments.
For promoter element identification: the αAtp2b2 promoter constructs in the pGL3 vector were transfected at 400 ng/well into OC-1 or OC-2 mammalian cells. The renilla luciferase standard vector (Promega) was included at 5 ng/well and a DNA:transfection reagent ratio of 1:4 was used (Fugene HD, Promega: E2311).
For luciferase and transcription factor co-transfection assays: The αAtp2b2 promoter construct [+572/−2133] was transfected (200 ng/well) with each pcDNA transcription factor vector: Atoh1, Egr1, Gata3, Pou4f3, OR Usf1 (200 ng/well). For Egr1 only, the pcDNA Egr1 transcription factor vector was transfected at 200 ng/well with three of the αAtp2b2 promoter construct truncations [+572/−287], [+572/−855], OR [+572/−2133] (200 ng/well). For all co-transfection assays: the renilla luciferase standard vector was included at 2.5 ng/well and a DNA:transfection reagent ration of 1:4 was used.
OC-1 and N2A cells were plated 24 h before transfection in 12 well plates. OC-1 cells were plated at a density of 8 × 104 cells per well and N2A cells were plated at a density of 1 × 105 cells per well. Transcription factor constructs were transfected at 800 ng/well using a 1:3 ratio of DNA:transfection reagent. Cells were incubated with constructs for 48 h and then harvested in 400 µl of QiAzol™ lysis reagent. mRNA was extracted using RNeasy Plus Universal Mini Kit (Qiagen: 73404) and subsequently converted into cDNA using random hexamer (ThermoFisher Scientific: N8080127) and reverse transcriptase (Clontech SMARTscribe: 639537).
Baseline expression experiments utilized primers designed to recognize specific cDNAs (Table 2). IQ™ SYBR® Green Supermix (Bio-Rad: 170-8886) reagent was used. Endogenous levels of Atp2b2 and Atp2b4 transcript were quantified using premade gene expression assays from Applied Biosystems. (Atp2b2: Mm00437640_m1 and Atp2b4: Mm01285597_m1). TaqMan reactions were run using the SSOAdvanced™ Universal Probes Supermix (Bio-Rad: 172-5281). Expression levels were normalized to housekeeping genes, hydroxymethylbilane synthase (Hmbs) and γ-Actin (Actg1) (Applied Biosystems: Hmbs: Mm01143545_m1 and Actg1: Mm01963702_s1). All experiments were done on a Bio-Rad Q5 iCycler® or CFX96 Touch™ (Bio-Rad). For all experiments, three biological replicates were assayed, data is the average of at least three runs where replicates within a run had a relative standard deviation of <0.03 (replicates with standard deviation/cycle threshold >0.03 were excluded).
HeLa cells were grown to confluency in two T-75 cm2 flasks and extracted using trypsin. Cells were thoroughly washed and nuclear extracts were isolated using cell extraction buffer (ThermoFisher Scientific: FNN0011) supplemented with protease inhibitor (Sigma–Aldrich). Protein concentration was determined using a bicinchoninic acid assay (ThermoFisher Scientific: 23227). Double stranded DNA probes were made by boiling complimentary single stranded probes at a concentration of 25 µm in annealing buffer for 5 min in 1 l of boiling water. The boiling water was removed from the heat and the tubes were allowed to slow cool to room temperature. The canonical EGR1 positive control was designed based on the canonical consensus site for EGR1 (GCGGGGGCG) [17, 18]. The Santa Cruz EGR1 control probe was purchased from Santa Cruz Biotechnology (SC-2529). All other probes were developed using binding site predictions and consensus information from MotifMap, MatInspector and TFBIND. Biotin labeled probes were purchased from Integrated DNA Technologies (Table 3). EMSAs were performed using the LightShift™ Chemiluminescent EMSA protocol (ThermoFisher Scientific: 21048) with 18 µg of HeLa nuclear extract per 20 µl reaction. The reactions were incubated at room temperature for 30 min without the addition of the biotin labeled probe. Upon addition of the biotin labeled probes the reactions were incubated for another 30 min. For the supershift assays: 2 µg of antibody were added to the reaction and incubated for 30 min prior to the addition of biotin labeled probes. The EGR1 antibody was purchased from Abcam® (Ab174509, Rabbit polyclonal to N-terminus) and the negative control antibody to ATP2B2 was purchased from ThermoFisher Scientific (PA1-915, Rabbit polyclonal). The reactions were electrophoresed in 0.5% TBE on 5% TBE polyacrylamide gels. The gels were run at 4 °C for 2.5–3 h at 60 V. Gels were transferred at 350 mA for 1 h at 4 °C to Biodyne™ B Nylon Membrane (Pierce: 77016), cross-linked for 15 min with 312 nm bulbs. Gels were visualized using the Chemiluminescent Nucleic Acid Detection Module (ThermoFisher Scientific: 89880). Each experiment included the standard shift assay and the test assay. Test assays include: self-competition assay (n = 5), modified EGR1 probe assay (n = 2) and the supershift assay (n = 4). Quantification was performed on blots with exposures between 5 and 10 min where none of the bands were saturated. For each blot, background intensity was subtracted from shift bands. In the self-competition assays, competition probe concentrations ranged from 16× to 600× the standard probe. Self-competition experiments including more than one competition probe concentration were averaged and counted as a single experiment.
Genomic modifications at Atp2b2 indicate transcriptional hot-spots in mouse cerebellum
The genomic region surrounding Atp2b2 is vast and complex. With 22 exons and four different transcriptional start site (TSS) variants, there are many opportunities for gene regulation. Of the four distinct Atp2b2 TSS variants, α and β are the variants expressed in neurons. αAtp2b2 is the primary transcript in the hair cells of young mice (postnatal day 9) . The α and β transcripts share the same coding exons but differ in the sequence of their 5′ untranslated regions. Utilizing the UCSC genome browser, we compared the genomic regions of α and β Atp2b2 in two distinct tissues to identify transcriptionally active domains. Chromatin features of interest included DNase I hypersensitivity (DHS), an indication of uncoiled, accessible DNA; RNA polymerase II (RNA Pol II) occupancy, the enzyme that transcribes DNA into RNA; and trimethylation of histone H3 at lysine 4 (H3K4me3), a hallmark of transcriptionally active promoters. The heat maps for these marks in the region surrounding the 5′ end of Atp2b2 reveal several distinct signals in 8 week mouse cerebellum (indicative of active transcription) not detected in 8 week old heart  (Fig. 1). In particular, the regions directly surrounding the noncoding Iα and Iβ exons contain significant levels of H3K4me3 and RNA Pol II occupancy and coincide with predicted CpG islands, known initiation sites for ubiquitously expressed genes . Together, these data suggest that αAtp2b2 and βAtp2b2 transcripts are both transcribed in the cerebellum and have unique start sites. This allows for differential regulation and tissue specific expression. We have elected to focus on the αAtp2b2 transcript because of its unique expression profile in the brain and sensory epithelium.
Identifying the promoter of the αAtp2b2 transcript
The location of the CpG island and pattern of histone methylation in the cerebellum suggest that promoter elements are likely located directly upstream of the αAtp2b2 TSS. To test this empirically, the 5′ UTR and 5 kbp region containing the CpG island were cloned into a luciferase reporter vector (Fig. 2a). All constructs are numbered relative to the TSS, which maps to the negative strand of chromosome 6 at base 114,042,026 (Ensembl version 82, GRCm38.p4). The “full-length” αAtp2b2 promoter construct contains 5399 bases located directly upstream of the TSS from 114,042,027 [closest to TSS] to 114,047,365 [farthest from TSS]. Note that only the exons of the 5′ UTR were included in the construct, the 5′ UTR spans approximately 200 kbp of the genome (Table 4).
Ten predominately neuronal or sensory-epithelium derived cell lines were screened for expression of Atp2b2 transcripts (OC-K3, HEI-OCI, E-36, N2A, NIH/3T3, C2C12, OC-2, OC-1, ARPE-19, HeLa). Atp2b2 gene expression was measured via quantitative PCR (qPCR) with a cutoff at 40 cycles. Under these conditions only two cell lines expressed detectable levels of Atp2b2, OC-1 cells from the immortomouse organ of Corti and the N2A neuroblastoma cell line. Given the origins of the other cell lines it is surprising that Atp2b2 is not expressed; we speculate that this is an artifact of immortalization.
For the luciferase assay, we chose to use OC-1 cells and related OC-2 cells. These cell lines were immortalized during hair cell differentiation and, although OC-2 cells do not consistently express Atp2b2, OC-2 cells express hair cell markers that OC-1 lacks . For this reason it was included in our luciferase assays. A series of αAtp2b2 promoter luciferase constructs were transiently transfected into OC-1 and OC-2 cells and the fold change in promoter activity over empty pGL3 vector was determined (Fig. 2b). The greatest amount of luciferase activity above background was observed with the shortest promoter construct [+572/−287]. The inclusion of an additional 500 bp of 5′ sequence [+572/−855] had no detectable effect on promoter activity. When the length of the promoter fragment was extended further in the 5′ direction, the amount of luciferase activity began to decrease until it reached baseline levels. Reversing the orientation of the promoter or removal of the 5′ UTR and CpG islands completely abolished any detectable activity. These data demonstrate that the 5′ UTR and CpG islands are necessary and sufficient for promoter activity, and maps the minimal promoter between bases +572 and −287. The remaining proximal promoter region most likely contains important tissue specific enhancer and repressor elements.
ATOH1, EGR1, GATA3, POU4F3, and USF1 are potential transcriptional modulators of αAtp2b2
Binding prediction software, MatInspector and TFBIND [25, 26] were used to identify potential transcription factor binding sites in the proximal promoter of the αAtp2b2 transcript. The entire region (+572 to −5399 bps) was analyzed and contained putative sites for thousands of transcription factors. To narrow down likely candidates we compared factors from the αAtp2b2 promoter to the βAtp2b2 proximal promoter and identified factors that were enriched in or unique to the αAtp2b2 promoter. Using the Shared Harvard Inner Ear Laboratory Database (SHIELD) we determined expression of candidate genes in auditory hair cells  (Fig. 3a). We cross referenced these lists with transcription factor expression data for the OC-1 and OC-2 cells used in the luciferase assays. Published data for transcript expression in OC-1 cells was verified for select genes via qPCR (Fig. 3b). Two important hair cell factors, ATOH1 and POU4F3, were included based on manually identified binding sites in the proximal promoter (Fig. 3c). Candidates were narrowed down to the following transcription factors: ATOH1, EGR1, GATA3, POU4F3, and USF1. Together these transcription factors are involved in hair cell development and maintenance, synaptic plasticity, regeneration and tissue repair [28–38]. The predicted binding sites for these genes can be found in Fig. 3c.
EGR1 and ATOH1 modulate the isolated αAtp2b2 promoter in OC-1 cells
To determine which of these transcription factors modulates the αAtp2b2 transcript promoter, cDNA of the candidate transcription factors was cloned into a pcDNA vector upstream of IRES + EGFP. These vectors were co-transfected into OC-1 cells along with the [+572/−2133] αAtp2b2 luciferase construct. Gene expression levels were similar for all of the transcription factors post-transfection (Additional file 1: Table S1). The [+572/−2133] construct was used because it is the longest promoter construct that has activity over baseline in the OC-1 cell line (Fig. 2b). Figure 3c shows a cartoon of the sequence from base +101 to base −2133 in detail with the location of transcription factor binding sites and the CpG island. This experiment reveals that ATOH1 decreases αAtp2b2 promoter activity by 30% while EGR1 produces a significant twofold increase (Fig. 4a). Overexpression of GATA3, POU4F3 or USF1 had no significant effect on αAtp2b2 promoter activity (Fig. 4a).
Identifying the region of EGR1 activity in the αAtp2b2 promoter
ATOH1 and EGR1 are capable of modulating the [+572/−2133] αAtp2b2 luciferase construct in the OC-1 cell line, indicating functional significance in a mammalian cell model. The potential binding site of ATOH1 can be narrowed to a single site in the CpG island of the αAtp2b2 promoter (Fig. 3c). Conversely, EGR1 has numerous predicted binding sites throughout the αAtp2b2 promoter region. To map the location of functional EGR1 binding sites within the promoter, three different luciferase reporter constructs: [+572/−287], [572/−855], and [+572/−2133] were co-transfected with EGR1 or empty vector into OC-1 cells. All three constructs showed a significant increase in luciferase activity with co-expression of EGR1 (Fig. 4b). This data suggests that EGR1 binds the αAtp2b2 promoter between bases +572 and −287. Due to limitations of restriction sites in this region, putative binding sites for EGR1 were not narrowed further. The putative binding sites for EGR1 are contained in the CpG island which is 80% GC. Due to limitations of PCR and sequencing under these conditions, mutagenesis was not performed .
The effect of transcription factor overexpression on endogenous Atp2b2 and Atp2b4 transcript levels
The luciferase assay results indicate that EGR1 and ATOH1 have a functional effect on the isolated αAtp2b2 transcript promoter. However, this observation does not indicate that ATOH1 and EGR1 effect endogenous levels of Atp2b2 transcript expression. To test this, EGR1 and ATOH1 were overexpressed in N2A and OC-1 cells that transcribe Atp2b2 endogenously. Gene expression levels were similar for all of the transcription factors after transfection (Additional file 1: Table S1). Due to the low basal levels of Atp2b2 expression, TaqMan assays with high target specificity and accuracy were utilized for this experiment. The TaqMan probes do not discriminate between the different alternate start site variants of Atp2b2 so total Atp2b2 transcript levels were quantified. Atp2b4 transcript expression was also measured because there is evidence that Atp2b4 compensates for decreased expression of Atp2b2 in the stereocilia of Atp2b2 mutant mice .
In OC-1 and N2A cells, EGR1 overexpression increases transcript levels of Atp2b2 recapitulating our luciferase results (Fig. 5a, b). When ATOH1 is overexpressed in OC-1 and N2A cells, Atp2b2 transcript levels are unaffected (Fig. 5a, b). Although ATOH1 overexpression did not significantly inhibit Atp2b2 transcript, Atp2b4 levels were found to increase 3.5 fold over baseline in OC-1 cells (Fig. 5c). An ATOH1 driven increase in Atp2b4 transcript was not observed in N2A cells (Fig. 5d). Given that ATOH1 is specific to developing hair cells it is unsurprising that there is a cell-type specific effect. EGR1 overexpression has no significant effects on Atp2b4 expression in OC-1 or N2A cells (Fig. 5c, d). The luciferase data taken together with Atp2b2 and Atp2b4 expression data suggest interplay between EGR1 and ATOH1 in OC-1 cells. EGR1 increases Atp2b2 expression while ATOH1 inhibits the αAtp2b2 promoter and increases Atp2b4 expression (Figs. 4a, 5a, c).
Shift and supershift assays indicate that EGR1 binds to the αAtp2b2 promoter
It is likely that both EGR1 and ATOH1 play a role in modulation of Atp2bs. However, ATOH1 did not significantly affect endogenous levels Atp2b2 in OC-1 and N2A cells. Given this negative data, only EGR1 was included in the binding assays. To identify potential EGR1 binding sites, data from MatInspector was cross-referenced with TFBIND software predictions. Sites with matrix similarity of 0.75 or higher (highest predicted binding being 1.0) were included. This identified 12 overlapping binding sites that could be divided into seven different regions of the promoter (Fig. 3). Of particular interest are three sites in the CpG island near the TSS between base −71 and base −98 (Fig. 3). A probe targeting this region was designed to experimentally test the protein binding potential of this region. Two positive control probes were used to verify the shift assay, one derived from the EGR1 consensus sequence (GCGGGGGCG) [17, 18] and a gel shift probe commercially available from Santa Cruz Biotechnology (Table 3). EGR1 protein levels were measured in OC-1, HeLa and NIH/3T3 cells by immunoblotting. Protein expression of EGR1 was not detectable in OC-1 cells (data not shown) or in NIH/3T3 cells (Fig. 6a). EGR1 was found in the nucleus and cytoplasm of HeLa cells (Fig. 6a).
All EMSAs were performed using HeLa cell nuclear extracts (Fig. 6b, c). No protein binding activity was detected with two test probes for sequences either downstream or at the TSS (data not shown). A mobility shift is only observed with the DNA probe derived from the site upstream of the TSS [−71 to −98] (Fig. 6b). Formation of this shifted DNA–protein complex is inhibited when increasing concentrations of unlabeled [−71 to −98] DNA probe are added as competitor. Four nucleotides in EGR1 binding cluster of the [−71 to −98] DNA probe were modified to investigate sequence specificity of binding (Table 3). This modified probe (Mod) is unable to produce a shift when compared to the [−71 to −98] probe (Test) derived from the αAtp2b2 promoter and two positive control probes (Can and SC) (Fig. 6c). These experiments indicate that the shift is sequence specific but does not indicate which protein is binding to the probe.
To determine if EGR1 is the protein responsible for the mobility shift in Fig. 6a and b, antibody to EGR1 was added to the binding reactions. For these assays, the shifted complexes resolved into three distinct populations (Fig. 6d). The formation of complex II was dramatically inhibited upon addition of the EGR1 specific antibody but not the negative control antibody targeting ATP2B2 (Fig. 6d, Lane 4). These results suggest that EGR1 is binding to the [−71 to −98] probe in a sequence specific manner to stimulate αAtp2b2 promoter activity. Although there are different populations (I and III) that remain bound to the probe in the presence of EGR1 antibody it is likely that this is non-specific binding or it is binding of the probe to a protein besides EGR1. There are a number of predicted SP1 binding sites in the CpG island. Previously published experiments have exhibited the similarities between binding sites for SP1 and EGR1 [18, 42]. It is plausible that this is the identity of the protein that remains bound to the probe in the supershift assay (Fig. 6d, Lane 3).
Densitometry was employed to compare the standard shift assay against the test shift assay for each of the above experiments: self-competition, modified EGR1 probe and supershift. The average pixel density of the standard shift was significantly darker than the test shift in all three experiments (Fig. 6e). This quantification supports the conclusion that EGR1 binds to the [−71 to −98] probe in a sequence specific manner.
Initiation of gene transcription in Atp2bs and αAtp2b2
The Ca2+ ATPases play an important role in maintaining cell homeostasis but little is known about initiation of transcription of these genes. Atp2b2 has the fastest activation and is one of the primary ATP pumps found in the brain and sensory epithelium. Investigations of this gene have highlighted its pivotal role in normal function of the auditory system as well as the vestibular system [19, 6, 7, 12, 43, 44]. Of the two transcripts found in the brain (α and β) only the α transcript is found in the auditory hair cells of young mice (postnatal day 9). This study focuses on the α transcript because of its pivotal role in the normal function of the auditory system.
Identification of the αAtp2b2 minimal promoter and repressor elements
Studies of the genomic region surrounding the Atp2b2 gene indicate that the minimal promoter of the α transcript is likely located directly adjacent to the TSS. Our luciferase assay data narrows this to the region 2133 bp upstream of the TSS. The abolishment of promoter activity upon removal of the CpG island and the 5′ UTR narrows the minimal promoter to the 287 bps upstream of the αAtp2b2 TSS. In both OC-1 and OC-2 cell lines, activity decreases as the promoter is extended in the 5′ direction. This is likely due to repressor elements upstream of the minimal promoter. Interestingly, the activity of the promoter falls off faster in OC-2 cells suggesting cell type specific repression and modulation of the promoter.
The effect of ATOH1 and EGR1 on Atp2b2
In OC-1 cells, we find that ATOH1 decreases αAtp2b2 proximal promoter activity by ~30%. However, we did not detect a significant effect of ATOH1 on endogenous levels of Atp2b2 in N2A or OC-1 cells. This contradicts the luciferase assay which predicts that ATOH1 would inhibit endogenous levels of Atp2b2 in OC-1 cells. A caveat of this result is that Atp2b2 is very close to the limit of detection in OC-1 cells making any decrease in expression challenging to measure. It is therefore possible that ATOH1 does inhibit endogenous levels of Atp2b2 in OC-1 cells but we were unable to detect it. This caveat does not apply to the EGR1 overexpression assay because EGR1 increases expression of Atp2b2 in OC-1 cells (and in N2A cells). Any increase in Atp2b2 expression is in the detectable range of the assay. In addition to increasing endogenous levels of Atp2b2 in OC-1 and N2A, EGR1 also increases the activity of the αAtp2b2 promoter twofold over empty vector. We conclude that EGR1 increases promoter activity of αAtp2b2 and expression of Atp2b2 transcript in OC-1 and N2A cells while ATOH1 represses promoter activity of αAtp2b2 in OC-1 cells.
The interplay of EGR1, ATOH1, Atp2b2 and Atp2b4
In OC-1 cells, overexpression of ATOH1 increased Atp2b4 transcript 3.5 times over baseline. However, the effect was not recapitulated in N2A cells. Our data suggests that in organ of Corti derived cells ATOH1 causes increased expression of Atp2b4. SHIELD data indicates that ATOH1 is expressed throughout development in hair cells with a peak around postnatal day 0  (Table 5). Although SHIELD data shows only low Atp2b4 expression, other studies have exhibited expression of Atp2b4 in the hair cells at postnatal day 12 . This expression may be driven by ATOH1.
Our data shows that EGR1 increases expression of Atp2b2 while ATOH1 represses the αAtp2b2 promoter. SHIELD data indicates that EGR1 levels are high in the hair cells of mice at embryonic day 16 and remain high with a peak at postnatal day 0. Atp2b2 expression steadily rises through all of the ages assayed and it is possible that Atp2b2 levels are inhibited in early development by ATOH1 and rises as ATOH1 turns off (Table 5).
This exhibits an interesting interplay between Atp2b2 and Atp2b4. It is known that ATPases compensate for each other in knockout animals . Notably, in Atp2b2 null mice, Atp2b4 is expressed abnormally in stereocilia into adulthood . In OC-1 cells, ATOH1 increases expression of Atp2b4 and represses the promoter of Atp2b2. This suggests that ATOH1 may be a developmental switch for Atp2bs in hair cells. For example, as ATOH1 turns off it may allow for Atp2b2 to be turned on.
The significance of EGR1 and Atp2b2
Although EGR1 increases expression of Atp2b2 in both OC-1 and N2A cells, the low level of Atp2b2 expression in both cell lines suggests there are other factors involved in transcriptional activation of this gene. This is unsurprising given the complexity of this gene and the necessity for tight control of Atp2b2 for normal physiology. The effect of EGR1 seems specific to Atp2b2 as overexpression did not increase expression of the closely related Atp2b4 gene in either OC-1 or N2A cells. Investigation of the proximal promoter of Atp2b4 on the UCSC genome browser suggests this is due to the lack of a CpG island. EGR1 is important for neuronal and sensory systems through involvement in synaptic plasticity, retinal formation and acoustic trauma. Additionally, it is implicated in auditory hair cell regeneration pathways [33,34,35,36,37,38].
This is the first investigation of the Atp2b2 promoter and is the second investigation of a promoter in the vast Atp2b family. This study identified the minimal promoter elements of αAtp2b2 including a CpG island. Additionally, our investigations indicate that ATOH1 and EGR1 modulate αAtp2b2 promoter activity. There appears to be interplay between EGR1, ATOH1, Atp2b2 and Atp2b4 in OC-1 cells; further studies are necessary to confirm these interactions and determine their physiological relevance in vivo. EGR1 was also found to bind to a specific sequence 71 bases upstream of the TSS in the CpG island of the αAtp2b2 promoter. In conclusion, we have elucidated transcriptional regulatory mechanisms for Atp2b2, identified transcription factors that modulate expression of Atp2b2 in the brain and peripheral auditory system, and established a foundation for future studies investigating expression of αAtp2b2.
transcriptional start site
- Atp2b :
plasma membrane calcium ATPase
DNase I hypersensitivity
- RNA Pol II:
RNA polymerase II
histone H3 lysine 4 trimethylation
electrophoretic mobility shift assay
atonal bHLH transcription factor 1
early growth response protein 1
POU domain class 4 transcription factor 3
trans-acting T cell-specific transcription factor GATA3
upstream stimulatory factor 1
- Ca2+ :
- dfw :
internal ribosome entry site
- Hmbs :
- Gamma (γ) actin:
Tempel BL, Shilling DJ. The plasma membrane calcium ATPase and disease. Subcell Biochem. 2007;45:365–83.
Strehler EE. Plasma membrane calcium ATPases: from generic Ca2+ sump pumps to versatile systems for fine-tuning cellular Ca2+. Biochem Biophys Res Commun. 2015;460:26–33.
Prasad V, Okunade GW, Miller ML, Shull GE. Phenotypes of SERCA and PMCA knockout mice. Biochem Biophys Res Commun. 2004;322:1192–203.
Hirawa N, Fujiwara A, Umemura S. ATP2B1 and blood pressure: from associations to pathophysiology. Curr Opin Nephrol Hypertens. 2013;22:177–84.
Ho PW-L, Pang SY-Y, Li M, Tse ZH-M, Kung MH-W, Sham P-C, Ho S-L. PMCA4 (ATP2B4) mutation in familial spastic paraplegia causes delay in intracellular calcium extrusion. Brain Behav. 2015;5:1–2.
Street VA, McKee-Johnson JW, Fonseca RC, Tempel BL, Noben-Trauth K. Mutations in a plasma membrane Ca2+-ATPase gene cause deafness in deafwaddler mice. Nat Genet. 1998;19:390–4.
Konrad-Martin D, Norton JS, Mascher KE, Tempel BL. Effects of PMCA2 mutation on DPOAE amplitudes and latencies in deafwaddler mice. Hear Res. 2001;151:205–20.
Schultz JM, Yang Y, Caride AJ, Filoteo AG, Penheiter AR, Lagziel A, Morell RJ, Mohiddin SA, Fananapazir L, Madeo AC. Modification of human hearing loss by plasma-membrane calcium pump PMCA2. N Engl J Med. 2005;352:1557–64.
Ficarella R, Di Leva F, Bortolozzi M, Ortolano S, Donaudy F, Petrillo M, Melchionda S, Lelli A, Domi T, Fedrizzi L, Lim D, Shull GE, Gasparini P, Brini M, Mammano F, Carafoli E. A functional study of plasma-membrane calcium-pump isoform 2 mutants causing digenic deafness. Proc Natl Acad Sci USA. 2007;104:1516–21.
Roberts-Thomson SJ. Plasma membrane calcium pumps and their emerging roles in cancer. World J Biol Chem. 2010;1:248.
Yang W, Liu J, Zheng F, Jia M, Zhao L, Lu T, Ruan Y, Zhang J, Yue W, Zhang D, Wang L. The evidence for association of ATP2B2 polymorphisms with autism in Chinese Han population. PLoS ONE. 2013;8:e61021.
Watson CJ, Tempel BL. A new Atp2b2 deafwaddler allele, dfwi5, interacts strongly with Cdh23 and other auditory modifiers. Hear Res. 2013;304:41–8.
Penheiter AR, Filoteo AG, Croy CL, Penniston JT. Characterization of the deafwaddler mutant of the rat plasma membrane calcium-ATPase 2. Hear Res. 2001;162:19–28.
McCullough BJ, Tempel BL. Haplo-insufficiency revealed in deafwaddler mice when tested for hearing loss and ataxia. Hear Res. 2004;195:90–102.
Du Y, Carlock L, Kuo TH. The mouse plasma membrane Ca2+ pump isoform 1 promoter: cloning and characterization. Arch Biochem Biophys. 1995;316:302–10.
Rivolta MN, Grix N, Lawlor P, Ashmore JF, Jagger DJ, Holley MC. Auditory hair cell precursors immortalized from the mammalian inner ear. Proc R Soc Lond B Biol Sci. 1998;265:1595–603.
Daily K, Patel VR, Rigor P, Xie X, Baldi P. MotifMap: integrative genome-wide maps of regulatory motif sites for model species. BMC Bioinform. 2011;12:495.
Cao X, Mahendran R, Guy GR, Tan YH. Detection and characterization of cellular EGR-1 binding to its recognition site. J Biol Chem. 1993;268:16949–57.
Silverstein RS, Tempel BL. Atp2b2, encoding plasma membrane Ca2+-ATPase type 2, (PMCA2) exhibits tissue-specific first exon usage in hair cells, neurons, and mammary glands of mice. Neuroscience. 2006;141:245–57.
Thurman RE, Rynes E, Humbert R, Vierstra J, Maurano MT, Haugen E, Sheffield NC, Stergachis AB, Wang H, Vernot B, Garg K, John S, Sandstrom R, Bates D, Boatman L, Canfield TK, Diegel M, Dunn D, Ebersol AK, Frum T, Giste E, Johnson AK, Johnson EM, Kutyavin T, Lajoie B, Lee B-K, Lee K, London D, Lotakis D, Neph S, et al. The accessible chromatin landscape of the human genome. Nature. 2012;489:75–82.
Lenhard B, Sandelin A, Carninci P. Metazoan promoters: emerging characteristics and insights into transcriptional regulation. Nat Rev Genet. 2012;13(4):233-45. doi:10.1038/nrg3163.
ENCODE Project Consortium. A user’s guide to the encyclopedia of DNA elements (ENCODE). PLoS Biol. 2011;9:e1001046.
Consortium TEP. An integrated encyclopedia of DNA elements in the human genome. Nature. 2012;489:57–74.
Rivolta MN, Holley MC. Cell lines in inner ear research. J Neurobiol. 2002;53:306–18.
Cartharius K, Frech K, Grote K, Klocke B, Haltmeier M, Klingenhoff A, Frisch M, Bayerlein M, Werner T. MatInspector and beyond: promoter analysis based on transcription factor binding sites. Bioinform Oxf Engl. 2005;21:2933–42.
Tsunoda T, Takagi T. Estimating transcription factor bindability on DNA. Bioinform Oxf Engl. 1999;15:622–30.
Shen J, Scheffer DI, Kwan KY, Corey DP. SHIELD: an integrative gene expression database for inner ear research. Database J Biol Databases Curation. 2015;2015:bav071.
Pan N, Jahan I, Kersigo J, Duncan JS, Kopecky B, Fritzsch B. A novel Atoh1 “self-terminating” mouse model reveals the necessity of proper Atoh1 level and duration for hair cell differentiation and viability. PLoS ONE. 2012;7:e30358.
Alvarado DM, Veile R, Speck J, Warchol M, Lovett M. Downstream targets of GATA3 in the vestibular sensory organs of the inner ear. Dev Dyn. 2009;238:3093–102.
van der Wees J, van Looij MA, de Ruiter MM, Elias H, van der Burg H, Liem S-S, Kurek D, Engel JD, Karis A, van Zanten BG, De Zeeuw CI, Grosveld FG, van Doorninck JH. Hearing loss following Gata3 haploinsufficiency is caused by cochlear disorder. Neurobiol Dis. 2004;16:169–78.
Hertzano R. Transcription profiling of inner ears from Pou4f3ddl/ddl identifies Gfi1 as a target of the Pou4f3 deafness gene. Hum Mol Genet. 2004;13:2143–53.
Horbach T, Götz C, Kietzmann T, Dimova EY. Protein kinases as switches for the function of upstream stimulatory factors: implications for tissue injury and cancer. Front Pharmacol. 2015;6:3.
Lomax MI, Gong T-W, Cho Y, Huang L, Oh S-H, Adler HJ, Raphael Y, Altschuler RA. Differential gene expression following noise trauma in birds and mammals. Noise Health. 2001;3:19–35.
Schippert R, Schaeffel F, Feldkaemper MP. Microarray analysis of retinal gene expression in Egr-1 knockout mice. Mol Vis. 2009;15:2720–39.
Wei F, Xu ZC, Qu Z, Milbrandt J, Zhuo M. Role of Egr1 in hippocampal synaptic enhancement induced by tetanic stimulation and amputation. J Cell Biol. 2000;149:1325–34.
Hu SS, Mei L, Chen JY, Huang ZW, Wu H. Expression of immediate-early genes in the inferior colliculus and auditory cortex in salicylate-induced tinnitus in rat. Eur J Histochem EJH. 2014;58:2294.
Hawkins RD, Bashiardes S, Powder KE, Sajan SA, Bhonagiri V, Alvarado DM, Speck J, Warchol ME, Lovett M. Large scale gene expression profiles of regenerating inner ear sensory epithelia. PLoS ONE. 2007;2:e525.
Smith ME, Rajadinakaran G. The transcriptomics to proteomics of hair cell regeneration: looking for a hair cell in a haystack. Microarrays Basel Switz. 2013;2:186–207.
Scheffer DI, Shen J, Corey DP, Chen Z-Y. Gene expression by mouse inner ear hair cells during development. J Neurosci. 2015;35:6366–80.
Rivolta MN. Transcript profiling of functionally related groups of genes during conditional differentiation of a mammalian cochlear hair cell line. Genome Res. 2002;12:1091–9.
McDowell DG, Burns NA, Parkes HC. Localised sequence regions possessing high melting temperatures prevent the amplification of a DNA mimic in competitive PCR. Nucleic Acids Res. 1998;26:3340–7.
Deckmann K, Rörsch F, Geisslinger G, Grösch S. Identification of DNA–protein complexes using an improved, combined western blotting-electrophoretic mobility shift assay (WEMSA) with a fluorescence imaging system. Mol Biosyst. 2012;8:1389.
Bortolozzi M, Brini M, Parkinson N, Crispino G, Scimemi P, De Siati RD, Di Leva F, Parker A, Ortolano S, Arslan E, Brown SD, Carafoli E, Mammano F. The novel PMCA2 pump mutation Tommy impairs cytosolic calcium clearance in hair cells and links to deafness in mice. J Biol Chem. 2010;285:37693–703.
Wood JD, Muchinsky S, Filoteo A, Penniston J, Tempel B. Low endolymph calcium concentrations in deafwaddler2J mice suggest that PMCA2 contributes to endolymph calcium maintenance. J Assoc Res Otolaryngol. 2004;5:99–110.
Furuta H, Luo L, Hepler K, Ryan AF. Evidence for differential regulation of calcium by outer versus inner hair cells: plasma membrane Ca-ATPase gene expression. Hear Res. 1998;123:10–26.
RRM conducted the majority of the experiments, analyzed the results, and wrote the paper. JL conceived the original cloning strategies for the αAtp2b2 promoter and worked closely with RRM on the design and execution of luciferase assays and cloning. BLT contributed fundamental ideas for the direction of the project and wrote the paper with RRM. All authors read and approved the final manuscript.
We thank Prof. M. Holley for the gift of OC-1 and OC-2 cell lines. We also acknowledge Dr. B. Peguero, Profs. E. Wang and J. Stone for their helpful comments on this manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
Raw data for luciferase assays, qPCR and EMSAs will be made available upon request.
Ethics approval and consent to participate
No human or animal subjects were used in this study. Neuro-2A and HeLa cells were originally obtained from ATCC® (ATCC®: Neuro-2a/CCL-131™ and HeLa/CCL-2™). OC-1 and OC-2 cells were originally obtained from Prof. Matthew Holley (University of Sheffield). Cell lines were subsequently passaged and maintained in the Tempel lab.
This work was supported in whole or part by National Institute of Health Grants RO1 DC02739 (BLT), R21 DC014864 (BLT), P30 DC04661 (BLT), T32 DC005361 (RRM). Reviewers for the funding body gave helpful feedback on the original grant outlining this study. The funding body had no role in the study design, data collection, analysis or interpretation or the writing of the manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1. Change in Transcription Factor Expression in OC-1 and N2A cells After Transfection. Normalized expression of transcription factors (AU) in OC-1 and N2A cells at baseline and after transfection with transcription factor constructs.
About this article
Cite this article
Minich, R.R., Li, J. & Tempel, B.L. Early growth response protein 1 regulates promoter activity of α-plasma membrane calcium ATPase 2, a major calcium pump in the brain and auditory system. BMC Molecular Biol 18, 14 (2017). https://doi.org/10.1186/s12867-017-0092-1