Next Article in Journal
Acknowledgment to Reviewers of Non-Coding RNA in 2020
Next Article in Special Issue
HiMoRNA: A Comprehensive Database of Human lncRNAs Involved in Genome-Wide Epigenetic Regulation
Previous Article in Journal
Evaluation of the Interplay between the ADAR Editome and Immunotherapy in Melanoma
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
ncRNA
Volume 7
Issue 1
10.3390/ncrna7010006
ncrna-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Hypothesis

Secondary Structural Model of MALAT1 Becomes Unstructured in Chronic Myeloid Leukemia and Undergoes Structural Rearrangement in Cervical Cancer

by
Matthew C. Wang
,
Phillip J. McCown
†,
Grace E. Schiefelbein
and
Jessica A. Brown
*
Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, IN 46556, USA
*
Author to whom correspondence should be addressed.
†
Present address: Department of Internal Medicine, Division of Nephrology, University of Michigan, Ann Arbor, MI 48109, USA.
Non-Coding RNA 2021, 7(1), 6; https://doi.org/10.3390/ncrna7010006
Submission received: 30 November 2020 / Revised: 11 January 2021 / Accepted: 11 January 2021 / Published: 13 January 2021
(This article belongs to the Special Issue Systematic Analysis of lncRNA Structures and Functions)
Browse Figures
Versions Notes

Abstract

:
Long noncoding RNAs (lncRNAs) influence cellular function through binding events that often depend on the lncRNA secondary structure. One such lncRNA, metastasis-associated lung adenocarcinoma transcript 1 (MALAT1), is upregulated in many cancer types and has a myriad of protein- and miRNA-binding sites. Recently, a secondary structural model of MALAT1 in noncancerous cells was proposed to form 194 hairpins and 13 pseudoknots. That study postulated that, in cancer cells, the MALAT1 structure likely varies, thereby influencing cancer progression. This work analyzes how that structural model is expected to change in K562 cells, which originated from a patient with chronic myeloid leukemia (CML), and in HeLa cells, which originated from a patient with cervical cancer. Dimethyl sulfate-sequencing (DMS-Seq) data from K562 cells and psoralen analysis of RNA interactions and structure (PARIS) data from HeLa cells were compared to the working structural model of MALAT1 in noncancerous cells to identify sites that likely undergo structural alterations. MALAT1 in K562 cells is predicted to become more unstructured, with almost 60% of examined hairpins in noncancerous cells losing at least half of their base pairings. Conversely, MALAT1 in HeLa cells is predicted to largely maintain its structure, undergoing 18 novel structural rearrangements. Moreover, 50 validated miRNA-binding sites are affected by putative secondary structural changes in both cancer types, such as miR-217 in K562 cells and miR-20a in HeLa cells. Structural changes unique to K562 cells and HeLa cells provide new mechanistic leads into how the structure of MALAT1 may mediate cancer in a cell-type specific manner.

1. Introduction

Long noncoding RNAs (lncRNAs; all abbreviations henceforth can be found in Supplemental Table S1, “Abbreviations” tab) are involved in a wide array of biological functions in humans, including imprinting, cell differentiation, and disease proliferation [1,2]. lncRNAs regulate gene expression through multiple mechanisms, including alternative splicing [3], binding transcription factors [4,5], and binding microRNAs (miRNAs) [6]. Specifically, lncRNAs can act as competing endogenous RNAs (ceRNAs) and sponge miRNAs, thus hindering those miRNAs from promoting degradation of their intended mRNA targets [7,8,9,10]. In some situations, miRNA binding can even promote lncRNA degradation [9]. The ability of lncRNAs to interact with other RNAs and proteins is largely mediated by secondary structure. For instance, lncRNAs can sponge miRNAs when the binding site is free of secondary structures (e.g., hairpins and pseudoknots) and/or proteins. In cancer and other human diseases, the proper functioning of lncRNAs is directly impacted by miRNA-binding events, where miRNA sponging can promote mRNA dysregulation and aberrant gene expression [11]. lncRNAs are often associated with oncogenic or tumor-suppressing effects and have high potential for use as cancer biomarkers [12]. While the structure dictates function of RNA and proteins, the distinct functional roles of structural elements in many lncRNAs are unclear.
One notable cancer-associated lncRNA is human metastasis-associated lung adenocarcinoma transcript 1 (MALAT1), whose longest isoform has over 8700 nucleotides [13] and its mature form terminates in a triple helix at its 3′ end [14]. Additionally, MALAT1 has hnRNPC- and hnRNPG-binding sites which are made available when N6-methyladenosine (m6A) modifications destabilize characterized hairpins [15,16,17]. Numerous roles for MALAT1 have been proposed, including regulation of pre-mRNA splicing, transcriptional regulation through binding of transcription factors, and acting as a ceRNA [18]. MALAT1 has previously been linked to SR protein phosphorylation and alternative splicing of pre-mRNAs [3]. Aberrant MALAT1 upregulation is considered oncogenic and occurs in breast, cervical, liver, and lung cancers, among others [18]. Moreover, MALAT1 has numerous known miRNA- and protein-binding partners, some of which mediate cancer. For example, the sponging of miR-217 by MALAT1 promotes dasatinib resistance in K562 cells by upregulating AGR2 [19]. Similarly, MALAT1 sponges the tumor suppressors miR-145 [20] and miR-202-3p [21], prompting invasiveness and epithelial-mesenchymal transitioning in HeLa cells [22]. The strong correlation between MALAT1 expression and cancer [18], as well as numerous known interactions with miRNAs and proteins that can bind to MALAT1, has made MALAT1 a promising biomarker and anticancer therapeutic target [23]. While increased expression of MALAT1 has been correlated with cancer or severe cancer phenotypes, how the roles of MALAT1 are influenced by its myriad structural elements is not yet clear.
Previously, a secondary structural model of human MALAT1 in noncancerous cells was proposed [15]. The model posits 194 helices and 13 pseudoknots covering 8425 nucleotides of MALAT1 and identified several unique and dynamic MALAT1 secondary structures, including a putative m6A switch that regulates miRNA binding [15]. According to this model, when there is a lack of methylation at A5044 in HeLa cells, a short hairpin is favored over a pseudoknot, thereby possibly increasing the accessibility of cancer-associated miR-101-3p, miR-217-5p, and miR-383-5p to their binding sites in MALAT1 [15,24,25,26]. Based on this structural switch and other cell-dependent structural predictions, we decided to examine structural data for MALAT1 in cancerous contexts. Herein, the working noncancer-derived secondary structural model of MALAT1 was compared to previously published RNA structural probing data in chronic myeloid leukemia (CML)-derived K562 cells [27] and in cervical cancer-derived HeLa cells [28,29] to identify putative differences between cancerous and noncancerous MALAT1 structural models [15]. DMS-Seq data from K562 cells [27] and PARIS data from HeLa cells [28,29] were used to determine how the structural model of MALAT1 changes in cancer. Both datasets were analyzed to generate a more complete picture of MALAT1, as the output for DMS-Seq provides structural information for only adenosine and cytidine and PARIS detects only psoralen-crosslinked duplexes [27,28,29]. From these comparative analyses, MALAT1 is predicted to become unstructured in K562 cells, with 59 of 101 analyzed hairpins losing at least half of the base pairs predicted in MALAT1 in a noncancerous context. Conversely, the MALAT1 model is largely maintained in HeLa cells, but 18 novel RNA-RNA interactions suggest structural rearrangements predominate. These structural changes would subsequently impact accessibility of 50 miRNA-binding sites. For instance, the binding site for miR-217, which is associated with tumor suppression [30], is subject to linearization in K562 cells, which would allow MALAT1 to sponge miR-217. Overall, this work highlights the potential, novel dynamic secondary structural changes in K562 and HeLa cells, whereby differential effects on miRNA-binding sites, protein-binding sites, RNA modifications, single nucleotide polymorphisms (SNPs), and cancer-associated mutations suggest mechanisms by which MALAT1 plays different roles in K562 and HeLa cells.

2. Results and Discussion

2.1. DMS-Seq Data Suggest Unfolding of MALAT1 Structure in K562 Cells

DMS-Seq involves chemically labeling RNA with DMS on unstructured adenosine and cytidine residues, which stops reverse transcriptase in a manner that can be detected by sequencing [27]. K562 DMS-Seq data [27] for MALAT1 (herein referred to as K562-MALAT1) were first analyzed to determine which adenosine and cytidine nucleotides are unstructured or structured (Supplemental Table S1). Of the 8425 nucleotides within the human MALAT1 transcript that were examined, 3951 nucleotides are either adenosine or cytidine. Of these, DMS-Seq data determined using MALAT1 isolated from K562 cells were available for 2554 adenosine and cytidine nucleotides, whereby 1835 DMS-Seq (71.8%) datapoints (i.e., number of DMS-Seq reads corresponding to a single adenosine or cytidine) were classified as unstructured and 719 DMS-Seq datapoints (28.2%) were classified as structured pursuant to the 250-read threshold (see Section 3 and Figure 1A). When K562-MALAT1 was compared to the working noncancerous MALAT1 model (Supplemental Table S1), 1504 datapoints (58.9%) agreed with the consensus model while 733 datapoints (28.7%) corresponded to loss of structure and 317 datapoints (12.4%) corresponded to gain of structure (Figure 1B). While the majority of K562-MALAT1 agreed with the MALAT1 consensus model, 41.1% of datapoints diverged from the noncancerous model, suggesting that wide-sweeping changes in MALAT1 secondary structure may occur in K562 cells.
K562-MALAT1 data that corresponded to the loss or gain of structure were mapped onto the noncancerous MALAT1 model to examine how the secondary structure of MALAT1 may vary in K562 cells (Figure 2 and Supplemental Figure S1). Of the 194 hairpins predicted in the noncancerous MALAT1 model, 101 hairpins from the noncancerous model are supported by DMS-Seq data and 59 of these hairpins (58.4%) lose at least half of their base pairs in K562 cells. Among these hairpins are: H44, losing 13 of 25 base pairs (52.0%); H45, losing 5 of 8 base pairs (62.5%); H98, losing 8 of 12 base pairs (66.7%); and H155, losing 4 of 6 base pairs (66.7%) (Figure 3, “K562-MALAT1”). These correspond to novel losses of eight base pairs in H44, four base pairs in H45, six base pairs in H98, and two base pairs in H155 (Figure 3, “Change”). The comparative analysis suggests widespread loss of structural features within the context of K562-MALAT1. While individual nucleotides occasionally appear to gain structure in K562-MALAT1, the sporadic occurrences do not suggest the development of any unambiguous novel secondary structures in MALAT1 in K562 cells (Figure 2). It is worth noting that refolding of the MALAT1 structure based on the K562 DMS-Seq data is expected to produce a novel secondary structure of K562-MALAT1. However, DMS-Seq is the only major RNA structural probing dataset available for K562 cells and it lacks data for about 1600 nts from a central region of MALAT1; therefore, a novel model cannot be constructed and the analysis herein is restricted to identifying regions of MALAT1 that potentially change in K562 cells. Overall, K562-MALAT1 results suggest many hairpins in the working MALAT1 model lose structure, and unstructured regions remain unstructured. This result is in agreement with prior work using DMS-Seq data, which found loss of structure in mRNAs [27]. Cumulatively, the K562-MALAT1 data indicate general loss of structure in MALAT1, thereby suggesting possible functional ramifications within the context of K562 cells.

2.2. Predicted Secondary Structural Changes in K562-MALAT1 Would Impact Multiple RNA- and Protein-Binding Sites

Loss of secondary structure in K562-MALAT1 may signal that certain RNA- and protein-binding sites are now available in MALAT1, especially for single-stranded RNA-binding proteins. As such, this possibility was examined to identify the aberrant binding events in K562-MALAT1 that are different from binding and interaction events for MALAT1 in noncancerous conditions. miRNAs, ncRNAs, proteins, RNA modifications, SNPs, and cancer-associated mutations were re-aligned to MALAT1 to identify structure-function relationships that provide a starting point to examine their possible roles in CML (Supplemental Table S1).
miRNAs are currently known to play a pivotal role in the development and progression of CML [31]. Of the 98 validated miRNA-binding sites in MALAT1, 28 sites occur in hairpins predicted to lose structure in K562-MALAT1, thereby increasing accessibility of binding site and potential for sponging (Supplemental Table S1). Examples include miR-320, which overlaps with H101; miR-217, which overlaps with H160; and miR-140-5p, which overlaps with H168 (Figure 4). miR-320 is considered a tumor suppressor in K562 cells, but K562 cells often bypass its action by transporting miR-320 to exosomes via hnRNPA1 [32]. Sponging of miR-320 by MALAT1 in K562 cells could also dampen the tumor-suppressive effects of miR-320, as is the case with lncRNA SNHG12 sponging miR-320 in gastric cancer [33]. miR-217 reportedly targets the mRNA of oncogenic protein AGR2 in K562 cells [19]. As decreases in unbound miR-217 accompany AGR2 upregulation and subsequent dasatinib resistance in K562 cells [19], sponging of miR-217 by MALAT1 may have similar effects in K562 cells. miR-140-5p has been linked to CML cell apoptosis via targeting of the SIX1 mRNA transcript [34]; therefore, possible sponging of miR-140-5p by MALAT1 in K562 cells may promote cell survival. These examples highlight how the novel availability of miRNA-binding sites in MALAT1 may aid in K562 cell progression via multiple miRNA-mediated mechanisms.
Eleven hairpins expected to lose structure in K562-MALAT1 (H49, H77, H79, H80, H147, H148, H155, H156, H164, H165, H168) overlap with eight of the ten U1 snRNA-binding sites (nts 1825–1925, 3015–3067, 3152–3185, 5924–6023, 6127–6277, 6850–6884, 6985–7045, and 7138–7206) (Supplemental Table S1). U1 snRNA is known to be mutated in multiple cancer types to promote aberrant gene splicing patterns [35]. Although the roles of U1-MALAT1 interactions have not been elucidated, it is conceivable that binding of U1 snRNA to MALAT1 may also contribute toward alternate, oncogenic splicing patterns that promote CML. Furthermore, one U1 snRNA-binding site (nts 3152–3185) overlaps with a HuR/ELAV1-binding site (nts 3158–3163) that may become accessible upon loss of H79 in K562-MALAT1. Although the effects of HuR-MALAT1 binding on CML have not been investigated, the HuR-MALAT1 RNP complex has been shown to stop breast cancer cells from undergoing epithelial-mesenchymal transitioning by decreasing the levels of CD133 [36]. Thus, increased HuR-MALAT1 binding in K562 cells is expected to hinder cancer progression, suggesting the existence of alternate pathways by which HuR-MALAT1 binding affects K562 cells. HuR typically binds to mRNAs in cancer in order to promote cancerous functions, such as metastasis and apoptosis resistance [37]. Thus, competition between HuR and U1 snRNA for a binding site around nt 3160 may point to a carefully controlled cancer-promoting mechanism mediated by MALAT1. In general, characterized protein-binding sites on MALAT1 are expected to become more available as a result of widespread structural loss and these changes in protein-MALAT1 binding predicted by K562 DMS-Seq data hint at novel pathways to explore further.
Besides RNA- and protein-binding sites, RNA structure can also be modulated by RNA modifications, SNPs, and cancer-associated mutations. RNA modifications on MALAT1 in K562 cells are undetermined, so modification sites from other cell lines were used. Of all the 82 m6A modifications mapped to MALAT1 at single-nucleotide resolution in HEK293, HEK293T, and HeLa cells [24,38,39,40,41], 61 m6A modifications either overlap with hairpins predicted to lose structure or, if not overlapping with hairpins, correspond to adenosines predicted to be unstructured (Supplemental Table S1). The METTL3/14 complex, which is responsible for about 80% of m6A marks in human mRNAs and ncRNAs, has no strong preference for ssRNA or dsRNA [42], suggesting secondary structural changes are insufficient to predict changes in m6A levels caused by METTL3/14. It is worth noting that the METTL3/14 complex is considered tumor suppressive and is downregulated in cancers like endometrial cancer, whereas m6A erasers like ALKBH5 and FTO are oncogenic and often upregulated in cancers like acute myeloid leukemia (AML) and breast cancer [43]. Correspondingly, ALKBH5 is associated with MALAT1 upregulation [44] and FTO regulates MALAT1 levels via demethylation [45]. ALKBH5 does not discriminate between ssRNA and dsRNA [46] and FTO targets ssRNA [47]; therefore, m6A marks in MALAT1 are potential substrates for both m6A erasers. Also notable, m6A2515 enables the binding of hnRNPG to MALAT1 at H63 [16], which is lost in K562-MALAT1 (Figure 2). As hnRNPG has stronger binding affinity for ssRNA, particularly in A-rich regions, [48,49] and has generally been associated with tumor suppressive effects [50,51,52], increased MALAT1-hnRNPG binding may decrease the tumor suppressive activity of hnRNPG and may promote K562 cell progression. While the functional effects of any aberrant methylation patterns are difficult to predict in CML, m6A modifications and their roles in RNA regulation and function, particularly with regard to mRNAs where such modifications are the most common, have been explored in attempts to develop novel cancer biomarkers and treatments [53]. Therefore, understanding how structural alterations in MALAT1 modulate m6A modification sites in K562 cells is of particular interest.
Seventeen SNPs have been identified in MALAT1 [54]. rs664589 (C4117G), rs115795653 (A6415G), and rs60151940 (C7151W) are three SNPs that correspond to nucleotides that are predicted to lose structure in K562-MALAT1 (Supplemental Table S1). SNPs in structured RNAs are generally believed to alter the local secondary structure [55], although the severity of alterations can vary and is difficult to predict [56]. Interestingly, SNP rs664589 has been characterized as aiding colorectal cancer progression by inhibiting MALAT1-miR-194-5p binding [57]. Besides SNPs, 655 somatic cancer-associated mutations have been identified in MALAT1 [58]. Fifty-nine mutations (9.0%) correspond to nucleotides predicted to lose structure in K562-MALAT1 (Supplemental Table S1). Such mutations could further weaken the hairpin structures or reduce miRNA binding, particularly if the mutation were to disrupt base pairing in the seed region of the miRNA-binding site. Within H71, the A2875U mutation alters seed-region base pairing for miR-92a-3p, miR-363-3p, and miR-25-3p (Figure 5 and Supplemental Table S1). Although miR-363-3p is associated with tumor suppression in other cancer types [59], miR-92a-3p and miR-25-3p promote progression in cancers like liposarcomas [60]. Moreover, miR-92a-3p was previously found to aid CML by downregulating C/EBPα and subsequently causing cachexia, i.e., severe weight and muscle loss associated with cancer [61]. As predicted previously, this proposed role of H71 in regulating MALAT1-miRNA interactions illustrates how H71 can modulate the different outputs depending on the cellular context [15]. In total, 16 miRNA-binding sites have seed regions within hairpins predicted to become unstructured in K562-MALAT1 (Supplemental Table S1, “SeqMarkup” tab), and may experience reduced binding affinity due to somatic cancer-associated mutations. Additionally, one METTL3/14-binding site, two HuR/ELAV-binding sites, and ten U1-binding sites face similar conditions because of mutations, pointing to complex regulatory pathways that may depend on the K562-MALAT1 structure. Together, the effects of secondary structural loss on miRNA-, U1 snRNA and HuR-, m6A-, SNP-, and cancer-associated mutation-related effects in K562 cells represent potential avenues for further characterization of MALAT1 activity in cancer.

2.3. PARIS Data Suggest Maintenance of Overall Structure of MALAT1 with Rearrangements of Select Long-Range Interactions

PARIS involves sequencing of RNA fragments that were once photocrosslinked-duplexes isolated from psoralen-treated cells [28]. PARIS data for MALAT1 in HeLa cells (herein referred to as HeLa-MALAT1) were compared and mapped to the noncancerous consensus model [15] to see how the MALAT1 secondary structural model is expected to change within the context of HeLa cells (Figure 6, Supplemental Figure S2 and Supplemental Tables S1 and S2) [28,29]. Eighty unique PARIS interactions (Supplemental Table S2) were aligned to the MALAT1 model. Of these, 18 PARIS interactions (22.5%) diverged from hairpins described in the consensus model while 62 local interactions (77.5%) agreed with the hairpins in the consensus model, suggesting that the secondary structure of MALAT1, as it is hypothesized to exist, may be largely maintained in HeLa cells (Figure 6). The 18 PARIS interactions that diverge from the model typically suggest that secondary structural elements undergo structural rearrangement. Of the 18 divergent PARIS interactions, six interactions (33.3%) were denoted as short-range or local interactions and 12 interactions (66.7%) were denoted as long-range interactions, most of which are separated by at least 80 nucleotides in their primary structure in the working noncancerous MALAT1 model (Supplemental Table S2) [15]. The 12 divergent long-range PARIS interactions typically signal the structural rearrangement of multiple hairpins whereas the six short-range PARIS interactions signal the formation of novel structures. Up to 42 hairpins out of 161 hairpins (26.1%) are expected to undergo rearrangement and 119 hairpins (73.9%), which are conserved among mammals [15], appear maintained in the HeLa-MALAT1 model (Supplemental Tables S1 and S2). Thus, with regard to structural alterations of hairpins, rearranging long-range interactions is preferred over novel short-range interactions.
Five of the aforementioned local PARIS interactions occur in predominantly unstructured regions of the working MALAT1 consensus model (Figure 6, dark red lines). Four of these interactions fall between nts 1897 and 1941 (i.e., between H49 and H50) and the fifth interaction falls between nts 7458 and 7461, preceding H174 (Figure 6). The noncancerous HEK293T PARIS data did not highlight any such structures at the corresponding locations (Supplemental Table S2). Together, these five interactions suggest distinct instances of dynamic, novel structures. In contrast, most of the 12 long-range PARIS interactions indicate distinct instances of structural rearrangement (Figure 6, purple lines). Curiously, five of the long-range PARIS interactions and one divergent short-range PARIS interaction start within 561 nucleotides of one another, spanning nts ~4950 to ~5600 (Figure 6 and Supplemental Table S2). This region is largely conserved among mammals as well as some vertebrates [15]. Thus, a core of MALAT1 undergoes structural rearrangement in HeLa cells. The 12 long-range PARIS interactions suggest rearrangement of 39 hairpins, such as H126, H134, H136, and H178. Long-range interactions suggest rearrangement of H105 (coordinates 6446,6564), which notably forms a 56-way junction, and H170 (coordinates 7631,8196), which notably forms a 20-way junction (Figure 6, 56WJ and 20WJ). The PARIS data suggest these long-range interactions are lost in favor of structural rearrangement in HeLa cells, as opposed to general structural loss in K562-MALAT1.
In addition to hairpins, the hypothetical consensus model predicts 13 pseudoknots in noncancerous cells [15]. As previously noted, m6A5044 is absent in HeLa cells [15,24,25,26]. This loss of methylation may result in the loss of PK7 as there is a lack of PARIS data for PK7 (coordinates 5038, 6642) in HeLa cells, as previously reported [15]. Instead, PARIS reads (coordinates 5038,5145) suggest formation of a local hairpin [15], as indicated by the sixth divergent short-range PARIS interaction (Figure 6). Additionally, long-range PARIS interactions suggest structural rearrangement of PK3 and PK9 while the structural rearrangement of PK10 and PK11 is supported by short-range interactions. PARIS data do not predict disruption of any other pseudoknots. Unlike hairpins, pseudoknots typically span long ranges of MALAT1 in the working noncancerous model [15]. As a result, loss of many pseudoknots would indicate widespread structural changes in MALAT1. While the loss of these pseudoknots signals some propensity for long-range structural changes, the maintenance of eight pseudoknots reaffirms the trend of structural maintenance within HeLa-MALAT1. Overall, most local secondary structural features are retained in HeLa-MALAT1, with rearrangement of select long-range secondary structures and formation of a small number of novel, local structures.

2.4. Predicted Structural Changes in HeLa-MALAT1 Would Impact RNA-Binding Sites and Modifications

Structural rearrangements or novel structures detected in HeLa-MALAT1 means that the structures of RNA- and protein-binding sites underwent changes that may potentially alter their function. MALAT1 has 98 experimentally verified miRNA-binding sites (Supplemental Table S1). Duplex formation in HeLa-MALAT1 suggests disruption of seed region-binding sites for 25 of these validated miRNAs (Supplemental Table S2), which means these binding sites would be less accessible in HeLa cells. The four local PARIS interactions spanning nts 1897–1941 potentially decrease binding site availability for miR-145-5p [62]. miR-145 has been shown to inhibit HeLa cell proliferation by targeting the FSCN1 mRNA transcript [63]. miR-145 is reportedly a tumor suppressor in HeLa cells via the regulation of several proteins, including CDKs and Cyclin D1 [64]. Although miR-145 is downregulated in HeLa cells [65], its function relative to expected changes in the HeLa-MALAT1 structure raise questions regarding the full role of miR-145 in HeLa cells. The remaining PARIS interactions indicate structural rearrangement of binding sites for 24 other supported miRNAs, including miR-200b-3p, miR-20a-5p, and miR-106b-5p [62]. miR-200b is upregulated in cervical cancer and aids cervical cancer metastasis by downregulating FOXG1 [66]. Likewise, miR-20a is upregulated in HeLa cells and leads to the upregulation of the oncogenic protein TNKS2 in HeLa cells [67]. miR-106b is also upregulated in HeLa cells [68] and inhibits HeLa cell proliferation by downregulating PTEN via sponging of miR-106b by the lncRNA PTENP1 [69]. The PARIS data suggest these three latter miRNAs will not be sponged by HeLa-MALAT1, thus possibly aiding HeLa cell growth and survival. Collectively, these studies suggest the presence of complex miRNA-lncRNA-mRNA networks that may be disrupted by changes to MALAT1 secondary structure in HeLa cells. Additional work is required to elucidate the full pathways governed by such miRNAs and to fully understand how structural changes in MALAT1 affect miRNA function in HeLa cells.
Besides miRNAs, MALAT1 has been described as forming intermolecular RNA–RNA interactions with rRNA and U1 snRNA [70,71]. The structural status of one of the five rRNA-binding sites (C2700) and two of the ten U1 snRNA-binding sites (nts 1825–1925 and 6985–7045) is changed in HeLa-MALAT1 (Supplemental Table S1). Because few sites are affected, little to no significant alteration to MALAT1-mediated U1 snRNA and rRNA function is expected in HeLa cells. Additionally, protein-binding sites on MALAT1 are expected to become less available as a result of structural rearrangement throughout HeLa-MALAT1. Three METTL3/14-binding sites (nts 2412–2416, 5042–5046 and 8179–8184) and one HuR/ELAV1-binding site (nts 3248–3258) are hypothesized to undergo structural rearrangement, as indicated by the PARIS data. As previously described, the lack of affinity for ssRNA or dsRNA makes analysis of novel METTL3/14 function with regard to MALAT1 difficult [42]. However, the METTL3/14 complex shows some increased affinity for single-stranded nucleic acids [72], so aberrant m6A levels are possible under such circumstances as a result of MALAT1 rearrangement. Unlike the K562 cells, loss of HuR-MALAT1 binding is expected in HeLa cells as structural rearrangement will make the HuR-binding site less available. This loss mirrors the aforementioned functions of HuR-MALAT1 binding in breast cancer [36], suggesting HuR-MALAT1 binding may be decreased in HeLa cells in order to target CD133 expression and subsequently promote cancer progression. Because only one HuR-binding site is expected to undergo structural rearrangement, the repercussions on HuR function may be muted. Although probing of this particular pathway is needed to confirm such a hypothesis, a possible role of HuR-MALAT1 binding is more compelling in HeLa cells than in K562 cells.
Interestingly, several RNA modification sites identified in MALAT1 isolated from HeLa cells occur in structurally rearranged regions: 19 m6A sites, five m5C sites (C4834, C5518, C5520, C5538, and C5539), and one Am (2′-O-methyladenosine) modification site (A1909) (Supplemental Tables S1 and S2). m5C modifications have been found to regulate chromatin-related roles in other lncRNAs, such as HOTAIR and Xist, for this modification often occurs specifically in regions of the lncRNA that interact with chromatin-associated protein complexes [73]. The five aforementioned m5C sites were specifically identified in HeLa cells (see Supplemental Table S1). All five m5C sites in MALAT1 are clustered within 705 nucleotides of each other, with four of them clustered within 21 nucleotides (Supplemental Table S1). Thus, because MALAT1 binds active chromatin, a novel structure in HeLa-MALAT1 may promote a distinct and cancer-specific chromatin-associated complex via m5C [74]. Moreover, the existence of modified nucleotides in MALAT1 and the diversity of RNA modifications, along with advances in modification detection, may result in the discovery of novel MALAT1 modifications that can be implemented as biomarkers [75]. Thus, integrating PARIS and RNA modification data yielded insights into how RNA modifications, particularly m5C, may influence MALAT1 function in HeLa cells.
MALAT1 has 17 SNPs [54]. The HeLa-MALAT1 data suggest structural rearrangements for three SNP sites: rs11540782 (U1876C), rs1056816 (A4872K), and rs79910129 (G3247W) (Supplemental Table S1). As previously stated, the exact effects of a given SNP on secondary structure can vary but often result in the disruption of duplexes and loss of secondary structure [56]. As such, based on the PARIS data, no major cellular changes are expected in HeLa cells related to the MALAT1 SNPs. Of the 655 somatic cancer-associated mutations that have been identified in MALAT1 [58], 102 mutations (15.6%) occur in regions predicted to undergo structural rearrangement in HeLa-MALAT1. Mutations within PARIS interactions are liable to destabilize the corresponding RNA duplexes but are also likely to disrupt the binding sites, thus decreasing binding of molecules like miRNAs. Eight miRNA seed-region binding sites in HeLa-specific PARIS interactions are altered by mutations, as are one METTL3/14-binding site and one HuR/ELAV1-binding site. A U5520 insertion alters the seed-region binding sites of three miRNAs within a long-range PARIS interaction (coordinates 5503, 5708): miR-17-5p, miR-20ab-5p, and miR-106b-5p. As discussed previously, free miR-20a is expected to aid HeLa cells via TNKS2 expression [67]. Both miR-17-5p and miR-106b-5p are described as oncogenic in cervical cancer [76,77]. miR-17-5p targets TGFBR2 and stimulates proliferation, and miR-106b-5p promotes PTEN downregulation to achieve similar effects [69]. Hence, there is the potential for somatic cancer-associated mutations to regulate MALAT1 function through structural changes in HeLa cells.

3. Materials and Methods

3.1. Dataset Acquisition

All sequence datasets used in this analysis were accessed using the Gene Expression Omnibus. DMS-Seq data from K562 cells were extracted from data file GSM1297493-GSM1297494 (GSE45803) [27]. PARIS data from HeLa cells were extracted from data file GSM1917754 (GSE74353), representing high RNase data [28,29]. K562 DMS-Seq and HeLa PARIS data were mapped to human MALAT1 relative to hg38 using the open source platform Galaxy [78] and the UCSC Genome Browser [79]. VARNA was used to visualize the MALAT1 secondary structure [80]. Human MALAT1 nucleotide positions correspond to accession NR_002819.2 and ENST00000534336.1.
The data for miRNA-binding sites, U1 snRNA-binding sites, rRNA-binding sites, protein-binding sites, RNA modifications, SNPs, and somatic mutations associated with cancer were previously curated by McCown et al. [15] (Supplemental Table S1) and updated as follows: miRNA data were downloaded from ENCORI (Sun Yat-sen University, Guangzhou, China) on 19 August 2020 [81]. Only experimentally verified miRNAs were considered in our analysis. Additional validated U1 snRNA sites were reported by Cai et al. [82]. Additional m6A modification data were downloaded from m6AVAR (Sun Yat-sen University, Guangzhou, China) and are accurate as of 1 November 2020 [26]. MALAT1 somatic mutations associated with cancer were accessed from the National Cancer Institute Genomic Data Commons (Bethesda, MD, USA) and were accurate as of 1 April 2020 [58].

3.2. Dataset Analysis and Comparison to MALAT1 Secondary Structural Model

K562 DMS-Seq data pertaining to MALAT1 were extracted and examined with respect to the previously determined hypothetical MALAT1 secondary structural model in noncancerous human foreskin fibroblasts [15]. As the ensuing analysis sought to investigate the hypothetical differences in structure as opposed to establish a new model or identify concrete structural changes, no data were excluded on the basis of ambiguity with respect to the noncancerous model [15]. Thresholds for determining whether DMS-Seq data indicated structure were established using several verified hairpins: H190 (triple helix-containing hairpin) [15] and H191- H194 (mascRNA) (Figure 2) [14,83]. Because MALAT1 expression is higher in cancer cells, including leukemia, than in noncancerous tissue [84], the threshold value of structure calls for DMS-Seq data was set at 250 reads rather than 20 [15]. This threshold of 250 reads was established in a manner similar to the threshold of 20 determined for the consensus MALAT1 structure [15]. Briefly, DMS-Seq labels only adenosine and cytidine in a statistically significant manner; therefore, only these residues were considered in our structural analysis. Adenosine and cytidine residues having no more than 250 DMS-Seq reads were classified as structured. Adenosine and cytidine residues with more than 250 DMS-Seq reads were considered unstructured. We chose 250 as the threshold because this value was approximately midway between nucleotides that were known to be unstructured (C8291 and A8292 at 314 and 322 reads, respectively) in H190 and nucleotides known to be structured (e.g., m1A8398 at 217 reads) in H190-H194 (Supplemental Table S1) [14,83,85].
Like the noncancerous fibroblast DMS-Seq data, K562 DMS-Seq data did not contain nts 1–1280 of MALAT1 because the major isoform of MALAT1 is nts 1281–8425 [83] nor nts 4258–5838 (see Supplemental Table S1). Of the 194 hairpins predicted by McCown et al., 124 hairpins were predicted between nts 1281–4258 and nts 5838–8425 [15]. Of these, 23 hairpins were predicted using PARIS data and were thus excluded from the analysis for this cell line. Hence, 101 hairpins were analyzed for structural alterations in K562 cells. Structural alterations within hairpins were determined by comparing the K562 DMS-Seq data to the noncancerous secondary structural model of human MALAT1 excluding PARIS-derived hairpins, as well as to human foreskin fibroblasts [27,29]. As all 13 pseudoknots were established using PARIS data, structural changes in pseudoknots were not considered in the K562 MALAT1 model.
To verify the hypothesis that the unstructured adenosines and cytidines have higher DMS reactivity values than structured adenosines and cytidines (Supplemental Table S1, see U Test of Noncancerous Data tab and U Test of K562-MALAT1 Data tab), we conducted a Mann–Whitney U test on the DMS-Seq data for the structured and unstructured adenosines and cytidines in noncancerous MALAT1 [15] and K562-MALAT1. The DMS-Seq values from all adenosines and cytidines classified as structured or unstructured were tabulated. For the noncancerous MALAT1 model, DMS-Seq data derived from fibroblasts [27] were sorted into structured and nonstructured categories based on their presence in structured or unstructured regions of the consensus noncancerous MALAT1 model. These DMS-Seq values were then ordered from smallest to largest, ranked, and subjected to a U test. This U test produced a statistically significant result, with a p value of ~0. For the K562-MALAT1 U test, structured nucleotides were all adenosine/cytidine nucleotides in stems of hairpins that maintained more than 50% of the base pairs (Supplemental Table S1, see Hairpin Coordinates tab). PARIS-derived hairpins were included. Ranks for structured and unstructured DMS-Seq values were determined and subjected to a U test. The worse performing U, corresponding to the unstructured portion of the K562-MALAT1 model, was subjected to a z statistic calculation and was determined to be significant with a p value of ~0, demonstrating that the mean values of structured DMS-Seq adenosines and cytidines do not differ randomly from the mean values of the unstructured DMS-Seq adenosines and cytidines and, as stated, the null hypothesis can be rejected. Although these U tests suggest that the unstructured DMS-Seq values are larger than the structured DMS-Seq values, please bear in mind that these DMS-Seq values are not truly independent observations due to nearest-neighbor effects of nucleotides within helices. Two additional U tests were also conducted to compare the DMS-Seq values for adenosine/cytidine nucleotides classified as unstructured in cancerous K562 and noncancerous fibroblasts cells and structured in cancerous and noncancerous cells. The DMS-Seq values for all structured and all unstructured nucleotides in both cell types were sorted into respective columns, ranked, and subjected to U tests. Both U tests produced statistically significant results, with p values of ~0. Please note that these statistical tests may be impacted by the different thresholds applied to the noncancerous and cancerous DMS-Seq data, thereby ensuring different means.
HeLa PARIS data were examined with respect to a previously determined hypothetical MALAT1 secondary structural model in noncancerous HEK293T cells [15]. Because the HEK293T PARIS samples were prepared under high RNase conditions, we likewise examined the high RNase dataset for HeLa cells [28,29]. First, we removed low quality RNA-seq reads, which were lower than Q30 on the Illumina quality score metrics for each read. Reads that formed duplexes outside the MALAT1 coordinates were excluded. Next, PCR duplicates were removed from our analysis by barcode matching in Linux. Finally, PARIS reads that overlapped by 20 or fewer nts were compressed into one region of at least 10 nts and defined as a double-stranded region for any region having at least three PARIS reads. Importantly, PARIS data from HEK293T and HeLa cells were not available for nts 1–1280 of MALAT1 because the major isoform of MALAT1 is nts 1281–8425 [83]; therefore, nts 1–1280 of MALAT1 were not examined for structural alterations in HeLa cells. Structural alterations within hairpins and pseudoknots were determined by comparing the HeLa PARIS data to the noncancerous secondary structural model of human MALAT1, as well as to HEK293T PARIS data [28,29]. When aligning HeLa PARIS data to MALAT1, the PARIS data were assigned specific ranges of nts (see Supplemental Table S2), occasionally creating an apparent overlap among the PARIS interactions (see Figure 6). Because these overlaps are likely artificial, they were not considered to be indicating dynamic structures but to be indicating PARIS interactions in close proximity.

3.3. Data and Software Availability

All data and software are freely available at their designated repositories as indicated above. To the best of our knowledge, there are no restrictions or embargoes in place on any of these data.

4. Conclusions

Analyzing the K562-MALAT1 and HeLa-MALAT1 models provides insights into the roles and mechanisms of MALAT1 in two different cancer cell lines: K562 and HeLa cells. Based on the putative structural changes with respect to a working model in noncancerous cells, we hypothesize that MALAT1 possesses different secondary structures in both K562 and HeLa cells. However, the nature of these structural changes is distinct to each cancer type and does not appear to have obvious overlaps. In K562 cells, we predict that at least 30% of all MALAT1 hairpins will lose at least 50% of their base pairs (Figure 1, Figure 2 and Figure 3). This putative loss of structure may increase the propensity for miRNA sponging by K562-MALAT1, causing profound effects on cancer cell function. In contrast, PARIS data in HeLa cells suggest that long-range interactions occur in HeLa-MALAT1 that are not seen in noncancerous cells, likely due to alternative structuring of some hairpins and select pseudoknots in MALAT1, although most local secondary structures are preserved. The novel long-range interactions that we predict in HeLa-MALAT1 are expected to decrease instances of miRNA sponging and perhaps alter the functional readout of m5C methylations on MALAT1, having widespread consequences in K562 cells. Toggling of structures near m5C marks may influence the chromatin restructuring and gene expression in HeLa cells. Whereas K562-MALAT1 is predicted to have increased miRNA-sponging capabilities in K562 cells, alterations to m5C-containing structures in HeLa-MALAT1 may point to novel MALAT1-chromatin interactions, hinting at diverse biological processes being coordinated by MALAT1. The full extent to which the proposed structural alterations affect cancer cell development and progression await experimental validation, from establishing a secondary structural model of MALAT1 to confirming isolated structure-function relationships. However, we have identified the possible structural differences between MALAT1 in different cellular contexts that could exacerbate K562 and HeLa cells, leading to mechanistic insights regarding the complex cancer-specific functions of MALAT1.

Supplementary Materials

The following are available online at https://www.mdpi.com/2311-553X/7/1/6/s1. Supplemental Figure S1: Poster-size image (approximately 19 in. by 31 in.) of Figure 2 showing K562 DMS-Seq structural changes on the MALAT1 secondary structural model. Supplemental Figure S2: Poster-size image (approximately 19 in. by 31 in.) of Figure 6 showing HeLa PARIS data on the MALAT1 secondary structural model. Supplemental Table S1: Sequence markup and analysis of K562-MALAT1 structural changes on various structure-dependent functions. Supplemental Table S2: Analysis of HeLa-MALAT1 structural changes on various structure-dependent functions.

Author Contributions

Conceptualization, M.C.W., P.J.M., and J.A.B.; methodology, M.C.W. and P.J.M.; software, P.J.M.; validation, M.C.W. and P.J.M.; formal analysis, M.C.W., P.J.M., and J.A.B.; investigation, M.C.W., P.J.M., and G.E.S.; resources, P.J.M.; data curation, M.C.W., P.J.M., and G.E.S.; writing—original draft preparation, M.C.W., P.J.M., and J.A.B.; writing—review and editing, M.C.W., P.J.M., G.E.S., and J.A.B.; visualization, M.C.W. and P.J.M.; supervision, P.J.M. and J.A.B.; project administration, M.C.W. and J.A.B.; funding acquisition, J.A.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by startup funds provided by the University of Notre Dame and the Clare Boothe Luce Program of the Henry Luce Foundation to J.A.B. and by the Sean O’Keefe Cancer Fund provided by the University of Notre Dame College of Science Summer Undergraduate Research Fellowship Program to M.C.W.

Data Availability Statement

The data that support the findings of this study are openly available at the Gene Expression Omnibus at http://www.ncbi.nlm.nih.gov/geo using GSE or GSM numbers as indicated within the article and/or its Supplementary Materials.

Acknowledgments

We thank Charlotte Kunkler and Mika Schievelbein for critical review of the manuscript and for all Brown laboratory members for thoughtful discussions. We acknowledge that some of the publicly available data used in this study were obtained using HEK cell lines; the authors do not agree with nor encourage research performed on HEK cell lines.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Quinn, J.J.; Chang, H.Y. Unique features of long non-coding RNA biogenesis and function. Nat. Rev. Genet. 2016, 17, 47–62. [Google Scholar] [CrossRef] [PubMed]
  2. Chillón, I.; Marcia, M. The molecular structure of long non-coding RNAs: Emerging patterns and functional implications. Crit. Rev. Biochem. Mol. Biol. 2020, 55, 662–690. [Google Scholar] [CrossRef] [PubMed]
  3. Tripathi, V.; Ellis, J.D.; Shen, Z.; Song, D.Y.; Pan, Q.; Watt, A.T.; Freier, S.M.; Bennett, C.F.; Sharma, A.; Bubulya, P.A.; et al. The nuclear-retained noncoding RNA MALAT1 regulates alternative splicing by modulating SR splicing factor phosphorylation. Mol. Cell 2010, 39, 925–938. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Jeon, Y.; Lee, J.T. YY1 Tethers Xist RNA to the inactive X nucleation center. Cell 2011, 146, 119–133. [Google Scholar] [CrossRef] [Green Version]
  5. Vance, K.W.; Ponting, C. Transcriptional regulatory functions of nuclear long noncoding RNAs. Trends Genet. 2014, 30, 348–355. [Google Scholar] [CrossRef] [Green Version]
  6. Ballantyne, M.D.; McDonald, R.A.; Baker, A.H. LncRNA/MicroRNA interactions in the vasculature. Clin. Pharmacol. Ther. 2016, 99, 494–501. [Google Scholar] [CrossRef] [Green Version]
  7. Salmena, L.; Poliseno, L.; Tay, Y.; Kats, L.; Pandolfi, P.P. A ceRNA hypothesis: The Rosetta Stone of a hidden RNA language? Cell 2011, 146, 353–358. [Google Scholar] [CrossRef] [Green Version]
  8. Paraskevopoulou, M.D.; Hatzigeorgiou, A.G. Analyzing miRNA-lncRNA interactions. In Long Non-Coding RNAs: Methods and Protocols; Feng, Y., Zhang, L., Eds.; Humana Press: New York, NY, USA, 2016; Volume 1402, pp. 271–286. [Google Scholar]
  9. Yoon, J.H.; Abdelmohsen, K.; Gorospe, M. Functional interactions among microRNAs and long noncoding RNAs. Semin. Cell Dev. Biol. 2014, 34, 9–14. [Google Scholar] [CrossRef] [Green Version]
  10. Ergun, S.; Oztuzcu, S. Oncocers: ceRNA-mediated cross-talk by sponging miRNAs in oncogenic pathways. Tumor Biol. 2015, 36, 3129–3136. [Google Scholar] [CrossRef]
  11. Bak, R.O.; Mikkelsen, J.G. miRNA sponges: Soaking up miRNAs for regulation of gene expression. Wiley Interdiscip. Rev. RNA 2013, 5, 317–333. [Google Scholar] [CrossRef]
  12. Bhan, A.; Soleimani, M.; Mandal, S.S. Long noncoding RNA and cancer: A new paradigm. Cancer Res. 2017, 77, 3965–3981. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Gutschner, T.; Hämmerle, M.; Diederichs, S. MALAT1—A paradigm for long noncoding RNA function in cancer. J. Mol. Med. 2013, 91, 791–801. [Google Scholar] [CrossRef] [PubMed]
  14. Brown, J.A.; Bulkley, D.; Wang, J.; Valenstein, M.L.; Yario, T.A.; Steitz, T.A.; Steitz, J.A. Structural insights into the stabilization of MALAT1 noncoding RNA by a bipartite triple helix. Nat. Struct. Mol. Biol. 2014, 21, 633–640. [Google Scholar] [CrossRef] [PubMed]
  15. McCown, P.J.; Wang, M.C.; Jaeger, L.; Brown, J.A. Secondary structural model of human MALAT1 reveals multiple structure–function relationships. Int. J. Mol. Sci. 2019, 20, 5610. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Liu, N.; Zhou, K.I.; Parisien, M.; Dai, Q.; Diatchenko, L.; Pan, T. N6-methyladenosine alters RNA structure to regulate binding of a low-complexity protein. Nucleic Acids Res. 2017, 45, 6051–6063. [Google Scholar] [CrossRef] [Green Version]
  17. Liu, N.; Dai, Q.; Zheng, G.; He, C.; Parisien, M.; Pan, T. N6-methyladenosine-dependent RNA structural switches regulate RNA–protein interactions. Nat. Cell Biol. 2015, 518, 560–564. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Arun, G.; Aggarwal, D.; Spector, D.L. MALAT1 long non-coding RNA: Functional implications. Non Coding RNA 2020, 6, 22. [Google Scholar] [CrossRef]
  19. Pan, B.; Yang, J.; Wang, X.; Xu, K.; Ikezoe, T. miR-217 sensitizes chronic myelogenous leukemia cells to tyrosine kinase inhibitors by targeting pro-oncogenic anterior gradient 2. Exp. Hematol. 2018, 68, 80–88. [Google Scholar] [CrossRef]
  20. Sachdeva, M.; Zhu, S.; Wu, F.; Wu, H.; Walia, V.; Kumar, S.; Elble, R.; Watabe, K.; Mo, Y.Y. p53 represses c-Myc through induction of the tumor suppressor miR-145. Proc. Natl. Acad. Sci. USA 2009, 106, 3207–3212. [Google Scholar] [CrossRef] [Green Version]
  21. Zhao, Y.; Li, C.; Wang, M.; Su, L.; Qu, Y.; Li, J.; Yu, B.; Yan, M.; Yu, Y.; Liu, B.; et al. Decrease of miR-202-3p expression, a novel tumor suppressor, in gastric cancer. PLoS ONE 2013, 8, e69756. [Google Scholar] [CrossRef] [Green Version]
  22. He, J.; Huang, B.; Zhang, K.; Liu, M.; Xu, T. Long non-coding RNA in cervical cancer: From biology to therapeutic opportunity. Biomed. Pharmacother. 2020, 127, 110209. [Google Scholar] [CrossRef]
  23. Pan, J.; Bian, Y.; Cao, Z.; Lei, L.; Pan, J.; Huang, J.; Cai, X.; Lan, X.; Zheng, H. Long noncoding RNA MALAT1 as a candidate serological biomarker for the diagnosis of non-small cell lung cancer: A meta-analysis. Thorac. Cancer 2019, 11, 329–335. [Google Scholar] [CrossRef]
  24. Linder, B.; Grozhik, A.V.; Olarerin-George, A.O.; Meydan, C.; Mason, C.E.; Jaffrey, S.R. Single-nucleotide-resolution mapping of m6A and m6Am throughout the transcriptome. Nat. Methods 2015, 12, 767–772. [Google Scholar] [CrossRef] [PubMed]
  25. Liu, N.; Parisien, M.; Dai, Q.; Zheng, G.; He, C.; Pan, T. Probing N6-methyladenosine RNA modification status at single nucleotide resolution in mRNA and long noncoding RNA. RNA 2013, 19, 1848–1856. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Zheng, Y.; Nie, P.; Peng, N.; He, Z.; Liu, M.; Xie, Y.; Miao, Y.; Zuo, Z.; Ren, J. m6AVar: A database of functional variants involved in m6A modification. Nucleic Acids Res. 2018, 46, D139–D145. [Google Scholar] [CrossRef] [PubMed]
  27. Rouskin, S.; Zubradt, M.; Washietl, S.; Kellis, M.; Weissman, J.S. Genome-wide probing of RNA structure reveals active unfolding of mRNA structures in vivo. Nat. Cell Biol. 2014, 505, 701–705. [Google Scholar] [CrossRef] [Green Version]
  28. Lu, Z.; Gong, J.; Zhang, Q.C. PARIS: Psoralen analysis of RNA interactions and structures with high throughput and resolution. Comput. Biol. 2018, 1649, 59–84. [Google Scholar] [CrossRef] [Green Version]
  29. Lu, Z.; Zhang, Q.C.; Lee, B.; Flynn, R.A.; Smith, M.A.; Robinson, J.T.; Davidovich, C.; Gooding, A.R.; Goodrich, K.J.; Mattick, J.S.; et al. RNA duplex map in living cells reveals higher-order transcriptome structure. Cell 2016, 165, 1267–1279. [Google Scholar] [CrossRef] [Green Version]
  30. Zhao, W.G.; Yu, S.N.; Lu, Z.H.; Ma, Y.H.; Gu, Y.M.; Chen, J. The miR-217 microRNA functions as a potential tumor suppressor in pancreatic ductal adenocarcinoma by targeting KRAS. Carcinogenesis 2010, 31, 1726–1733. [Google Scholar] [CrossRef] [Green Version]
  31. Polakova, K.M.; Lopotová, T.; Klamova, H.; Burda, P.; Trněný, M.; Stopka, T.; Moravcová, J. Expression patterns of microRNAs associated with CML phases and their disease related targets. Mol. Cancer 2011, 10, 41. [Google Scholar] [CrossRef] [Green Version]
  32. Gao, X.; Wan, Z.; Wei, M.; Dong, Y.; Zhao, Y.; Chen, X.; Li, Z.; Qin, W.; Yang, G.; Liu, L. Chronic myelogenous leukemia cells remodel the bone marrow niche via exosome-mediated transfer of miR-320. Theranostics 2019, 9, 5642–5656. [Google Scholar] [CrossRef] [PubMed]
  33. Zhang, H.; Lu, W. LncRNA SNHG12 regulates gastric cancer progression by acting as a molecular sponge of miR-320. Mol. Med. Rep. 2017, 17, 2743–2749. [Google Scholar] [CrossRef] [Green Version]
  34. Nie, Z.Y.; Liu, X.J.; Zhan, Y.; Liu, M.H.; Zhang, X.Y.; Li, Z.Y.; Lu, Y.Q.; Luo, J.M.; Yang, L. miR-140-5p induces cell apoptosis and decreases Warburg effect in chronic myeloid leukemia by targeting SIX1. Biosci. Rep. 2019, 39, 1–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Shuai, S.; Suzuki, H.; Diaz-Navarro, A.; Nadeu, F.; Kumar, S.A.; Gutierrez-Fernandez, A.; Delgado, J.; Pinyol, M.; López-Otín, C.; Puente, X.S.; et al. The U1 spliceosomal RNA is recurrently mutated in multiple cancers. Nat. Cell Biol. 2019, 574, 712–716. [Google Scholar] [CrossRef] [PubMed]
  36. Latorre, E.; Carelli, S.; Raimondi, I.; D’Agostino, V.; Castiglioni, I.; Zucal, C.; Moro, G.; Luciani, A.; Ghiraldi, G.; Monti, E.; et al. The ribonucleic complex HuR-MALAT1 represses CD133 expression and suppresses epithelial-mesenchymal transition in breast cancer. Cancer Res. 2016, 76, 2626–2636. [Google Scholar] [CrossRef] [Green Version]
  37. Srikantan, S.; Gorospe, M. HuR function in disease. Front. Biosci. 2012, 17, 189–205. [Google Scholar] [CrossRef] [Green Version]
  38. Dai, Q.; Moshitch-Moshkovitz, S.; Han, D.; Kol, N.; Amariglio, N.; Rechavi, G.; Dominissini, D.; He, C. Nm-seq maps 2′-O-methylation sites in human mRNA with base precision. Nat. Methods 2017, 14, 695–698. [Google Scholar] [CrossRef]
  39. Zhang, L.S.; Liu, C.; Ma, H.; Dai, Q.; Sun, H.L.; Luo, G.; Zhang, Z.; Zhang, L.; Hu, L.; Dong, X.; et al. Transcriptome-wide mapping of internal N7-methylguanosine methylome in mammalian mRNA. Mol. Cell 2019, 74, 1304–1316.e8. [Google Scholar] [CrossRef]
  40. Jacob, R.; Zander, S.; Gutschner, T. The dark side of the epitranscriptome: Chemical modifications in long non-coding RNAss. Int. J. Mol. Sci. 2017, 18, 2387. [Google Scholar] [CrossRef] [Green Version]
  41. Squires, J.E.; Patel, H.R.; Nousch, M.; Sibbritt, T.; Humphreys, D.T.; Parker, B.J.; Suter, C.M.; Preiss, T. Widespread occurrence of 5-methylcytosine in human coding and non-coding RNA. Nucleic Acids Res. 2012, 40, 5023–5033. [Google Scholar] [CrossRef]
  42. Liu, J.; Yue, Y.; Han, D.; Wang, X.; Fu, Y.; Zhang, L.; Jia, G.; Yu, M.; Lu, Z.; Deng, X.; et al. A METTL3–METTL14 complex mediates mammalian nuclear RNA N6-adenosine methylation. Nat. Chem. Biol. 2014, 10, 93–95. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Chen, X.; Zhang, J.; Zhu, J.S. The role of m6A RNA methylation in human cancer. Mol. Cancer 2019, 18, 1–9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Zhu, S.; Lu, Y. Dexmedetomidine suppressed the biological behavior of HK-2 cells treated with LPS by down-regulating ALKBH5. Inflammation 2020, 43, 2256–2263. [Google Scholar] [CrossRef]
  45. Berulava, T.; Ziehe, M.; Klein-Hitpass, L.; Mladenov, E.; Thomale, J.; Ruther, U.; Horsthemke, B. FTO levels affect RNA modification and the transcriptome. Eur. J. Hum. Genet. 2012, 21, 317–323. [Google Scholar] [CrossRef] [Green Version]
  46. Xu, C.; Liu, K.; Tempel, W.; Demetriades, M.; Aik, W.; Schofield, C.J.; Min, J. Structures of human ALKBH5 demethylase reveal a unique binding mode for specific single-stranded N6-methyladenosine RNA demethylation. J. Biol. Chem. 2014, 289, 17299–17311. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Yeo, G. The role of the FTO (Fat Mass and Obesity Related) locus in regulating body size and composition. Mol. Cell. Endocrinol. 2014, 397, 34–41. [Google Scholar] [CrossRef]
  48. Moursy, A.; Allain, F.H.T.; Cléry, A. Characterization of the RNA recognition mode of hnRNPG extends its role in SMN2 splicing regulation. Nucleic Acids Res. 2014, 42, 6659–6672. [Google Scholar] [CrossRef] [Green Version]
  49. Pal, A.; Levy, Y. Structure, stability and specificity of the binding of ssDNA and ssRNA with proteins. PLoS Comput. Biol. 2019, 15, e1006768. [Google Scholar] [CrossRef]
  50. Shin, K.H.; Kim, R.H.; Kim, R.H.; Kang, M.K.; Park, N.H. hnRNPG elicits tumor-suppressive activity in part by upregulating the expression of Txnip. Biochem. Biophys. Res. Commun. 2008, 372, 880–885. [Google Scholar] [CrossRef] [Green Version]
  51. Shin, K.H.; Kang, M.K.; Kim, R.H.; Christensen, R.; Park, N.H. Heterogeneous nuclear ribonucleoprotein G shows tumor suppressive effect against oral squamous cell carcinoma cells. Clin. Cancer Res. 2006, 12, 3222–3228. [Google Scholar] [CrossRef] [Green Version]
  52. Shin, K.; Kim, R.H.; Yu, B.; Kang, M.K.; Elashoff, D.; Christensen, R.; Pucar, A.; Park, N.H. Expression and mutation analysis of heterogeneous nuclear ribonucleoprotein G in human oral cancer. Oral Oncol. 2011, 47, 1011–1016. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Dai, D.; Wang, H.; Zhu, L.; Jin, H.; Wang, X. N6-methyladenosine links RNA metabolism to cancer progression. Cell Death Dis. 2018, 9, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Wang, J.Z.; Xiang, J.J.; Wu, L.G.; Bai, Y.S.; Chen, Z.W.; Yin, X.Q.; Wang, Q.; Guo, W.H.; Peng, Y.; Guo, H.; et al. A genetic variant in long non-coding RNA MALAT1 associated with survival outcome among patients with advanced lung adenocarcinoma: A survival cohort analysis. BMC Cancer 2017, 17, 1–8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Ritz, J.; Martin, J.S.; Laederach, A. Evaluating our ability to predict the structural disruption of RNA by SNPs. BMC Genom. 2012, 13 (Suppl. 4), 1–11. [Google Scholar] [CrossRef] [Green Version]
  56. Sabarinathan, R.; Tafer, H.; Seemann, S.E.; Hofacker, I.L.; Stadler, P.F.; Gorodkin, J. RNAsnp: Efficient detection of local RNA secondary structure changes induced by SNPs. Hum. Mutat. 2013, 34, 925. [Google Scholar] [CrossRef] [Green Version]
  57. Wu, S.; Sun, H.; Wang, Y.; Yang, X.; Meng, Q.; Yang, H.; Zhu, H.; Tang, W.; Li, X.; Aschner, M.; et al. MALAT1 rs664589 polymorphism inhibits binding to miR-194-5p, contributing to colorectal cancer risk, growth, and metastasis. Cancer Res. 2019, 79, 5432–5441. [Google Scholar] [CrossRef] [Green Version]
  58. Grossman, R.L.; Heath, A.P.; Ferretti, V.; Varmus, H.E.; Lowy, D.R.; Kibbe, W.A.; Staudt, L.M. Toward a shared vision for cancer genomic data. N. Engl. J. Med. 2016, 375, 1109–1112. [Google Scholar] [CrossRef]
  59. Song, B.; Yan, J.; Liu, C.; Zhou, H.; Zheng, Y. Tumor suppressor role of miR-363-3p in gastric cancer. Med. Sci. Monit. 2015, 21, 4074–4080. [Google Scholar] [CrossRef] [Green Version]
  60. Casadei, L.; Calore, F.; Creighton, C.J.; Guescini, M.; Batte, K.; Iwenofu, O.H.; Zewdu, A.; Braggio, D.A.; Bill, K.L.; Fadda, P.; et al. Exosome-derived miR-25-3p and miR-92a-3p stimulate liposarcoma progression. Cancer Res. 2017, 77, 3846–3856. [Google Scholar] [CrossRef] [Green Version]
  61. Wan, Z.; Chen, X.; Gao, X.; Dong, Y.; Zhao, Y.; Wei, M.; Fan, W.; Yang, G.; Liu, L. Chronic myeloid leukemia-derived exosomes attenuate adipogenesis of adipose derived mesenchymal stem cells via transporting miR-92a-3p. J. Cell. Physiol. 2019, 234, 21274–21283. [Google Scholar] [CrossRef]
  62. Willmann, M.R.; Berkowitz, N.D.; Gregory, B.D. Improved genome-wide mapping of uncapped and cleaved transcripts in eukaryotes—GMUCT 2.0. Methods 2014, 67, 64–73. [Google Scholar] [CrossRef] [PubMed]
  63. Ma, L.; Li, L.L. miR-145 Contributes to the progression of cervical carcinoma by directly regulating FSCN1. Cell Transplant. 2019, 28, 1299–1305. [Google Scholar] [CrossRef] [Green Version]
  64. Lu, H.; He, Y.; Lin, L.; Qi, Z.; Ma, L.; Li, L.; Su, Y. Long non-coding RNA MALAT1 modulates radiosensitivity of HR-HPV+ cervical cancer via sponging miR-145. Tumor Biol. 2015, 37, 1683–1691. [Google Scholar] [CrossRef] [PubMed]
  65. Wang, X.; Tang, S.; Le, S.Y.; Lu, R.; Rader, J.S.; Meyers, C.; Zheng, Z.M. Aberrant expression of oncogenic and tumor-suppressive microRNAs in cervical cancer is required for cancer cell growth. PLoS ONE 2008, 3, e2557. [Google Scholar] [CrossRef] [Green Version]
  66. Zeng, F.; Xue, M.; Xiao, T.; Li, Y.; Xiao, S.; Jiang, B.; Ren, C. MiR-200b promotes the cell proliferation and metastasis of cervical cancer by inhibiting FOXG1. Biomed. Pharmacother. 2016, 79, 294–301. [Google Scholar] [CrossRef]
  67. Kang, H.W.; Wang, F.; Wei, Q.; Zhao, Y.F.; Liu, M.; Li, X.; Tang, H. miR-20a promotes migration and invasion by regulating TNKS2 in human cervical cancer cells. FEBS Lett. 2012, 586, 897–904. [Google Scholar] [CrossRef] [Green Version]
  68. Piao, J.; You, K.; Guo, Y.; Zhang, Y.Y.; Li, Z.; Geng, L. Substrate stiffness affects epithelial-mesenchymal transition of cervical cancer cells through miR-106b and its target protein DAB2. Int. J. Oncol. 2017, 50, 2033–2042. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Fan, Y.; Sheng, W.; Meng, Y.; Cao, Y.; Li, R. LncRNA PTENP1 inhibits cervical cancer progression by suppressing miR-106b. Artif. Cells Nanomed. Biotechnol. 2020, 48, 393–407. [Google Scholar] [CrossRef] [PubMed]
  70. Engreitz, J.M.; Sirokman, K.; McDonel, P.; Shishkin, A.A.; Surka, C.; Russell, P.; Grossman, S.R.; Chow, A.Y.; Guttman, M.; Lander, E.S. RNA-RNA interactions enable specific targeting of noncoding RNAs to nascent pre-mRNAs and chromatin sites. Cell 2014, 159, 188–199. [Google Scholar] [CrossRef] [Green Version]
  71. Aw, J.G.A.; Shen, Y.; Wilm, A.; Sun, M.; Ni Lim, X.; Boon, K.L.; Tapsin, S.; Chan, Y.S.; Tan, C.P.; Sim, A.Y.; et al. In vivo mapping of eukaryotic RNA interactomes reveals principles of higher-order organization and regulation. Mol. Cell 2016, 62, 603–617. [Google Scholar] [CrossRef] [Green Version]
  72. Woodcock, C.B.; Yu, D.; Hajian, T.; Li, J.; Huang, Y.; Dai, N.; Correa, I.R., Jr.; Wu, T.; Vedadi, M.; Zhang, X.; et al. Human MettL3–MettL14 complex is a sequence-specific DNA adenine methyltransferase active on single-strand and unpaired DNA in vitro. Cell Discov. 2019, 5, 1–4. [Google Scholar]
  73. Amort, T.; Soulière, M.F.; Wille, A.; Jia, X.Y.; Fiegl, H.; Wörle, H.; Micura, R.; Lusser, A. Long non-coding RNAs as targets for cytosine methylation. RNA Biol. 2013, 10, 1002–1008. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. West, J.A.; Davis, C.P.; Sunwoo, H.; Simon, M.D.; Sadreyev, R.I.; Wang, P.I.; Tolstorukov, M.Y.; Kingston, R.E. The long noncoding RNAs NEAT1 and MALAT1 bind active chromatin sites. Mol. Cell 2014, 55, 791–802. [Google Scholar] [CrossRef] [Green Version]
  75. McCown, P.J.; Ruszkowska, A.; Kunkler, C.N.; Breger, K.; Hulewicz, J.P.; Wang, M.C.; Springer, N.A.; Brown, J.A. Naturally occurring modified ribonucleosides. Wiley Interdiscip. Rev. RNA 2020, 11, e1595. [Google Scholar] [CrossRef] [PubMed]
  76. Yang, C.; Dou, R.; Yin, T.L.; Ding, J. MiRNA-106b-5p in human cancers: Diverse functions and promising biomarker. Biomed. Pharmacother. 2020, 127, 110211. [Google Scholar] [CrossRef]
  77. Cai, N.; Hu, L.; Xie, Y.; Gao, J.H.; Zhai, W.; Wang, L.; Jin, Q.J.; Qin, C.Y.; Qiang, R. MiR-17-5p promotes cervical cancer cell proliferation and metastasis by targeting transforming growth factor-β receptor 2. Eur. Rev. Med. Pharmacol. Sci. 2018, 22, 1899–1906. [Google Scholar]
  78. Afgan, E.; Baker, D.; Batut, B.; van den Beek, M.; Bouvier, D.; Čech, M.; Chilton, J.; Clements, D.; Coraor, N.; Gruning, B.A.; et al. The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2018 update. Nucleic Acids Res. 2018, 46, W537–W544. [Google Scholar] [CrossRef] [Green Version]
  79. Kent, W.J.; Sugnet, C.W.; Furey, T.S.; Roskin, K.M.; Pringle, T.H.; Zahler, A.M.; Haussler, A.D. The Human Genome Browser at UCSC. Genome Res. 2002, 12, 996–1006. [Google Scholar] [CrossRef] [Green Version]
  80. Darty, K.; Denise, A.; Ponty, Y. VARNA: Interactive drawing and editing of the RNA secondary structure. Bioinformatics 2009, 25, 1974–1975. [Google Scholar] [CrossRef] [Green Version]
  81. Li, J.H.; Liu, S.; Zhou, H.; Qu, L.H.; Yang, J.H. StarBase v2.0: Decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. 2014, 42, D92–D97. [Google Scholar] [CrossRef] [Green Version]
  82. Cai, Z.; Cao, C.; Ji, L.; Ye, R.; Wang, D.; Xia, C.; Wang, S.; Du, Z.; Hu, N.; Yu, X.; et al. RIC-seq for global in situ profiling of RNA–RNA spatial interactions. Nat. Cell Biol. 2020, 582, 432–437. [Google Scholar] [CrossRef] [PubMed]
  83. Wilusz, J.E.; Freier, S.M.; Spector, D.L. 3′ End processing of a long nuclear-retained noncoding RNA yields a tRNA-like cytoplasmic RNA. Cell 2008, 135, 919–932. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Lv, D.; Xu, K.; Jin, X.; Li, J.; Shi, Y.; Zhang, M.; Jin, X.; Li, Y.; Xu, J.; Li, X. LncSpA: LncRNA spatial atlas of expression across normal and cancer tissues. Cancer Res. 2020, 80, 2067–2071. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Safra, M.; Sas-Chen, A.; Nir, R.; Winkler, R.; Nachshon, A.; Bar-Yaacov, D.; Erlacher, M.; Rossmanith, W.; Stern-Ginossar, N.; Schwartz, S. The m1A landscape on cytosolic and mitochondrial mRNA at single-base resolution. Nat. Cell Biol. 2017, 551, 251–255. [Google Scholar] [CrossRef]
Ncrna 07 00006 g001 550
Figure 1. Summary of DMS-Seq data analysis. (A) 2311 available DMS-Seq datapoints (i.e., number of DMS-Seq reads corresponding to a single adenosine or cytidine) in human foreskin fibroblasts and 2554 DMS-Seq datapoints in K562 cells were classified as unstructured or structured based on whether they exceeded the 20-read or 250-read threshold, respectively (see Section 3 and [15]). (B) K562 DMS-Seq datapoints were compared to the working MALAT1 consensus model and further classified as losing structure, gaining structure, or agreeing with the model [15].
Figure 1. Summary of DMS-Seq data analysis. (A) 2311 available DMS-Seq datapoints (i.e., number of DMS-Seq reads corresponding to a single adenosine or cytidine) in human foreskin fibroblasts and 2554 DMS-Seq datapoints in K562 cells were classified as unstructured or structured based on whether they exceeded the 20-read or 250-read threshold, respectively (see Section 3 and [15]). (B) K562 DMS-Seq datapoints were compared to the working MALAT1 consensus model and further classified as losing structure, gaining structure, or agreeing with the model [15].
Ncrna 07 00006 g001
Ncrna 07 00006 g002 550
Figure 2. Map of K562 DMS-Seq structural changes onto the MALAT1 secondary structural model. K562 DMS-Seq data were available for nts 1284–4257 and nts 5839–8425. Of those indicating change in structure relative to the noncancerous consensus model, adenosine and cytidine residues that are predicted to lose structure are depicted in orange (733 datapoints) and those residues that are predicted to gain structure are depicted in blue (317 datapoints). Orange marks in hairpins have been extended to cover both nucleotides in a given base pair, and base pairs marked as unstructured have the two nucleotides moved apart to emphasize loss of structure. Labels for secondary structures correspond to those established for the working MALAT1 model in noncancerous cells [15]. Please note that PARIS-derived hairpins appear in the consensus model, although these hairpins were not considered in the differential structural analysis of K562-MALAT1.
Figure 2. Map of K562 DMS-Seq structural changes onto the MALAT1 secondary structural model. K562 DMS-Seq data were available for nts 1284–4257 and nts 5839–8425. Of those indicating change in structure relative to the noncancerous consensus model, adenosine and cytidine residues that are predicted to lose structure are depicted in orange (733 datapoints) and those residues that are predicted to gain structure are depicted in blue (317 datapoints). Orange marks in hairpins have been extended to cover both nucleotides in a given base pair, and base pairs marked as unstructured have the two nucleotides moved apart to emphasize loss of structure. Labels for secondary structures correspond to those established for the working MALAT1 model in noncancerous cells [15]. Please note that PARIS-derived hairpins appear in the consensus model, although these hairpins were not considered in the differential structural analysis of K562-MALAT1.
Ncrna 07 00006 g002
Ncrna 07 00006 g003 550
Figure 3. Hairpins predicted to be disrupted in K562-MALAT1. Loss of structure is shown for four hairpins: (A) H44, (B) H45, (C) H98, and (D) H155. Canonical base pairs are shown with a line and noncanonical base pairs are shown with a dot. Nucleotides with an orange background are classified as unstructured and nucleotides with a blue background are classified as structured. For each hairpin, the following schematics are shown from top to bottom: original secondary structure as predicted in the working MALAT1 model (“Model”) using fibroblast DMS-Seq data, the K562 DMS-Seq data (“K562-MALAT1”), and the K562 DMS-Seq datapoints that differ from their noncancerous counterparts (“Change”).
Figure 3. Hairpins predicted to be disrupted in K562-MALAT1. Loss of structure is shown for four hairpins: (A) H44, (B) H45, (C) H98, and (D) H155. Canonical base pairs are shown with a line and noncanonical base pairs are shown with a dot. Nucleotides with an orange background are classified as unstructured and nucleotides with a blue background are classified as structured. For each hairpin, the following schematics are shown from top to bottom: original secondary structure as predicted in the working MALAT1 model (“Model”) using fibroblast DMS-Seq data, the K562 DMS-Seq data (“K562-MALAT1”), and the K562 DMS-Seq datapoints that differ from their noncancerous counterparts (“Change”).
Ncrna 07 00006 g003
Ncrna 07 00006 g004 550
Figure 4. Structural loss in K562-MALAT1 improves accessibility of miRNA-binding sites. Three hairpin-miRNA interactions are shown: (A) miR-320 seed binding to H101, (B) miR-217 seed binding to H160, and (C) miR-140-5p seed binding to H168. Depicted is the predicted structure of each hairpin in noncancerous cells (“Model”) adjacent to the hairpin structure in K562-MALAT1 (“K562”), where the relevant base pairs are removed according to our analysis of DMS-Seq data. Nucleotides with an orange background are classified as unstructured and nucleotides with a blue background are classified as structured. The binding of a given miRNA (red) seed region to the K562-MALAT1 hairpin is shown. The nucleotide R denotes a G and A residue in miR-320a and miR-320b, respectively.
Figure 4. Structural loss in K562-MALAT1 improves accessibility of miRNA-binding sites. Three hairpin-miRNA interactions are shown: (A) miR-320 seed binding to H101, (B) miR-217 seed binding to H160, and (C) miR-140-5p seed binding to H168. Depicted is the predicted structure of each hairpin in noncancerous cells (“Model”) adjacent to the hairpin structure in K562-MALAT1 (“K562”), where the relevant base pairs are removed according to our analysis of DMS-Seq data. Nucleotides with an orange background are classified as unstructured and nucleotides with a blue background are classified as structured. The binding of a given miRNA (red) seed region to the K562-MALAT1 hairpin is shown. The nucleotide R denotes a G and A residue in miR-320a and miR-320b, respectively.
Ncrna 07 00006 g004
Ncrna 07 00006 g005 550
Figure 5. A2875U disrupts the miR-92-3p binding site on H71. (A) A schematic is shown for the predicted structure of wild type (WT) H71 in noncancerous cells [15]. Nucleotides with an orange background are classified as unstructured and nucleotides with a blue background are classified as structured. The expected binding of miR-92a-3p (red) to H71 in K562-MALAT1 is shown for (B) wild type H71 and (C) H71 with the A2875U mutation, which is outlined in yellow. The same disruption to seed binding is expected for miR-363-3p and miR-25-3p (Supplemental Table S1, “Seq-Markup” tab).
Figure 5. A2875U disrupts the miR-92-3p binding site on H71. (A) A schematic is shown for the predicted structure of wild type (WT) H71 in noncancerous cells [15]. Nucleotides with an orange background are classified as unstructured and nucleotides with a blue background are classified as structured. The expected binding of miR-92a-3p (red) to H71 in K562-MALAT1 is shown for (B) wild type H71 and (C) H71 with the A2875U mutation, which is outlined in yellow. The same disruption to seed binding is expected for miR-363-3p and miR-25-3p (Supplemental Table S1, “Seq-Markup” tab).
Ncrna 07 00006 g005
Ncrna 07 00006 g006 550
Figure 6. Map of HeLa PARIS data onto the MALAT1 secondary structural model. Short-range (or local) PARIS interactions are depicted with 25 nucleotides flanking the datapoint. Short-range PARIS interactions diverging from the model are dark blue lines, whereas those that are consistent with the consensus model are light blue lines. Long-range PARIS interactions diverging from the model are depicted as red lines and show 50 nucleotides flanking the datapoints. All datapoint coordinates are in Supplemental Table S2. PKs expected-to-be-lost are depicted with dashed orange lines. Labels for secondary structures correspond to those established for the working noncancerous MALAT1 model [15].
Figure 6. Map of HeLa PARIS data onto the MALAT1 secondary structural model. Short-range (or local) PARIS interactions are depicted with 25 nucleotides flanking the datapoint. Short-range PARIS interactions diverging from the model are dark blue lines, whereas those that are consistent with the consensus model are light blue lines. Long-range PARIS interactions diverging from the model are depicted as red lines and show 50 nucleotides flanking the datapoints. All datapoint coordinates are in Supplemental Table S2. PKs expected-to-be-lost are depicted with dashed orange lines. Labels for secondary structures correspond to those established for the working noncancerous MALAT1 model [15].
Ncrna 07 00006 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, M.C.; McCown, P.J.; Schiefelbein, G.E.; Brown, J.A. Secondary Structural Model of MALAT1 Becomes Unstructured in Chronic Myeloid Leukemia and Undergoes Structural Rearrangement in Cervical Cancer. Non-Coding RNA 2021, 7, 6. https://doi.org/10.3390/ncrna7010006

AMA Style

Wang MC, McCown PJ, Schiefelbein GE, Brown JA. Secondary Structural Model of MALAT1 Becomes Unstructured in Chronic Myeloid Leukemia and Undergoes Structural Rearrangement in Cervical Cancer. Non-Coding RNA. 2021; 7(1):6. https://doi.org/10.3390/ncrna7010006

Chicago/Turabian Style

Wang, Matthew C., Phillip J. McCown, Grace E. Schiefelbein, and Jessica A. Brown. 2021. "Secondary Structural Model of MALAT1 Becomes Unstructured in Chronic Myeloid Leukemia and Undergoes Structural Rearrangement in Cervical Cancer" Non-Coding RNA 7, no. 1: 6. https://doi.org/10.3390/ncrna7010006

APA Style

Wang, M. C., McCown, P. J., Schiefelbein, G. E., & Brown, J. A. (2021). Secondary Structural Model of MALAT1 Becomes Unstructured in Chronic Myeloid Leukemia and Undergoes Structural Rearrangement in Cervical Cancer. Non-Coding RNA, 7(1), 6. https://doi.org/10.3390/ncrna7010006

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop