- Open Access
Selective activation of miRNAs of the primate-specific chromosome 19 miRNA cluster (C19MC) in cancer and stem cells and possible contribution to regulation of apoptosis
Journal of Biomedical Science volume 24, Article number: 20 (2017)
The human chromosome 19 miRNA cluster (C19MC) of 43 genes is a primate-specific miRNA cluster that may have biological significance in the genetic complexity of the primate. Despite previous reports on individual C19MC miRNA expression in cancer and stem cells, systematic studies on C19MC miRNA expression and biological functions are lacking.
Cluster-wide C19MC miRNA expression profiling by microarray analysis showed wholesome C19MC activation in embryonic stem cells (ESCs) and induced pluripotent stem cells (iPSCs). However, in multipotent adipose-derived mesenchymal stem cells (MSCs) and a unipotent human white pre-adipocyte cell line, only selected C19MC miRNAs were expressed. MiRNA copy number analysis also showed selective C19MC expression in cancer cells with expression patterns highly similar to those in MSCs, suggesting similar miRNA regulatory mechanisms in these cells. Selective miRNA expression also suggests complex transcriptional mechanism(s) regulating C19MC expression under specific cellular and pathological conditions. Bioinformatics analysis showed that sixteen of the C19MC miRNAs share the same “AAGUGC” seed sequence with members of the miR-302/-372 family, which are known cellular reprogramming factors. In particular, C19MC-AAGUGC-miRNAs with the nucleotides 2-7 canonical seed position as in miR-302/-372 miRNAs, may play similar roles as miR-302/-372 in induced pluripotency. A biased 3p-arm selection of the C19MC-AAGUGC-miRNAs was observed indicating that targets of the 3p species of these miRNAs may be biologically significant in regulating stemness. Furthermore, bioinformatics analysis of the putative targets of the C19MC-AAGUGC-miRNAs predicted significant involvement of signaling pathways in reprogramming, many of which contribute to promoting apoptosis by indirect activation of the pro-apoptotic proteins BAK/BAX via suppression of genes of the cell survival pathways, or by enhancing caspase-8 activation through targeting inhibitors of TRAIL-inducing apoptosis.
This work demonstrated selective C19MC expression in MSCs and cancer cells, and, through miRNA profiling and bioinformatics analysis, predicted C19MC modulation of apoptosis in induced pluripotency and tumorigenesis.
MicroRNAs (miRNAs) are short noncoding single-stranded RNAs that act post-transcriptionally as negative regulators of gene expression. In the miRNA biogenesis process, either or both 5p and 3p miRNA species may be generated from the pre-miRNA precursor arms [1, 2]. Selective maturation or co-existence of the 5p and 3p species is biologically important since they target different gene sets. MiRNAs also form family groups defined by short homologous sequences, called the seed sequence, which is located at the 5’- end of the mature miRNAs. MiRNAs of the same family often form functionally-related groups that cross-regulate targets to ensure conservation of biological functions [3, 4].
miRNAs are involved in regulating developmental patterning, maintenance of stem cell self-renewal and cancer progression [5, 6]. Specific miRNAs are up-regulated in pluripotent stem cell population but not in mature differentiated cell types in early embryonic development . Specific miRNAs have also been shown to be able to reprogram somatic cells to induced pluripotent stem cells (iPSCs) [8, 9]. MiR-302-driven cellular reprogramming coordinates stem cell division by regulating targets in the cell cycle, particularly at the G1/S restriction point . Besides the roles in stem cell biology, miRNAs also act as tumor suppressors or oncogenes in the tumorigenesis process [10,11,12,13,14]. Thus, aberrant expression of miRNA affects crucial processes in the development and progression of tumors, including induction of anti-apoptosis, development of drug resistance and cancer invasion and metastasis [15,16,17,18,19,20,21,22].
For regulatory advantages, miRNAs, particular those from the same family, are often clustered in specific chromosomal locations . One such human miRNA cluster is mapped on chromosome 19, and is called the chromosome 19 miRNA cluster, or C19MC . C19MC, one of the largest miRNA gene clusters in the human genome, contains 46 highly homologous miRNA genes, including 7 duplicated pairs of the same genes, within a ~100-kb genomic region. Importantly, C19MC is a primate-specific miRNA cluster that appeared late in the evolution of the primate lineage; bioinformatics analysis has predicted that C19MC miRNAs play critical roles in reproduction, development and differentiation in the primate compared to the lower vertebrates . Reproduction-related role of C19MC is further reflected in its restrictive expression in reproductive tissues, but not in other adult organs and tissues; C19MC expression has also been shown in pluripotent embryonic stem cells (ESCs) [7, 11, 25]. In the human placenta, C19MC is expressed en bloc from the paternal allele governed by a major promoter located 17.6 kb upstream of the first miRNA gene in the cluster . Transcription of the C19MC cluster was further suggested to be mediated by demethylation of the upstream CpG-rich master promoter region to first generate a primary transcript encompassing the entire C19MC gene cluster, followed by splicing to generate the individual precursor miRNA species, and subsequently processed by the DGCR8-Drosha microprocessor complex to generate individual mature miRNAs . The implication of a master promoter mode of transcription is the all-or-none presence of C19MC miRNAs in the expressed cells. At present, the biological functions and expression patterns of C19MC members in other stem cell types and in cancer cells have not been systematically examined in a cluster-wide manner.
In a previous study, we reported genome-wide miRNA profiling analysis of ESCs, iPSCs and mesenchymal stem cells (MSCs), proposing cross- and co-regulation by 5p and 3p paired miRNA species during reprogramming . Using the same miRNA microarray profiling dataset, we focused in this work on the expression profiles of C19MC miRNAs in various stem cell types and in cancer cells. Possible biological functions of a subset of miR-302-like C19MC miRNAs, were further investigated by bioinformatics analysis, which predicted targeting at the apoptosis pathway in the tumorigenesis of cancer cells and induced pluripotency in stem cells.
We have previously described iPSC lines derived from two MSCs, namely adipose stem cell (ASC; Invitrogen, Carlsbad, CA, USA) and human adipose-derived MSC (MSC-AT; PromoCell, Heidelberg, Germany), and from a human white pre-adipocyte (HWP) cell line [4, 27]. In this work, human adipose-derived MSC, designated ASC Lonza, was purchased from Lonza, Lonza, Verviers, Belgium. MH#1 was an iPSC cell lined established from ASC Lonza in our lab (S. Sugii, unpublished data). WJ0706 is a human MSC cell line derived from Wharton’s Jelly (WJ) obtained from Cytopeutics Sdn. Bhd, Selangor, Malaysia (http://www.cytopeutics.com). The MSC cell lines were isolated and characterized at Cytopeutics according to standard procedures and with ethical clearance . Human placenta choriocarcinoma cell line JEG-3 (ATCC HTB-36), human normal placental cell line HS 799. PI (ATCC CRL-7530) and human normal colon cell line CRL-1790 (ATCC CRL-1790) were purchased from ATCC (Manassas, VA, USA). Cancer cell lines were kindly provided by Professor Y.M. Lim, Cancer Research Center, Universiti Tunku Abdul Rahman.
miRNA microarray profiling
Total RNA was isolated from the cell lines by using the MiRNeasy Mini Kit (Qiagen, USA) according to the manufacturer’s manual. Microarray analysis was performed using the SmartChip Human MicroRNA Panel version 3.0 (WaferGen Biosystems, Fremont, CA, USA) containing 1036 unique real-time PCR reactions in quadruplicates as previously described . To identify differentially expressed miRNAs, the iPSC data were compared with data of the parental MSC or HWP cells from which they were derived. For calculation of expression levels, the All-Mean Normalization method was employed, where mean Ct values of all expressed genes were used. To compute the expression levels of expressed miRNAs, the Ct values of each sample were compared to its average Ct (All-Mean) to obtain the ΔCt values. ΔΔCt was then calculated by the two ΔCt values between the iPSC and its parental cells. Log2(fold change), or log2(FC), was computed as log2[FC (2-ΔΔCt)]. The selection criteria for differentially expressed miRNA was log2(FC) > 1.5 or < -1.5, with p < 0.05 as determined by the Student’s t test.
Reverse transcription and miRNA real-time PCR assays
RNA was reverse-transcribed using the TaqMan MicroRNA Reverse Transcription Kit and miRNA-specific stem-loop primers (Applied Biosystems, USA). The reverse transcription reaction was as follows: 0.075 μl of 100 mM dNTP, 0.5 μl RT enzyme (50 U/μl), 0.75 μl 10X RT buffer, 0.094 μl RNAse inhibitor (20 U/μl), 50 ng RNA. Reaction conditions were 16 °C for 30 min, 42 °C for 30 min, 85 °C for 5 min, and held at 4 °C. Real-time PCR assays of the transcribed cDNA were performed using the TaqMan MicroRNA assays (Applied Biosystems, USA). The reaction mixture was as follows: 10 μl TaqMan 2X universal PCR master mix, 1 μl 20X TaqMan MicroRNA Assay, 1.33 μl cDNA and 7.67 μl nuclease-free water. Reaction conditions were: 95 °C for 10 min, followed by 40 cycles at 95 °C for 15 s, and 60 °C for 1 min.
Determination of absolute copy number of mature miRNAs
Synthetic mature miRNAs (Integrated DNA Technologies IDT, Coralville, IA, USA) were serially diluted to final concentrations of 200 nM, 20 nM, 2 nM, 0.2 nM, 0.02 nM, 2 pM, 0.2 pM and 0.02 pM. Serially-diluted synthetic RNAs were reverse-transcribed and subjected to real-time PCR analysis concurrently with the sample RNAs. Standard curves were included on each plate of the miRNA TaqMan assays to convert the cycle threshold (Ct) values of each sample into the corresponding number of miRNA copies in each cell, assuming that each cell contains 15 pg total RNA, as previously described . Ct values ≥ 35 indicated that their expression levels were too low for accurate analysis, and were considered no detectable expression. The cut-off threshold of miRNA expression was, therefore, standardized at Ct < 35.
Construction of phylogeny tree
The stem-loop sequences for C19MC were downloaded from miRBase ver. 21. A phylogenetic tree was generated by multiple sequence alignment using the Clustal method of the Megalign project provided by DNAstar® (Salt Lake City, Utah, USA).
Prediction of miRNA target genes
miRNA:mRNA interactions were predicted using the major miRNA databases TargetScan and microRNA.org. To identify genes and pathways targeted specifically by selected C19MC-AAGUGC-miRNAs, overlapping target gene sets of the selected miRNAs were used for the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and Gene Ontology (GO) annotation analysis based on the web-based DAVID (Database for Annotation, Visualization and Integrated) algorithm. The criteria of analysis was EASE score ≤ 0.05, in which EASE score is a modified Fisher Exact P value in the DAVID system used for gene-enrichment analysis. An EASE score P value = 0 represents perfect enrichment; P value ≤ 0.05 was considered as significant gene-enrichment in a specific annotation category.
Selective activation of C19MC miRNAs in mesenchymal stem cells
We have previously reported genome-wide miRNA expression profiling of two ESCs, two multipotent adipose-derived MSCs, a unipotent HWP cell line, and the three iPSC lines derived from the MSCs and HWP . In this study, the expression data of C19MC miRNAs in the MSCs and HWP, and the MSC- and HWP-derived iPSCs were extracted for further analysis. Our results showed that all the forty-five C19MC miRNAs included in the microarray in either the 5p or 3p or in both 5p/3p configurations were expressed, albeit to different extents, in all the three pluripotent iPSC cell lines tested (Table 1), and in the hESC controls (data not shown). The miR-372 family that lies adjacent to the C19MC cluster (Fig. 1a) was also included in the analysis since they have been reported to be expressed in pluripotent stem cells . Of the forty-five C19MC miRNAs, thirty-nine were significantly (p < 0.05) expressed, as previously reported . Expression of the C19MC miRNAs in the iPSCs was generally two-fold or greater than that in the parental cell lines; the highest level of expression was 8.375 log2(fold change) in miR-520b (Table 1). Notably, both the 5p and 3p miRNA species were expressed in most cases; otherwise, the 3p species was the favored precursor arm selected for the mature miRNAs, as opposed to frequent 5p arm expression in most other miRNA genes .
On the other hand, only selected C19MC miRNAs were found to be expressed in MSC and HWP (Tables 1 and 2). Many of the expressed miRNAs share the “AAGUGC” seed sequence of the known reprogramming miR-302 miRNA family; these miRNAs are called the C19MC-AAGUGC-miRNAs in this work (see Fig. 2a and depiction below). Twenty-two (48.9%) of the forty-five C19MC miRNAs were activated in one or both MSC cell lines. Only eight miRNAs were expressed in HWP, which were, interestingly, also all expressed in the two MSC and all pluripotent cells (Table 2). This may suggest that these eight miRNAs constitute the minimal miRNA set require for minimal potency in the unipotent HWP. Thus, the cluster-wide microarray results indicated selective activation of twenty-two C19MC miRNAs in multipotent mesenchymal stem cells.
It has been reported that C19MC miRNAs are not expressed in adult tissues except in tissues of the reproductive system . To obtain further supporting evidences on selective activation, expression of eight miRNAs spanning the C19MC cluster (Fig. 1a), but with different genomic structures, was selected for further experimentally verification; amongst the selected miRNAs, miR-512-3p is transcribed by the two miR-512-1 and-512-2 genes located at the 5’-end of the C19MC miRNA gene cluster; miR-520c-3p, -519b-3p and -520f-3p are single miRNA genes located between previously proposed exons; miR-524-5p and -517a-3p are two of three miRNA genes mapped on intron 18 and miR520d-5p and -520g-3p are two of four miRNAs mapped on intron 20 (Fig. 1a) . Verification was done in three other different MSC cell lines, namely the MSC cell line WJ0706 derived from the Wharton’s Jelly , and two other adipose-derived MSC cell lines, ASC-Inv and ASC Lonza (Fig. 1b). In the experiments, two other adipose MSC-derived iPSCs, ASC-iPSC and MH#1, and two hESCs, H6 and H9, were included. The miRNA expression levels were determined as the absolute miRNA transcript copy number per cell, which ranged from 0 copy, at a real-time RT-PCR Ct value ≥ 35 (see Methods), to 377,200 copies per cell at a Ct value of 25.7 in miR-520g-3p in ESC H9 cells (Fig. 1b).
Consistent with the miRNA microarray results, the selected miRNAs were all expressed to different levels in all four iPSCs and ESCs (Fig. 1b). In contrast, the tested miRNAs were either not expressed, or expressed to different but lower levels in the MSCs tested. MSC expression of the eight C19MC miRNAs could be grouped in three expression patterns: group A, which included miR-512-3p and -520c-3p, showed very low or undetectable expression in the MSCs; expression of the group B miR-520d-5p, 519b-3p and -524-5p was detected in at least one or both MSC cell lines, whereas miR-520f-3p, -517a-3p and -520g-3p in group C were all expressed all three MSC cell lines (Fig. 1b and Table 2). The collective results obtained from the microarray and real-time RT-PCR experiments, therefore, confirmed selective C19MC expression in multipotent MSCs, and en bloc expression in pluripotent iPSCs. Furthermore, there seemed to be no correlation between the expression pattern and the physical location of the miRNA genes tested (Fig. 1a). Notably, the miR-524-5p and -517a-3p and the miR520d-5p and -520g-3p couples are flanked by two proposed exons but belong to different expression groups B and C (Fig. 1a). The data suggest regulation by different promoters or transcriptional regulatory mechanism(s) other than simple splicing of the two flanking exon and co-processing of the spliced intron sequence as previously proposed for C19MC expression in a choriocarcinoma JEG-3 cell line . The observation further suggests a critical biological role of the expressed C19MC miRNA in conferring different degrees of stemness to the stem cells, particularly in MSCs.
Selective activation of C19MC miRNAs in cancer cells
Previous reports have indicated frequent activation of C19MC miRNAs in different cancer types, including colorectal cancer, breast cancer and primitive neuroectodermal brain tumor [13, 30, 31] (see below). To investigate C19MC expression in cancer cells, the expression of the same set of eight C19MC miRNAs was also quantified as gene copy number per cell in two colorectal cancer (HCT15 and SKCO1), one breast cancer (MCF-7) and one hepatocellular carcinoma (HepG2) cell lines; the choriocarcinoma (JEG-3) cell line, which was derived from the reproductive system, was included a positive control since JEG-3 cells have been shown to express all C19MC miRNAs in high levels  (Fig. 1c). Two cell lines CRL-1790 and HS799. PI, derived from normal colon and placenta tissues, respectively, were also included in the analysis. Despite en bloc and high-level C19MC expression in JEG-3 cells, only four of the eight miRNAs, namely miR-520d-5p of Group B as defined above for stem cell expression, and all three Group C miRNAs, miR-520f-3p, -517a-3p and -520g-3p, were shown to be expressed in the normal placenta cell line Hs799. PI. Furthermore, expression of the Group B miR-524-5p, and all three Group C miRNAs was detected in CRL-1790, which was derived from normal fetal colon epithelium (Fig. 1c). The observed expression of selective C19MC in fetal colon epithelium and in the placenta is consistent with previous conclusions that C19MC is specifically expressed in reproduction and developmental process-related tissues and is silenced in normal tissues [7, 11, 25]. Interestingly, in the five cancer cell lines examined, the selective expression patterns of the eight miRNAs was similar those shown in MSCs above (Fig. 1b). Group A miRNAs also showed very low or undetectable expression in normal and cancer cells, except in JEG-3, whereas the Group B miRNAs were detected in one or more cancer cell lines; all three Group C miRNAs were expressed all four cancer cell lines (Fig. 1c). Taken together, quantitative expression analysis showed highly similar C19MC miRNA expression profiles found in MSCs and cancer cells, suggesting that the C19MC miRNAs may share some similar molecular and biological features in transcriptional regulation and in the etiological pathways in acquiring multipotency and cancer phenotype.
Identification of C19MC miRNAs harboring the “AAGUGC” seed sequence
miRNA-mRNA interactions involve the seed region at the 5’ end of the miRNA; hence, seed sequences are important predictors for the identification of miRNA-targeted transcripts . MiRNAs that share a common seed sequence also might share target specificity and possibly biological functions. On sequence alignment, sixteen C19MC miRNAs were found to share the same seed sequence, 5’-AAGUGC-3’, with the reported reprogramming-able miR-302 and miR-372 miRNA families [8, 9] (Fig. 2a). These miRNAs are designated as “C19MC-AAGUGC-miRNAs”. Furthermore, it is noted that the AAGUGC seed position at 5’ end is variable among the C19MC-AAGUGC-miRNAs: subgroup I miRNAs, which includes eight miR-519 and -520 subfamilies, have the seed sequence located at the canonical and optimal 5’-nucleotide positions (nts) 2-7, as in the miR-302/-372 families; the seed sequence of the four subgroup IIa miRNAs is at location nts 1-6, and that of the remaining subgroup IIb miRNAs is at nts 3-8 and 4-9 (Fig. 2a). Hence, despite the presence of the AAGUGC seed sequence, it is more likely that the nts 2-7 canonical subgroup of the C19MC-AAGUGC-miRNAs may target genes that share similar functions as the miR-302/-372 miRNAs.
While the 5p arm of a pre-miRNA precursor is normally selected for maturation , it is noted that the C19MC-AAGUGC-miRNAs are predominantly derived from the 3p arm of the precursor miRNAs, hinting at an evolutionary bias in 3p selection with possible biological implications. Further supporting evidence of conservation of the C19MC-AAGUGC-miRNAs was derived from the construction of a phylogenetic tree of all precursor sequences of the C19MC miRNAs (Fig. 2b). Most C19MC-AAGUGC-miRNAs are grouped into the same cluster in the top half of the phylogenetic tree. Four of the remaining C19MC-AAGUGC-miRNAs form another cluster in the middle of the tree and the remaining two miRNAs are scattered in different branches in the lower half of the tree (Fig. 2b).
Bioinformatics predictions of possible biological functions of group I C19MC-AAGUGC-miRNAs
It is noted that the C19MC-AAGUGC-miRNAs with the canonical nts 2-7 seed position, defined here as Group I (Fig. 2), contributed more significantly in gene targeting. Thus, in this study, we focused on analysis of potential biological functions of C19MC-AAGUGC-miRNAs in group I. Bioinformatics searches showed a total of 2058 putative target genes targeted by group I C19MC-AAGUGC-miRNAs (Fig. 3a and Additional file 1: Table S1). However, construction of a Venn diagram showed that only 262 putative target genes are common between the miR-519 and miR-520 subfamilies in group I, indicating that the miR-519 and -520 subfamilies target different sets of genes. The overlapping gene sets among miR-302/372 and the miR-519 and miR-520 subfamilies in group I were further compared (Fig. 3a). The results showed that 1185 putative shared genes were obtained between the miR-520 and -302/372 families (Fig. 3a, blue box and Additional file 1: Table S1), suggesting that the miR-520 subfamily might share similar biological functions with the miR-302/372 family. The group I miR-519 subfamily also shares 262 putative target genes with the miR-302/-372 families, far fewer than the miR-520 subfamily (Fig. 3a, red box). Consistent with the bioinformatics prediction, a literature review showed that a number of validated targets have indeed been reported to be shared between the miR-302/372 and the group I C19MC-AAGCGU-miRNA families (Table 3).
The 2058 putative target genes were further subjected to GO analysis and KEGG pathway annotation (Fig. 3b-d). Of the 828 predicted targets in the top 10 GO terms in biological functions, 616 (74.4%) putative genes are associated with transcriptional and translational regulation of gene expression (Fig. 3b, GO terms 1, 2, 4, 5, 7 & 8). The remaining predicted targets regulate apoptosis, nervous system development, cellular response to DNA damage stimulus and cell cycle. The majority of the 2058 predicted genes in GO terms in molecular functions is likewise associated with transcriptional and translational regulation in some way (Fig. 3c), and in epigenetic regulation (Fig. 3c, GO term 9). Four hundred eleven genes (20.0%) are related to metal or zinc ion binding (Fig. 3c, GO terms 2 & 4). which may also be components of signaling pathways. Taken together, the GO analysis data suggested that the group I C19MC-AAGUGC-miRNAs are mainly associated with the regulation of gene expression, cell proliferation and apoptosis via various signaling pathways.
The regulatory pathways were further annotated by interrogation of the KEGG database, which yielded 24 pathways, which included 568 genes in total (Additional file 2: Table S2); 14 of the 24 KEGG pathways which may be related to pluripotency and cancer are shown in Fig. 3d. Ten of the 24 pathways, which included 260 (45.8%) genes, are different signaling pathways that are known to be involved in the growth and development processes [9, 10, 12, 21, 32,33,34,35,36]. Notably, 129 (22.7%) genes are associated with pathways regulating apoptosis including PI3K-AKT, MAPK, HIF-1 and TNF (Fig. 3d; see also Fig. 4 and Discussion below). The highest-enriched PI3K-Akt signaling pathway (56 genes) regulates cell survival by reducing apoptosis, stimulating cell growth and increasing proliferation . Furthermore, many of the genes are related to pathways that regulate the cell cycle (22 genes) and apoptosis (13 genes) (Fig. 3d), which are important cellular events in the initiation and maintenance of stem cell pluripotency and tumorigenesis.
Possible group I C19MC-AAGUGC-miRNAs targeting of the pro-apoptosis functions in the survival pathway
Suppression of apoptosis is an important feature of the initiate phase of the reprogramming process . On the other hand, apoptosis dysregulation is associated with the different stages of tumorigenesis, including initiation, progression and metastasis . A database search showed that the group I C19MC-AAGUGC-miRNAs target 179 apoptosis-associated genes (Additional file 3: Table S3). On the other hand, the KEGG pathway analysis above (Fig. 3d) has also revealed that the highest number of putative target genes of group I miRNAs are associated with PI3K-Akt, a survival pathway. Hence, we hypothesized that the group I miRNAs acted more specifically to inhibit apoptosis by targeting survival-related genes. Fifteen survival-related genes were predicted targets of the group I miRNAs (Table 4). Two out of the fifteen genes, viz. NIK and RelA, have been experimentally validated as direct targets miR-520e and miR-520c-3p [10, 12]. Importantly, the group I miRNAs may promote apoptosis either by indirectly activating pro-apoptotic proteins BAK/BAX through suppression of the cell survival-related genes , or by enhancing caspase-8 activation through targeting inhibitors of TRAIL-inducing apoptosis [40, 41]. Taken together, the group I C19MC-AAGUGC-miRNAs were predicted by bioinformatics analysis to regulate apoptosis, which is important in the initial phase of cellular reprogramming, and in particular the cell survival pathways, which are directly relevant to tumorigenesis processes.
Selective C19MC miRNA expression in MSC and in cancer cells suggests a complex transcriptional regulatory mechanism
In the present study and in the literature, data showed similar and disperse expression patterns of eight tested C19MC miRNAs in both mesenchymal stem and cancer cells (Fig. 1), in contrast to the previous model of en bloc expression in the choriocarcinoma JEG-3 cell line regulated by a master promoter . Another studies have shown that the highly abundant Alu repetitive sequences embedded within the C19MC genomic region may function as independent RNA polymerase II promoters [42, 43]. Our study clearly showed selective C19MC miRNA activation in MSCs and HWP, and in cancer cells, suggesting that C19MC transcripts are more likely regulated by multiple promoters, which may in turn be active by condition-specific transcription factors. Furthermore in cancer cells, chromosomal rearrangements, amplification and modification of the promoter(s) or specific transcription factors could further regulate the selective C19MC miRNA expression. Previous reports have, indeed, shown that translocation of chromosomal band 19q13.4 selectively activated C19MC miRNAs in thyroid adenomas, and that C19MC genomic amplifications in an aggressive primitive neuroectodermal brain tumors were associated with specific and abundant expression of miR-517c and -520g . Moreover, epigenetic alterations in the C19MC genomic region may also play important role in regulating C19MC expression, particularly in cancer and the dynamic stem cells. Promoter silencing of C19MC miRNAs by the DNA methylation inhibitor, 5-azacytidine, activated sixteen C19MC miRNAs . Furthermore, placenta-derived mesenchymal stem cells were reported to escape epigenetic silencing of the paternal allele resulting in a number of the C19MC miRNAs being abundantly expressed . Specific activation of the C19MC miR-512-5p by histone deacetylase inhibitors was also reported in human gastric cancer cells . Transcription factors acting in trans are essential regulators of C19MC miRNA expression as shown by direct binding of p53 and the estrogen receptor α (ERα) to presumptive promoters of C19MC miR-519d and miR-515-5p, respectively, in chromatin immunoprecipitation assays [47, 48]. As a result of the specific DNA binding, miR-519d is up-regulated by p53, whereas ERα mediates both down- and up-regulated expression of miR-515-5p induced by estrogens and tamoxifen, respectively. Thus, transcription of specific C19MC miRNAs in MSC and cancer cells is likely highly complex, and may be dependent on the cellular and pathological state of the cells.
It was previously reported that C19MC is silenced in normal tissues [11, 25] due to hypermethylation of both the paternal and maternal alleles . However, placenta is able to escape epigenetic silencing by maintaining paternal allele-expression [26, 45]. Moreover, the expression of miR-498, a member of C19MC, was reported in the fetal brain , echoing our report of Group C miRNAs being expressed in a fetal colon epithelium-derived cell line, CRL-1790 and placental Hs799. PI (Fig. 1c), consistent with C19MC expression in reproductive and developmental process-related tissues, relevant to the primate-specificity of the C19MC cluster.
In this study, we found that the 3p arms of the C19MC miRNA precursors were predominantly selected in ESCs and iPSCs (Table 2 and Fig. 2). Several studies have previously demonstrated that preferred arm selection is temporal- and spatial-dependent [50, 51]. Indeed, the 3p miRNA species have been shown to be more abundantly expressed in tumor tissues as opposed to preferred 5p selection in normal tissues . Echoing these findings, the miR-302-like C19MC are also predominantly 3p-biased, possibly targeting genes which are biologically significant in regulating the stemness of stem cells and the tumor phenotype in cancers.
Structural and function significance of the group I C19MC-AAGUGC-miRNAs
Our results also showed divergence in the positions of the AAGUGC seed sequence among C19MC miRNAs carrying the hexameric sequence. The canonical seed region situated at nucleotides 2-7 is a perfect seed match which markedly decreases the presence of false-positive bioinformatics predictions, thus improving prediction reliability . Furthermore, the canonical seed region is crucial and sufficient to trigger target silencing . The hexamer of half of the C19MC-AAGUGC-miRNAs reported here are located at nts 2-7, designated as group I C19MC-AAGUGC-miRNAs in this report (Fig. 2), suggesting high possibility that the predicted genes are the putative targets. Other non-canonical C19MC-AAGUGC-miRNAs are likely to have lower affinity and specificity and may be limited in mediating repression without the 3’-compensatory binding .
The group I miRNAs are composed of the miR-519 and -520 subfamilies. Despite their similar seed location at nts 2-7, Venn diagram analysis shows that these two subfamilies share only a small number of putative target genes (Fig. 3a). Common prediction algorithms that use identical powerful prediction characteristics, such as the mandatory stringent seed base-pairing produce different prediction results properly due to usage of various UTR databases as well as different internal criteria . In this study, the putative target sets of the miR-519 and -520 subfamilies are overlapping gene sets predicted by two different prediction algorithms. Furthermore, one of the characteristics of target prediction, the sequence context surrounding the seed binding site of the target transcript , between the miR-519 and -520 subfamilies are also dissimilar (Fig. 2a). This may explain the different target gene sets of these two subfamilies. It has been reported that miR-96 and -182 that have identical seed region (UUGGCA, nucleotides 2-7) regulate different targets . However, the miR-520 and miR-302/372 families share a significant number of target genes (Fig. 3a) suggesting common biological functions. Hence, it is highly likely that the group I miR-520 miRNAs may also contribute to reprogramming, as supported by the predicted involvement of miR-520 miRNAs in the reprogramming-related apoptosis and cell proliferation pathways (see Fig. 4 and Discussion below).
Regulation of C19MC miRNAs in tumorigenesis and stemness
Selective activation of C19MC miRNAs in MSCs and cancer cells reported here suggests functional involvement of the activated miRNAs in maintaining the stemness and promoting cancer development. Frequent aberrant C19MC miRNA expression in cancers has been reported [20, 44, 47]. Activation of the C19MC miR-519d was shown to target CDKN1A/p21, PTEN, AKT3 and TIMP2, and is closely associated with the pathogenesis of hepatocellular carcinoma by promoting cell proliferation and invasion, and in inhibiting apoptosis . In breast cancer, high expression levels of plasma miR-520g is correlated with patients with lymph node metastasis and mammary gland invasion, and suppressed p53 expression .
On the other hands, C19MC miRNAs have also been shown to play important role in cellular stemness state. In normal embryonic development, many C19MC miRNAs have been shown to be expressed only in undifferentiated or germinal tissues, and C19MC expression inhibits differentiation of human embryonic stem cells [7, 26, 55]. The observation that the cellular reprogramming-able transcription factors OCT4 and NANOG regulate C19MC miRNA expression in human embryonic stem cells further supports close association of C19MC with induced pluripotency . Moreover, the identification of sixteen miR-302-like C19MC miRNAs predicts functions in promoting “stemness” as the miR-302 and miR-372 families. Similarly, eight miR-302-like C19MC miRNAs were previously shown to promote cell proliferation and cell-cycle progression by targeting p21, an inhibitor of the G1/S transition, as for the miR-302 and -372 families [5, 6].
Possible involvement of group I C19MC-AAGUGC-miRNAs in regulating the apoptosis pathway common to stemness and cancer phenotype
Suppressed apoptosis is important to both the initial phase of acquiring pluripotency and in cancer progression [57, 58]. A combined expression profile and bioinformatics analysis reported in this work has, indeed, shown that the group I C19MC-AAGUGC-miRNAs, target genes related to the survival pathways (Table 4). Based on the predicted target genes, a scheme that correlates the group I C19MC-AAGUGC-miRNAs to stemness and cancer phenotype is proposed (Fig. 4). In general, group I miRNAs may enhance apoptosis through the PIK3/ATK, TNFs/NF-κB and TRAIL pathways, as predicted by KEGG pathway analysis (Fig. 4) [40,59,, 58–60]. The PIK3 pathway is activated by a wide range of extracellular signals, including cytokines, e.g. IL-2 , growth factors, e.g. IGF1  and components of the extracellular matrix (ECM) such as TSP-1 , all of which are the predicted targets of the group I miRNAs (Fig. 4). It is proposed here that the miRNAs target and inactivate the PIK3/AKT3 pathway by inhibition of the PIK3-related upstream genes TSP-1, IL-2, IGF1, KIT, SOS1 and PIK3CA, and the downstream AKT1 gene. The second important mechanism of cell survival is tumor necrosis factors (TNFs) activation of anti-apoptotic proteins via the nuclear factors of kappa B (NF-κB) signaling cascade (Fig. 4). Similar to the PIK3/ATK pathway, group I C19MC-AAGUGC- miRNAs may enhance apoptosis by the predicted targeting of the TNFα, TLR4, TRAF6, TAK1, NIK, MALT1 and RelA genes. Thirdly, group I miRNAs are also predicted to silence genes, such as DcR2, that are inhibitory to the TRAIL-induced apoptosis pathway, resulting in pro-apoptosis . The group I miRNAs-modulated pathways subsequently suppress the activation of downstream effector caspase-3, -6, and -7, thus inhibiting apoptosis and promoting proliferation .
In the present study, the data show selective expression of C19MC miRNAs in cancer and stem cells, offering insights into possible involvement of selective C19MC miRNAs in regulation of “stemness” and tumorigenesis, possibly via the cell survival pathways. More specifically, a subgroup of sixteen C19MC miRNAs has been identified that shares the same AAGUGC seed sequence as the reprogramming miR-302/372 family, predicting contribution of the C19MC-AAGUGC-miRNAs to the reprograming process. Further elucidation of the biological functions of C19MC miRNAs, particular the miRNA-302-like subclass, may lead to potential applications in more efficient cellular reprogramming and in cancer therapy.
The human chromosomal 19 miRNA cluster
Estrogen receptor α
Embryonic stem cells
- HCT15 and SKCO1:
Human white pre-adipocyte
Induced pluripotent stem cells
Human placental choriocarcinoma cell line
Kyoto Encyclopedia of Genes and Genomes
Mesenchymal stem cells
Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136:215–33.
Meijer HA, Smith EM, Bushell M. Regulation of miRNA strand selection: follow the leader? Biochem Soc Trans. 2014;42:1135–40.
Choo KB, Soon YL, Nguyen PNN, Hiew MSY, Huang C-J. MicroRNA-5p and -3p co-expression and cross-targeting in colon cancer cells. J Biomed Sci. 2014;21:95.
Huang CJ, Nhi Nguyen PN, Choo KB, Sugii S, Wee K, Cheong SK, et al. Frequent co-expression of miRNA-5p and -3p species and cross-targeting in induced pluripotent stem cells. Int J Med Sci. 2014;11:824–33.
Wang Y, Baskerville S, Shenoy A, Babiarz JE, Baehner L, Blelloch R. Embryonic stem cell-specific microRNAs regulate the G1-S transition and promote rapid proliferation. Nat Genet. 2008;40:1478–83.
Wu S, Huang S, Ding J, Zhao Y, Liang L, Liu T, et al. Multiple microRNAs modulate p21Cip1/Waf1 expression by directly targeting its 3’ untranslated region. Oncogene. 2010;29:2302–8.
Razak SRA, Ueno K, Takayama N, Nariai N, Nagasaki M, Saito R, et al. Profiling of microRNA in human and mouse ES and iPS cells reveals overlapping but distinct microRNA expression patterns. PLoS One. 2013;8:e73532.
Anokye-Danso F, Trivedi CM, Juhr D, Gupta M, Cui Z, Tian Y, et al. Highly efficient miRNA-mediated reprogramming of mouse and human somatic cells to pluripotency. Cell Stem Cell. 2011;8:376–88.
Subramanyam D, Lamouille S, Judson RL, Liu JY, Bucay N, Derynck R, et al. Multiple targets of miR-302 and miR-372 promote reprogramming of human fibroblasts to induced pluripotent stem cells. Nat Biotechnol. 2011;29:443–8.
Zhang S, Shan C, Kong G, Du Y, Ye L, Zhang X. MicroRNA-520e suppresses growth of hepatoma cells by targeting the NF-κB-inducing kinase (NIK). Oncogene. 2012;31:3607–20.
Liang Y, Ridzon D, Wong L, Chen C, Du T, Zamore P, et al. Characterization of microRNA expression profiles in normal human tissues. BMC Genomics. 2007;8:166.
Keklikoglou I, Koerner C, Schmidt C, Zhang JD, Heckmann D, Shavinskaya A, et al. MicroRNA-520/373 family functions as a tumor suppressor in estrogen receptor negative breast cancer by targeting NF-κB and TGF-β signaling pathways. Oncogene. 2012;31:4150–63.
Kleinman CL, Gerges N, Papillon-Cavanagh S, Sin-Chan P, Pramatarova A, Quang D-AK, et al. Fusion of TTYH1 with the C19MC microRNA cluster drives expression of a brain-specific DNMT3B isoform in the embryonal brain tumor ETMR. Nat Genet. 2014;46:39–44.
Lu S, Zhu Q, Zhang Y, Song W, Wilson MJ, Liu P. Dual-Functions of miR-373 and miR-520c by Differently Regulating the Activities of MMP2 and MMP9. J Cell Physiol. 2015;230:1862–70.
Gui S, Chen X, Zhang M, Zhao F, Wan Y, Wang L, et al. Mir-302c mediates influenza A virus-induced IFNβ expression by targeting NF-κB inducing kinase. FEBS Lett. 2015;589:4112–8.
Cui W, Zhang Y, Hu N, Shan C, Zhang S, Zhang W, et al. miRNA-520b and miR-520e sensitize breast cancer cells to complement attack via directly targeting 3’UTR of CD46. Cancer Biol Ther. 2010;10:232–41.
Min D, Lv X, Wang X, Zhang B, Meng W, Yu F, et al. Downregulation of miR-302c and miR-520c by 1,25(OH)2D3 treatment enhances the susceptibility of tumour cells to natural killer cell-mediated cytotoxicity. Br J Cancer. 2013;109:723–30.
Su C-M, Wang M-Y, Hong C-C, Chen H-A, Su Y-H, Wu C-H, et al. miR-520h is crucial for DAPK2 regulation and breast cancer progression. Oncogene. 2016;35(9):1134–42.
Zhou J-Y, Zheng S-R, Liu J, Shi R, Yu H-L, Wei M. MiR-519d facilitates the progression and metastasis of cervical cancer through direct targeting Smad7. Cancer Cell Int. 2016;16:21.
Kan H, Guo W, Huang Y, Liu D. MicroRNA-520 g induces epithelial-mesenchymal transition and promotes metastasis of hepatocellular carcinoma by targeting SMAD7. FEBS Lett. 2015;589:102–9.
Liu P, Wilson MJ. miR-520c and miR-373 upregulate MMP9 expression by targeting mTOR and SIRT1, and activate the Ras/Raf/MEK/Erk signaling pathway and NF-κB factor in human fibrosarcoma cells. J Cell Physiol. 2012;227:867–76.
Zhang J, Liu L, Sun Y, Xiang J, Zhou D, Wang L, et al. MicroRNA-520 g promotes epithelial ovarian cancer progression and chemoresistance via DAPK2 repression. Oncotarget. 2016;7:26516–34.
Wang J, Haubrock M, Cao K-M, Hua X, Zhang C-Y, Wingender E, et al. Regulatory coordination of clustered microRNAs based on microRNA-transcription factor regulatory network. BMC Syst Biol. 2011;5:199.
Bortolin-Cavaille ML, Dance M, Weber M, Cavaille J. C19MC microRNAs are processed from introns of large Pol-II, non-protein-coding transcripts. Nucleic Acids Res. 2009;37:3464–73.
Lin S, Cheung WKC, Chen S, Lu G, Wang Z, Xie D, et al. Computational identification and characterization of primate-specific microRNAs in human genome. Comput Biol Chem. 2010;34:232–41.
Noguer-dance M, Abu-amero S, Al-khtib M, Lefe A, Coullin P, Moore GE. The primate-specific microRNA gene cluster (C19MC) is imprinted in the placenta. Hum Mol Genet. 2010;19:3566–82.
Sugii S, Kida Y, Berggren WT, Evans RM. Feeder-dependent and feeder-independent iPS cell derivation from human and mouse adipose stem cells. Nat Protoc. 2011;6:346–58.
Choo KB, Tai L, Hymavathee KS, Wong CY, Nguyen PNN, Huang CJ, et al. Oxidative stress-induced premature senescence in Wharton’s jelly-derived mesenchymal stem cells. Int J Med Sci. 2014;11:1201–7.
Chen C, Ridzon DA, Broomer AJ, Zhou Z, Lee DH, Nguyen JT, et al. Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res. 2005;33:e179.
Ma W, Yu Q, Jiang J, Du X, Huang L, Zhao L, et al. miR-517a is an independent prognostic marker and contributes to cell migration and invasion in human colorectal cancer. Oncol Lett. 2016;11:2583–9.
Ren G, Wang L. Study on the relationship between miR-520 g and the development of breast cancer. Eur Rev Med Pharmacol Sci. 2016;20:657–63.
Zhang X, Tang N, Hadden TJ, Rishi AK. Akt, FoxO and regulation of apoptosis. Biochim Biophys Acta. 1813;2011:1978–86.
Tsai C-H, Tsai H-C, Huang H-N, Hung C-H, Hsu C-J, Fong Y-C, et al. Resistin promotes tumor metastasis by down-regulation of miR-519d through the AMPK/p38 signaling pathway in human chondrosarcoma cells. Oncotarget. 2014;6:258–70.
Brownlie RJ, Zamoyska R. T cell receptor signalling networks: branched, diversified and bounded. Nat Rev Immunol. 2013;13:257–69.
Kurosaki T, Shinohara H, Baba Y. B cell signaling and fate decision. Annu Rev Immunol. 2010;28:21–55.
Danielsen SA, Eide PW, Nesbakken A, Guren T, Leithe E, Lothe RA. Portrait of the PI3K/AKT pathway in colorectal cancer. Biochim Biophys Acta. 1855;2015:104–21.
Smith ZD, Nachman I, Regev A, Meissner A. Dynamic single-cell imaging of direct reprogramming reveals an early specifying event. Nat Biotechnol. 2010;28:521–6.
Fulda S. Tumor resistance to apoptosis. Int J Cancer. 2009;124:511–5.
Westphal D, Dewson G, Czabotar PE, Kluck RM. Molecular biology of Bax and Bak activation and action. Biochim Biophys Acta. 1813;2011:521–31.
Wang S, El-Deiry WS. TRAIL and apoptosis induction by TNF-family death receptors. Oncogene. 2003;22:8628–33.
Garofalo M, Croce CM. MicroRNAs as therapeutic targets in chemoresistance. Drug Resist Updat. 2013;16:47–59.
Borchert GM, Lanier W, Davidson BL. RNA polymerase III transcribes human microRNAs. Nat Struct Mol Biol. 2006;13:1097–101.
Saito Y, Suzuki H, Tsugawa H, Nakagawa I, Matsuzaki J, Kanai Y, et al. Chromatin remodeling at Alu repeats by epigenetic treatment activates silenced microRNA-512-5p with downregulation of Mcl-1 in human gastric cancer cells. Oncogene. 2009;28:2738–44.
Li M, Lee KF, Lu Y, Clarke I, Shih D, Eberhart C, et al. Frequent amplification of a chr19q13.41 microRNA polycistron in aggressive primitive neuroectodermal brain tumors. Cancer Cell. 2009;16:533–46.
Tsai K, Kao H, Chen H, Chen S, Lin W. Epigenetic control of the expression of a primate-specific microRNA cluster in human cancer cells. Epigenetics. 2009;4:587–92.
Flor I, Neumann A, Freter C, Helmke BM, Langenbuch M, Rippe V, et al. Abundant expression and hemimethylation of C19MC in cell cultures from placenta-derived stromal cells. Biochem Biophys Res Commun. 2012;422:411–6.
Fornari F, Milazzo M, Chieco P, Negrini M, Marasco E, Capranico G, et al. In hepatocellular carcinoma miR-519d is up-regulated by p53 and DNA hypomethylation and targets CDKN1A/p21, PTEN, AKT3 and TIMP2. J Pathol. 2012;227:275–85.
Pinho FG, Frampton AE, Nunes J. Downregulation of microRNA-515-5p by the estrogen receptor modulates sphingosine kinase 1 and breast cancer cell proliferation. Cancer Res. 2013;73:5936–48.
Zhang R, Wang YQ, Su B. Molecular evolution of a primate-specific microRNA family. Mol Biol Evol. 2008;25:1493–502.
Jagadeeswaran G, Zheng Y, Sumathipala N, Jiang H, Arrese EL, Soulages JL, et al. Deep sequencing of small RNA libraries reveals dynamic regulation of conserved and novel microRNAs and microRNA-stars during silkworm development. BMC Genomics. 2010;11:1–18.
Li S-C, Liao Y-L, Ho M-R, Tsai K-W, Lai C-H, Lin W. miRNA arm selection and isomiR distribution in gastric cancer. BMC Genomics. 2012;13(Suppl 1):S13.
Carroll AP, Goodall GJ, Liu B. Understanding principles of miRNA target recognition and function through integrated biological and bioinformatics approaches. WIREs RNA. 2014;5:361–79.
Akbari Moqadam F, Pieters R, den Boer ML. The hunting of targets: challenge in miRNA research. Leukemia. 2013;27:16–23.
Jalvy-Delvaille S, Maurel M, Majo V, Pierre N, Chabas S, Combe C, et al. Molecular basis of differential target regulation by miR-96 and miR-182: the Glypican-3 as a model. Nucleic Acids Res. 2012;40:1356–65.
Ren J, Jin P, Wang E, Marincola FM, Stroncek DF. MicroRNA and gene expression patterns in the differentiation of human embryonic stem cells. J Transl Med. 2009;7:20.
Bar M, Wyman SK, Fritz BR, Qi J, Garg KS, Parkin RK, et al. MicroRNA discovery and profiling in human embryonic stem cells by deep sequencing of small RNA libraries. Stem Cells. 2008;26:2496–505.
David L, Polo JM. Phases of reprogramming. Stem Cell Res. 2014;12:754–61.
Portt L, Norman G, Clapp C, Greenwood M, Greenwood MT. Anti-apoptosis and cell survival: a review. Biochim Biophys Acta. 1813;2011:238–59.
Khandelwal N, Simpson J, Taylor G, Rafique S, Whitehouse A, Hiscox J, et al. Nucleolar NF-κB/RelA mediates apoptosis by causing cytoplasmic relocalization of nucleophosmin. Cell Death Differ. 2011;18:1889–903.
Schiaffino S, Mammucari C. Regulation of skeletal muscle growth by the IGF1-Akt/PKB pathway: insights from genetic models. Skelet Muscle. 2011;1:4.
Fung MM, Rohwer F, Mcguire KL. IL-2 activation of a PI3K-dependent STAT3 serine phosphorylation pathway in primary human T cells. Cell Signal. 2003;15:625–36.
Pallero MA, Elzie CA, Chen J, Mosher DF, Murphy-Ullrich JE. Thrombospondin 1 binding to calreticulin-LRP1 signals resistance to anoikis. FASEB J. 2008;22:3968–79.
Indran IR, Tufo G, Pervaiz S, Brenner C. Recent advances in apoptosis, mitochondria and drug resistance in cancer cells. Biochim Biophys Acta. 1807;2011:735–45.
We thank Professor Yang-Mooi Lim, UTR, Selangor, Malaysia, for providing the cancer cell lines used in this study.
This work was supported by a Universiti Tunku Abdul Rahman UTARRF grant (6200/CB7), Selangor, Malaysia, to KBC and the Ministry of Higher Education Malaysia (FRGS/2/2014/SKK01/UTAR/01/2).
Availability of data and material
The supplementary materials of the manuscript are as followed:
Additional file 1: Table S1. List of putative target genes
Additional file 2: Table S2. KEGG pathways of the Group I C19MC miRNAs
Additional file 3: Table S3. Predicted target genes of group I C19MC miRNAs related to apoptosis
PNNN, CJH and KBC conceived and designed all the experiments and did the data analysis in the study. PNNN performed the experiments with contribution from CJH in bioinformatics analysis. PNNN and KBC wrote the paper. SS provided the iPSC and hESC cell lines and advice on culture. SKC provided advice and resources to conduct the experiments. All authors have read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate