- Open Access
Discovery, characterization and potential roles of a novel NF-YAx splice variant in human neuroblastoma
Journal of Experimental & Clinical Cancer Research volume 38, Article number: 482 (2019)
Identification of novel cancer-associated splice variants is of potential diagnostic, prognostic and therapeutic importance. NF-Y transcription factor is comprised of NF-YA, NF-YB and NF-YC subunits, binds inverted CCAAT-boxes in ≈70% of gene promoters, regulates > 1000 cancer-associated genes and proteins involved in proliferation, staminality, differentiation, apoptosis, metabolism and is subject to component alternative splicing. RT-PCR evaluation of alternative NF-YA splicing in primary human neuroblastomas (NBs), led to discovery of a novel NF-YAx splice variant, also expressed during mouse embryo development and induced by doxorubicin in NB cells. Here, we report the discovery and characterisation of NF-YAx and discus its potential roles in NB.
NF-YAx cDNA was RT-PCR-cloned from a stage 3 NB (provided by the Italian Association of Haematology and Paediatric Oncology, Genova, IT), sequenced and expressed as a protein using standard methods and compared to known fully-spliced NF-YAl and exon B-skipped NF-YAs isoforms in: EMSAs for capacity to form NF-Y complexes; by co-transfection, co-immunoprecipitation and Western blotting for capacity to bind Sp1; by IF for localisation; in AO/EtBr cell-death and colony formation assays for relative cytotoxicity, and by siRNA knockdown, use of inhibitors and Western blotting for potential mechanisms of action. Stable SH-SY5Y transfectants of all three NF-YA isoforms were also propagated and compared by RT-PCR and Western blotting for differences in cell-death and stem cell (SC)-associated gene expression, in cell-death assays for sensitivity to doxorubicin and in in vitro proliferation, substrate-independent growth and in vivo tumour xenograft assays for differences in growth and tumourigenic capacity.
NF-YAx was characterized as a novel variant with NF-YA exons B, D and partial F skipping, detected in 20% of NF-YA positive NBs, was the exclusive isoform in a stage 3 NB, expressed in mouse stage E11.5–14 embryos and induced by doxorubicin in SH-SY5Y NB cells. The NF-YAx protein exhibited nuclear localisation, competed with other isoforms in CCAAT box-binding NF-Y complexes but, in contrast to other isoforms, did not bind Sp1. NF-YAx expression in neural-related progenitor and NB cells repressed Bmi1 expression, induced KIF1Bβ expression and promoted KIF1Bβ-dependent necroptosis but in NB cells also selected tumourigenic, doxorubicin-resistant, CSC-like sub-populations, resistant to NF-YAx cytotoxicity.
The discovery of NF-YAx in NBs, its expression in mouse embryos and induction by doxorubicin in NB cells, unveils a novel NF-YA splice mechanism and variant, regulated by and involved in development, genotoxic-stress and NB. NF-YAx substitution of other isoforms in NF-Y complexes and loss of capacity to bind Sp1, characterises this novel isoform as a functional modifier of NF-Y and its promotion of KIF1Bβ-dependent neural-lineage progenitor and NB cell necroptosis, association with doxorubicin-induced necroptosis and expression in mouse embryos coinciding with KIF1Bβ-dependent sympathetic neuroblast-culling, confirm a cytotoxic function and potential role in suppressing NB initiation. On the other hand, the in vitro selection of CSC-like NB subpopulations resistant to NF-YAx cytotoxicity not only helps to explain high-level exclusive NF-YAx expression in a stage 3 NB but also supports a role for NF-YAx in disease progression and identifies a potential doxorubicin-inducible mechanism for post-therapeutic relapse.
Alternative gene splicing is a fundamental physiological mechanism for the differential expression of proteins from the same gene coding sequence and is largely responsible for the increased proteomic complexity of higher organisms that cannot be explained by differences in individual gene numbers alone . Aberrant alternative splicing has been reported in cancer and is promoted by stress within the tumour microenvironment, oncogenic viral infection, gene translocation and oncogene de-regulation of splice factor expression. Cancer-associated alternative splicing has been shown to inactivate onco-suppressors and activate oncogenes, making the identification of novel cancer-associated splice isoforms of potential diagnostic, prognostic and therapeutic importance [1,2,3,4,5].
The ubiquitous transcription factor NF-Y binds inverted CCAAT-boxes (5′-ATTGG-3′) in ≈70% of gene promoters, recruits other transcription factors and proteins to promoters and regulates > 1000 cancer-associated genes involved in proliferation, stemness, differentiation, apoptosis, miRNA expression and metabolism [6,7,8,9]. NF-Y also regulates chromatin de-condensation during mitotic bookmarking and post-mitotic transcriptional re-activation, and acts as a bi-directional histone-like transcription factor, switching promoter histone methyl-marks from positive to negative. NF-Y is necessary for development, its inactivation is embryonically lethal and alterations in NF-Y activity have wide-ranging effects on cell behaviour [6,7,8,9,10]. The involvement of NF-Y in cancer is underpinned by its role in proliferation, the presence of inverted CCAAT-boxes in the promoters of many cancer-associated genes and interaction with and regulation of cancer-related proteins, including c-fos, wild type and mutant p53, p73, ΔNp63α, p21, Ash2L, BRCA-1, cMyc, Sox9, lamin, ZHX1/2 and in particular the ubiquitous cancer-related transcription factor Sp1 [6,7,8,9,10,11].
NF-Y is a hetero-trimer composed of NF-YA, NF-YB and NF-YC subunits, all of which are required for DNA binding and transcriptional activity. NF-YB complexes with NF-YC prior to binding NF-YA, which confers DNA binding and transcriptional activity to the NF-Y complex [6,7,8, 12,13,14]. The 29.44 kb NF-YA gene localises to chromosome 6p21, is organized into 9 exons  and is predominantly expressed as a fully-spliced 42 kDa, 347 amino acid (aa) long-form NF-YAl with glutamine-rich, S/T-rich transactivation, subunit-interaction and DNA-binding domains or an alternative exon B-spliced 40 kDa, 318 aa short-form NF-YAs, deleted of glutamine-rich aa’s 26–54 [15, 16]. Minor NF-YA isoforms also include a 3 bp aa 27 deletion-variant, a 18 bp aa 548–565 deletion-variant and L2-L6 alternative splice variants deleted in glutamine-rich and S/T domains, driven by nucleotide 79, 154 and 548 deletions at A/B, B/C and E/F splice junctions [15, 17].
Alternative splicing of the NF-YA gene has been implicated in the regulation of cell staminality, differentiation, apoptosis and transformation. NF-YAs forms part of the stem cell (SC) transcriptional circuitry, predominates in embryonic SCs and is lost upon SC differentiation. In contrast, NF-YAl promotes differentiation and loss of NF-YA expression induces senescence or apoptosis. Alternative NF-YAs splicing is promoted by the oncogenic polyomavirus SV40 and by v-ras oncogene and converts tumor-suppressing, differentiation-promoting NF-Y complexes predominated by NF-YAl into tumor and CSC promoting complexes predominated by NF-YAs [8, 18,19,20,21,22,23].
Neuroblastomas (NB) are aggressive embryonic tumours of neural crest origin, derived from immature sympathetic neuroblasts . These primitive tumours initiate under conditions that impair sympathetic neuroblast culling during development, reported to depend upon either loss of the KIF1B gene associated with chromosome 1p36-deletion, germline KIF1B mutations or Nmyc amplification [25,26,27,28,29,30,31,32,33]. NF-Y involvement in NB pathogenesis and progression, however, has received scant attention. In the few existing reports, NF-Y has been shown to be critical for expression of soluble guanyl cyclase in NB cells required for cGMP production and differentiation  and is involved in elevated glypican 3 expression in NBs . NF-Y and Sp1 transcription factors combine to promote tetramethylpyrazine-induced neuronal differentiation of NB cells  and regulate expression of the α3 Na+, K + -ATPase subunit, essential for maintaining electrochemical gradients across cell membranes . Suboptimal NF-Y function in NB cells has also been implicated in de-regulating the matrix metalloproteinase and tissue inhibitor of metalloproteinase equilibrium, resulting in invasion  and increased expression of the NF-YA subunit has been reported to differentiate between aggressive stage 4 NBs and stage 4S NBs that exhibit spontaneous regression .
Considering the relative absence of studies of NF-Y expression in NB, combined with reports associating fully spliced NF-YAl with cellular differentiation and reduced malignancy and associating alternative exon B spliced NF-YAs with cellular staminality and increased malignancy [18,19,20,21,22,23], we initiated a study of NF-YAl and NF-YAs expression in human primary NBs. This led to the unexpected discovery of a novel NF-YAx splice variant, with NF-Y functional modifying activity, which forms the subject of this report.
Materials and methods
Aim, design and setting
The aim of this study was to report the discovery and characterization of the novel alternative NF-YAx splice variant, discovered as the exclusive NF-YA isoform expressed in an advanced stage 3 NB, expressed during mouse embryo development and induced in NB cells by doxorubicin, and to provide insights into its potential function in NB pathogenesis and progression. Experimental design included cloning and sequence characterization, expression vector construction, protein expression and characterization, transient and stable transfection and biological characterizations in terms of growth, cytotoxicity and tumorigenicity.
Characteristics of participants and materials
RNAs from primary human stage 1, 2, 3 and 4 NBs (International NB Staging System) were kindly provided by the tissue bank of the Italian Association of Haematology and Pediatric Oncology (Genova, Italy). Patient age, sex or gender was not provided. RNAs from stage E8.5-E18.5 mouse embryos were provided by Dr. Rita Gallo (University of L’Aquila, L’Aquila, Italy). The following cell lines were used: stable transfected pcDNA3.1 control, TrkAIII and TrkT3 SH-SY5Y cell lines ; SV40 large T-antigen immortalized ST14A rat neuronal progenitor cells (kindly provided by Dr. Elena Cattaneo) ; NB cell lines SK-N-SH (HTB-11); SH-SY5Y (CRL-2266); IMR32 (CCL-127); SK-N-BE (CRL-2271); LAN-1 (HB-8568), HEK-293 embryonic neural lineage kidney cells [41, 42] (CRL-1573) and SL2 Sneider’s Drosphila Line 2 insect (CRL-1963) cell lines from ATCC (Rockville, MD); LAN-5 NB cell line from DSMZ (ACC673) (Braunschweig, Germany), and SMS-KCNR and SH-EP NB cell lines from Dr. U.P. Thorgiersson (NCI, NIH, Bethesda, MD). Cells were routinely grown at 37 °C and 5% CO2, in ATCC or DSMZ recommended culture medium (RPMI, DMEM), supplemented with appropriate antibiotics (Zeocin for stable-transfectants, penicillin and streptomycin) and 10% foetal calf serum. SL2 cells were cultured in Schneider’s Drosophila medium supplemented with 10% foetal bovine serum. Human foetal brain-derived neural stem cells were maintained as neurospheres in NPMM maintenance medium supplemented with Growth and Supplement Singlequotes, as directed (Cambrex Bio Science, Milan, Italy). All cell lines were authenticated at source. The following antibodies were used: KIF1Bβ (A301-055A-M, Bethyl Laboratories, Montgomery, TX); EglN3 (BS0532R-TR, Bioss, Woburn, MA); β-Actin (sc-47,778), Caspase 9 (sc-819), Mcl1 (sc56077), Bcl2 (sc-509), NF-YA (sc-10,779), NF-YB (sc-13,045), Sp1 (sc-420), Ubiquitinated proteins (sc-9133, P62 (sc48402), Bcl-xL (sc-136,207) (Santa Cruz, Dallas, TX); Alix (92880), Caspase 3 (9662), PARP (9642), pJNK (4668) and JNK (9258) (Cell Signaling Technology, Danvers, MA).
RNA purification and RT-PCR
Total RNAs from cell lines and xenograft tumor tissues were purified using an RNA-easy Plus RNA purification kit, as directed (Qiagen, Milan, Italy). All RNAs (1 μg) were reverse transcribed using a Maxima H Minus First Strand cDNA synthesis kit, as directed (Thermo-Fisher Scientific, Waltham, MA) and subjected to semi-quantitative PCR (35 cycles of 40 s at 95 °C, 30 s at 1° below melting temperature and 42 s at 72 °C), using the following primer sets: human NF-YA (5′-AATAGTTCGACAGAGCAGATTG-3′ and 5′-TCCTGCCAAAC TGGCTGCTGGGAT-3′, generating products of 575 bp for NF-YAl, 488 bp for NF-YAs and 338 bp for NF-YAx); mouse NF-YA (5′-AAACAGCAATAGTTCCACAGAGCAGATCG-3′ and 5′-CGATCTGCTCTGTGGAACTATTGCTG TT T-3′; products of 582 bp for NF-YAl, 494 bp for NF-YAs and 361 bp for NF-YAx); CD133 (5′-GCTGATGTTGAAACTGCTTGAG-3′ and 5′-GGTGCCGTTGCCTTGG-3′); SOX2 (5′-CCAAGACGCTCATGAAGAAG-3′ and 5′-TGGTCATGGAGTTGTACTGC-3′;); CD117 (5′-AACGCTCGAGTACCTGTGAA-3′ and 5′-GACAGAATTGATCCGCACAG-3′); Nestin (5′-GGATCAGATGACATTAAGACCC-3′ and 5′-TCCAGTGGTTCTTGAATTTCC-3′); Nanog (5′-TCTCTCCTCTTCCTTCCTCCATG-3′ and 5′-CTGTTTGTAGCTGAGGTTCAGGATG-3′); p75NGFR (5′-GTGGGACAGAGTCTGGGTGT-3′ and 5′-AAGGAGGGGAGGTGATAGGA-3′); Bmi1 (5′-GCTGCCAATGGCTCTAATGA-3′ and 5′-TGCTGGGCATCGTAAGTATCTT-3′); Bcl2 (5′-GACT TCGCCGAGATGTCCAGC-3′ and 5′-CAAGCT CCCACCAGGGCCAAAC-3′); Mcl1 (5′-AAGCCAATGGGCAGGTCT-3′ and 5′-TGTCCAGTTTCCGAAGCAT-3′); BclxL (5′-GAACGGCGGCTGGGATACTTTTG-3′ and 5′-GT GAATTCTGAGGCCAAGGGAAC-3′); PUMA (5′-GCATGGGGTCTGCCCAGG-3′ and 5′-CCGCCGCTCGTACTGTGC-3′); BAD (5′-CAGTCACCAGCAGGAGCAG-3′ and 5′-CCCTCCCTCCAAAGGAGACAG-3′); BAX (5′-ATGGACGGGTCCGGGGAGCAG-3′ and 5′-GCACCAGTTTGCTGGCAAAGTAG-3′); EglN3 (5′-GCGTCTCCAAGCGACA-3′ and 5′-GTCTTCAGTGAGGGCAGA-3′); KIF1Bβ (5′-CAGTGACTGTAAGTTGTCTGATATA-3′ and 5′-GTAAAGAGGCTCCTTGAAAT-3′); GAPDH (5′-CTGCACCACCAACTGCTTAG-3′ and 5′-AGGTCCACCACTGACACGTT-3′) and 18S rRNA (5′-AAACGGCTACCA CATCCAAG-3′ and 5′-CGCTCCCAAGATCCAACTAC-3′).
NF-YA variant cDNA cloning and sub-cloning
NF-YAx coding cDNA was cloned from the primary human stage 3 NB, exhibiting exclusive NF-YAx expression. NF-YAl, NF-YAs, NF-YB and NF-YC cDNAs were cloned from SH-SY5Y mRNA, using primers 5′-CCTCCTGATTGGGTTTCGGAGT-3′ and 5′-GGGGTTAGGACACTCGGATGAT-3′ spanning fully-spliced human NF-YAl cDNA (NM_002505.5); primers 5′-GGTTCTGACAGTATTTCATGACAA-3′ and 5′-CAGATCATGAAAACTGAATTTGCT-3′ spanning fully spliced human NF-YB cDNA (NM_006166) and primers 5′-GGACTCCTGAGCACAGTTGTCGAG-3′ and 5′-CCTTGGCCTTGCCAGCTCAGGCCC-3′ spanning fully spliced NF-YC cDNA (HSU78774). RT-PCR products were agarose gel purified, sub-cloned into the TA cloning vectors (Invitrogen, Thermo-Fisher Scientific, Waltham, MA) and sequenced in an ABI prism automated DNA-sequencer. Full length NF-YAl, NF-YAs, NF-YAx, NF-YB and NF-YC cDNA clones were then sub-cloned into mammalian pcDNA3.1Zeo (+) and pAc insect expression vectors (Thermo-Fisher Scientific, Waltham, MA), together with cDNA for full-length Sp1, available in our laboratory. NF-YAdn expression vector was from Dr. R. Mantovani (University of Milan, Italy).
SiRNA knockdown was achieved using a TriFECTa Dicer-Substrate RNAi kit and three KIF1B-specific Dicer-Substrate siRNA duplexes (hs.Ri.KIF1B.13.1: 5′-UCCACUGAGAAGGUCAGUAAAAUCA-3′ and 5′-UGAUUUUACUGACCUUCUCAGUGGAAA-3′; hs.Ri.KIF1B.13.2: 5′-CAAGGAAUCCAAAUGCAUCAUUCAG-3′ and 5′-CUGAAUGAUGCAUUUGGAUUCCUUGCU-3′; hs.Ri.KIF1B.13.3: 5′-AGCUAUUGAACGUUUAAAG GAAUCA-3′), as directed (Integrated DNA Technologies, Bologna, Italy). Briefly, 1 × 105 cells/ml HEK-293 cells on 24 well plates were grown overnight to ≈ 80% confluence and co-transfected with either 50 nM of negative control siRNA duplex (5′-CGUUAAUCGCGUAUAAUACGCGUAT-3′ and 5′-AUACGCGUAUUAUACGCGAUUAACGA C-3′) or 50 nM of a mix of KIF1B-specific siRNA duplexes, plus 1 μg/ml of pcDNA3.1 NF-YAx expression plasmid or empty pcDNA control plasmid, using TransIT-X2 Dynamic Delivery System, as directed (Mirus Bio, www.mirusbio.com/6000). Sham-transfected controls received transfection reagent alone. Knockdown of NF-YAx-induced KIF1Bβ expression was confirmed at 24 h by RT-PCR and the effect of knockdown assessed on NF-YAx-induced cell death at 48 h, by phase contrast and AO/EtBr cell death assays. Knockdown experiments were performed in duplicate and repeated (n = 4). Transfection efficiency was confirmed using HPRT-S1 DS positive control and validated using a negative control duplex (NC1), as directed (Integrated DNA Technologies, Bologna, Italy).
Immunoprecipitations and Western blotting
Cell proteins, extracted in lysis buffer (PBS containing 0.5% sodium deoxycholate, 1% NP40, 0.1% SDS, 1 mM sodium orthovanadate, 1 mM PMSF, 1 μg/ml of pepstatin A and Aprotinin), were analysed by reducing SDS-PAGE/Western blotting, as previously described . For immunoprecipitations, cell extracts (200–500 μg), pre-cleared with IgG and Protein A Sepharose (Fast flow, Sigma-Aldrich, St Louis, MI), were incubated overnight with primary antibody (0.1–1.0 μg), with rotation at 4 °C, then 20 μl of Protein A Sepharose in lysis buffer was added and incubated for 30 min at 4 °C. Protein A Sepharose/IgG conjugates, collected by centrifugation (10,000 x g for 5 min), were washed in lysis buffer and examined by reducing SDS-PAGE/Western blotting.
Stable SH-SY5Y transfection
SH-SY5Y cells were stable-transfected with pcDNA 3.1 NF-YAl, NF-YAs or NF-YAx or empty pcDNA3.1Zeo expression vectors, as previously described . Zeocin-resistant colonies were isolated upon appearance, clonally expanded and characterized by RT-PCR and Western blot.
Cells at 1 × 105/ml were grown overnight to ≈80% confluence on 6 or 24 well plates, transfected with NF-YA variant or empty plasmid DNA (1 μg/ml) in Fugene HD, as directed (Promega, Madison, WI), washed at 6 h, grown in complete medium, monitored at 6, 12, 24 and 48 h for evidence of cytotoxicity and photographed. At 48 h, suspension and adherent cell populations were separated and either analysed in AO/EtBr cell-death assay, used for RNA purification and RT-PCR analysis or extracted for Western blotting. Transfection efficiency was estimated using pEGFP-N1 green fluorescent protein reporter plasmid (Clontech, Mountain View, CA).
Southern blots of NF-YA and GAPDH RT-PCR products were prepared, as previously described . Briefly, RT-PCR products from primary human NB RNAs (1 μg) and doxorubicin-treated SH-SY5Y mRNAs (100 ng), resolved by 1.5% agarose gel electrophoresis, were in-gel denatured with 0.5 M NaOH, transferred to Hybond N+ nylon membranes (Amersham International, Little Chalfont, UK) and hybridized with P32-labelled probes for either NF-YA exon C (5′-CAAGGGCAGCCATTAATGGTGCAGGTCAGTGGAGGCCAGCTAATCACATCAACTGGCCAACCCATCATG-3′), NF-YA exon D (5′-AGGGCCAGCAGGGCCAGACCCAGCAGATCATCATCCAGCAGCCCCAGACGGCT GTCACTGCTGGCCAGA-3′) or GAPDH and visualized by autoradiography.
MTS proliferation and thymidine incorporation assays
In MTS proliferation assays (Promega, Madison, WI), 5 × 103 cells in 96-well microtiter plates were grown for 0–4 days. MTS reagent (20 μl) was added at 24-h intervals to individual wells and colorimetric conversion monitored at 492 nm every 30 min for 4 h in a microtiter plate reader. Assays were performed in duplicate and repeated three times (n = 6). For thymidine incorporation assays, 5 × 103 cells seeded on 96-well microtiter plates were grown for 0–4 days. At 24-h intervals, 1 μCi of tritiated [3H] thymidine (Amersham International, Little Chalfont, UK) was added and incubated for 4 h. Cultures were then washed in pre-warmed PBS, trypsinized, vacuum transferred onto 3MM Whatman paper filters and counted in a β-Scintillation counter (Becton, Dickinson and Company, Franklyn Lakes, NJ). Assays were performed in duplicate and repeated 3 times (n = 6).
In vitro protein translation
NF-YAl, NF-YAs, NF-YAx, NF-YB and NF-YC proteins were translated in vitro from linearized cDNAs (1 μg), in the presence of 35S labelled methionine, for 90 min at 37 °C, in a rabbit reticulocyte TNT T7 Quick Coupled Transcription/translation System, as directed (Promega, Madison, WI). Translated proteins were resolved by reducing SDS-PAGE autoradiography. For EMSAs, 35S-labelled methionine was omitted from in vitro translation reactions.
In vitro Colony formation assays
Cell suspensions (1 × 105/ml) in 6 well culture plates were grown overnight to ≈80% confluence, transfected with empty pcDNA3.1 or pcDNA3.1 NF-YAl, NF-YAs or NF-YAx expression vectors (2 μg/ml) in Fugene, as directed (Promega, Madison, WI), and grown for a further 48 h. Zeocin (200 μg/ml) was then added, replaced every 3 days and colonies grown for 21 days, fixed in 100% Methanol, stained in 0.5% w/v crystal violet, 1% Formaldehyde, 1x PBS, 1% methanol, washed in tap water, air dried, photographed and counted. Assays were performed in triplicate and repeated 3 times (n = 9).
Neural stem cell sphere growth assays
Cells seeded at 1 × 105/ml in vertical T75 culture flasks (NUNC) were grown in DMEM/F1, containing 0.6% glucose, 1 x B27 supplement (Gibco, Thermo-Fisher Scientific, Waltham, MA), 1 x N2 supplement (Gibco), antibiotics (1 x Pen/Strep), glutamine, EGF (40 ng/ml) and FGF (40 ng/ml), and monitored for formation of neuro-spheres consisting of > 50 cells. Assays were performed in duplicate and repeated (n = 4).
Direct RT-PCR sequencing
Agarose gel-purified RT-PCR product sequencing was performed using BigDye Direct Cycle Sequencing kit, as directed (Thermo-Fisher Scientific, Waltham, MA). Briefly, cDNAs (final concentration 4 ng/ml) were added to reaction mixtures containing forward or reverse primer (0.8 μM), BigDye direct PCR master mix, deionized water in a final volume of 20 μl and subjected to 35 PCR cycles (94 °C and for 10 min and 96 °C for 10 s, 55 °C for 5 min and 60 °C for 4 min). PCR reactions were mixed with 2 ul of NaAc and 50 μl absolute ethanol, incubated for 10 mins, spin-dried at 120000 rpm for 20 mins and sequenced in a mono-capillary sequencer, as directed (ABI PRISM 310, Life Technologies, Monza, IT).
Tumor growth in soft agar
Single-cell suspensions (passed through a gauge × 18 syringe needle) of 5 × 104 cells were mixed in a 33% solution of agar (BiTec; Difco) in RPMI containing 5% FCS at 37 °C and layered onto a solid 0.6% agarose substrate prepared in the same growth medium. Following agar solidification, complete medium was added, replaced every 2 days and growth monitored over 14 days. Assays were performed in duplicate and repeated (n = 4).
Tumor growth in NGS mice
Tumorigenesis in vivo was performed as previously described . Stable-transfected SH-SY5Y cells, prepared as a single-cell suspension in PBS without Ca2+ and Mg2+, were injected subcutaneously into the flanks of anesthetized, 6-week-old female NGS mice (Charles River, Calco, Italy), at a concentration of 1 × 107 cells in 200 μl per site. Tumor initiation was recorded at a minimum volume (tumor length × [tumor width]2 × 0.44) of 12 mm3 and animals sacrificed at 21 days. Groups consisted of 5 animals and the assay performed once (n = 5). All experiments were performed in accordance with Italian national and Rome University guidelines.
Nuclear extracts and EMSA
Nuclear extracts and EMSAs were performed as previously described . Briefly, binding reactions, performed at room temperature for 20 min, contained 32P-5′ end-labelled, double-stranded oligonucleotide probes, 2 μg of poly-deoxyinosinic-deoxy-cytidilic acid, 5 μg of nuclear extract and additional competitor DNAs or antibodies as specified in the figure legends. The oligonucleotides used were as follows: inverted CCAAT-box: 5′-GGGAGACCGTACGTGATTGGTTAATCTCTT-3′ and non-specific: 5′-GGTCGATAGGGAATTTACACGC-3′ oligonucleotides. All oligonucleotides were double stranded, the complementary affinity strands are not indicated. Assays were repeated 3 times (n = 3).
Cell death assays
Cell death assays were as previously described . Suspension or adherent cells, detached in ice cold PBS containing 1 mM EDTA, were transferred into sterile tubes, pelleted at 1000 x g at 4 °C, washed in ice cold PBS, re-pelleted, re-suspended in 25 μl of PBS containing 2 μl of acridine orange/ethidium bromide solution (100 μg/ml acridine orange and 100 μg/ml ethidium bromide in PBS), mounted onto glass slides, examined immediately under a Zeiss “Axioplan-2” fluorescence microscope, digitally photographed and dead (orange/red nuclei) and live cells (green nuclei) counted. Assays were performed in duplicate and repeated 3 times (n = 6).
Cells grown on Nunc glass chamber slides (Sigma-Aldrich, St Louis, MI) were washed in PBS, fixed and permeabilized in 100% ice cold methanol (− 20 °C), incubated for 1 h in blocking solution (1% bovine serum albumin in PBS-0.03% TX100), incubated for 2 h with primary antibody diluted in blocking solution at room temperature, washed, incubated with secondary fluorochrome-conjugated antibody diluted in blocking solution, for 1 h at room temperature, mounted with VectaMount (Vector Laboratories, Berlingame, CA) and observed using a Zeiss Axioplan 2 fluorescence microscope with digital camera and Leica M500 Image Manager software.
Data were analysed by Student’s t-test (https://www.graphpad.com/quickcalcs/ttest1.cfm) and statistical significance associated with probabilities of ≤0.05.
This study, originally designed to evaluate fully-spliced NF-YAl and alternatively exon B spliced NF-YAs mRNA expression in primary human NBs, resulted in the unexpected discovery of a novel NF-YAx splice variant.
NF-YAx cloning and sequencing
NF-YAl (1284 bp) and NF-YAs (1197 bp) coding cDNAs from SH-SY5Y cells and full-length NF-YAx cDNA (947 bp) from a stage 3 NB were sequenced, and NF-YAx characterized as a novel splice variant exhibiting in-frame exons B (87 bp, amino acids (aa) 26–54), D (132 bp, aa 103–146) and partial F (18 bp, aa 182–188) skipping (Fig. 1a) .
NF-YAx is expressed in NBs and mouse embryos
RT-PCR of NB RNAs using primers spanning NF-YA exons A-F, generated 575 bp (NF-YAl), 488 bp (NF-YAs) and 338 bp (NF-YAx) products that hybridized to NF-YA exon C probe in Southern blots (Fig. 1b). The 338 bp product, characterized as NF-YAx by direct PCR sequencing (not shown), was the exclusive isoform expressed in a stage 3 NB and a minor isoform, together with NF-YAl and NF-YAs, in two stage 2 NBs but was not detected in either stage 1 or 4 NBs. RNAs from SK-N-MC, SK-N-SH, SH-SY5Y, KCNR, IMR-32 and SK-N-BE NB cell lines, HEK293 embryonic neuronal-lineage kidney cells and human neonatal brain stem cells (SCs) exhibited predominant NF-YAs but not constitutive NF-YAx expression; CHP-126, LAN-1 and LAN-5 NB cell lines exhibited equivalent NF-YAs and NF-YAl but not NF-YAx expression and SHEP NB cells exhibited predominant NF-YAl but not NF-YAx expression. Predominant NF-YAs with minor NF-YAx expression characterized ST14A embryonic striatal neuronal progenitors under progenitor-maintaining (33 °C) and differentiation-inducing (39 °C) conditions  (Fig. 1c and d). Western blots confirmed predominant NF-YAs protein expression in SH-SY5Y cells, predominant NF-YAl expression in SHEP cell extracts (Fig. 1e). NF-YA isoforms translated in vitro were characterized by SDS-PAGE with approximate molecular masses of 46 kDa (NF-YAl), 42 kDa (NF-YAs) and 35 kDa (NF-YAx) (Fig. 1f). Samples were not available for NF-YAx protein analysis in primary NBs.
In mouse embryo RNAs, NF-YAs mRNA expression was predominant in stage E8.5 embryos, NF-YAl mRNA expression predominant in stage E9.5-E18.5 embryos and NF-YAx mRNA expression, confirmed by direct sequencing (data not shown), was detected in main body (M) but not head (U) or lower limbs (L) RNAs from stage E12.5, 13.5 and 14.5 embryos (Fig. 2a).
Doxorubicin induces NF-YAx expression in SH-SY5Y cells
Doxorubicin (10 μM) induced NF-YAx mRNA expression in SH-SY5Y cells at times of maximal > 95% cell death (Fig. 2b), which in Southern blots exhibited time and concentration-dependent induction kinetics at concentrations of 5 and 10 μM, with NF-YAx identified by hybridization to the NF-YA exon C (Exon C) but not exon D (Exon D) probe (Fig. 2c). NF-YAx protein expression was also detected in SH-SY5Y nuclear extracts at 24 and 48 h, following Doxorubicin (10 μM) treatment (Fig. 2d). Necrostatin-1 (100 μM) significantly reduced cell death induced by 10 μM doxorubicin from 94 ± 7.5% to 74 ± 8.7% (p < 0.002, n = 12), at 48 h, (Fig. 2e) and doxorubicin-induced cell death was also characterised by vacuolation, swelling and lysis, consistent with a proportion of necroptotic cell-death (Fig. 2e). NF-YAx expression in SH-SY5Y cells was not induced by treatment with either DTT (5 mM), thapsigargin (10 ng/ml), tunicamycin (1 μM), TPA (10 ng/ml), dibutyryl cAMP (0.5 mM), retinoic acid (10 μM) or CoCl2 (150 μM) over 48 h and was not detected in either stable TrkT3 or TrkAIII SH-SY5Y-transfectants  (Fig. 2f).
NF-YAx forms CCAAT-box binding NF-Y complexes
In EMSAs, in vitro-translated NF-YAl, NF-YAs and NF-YAx all formed specific CCAAT-box binding complexes with NF-YB and NF-YC, confirmed by competition EMSAs in which excess un-labelled CCAAT-box oligonucleotide (S-oligo), anti-NF-YA and anti-NF-YB antibodies but not non-specific oligonucleotide (NS-oligo) or pre-immune IgG abrogated binding (Figs. 1f and 3a and b). EMSAs performed with constant NF-YAl, NF-YB and NF-YC and increasing NF-YAx levels (0–4 μl), or constant NF-YAx, NF-YB and NF-YC and increasing NF-YAl levels (0–4 μl), revealed reciprocal NF-YAx/NF-YAl substitution, confirmed by changes in electrophoretic mobility (Fig. 3b). EMSAs also detected specific CCAAT-box binding NF-Y complexes, containing NF-YAl, NF-YAs and NF-YAx, in nuclear extracts from stable NF-YA isoform SH-SY5Y-transfectants, exhibiting similar levels of nuclear-localized NF-YAl, NF-YAs and NF-YAx expression (Fig. 3c, d and e). NF-YAx-specific antibodies, however, are not available to confirm NF-YAx recruitment to gene promoters in vivo by chromatin immunoprecipitation assay.
NF-YAx does not bind Sp1
NF-YAl and NF-YAs bind Sp1 and regulate Sp1 recruitment to promoters [11, 44]. IP/Western blots confirmed Sp1 pull-down of NF-YAl and NF-YAs but not NF-YAx in Sp1-negative SL2 cells , co-transfected with pAc Sp1 and NF-YAl, NF-YAs and NF-YAx insect expression vectors (Fig. 3f). In stable SH-SY5Y-transfectants evidence for Sp1 pull down of NF-YAl and NF-YAs but not NF-YAx was also supported by augmented immunoreactivity overlapping non-specific (NS1 and NS2) bands, consistent with Sp1 binding of NF-YAl and NF-YAs but not NF-YAx (Fig. 3g).
NF-Yx is cytotoxic to embryonal neural-lineage progenitors and NB cells and inhibits colony formation
Transient 48-h transfection of ST14A, HEK-293 and SH-SY5Y with pEGFP-N1 reporter gene, resulted in transfection-efficiencies of 28.3 ± 8.2%, 55.9 ± 4.7% and 30.3 ± 6.5% (mean ± SD) (n = 6, for all three), respectively (not shown). Transient 48-h ST14A, HEK-293 and SH-SY5Y transfection with empty pcDNA, NF-YAl, NF-YAs, NF-YAx or NF-YAdn vectors, resulted in similar levels of NF-YA isoform expression (Fig. 4a, shown for HEK-293 cells only).
In ST14A cells, NF-YAx and NF-YAdn transfection significantly increased the percentage of suspension cells at 48 h, from 7.5 ± 2.25% (n = 6) in pcDNA, 8.5 ± 2.1% (n = 6) in NF-YAl and 8 ± 5.2% (n = 6) in NF-YAs-transfectants, to 25 ± 3.1% (n = 6) in NF-YAx and 27 ± 4.2% (n = 6) (means ± SD) in NF-YAdn-transfectants (p < 0.0001 for both NF-YAx and NF-YAdn vs controls), with 98 ± 5.5% (n = 6) of NF-YAx-induced suspension cells and 95 ± 7.6% (n = 6) (means ± SD) of NF-YAdn-induced suspension cells confirmed dead by AO/EtBr assay (Fig. 4b and c). In SH-SY5Y cells, NF-YAx and NF-YAdn also significantly increased the percentage of suspension cells at 48 h, from 4.25 ± 2.25% (n = 6) in pcDNA, 6.25 ± 2.1% (n = 6) in NF-YAl and 5 ± 5.2% (n = 6) in NF-YAs-transfectants, to 26 ± 3.1% (n = 6) in NF-YAx and 26.75 ± 4.2% (n = 6) in NF-YAdn-transfectants (means ± SD) (p < 0.0001 for NF-YAx and NF-YAdn vs controls), with 88 ± 5.5% (n = 6) of NF-YAx-induced suspension cells and 93 ± 5.6% (n = 6) (means ± SD) of NF-YAdn-induced suspension cells confirmed dead by AO/EtBr assay (Fig. 4b and c). In HEK-293 cells, NF-YAx and NF-YAdn also significantly increased the percent of suspension cells, from 5.25 ± 2.22% (n = 6) in pcDNA, 8.7 5 ± 1.9% (n = 6) in NF-YAl and 11.25 ± 4.8% (n = 6) in NF-YAs-transfectants, to 40.75 ± 3.1% (n = 6) in NF-YAx and 39.6 ± 4.3% (n = 6) in NF-YAdn-transfectants (means ± SD) (p < 0.0001 for NF-YAx and NF-YAdn vs controls), with 92 ± 6.5% (n = 6) of NF-YAx-induced suspension cells and 94 ± 7.8% (n = 6) (means ± SD) of NF-YAdn-induced suspension cells confirmed dead by AO/EtBr assay (Fig. 4b and c). In contrast, NF-YAl, NF-YAs or NF-YAx did not induce the death of SH-SY5Y cells expressing TrkAIII oncoprotein (Fig. 4d). NF-YAx-induced death was not detected prior to 16 h, was maximal at 48 h, was not prevented by z-VAD-fmk (20 μM) but was significantly inhibited by necrostatin-1 (100 μM) , from 29 ± 3.1% (n = 6) to 12 ± 1.6% (n = 6) in SH-SY5Y cells and from 30 ± 4.5% (n = 6) to 15 ± 1.5% (n = 6) in HEK-293 cells (means ± SD) (P < 0.001, for both) (Fig. 5a).
In colony formation assays, sham-transfected SH-SY5Y, HEK-293 and ST14A cells did not form colonies and pcDNA-transfected counterparts formed maximum colony numbers. Compared to respective pcDNA controls, in SH-SY5Y cells, NF-YAl-transfectants formed 22.6 ± 7.5% fewer colonies (p < 0.001, n = 9), NF-YAs-transfectants 39.5 ± 2.1% fewer colonies (p < 0.0001, n = 9), NF-YAx-transfectants 91.1 ± 5.5% fewer colonies (means ± SD) (p < 0.0001, n = 18) and NF-YAdn-transfectants did not form any colonies (p < 0.0001, n = 18). In ST14A cells, NF-YAl transfectants formed 23.3 ± 7.5% fewer colonies (P < 0.0001, n = 9), NF-YAs 27.3 ± 2.1% fewer colonies (P < 0.001, n = 9) whereas NF-YAx and NF-YAdn-transfectants did not form colonies (P < 0.001, n = 9 for both). In HEK-293 cells, NF-YAl-transfectants formed 28 ± 9.1% fewer colonies and NF-YAs-transfectants 38.6 ± 8.5% fewer colonies (p < 0.0001 for both, n = 9 for each), whereas NF-YAx and NF-YAdn-transfectants did not form any colonies at all (Fig. 5b).
NF-YAx induces Necroptosis
NF-YAx-induced ST14A, SH-SY5Y and HEK293 cell-death was characterized by vacuolation, swelling and necrotic-lysis but not by apoptotic-body formation (Fig. 6a) or chromatin condensation (Fig. 4f). Western blots confirmed NF-YAx expression in dead suspension SH-SY5Y and HEK-293 cells (Fig. 6b) but did not detect caspase-3, caspase-9 or PARP-cleavage, enhanced JNK phosphorylation, reduced AIP/Alix expression or altered Bcl-xL, Bcl2, Mcl-1 and EglN3 expression (Fig. 6c and d). In both SH-SY5Y and HEK-293 cells, NF-YAx but not NF-YAl reduced Bmi1 expression, induced KIF1Bβ expression but did not alter EglN3 expression (Fig. 6c and d). In HEK-293 cells, siRNA knockdown of NF-YAx-induced KIF1Bβ expression significantly reduced NF-YAx-induced death from 24.6 ± 0.75% (n = 6) without siRNAs and 27.2 ± 1.2% (n = 6) with scrambled siRNAs to 5.4 ± 0.61% (n = 6) with KIF1B-specific siRNAs (means ± SD) (p < 0.0001 versus siRNA controls) (Fig. 6e). Western blots also detected a significant 187.08 ± 10.36% (n = 3) increase in ratio of ubiquitinated proteins to β-actin and significant 148.2 ± 8.6% (n = 3) increase in the ratio p62 to β-actin (means ± SD) (p < 0.036 for all comparisons) in extracts from HEK-293 cells transfected with NF-YAx compared to pcDNA, NF-YAl or NF-YAs transfected counterparts (Fig. 6f).
NF-YAx selects tumorigenic, doxorubicin-resistant cancer SCs
Duplicate stable SH-SY5Y NF-YAx-transfectants were established and compared to stable control, NF-YAl and NF-YAs-transfectants. Stable NF-YAl, NF-YAs, NF-YAx1 and x2-transfectants expressed similar levels of NF-YAl, NF-YAs and NF-YAx (Fig. 3c) and did not differ significantly in either proliferation or mitotic rates, assessed by MTS and thymidine-incorporation assays (Fig. 7a). NF-YAl, NF-YAs and NF-YAx stable transfectants all formed neuro-spheres in neural stem cell assays in vitro (Fig. 7b) and similar numbers of similar sized tumor spheroids in soft-agar tumorigenesis assays in vitro (Fig. 7c). In xenograft tumorigenesis assays in NSG mice, pcDNA, NF-YAs, NF-YAx1 and x2-transfectants formed sub-cutaneous tumours of similar dimensions, significantly larger than tumours formed by stable NF-YAl -transfectants (> 2 fold smaller, p < 0.049 compared to all other stable-transfectants, n = 5 per group) (Fig. 7d). In cell death assays, doxorubicin induced similar levels of stable pcDNA, NF-YAl and NF-YAs-transfectant cell death at all concentrations at 6, 24 and 48 h. In contrast, stable NF-YAx1 and x2-transfectants exhibited significantly enhanced survival at 6 but not 24 or 48 h (p < 0.0001, n = 6 for both) in the presence of 10 μM doxorubicin, at 6 and 24 but not 48 h (p < 0.0001, n = 6 for both) in the presence of 5 μM doxorubicin and at all time points in the presence of 0.01, 0.1 and 1 μM doxorubicin (p < 0.0001, n = 6 for both at all time points and doses) (Fig. 7e).
In RT-PCR assays NF-YAx1 and x2 transfectants exhibited higher levels of p75NTR, Nanog, Nestin and EglN3 expression than NF-YAl and NF-YAs-transfectants, similar levels of Sox-2, CD133 and CD117 expression to NF-YAs-transfectants, elevated above control and NF-YAl-transfectants, high-level Bmi1 expression but did not express KIF1Bβ. In contrast, NF-YAl-transfectants expressed lower levels of p75NTR and Nanog than the other transfectants. Bcl2, Mcl1, Bcl-xL, PUMA, BAD and Bax expression levels did not differ between stable transfectants (Fig. 8a, b and c). RNAs purified from xenograft tumours exhibited a similar pattern of CSC gene expression to corresponding cell cultures (Fig. 8d). PTC-209 (10 μM for 24 h) abrogated proliferation of all stable transfectants (Fig. 9a) but did not reduce Bmi1 expression nor induce KIF1Bβ mRNA expression in NF-YAx1 or x2-transfectants (Fig. 9b).
We report a novel development and genotoxic stress-regulated alternative splice mechanism for promoting embryonic neural-lineage progenitor and NB cell-death, characterized by a switch to alternative NF-YA splicing and expression of a novel cytotoxic NF-YAx extra short-form variant. This novel isoform, originally discovered in human primary stage 2 and stage 3 NB RNAs, was the exclusive NF-YA isoform expressed at a high-level in an advanced stage 3 NB and was characterised as a novel NF-YA splice variant, exhibiting in-frame exon B, D and partial F skipping, responsible for truncating NF-YA transactivation domain sequence. NF-YAx readily competed with fully-spliced NF-YAl in CCAAT-box binding NF-Y complex formation but in contrast to NF-YAl and NF-YAs isoforms did not bind Sp1 and, therefore, represents a functional modifier of one of the more important physiological and cancer-associated transcription factors. In mouse embryos, NF-YAx expression coincided with the reported phase of neurotrophin-regulated KIF1Bβ-dependent sympathetic neuroblast-culling, unrestrained NF-YAx expression induced KIF1Bβ-dependent necroptosis in neural-lineage progenitors and NB cells and association between doxorubicin-induced NF-YAx expression and necroptosis in NB cells, supports a pro-necroptotic cytotoxic function for NF-YAx and a potential role in KIF1Bβ-dependent suppression of NB initiation, during development. On the other hand, propagation through selection of tumorigenic, doxorubicin-resistant CSC-like stable NF-YAx SH-SY5Y transfectants, resistant to NF-YAx cytotoxicity, not only helps to explain the high-level exclusive NF-YAx expression detected in an advanced stage 3 NB but also supports an additional potential role for NF-YAx, within the tumour context, in disease progression and identifies a potential mechanism for doxorubicin-induced post-therapeutic relapse, through CSC selection.
NF-YAx cDNA was cloned from a stage 3 NB and sequence characterised as a novel NF-YA splice variant with in-frame exons B, D and partial F sequence skipping, adding to existing NF-YAl, NF-YAs and NF-YA L2–6 variants [8, 12, 13, 17], with potential implications for NB pathogenesis and progression. NF-YAx expression, detected in 20% (3/15) NF-YA-positive NB RNA samples, was the exclusive high-level NF-YA mRNA isoform expressed in a stage 3 NB but was not detected in stage 1 and 4 NBs, human neonatal neural stem cells, HEK-293 human neural-lineage embryonic kidney cells and 14 human NB cell lines. Although samples were not available for NF-YAx protein detection in primary NBs, expression of the endogenous NF-YAx protein was confirmed in SH-SY5Y cells, following treatment with doxorubicin, corroborating the characterisation of NF-YAx cDNA as a complete, in-frame, non-mutated NF-YA splice variant that is readily translated into the NF-YAx protein in vitro and in vivo. Considering the relatively small number of NB samples analysed in this study, however, the possibility that NF-YAx mRNA expression is restricted to localized NB disseminated at most to local lymph nodes (Stages 2 and 3)  and may differentiate stage 3 from other disease stages, must await confirmation in a future larger NB cohort study.
NF-YAx mRNA was also expressed during mouse stage E12.5-E14.5 embryo development, consistent with identical exon organization of human (NCBI: NM_002505.5) and mouse (NCBI: NM_001110832.1) NF-YA genes. NF-YAx expression during stage E12.5-E14.5 embryo development coincides with the reported phase of neurotrophin-regulated EglN3/KIF1Bβ-dependent sympathetic neuroblast-culling, required for sympathetic nervous system development and suppressing NB initiation [25,26,27,28, 48,49,50,51,52,53,54,55]. Furthermore, unrestrained NF-YAx expression in neural-related progenitors and NB cells reduced expression of the KIF1B repressor Bmi1 , induced KIF1Bβ expression, promoted KIF1Bβ-dependent necroptosis and abrogated colony forming capacity in vitro. This suggests that NF-YAx expression during mouse embryo development may be involved in KIF1Bβ-dependent sympathetic neuroblast culling and subsequent suppression of NB initiation. In contrast, unrestrained NF-YAx expression did not kill NB cells expressing TrkAIII oncoprotein . This bears similarity to NGF-activated TrkA protection of sympathetic neuroblast against KIF1Bβ-dependent death [26, 27] and suggests that TrkAIII, which associates with NB, may facilitate NB initiation by preventing KIF1Bβ-induced death in a manner analogous to 1p36-deletion of KIF1B, germline KIF1B mutations or Nmyc amplification [25,26,27,28,29,30,31,32]. NF-YAx induction of KIF1Bβ expression, however, did not associate with changes in EglN3 expression, suggesting that this novel death-mechanism may act down stream of EglN3 through NF-YAx reduction of Bmi1 expression. Whether EglN3 promotion of alternative splicing , extends to NF-YAx remains to be determined.
NF-YAx cytotoxicity to neural-related progenitors and NB cells was not associated with caspase 3, caspase 9 or PARP cleavage nor with changes in Bcl-2, BclXL, Mcl-1, NOXA, PUMA, BAX, BAD and AIP1/Alix expression or JNK-phosphorylation, was not inhibited by the pan-caspase inhibitor z-VAD-fmk but was characterized by cell vacuolation, swelling and necrotic-lysis and inhibited by the RIPK-1 inhibitor necrostatin-1 , consistent with necroptosis rather than apoptosis [57, 58]. NF-YAx also increased the levels of ubiquitinated and p62 proteins in HEK-293 cells, suggesting reduced autophagic flux, consistent with reports that both Bmi1 inhibition and KIF1Bβ expression promote lysosomal/autophagosome uncoupling and necroptosis [59, 60], which may explain the NB suppressing function of KIF1Bβ [25, 32, 33, 54, 61, 62]. Although KIF1Bβ-induced apoptosis [32, 33, 57, 62, 63] was not detected in this study, we do not exclude potential roles for KIF1Bβ-induced ROS-mediated mitophagy or Drp1-dependent mitochondrial-fission in this death mechanism [33, 60, 64,65,66].
NF-YAx was significantly more cytotoxic to neural-related progenitor and NB cells than either NF-YAl or NF-YAs, contrasting with a previous report that unrestrained NF-YAl expression is highly cytotoxic and induces p53-dependent apoptosis . This discrepancy can be explained by fact that p53 is compromised by SV40 large T-antigen in ST14A cells , by adenovirus-5 in HEK-293 cells [41, 42] and by different mechanisms in SH-SY5Y cells [67,68,69], suggesting that this novel NF-YAx-dependent necroptosis mechanism may be restricted to genotoxic-stress under p53 compromised conditions. Furthermore, the sensitivity of ST14A neural progenitors, which exhibit low-level constitutive NF-YAx expression, to NF-YAx-induced necroptosis, confirms that this mechanism depends upon predominant NF-YAx expression. Unrestrained NF-YAdn expression, which contains a DNA binding-domain mutation that prevents NF-Y binding and transcription , also induced HEK-293 and SH-SY5Y necroptosis. This indicates that survival can switch to necroptosis in these cell types when either NF-Y binding is prevented (i.e. with NF-YAdn) or when a change in NF-Y function reaches a critical level (i.e. with NF-YAx ). Similar necroptotic-like neuronal death has also been reported during development [25, 48, 54, 62, 71, 72]. A pro-necroptotic role for NF-YAx was also supported by the association between doxorubicin-induced NF-YAx expression and SH-SY5Y cell-death, which was also characterised by vacuolation, swelling and cell lysis and significantly inhibited by the necroptosis inhibitor necrostatin-1 , consistent with a percentage of doxorubicin-induced necroptosis. Confirmation of the role of NF-YAx in this necroptotic proportion of doxorubicin-induced death, however, must await development of reagents specific for NF-YAx depletion.
Doxorubicin induction of alternative NF-YAx splicing not only implicates the DNA damage-associated alternative splice mechanism  and a role for NF-YAx in the response to genotoxic stress but also identifies NF-YAx as a potential biomarker of response to genotoxic therapy, suggesting that NF-YAx expression in primary NBs could also reflect neo-adjuvant genotoxic chemotherapy . Furthermore, the propagation of tumorigenic, doxorubicin-resistant, CSC-like stable NF-YAx SH-SY5Y transfectants, resistant to NF-YAx cytotoxicity, unveils an additional potential role for NF-YAx in genotoxic drug-induced, post-therapeutic relapse through CSC selection and maintenance. Moreover, the fact that NF-YAx expression was not induced by agents that promote ER-stress, hypoxic-stress, differentiation or malignant behaviour, suggests that this alternative splice mechanism may be relatively restricted to conditions of genotoxic stress in NB cells.
NF-YAx substitution of NF-YAl in DNA-binding NF-Y complexes and the capacity of NF-YAl and NF-YAs but not NF-YAx to bind Sp1 (this study, [11, 44, 75]), characterises NF-YAx as a potential modifier of both NF-Y and NF-Y-dependent Sp1 function. NF-YAx loss of Sp1 binding results from truncation of the aa 55–139 binding site for Sp1 in NF-YAl and NF-YAs to aa 55–103 in NF-YAx (this study, [44, 76]). This truncation may also compromise NF-Y interaction with ZHX transcriptional repressors that bind NF-YA aa’s 31–140, with potential to also alter NF-Y/ZHX-regulated genes expression, including MDR-1 chemotherapeutic-cytotoxicity-regulator and polo kinase-1 mitosis-regulator [77,78,79,80,81]. Transactivation-domain truncation may also weaken NF-Y function and interaction with other factors and the reduced size of NF-YAx, particularly in complexes with smaller NF-YB and NF-YC variants , may de-regulate transcription from promoters with precisely-spaced NF-Y-dependent transcriptional domains, e.g. ER-stress response gene promoters . Although, NF-YAx-specific antibodies are not yet available for confirmation, by chromatin immunoprecipitation assay, that NF-YAx is recruited to gene promoters in vivo, the fact that NF-YAx contains an intact DNA binding domain, competes with other isoforms in NF-Y complex formation and binds double stranded inverted CCAAT-box oligonucleotides in NF-Y complexed-form in vitro with similar kinetics to other NF-YA isoforms, strongly supports this probability. Although, we are only beginning to understand how NF-YAx may influence transcription, transient NF-YAx expression reduced Bmi1 and induced KIF1Bβ expression in neural-lineage progenitors and NB cells, and stable NF-YAx expression was associated with enhanced p75NTR, SOX2, Nestin, Nanog, CD117, CD133 and EglN3 expression.
Despite the acute cytotoxicity of NF-YAx to SH-SY5Y cells, stable SH-SY5Y transfectants exhibiting predominant NF-YAx expression were eventually propagated. These transfectants did not exhibit altered Bcl2, Bcl-xL or Mcl1 expression nor constitutive KIF1Bβ expression, suggesting selection based upon KIF1Bβ repression but not Bcl2-regulated mitochondrial apoptosis. Stable NF-YAx transfectants also exhibited high constitutive Bmi1 expression, consistent with resistance to NF-YAx-induced KIF1Bβ-dependent cytotoxicity. This indicates that the SH-SY5Y cell line contains sub-populations that are sensitive and resistant to acute NF-YAx cytotoxicity, characterised by differences in NF-YAx regulation of Bmi1 and KIF1Bβ expression. The Bmi1 inhibitor PTC-209 , however, failed to induce KIF1Bβ expression in stable NF-YAx-transfectants, indicating that Bmi1 alone is not responsible for KIF1Bβ repression in the resistant subpopulations. PTC-209 did, however, abrogate the proliferation of all SH-SY5Y-transfectants, consistent with the role of Bmi1 in proliferation and reports that PTC-209 induces G1/S checkpoint cell-cycle arrest [84,85,86,87,88,89]. We are currently investigating whether KIF1B repression in stable NF-YAx SH-SY5Y transfectants, as a potential selection mechanism, may involve inhibition EglN3 prolyl hydroxylase activity, required for KIF1Bβ expression and regulated by oxygen and high/low-order complexing [55, 90,91,92] and/or changes in NB-associated NF-Y and Sp1-regulated microRNA expression [93,94,95].
Stable NF-YAx SH-SY5Y transfectants were more CSC-like than stable control, NF-YAl or NF-YAs transfectants and were characterized by higher levels of p75NTR, Nestin, Nanog and EglN3 expression compared to all other stable transfectants, similar levels of SOX2, CD133 and CD117 expression to NF-YAs-transfectants elevated over other transfectants, high-level Bmi1 expression and KIF1Bβ repression, consistent with reports that elevated Bmi1 and EglN3 expression associated with embryonic staminality [84,85,86,87,88,89, 92, 95, 96], elevated p75NTR expression characterises neural crest and SH-SY5Y SCs [97,98,99] and KIF1Bβ-repression characterises undifferentiated aggressive NBs with 1p36 deletion or Nmyc amplification or that NBs that initiate in KIF1Bβ −/− mice [25, 26, 29,30,31,32,33, 100]. Stable NF-YAx SH-SY5Y transfectants were tumorigenic both in vitro and in vivo, confirming that acute NF-YAx cytotoxicity selects a tumorigenic NB CSC-like subpopulations resistant to NF-YAx-induced KIF1Bβ-dependent necroptosis. This provides an explanation for high-level exclusive NF-YAx expression in the advanced stage 3 NB and supports a role for NF-YAx, within the tumour context, in disease progression and suggests that doxorubicin induction of NF-YAx expression may also promote CSC selection in a potential mechanisms for post-therapeutic relapse, enhanced drug-resistance, increase adaptive plasticity and metastatic behaviour (this study, [101, 102]).
The discovery of NF-YAx mRNA in primary stage 2 and 3 NBs, its expression as the exclusive isoform in an advanced stage 3 NB, expression in stage E12.5 to E14.5 mouse embryos and induction by doxorubicin in NB cells, unveils a novel NF-YA splice mechanism and variant that is both regulated by and involved in NB, development and conditions of genotoxic-stress. NF-YAx substitution of other isoforms in NF-Y complexes and loss of Sp1 binding capacity characterises this novel isoform as a functional modifier of one of the more important physiological and cancer-associated transcription factors. NF-YAx induction of KIF1Bβ-dependent embryonic neural-related progenitor and NB necroptosis, association with doxorubicin-induced necroptosis and expression in murine embryos at times corresponding to the phase of neurotrophin-regulated KIF1Bβ-dependent sympathetic neuroblast culling, supports a predominant pro-necroptotic cytotoxic function for NF-YAx in neural progenitors and NB cells, with potential to suppress NB initiation during development. On the other hand, propagation by selection of tumorigenic, doxorubicin-resistant CSC-like stable NF-YAx expressing SH-SY5Y transfectants, resistant to NF-YAx cytotoxicity, not only helps to explain the exclusive high-level NF-YAx expression in an advanced stage 3 primary NB but also suggests that NF-YAx expression within the tumour context may promote disease progression and provides a possible doxorubicin-inducible mechanism for post-therapeutic relapse, through CSC selection and maintenance, with potential to enhance survival within stressful chemotherapeutic tumour microenvironments [62, 103].
Availability of data and materials
The data sets used and/or analysed during this study are either included in this published article or are available from the corresponding author on reasonable request.
B cell-specific Moloney murine leukaemia virus integration site 1
cancer stem cell
kinesin family member 1B
- NF-YAl :
long-form nuclear transcription factor YA
- NF-YAs :
short-form nuclear transcription factor YA
- NF-YAx :
extra short-form nuclear transcription factor YA
- NSG mice:
NOD severe combined immunodeficient gamma mice
specificity protein 1
Modrek B, Lee C. A genomic view of alternative splicing. Nucleic Acids Res. 2001;30:13–9.
Kalnina Z, Zayakin P, Silina K, Line A. Alterations in pre-mRNA splicing in cancer. Genes Chromosomes Cancer. 2005;42:342–57.
Tacconelli A, Farina AR, Cappabianca L, DeSantis G, Tessitore A, Vetuschi A, et al. TrkA alternative splicing: a regulated tumor-promoting switch in human neuroblastoma. Cancer Cell. 2004;6:347–60.
Cappabianca L, Guadagni S, Maccarone R, Sebastiano M, Chiominto A, Farina AR, Mackay AR. A pilot study of alternative TrkAIII splicing in Merkel cell carcinoma: a potential oncogenic mechanism and novel therapeutic target. J Exp Clin Cancer Res. 2019;38:424. https://doi.org/10.1186/s13046-19-1425-3.
Wang B-D, Lee NH. Aberrant splicing in Cancer and drug resistance. Cancers (Basel). 2018. https://doi.org/10.3390/cancers10110458.
Mantovani R. Nuclear factor Y in development and disease. Biochim Biophys Acta Gene Regul Mech. 2017;1860:523–4.
Maity SN. NF-Y (CBF) regulation in specific cell types and mouse models. Biochim Biophys Acta Gene Regul Mech. 2017;1860:598–603.
Li G, Zhao H, Wang L, Wang Y, Guo X, Xu B. The animal nuclear factor Y: an enigmatic heterotrimeric transcription factor. Am J Cancer Res. 2018;8:1106–25.
Gurtner A, Manni I, Piaggio G. NF-Y in cancer: impact on cell transformation of a gene essential for proliferation. Biochim Biophys Acta Gene Regul Mech. 1860;2017:604–16.
Zambelli F, Pavesi G. Genome wide features, distribution and correlations of NF-Y binding sites. Biochim Biophys Acta Gene Regul Mech. 2017;1860:581–9.
Suske G. NF-Y and SP transcription factors-New insights in a long-standing liaison. Biochim Biophys Acta Gene Regul Mech. 2017;1860:590–7.
Maity SN, de Crombrugghe B. Biochemical analysis of the B subunit of the heterotrimeric CCAAT-binding factor. J Biol Chem. 1992;267:8286–92.
Coustry F, Maity SN, de Crombrugghe B. Studies on transcription activation by the multimeric CCAAT-binding factor CBF. J Biol Chem. 1995;270:468–75.
Romier C, Cocchiarella F, Mantovani R, Moras D. The NF-YB/NF-YC structure gives insight into DNA binding and transcription regulation by CAAT factor NF-Y. J Biol Chem. 2003;278:1336–45.
Li X-Y, Hooft van Huijsduijnen R, Mantovani R, Benoist C, Mathis D. Intron-exon organisation of the NF-Y genes. J Biol Chem. 1992;267:8984–90.
Dolfini D, Gatta R, Mantovani R. NF-Y and the transcriptional activation of CCAAT promoters. Crit Rev Biochem Mol Biol. 2012;47:29–49.
Ge Y, Jensen TL, Matherly LH, Taub JW. Synergistic regulation of cystathionine-B-synthase-1b promoter by transcription factors NF-YA isoforms and Sp1. Biochim Biophys Acta. 2002;1579:73–80.
Zhu J, Zhang Y, Joe GJ, Pompetti R, Emerson SG. NF-YA activates multiple hematopoietic stem cell (HSC) regulatory genes and promotes HSC self-renewal. Proc Natl Acad Sci U S A. 2005;102:11728–33.
Grskovic M, Chaivorapol C, Gaspar-Maia A, Li H, Ramalho-Santos M. Systemic identification of cis-regulatory sequences active in mouse and human embryonic stem cells. PLOS ONE. 2007; e145, doi:https://doi.org/10.1371/jounal.pgen.0030145.
Dolfini D, Minuzzo M, Pavesi G, Mantovani R. The short isoform of NF-YA belongs to the embryonic stem cell transcription factor circuitry. Stem Cells. 2012;30:2450–9.
Basile V, Baruffaldi F, Dolfini D, Belluti S, Benatti P, Ricci L, et al. NF-YA splice variants have different roles on muscle differentiation. Biochim Biophys Acta. 1859;2016:627–38.
Bevilacqua MA, Faniello MM, Iovine B, Russo T, Cimino F, Costanzo F. Transcription factor NF-Y regulates differentiation of CaCo-2 cells. Arch Biochem Biophys. 2002;407:39–44.
Gu Z, Kuntz-Simon G, Rommelaere J, Cornelius J. Oncogenic transformation-dependent expression of a transcription factor NF-Y subunit. Mol Carcinogen. 1999;24:294–9.
Ruggeri P, Farina AR, Cappabianca L, Di Ianni N, Ragone M, Merolla S, et al. Neurotrophin and neurotrophin involvement in human neuroblastoma. In “Neuroblastoma” (ed: Hiroyuki Shimada, Intech, Rijeka, Cro) 2013: pp47–108.
Schlisio S, Kenchappa RS, Vredeveld LCW, George RE, Stewart R, Greulich H, et al. The kinesin KIF1Bβ acts downstream of EglN3 to induce apoptosis and is a potential 1p36 tumor suppressor. Genes Dev. 2008;22:884–93.
Kristiansen M, Ham J. Programmed cell death during neuronal development: the sympathetic neuron model. Cell Death Diff. 2014;21:1025–35.
Fell SM, Li S, Wallis K, Kock A, Surova O, Rakili V, et al. Neuroblast differentiation during development and in neuroblastoma requires KIF1Bβ-mediated transport of TrkA. Genes Dev. 2017;31:1036–53.
Jiang M, Stanke J, Lahti JM. The connections between neural crest development and neuroblastoma. Curr Top Dev Biol. 2011;94:77–127.
Ochiai H, Takenobu H, Nakagawara A, Yamaguchi Y, Kimura M, Ohira M, et al. Bmi1 is a MYCN target gene that regulates tumorigenesis through repression of KIF1Bβ and TSLC1 in neuroblastoma. Oncogene. 2010;29:2681–90.
Ohira M, Kageyama H, Mihara M, Furuta S, Machida T, Shishikura T, et al. Identification and characterization of a 500-kb homozygous deleted region in a neuroblastoma cell line. Oncogene. 2000;19:4302–7.
Bauters C, Pigny P, Dahia PLM. A germline mutation of the KIF1Bβ gene on 1p36 in a family with neural and non-neural tumors. Hum Genet. 2008;124:279–85.
Munirajan A, Ando K, Mukai A, Takahashi M, Suenaga Y, Ohira M, et al. KIF1Bβ functions as a haploinsufficient tumor suppressor gene mapped to chromosome 1p36.2 by inducing apoptotic cell death. J Biol Chem. 2008;283:24426–34.
Chen ZX, Wallis K, Fell SM, Sobrado VR, Hemmer MC, et al. RNA helicase a is a downstream mediator of KIF1Bβ tumor-suppressor function in neuroblastoma. Cancer Discov. 2014;4:434–51.
Sharina IG, Martin E, Thomas A, Uray KL, Murad F. CCAAT-binding factor regulates expression of the beta-1 subunit of soluble guanylyl cyclase gene in the BE2 human neuroblastoma cell lines. Proc Natl Acad Sci U S A. 2003;100:11523–8.
Boily G, Ouellet S, Langlois S, Lariviere M, Drouin R, Sinnett D. In vivo footprinting of the Glypican 3 (GPC3) promoter region in neuroblastoma cells. Biochim Biophys Acta. 2007;1769:182–93.
Yan YX, Zhao JX, Han S, Zhou NJ, Jia ZQ, Yao SJ, et al. Tetramethylpyrazine induces SH-SY5Y cell differentiation toward the neuronal phenotype through activation of the PI3K/Akt/Sp1/Topo11β pathway. Eur J Cell Biol. 2015;94:626–41.
Benfante R, Antonini RA, Vaccari M, Flora A, Chen F, Clementi F, et al. The expression of the human neuronal α3 Na+, K+-ATPase subunit gene is regulated by the activity of the Sp1 and NF-Y transcription factors. Biochem J. 2005;386:63–72.
Cappabianca L, Farina AR, Tacconelli A, Mantovani R, Gulino A, Mackay AR. Reconstitution of TIMP-2 expression in SH-SY5Y neuroblastoma cells by 5-azacytidine is mediated transcriptionally by NF-Y through an inverted CCAAT site. Exp Cell Res. 2003;286:209–18.
Bernard J, Raguenez G, Kauffmann A, Valent A, Ripoche H, Joulin V, et al. MYCN-non-amplified metastatic neuroblastoma with good prognosis and regression: a molecular portrait of stage 4S. Mol Oncol. 2008;2:261–71.
Cattaneo E, Conti L. Generation and characterization of embryonic striatal conditionally immortalized ST14A cells. J Neurosci Res. 1998;53:223–34.
Graham FL, Smiley J, Russel WC, Nairn R. Characteristics of a human cell line transformed by DNA from human adenovirus 5. J Gen Virol. 1977;36:59–74.
Stepanenko AA, Dmitrenko VV. HEK293 in cell biology and cancer research: phenotype, karyotype, tumorigenicity, and stress-induced genome-phenotype evolution. Gene. 2015;569:182–90.
Ruggeri P, Farina AR, Cappabianca L, Di Ianni N, Ragone M, Merolla S, et al. The TrkAIII oncoprotein inhibits mitochondrial free radical ROS-induced death of SH-SY5Y neuroblastoma cells by augmenting SOD2 expression and activity at the mitochondria, within the context of a tumor stem cell-like phenotype. PLOS ONE. 2014;15:e94568.
Roder K, Wolf SS, Larkin KJ, Schweizer M. Interaction between the two ubiquitously expressed transcription factors NF-Y and Sp1. Gene. 1999;234:61–9.
Emami KH, Smale ST, Burke TW. Sp1 activation of TATA-less promoter requires a species-specific interaction involving transcription factor IID. Nucl Acids Res. 1998;26:839–46.
Degterev A, Hitomi J, Germscheid M, Ch’en IL, Korkina O, Teng X, et al. Identification of RIP1 kinase as a specific cellular target of necrostatins. Nat Chem Biol. 2008;4:313–21.
PDQ Pediatric Treatment Editorial Board. Neuroblastoma Treatment (PDQ®): Health Professional Version. 2018 Aug 17. In: PDQ Cancer Information Summaries [Internet]. Bethesda (MD): National Cancer Institute (US); 2002-.https://www.ncbi.nlm.nih.gov/books/NBK65747/.
Fernholm M. On the development of the sympathetic chain and the adrenal medulla in the mouse. Z Anat Entwicki Gesch. 1971;133:305–17.
Mesner PW, Winters TR, Green SH. Nerve growth factor withdrawal-induced cell death in neuronal PC12 cells resembles that in sympathetic neurons. J Cell Biol. 1992;119:1669–80.
El Shamy WM, Linnarsson S, Lee K-F, Jaenisch R, Ernfors P. Prenatal and postnatal requirements for sympathetic neuroblast survival and innervation of specific targets. Development. 1996;122:491–500.
Pinon LGP, Minichiello L, Klein R, Davies AM. Timing of neuronal death in trkA, trkB and trkC mutant embryos reveals developmental changes in sensory neuron dependence on Trk sigaling. Development. 1996;122:3255–61.
Lipscomb EA, Sarmiere PD, Crowder RJ, Freeman S. Expression of the SM-20 gene promotes death in nerve growth factor-dependent sympathetic neurons. J Neurochem. 1999;73:429–32.
Zhao C, Takita J, Tanaka Y, Setou M, Nakagawara T, Tekada S, et al. Charcot-Marie-tooth disease type 2A caused by mutation in a microtubule motor KIF1Bbeta. Cell. 2001;105:587–97.
Sommer L, Rao M. Neural stem cells and regulation of cell number. Prog Neurobiol. 2002;66:1–18.
Lee S, Nakamura E, Yang H, Wei W, Linggi MS, Sajan MP, et al. Neuronal apoptosis linked to EglN3 prolyl hydroxylase and familial pheochromocytoma genes: developmental culling and cancer. Cancer Cell. 2005;8:155–67.
Stoer A, Yang Y, Patel S, Evangelista AM, Aponte A, Wang G, et al. Prolyl hydroxylation regulates protein degradation, synthesis, and splicing in human induced pluripotent stem cell-derived cardiomyocytes. Cardiovasc Res. 2016;110:346–58.
Sperandio S, de Belle I, Bredesen DE. An alternative, nonapoptotic form of programmed cell death. Proc Natl Acad Sci U S A. 2000;97:14376–81.
Sperandio S, Poksay K, de Belle I, Lafuente MJ, Liu B, Nasir J, et al. Paraptosis: mediation by MAP kinases and inhibition by AIP-1/Alix. Cell Death Differ. 2004;11:1066–75.
Dey A, Mustafi SB, Saha S, Kumar Dhar Dwivedi S, Mukherjee P, Battacharya R. Inhibition of Bmi1 induces autophagy-mediated necroptosis. Autophagy. 2016;12:659–70.
Fricker M, Tolkovsky AM, Borutaite V, Coleman M, Brown GC. Neuronal cell death. Physiol Rev. 2018;98:813–80.
Schlisio S. Neuronal apoptosis by prolyl hydroxylation: implication in nervous system tumours and the Warburg conundrum. J Cell Mol Med. 2009;3:4104–12.
Nazio F, Bordi M, Cianfanelli V, Locatelli F, Cecconi F. Autophagy and cancer stem cells: molecular mechanisms and therapeutic applications. Cell Death Differ. 2019;26:690–702.
Hoa N, Myers MP, Douglass TG, Zhang JG, Delgado C, Driggers L, et al. Molecular mechanisms of paraptosis induction: implications for genetically modified tumor vaccine. PLoS One. 2009;4:e4631.
Blackstone C. KIF1Bβ and neuroblastoma; failure to divide and cull. Dev Cell. 2016;36:127–8.
Angelina C, Tan ISY, Choo Z, Lee OZJ, Perviaz S, Chen ZX. KIF1Bβ increases ROS to mediate apoptosis and reinforces its protein expression through O2- in a positive feedback mechanism in neuroblastoma. Sci Reports. 2017;7:16867.
Han H, Chou C-C, Li R, Liu J, Zhang L, Zhu W, et al. Chalcomoracin is a potent anticancer agent acting through triggering oxidative stress via a mitophagy- and paraptosis-dependent mechanism. Sci Reports. 2018;8:9566. https://doi.org/10.1038/s41598-018-27724-3.
Toh WH, Logette E, Corcos L, Sabapathy K. Tap73β and DNp73β activate the expression of the prosurvival cappase-2s. Nuc Acids Res. 2006;36:4498–509.
Van Maerken T, Rihani A, Dreidax D, De Clercq S, Yigit N, Marine J-C, et al. Functional analysis of the p53 pathway in neuroblastoma cells using a small-molecule MDM2 antagonist nutlin-3. Mol Cancer Ther. 2011;10:983–93.
Miyazaki M, Otomo R, Matsushima-Hibiya Y, Suzuki H, Nakajima A, Abe N, et al. The p53 activator overcomes resistance to ALK inhibitors by regulating p53-target selectivity in ALK-driven neuroblastomas. Cell Death Disc. 2018;5:56. https://doi.org/10.1038/s41420-018-0059-0.
Mantovani R, Li XY, Pessara U, Hooft van Huisjduijnen R, Benoist C, Mathis D. Dominant negative analogues of NF-YA. J Biol Chem. 1994;269:20340–6.
Pilar G, Landmesser L. Ultrastructural differences during embryonic cell death in normal and peripherally deprived ciliary ganglia. J Cell Biol. 1976;68:339–56.
Kutcher LM, Shaham S. Non-apoptotic cell death in animal development. Cell Death Differ. 2017;24:1326–36.
Chen J, Crutchley J, Zhang D, Owzar K, Kastan MB. Identification of a DNA damage-induced alternative splicing pathway that regulates p53 and cellular senescence markers. Cancer Disc. 2017;7:766–81.
Genesha VR, Schor NF. Pharmacological management of high-risk neuroblastoma in children. Paediatr Drugs. 2011;13:245–55.
Liang F, Schaufele F, Gardner DG. Functional interaction of NF-Y and Sp-1 is required for type a natriuretic peptide receptor gene transcription. J Biol Chem. 2001;276:1516–22.
Yamada K, Tanaka T, Miyamoto K, Noguchi T. Sp family members and nuclear factor-Y cooperatively stimulate transcription from the rat pyruvate kinase M gene distal promoter region via their direct interactions. J Biol Chem. 2000;275:18129–37.
Wang Y, Wan C, Yu S, Yang L, Li B, Lu T, et al. Upregulating expression of NF-YC contributes to neuronal apoptosis via proapoptotic protein bim in rats’ brain hippocampus following middle cerebral artery occlusion (MCAO). J Mol Neurosci. 2014;52:552–65.
Ma H, Yue X, Gao L, Liang X, Yan W, Zhang Z, et al. ZHX2 enhances the cytotoxicity of chemotherapeutic drugs in liver tumor cells by repressing MDR1 via interfering with NF-YA. Oncotarget. 2015;6:1049–63.
Kawata H, Yamada K, Shou Z, Mizutani T, Yazawa T, Yoshino M, et al. Zinc-fingers and homeoboxes (ZHX) 2, a novel member of the ZHX family functions as a transcription factor. Biochem J. 2003;373:747–57.
Hirano S, Yamada K, Kawata H, Shou Z, Mizutani T, Yazawa T, et al. Rat zinc-fingers and homeoboxes 1 (ZHX1), a nuclear factor-YA-interacting nuclear protein, forms a homodimer. Gene. 2002;290:107–14.
Yamada K, Osawa H, Granner DK. Identification of proteins that interact with NF-YA. FEBS Lett. 1999;460:41–5.
Ceribelli M, Benatti P, Imbriano C, Mantovani R. NF-YC complexity is generated by dual promoters and alternative splicing. J Biol Chem. 2009;284:34189–200.
Dolfini D, Zambelli F, Pedrazzoli M, Mantovani R, Pavesi G. A high definition look at the NF-Y regulome reveals genome-wide associations with selected transcription factors. Nucl Acid Res. 2016;44:4684–702.
Valk-Lingbeek ME, Bruggeman SWM, van Lohuizen M. Stem cells and cancer: the polycomb connection. Cell. 2004;118:409–18.
Molofsky AV, Pardal R, Iwashita T, Park IK, Clarke MF, Morrison SJ. Bmi-1 dependence distinguishes neural stem cell self-renewal from progenitor proliferation. Nature. 2003;425:962–7.
Molofsky AV, He S, Bydon M, Morrison SJ, Pardal R. Bmi-1 promotes neural stem cell self-renewal and neural development but not mouse growth and survival by repressing the p16Ink4a and p19Arf senescence pathways. Genes Dev. 2005;19:1432–7.
Fasano CA, Dimos JT, Ivanova NB, Lowry N, Lemischka IR, Temple S. shRNA knockdown of Bmi-1 reveals a critical role for p21-Rb pathway in NSC self-renewal during development. Cell Stem Cell. 2007;1:87–99.
Lessard J, Baban S, Sauvageau G. Stage-specific expression of polycomb group genes in human bone marrow cells. Blood. 1998;91:1216–24.
Srinivasan M, Bharali DJ, Sudha T, Khedr M, Guest I, Sell S, et al. Downregulation of Bmi1 in breast cancer stem cells suppresses tumor growth and proliferation. Oncotarget. 2017;8:38731–42.
Nakayama K, Gazdoiu S, Abraham R, Pan ZQ, Ronai Z. Hypoxia-induced assembly of prolyly hydroxylase PHD3 into complexes: implications for its activity and susceptibility for degradation by the E3 ligase Siah2. Biochem J. 2007;401:217–26.
Zhou T, Erber L, Liu B, Gao Y, Ruan H-B, Chen Y. Proteomic analysis reveals diverse proline hydroxylation-mediated oxygen-sensing cellular pathways in cancer cells. Oncotarget. 2016;7:79154–69.
Appelhoff RJ, Tian YM, Raval RR, Turley H, Harris AL, Pugh CW, et al. Differential function of the prolyl hydroxylases PHD1, PHD2, and PHD3 in the regulation of hypoxia inducible factor. J Biol Chem. 2004;279:38458–65.
Mestdagh P, Bostrom A-K, Impens F, Fredlung E, Van Peer G, De Antonellis P, et al. The miR-17-92 microRNA cluster regulates multiple components of the TGFβ-pathway in neuroblastoma. Mol Cell. 2010;40:762–73.
Ji M, Rao E, Ramachandrareddy H, Shen Y, Jiang C, Chen J, et al. The miR-17-92 MicroRNA cluster is regulated by multiple mechanisms in B-cell malignancies. J Am Pathol. 2011;179:1645–56.
Narva E, Pursiheimo J-P, Laiho A, Rahkonen N, Emani MR, Viitala M, et al. Continuous hypoxic culturing of human embryonic stem cells enhances SSEA-3 and MYC levels. PLoS One. 2013;8:e78847.
Tuckerman JR, Zhao Y, Hewitson KS, Tian YM, Pugh CW, Ratcliff PJ, et al. Determination and comparison of specific activity of the HIF-prolyl hydroxylases. FEBS Lett. 2004;576:145–50.
Achilles A, Trainor PA. Neural crest stem cells: discovery, properties and potential for therapy. Cell Res. 2012;22:288–304.
Biagiotti T, D’Amico M, Marzi I, Di Gennaro P, Arcangeli A, Wanke E, Oliviotto M. Cell renewing in neuroblastoma: electrophysiological and immunocytochemical characterization of stem cells and derivatives. Stem Cells. 2006;24:443–53.
Ross RA, Spengler BA. Human neuroblastoma stem cells. Semin Cancer Biol. 2007;17:241–7.
Li S, Fell SM, Surova O, Smedler E, Wallis K, Chen ZX, et al. The 1p36 tumor suppressor KIF1Bβ is required for calcineurin activation controlling mitochondrial fission and apoptosis. Dev Cell. 2016;36:164–78.
Pandian V, Ramraj S, Khan FH, Azim T, Aravindan N. Metastatic neuroblastoma cancer stem cells exhibit flexible plasticity and adaptive stemness signaling. Stem Cell Res Ther. 2015;6:2. https://doi.org/10.1186/s13287-015-0002-8.
Garner EF, Beierle EA. Cancer stem cells and their interaction with the tumour microenvironment in neuroblastoma. Cancers. 2016; 8:doi:https://doi.org/10.3390/cancers8010005.
Zheng X, Naiditch J, Czurylo M, Jie C, Lautz T, Clark S, et al. Differential effect of long term drug selection with doxorubicin and virinostat on neuroblastoma cells with cancer cell characteristics. Cell Death Dis. 2013;4:e740.
Lucia Cappabianca and Antonietta R. Farina are co-first authors. We wish to thank Professor R. Mantovani (University of Milan, Italy) for the NF-YAdn expression vector and Professor E. Cattaneo (University of Milan, Italy) for the ST14A cell line.
This work was funded by Fondazione Salvatore Maugeri (CTERZMAUGERI11) and AIRC grant (IG20801).
Ethics approval and consent to participate
Written consent was obtained for the use of all patient NB RNA samples by the tissue bank of the Italian Association of Haematology and Pediatric Oncology (Genova, Italy). Animals experiments were performed at University of Rome (La Sapienza) in accordance with Italian National and Rome university regulations.
Consent for publication
No identifying patient details are contained within this manuscript.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Cappabianca, L., Farina, A.R., Di Marcotullio, L. et al. Discovery, characterization and potential roles of a novel NF-YAx splice variant in human neuroblastoma. J Exp Clin Cancer Res 38, 482 (2019) doi:10.1186/s13046-019-1481-8
- Alternative splicing
- Cancer stem cells
- Genotoxic stress