IRF-1 regulates alternative mRNA splicing of carcinoembryonic antigen-related cell adhesion molecule 1 (CEACAM1) in breast epithelial cells generating an immunoreceptor tyrosine-based inhibition motif (ITIM) containing isoform
Molecular Cancer volume 13, Article number: 64 (2014)
Interferon regulatory factor-1 (IRF-1) is a master regulator of IFN-γ induced gene transcription. Previously we have shown that IRF-1 transcriptionally induces CEACAM1 via an ISRE (Interferon-Stimulated Response Element) in its promoter. CEACAM1 pre-mRNA undergoes extensive alternative splicing (AS) generating isoforms to produce either a short (S) cytoplasmic domain expressed primarily in epithelial cells or as an ITIM-containing long (L) isoform in immune cells.
The transcriptional and molecular mechanism of CEACAM1 minigenes AS containing promoter ISREs mutations in the breast epithelial, MDA-MB-468, cell line was detected using flow cytometry. In addition, transcriptome sequencing was utilized to determine whether IRF-1 could direct the AS of other genes as well. Tumor xenografts were used to evaluate CEACAM1 isoform expression on the leading edge of breast tumor cells.
In the present study, we provide evidence that CEACAM1’s promoter and variable exon 7 cross-talk allowing IRF-1 to direct AS events. Transcriptome sequencing shows that IRF-1 can also induce the global AS of genes involved in regulation of growth and differentiation as well as genes of the cytokine family. Furthermore, MDA-MB-468 cells grown as tumor xenografts exhibit an AS switch to the L-isoform of CEACAM1, demonstrating that an in vivo inflammatory milieu is also capable of generating the AS switch, similar to that found in human breast cancers Mol Cancer 7:46, 2008.
The novel AS regulatory activities attributed to IRF-1 indicate that the IFN-γ response involves a global change in both gene transcription and AS in breast epithelial cells.
Interferons (IFNs) are a group of immunomodulatory proteins made and released by host cells in response to the presence of pathogens such as viruses, bacteria, parasites or tumor cells. IFNs are primarily known for their role in inhibiting viral infections and in stimulating the immune system . IFN-γ, a type II interferon, is a pleiotropic cytokine involved in the regulation of immune and inflammatory responses, including the activation, growth and differentiation of B-cells, T-cells, natural killer (NK) cells and macrophages. Aberrant IFN-γ expression, associated with autoinflammatory and autoimmune diseases, acts as a pro-inflammatory cytokine . IFN-γ is produced predominantly by natural killer and natural killer T cells as part of the innate immune response, and by CD4 Th1 and CD8 cytotoxic T lymphocyte effector T cells once antigen-specific immunity develops .
IFN-γ induces the JAK/STAT pathway leading to the activation and binding of transcriptional activators that induce expression of IFN-stimulated genes, notably IRF-1. IRF-1 is expressed at low levels in unstimulated cells and recognizes the ISRE found in most IFN-inducible gene promoters containing repeating short GAAA core sequences [4, 5]. We have previously shown that CEACAM1 transcription can be induced by IFN-γ via an ISRE in the CEACAM1 promoter [6–8]. CEACAM1 is a member of the immunoglobulin super family of glycoproteins . It is expressed on the surface of epithelial and activated endothelial cells, as well as on cells from the immune system and plays a role in a variety of cellular processes including cell-cell adhesion, differentiation, apoptosis and immune response .
In humans and rat, CEACAM1 pre-mRNA undergoes extensive AS generating isoforms consisting of an N-terminal and a variable number of multiple extracellular Ig-like domains, a transmembrane domain, and either short (CEACAM1-S) or long (CEACAM1-L) cytoplasmic domains [11–13]. The short cytoplasmic domain isoform, the predominant isoform in epithelial cells, has been shown to bind actin, tropomyosin, calmodulin, and annexin II and is involved in lumen formation [14–16]. In normal breast, it is the predominant isoform, while in breast cancer the S/L ratio is greatly reduced . The long cytoplasmic domain isoform with its two cytosolic phosphotyrosine residues and immunoreceptor tyrosine-based inhibitory motifs (ITIM), predominant in immune cells, bind SHP-1 when phosphorylated and convey inhibitory activities to CEACAM1-L . Furthermore, in T lymphocytes CEACAM1-L isoforms are significantly upregulated at the cell surface in response to IL-2 treatment, mediating their cell adhesion to other lymphocytes or tumor cells [19, 20].
The observation that CEACAM1 functions in recognition and binding between cells with adhesion properties is integral to its signal transduction leading to growth suppressive activities . Indeed, CEACAM1 is a tumor suppressor gene that is down-regulated in most tumor cells such as colon, breast and prostrate cancer [21–23]. Introduction of CEACAM1-L but not CEACAM1-S cDNA by adenovirus as a gene delivery vehicle reduces tumorigenicity of breast cancer cells . These findings are consistent with the idea that the expression of CEACAM1 plays a role in keeping cells in a differentiated state, while its down-regulation may be an essential event in the development of epithelial malignancy . However, the implication of the regulation of its AS generated isoforms is less well understood.
AS is a genomic evolutionary strategy that generates proteome expansion in metazoans [25, 26]. cis- Regulatory elements known as exonic splicing enhancers and exonic splicing silencers (ESS ) provide the binding platform for splicing regulators to promote or inhibit the use of nearby splice-sites (ss). AS effectors include members of the heterogeneous nuclear ribonucleoprotein (hnRNP) family and serine/arginine-rich (SR) protein families are among the most common splicing regulators [27, 28]. hnRNPs rapidly associate with nascent transcripts and repress certain AS events leading to exon skipping events . Three abundant hnRNPs (hnRNP L, hnRNP A1, and hnRNP M) are implicated in many AS events in human and several other eukaryotes and play a role in regulating splicing events by inhibiting the use of splice-sites from the spliceosome . Previously we showed that hnRNP L and hnRNP A1 interact specifically with exon 7 to control splice-site recognition in the formation of CEACAM1-S mRNA. hnRNP M is important for the AS of CEACAM1-L mRNA .
We recently showed that IFN-γ altered the ratio of splice isoforms of CEACAM1 through the presumable induction of its down-stream effector IRF-1 . The goal of the present study was to determine whether forced expression of IRF-1 directly plays a role in regulating CEACAM1 AS. Using GFP-driven reporter constructs containing the CEACAM1 genomic promoter and its exons 6-7-8 (heretofore called CAM1 6-7-8) cassettes, we show that AS depends on IRF-1 to couple transcriptional complex communication with the splicing machinery. Transcriptome sequencing analysis shows that IRF-1 can induce the AS of other genes as well, including those involved in membrane regeneration and repair, microtubule assembly, and tumor suppressor functions. We also report extensive changes in transcriptional and AS regulation of cytokines in breast cancer epithelial cells in response to IRF-1 treatment. The AS regulatory activities of IRF-1 demonstrate that the mechanism of the interferon response includes a global change in both gene transcription and AS in breast epithelial cells.
Ad-IRF-1 treatment induces the CEACAM1-L isoform in a breast cancer line
Although normal breast epithelial cells predominantly express the CEACAM1-S isoform, we have observed a switch of the ratio of S/L isoforms in favor of the L isoform in breast cancer by an unknown mechanism . More recently, while investigating changes in transcription factor binding to the CEACAM1 promoter in MDA-MB-468 breast cells after treatment with IFN-γ, we observed the induction of CEACAM1-L isoform suggesting that the AS mechanism involves the ISRE of CEACAM1’s promoter . Since the ISRE involves binding of IRF-1, we investigated if forced expression of IRF-1 could directly generate the L isoform.
The AS of CEACAM1 is dependent on splicing regulatory domains in exon 7 . These cis-elements are capable of recruiting various trans-regulatory proteins to direct splice-site recognition involved in exon 7 definition. Specifically, exon 7 contains two distinct ESSs that recruit hnRNP L and hnRNP A1 to direct exon-skipping by the spliceosome whereas hnRNP M leads to exon-inclusion, in the production of the L-isoform (Figure 1A).
Ectopically expressed IRF-1 in the breast carcinoma cell line MDA-MB-468 has been previously described using a recombinant adenovirus (Ad-IRF-1) as a gene transfer vehicle containing the mouse IRF-1 gene under a constitutive, cytomegalovirus promoter [29, 30]. Using this system, we observed the expression of IRF-1 but not IRF-2, a repressor of CEACAM1 transcription and other target genes , 12 h after infection with an apparent peak at 24 h (Figure 1B and Additional file 1: Figure S1A). IRF-1 expression was absent in the Ad-null control-infected cells. The induction of the cyclin-dependent kinase inhibitor p21 (WAF1/CIP1) was included to demonstrate that effector IRF-1 was functionally competent .
We next addressed whether Ad-IRF-1 in MDA-MB-468 cells could replicate IFN-γ-directed AS of CEACAM1. CEACAM1-4L and –S mRNA isoforms, but not CEACAM1-3L and -3S isoforms, are strongly up-regulated by Ad-IRF-1 as compared to virus vector control (Figure 1C, lane 4 versus lane 3). CEACAM1-4L and CEACAM1-4S contain four extracellular immunoglobulin-like domains while CEACAM1-3L and CEACAM1-3S only three. Notably, in both isoforms we detected strong activation of spliceosome processing (more than 2-fold) of exon 7 containing mRNA transcripts between 12-24 h post infection (CEACAM1-4L, p < 0.001 and CEACAM1-3L, p < 0.01; Figure 1C). We next tested MCF7 cells, derived from breast adenocarcinoma, with Ad-IRF-1 treatment to see if we could detect a similar AS phenotype (Additional file 2: Figure S2A-B). We did not detect evidence of CEACAM1 mRNA using quantitative real-time PCR analyses in untreated cells consistent with previous observations . By contrast, stimulation of cells by Ad-null treatment caused a significant and strong up-regulation of CEACAM1-S but not CEACAM1-L mRNA as early as 6 h post-treatment. We also observed that IRF-1 transcriptionally induces more CEACAM1-S but not CEACAM1-L mRNA at all time points tested. Only at 24 h we observed the emergence of CEACAM1-L mRNA, albeit approximately 1000-fold less. This is in contrast to our observations in MDA-MB-468 cells that production of the L-isoform starts as early as 6 hrs after Ad-IRF-1 treamment (Additional file 1: Figure S1B). In the case of the MCF7 cell line response to Ad-null viral infection, it is likely that an additional immune response is generated unrelated to splicing; e.g., a TLR9 response could account for this.
Total lysates from treated and untreated Ad-IRF-1 cells were next subjected to Western blots to establish correlative protein to RNA expression levels (Figure 1D and Additional file 3: Figure S3A-B). Using antibody 229 (which detects CEACAM1-L; green channel) and T84.1 (which detects both CEACAM1 isoforms, red channel) and fluorescent secondary antibodies we were able to distinguish between low levels of CEACAM1-S in untreated and viral control cells and stimulated production of CEACAM1-L after Ad-IRF-1 treatment (yellow band, Figure 1D, lane 4). After Ad-IRF-1 treatment, protein levels of CEACAM1-S increased over basal levels by Western blot and FACS analysis (Figure 1D, lane 4 and Additional file 4: Figure S4). We conclude that a novel mechanism for IRF-1 function has been uncovered and it involves post-transcriptional AS.
An ISRE with a GG->CC substitution disrupts IRF-1 binding
To elucidate the molecular mechanism of IRF-1 dependent AS, we created ISRE mutants that contained purine to pyrimidine substitutions. Previously, it was shown that the major groove of IRF-1 contacts bases that localize to a well-established consensus GAAA sequence within the 13-bp PRD I element . Examination of CEACAM1’s core ISRE shows three such motifs flanking a central 223-GG-224 dinucleotide (Figure 2B). We created mutant ISRE A->T by introducing AAA->TTT or TTC substitutions along the ISRE. Additionally, we previously reported an in vivo footprint in the ISRE (-223/4 nt relative to the TSS) that indicated 223-GG-224 bases interact directly with IRF-1 on the sense but not antisense strand . Therefore, mutant ISRE GG->CC was created to test this target sequence.
We hypothesized that weakening IRF-1’s association with CEACAM1’s ISRE would inhibit AS events. To test this hypothesis, we first cloned WT and mutant ISREs (from a 1.135-kb promoter CEACAM1 DNA fragment relative to the TSS) in translational fusion to a GFP reporter (Figure 2A). We first tested whether our system was responsive to infection with Ad-IRF-1 using the ISRE WT construct and performed RT-PCR to determine whether we could detect a stable GFP-containing message (Figure 2C). Our data of IRF-1 treatment in ISRE WT compared to Vector alone demonstrates responsiveness to Ad-IRF-1 delivery [lane 6 (0.87) compared to lane 4 (0.58) GFP:GAPDH ratio]. These ISRE promoter-GFP constructs were next transiently transfected into MDA-MB-468 cells. To avoid misinterpretation of weak reporter activity resulting from potential loss-of-function phenotypes, we also co-transfected WT CEACAM1 ISRE fused to mCherry, as a transfection control (Figure 2A). Our dual reporter system provides in effect an in vivo flow-cytometry read-out of the efficacy of promoter mutations as a function of transfection efficiency, as previously shown . As shown in Figure 2D, one of three independent experiments, and summarized in Figure 2E, mock-treated ISRE WT samples but not vector controls produced double positive GFP/mCherry cells (21.1% compared to 3.6%). Interestingly, we note that the percent double positive cells observed for the ISRE A->T mutant decreased by only 28.9% (15.0/21.1), whereas by comparison, the mutant ISRE GG->CC decreased by a potent 94.3% (1.2/21.1; p < 0.01). We found that the addition of exogenous IRF-1 increased reporter activity by almost 1.5-fold in the ISRE WT reporter as compared to the vector control which showed no significant difference. Notably, the target mutant ISRE GG->CC predicted by footprinting showed more than a 5-fold increase in GFP expression [6.3/1.2; p < 0.01] when presented with Ad-IRF-1 compared to the weaker result by the consensus mutant ISRE A->T. We next harvested and collected RNA from these cells and tested by RT-PCR whether we could detect changes in GFP expression at the level of RNA (data not shown). Consistent with our flow cytometry observations, we conclude that transcriptional activation of GFP by ISRE WT and ISRE GG->CC promoter is enhanced by the addition of Ad-IRF-1 (Figure 2E). We also conclude that no other transcriptional factors (SP1, USF1) substitute for IRF-1 when ISRE promoter association has been dysregulated (Figure 2A; see genetic organization of CEACAM1 promoter).
To characterize the high inhibition observed for the ISRE GG->CC promoter mutant, we next performed a fluorescent Electromobility Shift Assay (fEMSA) using IRDye (infrared dye)-labeled 60-mer oligonucleotides (Figure 3A). We chose the Epstein-Barr virus internal ribosome entry site (EBNA) sequence to serve as our positive control for EMSA binding. From the data presented we show that the binding affinity for IRF-1 by the ISRE GG->CC mutant (green channel) is significantly impaired when compared to the ISRE WT sequence (Figure 3B, lanes 4-7 versus lanes 9-12; red channel). At the lowest ratio of IRF-1 to ISRE DNA (17:0.3), we observed the majority of ISRE WT DNA formed a slow migrating IRF-1 containing complex (denoted C1 in the figure) and a faster migrating complex (denoted C2; lane 9). By contrast, the ISRE GG->CC mutant made no C1 complex, while 50% of the bound DNA was in the C2 complex and another 50% remained in an unbound state, suggesting that the dinucleotide mutation caused a weakening of IRF-1 binding (lane 4). Since complex formation still occurs in the mutant, it suggests that CEACAM1’s ISRE containing repeat-short GAAA sequences still play an important role in binding .
We also tested whether IRF-2 binding is affected by the target mutant ISRE GG->CC. Previous studies have shown that IRF-2 competitively inhibits IRF-1 by binding to the same regulatory elements causing inhibition of IRF-1-mediated transcriptional activation of IFN-inducible genes . As shown in Figure 3C, fEMSAs performed with IRF-2 and the ISRE GG->CC mutant caused a similar disruption of protein–DNA interaction with the notable exception that at the lowest ratio of IRF-2 to ISRE DNA tested (17:0.3), some DNA remains unbound for both the ISRE WT and ISRE GG->CC DNA duplex (Figure 3C, lanes 9 versus 4). This suggests IRF-1 has a greater affinity for CEACAM1’s ISRE than IRF-2 and under natural conditions this difference could be important for modulating transcriptional activity at the promoter. Competition assays were performed with increasing amounts of competitor ISRE-WT DNA in the presence of 0.3 pmol ISRE GG->CC DNA and either 67 pmol IRF-1 or IRF-2 protein (Figure 3C). We observed competitor DNA caused accumulation of complex C2 in ISRE GG->CC DNA (see lanes 11 and 12 versus 5 and 6) which agrees with the predicted enhanced affinity of ISRE WT for both IRF-1 and IRF-2. Addition of 3.2 pmol competitor WT DNA, which produced all three binding states (unbound, yellow C2 and C1 complexes which is composite of red and green channel, lanes 6 and 12), indicates complete saturation of exogenous IRF protein. Taken together, these data validate the importance of the dinucleotide 223-GG-224 in the ISRE of CEACAM1.
Splicing reporter constructs show production of L-isoform requires only IRF-1
Having demonstrated GFP fluorescence correlates with induction of transcriptional regulation in our Ad-IRF-1 system, we next generated a tripartite splicing reporter construct consisting of a WT or mutant ISRE fused to the CAM1 6-7-8 minigene immediately upstream of the GFP reading frame such that fluorescence activity occurs from a single mRNA (Figure 4A). Essentially, these splicing reporters utilize AS mechanisms to modulate on–off expression of GFP such that exon inclusion but not exon skipping of a variable region puts GFP in frame. As shown in Figure 4B, one of three independent experiments, and summarized in Figure 4C, mock-treated ISRE WT samples but not vector controls produced double positive GFP/mCherry cells (6.3% compared to 4.9%). This correlated well with the fact that the default splicing program for CEACAM1 in MDA-MB-468 cells is the short isoform (see Figure 1C, lanes 2-3). Intriguingly, the percent double positive cells observed for the ISRE A->T mutant did not decrease significantly in comparison to the ISRE WT (6.4% compared to 6.3%) and the ISRE GG->CC mutant was now similar to background vector control levels (4.9% compared to 4.9%). In direct contrast to our promoter-only driven GFP reporters (Figure 2), the addition of exogenous IRF-1 strongly increased reporter activity by a potent 2.9-fold in the ISRE WT reporter [18.2/6.3; p < 0.01] as compared to the vector control that showed no significant difference. Flow cytometry of the ISRE GG->CC mutant revealed similar trends in its ability to stimulate GFP expression (Figure 4C; 13.75% compared to 4.9%; p < 0.05). From these data, we conclude that IRF-1 can activate spliceosome processing leading to the generation of CEACAM1-L. Notably, recovery of double positive GFP/mCherry cells in the ISRE A->T and ISRE GG->CC mutants to levels approaching the ISRE WT strongly indicated the possibility that regulatory elements were compensating for the strong differential potencies of these mutants.
Exon 7 directs AS through promoter activation of CEACAM1
We next tested if the promoter itself was responsible for recruiting splicing factors, such as hnRNPs, to the site of transcription. Proteins important for initiation of the transcription machinery, such as Med23 of the Mediator Complex, have been shown to interact with hnRNP L in the recruitment function for the preinitiation complex formation of the spliceosome . In previous studies, we have demonstrated that binding of hnRNP L and hnRNP A1 bind to exonic silencer sequences within exon 7 of CEACAM1 to repress exon inclusion . Additionally, we further demonstrated an extensive spliceosome regulatory complex specifically assembles on exon 7, presumably containing a multiplicity of splicing factors. We therefore proceeded to test the hypothesis that exon 7 RNA acts as a regulatory element in the promoter activation of CEACAM1.
Our experimental approach consisted of coupling ISRE promoter mutations with substitution of the cis-elements involved in hnRNP L and hnRNP A1 binding in exon 7 and studying the effects of splicing regulation on the production of CEACAM1-L:GFP. To effectively delete the binding sites of these splicing factors, we targeted two regions for further mutagenesis, namely E7-5 (nt 21-30) and E7-9 (nt 41-50) of exon 7 (total 53 nt). Previously we had shown potent exon-inclusion regulatory responses during our scanning mutagenesis study in the areas of exon 7 [see Figure 1B in ]. For simplicity, we refer to these RNA elements from this point forwards as ΔL (hnRNP L) and ΔA1 (hnRNP A1) (Figure 5A). As shown in Figure 5B, one of three independent experiments, and summarized in Figure 5C, mock-treated ISRE WT samples but not vector controls nor the ISRE GG->CC mutant produced double positive GFP/mCherry cells (8.9% compared to 3.7% and 4.6%). Consistent with our hypothesis of how hnRNP A1 and hnRNP L coordinate spliceosome functioning to generate the CEACAM1-S isoform, WT ISRE coupled-exon 7 ΔLΔA1 double mutants revert double positive GFP/mCherry cells towards WT levels (7.0% compared to 8.9%). Most importantly, the ISRE GG->CC-6-7-8 (ΔLΔA1) mutant creates a striking interaction that decreases double positive GFP/mCherry cells compared to WT ISRE-6-7-8 (2.6% compared to 8.9%; p < 0.01). When IRF-1 is reintroduced to these double mutant promoter-exon 7 cells, we observed an impressive three-fold recovery of double positive GFP/mCherry cells strongly suggesting mutations that disrupt proper folding of exon 7 RNA can be rescued by overexpression of IRF-1. Taken together, these data support our hypothesis that interactions between the spliceosome and transcriptosome occur through IRF-1 (9.7% compared to 2.6%; p < 0.01) and a disruption of either complex dramatically alters AS.
We next performed a series of RNA-protein binding experiments to probe whether the mechanism of IRF-1 associated CEACAM1 binding occurs through direct or indirect interactions. We hypothesized that slight changes in gene expression of hnRNP L and hnRNP A1 by IRF-1 might contribute to altered AS of CEACAM1. Ad-IRF-1 treated cells at 24 h were subjected to RIP assays, which can be used to identify specific RNA molecules associated with specific nuclear or cytoplasmic binding proteins in endogenously formed complexes. Using antibodies specific to hnRNP L, hnRNP M, hnRNP A1 and IRF-1, we tested for the ability to pull-down CEACAM1 transcripts that included or lacked exon 7 (Figure 6). We observed Ad-IRF-1 treatment led to a greater than 45-fold and 66-fold enrichment of CEACAM1-L mRNA by antibodies to hnRNP L and hnRNP A1, respectively versus Ad-null alone (Figure 6A). By contrast, we observed negligible association of CEACAM1-L mRNA by antibodies to hnRNP M. Similar trends were observed for CEACAM1-S mRNA, though more enrichment of short isoform mRNA was pulled down by antibodies to hnRNP L and hnRNP A1 (Figure 6B). This phenotype is reasonable given the translational induction of CEACAM1-S exceeds that compared to CEACAM1-L (see Figure 1D, lane 4 versus 3). hnRNP M did not strongly associate with CEACAM1 mRNA (Figures 6A-B), as expected since our RNA co-immunoprecipitation experiments never showed a direct interaction of hnRNP M with exon 7 RNA . This supports the hypothesis that hnRNP M transiently associates with the spliceosomal complex that regulates exon 7 AS (Figure 1A). Taken together, these results also validate our previous finding that hnRNP A1 influences exon 7 despite not having a typical consensus sequence found in several eukaryotic genes suggesting that RNA species require interrogation beyond simple identification of local splicing factors and sequences . Since we could not detect precipitation of CEACAM1 mRNA using IRF-1, we conclude that the effect of the CEACAM1 promoter on AS of CEACAM1 occurs through indirect interactions (Figure 6C).
IRF-1 mediates Global AS
We next asked whether this ability of IRF-1 to participate in splice-site selection is representative of a generalized pathway of induction of molecular networks conferring putative host defense to the cell. We evaluated gene expression and AS levels from RNA-seq data following a workflow that was designed to identify skipped exon events in the transcriptome (Figure 7A, Additional file 5: Table S1). We evaluated our dataset for genes that were up (2.3%) or down-regulated (1.7%) and that had exons exclusively expressed (1.4%) in only the Ad-IRF-1 treatment (Figure 7B). When we focused our attention on genes that were either up-regulated with exons that were exclusively expressed or down-regulated with exons exclusively expressed, we observed subsets that included similar numbers of genes (12 each; Figure 7B). To provide evidence for the inclusion of the exclusively expressed exon, we performed quantitative reverse transcription PCR (using primer pairs as shown next to each gene studied; Figure 7C). The relative abundance of IRF-1-dependent transcripts were characterized for genes DYSF, MYOF, DDC, MAP2 and ST7L in comparison to housekeeping β-actin and a positive control of CEACAM1-L mRNA using isoform-specific primers (Figure 7C). We observed an approximately 5-fold (DDC; p < 0.001) to 20-fold (MAP2; p < 0.001) increase in the number of exons exclusively expressed in transcripts induced in IRF-1 treated samples as compared to the control Ad-null. Based on our customized bioinformatics approach and the small number of genes tested as a proof-in-principle, we predict many more genes may have their AS influenced by IRF-1 (Additional file 5: Table S1). We conclude that besides its role as a transcription factor, IRF-1 plays a significant role in the regulation of AS of human genes involved in a variety of cellular functions. This newly discovered role may augment the immune response by global changes in AS.
Since IRF-1 is regulated by IFNs, and since cytokine induction is associated with an increase in alternatively spliced mRNA transcripts [36, 37], we next considered the possibility that cytokines that undergo AS are directly mediated by the expression of IRF-1. Using a human cytokine array, mRNA transcripts were analyzed by quantitative real-time PCR for cytokines representing Type I (IFN-α genes and IFN-β), Type II (IFN-γ), Type III (interleukins) and TNF superfamily (LTA) genes in comparison to housekeeping genes, GUSB and HPRT (Figure 8). We observed uniform and statistically significant transcriptional differences between untreated and treated Ad-IRF-1 infected cells for the Type I interferons, as expected. Most notably, this included IFNA1 (p < 0.05), IFNB1, IFNA2, IFNA6, IFNA7 (all p < 0.01) and IFNA8, IFNA16, IFNA17 (all p < 0.001). We also observed highly responsive changes in expression in the interleukins, most notably IL1β, IL3, IL4, IL12B, IL13 and IL15 which showed several orders of magnitude difference (all p < 0.001). Taken together, these data extend the number of genes that are under the control of AS master regulator, IRF-1.
Since the CEACAM1 long isoform contains an ITIM, a characteristic motif found in immune cells, we next addressed the potential impact of the isoform switch in the connection to tumorgenesis where an imflammatory environment is considered pro-tumorigenic. For this study, we chose to continue with our model breast cancer cell line MDA-MB-468 and asked how the CEACAM1 isoform ratio would behave in a tumor. MDA-MB-468 cells were orthotopically implanted into mammary fat pads of NOD/SCID mice (n = 6) and tumor volumes were measured until they reached an average size of 150 mm3 (Figure 9A). Total RNAs were isolated from 6 tumors and two parental cell populations and subjected to quantitative PCR using human exon-junction specific primers for CEACAM1. Our data shows that the cytoplasmic domain of CEACAM1, produced primarily as a short isoform in the parental cells grown on plastic, undergoes a significant 4-fold induction (12.03/3.15) of AS to generate the cytoplasmic L-isoform (Figure 9B). Immunofluorescence staining of human CEACAM1 expressing cells from MDA-MB-468 tumor tissue using antibodies directed to the ectodomain shown in red (5F4) or the -L cytoplasmic tail of CEACAM1 (229) shown in green, showed strong expression along the invasive front of the breast tumor (Figure 9C; Upper) as compared to control secondary antibodies (Figure 9C; Lower). When we investigated whether invasive CEACAM1 positive mouse cells on the tumor front might be responsible for the AS we observed for human CEACAM1 MDA-MB-468 cells (Additional file 6: Figure S5), we observed little staining with antibody CC1, which recognizes mouse CEACAM1 exclusively. Taken together, our results with human CEACAM1 are consistent with recent findings that showed increased expression of the L-isoform associates with the invasive front in colorectal cancer  and that the S/L isoform ratio of CEACAM1 is altered in breast cancer specimens . We conclude that differences between long and short cytoplasmic isoforms occur in human breast tumor CEACAM1 expressing cells, and are influenced by the tumor microenvironment, and specifically, favor the long isoform characteristic of immune cells rather than the normal epithelial short isoform.
CEACAM1 AS is dependent upon IRF-1 activity
In this study, we have shown that promoter regulation by IRF-1 determines CEACAM1 AS, specifically to generate the cytoplasmic L-isoform. Previously we proposed that regulatory cis-acting elements in exon 7 in CEACAM1 in combination with the relative abundance of splicing factors hnRNP L, hnRNP A1, and hnRNP M play the critical role in the AS of CEACAM1 (Figure 1; ). Here we extend our findings to include functional interactions between the promoter of CEACAM1 and exon 7 RNA that play a role in isoform-specific splicing.
The critical observation that led to this study was that type II IFN-γ led to the AS of CEACAM1 in breast carcinoma cells, MDA-MB-468 . Previously, it was reported that IFN-γ can cause the AS of tryptophanyl-tRNA synthetase, though the direct mechanism causing this splicing switch was not explored . IFN-γ is an essential mediator of the immune response to viruses and cancer and is one of the strongest inducers of IRF-1 expression in cancer cells. While IRF-1 is itself subject to AS [40, 41], this is the first study to show that IRF-1 coordinates promoter activity and AS to control global changes in the cell. Previous studies have shown that promoter identity influences alternative RNA-processing decisions, through kinetic Pol II elongation and chromatin structure . For example, the catalytic subunit Brm of the SWI/SNF chromatin-remodeling complex was shown to favor inclusion of variant exons in the mRNA of several genes, including E-cadherin, BIM, cyclin D1 and CD44 . CoAA, an hnRNP-like protein, was shown to mediate transcriptional and splicing effects through interactions with transcriptional coregulator TRBP, a protein recruited to target promoters through interactions with activated nuclear receptors . From the point of view of the spliceosome, there are few studies that show how splicing factors mechanistically control transcription.
The novelty of this work lies in the following. First, we show that IRF-1 can control AS of CEACAM1 and a number of other genes involved in cellular maintenance and differentiation (Figures 1, 7, 8). Second, CAM1 6-7-8 minigenes were genetically fused to CEACAM1 promoter and GFP, where only the L-isoform indicates an authentic AS event. Third, to prove that pre-mRNA processing events were temporally coupled to IRF-1, we show that mutations in the ISRE that caused a loss-of-function phenotype in promoter-only studies (Figure 2) in Ad-null/MDA-MB-468 transfected cells were no longer effective (Figure 4; see explanation to follow). When Ad-IRF-1/MDA-MB-468 transfected cells were analyzed, we observed 3-fold increases in GFP expression in our mutant constructs strongly validating a role for IRF-1 in the AS of CEACAM1. Fourth, to prove that compensation of promoter activity was coming from the CEACAM1 minigene RNA, we created double mutations in the promoter and exon 7 RNA and showed pre-mRNA processing events were temporally coupled (Figure 5). Indeed, we show that cellular-driven AS fates in CEACAM1 expression of the L-cytoplasmic tail can secondarily co-regulate its own gene transcription, in a mechanism similar to how Drosophila melanogaster embryonic transcription units show intron looping and associated ribonucleoprotein complexes on transcripts tethered to DNA [45, 46] (Figure 6C). Our data suggests that in breast carcinoma epithelial cells the default splicing maintains production of CEACAM1-S. When IRF-1 binds to its cognate ISRE sequence, the transcriptional complex communicates with the splicing machinery via uncharacterized interactions through splicing regulators leading to an AS switch and production of CEACAM1-L. Fifth, we infer that SP1 and USF1, known transcription factors associated with the CEACAM1 promoter , do not rescue loss-of-function mutations in Ad-null/MDA-MB-468 transfected cells with our CEACAM1 ISRE:CAM1 6-7-8 reporter constructs (Figures 2, 4-5) strongly suggesting that IRF-1 is unique among the transcription factors for generating the CEACAM1 long ITIM containing isoform.
A Global role for IRF-1 and AS
While global characterization of IRF-1 genes in vertebrates as a mechanism for host-defense has been studied extensively , the study of how IFN-γ via IRF-1 increases the complexity and diversity of the immune system through AS is rather limited. Host defense serves two main functions: the generation of immune responses to invading pathogens and the suppression of tumor development. These are achieved by the efficient coordination of prompt and controlled cellular responses. Additionally, tight control of AS can be influenced by external conditions (e.g., the circadian clock , nutritional state , and cell-cycle controls , for example) to function as specific regulators of gene expression. In this respect, our work adds to the growing body of evidence that links AS to environmental cues that activate the spliceosome, in this case immune activation through IRF-1. Our RNA-sequencing and cytokine studies (Figures 7, 8) provide new evidence for a function beyond the known IRF-1 effector responses by showing that multiple genes involved in a wide-array of cellular functions are subject to IRF-1 AS (Figure 7; Additional file 5: Table S1).
In vivo generation of an ITIM-containing isoform
Breast cancers are phenotypically and genotypically heterogeneous tumors containing multiple cancer cell populations with various metastatic potential. In normal breast epithelial cells, where CEACAM1 is predominantly expressed as the S-isoform, but the S/L ratio is altered in breast cancer , one can ask if this change is a consequence of the inflammatory tumor microenvironment, and furthermore, if this change is advantageous to tumor growth and/or invasion. First, it is noteworthy that immune cells express the long and not the short isoform, perhaps due to their own response to inflammatory signals that include IFN-γ produced by activated T- and NK cells . Second, the change in microenvironment during tumorigenesis may cause this change in AS, again via IFN-γ and IRF-1. This phenotype has been reported previously in colon cancer, where the long cytoplasmic dominance of CEACAM1 was associated with invasion and migration of colorectal cancer cells . Indeed, in our study we demonstrate that a similar AS change of CEACAM1 was observed at the invasion front of breast cancer cells (Figure 9). This finding is consistent with recent studies that delineate the role of the microenvironment in determining tissue-specific CEACAM1 S or L dominance in immune cells  and in creating a pro-angiogenic tumor microenvironment . Finally, it is worth considering if expression of an inhibitory ITIM containing isoform is pro-tumorigenic. The answer to this question will require further delineation of the signaling pathways in these cells, a complex issue to say the least.
Taken together, our results reveal a novel role for IRF-1 in the cotranscriptional processing of RNA. In this study, we have shown that promoter regulation by IRF-1 determines CEACAM1 AS, specifically to generate the cytoplasmic L-isoform. Effector responses by IRF-1 in the context of AS should reveal novel cellular functions. Our work is consistent with the idea that immune cells reacting to invading areas of a growing tumor generate IRF-1 in the tumor perhaps through release of IFN-γ, which eventually leads to AS to generate CEACAM1 L-isoform expression not only on T cells but in non-inflamed epithelial cells as well. In cases of acute infection this sort of temporary change in phenotype may be advantageous, while a chronic pro-inflammatory phenotype associated with tumor growth may favor further tumor growth.
Cell culture, reagents, and treatments
MDAMB468 (ATCC® HTB132™) and MCF7 (ATCC® HTB-22) cells were purchased from ATCC and grown in a 5% CO2 incubator at 37°C in MEM supplemented with 1% Sodium Pyruvate, 1.5% Sodium Bicarbonate, 1 × Non-essential Amino Acids, 1 × Penicillin/Streptomycin/Amphotericin B and 10% heat-inactivated FBS. For induction of IRF-1, mouse IRF-1 was used essentially as described previously except for the following modifications . Cells were seeded at a density of 0.5 × 106 cells in 6-well plates 24 h prior to treatment. Recombinant mouse IRF-1 adenovirus (hereafter referred to as Ad-IRF-1) or the Ad-null empty vector control was added to cultures at PFU estimated at 100 particles per PFU. Cells were incubated for 24 h at 37°C and 5% CO2 unless otherwise specified. After incubation, RNA and proteins were isolated as described previously . For Western blots, FACS and immunofluorescent staining identifying CEACAM1 isoforms anti-CEACAM1-L (229; ), T84.1 and 5F4 were used. For CEACAM1 specific antibody, 5F4 was a kind gift of R. Blumberg (Harvard University). Antibodies used for Western blots identifying IRF-1 (ab26109), IRF-2 (C-19, sc-498), GAPDH (FL-335, sc-25778) and p21 (EA10, ab16767) were purchased from Santa Cruz Biotechnology or from Abcam.
RNA isolation semi and quantitative RT-PCR
Total RNA was isolated by the RNeasy mini kit (Qiagen). The RNA was treated with RNase-free DNase set (Qiagen), and RNA (1 μg) was reverse-transcribed in a 20 μl reaction using Maxima First Strand cDNA Synthesis Kit for RT-PCR (K1641, Thermo Scientific) following manufacturer instructions. Reactions were diluted 5-fold and 1/50 of cDNA were used for semiquantitative PCR with gene-specific primers and 2x Taq PCR Premix (Bioland Scientific) for 30 cycles for standard RT-PCR. Primers designed to amplify CEACAM1-3 and CEACAM1-4 isoform were described previously . Amplification of CEACAM1 mRNA for RT-PCR was performed by using sense exon 6 primers (5’-TTCTGCATTTCGGGAAGACCGGCAG-3’) coupled with antisense 3’ UTR primers (5’-AGCCTGGAGATGCCTATTAG-3’). Primer amplification regions for GFP are indicated in Figure 2C. The products were resolved on 2.25% agarose gels and visualized by staining with SYBR Green I (Invitrogen). Gels were photographed on a GelLogic 200 Imaging System. For real time PCR, 1/50 of cDNA reverse transcription reactions were used to amplify constitutively expressed housekeeping gene β-Actin and CEACAM1 mRNA levels. The reactions were performed on a CFX96 Touch Real-Time PCR Detection System (Bio-Rad) in a 20 μl volume with SsoFast EvaGreen Supermix (Bio-Rad). Primers BP56 (5’- GGTTGCTCTGATAGCAGTAGCCCTGGCATGTTTTCTGCATTTCGGGAAGACCGGCAGCTC-3’) and BP60 with a common sense primer to exon 6 amplify differentially spliced isoforms, CEACAM1-S and CEACAM1-L, respectively . TaqMan Array Human Cytokine Network (Life Technologies) plates were carried out according to manufacturer instructions. After denaturation for 30 sec at 95°C, 40 cycles were performed (5 sec at 95°C, 5 sec at 60°C). At the end of the reaction, melting curves were generated between 65°C and 95°C, for every 0.5°C. CEACAM1 mRNA and genes identified in our RNA Seq study were calculated by using β-Actin for normalisation. Quantification of gene expression was calculated by using the ΔCT method .
Promoter and splicing reporters
All DNA constructs were generated by standard cloning procedures and confirmed by sequencing. Construction of promoter CEACAM1 WT ISRE:GFP plasmid proceeded by amplifying a genomic DNA fragment corresponding to -1,135 bp relative to the transcription start site (TSS) from MDA-MB-468 cells using Bgl II and Hnd III restriction sites cloned into pZsGreen-1 (Clontech). To construct CEACAM1 WT ISRE:CAM1 6-7-8:GFP plasmid, the Bam HI restriction site in the CEACAM1 promoter was disrupted to allow ligation of our minigene splicing cassette containing genomic CEACAM1 exons 6 through exon 8, described previously , using Sac II and Bam HI restriction sites. AS events that favored only L-isoform led to the excitation of GFP fluorescence. To control for external validity and allow for quantitation of transfection efficiency, promoter and splicing reporter experiments were conducted using a second co-transfected promoter CEACAM1 WT ISRE:mCherry plasmid. Mutations for all constructs were introduced using QuikChange Lightning Site-Directed Mutagenesis Kit (Agilent Technologies) following manufacturer’s instructions. To delete both hnRNP L and hnRNP A1 binding sites in exon 7, the sequence AAGCATGCAA was introduced using GeneArt mutagenesis (Life Technologies). This sequence does not associate with any known cellular factor and therefore does not influence AS .
Fluorescence-based electrophoretic mobility shift assay (fEMSA)
The fEMSA was operated using a Light Shift Chemiluminescent EMSA kit (Thermo Scientific) according to the manufacturer’s instructions with the only exception that IRDye fluorescence-labeled oligonucleotides for ISRE WT or ISRE GG->CC were used to visualize DNA-protein binding interactions as described previously . EBNA was used as a positive control for the assay reagents. Sequence information for each set of oligos is contained in Figure 3A. The DNA-protein binding reaction buffer was composed of 1X binding buffer supplemented with 5% Glycerol, 10 mM MgCl2, 200 ng Poly (dI-dC) final concentration and 0.3 pmol IRDye-labeled oligonucleotides in a total volume of 10 μl. Recombinant IRF-1 and IRF-2 (ProSpec) each containing a His-tag were used in binding reactions and the reactions were incubated at room temperature for 30 min before separation of complexes on a 6% DNA Retardation Gel (Invitrogen) at 100 V for 60 min. Detection of DNA binding was performed using Licor Odyssey infrared scanning using the following parameters: 700 and 800 nm channels, focus height 3.5 mm, and intensity set to 5.
Flow cytometric analysis
Promoter and splicing reporters were transfected using cells at a density of 1 × 106 cells using Amaxa Cell Line Nucleofector Kit V (Lonza) and 6 μg DNA from each construct. To pre-heated 6-well plates containing growth media, cells were electroporated and immediately incubated for 24 h before flow cytometric analyses using a BD Fortessa (Becton Dickinson, San Jose, CA). GFP and mCherry fluorescence intensities were detected with 488 nm and 561 nm lasers, respectively. Data were analyzed with FlowJo Software version 8.1.1 (Treestar).
We performed RNA immunoprecipitation (RIP) experiments using the Magna RIP RNA-Binding Protein Immunoprecipitation Kit (Millipore, Bedford, MA), according to the manufacturer’s instructions. Antibodies to hnRNP L (4D11, sc-32317), hnRNP A1 (4B10, sc-32301), and hnRNP M (3C181, sc-56702) were purchased from Santa Cruz Biotechnology. The co-precipitated RNAs were detected by RT-PCR. Total RNAs (input controls) and isotype controls were assayed simultaneously to demonstrate that the detected signals were from RNAs specifically binding to CEACAM1 mRNA (data not shown). snRNP 70 antibody was used as a negative control for CEACAM1 binding.
Sequences obtained from Illumina’s CASAVA pipeline were aligned to human genome assembly hg19 using Tophat v1.2.0 . The expression levels of Refseq transcripts were calculated by the number of reads falling into the exon region of each transcript. The expression data were then quantile normalized and log2 transformed. Differentially expressed transcripts were identified by the following criteria: 1) Chi-square p value < 0.05; 2) Log2 ratio > 1 or < -1; 3) at least one sample has coverage > 10, where the coverage was calculated by the expression level of the transcript divided by the transcript length. The up-regulated and down-regulated genes were then subject to functional annotation analysis with DAVID (http://david.abcc.ncifcrf.gov/). A customized bioinformatics pipeline was developed to detect exon skipping event in the transcriptome. Briefly, for each Refseq transcript with 3 or more exons, the middle exons excluding the first and last exon were kept for further analysis. When exons were shared by multiple transcripts, only one exon was kept. We defined inclusion reads and exclusion reads as described previously . The inclusion and exclusion reads were counted in both control and IRF-1 treated samples. The inclusion/exclusion ratio for each sample was then calculated. The potential skipped exons in control sample were identified by the following criteria: 1) inclusion reads in control = 0; 2) exclusion reads in control > 1; 3) inclusion reads in Ad-null > =5; 4) exclusion reads in Ad-null =0; 5) Difference of the inclusion/exclusion ratio Ad-null vs. control sample > 8-fold. Similar approach was used to identify skipped exons in the Ad-IRF-1 sample. These combinations of criteria were used to reduce the chance of false positives. The identified exons were then annotated back to Refseq transcripts and compared to differentially expressed transcripts identified above.
Tumor xenografts in NOD/SCID mice
NOD/SCID transgenic mice were purchased from the Jackson Laboratory. MDA-MB-468 tumor cells (1 × 106) were mixed with Matrigel (BD Biosciences) and orthotopically injected into mammary fat pads. The experiment was terminated after 70 days when tumors reached an average size of 150 mm3. Mouse care and experimental procedures were performed under pathogen-free conditions in accordance with established institutional guidance and approved protocols from the Institutional Animal Care and Use Committee of the Beckman Research Institute of City of Hope.
Frozen tumor tissue sections were fixated in 2% paraformaldehyde and permeabilized in methanol. After blocking first using Mouse on Mouse (M.O.M.) Blocking Reagent (Vector Lab) and then in PBS containing 10% goat serum for 1 h each, slides were incubated overnight with antibodies: mouse anti-human CEACAM1 ectodomain (5F4, 1:200), rabbit anti-human -L cytoplasmic tail of CEACAM1 (229, 1:200)  or mouse anti-mouse CC1, a kind gift of K. Holmes (University of Colorado). The following day the slides were washed and secondary antibodies were applied for 1 h (goat anti-rabbit, Alexa Fluor 488 labeled and goat anti-mouse, Alexa Fluor 555, both 1:200, Invitrogen) were visualized using an Olympus IX81 Inverted microscope with Q Imaging Retiga 2000R cooled CCD camera with Image Pro Plus version 6.3 imaging software and 40X magnification. For nuclei staining, Hoechst 33342 (blue) was used at final concentration of 1 μg/ml.
Bonjardim CA: Interferons (IFNs) are key cytokines in both innate and adaptive antiviral immune responses–and viruses counteract IFN action. Microbes Infect. 2005, 7: 569-578. 10.1016/j.micinf.2005.02.001
Li J, Yu B, Song L, Eschrich S, Haura EB: Effects of IFN-gamma and Stat1 on gene expression, growth, and survival in non-small cell lung cancer cells. J Interferon Cytokine Res. 2007, 27: 209-220. 10.1089/jir.2006.0111
Schoenborn JR, Wilson CB: Regulation of interferon-gamma during innate and adaptive immune responses. Adv Immunol. 2007, 96: 41-101.
Savitsky D, Tamura T, Yanai H, Taniguchi T: Regulation of immunity and oncogenesis by the IRF transcription factor family. Cancer Immunol Immunother. 2010, 59: 489-510.
Escalante CR, Yie J, Thanos D, Aggarwal AK: Structure of IRF-1 with bound DNA reveals determinants of interferon regulation. Nature. 1998, 391: 103-106. 10.1038/34224
Chen CJ, Lin TT, Shively JE: Role of interferon regulatory factor-1 in the induction of biliary glycoprotein (cell CAM-1) by interferon-gamma. J Biol Chem. 1996, 271: 28181-28188. 10.1074/jbc.271.45.28181
Takahashi H, Okai Y, Paxton RJ, Hefta LJ, Shively JE: Differential regulation of carcinoembryonic antigen and biliary glycoprotein by gamma-interferon. Cancer Res. 1993, 53: 1612-1619.
Gencheva M, Chen CJ, Nguyen T, Shively JE: Regulation of CEACAM1 transcription in human breast epithelial cells. BMC Mol Biol. 2010, 11: 79- 10.1186/1471-2199-11-79
Kuespert K, Pils S, Hauck CR: CEACAMs: their role in physiology and pathophysiology. Curr Opin Cell Biol. 2006, 18: 565-571. 10.1016/j.ceb.2006.08.008
Lu R, Pan H, Shively JE: CEACAM1 negatively regulates IL-1beta production in LPS activated neutrophils by recruiting SHP-1 to a SYK-TLR4-CEACAM1 complex. PLoS Pathog. 2012, 8: e1002597- 10.1371/journal.ppat.1002597
Dery KJ, Gaur S, Gencheva M, Yen Y, Shively JE, Gaur RK: Mechanistic control of carcinoembryonic antigen-related cell adhesion molecule-1 (CEACAM1) splice isoforms by the heterogeneous nuclear ribonuclear proteins hnRNP L, hnRNP A1, and hnRNP M. J Biol Chem. 2011, 286: 16039-16051. 10.1074/jbc.M110.204057
Najjar SM, Accili D, Philippe N, Jernberg J, Margolis R: Taylor SI: pp 120/ecto-ATPase, an endogenous substrate of the insulin receptor tyrosine kinase, is expressed as two variably spliced isoforms. J Biol Chem. 1993, 268: 1201-1206.
Najjar SM, Boisclair YR, Nabih ZT, Philippe N, Imai Y, Suzuki Y, Suh DS, Ooi GT: Cloning and characterization of a functional promoter of the rat pp 120 gene, encoding a substrate of the insulin receptor tyrosine kinase. J Biol Chem. 1996, 271: 8809-8817. 10.1074/jbc.271.15.8809
Kirshner J, Chen CJ, Liu P, Huang J, Shively JE: CEACAM1-4S, a cell-cell adhesion molecule, mediates apoptosis and reverts mammary carcinoma cells to a normal morphogenic phenotype in a 3D culture. Proc Natl Acad Sci U S A. 2003, 100: 521-526. 10.1073/pnas.232711199
Kirshner J, Schumann D, Shively JE: CEACAM1, a cell-cell adhesion molecule, directly associates with annexin II in a three-dimensional model of mammary morphogenesis. J Biol Chem. 2003, 278: 50338-50345. 10.1074/jbc.M309115200
Lu R, Niesen MJ, Hu W, Vaidehi N, Shively JE: Interaction of actin with carcinoembryonic antigen-related cell adhesion molecule 1 (CEACAM1) receptor in liposomes is Ca2 + - and phospholipid-dependent. J Biol Chem. 2011, 286: 27528-27536. 10.1074/jbc.M111.235762
Gaur S, Shively JE, Yen Y, Gaur RK: Altered splicing of CEACAM1 in breast cancer: identification of regulatory sequences that control splicing of CEACAM1 into long or short cytoplasmic domain isoforms. Mol Cancer. 2008, 7: 46- 10.1186/1476-4598-7-46
Huber M, Izzi L, Grondin P, Houde C, Kunath T, Veillette A, Beauchemin N: The carboxyl-terminal region of biliary glycoprotein controls its tyrosine phosphorylation and association with protein-tyrosine phosphatases SHP-1 and SHP-2 in epithelial cells. J Biol Chem. 1999, 274: 335-344. 10.1074/jbc.274.1.335
Kammerer R, Hahn S, Singer BB, Luo JS, von Kleist S: Biliary glycoprotein (CD66a), a cell adhesion molecule of the immunoglobulin superfamily, on human lymphocytes: structure, expression and involvement in T cell activation. Eur J Immunol. 1998, 28: 3664-3674. 10.1002/(SICI)1521-4141(199811)28:11<3664::AID-IMMU3664>3.0.CO;2-D
Morales VM, Christ A, Watt SM, Kim HS, Johnson KW, Utku N, Texieira AM, Mizoguchi A, Mizoguchi E, Russell GJ, Russell SE, Bhan AK, Freeman GJ, Blumberg RS: Regulation of human intestinal intraepithelial lymphocyte cytolytic function by biliary glycoprotein (CD66a). J Immunol. 1999, 163: 1363-1370.
Luo W, Wood CG, Earley K, Hung MC, Lin SH: Suppression of tumorigenicity of breast cancer cells by an epithelial cell adhesion molecule (C-CAM1): the adhesion and growth suppression are mediated by different domains. Oncogene. 1997, 14: 1697-1704. 10.1038/sj.onc.1200999
Hsieh JT, Luo W, Song W, Wang Y, Kleinerman DI, Van NT, Lin SH: Tumor suppressive role of an androgen-regulated epithelial cell adhesion molecule (C-CAM) in prostate carcinoma cell revealed by sense and antisense approaches. Canc Res. 1995, 55: 190-197.
Rosenberg M, Nedellec P, Jothy S, Fleiszer D, Turbide C, Beauchemin N: The expression of mouse biliary glycoprotein, a carcinoembryonic antigen-related gene, is down-regulated in malignant mouse tissues. Canc Res. 1993, 53: 4938-4945.
Samineni S, Zhang Z, Shively JE: Carcinoembryonic antigen-related cell adhesion molecule 1 negatively regulates granulocyte colony-stimulating factor production by breast tumor-associated macrophages that mediate tumor angiogenesis. Int J Canc. 2013, 133: 394-407. 10.1002/ijc.28036.
Black DL: Mechanisms of alternative pre-messenger RNA splicing. Annu Rev Biochem. 2003, 72: 291-336. 10.1146/annurev.biochem.72.121801.161720
Dery KJ, Gusti V, Gaur S, Shively JE, Yen Y, Gaur RK: Alternative splicing as a therapeutic target for human diseases. Methods Mol Biol. 2009, 555: 127-144. 10.1007/978-1-60327-295-7_10
Jean-Philippe J, Paz S, Caputi M: hnRNP A1: The Swiss Army Knife of Gene Expression. Int J Mol Sci. 2013, 14: 18999-19024. 10.3390/ijms140918999
Lin S, Fu XD: SR proteins and related factors in alternative splicing. Adv Exp Med Biol. 2007, 623: 107-122. 10.1007/978-0-387-77374-2_7
Pizzoferrato E, Liu Y, Gambotto A, Armstrong MJ, Stang MT, Gooding WE, Alber SM, Shand SH, Watkins SC, Storkus WJ, Yim JH: Ectopic expression of interferon regulatory factor-1 promotes human breast cancer cell death and results in reduced expression of survivin. Canc Res. 2004, 64: 8381-8388. 10.1158/0008-5472.CAN-04-2223.
Yim JH, Wu SJ, Casey MJ, Norton JA, Doherty GM: IFN regulatory factor-1 gene transfer into an aggressive, nonimmunogenic sarcoma suppresses the malignant phenotype and enhances immunogenicity in syngeneic mice. J Immunol. 1997, 158: 1284-1292.
Harada H, Willison K, Sakakibara J, Miyamoto M, Fujita T, Taniguchi T: Absence of the type I IFN system in EC cells: transcriptional activator (IRF-1) and repressor (IRF-2) genes are developmentally regulated. Cell. 1990, 63: 303-312. 10.1016/0092-8674(90)90163-9
Armstrong MJ, Stang MT, Liu Y, Gao J, Ren B, Zuckerbraun BS, Mahidhara RS, Xing Q, Pizzoferrato E, Yim JH: Interferon Regulatory Factor 1 (IRF-1) induces p21(WAF1/CIP1) dependent cell cycle arrest and p21(WAF1/CIP1) independent modulation of survivin in cancer cells. Canc Lett. 2012, 319: 56-65. 10.1016/j.canlet.2011.12.027.
Newman EA, Muh SJ, Hovhannisyan RH, Warzecha CC, Jones RB, McKeehan WL, Carstens RP: Identification of RNA-binding proteins that regulate FGFR2 splicing through the use of sensitive and specific dual color fluorescence minigene assays. RNA. 2006, 12: 1129-1141. 10.1261/rna.34906
Harada H, Fujita T, Miyamoto M, Kimura Y, Maruyama M, Furia A, Miyata T, Taniguchi T: Structurally similar but functionally distinct factors, IRF-1 and IRF-2, bind to the same regulatory elements of IFN and IFN-inducible genes. Cell. 1989, 58: 729-739. 10.1016/0092-8674(89)90107-4
Huang Y, Li W, Yao X, Lin QJ, Yin JW, Liang Y, Heiner M, Tian B, Hui J, Wang G: Mediator complex regulates alternative mRNA processing via the MED23 subunit. Mol Cell. 2012, 45: 459-469. 10.1016/j.molcel.2011.12.022
Eissa NT, Strauss AJ, Haggerty CM, Choo EK, Chu SC, Moss J: Alternative splicing of human inducible nitric-oxide synthase mRNA. tissue-specific regulation and induction by cytokines. J Biol Chem. 1996, 271: 27184-27187. 10.1074/jbc.271.43.27184
Sellers RS, Luchin AI, Richard V, Brena RM, Lima D, Rosol TJ: Alternative splicing of parathyroid hormone-related protein mRNA: expression and stability. J Mol Endocrinol. 2004, 33: 227-241. 10.1677/jme.0.0330227
Ieda J, Yokoyama S, Tamura K, Takifuji K, Hotta T, Matsuda K, Oku Y, Nasu T, Kiriyama S, Yamamoto N, Nakamura Y, Shively JE, Yamaue H: Re-expression of CEACAM1 long cytoplasmic domain isoform is associated with invasion and migration of colorectal cancer. Int J Cancer. 2011, 129: 1351-1361. 10.1002/ijc.26072
Liu J, Shue E, Ewalt KL, Schimmel P: A new gamma-interferon-inducible promoter and splice variants of an anti-angiogenic human tRNA synthetase. Nucleic Acids Res. 2004, 32: 719-727. 10.1093/nar/gkh240
Tzoanopoulos D, Speletas M, Arvanitidis K, Veiopoulou C, Kyriaki S, Thyphronitis G, Sideras P, Kartalis G, Ritis K: Low expression of interferon regulatory factor-1 and identification of novel exons skipping in patients with chronic myeloid leukaemia. Br J Haematol. 2002, 119: 46-53. 10.1046/j.1365-2141.2002.03829.x
Lee EJ, Jo M, Park J, Zhang W, Lee JH: Alternative splicing variants of IRF-1 lacking exons 7, 8, and 9 in cervical cancer. Biochem Biophys Res Commun. 2006, 347: 882-888. 10.1016/j.bbrc.2006.06.145
Shukla S, Oberdoerffer S: Co-transcriptional regulation of alternative pre-mRNA splicing. Biochim Biophys Acta. 1819, 2012: 673-683.
Batsche E, Yaniv M, Muchardt C: The human SWI/SNF subunit Brm is a regulator of alternative splicing. Nat Struct Mol Biol. 2006, 13: 22-29. 10.1038/nsmb1030
Auboeuf D, Dowhan DH, Li X, Larkin K, Ko L, Berget SM, O’Malley BW: CoAA, a nuclear receptor coactivator protein at the interface of transcriptional coactivation and RNA splicing. Mol Cell Biol. 2004, 24: 442-453. 10.1128/MCB.24.1.442-453.2004
Beyer AL, Bouton AH, Miller OL: Correlation of hnRNP structure and nascent transcript cleavage. Cell. 1981, 26: 155-165. 10.1016/0092-8674(81)90299-3
Beyer AL, Osheim YN: Splice site selection, rate of splicing, and alternative splicing on nascent transcripts. Genes Dev. 1988, 2: 754-765. 10.1101/gad.2.6.754
Moore MJ, Silver PA: Global analysis of mRNA splicing. RNA. 2008, 14: 197-203.
McGlincy NJ, Valomon A, Chesham JE, Maywood ES, Hastings MH, Ule J: Regulation of alternative splicing by the circadian clock and food related cues. Genome Biol. 2012, 13: R54- 10.1186/gb-2012-13-6-r54
Marden JH, Fescemyer HW, Saastamoinen M, MacFarland SP, Vera JC, Frilander MJ, Hanski I: Weight and nutrition affect pre-mRNA splicing of a muscle gene associated with performance, energetics and life history. J Exp Biol. 2008, 211: 3653-3660. 10.1242/jeb.023903
Moore MJ, Wang Q, Kennedy CJ, Silver PA: An alternative splicing network links cell-cycle control to apoptosis. Cell. 2010, 142: 625-636. 10.1016/j.cell.2010.07.019
Nakano M, Saeki C, Takahashi H, Homma S, Tajiri H, Zeniya M: Activated natural killer T cells producing interferon-gamma elicit promoting activity to murine dendritic cell-based autoimmune hepatic inflammation. Clin Exp Immunol. 2012, 170: 274-282. 10.1111/j.1365-2249.2012.04664.x
Chen L, Chen Z, Baker K, Halvorsen EM, da Cunha AP, Flak MB, Gerber G, Huang YH, Hosomi S, Arthur JC, Dery KJ, Nagaishi T, Beauchemin N, Holmes KV, Ho JW, Shively JE, Jobin C, Onderdonk AB, Bry L, Weiner HL, Higgins DE, Blumberg RS: The short isoform of the CEACAM1 receptor in intestinal T cells regulates mucosal immunity and homeostasis via Tfh cell induction. Immunity. 2012, 37: 930-946. 10.1016/j.immuni.2012.07.016
Gerstel D, Wegwitz F, Jannasch K, Ludewig P, Scheike K, Alves F, Beauchemin N, Deppert W, Wagener C, Horst AK: CEACAM1 creates a pro-angiogenic tumor microenvironment that supports tumor vessel maturation. Oncogene. 2011, 30: 4275-4288. 10.1038/onc.2011.146
Chen CJ, Shively JE: The cell-cell adhesion molecule carcinoembryonic antigen-related cellular adhesion molecule 1 inhibits IL-2 production and proliferation in human T cells by association with Src homology protein-1 and down-regulates IL-2 receptor. J Immunol. 2004, 172: 3544-3552.
Yu Q, Chow EM, Wong H, Gu J, Mandelboim O, Gray-Owen SD, Ostrowski MA: CEACAM1 (CD66a) promotes human monocyte survival via a phosphatidylinositol 3-kinase- and AKT-dependent pathway. J Biol Chem. 2006, 281: 39179-39193. 10.1074/jbc.M608864200
Schmittgen TD, Livak KJ: Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc. 2008, 3: 1101-1108. 10.1038/nprot.2008.73
Szeszel-Fedorowicz W, Talukdar I, Griffith BN, Walsh CM, Salati LM: An exonic splicing silencer is involved in the regulated splicing of glucose 6-phosphate dehydrogenase mRNA. J Biol Chem. 2006, 281: 34146-34158. 10.1074/jbc.M603825200
Steiner S, Pfannschmidt T: Fluorescence-based electrophoretic mobility shift assay in the analysis of DNA-binding proteins. Methods Mol Biol. 2009, 479: 273-289. 10.1007/978-1-59745-289-2_18
Trapnell C, Pachter L, Salzberg SL: TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009, 25: 1105-1111. 10.1093/bioinformatics/btp120
Wang ET, Sandberg R, Luo S, Khrebtukova I, Zhang L, Mayr C, Kingsmore SF, Schroth GP, Burge CB: Alternative isoform regulation in human tissue transcriptomes. Nature. 2008, 456: 470-476. 10.1038/nature07509
We thank Lucy Brown in the City of Hope Analytical Cytometry Core and Jennifer Chean for help with flow cytometry experiments and Bing Mu within the City of Hope Bioinformatics core for Illumina base calling process. The authors are grateful to Ren-Jang Lin for critical discussions of the manuscript.
The authors declare that they have no competing financial interests.
Experiments, study design and data analysis were performed by KJD. Animal studies were performed by MK. RIP and splicing assays were performed by DG. XW performed bioinformatics analyses. TN performed western blot and cytokine array analyses. CC helped with transient transfection experiments for reporter assays. JHY provided Ad-IRF-1. The manuscript was written by KJD and edited by JES. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Figure S1: Time-course of MDA-MB-468 breast carcinoma cells after induction with Ad-IRF-1. (A) Western blotting of cell lysates after no treatment and induction of Ad-IRF-1 after 6, 12 and 24 h with antibodies that recognize IRF-1 and GAPDH. (B) Analyses of mRNA expression by RT-PCR (Upper) and histogram representing the quantification of CEACAM1-L/S ratio as determined in Figure 1 (Lower). All assays were performed in triplicate. *, p < 0.05; ***, p < 0.001 versus Ad-null control. (TIFF 382 KB)
Additional file 2: Figure S2: Time-course of MCF7 breast carcinoma cells after induction with Ad-IRF-1. Analyses of CEACAM1–S (A) and CEACAM1–L (B) mRNA expression by qRT-PCR and histogram representing the quantification of CEACAM1 relative to β-Actin mRNA levels. All assays were performed in triplicate. **, p < 0.01; ***, p < 0.001 versus Ad-null control. (TIFF 110 KB)
Additional file 3: Figure S3: Western blotting of cell lysates after induction of Ad-IRF-1 after 24 h. Shown are fluorescent (700 nm and 800 nm) images probed with antibody 229 (A) and antibody T84.1 (B) to detect CEACAM1 isoforms (companion to composite figure shown in Figure 1D). (TIFF 1 MB)
Additional file 4: Figure S4: Flow-cytometry analyses of MDA-MB-468 breast carcinoma cells after induction with Ad-IRF-1. Antibody 5F4, which recognizes both isoforms of CEACAM1, was added to cells that were either mock treated or induced with Ad-IRF-1 after 24 h. Treatments are: Blank, cells not treated with antibody; 2nd only, cells treated with Alexa Fluor 488 goat anti-mouse antibody only; Isotype, cells treated with mouse clone MG1-45 mIgG1 isotype control; MDA-MB-468, cells not treated with virus stimulation; Ad-null, cells treated with virus vector control; Ad-IRF-1, cells treated with virus expressing IRF-1. (TIFF 397 KB)
Additional file 5: Table S1: List of genes differentially expressed in the transcriptome analyses after IRF-1 treatment. (XLS 34 KB)
Additional file 6: Figure S5: Expression pattern of mouse CEACAM1 in vivo. Immunofluorescence staining of mouse CEACAM1 expressing cells from MDA-MB-468 tumor tissue. CEACAM1 expressing cells from MDA-MB-468 tumor tissue were treated with antibodies directed to the ectodomain shown in red (CC1, mouse antibody 1:200; shown with white arrowheads) or the -L cytoplasmic tail of human CEACAM1 by comparison and shown in green (229 rabbit antibody 1:200; shown with white arrows) followed by secondary antibodies (goat anti-rabbit, Alexa Fluor 488 labeled and goat anti-mouse, Alexa Fluor 555 labeled; 1:200). For nuclei staining, Hoechst 33342 (blue) was used at final concentration of 1 μg/ml. Scale bar, 100 μm. (TIFF 1 MB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Dery, K.J., Kujawski, M., Grunert, D. et al. IRF-1 regulates alternative mRNA splicing of carcinoembryonic antigen-related cell adhesion molecule 1 (CEACAM1) in breast epithelial cells generating an immunoreceptor tyrosine-based inhibition motif (ITIM) containing isoform. Mol Cancer 13, 64 (2014). https://doi.org/10.1186/1476-4598-13-64