Ago HITS-CLIP expands microRNA-mRNA interactions in nucleus and cytoplasm of gastric cancer cells
BMC Cancer volume 19, Article number: 29 (2019)
Intensive investigations have identified a collection of microRNAs (miRNAs) and their functional machineries in cytoplasm. However, a comprehensive view of miRNAs and mRNAs in cytoplasm and nucleus has not been explored. This study aims to reveal the mechanisms of miRNA-RNA interactions in nucleus and cytoplasm.
In this study, the miRNAs and their target mRNAs in the Argonaute2 (Ago2) complex of nucleus and cytoplasm of gastric cancer cells were characterized using high-throughput sequencing of RNAs isolated by crosslinking immunoprecipitation (HITS-CLIP). Then, the selected miRNAs were verified by Northern blot. The target mRNAs in the Argonaute2 (Ago2) complex of nucleus and cytoplasm of gastric cancer cells were analyzed through Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) analysis.
The results revealed that there were 243 miRNAs and 265 miRNAs in the Ago2 complexes of nucleus and cytoplasm, respectively. The majority of mature miRNAs existed in cytoplasm. The analysis of miRNA targetome from the Ago2 complexes indicated that a lot of mRNAs with high expression level existed in nucleus. The target genes of miRNAs in the Ago2 complexes of nucleus and cytoplasm played important roles in cell proliferation, cell differentiation, innate immune response and tumorigenesis.
microRNA-mRNA interactions occur in nucleus and cytoplasm of gastric cancer cells.
Therefore, our study demonstrated that miRNA-mRNA interactions not only took place in cytoplasm but also in nucleus.
MicroRNA (miRNA) as one kind of non-coding small RNA (~ 22 nucleotide) plays an important role in the regulation of gene expression [1, 2]. After the primary miRNA (pri-miRNA) is transcribed by RNA Pol II, it needs a series of processing to form the mature miRNA [1, 2]. The Drosha protein (~ 160 kDa) cleaves pri-miRNAs to initiate miRNA processing and maturation, and forms ~ 65 nt (nucleotide) precursor miRNAs with hairpin structures (pre-miRNAs) . Then, the Dicer protein can act on the pre-miRNA and cleave the end loop of the pre-miRNA to generate a double-stranded RNA with a length of approximately 22 base pairs (bp) [1, 2]. The double-stranded RNA duplex leads RNA-induced silencing complex (RISC), containing a member of the conserved Argonaute (Ago) protein family, to target sites of mRNAs, resulting in the destabilization of the mRNAs and/or inhibition of translation [4, 5]. Generally, miRNA function as gene regulator mainly by targeting 3’UTR of mRNAs through the “seed” sequence (nucleotides 2–7) of miRNA [4,5,6]. Intensive investigations during the past decades have established miRNAs as key regulators of various biological processes including organ development, cell differentiation, cell proliferation, virus-host interactions and tumorigenesis [1, 2, 7]. However, more and more studies demonstrated that the function of miRNA needs to be further elucidated. Some researches reported that mature miRNAs could also be detected in nucleus . Meister G found that miR-21 existed in both cytoplasm and nucleus . As reported, miR-25 and miR-92a show clear nuclear enrichment in rat primary cortical neurons by using a combination of microarray analysis and small RNA deep sequencing . Some reports have pointed out the existence of active RISC in nucleus because of the findings of human Ago proteins and other RNAi factors (DICER1, TARBP2 and GW182) in nucleus of cells [11,12,13,14]. However, the miRNAs in nucleus have not been intensively explored.
Generally, the interaction between a miRNA and its target requires base pairing with only 6–8 nucleotides. Therefore, it is a major challenge to predict target mRNAs. Based on evolutionarily conserved seed sequence and free binding energy, computational predictions have greatly improved the potentials to find bona fide miRNA target sites . Nonetheless, different prediction algorithms produce a large number of divergent results for one given miRNA [16, 17]. In addition, computer programs usually predict targets on a genome-wide scale irrespective of mRNA and miRNA that are specific to tissues or cell lines. Thus, it has become increasingly important to develop biochemical tools based on physical interactions between miRNAs and their specific target RNAs. As miRNA target sites are parts of RISC that contain, except for mRNAs, Argonaute (Ago) protein and guiding miRNAs, co-immunoprecipitation (Co-IP) of Ago protein together with associated target mRNAs is an effective tool for the identification or validation of miRNA target candidates [18, 19].
In order to explore the miRNA-mRNA interaction in nucleus and cytoplasm, the miRNAs and mRNAs in the Ago complexes of nucleus and cytoplasm of gastric cancer cells (MGC-803) were characterized using high-throughput sequencing of RNAs isolated by crosslinking immunoprecipitation (HITS-CLIP) in this study. The results revealed that there existed 243 miRNAs and 265 miRNAs in Ago2 complex of nucleus and cytoplasm, respectively. The targetome analysis indicated that these miRNAs could function in nucleus or cytoplasm.
Human gastric cancer cells (MGC-803) were cultured in RPIM-1640 Medium (Gibco, USA) supplemented with 10% fetal bovine serum (FBS) at 37 °C. The cell line was purchased from The Cell Bank of Type Culture Collection of Chinese Academy of Sciences, Shanghai, China in the year of 2014.
Isolation of nucleus and cytoplasm
MGC-803 cells were cross-linked at 254 nm and then washed twice with ice-cold phosphate-buffered saline (PBS). Subsequently, the cells were warmed at room for 45 s and lysed by TM-2 buffer [10 mM Tris-HCl, 2 mM MgCl2, phenylmethanesulfonyl fluoride (PMSF), pH 7.4] at room temperature for 1 min. To isolate nuclei and cytoplasm of cells, the lysate was centrifugated at 1000×g at 4 °C for 6 min.
Western blot analysis
8% SDS-PAGE was used to separate the protein samples. Then the protein samples were transferred to a nitrocellulose membrane (Bio-Rad, USA) at 80 V for 1.5 h. The membrane was blocked with 3% milk at room temperature for 2 h. Subsequently, the membrane was incubated with a primary antibody (anti-tubulin IgG, anti-Lamin A/C IgG or anti-Ago2 IgG) at 4 °C overnight. Tubulin and Lamin A/C antibodies were purchased from Abcam (USA). The anti-Ago2 IgG was prepared in our laboratory. After washes with TBST, the membrane was immersed in 2% milk containing AP-conjugated goat antimouse IgG (Sigma, USA) for 1 h. At last, NBT and BCIP solutions (BBI, Canada) were used to detect the protein.
Co-immunoprecipitation (co-IP) of Ago2 complex
The lysis buffer (20 mM Tris-HCl, 1 mM EDTA, 150 mM NaCl, 1% Triton X-100, pH 7.5) was prepared to lyse the isolated nuclei and cytoplasm. After treatment with RNAsin (Promega, USA), the lysate was incubated with RQ1 DNAase (Promega) for 5 min at 37°. Subsequently the Ago2 complex was immunoprecipitated using the polyclonal antibody against Ago2 for 3 h at 4 °C and then incubated with protein A-Sepharose (Bio-Rad) for 3 h at 4 °C. After washes with lysis buffer, the Ago2 complex was collected.
Sequencing and sequence analysis of RNAs of Ago2 complex
RNA was extracted from immunoprecipitated Ago2 complex and subjected to sequencing using GA-I genome analyzer (Illumina, San Diego, CA). To analyze small RNAs, the ACGT V3.1 program developed by LC Sciences (Houston, USA) was employed. Then the high-quality sequences were obtained by removing adaptor sequences, mRNA, rRNA, tRNA, snRNA, snoRNA and other noncoding RNA sequences available in Rfam (http://www.sanger.ac.uk/software/Rfam). The high quality sequences were then blasted with known human miRNAs in miRBase 19.0 to obtain miRNAs in this study. To analyze the longer RNAs, trinity software and the Illumina paired end method developed by LC Sciences (Houston, USA) were employed. Then the high-quality reads were aligned to human mRNA database (ftp://ftp.ensembl.org/pub/release-70/fasta/homo_sapiens/).
Total RNAs were extracted from the Ago2 complex of nucleus or cytoplasm using a mirVanaPTMP miRNA isolation kit according to the manufacturer’s protocol (Ambion, USA) and then separated on a denaturing 15% polyacrylamide gel containing 8 M urea. After transferred to a Hybond-N+ membrane (Amersham Biosciences, Buckinghamshire, UK), the RNAs were cross-linked under ultraviolet. The membrane was pre-incubated with DIG Easy Hyb granules buffer (Roche, Basel, Switzerland) for 1 h at 42 °C and then the miRNAs were probed against with complementary DIG-labeled miRNA probe at 42 °C overnight. The RNAs were detected using the DIG High Prime DNA Labeling and Detection Starter Kit II (Roche, USA). NBT and BCIP solutions (BBI, Canada) were used to detect the RNA.
Prediction of genes targeted by miRNAs
The transcriptome sequences of the Ago2 complexes of nucleus and cytoplasm served as a database for miRNA target prediction. MiRNA target genes were predicted by employing TargetScan 5.1 (http://www.targetscan.org) and miRanda (http://www.microrna.org/) as described in previous study .
Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) analysis
The coding sequences of transcripts obtained from the Ago2 complexes were aligned with the protein sequences in the GO database with the blast E value <1e− 5 . The best hit GO identities were described the properties of these genes and gene products. The coding sequences of transcripts were analyzed by KEGG Pathway (239 KEGG) to analyze the relationship between genes and enzyme annotations in KEGG, and the information mapped to pathway.
Acquirement of Ago2 complexes from cytoplasm and nucleus of gastric cancer cells
To identify miRNAs and mRNAs targeted by miRNAs in Ago2 complex of cytoplasm and nucleus, co-immunoprecipitation (Co-IP) assays with Ago2-specific antibody were performed using gastric cancer cells (MGC-803). The nuclear and cytoplasmic fractions of gastric cancer cells were extracted. Western blot analysis indicated that β-tubulin, a cytoplasmic marker, was detected in cell lysate and cytoplasm but not nucleus (Fig. 1a). At the same time, Lamin A/C, a nuclear marker, was found in cell lysate and nucleus (Fig. 1a). These data showed that cytoplasm and nucleus were successfully isolated.
The isolated cytoplasm and nucleus were subjected to Co-IP assays using Ago2-specific antibody. The results revealed that the Ago2 complexes of cytoplasm and nucleus were obtained (Fig. 1b).
Characterization of miRNAs in Ago2 complexes of cytoplasm and nucleus of gastric cancer cells
To characterize the miRNAs in Ago2 complexes of cytoplasm and nucleus, RNAs were extracted from Ago2 complexes and then subjected to deep sequencing. After removal of mRNA, rRNA, tRNA, snRNA and snoRNA sequences, high-throughput small RNA sequencing yielded an average of 7,813,944 high-quality reads. The results showed that most of small RNA reads were 20–25 nucleotides (nt) in length, which is typical for products generated by the enzyme Dicer (Fig. 2a). Of high-quality small RNA reads, 4,697,199 reads in nucleus and 4,545,958 in cytoplasm were mapped to known miRNAs in miRBase 19.0. In total, 243 miRNAs in the Ago2 complex of nucleus and 265 miRNAs in the Ago2 complex of cytoplasm were identified (Fig. 2b, Table 1), indicating that there was no difference of miRNA number between nucleus and cytoplasm of gastric cancer cells. The sequencing analysis revealed that there was a strong preference of U at the first nucleotide of small RNAs’ 5′ ends from the Ago2 complexes (Fig. 2c), showing that these small RNAs were in line with the characteristics of mature miRNA.
To confirm the existence of mature miRNAs in the Ago2 complexes of cytoplasm and nucleus, 5 miRNAs were randomly selected for Northern blot analysis. The results of Northern blots were consistent with those of miRNA sequencing (Fig. 2e).
Identification of miRNA targetome in Ago2 complexes of nucleus and cytoplasm of gastric cancer cells
To obtain an overview of gene expression profiles of miRNA targetome in cytoplasm and nucleus, cDNA libraries of target mRNA tags of miRNAs were prepared from the Ago2 complexes of cytoplasm and nucleus. After removal of repetitive and low-quality reads, an average of 12.2 million high-quality reads were aligned to human mRNA database. A total of 191,134 unigenes in nucleus and 191,134 unigenes in cytoplasm were obtained. The results showed that a large proportion of the target mRNA reads, 64.7% in nucleus and 55% in cytoplasm, were mapped to human mRNA database. Among the top 100 mRNAs with higher expression level, 80 genes were downregulated in cytoplasm and 20 upregulated (Fig. 3a), indicating that the main proportion of mRNAs existed in nucleus.
The analysis of cellular component by gene ontology (GO) revealed that the differentially expressed genes could be categorized into molecular function, cellular component and biological process (Fig. 3b). In the molecular function category, the most common functional annotations were the zinc ion binding, nucleotide binding, ATP binding and metal ion binding (Fig. 3b). A lot of differentially expressed genes could function as cellular components of cytoplasm and nucleus (Fig. 3b). In terms of biological process, most of genes were implicated in signal transduction, multicellular organismal development and translational elongation (Fig. 3b). The results of KEGG analysis indicated that the differentially expressed genes were classified into 239 KEGG pathways. Many genes were involved in the MAPK signaling pathway, the Wnt signaling pathway, focal adhesion, the TGF-beta signaling pathway and pathways in cancer (Fig. 3c). These data showed that the target genes of miRNAs in the Ago2 complexes of nucleus and cytoplasm played important roles in cell proliferation, cell differentiation, innate immune response and tumorigenesis.
Pathways mediated by miRNAs with higher expression level
In order to characterize the potential interactions between miRNAs and their targets in cytoplasm and nucleus, 10 miRNAs with higher expression level in nucleus and cytoplasm (miR-125b-5p, miR-92b-3p, miR-17-5p, miR-30b-5p, miR-20a-5p, miR-21-3p, let-7d-3p, miR-98-5p, miR-424-5p and miR-155-5p) and their predicted targets were subjected to GO and KEGG analyses. The results of GO analysis demonstrated that most of the miRNA targets were involved in single-organism process, positive regulation of biological activity and cellular component including cytoplasm, intracellular organelle and organelle (Fig. 4a). The data of KEGG analysis revealed that most target genes participated in axon guidance, natural killer cell mediated cytotoxity and T cell receptor signaling pathway.
Over a decade of years, investigations have identified a collection of miRNA functional machineries. Recently, more and more investigations focus on miRNA functions in nucleus. It is found that miR-9 can result in degradation of MALAT-1 mRNA by targeting the metastasis-associated lung adenocarcinoma transcript 1 (MALAT-1) in nucleus . MiR-671 could regulate the expression level of the circular non-coding (ncRNA) natural antisense transcript in nucleus . These studies demonstrate that miRNAs in nucleus can also regulate the expressions of genes. Consistent with previous studies, our study characterized 243 miRNAs and 265 miRNAs in the Ago2 complexes of nucleus and cytoplasm of gastric cancer cells in this study. Report has shown that miRNA could translocate from cytoplasm to nucleus, which has been demonstrated by using superquencher molecular probes complementary to the mature miR-122 . And RISC-associated proteins are involved in transport from the cytoplasm to the nucleus in human cells [24, 25]. For instance, Exportin-1 and importin-8 play positive roles in the nuclear import of miRNAs [24, 25]. However, miRNAs in nucleus and cytoplasm have not been extensively compared. In the present study, the results revealed that the expression pattern of miRNAs in cytoplasm of gastric cancer cells were obviously different from that in nucleus. The majority of mature miRNAs existed in cytoplasm. Therefore, our study contributed a comprehensive view of miRNAs and their targets in nucleus and cytoplasm.
MiRNAs, loaded into the miRISC, are thought to target multiple mRNAs, affecting the translation or stability of several target mRNAs [2, 4, 6, 26]. In this study, the miRNA-mRNA map in nucleus and cytoplasm was investigated using high-throughput sequencing of RNAs isolated by crosslinking immunoprecipitation (HITS-CLIP). The targetome of Ago2 complexes in nucleus and cytoplasm was obtained. The results showed that the main proportion of mRNAs existed in nucleus, while miRNAs function predominantly in cytoplasm. The existence of miRNAs in the Ago2 complexes of nucleus and cytoplasm revealed that miRNAs could function in nucleus and cytoplasm. Decades of studies have elucidated the function and mechanism of miRNAs in cytoplasm [1, 2, 4, 5]. In recent years, more and more investigations focus on the function and mechanism of miRNAs in nucleus. In C. elegans, it is reported that the intronic regions of pre-mRNA could be silenced by dsRNAs . In human HeLa cells, the U6 small nuclear RNA transcripts are also susceptible to RNAi-mediated post-transcriptional silencing . These findings support our miRNA-mRNA-Ago interactions in nucleus of gastric cancer cells. In nucleus, miR-320 can trigger the heterochromatinization and transcriptional gene silencing of POLR3D by binding the promoter locus of POLR3D . In addition, some splicing factors, such as U2 and U5 snRNP core subunits, SRSF1, SRSF3, SRSF7, SRSF10, PTBP1, PTBP2, KHDRBS1 and snRNAs, could be detected in immunoprecipitated Ago complex, indicating a cross-talk between splicing factors and nuclear RISC . However, the regulatory mechanism of nuclear miRNAs needs to be further investigated.
Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116(2):281–97.
Krol J, Loedige I, Filipowicz W. The widespread regulation of microRNA biogenesis, function and decay. Nat Rev Genet. 2010;11(9):597–610.
Lee Y, Ahn C, Han J, Choi H, Kim J, Yim J, Lee J, Provost P, Radmark O, Kim S, et al. The nuclear RNase III Drosha initiates microRNA processing. Nature. 2003;425(6956):415–9.
Shu L, Zhang X. Shrimp miR-12 suppresses white spot syndrome virus infection by synchronously triggering antiviral phagocytosis and apoptosis pathways. Front Immunol. 2017;8:855.
He Y, Ma T, Zhang X. The mechanism of synchronous precise regulation of two shrimp white spot syndrome virus targets by a viral microRNA. Front Immunol. 2017;8:1546.
Ren Q, Huang X, Cui Y, Sun J, Wang W, Zhang X. Two white spot syndrome virus microRNAs target the dorsal gene to promote virus infection in Marsupenaeus japonicas shrimp. J Virol. 2017;91:e02261–16.
Yang G, Yang L, Zhao Z, Wang J, Zhang X. Signature miRNAs involved in the innate immunity of invertebrates. PLoS One. 2012;7(6):e39015.
Roberts TC. The MicroRNA biology of the mammalian nucleus. Mol Ther Nucleic Acids. 2014;3:e188.
Meister G, Landthaler M, Patkaniowska A, Dorsett Y, Teng G, et al. Human Argonaute2 mediates RNA cleavage targeted by miRNAs and siRNAs. Mol Cell. 2004;15:185–97.
Khudayberdiev SA, Zampa F, Rajman M, Schratt G. A comprehensive characterization of the nuclear microRNA repertoire of post-mitotic neurons. Front Mol Neurosci. 2013;6:43.
Gagnon KT, Li L, Chu Y, Janowski BA, Corey DR. RNAi factors are present and active in human cell nuclei. Cell Rep. 2014;6:211–21.
Castel SE, Martienssen RA. RNA interference in the nucleus: roles for small RNAs in transcription, epigenetics and beyond. Nat Rev Genet. 2013;14:100–12.
Rudel S, Flatley A, Weinmann L, Kremmer E, Meister G. A multifunctional human Argonaute2-specific monoclonal antibody. RNA. 2008;14:1244–53.
Chu Y, Yue X, Younger ST, Janowski BA, Corey DR. Involvement of argonaute proteins in gene silencing and activation by RNAs complementary to a non-coding transcript at the progesterone receptor promoter. Nucleic Acids Res. 2010;38:7736–48.
Krek A, Grun D, Poy MN, Wolf R, Rosenberg L, et al. Combinatorial microRNA target predictions. Nat Genet. 2005;37:495–500.
Rajewsky N. microRNA target predictions in animals. Nat Genet. 2006;38(Suppl):S8–13.
Bentwich I. Prediction and validation of microRNAs and their targets. FEBS Lett. 2005;579:5904–10.
Cui Y, Huang T, Zhang X. RNA editing of microRNA prevents RNA-induced silencing complex recognition of target mRNA. Open Biol. 2015;5:150126.
Chi SW, Zang JB, Mele A, Darnell RB. Argonaute HITS-CLIP decodes microRNA-mRNA interaction maps. Nature. 2009;460:479–86.
Gene Ontology Consortium. The gene ontology (GO) project in 2006. Nucleic Acids Res. 2006;34:D322–6.
Leucci E, Patella F, Waage J, Holmstrom K, Lindow M, et al. microRNA-9 targets the long non-coding RNA MALAT1 for degradation in the nucleus. Sci Rep. 2013;3:2535.
Rossi JJ. A novel nuclear miRNA mediated modulation of a non-coding antisense RNA and its cognate sense coding mRNA. EMBO J. 2011;30:4340–1.
Foldes-Papp Z, Konig K, Studier H, Buckle R, Breunig HG, et al. Trafficking of mature miRNA-122 into the nucleus of live liver cells. Curr Pharm Biotechnol. 2009;10:569–78.
Castanotto D, Lingeman R, Riggs AD, Rossi JJ. CRM1 mediates nuclear-cytoplasmic shuttling of mature microRNAs. Proc Natl Acad Sci U S A. 2009;106:21655–9.
Weinmann L, Hock J, Ivacevic T, Ohrt T, Mutze J, et al. Importin 8 is a gene silencing factor that targets argonaute proteins to distinct mRNAs. Cell. 2009;136:496–507.
Gong Y, He T, Yang L, Yang G, Chen Y, Zhang X. The role of miR-100 in regulating apoptosis of breast cancer cells. Sci Rep. 2015;5:11650. https://doi.org/10.1038/srep11650.
Bosher JM, Dufourcq P, Sookhareea S, Labouesse M. RNA interference can target pre-mRNA: consequences for gene expression in a Caenorhabditis elegans operon. Genetics. 1999;153:1245–56.
Robb GB, Brown KM, Khurana J, Rana TM. Specific and potent RNAi in the nucleus of human cells. Nat Struct Mol Biol. 2005;12:133.
Kim DH, Saetrom P, Snove O Jr, Rossi JJ. MicroRNA-directed transcriptional gene silencing in mammalian cells. Proc Natl Acad Sci U S A. 2008;105:16230–5.
Ameyar-Zazoua M, Rachez C, Souidi M, Robin P, Fritsch L, Young R, et al. Argonaute proteins couple chromatin silencing to alternative splicing. Nat Struct Mol Biol. 2012;19:998–1004.
We are grateful for the funding of National Natural Science Foundation of China and National Program on the Key Basic Research Project.
This work was supported by National Natural Science Foundation of China (31430089) and National Program on the Key Basic Research Project (2015CB755903). No specific funding was received for this study.
Availability of data and materials
The data of this article are available in this manuscript.
Ethics approval and consent to participate
This study did not involve human subjects or human material. The cell lines in this study did not require ethics approval for their use.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Zhang, X., Shen, B. & Cui, Y. Ago HITS-CLIP expands microRNA-mRNA interactions in nucleus and cytoplasm of gastric cancer cells. BMC Cancer 19, 29 (2019). https://doi.org/10.1186/s12885-018-5246-0
- Ago HITS-CLIP
- miRNA-RNA interaction
- Nucleus and cytoplasm
- Gastric cancer cell