Abstract
Background/Aim: Overall survival for the high-risk group of neuroblastoma (NB) patients still remains at 40-50%, necessitating the establishment of a curable treatment. LIM domain only 1 (LMO1) gene encoding a transcriptional regulator is an NB-susceptibility gene with a tumor-promoting activity. Previously we conducted chromatin immunoprecipitation and DNA sequencing analyses on NB cell lines and identified 3 protein-coding genes regulated by LMO1. In this study, we extended our analyses to capture microRNA genes directly or indirectly regulated by LMO1. Materials and Methods: Using microarrays, we conducted a comparative gene expression analysis on an NB cell line SK-N-SH; between the cells with and without LMO1 suppression. Results: Overall, 18 microRNAs were identified to be indirectly down-regulated by LMO1 including 7 microRNAs of the let-7 family, whose cell proliferation inhibitory activity was observed. Conclusion: Target genes of the LMO1-regulated microRNAs and their relevant pathways may be a potential therapeutic target.
Neuroblastoma (NB), the most common extracranial solid tumor in childhood, is a cause for approximately 15% of all paediatric cancer deaths (1). With recent progress in NB treatment, patient prognosis has been significantly improved; Overall survival (OS) for low-risk and intermediate-risk groups is now >98% and 90-95%, respectively. On the contrary, OS for the high-risk (HR) group is still at 40-50%, and 50-60% of these patients have a relapse, for which no option of curable treatment remains (2, 3). To improve the OS and quality of life for the HR group, the identification of novel therapeutic target molecules is imperative (4-6).
Previously, the LIM domain only 1 (LMO1) gene encoding a transcriptional regulator was identified as an NB susceptibility gene by a genome-wide association study (7). In the study, it was also revealed that both the presence of an NB-related allele of single-nucleotide polymorphism (SNP) rs110419 in the LMO1 gene and the increased copy number of the gene contribute to augmented LMO1 expression, and that LMO1 promotes proliferation of NB cells. Importantly, it was shown that the relation of the variation at rs110419 in LMO1 as well as the increased copy number to NB is stronger in the HR group than in the non-HR group. Recently, the SNP rs110419 was found to be located in an enhancer region of the gene and have an influence on its function (8). Consequently, it is likely that the transcriptional regulator LMO1 is involved in the regulation of expression of a variety of genes important for NB progression in HR patients, and it is anticipated that identification of the LMO1's regulatory targets may lead us to identify novel tumor-promoting molecular pathways in NB, as well as novel therapeutic target molecules, especially for treatment of the HR group.
Previously, we conducted chromatin immunoprecipitation and DNA sequencing analyses on two NB cell lines, SK-N-SH and LA-N-5 and identified 3 protein-coding genes regulated by LMO1 in the 2 NB cells: LIM and senescent cell antigen-like domains 1 (LIMS1), Ras suppressor protein 1 (RSU1) and relaxin 2 (RLN2) (9). LIMS1 and RSU1 encode proteins functioning with integrin-linked kinase (ILK) in integrin adhesomes, and an ILK inhibiting compound Cpd22 inhibited proliferation of the NB cells in vitro (9-12). The findings suggest the possibility that in addition to the LIMS1/ILK pathway, more tumor-promoting pathways exist downstream of the LMO1-regulatory cascade.
As the previous study aiming at LMO1-regulated genes was limited to protein-coding genes, analyses were extended to microRNA genes in this study. MicroRNAs are known to have an important role in regulating gene function by binding to the 3’ untranslated region (UTR) of target messenger RNAs (mRNAs) and suppressing their translation (13, 14). Therefore, it is likely that the mRNAs targeted by LMO1-regulated microRNAs are significantly involved in NB oncogenesis.
Materials and Methods
Cell line. The NB cell line SK-N-SH was provided from the European Collection of Cell Cultures and maintained in D-MEM (Wako Pure Chemical Industries, Ltd., Osaka, Japan) supplemented with 10% fetal calf serum. LA-N-5 was provided by the RIKEN BRC through the National Bio-Resource Project of the MEXT, Japan, and maintained in RPMI1640 (Life Technologies Japan Ltd., Tokyo, Japan) supplemented with 10% FCS. SK-N-BE (2) was obtained from ECACC and maintained in a 1:1 mixture of EMEM (Wako) and Ham's F12 (Wako) supplemented with 0.5% non-essential amino acids (Life Technologies) and 15% FCS.
shRNA. Each shRNA construct was prepared by annealing the (a) strand to the (b) as follows: for shLMO1 (target sequence: 5’-GGCATTGGACAAGTACTGG-3’), (a) 5’-tcgagGGtATTGGAtAAGTAtTGGttcaagagaCCAGTACTTGTCCAATGCC-ttttttacgcgta-3’ and (b) 5’-agcttacgcgtaaaaaaGGCATTGGAGAAGTACTGGtctcttgaa-CCAaTACTTaTCCAATaCCc-3’, for shGFP (target sequence: 5’-GCTACCTGTTCCATGGCCAA-3’), (a) 5’-tcgagGCTAtCTGTTCCgTGGCCgAttcaagagaTTGGCCATGGAACAGGTAGttttttacgcgtg-3’ and (b) 5’-gatccacgcgtaaaaaaCTACCTGTTCCATGGCCAAtctcttgaaTcGGCCAcGGAACAGaTAGCc-3’. The shRNA was inserted into pLVSIN-CMV neo (Takara Bio, Shiga, Japan) whose CMV promoter was replaced by hU6 promoter. To obtain viral particles, the pLVSIN-hU6-shRNA constructs were introduced into Lenti-X™ 293T Cells (Takara Bio) with Lenti-X™ HTX Packaging System (Takara Bio), and after 72 h’ incubation, the medium was collected and its viral titer (infection units/ml) was determined by transduction to HT1080 cells. The NB cell lines were transducted with lentivirus (> 100,000 infection units) in the presence of polybrene (10 μg/ml in culture medium, Sigma-Aldrich, St. Louis, MO, USA).
Microarray expression analyses. Total RNA was extracted from SK-N-SH transduced with pLVSIN-shLMO1 or pLVSIN-shGFP, using Absolutely RNA miRNA Kit (Agilent Technologies, Santa Clara, CA, USA). Microarray expression analyses were performed with miRNA Complete Labeling and Hyb Kit (Agilent) and SurePrint G3 Human miRNA Microarray, Release 19.0 (Agilent), following a standard protocol recommended by the manufacturer. The data were analyzed using Agilent Feature Extraction Software (Agilent).
Real-time RT-PCR. For microRNA hsa-miR-3648, RT-PCR was conducted by converting about 10 ng of total RNA to the first strand cDNA with TaqMan MicroRNA Reverse Transcription Kit (Life Technologies), followed by performing TaqMan MicroRNA Assay (Life Technologies, Applied Biosystems Assay ID: 464401). For mRNA, RT-PCR was performed by converting approximately 5 μg of total RNA to the first strand cDNA with High Capacity cDNA Reverse Transcription Kit (Life Technologies), followed by conducting TaqMan Gene Expression Assay (Life Technologies, Applied Biosystems Assay ID: Hs00231133_m1 for LMO1 and Pre-Developed TaqMan Assay Reagent for glyceraldehyde-3-phosphate dehydrogenase (GAPDH)). The PCR was performed for 40 cycles that included 2 steps of temperature: 95°C for 15 sec and 60°C for 60 sec, using ABI PRISM 7900HT Sequence Detection System (Life Technologies). The relative transcript level was calculated using the Ct value of GAPDH transcript as reference.
Down-regulation of GAPDH and western blot analyses. Control siRNA and siRNA for GAPDH suppression were purchased from BIONEER Corp, Daejeon, Korea, and introduced into NB cells using Lipofectamine RNAiMAX (Life Technologies). Using CelLytic-M Mammalian Cell Lysis/Extraction Reagent (Sigma-Aldrich) and Protease Inhibitor Cocktail (Sigma-Aldrich), cell lysates were prepared and loaded onto a 15% SDS-polyacrylamide gel. Immune blotting was performed with an anti-GAPDH, an anti-Actin and a HRP-conjugated anti-goat antibody (Santa Cruz Biotechnology, Dallas, Texas, USA). The signals were detected by Pierce Western Blotting Substrate Plus (Thermo Fisher Scientific, Yokohama, Japan).
Cell growth assay. A mirVana miRNA inhibitor for miR-3648 and a negative control were obtained from Life Technologies. The microRNA mimics and Negative control #2 were purchased from BIONEER Corp. The NB cells were prepared in five wells of 96-well plates (10,000 cells/well). The microRNA inhibitor and mimics were added to the wells (6 pmol/well) using Lipofectamine RNAiMAX (Life Technologies). The cells were incubated for 72 h and the growth was evaluated by MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl tetrazolium bromide) assay using Cell Counting Kit-8 (Dojindo, Kumamoto, Japan).
Results
MicroRNA hsa-miR-3648 is a direct target of LMO1 regulation. In our previous study, we performed chromatin immunoprecipitation-DNA sequencing (ChIP-Seq) to identify LMO1-regulated genes using 2 NB cell lines SK-N-SH and LA-N-5 (9). From the sequenced reads, those mapped to a reference genome sequence were selected for peak detection (in other word, for detection of candidates for LMO1-binding site) and applied to 6 peak-detection programs with 9 different settings. A definite peak was defined as the peaks detected by more than 5 of the 9 settings of the programs in each cell line (9). In this study, we focused on the definite peaks located in or near known microRNA genes (within a 5-kb proximity of known genes), common to the two cell lines, and identified 5 microRNA genes; hsa-mir-548ap, hsa-mir-198, hsa-mir-4461, hsa-mir-3648 and hsa-mir-3687 (Table I). We examined the effect of LMO1 suppression on these genes by performing microRNA microarray expression analyses, comparing the SK-N-SH NB cells in which LMO1 was suppressed with shRNA, to its counterpart introduced with shRNA for a GFP gene as a reference (Figure 1). We listed the microRNAs that were down-regulated to half or less in NB cells following LMO1 suppression, compared to reference (Table II). The hsa-mir-3648 (Table I) was one of the listed microRNAs and its down-regulation was confirmed by real-time RT-PCR in both SK-N-SH and LA-N-5 cells where LMO1 expression was suppressed (Figure 2). Assuming that hsa-mir-3648 is an oncogenic microRNA up-regulated by LMO1, the consequence of inhibiting the microRNA in NB cells was examined. In addition to SK-N-SH and LA-N-5 cell lines, the SK-N-BE(2) cell line was used in the inhibition assay. The reason is that SK-N-SH and LA-N-5 have a low siRNA transfection efficiency and SK-N-BE(2) is more permissive for siRNA transfection (Figure 3). An oligonucleotide inhibitor for hsa-mir-3648 was introduced into the three NB cell lines and showed a very weak but statistically significant suppressive effect on LA-N-5, however, it did not suppress cell proliferation of the other two cell lines. This suggests that hsa-mir-3648 may not be related to the cell-proliferation promoting activity of LMO1 or may have just a limited role in the LMO1 activity (Figure 4).
MicroRNAs of the let-7 family are indirectly regulated by LMO1. Next, we focused on microRNAs which were up-regulated in SK-N-SH cells, in which LMO1 is suppressed. A total of 18 microRNAs were found to be up-regulated by 2- or more fold in SK-N-SH cells following suppression of LMO1 (Table III). As none of the 18 genes is a regulatory target of LMO1 identified by the ChIP-Seq analyses, all of them are likely to be under an indirect regulation by LMO1. Importantly, they included hsa-miR-34a-5p and let-7 family microRNAs which were previously reported to suppress the NB oncogene MYCN (15-17). We examined the effects of the 8 microRNA, hsa-miR-34a-5p, hsa-let-7a-5p, hsa-let-7b-5p, hsa-let-7c, hsa-let-7e-5p, hsa-let-7f-5p, hsa-let-7g-5p and hsa-let-7i-5p, on cell proliferation by introducing their oligonucleotide mimic into NB cell lines. As a result, all the mimics showed a statistically significant inhibitory effect on the proliferation of, at least, one of the three cell lines (Figure 5). It is noteworthy that the inhibitory effect of the mimics was clearly demonstrated especially in SK-N-BE(2), which display a much higher transfection efficiency for siRNA than SK-N-SH or LA-N-5.
Finally, the effect of the other 10 of the 18 microRNAs on cell proliferation was examined. Mimics for three microRNAs, hsa-miR-16-5p, hsa-miR-20b-5p and hsa-miR-409-3p significantly suppressed cell growth of SK-N-BE(2) but not that of SK-N-SH or LA-N-5, which may reflect their differences in transfection efficiency (Figure 6).
Discussion
In this study, hsa-miR-3648 was identified as a direct target of LMO1, however, no cell proliferation inhibitory activity of the microRNA was demonstrated in NB cells. There is no previous report on the functions of hsa-miR-3648. A recent report stated that hsa-miR-3648 is transcriptionally regulated in neural stem cells by the amyloid precursor protein (APP)-derived cleavage product, APP intracellular domain (AICD) (18). Taken together with our results indicating that hsa-miR-3648 is directly regulated by LMO1, we can suggest that AICD may cooperatively regulate the expression of hsa-miR-3648 with LMO1 in NB cells.
The MYCN gene is the best-known oncogene that functions in NB and promotes cell proliferation, cell invasion and angiogenesis by regulating a variety of molecules (17). In the let-7 family, let-7e was demonstrated to bind to the 3’ untranslated region of the MYCN transcript and inhibit its function, conferring a tumor suppressor activity in NB (15). Humans have 10 mature let-7 family members; let-7a, let-7b, let-7c, let-7d, let-7e, let-7f, let-7g, let-7i, mir-98 and mir-202 (19). Also, a tumor suppressive function of some let-7 family members was demonstrated in non-NB cancers. Stable transfection of let-7f reduced growth of thyroid cancer cells through inhibition of mitogen-activated protein kinase (20). It was reported that proliferation of lung cancer cells is suppressed by let-7a, let-7b, let-7c, let-7d and let-7g (21). Proliferation of nasopharyngeal carcinoma cell lines, HK1 and HONE1, was inhibited by let-7a, let-7b, let-7d, let-7e, let-7g and let-7i (22). By targeting MYH9 (Myosin, Heavy Chain 9, Non-Muscle) transcripts, let-7f inhibits invasion and metastasis of gastric cancer cells (23). A malignant phenotype of osteosarcoma cell lines, U2-OS and HOS, is changed by let-7g through inhibition of Aurora kinase B gene (24).
LIN28B is an RNA binding protein that binds to the terminal loop of let-7 microRNAs and inhibits processing of the microRNA (25). As silencing of LIN28B in NB cells up-regulates expression of let-7a, let-7e, let-7g, let-7i, and mir-98, it is likely that these microRNAs are regulated by LIN28B (26). A genome wide association study revealed LIN28B as an NB susceptibility gene, and Lin28b-transgenic mice develop NB (26). Intriguingly, Lin28 is down-regulated by let-7 and both molecules form a reciprocal autoregulatory loop that controls microRNA processing during embryonic neural stem cell commitment in mice (27).
In addition to the let-7 family, miR-34a was reported to target MYCN transcripts (16). It was also reported as a tumor suppressor that inhibits proliferation of NB cells through targeting mitogen-activated protein kinase kinase kinase 9 (MAP3K9) transcripts (28, 29).
Besides the 7 microRNAs of the let-7 family and miR-34a, this study detected 10 microRNAs indirectly suppressed by LMO1. Among these, miR-20b, miR-16 and miR-409 were clearly demonstrated to have a cell proliferation activity on SK-N-BE(2) (Table III, Figure 6). However, their targets in normal adrenal tissue or NB remain to be elucidated. Up-regulation of miR-20b expression was reported in gastric, breast and colorectal cancers and it is thought to act as an oncomir (30-32). However, an inhibitory activity of miR-20b on proliferation, migration and invasion of bladder cancer cells was also demonstrated (33). It may have a role in both tumor progression and suppression in a cancer type-dependent manner. It was reported that cisplatin down-regulates the expression of brain-derived neurotropic factor, a member of the neurotropin family associated with aggressive malignant behavior by targeting miR-16 and resulting in inhibition of cell proliferation of the NB cell line SH-SY5Y (34). miR-409 suppresses cancer cell invasion and metastasis by targeting a pro-metastatic gene radixin in gastric cancers (35).
The results of this study suggest that LMO1 executes its oncogenic function in NB, partly by down-regulating several tumor suppressive microRNAs.
Acknowledgements
This work was supported by the National Cancer Center Research and Development Fund (No. 23-B-18) and by JSPS KAKENHI Grant Number 26670510.
Footnotes
This article is freely accessible online.
- Received July 8, 2018.
- Revision received July 19, 2018.
- Accepted July 23, 2018.
- Copyright© 2018, International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved