BRIEF RESEARCH REPORT article
LncRNA H19 Suppresses Osteosarcomagenesis by Regulating snoRNAs and DNA Repair Protein Complexes
- 1Department of Integrative Biology and Pharmacology, McGovern Medical School, The University of Texas Health Science Center at Houston, Houston, TX, United States
- 2Department of Obstetrics and Gynecology and Women's Health, Einstein/Montefiore Medical Center, Bronx, NY, United States
- 3Department of Neurology, Renaissance School of Medicine at Stony Brook University, Stony Brook, NY, United States
- 4Department of Cell, Developmental and Regenerative Biology, Icahn School of Medicine at Mount Sinai, New York, NY, United States
- 5The Black Family Stem Cell Institute, Icahn School of Medicine at Mount Sinai, New York, NY, United States
- 6The Graduate School of Biomedical Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, United States
- 7Institute of Molecular Biology, National Chung Hsing University, Taichung, Taiwan
- 8Graduate Institute of Chinese Medical Science, China Medical University, Taichung, Taiwan
- 9The University of Texas MD Anderson Cancer Center UTHealth Graduate School of Biomedical Sciences, Houston, TX, United States
- 10Center for Stem Cell and Regenerative Medicine, The Brown Foundation Institute of Molecular Medicine for the Prevention of Human Diseases, The University of Texas Health Science Center at Houston, Houston, TX, United States
- 11Center for Precision Health, School of Biomedical Informatics and School of Public Health, The University of Texas Health Science Center at Houston, Houston, TX, United States
Osteosarcoma is one of the most frequent common primary malignant tumors in childhood and adolescence. Long non-coding RNAs (lncRNAs) have been reported to regulate the initiation and progression of tumors. However, the exact molecular mechanisms involving lncRNA in osteosarcomagenesis remain largely unknown. Li-Fraumeni syndrome (LFS) is a familial cancer syndrome caused by germline p53 mutation. We investigated the tumor suppressor function of lncRNA H19 in LFS-associated osteosarcoma. Analyzing H19-induced transcriptome alterations in LFS induced pluripotent stem cell (iPSC)-derived osteoblasts, we unexpectedly discovered a large group of snoRNAs whose expression was significantly affected by H19. We identified SNORA7A among the H19-suppressed snoRNAs. SNORA7A restoration impairs H19-mediated osteogenesis and tumor suppression, indicating an oncogenic role of SNORA7A. TCGA analysis indicated that SNORA7A expression is associated with activation of oncogenic signaling and poor survival in cancer patients. Using an optimized streptavidin-binding RNA aptamer designed from H19 lncRNA, we revealed that H19-tethered protein complexes include proteins critical for DNA damage response and repair, confirming H19's tumor suppressor role. In summary, our findings demonstrate a critical role of H19-modulated SNORA7A expression in LFS-associated osteosarcomas.
Li-Fraumeni syndrome (LFS) (OMIM #151623) is a rare familial cancer syndrome characterized by early onset of various tumors, soft-tissue sarcomas, osteosarcomas, breast cancers, brain tumors, adrenocortical carcinomas, and leukemia (Zhou et al., 2017). Germline mutations in the p53 tumor suppressor gene are responsible for LFS. Patient-derived iPSCs have been used to model various diseases (Lee et al., 2009b; Carvajal-Vergara et al., 2010; Itzhaki et al., 2011; Yagi et al., 2011; Mulero-Navarro et al., 2015; Gingold et al., 2016; Liu et al., 2018; Zhu et al., 2018). The development of various refined differentiation protocols utilizing induced pluripotent stem cells (iPSCs) has enabled the production of large quantities of differentiated cells from individual patients. In our previous studies, using a LFS patient-derived iPSC model, we demonstrated the tumor suppressor role of lncRNA H19 and the oncogenic role of SFRP2 during the formation of osteosarcoma in LFS patients (Lee et al., 2015; Lin et al., 2017b; Kim et al., 2018; Zhou et al., 2018). Although H19-mediated tumor suppression has been demonstrated in LFS-associated osteosarcomagenesis, the underlying mechanisms of its tumor suppressor activity remain unclear.
Long non-coding RNAs (lncRNAs) represent a class of nonprotein-coding RNAs longer than 200 nucleotides (Sanchez Calle et al., 2018). Despite not being translated into proteins, lncRNAs are important regulators of diverse biological processes and pathologies. LncRNAs are proposed to function via multiple mechanisms, including co-transcriptional regulation, bridging proteins and chromatin, offering cytoplasmic scaffolding, pairing with other RNAs, and serving as molecular decoys (Ulitsky and Bartel, 2013).
One of the most studied lncRNAs is the imprinted lncRNA H19, located on human chromosome 11 and expressed exclusively from the maternal allele (Gabory et al., 2009). Numerous functional studies have assessed the role of H19 in the pathogenesis of human cancers and yielded conflicting results. On one hand, H19 is a precursor of miR-675 (Keniry et al., 2012) as well as a “molecular sponge” for soaking up microRNAs let-7 (Kallen et al., 2013; Li et al., 2020) and miR-138 (Liang et al., 2015), supporting a role in promoting cancer cell proliferation and migration. On the other hand, in vivo mouse (Yoshimizu et al., 2008) and LFS iPSC-derived osteoblast (Lee et al., 2015) studies demonstrated that H19 displays a tumor-suppressive effect. Evidence of numerous human tumors displaying either overexpression or lack of expression of H19 (Hao et al., 1993; Lustig-Yariv et al., 1997; Yoshimizu et al., 2008; Lee et al., 2015; Li et al., 2020) suggests the possibility of context-dependent oncogenic and tumor-suppresive roles.
Small nucleolar RNAs (snoRNAs) are small non-coding RNAs 60–300 nucleotides in length, primarily located in the nucleolus. The main functions of snoRNAs are to assist with post-transcriptional modification and maturation of ribosomal RNAs (rRNAs), small nuclear RNA (snRNAs), and other cellular RNAs. SnoRNAs can interact with RNA binding proteins to form small nuclear ribonucleoproteins. SnoRNAs are divided into two classes according to their catalytic activity: C/D box snoRNAs, catalyzing 2-O-ribose methylation; and H/ACA box snoRNAs, catalyzing pseudouridylation (Bachellerie et al., 2002). Emerging evidence indicates that snoRNAs are widely involved in various cancer-related signaling pathways. For example, SNORA42 is overexpressed in non-small cell lung cancer (NSCLC) and plays an oncogenic role through suppressing p53 function and/or expression (Mei et al., 2012). SNORD76 leads to the activation of the WNT/β-Catenin pathway to promote hepatocellular carcinoma (HCC) tumorigenicity (Wu et al., 2018).
The DNA damage/repair response plays a key role in maintaining genome integrity and stability, and its dysfunction girds the development and progression of various cancer types. Recently, studies revealed that lncRNA also regulates DNA damage response and DNA repair through the ATM/ATR (Wan et al., 2013), p53 regulatory network (Hu et al., 2018), and DNA double-strand break repair pathway (Gazy et al., 2015). However, it remains unclear whether H19 is actively involved in DNA damage/repair response.
Our current study reveals that H19 functions as a tumor suppressor through negatively regulating the expression of oncogenic SNORA7A in the osteoblast context. Furthermore, by exploring its interactions with multiple essential DNA damage/repair response factors, we demonstrate how H19 executes its tumor-suppressive role.
Materials and Methods
U2OS cells were maintained in DMEM (Invitrogen, USA) supplemented with 10% FBS and 1% penicillin/streptomycin. iPSCs were maintained on Matrigel (Corning, USA)-coated plates in StemMACS™ iPS-Brew XF medium (MiltenyiBiotec, USA). All cells were cultured in a humidified incubator at 37°C and 5% CO2 and were tested to exclude mycoplasma contamination.
Differentiation of iPSCs to Mesenchymal Stem Cells (MSCs) and Then to Osteoblasts
In vitro differentiation of LFS iPSCs to MSCs was performed by a PDGF-AB-based method described previously (Lian et al., 2007; Zhou et al., 2018). Appropriately differentiated CD105+/CD166+/CD24- MSCs were plated in a 6-well plate at a density of 2 × 104 cells per well and cultured in an osteogenic differentiation medium (α-MEM supplemented with 10% FBS, 10 mM β-glycerol phosphate, 200 μM ascorbic acid, and 0.1 μM dexamethasone) to induce osteogenic differentiation as previously described (Lee et al., 2015).
RNA Isolation and RNA-Sequencing
LFS osteoblasts with ectopic H19 expression and vector control osteoblasts were collected and compared for H19-mediated tumor suppressor effects. Total mRNA was isolated using TRIzol reagent (Thermo Fisher Scientific, USA) according to the manufacturer's instructions. All RNA sample preparation and RNA-sequencing (RNA-seq) data analyses were performed as described previously (Lee et al., 2015). RNA was aggregated from biological triplicate experiments prior to running as a single sample for RNA-seq. The FPKM (fragments per kilobase of exon model per million reads mapped) of genes, lncRNAs, snoRNAs, and miRNAs were calculated and summarized in Supplementary Table 1.
GO Biological Process (GO_BP) and Wikipathway analyses were performed using Enrichr (https://amp.pharm.mssm.edu/Enrichr/) (Chen et al., 2013; Kuleshov et al., 2016) to identify enriched biological processes and pathways in H19-expressing LFS osteoblasts. Processes and pathways enriched with a p < 0.05 were considered significant. These genes were analyzed in Enrichr using Wikipathway and GO_BP pathway datasets to identify enriched GO terms, pathways, and functions. The TRANSFAC and JASPAR PWMs databases of transcription factors were used in Enrichr to identify transcription factors positively and negatively correlated with SNORA7A expression.
To determine pathological functions in osteosarcomagenesis significantly associated with SNORA7A, we identified genes positively and negatively correlated with SNORA7A expression in the TCGA-SARC dataset from SNORic (snoRNA in cancers; http://bioinfo.life.hust.edu.cn/SNORic/basic/) (Gong et al., 2017). Correlation of SNORA7A and mRNA expression of significant genes was run at default settings. Clinical analysis of 5-year survival in different cancers was analyzed from TCGA.
A modified 1x streptavidin-binding RNA aptamer (S1m) (Leppek and Stoecklin, 2014) was synthesized (Integrated DNA Technologies, USA) from DNA oligo pairs (S1m_Forward: AATTGgtagaaaATGCGGCCGCCGACCAGAATCATGCAAGTGCGTAAGATAGTCGCGGGTCGGCGGCCGCATctgctgggG; S1m_Reverse: AATTCcccagcagATGCGGCCGCCGACCCGCGACTATCTTACGCACTTGCATGATTCTGGTCGGCGGCCGCATtttctacC), annealed, and ligated into the multiple cloning sites of the tetracycline-inducible TetO-FUW vector (Lee et al., 2015) 4 times sequentially. The TetO-FUW vector containing 4×streptavidin-binding RNA aptamer was named TetO-FUW-4S1m. H19 was then cloned into the N terminal side of 4S1m to form TetO-FUW-4S1m-H19. For the SNORA7A expression construct, the synthesized full-length SNORA7A (Ensembl Transcript ID: ENST00000384765.1) was cloned into the pLKO.pig plasmid (Lee et al., 2012a,b) within EcoRI-AgeI cloning sites and confirmed by Sanger sequencing.
Total RNA was isolated using TRIzol reagent (Thermo Fisher Scientific, USA) and snoRNA was isolated using mirVana miRNA isolation kit (Thermo Fisher Scientific, USA) according to the manufacturer's instructions. The qRT-PCR primers for GAPDH (internal control), BGLAP, MEPE, FGF23, H19, and SNORA7A were described previously (Lee et al., 2015; Zhang et al., 2017).
Virus Packaging and Infection
TetO-FUW-4S1m or TetO-FUW-4S1m-H19 plasmids were co-transfected with packaging plasmids psPAX2 (Addgene, plasmid # 12260) and pMD2.G (Addgene, plasmid # 12259) into HEK-293T cells, and the virus collected from the cell culture medium 48 hours later. Osteosarcoma cell line U2OS was infected with viral particles together with the M2rtTA virus (Addgene, plasmid # 20342), which was similarly produced in HEK-293T cells, in the presence of 8ug/ml polybrene (Sigma-Aldrich, USA). Thirty-six hours post-infection, U2OS cells were treated with 1μg/ml doxycycline (Sigma-Aldrich, USA) for 24 hours to induce the expression of 4S1m or 4S1m-H9.
In vitro Anchorage-Independent Growth (AIG) Assay
AIG assay was performed as described previously (Lee et al., 2015). Briefly, LFS MSCs were transduced with H19 and/or SNORA7A. 10,000 LFS MSCs were mixed with 0.5% UltraPure low-melting-point agarose (Thermo Fisher Scientific, USA) and cultured in an osteogenic differentiation medium for 1 month. Colonies (≥50 μm) were counted under a Leica microscope DMi8.
H19-Interacting Protein Complex Purification and Mass Spectrometry Analysis
H19-interacting proteins were purified via the streptavidin-binding aptamers S1m as described previously (Leppek and Stoecklin, 2014). Briefly, S1m-H19 expressing U2OS cells were resuspended in 500 μl ice-cold lysis buffer (20 mM Tris-HCl (pH 7.5), 150 mM NaCl, 1.5 mM MgCl2, 2 mM DTT, 2 mM RNase inhibitor, and complete protease inhibitors cocktail). Lysates were subjected to centrifugation for 5 min at 12,000 rpm at 4°C. The supernatants were incubated with streptavidin agarose beads (Thermo Fisher Scientific, USA) overnight at 4°C under rotation. The streptavidin agarose beads were then washed five times for 5 min at 4°C with wash buffer [20 mM Tris-HCl (pH 7.5), 300 mM NaCl, 5 mMMgCl2, and 2 mM DTT]. The 4S1m-H19 associated protein complex samples were subjected to SDS-PAGE electrophoresis and micro-liquid chromatography/tandem mass spectrometry as described previously (Lee et al., 2007, 2009a, 2012a). Mass spectrometry analysis results are summarized in Supplementary Table 2.
Results were presented as mean ± standard error of the mean (SEM). Error bars in figures represent SEM. Differences between two groups were examined by the two-tailed unpaired Student t-test. *p < 0.05; **p < 0.01; and ***p < 0.001.
The data supporting the findings are available within the manuscript text, figures, and Supplementary Tables 1, 2. The RNA-seq data are available at the sequencing read archive (SRA) under accession number PRJNA673185.
H19 Modulates snoRNA Expression in LFS Osteoblasts
To explore the suppressive effects of H19 on osteosarcomagenesis, we analyzed the genome-wide transcriptomes of LFS iPSC-derived osteoblasts by RNA-seq following ectopic overexpression of H19. The RNA-seq experiment was performed with either vector- or H19-expressing LFS osteoblasts (n = 1) after pooling RNA from biological triplicate experiments. Transcriptome analysis confirmed ectopic H19 overexpression and identified a set of significantly altered genes [107 upregulated genes (FPKM ≥1 and fold change ≥5) and 81 downregulated genes (FPKM ≥1 and fold change ≥5)]. These H19-regulated genes were analyzed for enriched Gene Ontology biological processes (GO_BP) and Wikipathways using the comprehensive gene set enrichment analysis web tool Enrichr. GO_BP analyses showed that the GO_BPs enriched in H19-expressing LFS osteoblasts compared with vector-expressing (H19-depleted) LFS osteoblasts included ketone pathway, regulation of complement activation, skeletal system development, and positive regulation of cell death. In contrast, H19-depleted LFS osteoblasts demonstrated enrichment of GO_BPs including pre-mRNA cleavage required for polyadenylation, negative regulation of the cellular process, negative regulation of transport, and extracellular matrix organization (Figure 1A). Wikipathway analyses indicated that H19-overexpressing LFS osteoblasts were enriched for genes in the complement and coagulation cascade, complement activation, oxidative damage, prostaglandin synthesis, and regulation, TYROBP causal network, benzopyrene metabolism, apoptosis, and CCK2R signaling. In contrast, genes enriched in H19-depleted LFS osteoblasts were involving in methylation pathways, biogenic amine synthesis, estrogen metabolism, and hypertrophy model (Figure 1B).
Figure 1. H19 significantly affects miRNA and snoRNA expression. (A,B) Heatmaps show the GO_BP (A) and Wikipathway (B) categories enriched in either H19-expressing LFS osteoblasts (right column) or vector control LFS osteoblasts (left column) compared with each other based on mRNA expression of respective samples. Colors represent p values of enrichment. (C-E) H19 overexpression induces small changes in overall mRNA expression but large changes in miRNA and snoRNA expression in LFS osteoblasts. The percentage of mRNAs, lncRNAs, etc. (C), miRNAs (D), and snoRNAs (E) with different fold-changes (<2, 2-5, and >5 increase or decrease) are shown in bar plots. RNA expression between H19-overexpressing and vector control cells is presented in a scatter plot (right). Dot colors are used to represent the fold change of RNA expression. ns (not significant) indicates mRNAs, lncRNAs, miRNAs, snoRNAs, etc. with fold-change <2 and/or FPKM<1. The RNA-seq experiment was done with one sample pooled from 3 biological replicates. miRNAs and snoRNAs are defined based on the miRbase and snoDB database, respectively. “mRNA, lncRNA, etc.” indicates all transcripts except miRNAs and snoRNAs.
In analyzing the transcriptome data, we classified RNA alterations into three groups by their fold-change: less than 2, between 2 and 5, and higher than 5. We unexpectedly found that a much larger share of miRNAs and snoRNAs increased in expression following H19 expression in LFS osteoblasts compared with protein-coding mRNAs, lncRNAs, etc. (i.e all transcripts except miRNAs and snoRNAs). H19 restoration in LFS osteoblasts led to a more than 5-fold change (increased or decreased) in expression in 0.96% of protein-coding mRNAs, lncRNAs, etc. with the majority of mRNAs and lncRNAs changing <2-fold (FPKM <1) (Figure 1C). In contrast, 16% of miRNAs and 28.4% of snoRNAs changed greater than 5-fold in expression (Figures 1D,E). Scatter plots comparing miRNA and snoRNA expression between H19-expressing and vector control LFS osteoblasts showed that miRNAs (e.g., MIR570, MIR574, MIR943, MIR635, MIR7-1, and MIR145) and snoRNAs (e.g., SNORA6, SNORA7A, SNORA9, SNORA15B, SNORA19B, and SNORA36C) were significantly altered in H19-expressing LFS osteoblasts (Figures 1D,E). These findings suggest that H19 suppresses osteosarcomagenesis in LFS patients by primarily altering the expression of miRNAs and snoRNAs rather than mRNAs.
SnoRA7A Is an onco-snoRNA Involved in Tumor Progression and Associated With Poor Prognosis
Among the identified H19-regulated miRNAs and snoRNAs, SNORA7A was chosen for further study in light of its previously recognized role in controlling the self-renewal of MSCs (Zhang et al., 2017). To validate our transcriptome results, we ectopically expressed H19 in LFS iPSC-derived osteoblasts and found that H19 significantly inhibits SNORA7A expression (Figure 2A). H19-induced osteogenic gene expression (e.g., BGLAP, MEPE, and FGF23) was inhibited by SNORA7A (Figure 2B), indicating that SNORA7A inhibits osteogenesis. To investigate whether the inhibition of SNORA7A plays a role in H19-mediated tumor suppression, we performed an in vitro AIG assay and found retarded clonal growth of H19-transduced LFS osteoblasts upon SNORA7A ectopic expression (Figure 2C). These results suggest that SNORA7A is negatively regulated by H19, and that SNORA7A functions as an onco-snoRNA by antagonizing H19 tumor suppressor function.
Figure 2. Systems analyses of the oncogenic role of SNORA7A from TCGA datasets. (A) Ectopic expression of H19 downregulates SNORA7A expression in LFS osteoblasts. qRT-PCR data are represented as mean ±SEM; n = 3 biological replicates; statistical significance is determined using two-tailed Student's t-test; ***p < 0.001. (B) Left panel, ectopic expression of H19 increases osteogenic gene expression in LFS osteoblasts; in contrast, restoration of SNORA7A expression impairs H19-upregulated osteogenic gene expression in H19-transduced LFS osteoblasts. Right panel, qRT-PCR results demonstrate the expression of H19 and SNORA7A upon their ectopic expression for assays in the left panel. qRT-PCR data are represented as mean ± SEM; n = 3 biological replicates; statistical significance is determined using two-tailed Student's t-test; *p < 0.05; **p < 0.01; ***p < 0.001. (C) AIG assay for in vitro tumorigenicity demonstrates that H19 impairs the colony numbers of LFS osteoblasts and SNORA7A expression rescues the H19-suppressed tumorigenicity of LFS osteoblasts. H19 or H19/SNORA7A-transduced LFS osteoblasts were grown for 1 month and then assayed. Positive colonies are considered those larger than 50μm diameter. Data are represented as mean ± SEM; n = 6 biological replicates; statistical significance is determined using two-tailed Student's t-test; *p < 0.05; **p < 0.01; ***p < 0.001. (D,E) Genes positively (n = 780) and negatively (n = 945) correlated with SNORA7A expression in TCGA-SARC are identified. Pathway analysis by Wikipathway (A) and GO_BP (B) is performed on these gene sets using Enrichr to identify pathways significantly enriched or depleted (p ≤ 0.05) in association with SNORA7A expression. (F) Scatterplots of G6PD, RPS9, CDK4, CDKN1B, ADCY9, and ROCK1 mRNA expression correlation with SNORA7A expression. (G) Genes whose abnormal splicing (including exon skip, mutually exclusive splicing, or intron retention) by PSI was positively (n = 91, upper panel) or negatively (n = 182, lower panel) correlated with SNORA7A expression in sarcomas are identified by TRANSFAC and JASPAR PWWMs in Enrichr. PSI: percent spliced in index. (H) High SNORA7A expression is associated with poor cancer survival. Five-year overall survival is analyzed according to SNORA7A expression from TCGA in various cancers. SARC: sarcoma. BRCA: breast invasive carcinoma. ESCA: esophageal carcinoma. KIRC: kidney renal clear cell carcinoma.
To further explore the role of SNORA7A in tumorigenesis, we analyzed the TCGA-SARC dataset and identified 780 and 945 protein-coding genes whose mRNA levels were positively and negatively correlated with SNORA7A expression, respectively. Wikipathway analysis indicated that genes positively correlated with SNORA7A expression are mainly involved in cytoplasmic ribosomal protein function and cell metabolism (e.g., pentose phosphate metabolism, Cori cycle, pyrimidine metabolism, and glycolysis and gluconeogenesis), while genes negatively correlated with SNORA7A expression are associated with focal adhesion, integrin-mediated cell adhesion, and G protein signaling pathways, the latter of which is involved in osteoblast differentiation (Wu et al., 2010) (Figure 2D).
Similarly, GO_BP analysis demonstrated that genes positively correlated with SNORA7A expression are involved in rRNA metabolic process and processing, ribosome biogenesis, translation, glyceraldehyde-3-phosphate metabolic process, and ribose phosphate metabolic process, while genes negatively correlated with SNORA7A expression are associated with glycogen metabolic process and cytosolic calcium ion transport (Figure 2E). The molecules controlling pentose phosphate metabolism (G6PD), protein translation (RPS9), and cell cycle progression (CDK4) were positively correlated with SNORA7A expression. In contrast, cyclin-dependent kinase inhibitor (CDKN1B), adenylyl cyclase (ADCY9), and downstream effector of Rho (ROCK1) were negatively correlated with SNORA7A expression (Figure 2F). These systems analyses suggest that SNORA7A may transcriptionally and/or post-transcriptionally regulate multiple oncogenic features, including cellular metabolism, ribosome biogenesis, cell cycle, etc., culminating in osteosarcoma development in LFS patients.
The regulation of RNA splicing is increasingly recognized to be an essential mechanism underlying cancer development, and dysregulation of RNA splicing machinery has been found to contribute to tumorigenesis in various human cancers, including glioma, lymphoma, and breast cancer (David et al., 2010; Hsu et al., 2015; Koh et al., 2015). Given the well-recognized potential for snoRNAs to regulate pre-mRNA splicing (Falaleeva et al., 2016; Liang et al., 2019), we investigated the effects of SNORA7A levels toward potentially pathological mRNA splicing events including skipped exons, mutually exclusive alternative splicing, and intron retention. The percent spliced in index (PSI) was calculated for all identified genes in the TCGA-SARC dataset to identify genes whose alternative splicing was correlated with SNORA7A expression. These genes were then mapped to transcription factors using TRANSFAC and JASPAR analysis. Genes whose abnormal splicing was positively correlated to SNORA7A expression were mainly downstream targets of transcription factors such as the osteoblastic lineage regulator RUNX2 (Komori, 2019) and the tumor suppressors ZBTB7A (Liu et al., 2014) and ZFHX3 (Hu et al., 2019). Moreover, genes whose abnormal splicing was negatively correlated to SNORA7A expression were downstream targets of transcription factors such as the tumor angiogenesis regulator HIF1A (Koukourakis et al., 2002), the Hippo pathway transcription factor TEAD4 (Lin et al., 2017a; Shi et al., 2017) and the tumorigenic transcription factor NF-KB1 (Concetti and Wilson, 2018) (Figure 2G). These results imply that oncogenic effects of SNORA7A occur through regulation of RNA splicing of tumor suppressor genes and oncogenes.
We next analyzed the expression of SNORA7A in multiple cancers using TCGA datasets and correlated its expression with patient survival data. The five-year overall survival analyses showed that high SNORA7A levels were significantly associated with poorer survival in multiple cancer types including sarcoma (SARC, Log-rank p = 0.064, coxph = 0.017), breast cancer (BRCA, Log-rank p = 0.024, coxph = 0.026), esophageal carcinoma (ESCA, Log-rank p = 0.01, coxph = 0.0001), and kidney carcinoma (KIRC, Log-rank p ≤ 0.0001, coxph = 0.0001) (Figure 2H). Taken together, these results suggest a potential role of SNORA7A in promoting cancer progression and emphasize that H19 tumor suppressor functions may be mediated through the repression of SNORA7A expression.
H19 Potentially Interacts With DNA Damage/Repair Response Protein Complexes
Previous studies revealed that H19 is capable of interacting with EZH2 (Luo et al., 2013) to promote WNT/β-Catenin activation and subsequent downregulation of E-cadherin and that H19 also associates with KSRP (Giovarelli et al., 2014) to control mRNA decay. However, the protein complexes tethered by H19 in osteosarcoma remain largely unknown. To elucidate H19-associated protein complexes and gain insights into the potential biological functions controlled by the H19 interactome, we ligated H19 with an optimized 4 × streptavidin-binding RNA aptamer (S1m) and expressed H19-S1m in U2OS cells. H19-bound protein complexes could then be pulled down with streptavidin beads and identified by mass spectrometry (Supplementary Table 2). Compared with the S1m vector control, S1m-H19 pull-down identified 15 putatively interacting cellular proteins: TUBA1A, KHSRP, ALDH18A1, HIST1H1A, GANAB, TARS2, LRRFIP1, GPANK1, ZNF326, BLM, ABCC12, DMD, CLCN5, RBBP8, and OR8U1 (Figure 3A, left panel).
Figure 3. Potential H19 interacting proteins are primarily involved in DNA damage/repair response. (A) Left panel: Mass spectrometry analysis of TetO-FUW-4S1m-H19 streptavidin pull-down proteins identifies 15 unique H19-interacting proteins. TetO-FUW-4S1m vector is used as an S1m pull-down negative control (left). Right panel: The biological protein-protein interaction network based on H19-interacting proteins (orange, n = 15) and expanded to include genes (pink, n = 31) connected to at least two H19-interacting proteins is constructed in Cytoscape using the Biogrid human database. A solid line indicates physical interactions. A dashed line indicates genetic interactions. (B,C) Heatmaps representing p-values of enriched genesets (p ≤ 0.05) among the H19-interacting protein complex genes (n = 47) by Wikipathway (B) and GO_BP (C) analyses.
We applied the Biogrid human database to identify proteins previously demonstrated to bind to at least two of these 15 potential H19-binding proteins. This bioinformatics approach yielded 47 proteins strongly suspected to be connected to the 15 H19-associated proteins (Figure 3A, right panel). Wikipathway and GO_BP analyses revealed that these putative H19-associated protein complexes mainly function in DNA damage response, DNA repair, the ATM/ATR signaling pathway, the WNT signaling pathway, and the cell cycle (Figures 3B,C), emphasizing the potential role of H19 in regulating DNA damage response and repair. Particularly, BLM, an ATP-dependent DNA helicase, participates in DNA replication and DNA repair by modulating DNA double-strand break resection (Gravel et al., 2008). Mutations in the BLM gene are associated with Bloom syndrome, which carries a greatly increased risk of cancers including squamous cell carcinoma, leukemia, lymphoma, and gastrointestinal cancer (Lin et al., 2017b). Another H19-interacting protein is RBBP8 (also known as CtIP), an endonuclease that cooperates with the MRE11-RAD50-NBN (MRN) complex in DNA-end resection, the first step of homologous recombination (HR)-mediated double-strand break repair. RBBP8 mutations are commonly detected in various human cancer cell lines (You and Bailis, 2010). Inactivation of one RBBP8 allele predisposed mice to multiple types of cancers (e.g., lymphoma) suggesting that RBBP8 functions as a tumor suppressor (Chen et al., 2005). Based on these findings, we speculate that H19 tethers DNA repair proteins to other functional factors and helps maintain genomic integrity by suppressing oncogenic events such as DNA double-stranded breaks.
The essential biological functions of non-coding RNAs (e.g., miRNAs, snoRNAs, and lncRNAs) and their pathologic roles in tumorigenesis are becoming increasingly appreciated. In our previous work, we reported that H19, a lncRNA, is downregulated by mutant p53s (mutp53s) in LFS iPSC-derived osteoblasts and modulates osteoblastic differentiation and oncogenic repression through the imprinted gene network (IGN) (Gabory et al., 2009; Lee et al., 2015). While H19 has been demonstrated in various contexts to act as either a tumor suppressor or an oncogene, our current work further demonstrates the important regulatory role of H19 in osteosarcoma initiation and progression.
We discovered considerable H19-mediated regulation of snoRNAs. We confirmed that SNORA7A, a snoRNA known to promote MSC proliferation and self-renewal (Zhang et al., 2017), is suppressed by H19 and that restoration of SNORA7A expression impaired H19-mediated osteogenesis and tumor suppression. These findings suggest that SNORA7A functions as an onco-snoRNA in LFS-associated osteosarcoma. SNORNA7 expression in TCGA datasets is associated with numerous oncogenic pathways (ribosome biogenesis, pentose phosphate pathway, glycolysis, and HIF signaling) and RNA splicing events, suggesting a link between snoRNA expression and oncogenic events in clinical tumors.
p53 mutations have been widely linked to increased tumor oncogenesis through both gain-of-function interactions, for example by promoting HIF1A-regulated angiogenic genes (Amelio et al., 2018), and loss-of-function mutations, for example by abrogating wild-type p53 repression of FBL expression, culminating in increased ribosome biogenesis (Marcel et al., 2013). Upregulation of SNORA7A by mutp53-mediated H19 suppression provides yet another mechanism for mutp53 downstream pathways to cooperate in promoting bone malignancies.
In addition, mass spectrometry analysis of potential H19-interacting proteins revealed multiple DNA damage response and repair molecules directly associated with H19. H19 was previously proposed to regulate DNA damage response (Zheng et al., 2016; Zhu et al., 2017; Ma et al., 2018; Cheng et al., 2019), but the underlying mechanisms have yet to be defined. Our identification of potential interactions between H19 and BLM as well as RBBP8, both of which are recognized to regulate DNA damage response and repair, may provide another angle to elucidate these mechanisms. Importantly, gene mutations in both BLM and RBBP8 are associated with increased risks of tumor formation, including osteosarcoma. The evidence of physical interactions between H19 and DNA damage response genes leads us to speculate that H19 indirectly maintains genomic integrity, explaining its tumor-suppressor activity.
Finally, H19 expression was significantly increased in cancer cells treated with chemotherapeutic drugs (e.g., doxorubicin and cisplatin) (Zheng et al., 2016; Zhu et al., 2017). We speculate that upregulation of H19 following DNA damage events occurs in order to accelerate the repair process and prevent oncogenic effects of DNA damage.
In conclusion, our study indicated that lncRNA H19 plays a vital regulatory role in inhibiting osteosarcomagenesis and provides mechanistic insights for improving our understanding of H19-mediated tumor suppression in LFS patient-associated osteosarcomas.
Data Availability Statement
The transcriptome data has been deposited into the sequencing read archive (SRA; PRJNA673185).
AX, M-FH, DZ, BC, DW, C-CL, and RZ conducted all experiments. AX, M-FH, DZ, DW, C-CL, IL, RZ, and D-FL conceived and designed the study and interpreted results. AX, M-FH, DZ, JG, DB, RZ, and D-FL wrote the manuscript. All authors contributed to the article and approved the submitted version.
DZ was supported by DoD Horizon Award (W81XWH-20-1-0389). RZ was supported by the CPRIT UTHealth Cancer Genomic Core pilot grant (CGC-FY20-1). D-FL was a CPRIT Scholar in Cancer Research and is supported by CPRIT Award RR160019, the Rolanette and Berdon Lawrence Research Award, the Pablove Foundation Childhood Cancer Research Seed Grant Award, and CPRIT UTHealth Cancer Genomic Core pilot grant (CGC-FY20-1). We thank the Cancer Prevention and Research Institute of Texas (CPRIT RP180734) for supporting this work.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
In memory of our mentor and friend IL (1953-2017), who inspired us with his lifelong passion to understand stem cell biology.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2020.611823/full#supplementary-material
Amelio, I., Mancini, M., Petrova, V., Cairns, R. A., Vikhreva, P., Nicolai, S., et al. (2018). p53 mutants cooperate with HIF-1 in transcriptional regulation of extracellular matrix components to promote tumor progression. Proc. Natl. Acad. Sci. U.S.A 115, E10869–E10878. doi: 10.1073/pnas.1808314115
Carvajal-Vergara, X., Sevilla, A., D'Souza, S. L., Ang, Y. S., Schaniel, C., Lee, D. F., et al. (2010). Patient-specific induced pluripotent stem-cell-derived models of LEOPARD syndrome. Nature 465, 808–812. doi: 10.1038/nature09005
Chen, E. Y., Tan, C. M., Kou, Y., Duan, Q., Wang, Z., Meirelles, G. V., et al. (2013). Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinform. 14:128. doi: 10.1186/1471-2105-14-128
Chen, P. L., Liu, F., Cai, S., Lin, X., Li, A., Chen, Y., et al. (2005). Inactivation of CtIP leads to early embryonic lethality mediated by G1 restraint and to tumorigenesis by haploid insufficiency. Mol. Cell Biol. 25, 3535–3542. doi: 10.1128/MCB.25.9.3535-3542.2005
Cheng, T., Xu, M., Qin, B., Wu, J., Tu, Y., Kang, L., et al. (2019). lncRNA H19 contributes to oxidative damage repair in the early age-related cataract by regulating miR-29a/TDG axis. J. Cell Mol. Med. 23, 6131–6139. doi: 10.1111/jcmm.14489
David, C. J., Chen, M., Assanah, M., Canoll, P., and Manley, J. L. (2010). HnRNP proteins controlled by c-Myc deregulate pyruvate kinase mRNA splicing in cancer. Nature 463, 364–368. doi: 10.1038/nature08697
Falaleeva, M., Pages, A., Matuszek, Z., Hidmi, S., Agranat-Tamir, L., Korotkov, K., et al. (2016). Dual function of C/D box small nucleolar RNAs in rRNA modification and alternative pre-mRNA splicing. Proc. Natl. Acad. Sci. U.S.A. 113, E1625–1634. doi: 10.1073/pnas.1519292113
Gabory, A., Ripoche, M. A., Le Digarcher, A., Watrin, F., Ziyyat, A., Forne, T., et al. (2009). H19 acts as a trans regulator of the imprinted gene network controlling growth in mice. Development 136, 3413–3421. doi: 10.1242/dev.036061
Gazy, I., Zeevi, D. A., Renbaum, P., Zeligson, S., Eini, L., Bashari, D., et al. (2015). TODRA, a lncRNA at the RAD51 locus, is oppositely regulated to RAD51, and enhances RAD51-dependent DSB (double strand break) repair. PLoS ONE 10:e0134120. doi: 10.1371/journal.pone.0134120
Giovarelli, M., Bucci, G., Ramos, A., Bordo, D., Wilusz, C. J., Chen, C. Y., et al. (2014). H19 long noncoding RNA controls the mRNA decay promoting function of KSRP. Proc. Natl. Acad. Sci. U.S.A. 111, E5023–5028. doi: 10.1073/pnas.1415098111
Gong, J., Li, Y., Liu, C. J., Xiang, Y., Li, C., Ye, Y., et al. (2017). A pan-cancer analysis of the expression and clinical relevance of small nucleolar RNAs in human cancer. Cell Rep. 21, 1968–1981. doi: 10.1016/j.celrep.2017.10.070
Hsu, T. Y., Simon, L. M., Neill, N. J., Marcotte, R., Sayad, A., Bland, C. S., et al. (2015). The spliceosome is a therapeutic vulnerability in MYC-driven cancer. Nature 525, 384–388. doi: 10.1038/nature14985
Hu, Q., Zhang, B., Chen, R., Fu, C., Fu, A. J., et al. (2019). ZFHX3 is indispensable for ERbeta to inhibit cell proliferation via MYC downregulation in prostate cancer cells. Oncogenesis 8:28. doi: 10.1038/s41389-019-0138-y
Hu, W. L., Jin, L., Xu, A., Wang, Y. F., Thorne, R. F., Zhang, X. D., et al. (2018). GUARDIN is a p53-responsive long non-coding RNA that is essential for genomic stability. Nat. Cell Biol. 20, 492–502. doi: 10.1038/s41556-018-0066-7
Itzhaki, I., Maizels, L., Huber, I., Zwi-Dantsis, L., Caspi, O., Winterstern, A., et al. (2011). Modelling the long QT syndrome with induced pluripotent stem cells. Nature 471, 225–229. doi: 10.1038/nature09747
Keniry, A., Oxley, D., Monnier, P., Kyba, M., Dandolo, L., Smits, G., et al. (2012). The H19 lincRNA is a developmental reservoir of miR-675 that suppresses growth and Igf1r. Nat. Cell Biol. 14, 659–665. doi: 10.1038/ncb2521
Kim, H., Yoo, S., Zhou, R., Xu, A., Bernitz, J. M., Yuan, Y., et al. (2018). Oncogenic role of SFRP2 in p53-mutant osteosarcoma development via autocrine and paracrine mechanism. Proc. Natl. Acad. Sci. U.S.A. 115, E11128–E11137. doi: 10.1073/pnas.1814044115
Koh, C. M., Bezzi, M., Low, D. H., Ang, W. X., Teo, S. X., Gay, F. P., et al. (2015). MYC regulates the core pre-mRNA splicing machinery as an essential step in lymphomagenesis. Nature 523, 96–100. doi: 10.1038/nature14351
Koukourakis, M. I., Giatromanolaki, A., Sivridis, E., Simopoulos, C., Turley, H., Talks, K., et al. (2002). Hypoxia-inducible factor (HIF1A and HIF2A), angiogenesis, and chemoradiotherapy outcome of squamous cell head-and-neck cancer. Int. J. Radiat. Oncol. Biol. Phys. 53, 1192–1202. doi: 10.1016/s0360-3016(02)02848-1
Kuleshov, M. V., Jones, M. R., Rouillard, A. D., Fernandez, N. F., Duan, Q., Wang, Z., et al. (2016). Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 44, W90–97. doi: 10.1093/nar/gkw377
Lee, D. F., Kuo, H. P., Chen, C. T., Hsu, J. M., Chou, C. K., Wei, Y., et al. (2007). IKK beta suppression of TSC1 links inflammation and tumor angiogenesis via the mTOR pathway. Cell 130, 440–455. doi: 10.1016/j.cell.2007.05.058
Lee, D. F., Kuo, H. P., Liu, M., Chou, C. K., Xia, W., Du, Y., et al. (2009a). KEAP1 E3 ligase-mediated downregulation of NF-kappaB signaling by targeting IKKbeta. Mol. Cell 36, 131–140. doi: 10.1016/j.molcel.2009.07.025
Lee, D. F., Su, J., Ang, Y. S., Carvajal-Vergara, X., Mulero-Navarro, S., Pereira, C. F., et al. (2012a). Regulation of embryonic and induced pluripotency by aurora kinase-p53 signaling. Cell Stem Cell 11, 179–194. doi: 10.1016/j.stem.2012.05.020
Lee, D. F., Su, J., Sevilla, A., Gingold, J., Schaniel, C., and Lemischka, I. R. (2012b). Combining competition assays with genetic complementation strategies to dissect mouse embryonic stem cell self-renewal and pluripotency. Nat. Protoc. 7, 729–748. doi: 10.1038/nprot.2012.018
Lee, G., Papapetrou, E. P., Kim, H., Chambers, S. M., Tomishima, M. J., Fasano, C. A., et al. (2009b). Modelling pathogenesis and treatment of familial dysautonomia using patient-specific iPSCs. Nature 461, 402–406. doi: 10.1038/nature08320
Leppek, K., and Stoecklin, G. (2014). An optimized streptavidin-binding RNA aptamer for purification of ribonucleoprotein complexes identifies novel ARE-binding proteins. Nucleic Acids Res. 42:e13. doi: 10.1093/nar/gkt956
Li, A., Mallik, S., Luo, H., Jia, P., Lee, D. F., and Zhao, Z. (2020). H19, a long non-coding, RNA, mediates transcription factors and target genes through interference of MicroRNAs in pan-cancer. Mol. Ther. Nucleic Acids 21, 180–191. doi: 10.1016/j.omtn.2020.05.028
Lian, Q., Lye, E., Suan Yeo, K., Khia Way Tan, E., Salto-Tellez, M., Liu, T. M., et al. (2007). Derivation of clinically compliant MSCs from CD105+, CD24- differentiated human ESCs. Stem Cells 25, 425–436. doi: 10.1634/stemcells.2006-0420
Liang, W. C., Fu, W. M., Wong, C. W., Wang, Y., Wang, W. M., Hu, G. X., et al. (2015). The lncRNA H19 promotes epithelial to mesenchymal transition by functioning as miRNA sponges in colorectal cancer. Oncotarget 6, 22513–22525. doi: 10.18632/oncotarget.4154
Lin, Y. H., Jewell, B. E., Gingold, J., Lu, L., Zhao, R., Wang, L. L., et al. (2017b). Osteosarcoma: molecular pathogenesis and iPSC modeling. Trends Mol. Med. 23, 737–755. doi: 10.1016/j.molmed.2017.06.004
Liu, X. S., Haines, J. E., Mehanna, E. K., Genet, M. D., Ben-Sahra, I., Asara, J. M., et al. (2014). ZBTB7A acts as a tumor suppressor through the transcriptional repression of glycolysis. Genes Dev. 28, 1917–1928. doi: 10.1101/gad.245910.114
Luo, M., Li, Z., Wang, W., Zeng, Y., Liu, Z., and Qiu, J. (2013). Long non-coding RNA H19 increases bladder cancer metastasis by associating with EZH2 and inhibiting E-cadherin expression. Cancer Lett. 333, 213–221. doi: 10.1016/j.canlet.2013.01.033
Lustig-Yariv, O., Schulze, E., Komitowski, D., Erdmann, V., Schneider, T., de Groot, N., et al. (1997). The expression of the imprinted genes H19 and IGF-2 in choriocarcinoma cell lines. Is H19 a tumor suppressor gene? Oncogene 15, 169–177. doi: 10.1038/sj.onc.1201175
Ma, H., Yuan, L., Li, W., Xu, K., and Yang, L. (2018). The LncRNA H19/miR-193a-3p axis modifies the radio-resistance and chemotherapeutic tolerance of hepatocellular carcinoma cells by targeting PSEN1. J. Cell Biochem. 119, 8325–8335. doi: 10.1002/jcb.26883
Marcel, V., Ghayad, S. E., Belin, S., Therizols, G., Morel, A. P., Solano-Gonzalez, E., et al. (2013). p53 acts as a safeguard of translational control by regulating fibrillarin and rRNA methylation in cancer. Cancer Cell 24, 318–330. doi: 10.1016/j.ccr.2013.08.013
Mulero-Navarro, S., Sevilla, A., Roman, A. C., Lee, D. F., D'Souza, S. L., Pardo, S., et al. (2015). Myeloid dysregulation in a human induced pluripotent stem cell model of PTPN11-associated juvenile myelomonocytic leukemia. Cell Rep. 13, 504–515. doi: 10.1016/j.celrep.2015.09.019
Wan, G., Hu, X., Liu, Y., Han, C., Sood, A. K., Calin, G. A., et al. (2013). A novel non-coding RNA lncRNA-JADE connects DNA damage signalling to histone H4 acetylation. EMBO J 32, 2833–2847. doi: 10.1038/emboj.2013.221
Wu, L., Chang, L., Wang, H., Ma, W., Peng, Q., and Yuan, Y. (2018). Clinical significance of C/D box small nucleolar RNA U76 as an oncogene and a prognostic biomarker in hepatocellular carcinoma. Clin. Res. Hepatol. Gastr. 42, 82–91. doi: 10.1016/j.clinre.2017.04.018
Yagi, T., Ito, D., Okada, Y., Akamatsu, W., Nihei, Y., Yoshizaki, T., et al. (2011). Modeling familial Alzheimer's disease with induced pluripotent stem cells. Hum. Mol. Genet. 20, 4530–4539. doi: 10.1093/hmg/ddr394
Yoshimizu, T., Miroglio, A., Ripoche, M. A., Gabory, A., Vernucci, M., Riccio, A., et al. (2008). The H19 locus acts in vivo as a tumor suppressor. Proc. Natl. Acad. Sci. U.S.A. 105, 12417–12422. doi: 10.1073/pnas.0801540105
Zhang, Y., Xu, C., Gu, D., Wu, M., Yan, B., Xu, Z., et al. (2017). H/ACA box small nucleolar RNA 7A promotes the self-renewal of human umbilical cord mesenchymal stem cells. Stem Cells 35, 222–235. doi: 10.1002/stem.2490
Zheng, Z. G., Xu, H., Suo, S. S., Xu, X. L., Ni, M. W., Gu, L. H., et al. (2016). The essential role of H19 contributing to cisplatin resistance by regulating glutathione metabolism in high-grade serous ovarian cancer. Sci. Rep. 6:26093. doi: 10.1038/srep26093
Zhou, R., Xu, A., Gingold, J., Strong, L. C., Zhao, R., and Lee, D. F. (2017). Li-fraumeni syndrome disease model: a platform to develop precision cancer therapy targeting oncogenic p53. Trends Pharmacol. Sci. 38, 908–927. doi: 10.1016/j.tips.2017.07.004
Zhou, R., Xu, A., Tu, J., Liu, M., Gingold, J. A., Zhao, R., et al. (2018). Modeling osteosarcoma using li-fraumeni syndrome patient-derived induced pluripotent stem cells. J. Vis. Exp. 13:57664. doi: 10.3791/57664
Zhu, D., Kong, C. S. L., Gingold, J. A., Zhao, R., and Lee, D. F. (2018). Induced pluripotent stem cells and induced pluripotent cancer cells in cancer disease modeling. Adv. Exp. Med. Biol. 1119, 169–183. doi: 10.1007/5584_2018_257
Zhu, Q. N., Wang, G., Guo, Y., Peng, Y., Zhang, R., Deng, J. L., et al. (2017). LncRNA H19 is a major mediator of doxorubicin chemoresistance in breast cancer cells through a cullin4A-MDR1 pathway. Oncotarget 8, 91990–92003. doi: 10.18632/oncotarget.21121
Keywords: osteosarcoma, H19 lncRNA, iPSCs, Li-Fraumeni syndrome, snoRNA, p53
Citation: Xu A, Huang M-F, Zhu D, Gingold JA, Bazer DA, Chang B, Wang D, Lai C-C, Lemischka IR, Zhao R and Lee D-F (2021) LncRNA H19 Suppresses Osteosarcomagenesis by Regulating snoRNAs and DNA Repair Protein Complexes. Front. Genet. 11:611823. doi: 10.3389/fgene.2020.611823
Received: 29 September 2020; Accepted: 15 December 2020;
Published: 15 January 2021.
Edited by:Liang Chen, Wuhan University, China
Reviewed by:Matthias S. Leisegang, Goethe University Frankfurt, Germany
Dandan Li, Nanyang Normal University, China
Copyright © 2021 Xu, Huang, Zhu, Gingold, Bazer, Chang, Wang, Lai, Lemischka, Zhao and Lee. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
†These authors have contributed equally to this work