1. Introduction
Neuroblastoma (NB) is one of the most common solid tumors occurring in children. This tumor is divided into favorable and unfavorable biological categories, with the unfavorable tumors being more common in older children and characterized by rapid progress and chemo-resistance. Favorable tumors usually occur in infants, and some even regress or mature spontaneously. Since the malignant grade of NB biology greatly influences patient outcomes, biological factors related to the malignant grade of NB are of great clinical interest but are not yet fully understood.
In our department, we performed an exhaustive genetic analysis of biologically favorable and unfavorable NB tumors derived from 16 favorable cases who were alive with remaining residual tumors, including 4 spontaneous regressed-stage MS cases, and from clinical 16 unfavorable cases who died of tumor progression, and extracted genes with significantly higher expression levels in favorable tumors with oligo-microarrays [
1,
2,
3]. The top three highly expressed genes in favorable tumors were
DHRS3,
NROB1 (also known as
DAX1), and
CYP26A1, which are correlated with functions involving differentiation of the nervous system and cell division arrest [
4,
5,
6]. These genes are activated in neuronal development and differentiation. Therefore, expression levels of these genes might mediate the regression or maturation of NB cells.
NB cell lines were established from an aggressive neuroblastoma case (
Supplementary Table S1). In this study, we transfected each of these three genes into NB cell lines (SK-N-SH, NH12, and TGW, which have different characteristics) using the plasmid method of expression. To knockdown these genes, siRNA for each gene was also transfected into these cells. Subsequently, cell morphology was assessed, and assays for growth, cell cycle, migration, and colony formation to monitor anchorage-independent growth were performed. Gene expression profiles were further examined by RNA sequencing. Transfection of these highly expressed genes reduced cell growth, which consequently made it difficult to establish a stable cell-induced differentiation assay. To rule out the effect of the transfection process of these genes, we also used the retroviral gene stable expression system known as the Retro-X Tet-On Advanced System. To assess the correlation with NB differentiation induced by all-trans retinoic acid (ATRA), cells in which the Tet-On expression system was used were selected.
2. Materials and Methods
2.1. Cell Culture
TGW and SK-N-SH cells were obtained from the Japanese Cancer Resource Cell Bank (
http://cellbank.nihs.go.jp/, accessed on 7 March 2003) and the NH-12 cell line was established in our institution. TGW cells were cultured in Dulbecco’s Modified Eagle Medium, and SK-N-SH and NH12 cells were cultured in MEM-alpha (GIBCO BRL, Rockville, MD, USA) containing 2 mM L-glutamine. These media were supplemented with 10% fetal bovine serum (FBS) (Labome, Princeton, NJ, USA) and 100 µg/mL penicillin/streptomycin. Passage numbers of these cell lines were largely equal.
2.2. Plasmid Construction and Recombinant Retrovirus
To generate recombinant plasmid constructs, we amplified full-length Homo sapiens dehydrogenase/reductase 3 (DHRS3, RefSeqID, NM_004753, transcription variant 1) cDNA, Homo sapiens cytochrome P450 family 26 subfamily A member 1 (CYP26A1, RefSeqID, NM_000783 transcription variant 1) cDNA, and Homo sapiens nuclear receptor subfamily 0 group B member 1 (NROB1, RefSeqID, NM_000475) cDNA by PCR, using the following primer sets:
DHRS3-F: 5′- ATATGGTGTGGAAACGGCTG-3′, DHRS3-R: 5′-ATCCTATGTCCGCCCTTTGAA-3′; CYP26A1-F: 5′- ATATGGGGCTCCCGGCGCT-3′, CYP26A1-R: 5′- TCAGATTTCCCCATGGAAATGGGT-3′; NROB1-F: 5′-ATATGGCGGGCGAGAACCACCA-3′, NROB1-R: 5′- TTATATCTTTGTACAGAGCATTTCCAGCA-3′. These PCR products were cloned and inserted into a pIRES2 DsRed-Express2 vector (Clontech Lab, Takara bio, Mountain View, CA, USA), and then a 3xFLAG tag was inserted at the N-terminal of each gene.
These products were also used in the Retro-X™ Tet-On® Advanced Inducible Expression System (Clontech) which was applied to generate retrovirus constructs. Then, we established stable expression of each retroviral candidate gene for NB by expressing genes that activate transcription of downstream genes in the presence of doxycycline (Dox). We then confirmed expression following induction with Dox and examined changes following induction. The HEK-293-based retroviral packaging cell line GP2-293 was used for packaging retroviruses. Induction of gene expression was achieved by adding Dox (Clontech). Sequences were checked using the PRISM 3130 sequencer (ABI, Thermo Fisher Scientific, Waltham, MA, USA).
The day before transfection, culture cells were trypsinized and counted. The cells were seeded in a 24-well plate to achieve 50~70% confluency on the day of transfection in 100 µL of Serum free Opti-MEM® I Reduced Serum Medium with 500 ng of plasmid DNA and then 0.5 μL of PLUS Reagent was added. After incubation for 5 min at room temperature (r.t.), the cells were mixed thoroughly with Lipofectamine LTX and then incubated for 30 min at r.t. These cells were divided into each well and incubated for 6 hrs at 37 °C with 5% CO2. The selection process using G418 and Blasticidin, as described below, started from the next day. The selected clones were confirmed for expression of transfected genes using Western blotting and real-time quantitative RT-PCR (qRT-PCR), as described below.
2.3. In Vitro Cell-Migration Assay
An in vitro cell-migration assay was performed using a culture insert (ibidi, Gewerbehof, Germany). Cells (2 × 10
5) were seeded in the two wells of a culture insert, and after the cells adhered, the insert for producing the gap was removed. After 0 h, 24 h, and 48 h, the migration of cells into the gap area was analyzed by ImageJ software [
7].
2.4. Soft-Agar Colony-Formation Assay
We subjected cells to the soft-agar colony-formation assay in 35 mm dishes; 2000 cells were seeded into 1.5 mL of pre-autoclaved 0.8% low-melting agarose in filtrated 2× DMEM or a-Medium supplied with 20% FBS, which was overlaid with a 1.5 mL layer of 0.4% agar in the 2× appropriate medium. After 3 weeks, the colonies were analyzed and counted, and, after fixation, they were stained with crystal violet. Image-Pro software (Media Cybernetics, Rockville, MD, USA) was used for this assay.
2.5. Cell-Proliferation Assay
Cells were seeded (2 × 105) into 35 mm dishes and cultured in DMEM/alpha-MEM supplemented with 10% FBS and then harvested and counted using the LUNA automatic-cell-counting system (Logos Biosystems, Annandale, VA, USA) every two days. Proliferation was assessed in triplicate, and each experiment was repeated three times.
2.6. Cell-Cycle Analysis
Cell cycles were analyzed using a flow cytometer (LSRFortessa X-20TM, Becton Dickinson Biosciences, Franklin Lakes, NJ, USA). Cells were harvested and rinsed with PBS and suspended in PBS with 1% bovine serum albumin (BSA). They were then stained with Hochest 33342. The staining solution was incubated for 30 min in the dark and the cells were then washed with PBS three times. Nuclear DNA was analyzed by flow cytometry, and the population of cells in each cycle was estimated.
2.7. Western Blot Analysis, Immunocytochemistry, and Electron Microscopy
Cells were harvested and whole-cell lysates were extracted using a reagent kit (cOmplete Lysis-M, Roche, Basel, Switzerland), according to the manufacturer’s protocol, and acetone protein precipitation was then performed. After centrifugation at 15,000× g for 10 min at 4 °C, debris was discarded, and protein concentrations were measured. SDS–polyacrylamide gel electrophoresis (12%) was then carried out. For these analyses, 20 μg of protein was added to each lane. Anti-FLAG M2 (Sigma, St. Louis, MO, USA) was used as the primary antibody, and horseradish-peroxidase-conjugated goat anti-mouse antibodies were used as the secondary antibodies (GE Healthcare, Chicago, IL, USA). Western blot signals were detected with the chemiluminescence substrate of the ECL™ Prime Western Blotting detection system (GE Healthcare). Immunohistochemistry was assessed as follows: NB cells were seeded on type-1-collagen-coated glass cover slips. After incubation, cells were fixed by 4% paraformaldehyde in PBS and permeabilized by 0.2% Triton X-100, and they were then washed with PBS, blocked by 10% normal goat serum (Thermo Fisher Scientific, Waltham, MA, USA), and incubated with the primary antibody mixture in 1% normal goat serum for 1 hr at room temperature. After the incubation and PBS washing steps, samples were covered with secondary antibody staining solution and incubated for 1 h at room temperature. Samples were then washed several times with PBS and visualized by fluorescence microscopy (BZ-X710, Keyence, Osaka, Japan), confocal microscopy (Fluoview 1000-D, Olympus, Tokyo, Japan), and super-resolution microscopy (DeltaVision OMX 3D-SIM, GE Healthcare, Houston, TX, USA). Expression of exogenous proteins was detected with FLAG M2 antibodies and their respective endogenous antibodies, as described below. Forced expression strains were also fat-stained using Lipid Tox. For the detection of endogenous protein expression, we used the following primary antibodies: mouse anti-DHRS3 (H00009249-B01P, Abnova, Taipei, Taiwan), rabbit anti-DHRS3 (NBP1-80846, NOVUS, Tokyo, Japan), mouse anti-CYP26A1 (sc-53618, Santa Cruz Biotechnology, Inc., Dallas, TX, USA), mouse anti-DAX1 (MABD398, Millipore, Burlington, MA, USA), and rabbit anti-DAX-1 (sc-841, Santa Cruz Biotechnology, Inc.). Mouse anti-β-actin (ab8226, Abcam, Cambridge, UK) and rabbit anti-β-actin (Abcam, ab8227) were also used to confirm loading amounts.
For the analysis of intracellular changes using electron microscopy, the culture cells were fixed with 3% glutaric dialdehyde/0.1 M PBS (pH 7.2), followed by 1% osmium tetraoxide (OsO4/0.1 M PBS), and stained with 3% uranyl acetate. The cells were absolutely dehydrated with an ascending series of ethanol solutions and embedded in Epon 812 resin. After confirming the resin cure, thin slices were cut in the vertical direction to form cell layers. The sections were laid on a copper grid and stained with 2% uranil acetate and lead citrate. Electron-microscopical examination was performed using a transmission electron microscope (JEM-1230, JEOL, Tokyo, Japan) under the conditions of an 80 kv accelerating voltage and a 20,000-magnification ratio.
2.8. Establishment of NB Cell Lines Expressing Favorable Candidate Factors
Transient and stable transfection in NB cell lines were assessed using the lipofection method (Lipofectamine LTX, Invitrogen, Thermo Fisher Scientific), according to the manufacturer’s protocol. After three weeks of selection with G418 and Blasticidin, the BD FACS Aria II (Becton Dickinson) flow cytometer was used to sort DsRed and EmGFP expression in NB cells. Isolated colonies were expanded and established as stable polyclonal or clonal clones.
2.9. Real-Time Quantitative RT-PCR (qRT-PCR)
qRT-PCR was performed using the Power Up SYBR™ Green RT qPCR kit (Thermo Fisher Scientific). RT-qPCR was run on the Bio-Rad CFX96 Real-Time PCR Detection System. The thermal cycling conditions were 50 °C for 2 min and 95 °C for 2 min, followed by 40 cycles at 95 °C for 1 s and 60 °C for 30 s. A melting-curve analysis was recorded after each cycle. Cycle thresholds (Cts) were calculated using Bio-Rad CFX Manager software. Relative gene expressions were calculated using the 2−ΔΔCT method.
The primer sequences used were as follows (ACTB was used as the internal control): ACTB-F: 5′-TACCTCTATGCCAACACAGT-3′, ACTB-R: 5′- AGTACTTGCGCTCAGGAGGA-3′; DHRS3-F: 5′- TGGTCCATGGGAAGAGC-CTA-3′, DHRS3-R: 5′- TGTCCGCCCTTTGAAAGTGT-3′; CYP26A1-F: 5′-CTCCTCCAAATGGAATGAAG-3′, CYP26A1-R: 5′-CAGTCAGCATGGATGATATG-3′; NROB1-F: 5′- CAGTCAGCATGGATGATATG-3′, NROB1-R, 5′-CACAGCTCTTTATTCTTCCC-3′.
2.10. RNA Sequencing and Pathway Analysis
Comprehensive RNA-expression analysis was performed in DsRed- and EGFP-positive NB cells, as sorted by the FACS Aria cell sorter (Becton Dickinson Bioscience, Sparks, MD, USA). Total RNA was extracted from cell lines using the QIAwave RNA Mini Kit (Qiagen, Hilden, Germany). Total cellular RNA was extracted from tumor tissues by the acid-guanidinium–phenol–chloroform method [
8]. The extracted RNA was quantified using a 2100 Bioanalyzer (Thermo Fisher Scientific, Santa Clara, CA, USA) with an Agilent RNA 6000 Nano kit.
RNA sequencing was performed as follows. Sequencing libraries were generated using the TruSeq RNA LT Sample Prep Kit (Illumina, SanDiego, CA, USA), according to the manufacturer’s instructions, with RNA normalized to inputs of 1 μg of total RNA from each sample. mRNA-sequence libraries were prepared from 500 ng of cDNA, according to the manufacturer’s protocol. Sequencing was performed for 101 cycles using the HiSeq Reagent Kit v3 on a HiSeq
® 2500 instrument (Illumina). The next-generation sequencing platforms for analysis and data visualization applied belonged to the CLC Genomics Workbench (Qiagen). After read genome annotation from Fast-q files, quality trimming, PCR-duplicate-mapped-reads removal, and filtering of reads, the transcripts were processed. The BAM file data were exported to Strand NGS v2.1 (Agilent Scientific Instruments, Santa Clara, CA, USA), and raw transcript counts were normalized by cell line with the parameters from DESeq2 software [
9]. Then, expression analysis, gene ontology (GO) analysis, and pathway analysis were carried out. Ingenuity pathway analysis by IPA (TOMY digital biology, Tokyo, Japan) was applied for further investigation of biological pathway interactions.
2.11. Statistical Analysis
Unless otherwise specified, Fisher’s exact test was used for p-value calculations between two categorical variables. The Wilcoxon rank-sum and Kruskal–Wallis tests were used to analyze differences between two or more continuous variables. Log-rank tests were used to compare survival distributions in Kaplan–Meier plots. FDR corrections were performed for multiple testing corrections. Statistical analysis was carried out using the software programs RStudio, JMP v13, and Strand NGS. Statistical significance was defined at p < 0.05 by the Student’s t-test and Dunnett’s test. Results are presented as means ± SEs.
4. Discussion
The prognosis of patients with NB depends on the tumor biology and the stage of disease [
10,
11,
12]. To better understand NB biology, many studies have been conducted and several factors related to highly malignant NB have been identified, including
MYCN amplification and telomerase activation [
13,
14,
15,
16,
17]. However, few basic studies have assessed which factors characterize favorable tumors. Therefore, we conducted the present study, in which we have shown that
DHRS3,
CYP26A1, and
NROB1 are potential favorable prognostic factors for NB, as the expression of these genes in NB cell lines suppressed growth and reduced anchorage independency in plasmid-established lines.
In the cells transfected with one of these genes, morphological changes were largely similar in these NB cell lines, and our findings suggest that overexpression of these genes induces cellular senescence. Morphologically, senescence-like changes, such as enlargement and flattening of cells, were observed in the clones after transfection. In the growth assay, overexpression of these genes reduced growth in most cell lines. Finally, soft-agar formation suggested that overexpression of these genes reduced growth and cellular differentiation via a lower level of stemness activity.
Transfected DHRS3 accumulated in the ER and LDs in NH12 and SK-N-SH cells but was absent in TGW cells. In the growth assay (
Figure 4A), DHRS3-transfected SK-N-SH and NH12 cells showed reductions in growth, but growth was not changed in TGW cells. In the soft-agar colony-formation assay, TGWs showed little change in size and formation efficiency. Since TP53 promotes the accumulation of DHRS3 in LD and ER [
18], this was not seen in TGWs harboring a mutated TP53 gene. The TP53 mutation might thus be one of the genetic aberrations that are characteristic of highly malignant NB, and TP53 is necessary for this growth-reduction pathway through DHRS3. On the other hand, NROB1-transfected TGW cells showed severe gene dropout and low levels of protein expression, indicating the existence of a TP53-independent pathway of cellular senescence by NROB1.
Gene expression analysis by RNA sequencing revealed that pathways related to cellular senescence and cell adhesion were related. Pathway analysis of significantly up- and down-regulated genes found a reduction in several pathways, including Vitamin A and Carotenoid Metabolism, GPCR Ligand Binding, and WNT Ligand Biogenesis and Trafficking, implicating cellular senescence, cell adhesion and differentiation, and nervous system development.
DHRS3 is a member of the classical short-chain dehydrogenase reductase superfamily [
19], which was originally cloned from retinal cells. Retinol-like compounds affect intracellular processes mainly by regulating gene expression, including cell proliferation and differentiation [
20]. They are expressed in several non-ocular human tissues, including the liver, pancreas, heart, kidney, and lung. Although DHRS3 was found to be up-regulated in papillary thyroid cancer and neuroblasts, DHRS3 was negatively correlated with the metastasis of papillary thyroid carcinoma and was associated with a better prognosis for NB [
3,
21]. This is consistent with the data for DHRS3 overexpression in our study. Nonetheless, further research is necessary to understand the specific role of DHRS3 in NB. Notably, in a sarcoma study, it was found that a four-gene signature, including DHRS3, JRK, TARDBP, and TTC3, can predict patient survival based on a real-world cohort. These genes may be more feasible for clinical application and should be validated prospectively.
NROB1 (also known as DAX1) is a critical mediator of retinoic-acid (RA)-induced nNOS gene transcription in NB [
5]. Retinoids induce CYP26 enzyme production in the liver, enhancing their own rapid metabolic clearance, while retinoid resistance in tumor cells occurs in part due to increased CYP26 production. RA-metabolism-blocking agents, which inhibit CYP26 enzymes, can improve RA pharmacokinetics in preclinical NB models [
22]. The signaling networks through which RA mediates NB differentiation were inhibited upon
MYCN overexpression, which revealed opposing regulation of RA and for a number of differentiation-relevant genes, including LMO4, CYP26A1, ASCL1, RET, FZD7, and DKK1 [
23]. Furthermore, CYP26A1 was found to be a retinoid-responsive gene in retinoid-sensitive NB in
MYCN transgenic mice and lung and breast cancer cell lines [
24].
Based on the cellular phenomena and gene expression profiles we observed in the induction experiments on transfected NB cell lines, we conducted a differentiation-inducing experiment by administering ATRA to transfected cells. SK-N-SH and NH12 cells with DHRS3 expression became flattened and exhibited a giant morphology, and they strongly expressed senescence markers compared to cells without ATRA exposure. We observed little change in the expression of p16INK4a, a senescence-associated marker, and enhanced expression of p21Waf1/Cip1, suggesting that DHRS3 may be involved in senescence and differentiation.
Dropout of the transduced genes occurred during culture passage in transfected NB cell lines, making long-term stable gene transfer difficult. Thus, having a tetracycline-inducible expression system with retroviral vectors resulted in more pronounced cell-morphology changes. Though the expression of CYP26A1 did not affect the cell-growth rate or cell cycle in the Dox-expressing system, DHRS3 and NROB1 expression led to cell-morphology changes and inhibition and arrest of cell growth in a short period of time. Notably, DHRS3 expression arrested the cell cycle by engaging with the all-trans-retinol pathway and undergoing differentiation and senescence. Since p53 promotes LD accumulation in a manner consistent with DHRS3 enrichment in the ER [
18], TGW cells were resistant to DHRS3 induction due to the TP53 gene mutation they harbor. On the other hand, CYP26A1 induction in the Dox-expressing system showed morphological changes without growth inhibition (
Figure 4,
Supplementary Figure S3), suggesting that CYP26A1, as a retinoid-inducible enzyme, might play a role as a cell-differentiation-inducing factor instead of inhibiting growth arrest—a matter which requires further analysis.
Candidate genes for favorable prognostic factors may reduce the ability of NB cell lines to suppress tumor growth and anchorage dependency. To examine the usefulness of these findings as diagnostic markers of malignancy, we will next examine clinical samples from a bank of intractable tumors. These findings will offer insights into molecular targets and the feasibility of differentiation-induction therapies for the treatment of NB.