RNA Genes: Retroelements and Virally Retroposable microRNAs in Human Embryonic Stem Cells.

Embryonic stem cells (ESCs) are capable of undergoing self-renewal, and their developmental ability is known as the stemness. Recently, microRNAs (miRNAs) as regulators have been isolated from ESCs. Although Dicer and DiGeorge syndrome critical region gene 8 (DGCR8) are essential factors for the biogeneration of miRNA, Dicer-knockout (KO) ESCs have showed to fail to express differentiation markers and DGCR8-KO ESCs have showed to be arrest in the G1 phase. Furthermore, Dicer-KO ESCs lost the ability to epigenetically silence retroelemtns (REs). REs are expressed and transposed in ESCs, whose transcripts control expression of miRNAs, and their transposable retroelement (TE) expression is, therefore related to ESC proliferation and differentiation, suggesting that the interplay between miRNAs and REs may have a deep responsibility for the stemness including a short G1/S transition and for RE regulation in ESCs.

INTRODUCTION

Today we know that the human genome encodes non-coding RNA (ncRNA) from its non-coding DNA (ncDNA) sections. It is known that ncDNA is predominantly composed of multiple segments of transposable retroelements (TEs) such as, long interspersed elements (LINEs) and, short interspersed elements (SINEs, human Alu), retrotransposons, DNA transposons and satellites (Fig. ) [1]. The ncRNAs contain ribosomal RNA (rRNA), transfer RNA (tRNA), small nuclear RNA (snRNA), small nucleolar RNA (snoRNA), short interfering RNA (siRNA), microRNA (miRNA), other classes of small RNAs and longer RNA intermediates including complex patterns of interacting and overlapping sense and antisense transcripts from both coding and non-coding regions of genomic DNA [2, 3]. The small RNAs are functional and their sizes fall within a range of 19 to 31 nucleotides (nts). The most common small RNAs are miRNAs, which consist of genomic and resident miRNA genes, which are subsequently produced by digestions on two kinds of the nuclease III, Drosha and Dicer, and then directly target messenger RNAs (mRNAs) or possibly epigenetic machinery in human cells [4].

Embryonic stem cells (ESCs) including somatic and germline stem cells (GSCs), employ their capacity for unlimited self-renewal and differentiation for tissue development and homeostasis. In mammalian embryo, the first stage of this differentiation involves differentiation into either the cells which form the inner cell mass and those which form the trophoblast [5-7]. The cells of the inner cell mass, which are known as ESCs, are pluripotent, meaning that they can give rise to all the cells of the animal. This is quality known as stemness and this stemness retained by the ESCs themselves. Thus, the stem cells can produce differentiated daughter cells, as well as replicate themselves. The key transcriptional regulators Oct4, Sox2 and Nanog, etc. have been identified as the self-renewal and differentiation-related factors in ESCs [8-12]. Artificial pluripotent stem (iPS) cells demonstrate the roles played by these transcriptional regulator genes [13], but these cells developed into tumors (due to the operation of myc oncogene and oncoviral vectors). Recent cancer research indicates that alteration of miRNA profiles, such as those involved in oncogenic miRNA (oncomirs) upregulation or tumour suppressor miRNA downregulation could induce tumors even without an oncovirus or oncogenes [14]. This suggests that tumourigenesis of iPS cells might be due to dysregulation of miRNAs induced by transfection of transcriptional factors rather than transduction of the oncogene myc or oncovirus vectors themselves [15-17] as evidenced by tumourigenesis in iPS cells without oncovirus vectors and myc oncogenes [18]. Since miRNA genes are mobile and functional genetic elements [19], the miRNA profile of ESCs is likely to be different in depending on the niche cells in which they are located when it is self-renewal or differentiation conditions. The character of ESCs may be very different from that of cancer such as embryonic carcinoma and carcinoma-like cells including iPS cells [20]. Based on research with skin cell iPS cells and ESCs, it has become clear that miRNAs play an important role in the precise mechanisms regulating self-renewal and differentiation in these cells [21].

It appears that miRNAs may be the master regulators of stemness. Further, it also appears that tumourigenesis may be rooted in something other than the transcriptional and cellular signaling factors of oncoproteins. The ESC transcripts, including those of miRNAs, contain many retrovirus-like retrotransposon-derived repeats (about 50%) [22, 23]. De novo LINE-1 (L1) retrotransposition events occurred in the human brain [24]. Further, L1 retrotransposition events have been shown to occur at early stages in human embryogenesis [25] and L1 transcripts are related to the mobilization of human SINE Alu [26]. The miRNA genes are located within Alu and L1, and about 20% of known human miRNAs are derived from these repetitive-element sequences [27]. It has been reported that over 1,000 genes of the human genome may be controlled by an Alu-associated promoter and that a dual relationship exists between the miRNA gene clusters and the miRNA targets in Alu regions that may have evolved together [28, 29]. Furthermore, epigenetic chromatin remodeling of Alu sequences plays an important role in the regulation of miRNA expression [30]. Taken together, these findings suggest that miRNAs may be the master regulator of ESC self-renewal and differentiation and that the miRNA genes may control both abilities through expression of TEs. Therefore, it appears that the activity of Alu and L1 is more complex than previously documented in artificially-induced transcriptional factors because ESCs specific transcriptional factors do not control the epigenetic state except for that of polycomb-group proteins [31].

The seed sequences of 2- to 7 or 8-nts in the 5’-end of miRNAs recognize not only the 3’ untranslated regions (3’UTRs) of its mRNAs but also the coding regions of Oct4, Sox2 and Nanog [32] through an imperfect match to repress target mRNA and probable stability of the mRNA target. Computer analysis has shown that miRNAs target the 3’UTR as well as the 5’UTR and protein coding regions (CDS) [33]. Further, experimental evidence indicates that the miRNA genes could target the CDS and/or 3’UTR of OCT4, Sox2 and Nanog [34] and could regulate all multi-gene expression. The miRNA-induced silencing complex (miRISC) appears to regulate this process. The centre of miRISC is an Argonaute protein (Ago) and the seed region of miRNAs in the miRISC targets the sequences of mRNA with incomplete paring. However, a miRNA gene superposes to multiple-targets and a mRNA is targeted by multiple-miRNAs. Therefore, elucidation about how the targets of miRNA selected and how the miRNA genes switched on and off requires additional insights from theoretical ciphering of RNA codes even if suppose a quantum computing model independent of the nucleic acid pairing [19].

In this paper, we explore the significant expression of ESC miRNAs induced by TEs and we propose a model involving TEs forming ‘RNA wave’ in ESCs. To further understand an importance role of fine-tuning by the miRNA genes in ESC epigenetic silencing, we also explore an extended model of quantum theory for RNA coding operating without the operon involvement.

ALU TEs AS A SWITCH FOR STEMNESS

TEs in the human genome have been shown to occupy ~46% of its genomic sequences (Fig. ) [1]. Since TEs transcribe non-coding RNA and this RNA is reverse-transcribed and re-integrated into its own genome, if miRNA is encoded in TEs, the miRNA genes in TEs could also be duplicated according to the same process as TEs. It has been reported that human miRNAs may be derived from retroelements (REs) [29, 35, 36]. It has been suggested that around half of all human genomic information is derived from RNA. Further, human MADE1, miniature inverted-repeat TE (MITE) involves total 550 paralogs upon miR-548a-1 (Chromosome 6: Chr6), miR-548a-2 (Chr6), miR-548a-3 (Chr8), miR-548b (Chr6), miR-548c (Chr12), miR-548d-1 (Chr8) and miR-548d-2 (Chr17) as numbers of 24, 81, 82, 23, 124, 71 and 145 paralogs in individual genome of humans, respectively. MADE1 emerged along the primate evolutionary lineage and the miR-548 gene is an orthologs in the human, chimpanzee and rhesus macaque genomes [35]. Further, miR-548 targets MITEs of the human genome. In turn, MITEs including MADE1 are widely distributed among eukaryotes and miRNA and endogenous siRNA genes may exist due to the palindrome-like RNA structures form a hairpin, which can be catalyzed by Dicer. Thus, a self-regulation system may be involved in MADE1 and the miRNA genes. This system may be necessary for the maintenance of the small RNA-derived human genome. Subsequently, transposition of TEs seems to be deleterious to the genome. However, the genome itself may be constructed of and evolved by small RNAs in TEs. Previously, we proposed “RNA wave model,” which holds that retroviruses control their own amplification rate in order to minimize their deleterious effects on their hosts, as exemplified by the latent state of human immunodeficiency virus type 1 (HIV-1) [37]. According to RNA wave, RNA tags can transport miRNA genes, which are a mobile genetic element that allows cross-talk among cells, organs and individuals from vertical to horizontal. Further, it has been reported that human TE Alu evolved symbiotically with miRNAs [29, 38].

The relationship between Alu retroelements, one of the SINE TEs, and a miRNA cluster on chromosome 19 (C19MC) in humans is compatible with RNA wave [28, 39]. As mentioned above, members of the C19MC-gene family, the miR-515-1, miR-517a and miR-519a-1 genes, are transcribed by Pol III from the upstream section of Alu TE as miRNA promoters, suggesting that some interplay between miRNAs and TEs may act to the circuit on and off. The human miRNA gene cluster, including SINEs, miR-371, miR-372 and miR-373 in the down stream of C19MC, and the murine miR-290, miR-291a, miR-291b, miR-292, miR-293, miR-294 and miR-295 cluster are similar to ESCs in their self-renewal capability (Fig. , ). Furthermore, murine miR-290 cluster genes are transferred from ESCs to neighbouring cells by microvesicles (Fig. ) [40, 41]. Recently, synthetic miRNAs, viral miRNAs and endogenous miRNAs have been reported to be transferred between cells and to be able to spread their functionally silencing signals across cell boundaries [42, 43]. In addition, the vault particles implicated in intracellular transport contain small RNAs and these small RNAs may be related with the multidrug resistance mechanisms [44]. Thus, the direction of self-renewal or differentiation in ESCs and their neighbors may be commonly mediated through the information stored in the retroposable miRNA genes. Additionally, TEs is not a tag in this case, TE itself may be a RNA gene.

Alu transposition is required for LINE TEs but that is not dependent on LINE-encoded protein [26]. Two miRNAs, miR-95 and miR-151 in LINE2 target its LINE2 sequences [45]. The transcripts including miRNAs in ESCs contains many TE-derived repeats (~50%) [22, 23]. LINE1 retrotransposition events occur at early stages in human embryogenesis and human ESCs in vitro [24, 25]. Thus, the miRNA genes in TEs are completely disregarded in ESC investigations for human embryogenesis; this is a mistake at all and TEs may be a switch of self-renewal and differentiation of ESCs via miRNAs. Therefore, nobody has asked what is the code of RNA intermediates, such as miRNA, to elucidate the switch operation for ESC self-renewal and differentiation much shorter nanoRNA and TEs, furthermore intron and protein coding RNAs.

TRANSFORMATION AND REPROGRAMMING IN STEMNESS

Over 10 years ago, embryonic clone cells were produced from the nuclei of cultured cells [46] and then the first mammalian clone was produced [47]. DOLLY, the cloned ovine, proved that genome modification has been considered as irreversible event, but the genome can be reversible to able to reprogram by factors in the oocytes to make them pluripotent once again. The question is what is the factor(s). Recent research indicates that it may be miRNA genes [48, 49]. The initial profile of cloned bovine embryos was changed after somatic cell nuclear transfer as the reprogramming [49]. Further, the nuclear transfer with minimal contamination by the donor cytoplasm was accomplished by metaphase II G1-arrested G0 quiescent state cells [47]. This technological finding as reprogramming has forced us to accept that two factors may be in play here: resident miRNA genes and genomic miRNA genes.

‘Reprogramming’ was also used to create artificially induced adult somatic iPS cells. Recent papers by Werning et al. [50], Hanna et al. [51] or Takahashi et al. [52] gave the results that iPS cells from skin fibroblasts were isolated upon transformation with a combination of the transcription factor genes, Nanog, Sox2, Oct4 and/or Myc. These iPS cells were problematic in that Myc or vector retroviruses induced tumour formation, and iPS cells can essentially allow for isografting only, while ESCs can be considered for use of variable allografting [51-54]. Together with above nuclear transfer, experiments with iPS have caused biologist to completely reassess their thinking regarding the role of transcriptional factors, miRNA genes and other small nucleic acids (miRNA assessment).

It was originally thought that tumourigenesis might be avoided in iPS cells because it was believed that tumourigenesis was mediated by the oncogene Myc and other oncovirus vectors. However, it is now clear that it is not oncoviruses but oncomirs that can induce tumours and that oncogenes are controlled by the oncomirs or tumour suppressor miRNAs [17]. Further, iPS and ESCs up-regulate the expression of the miR-302 and 17-92 cluster, which is the oncomir in leukemia and lung cancer, and iPS cells show lower miR-371-cluster expression than ESCs, in which cluster expression is involved in self-renewal [55]. MiR-222 promotes continuous proliferation of human carcinoma cells [56-58]. In contrast, miR-222 is involved in helping human endometrial stromal cells to exit the cell cycle and enter the differentiation [59]. Therefore, at least miR-222 and likely other miRNA sets play an important role in reprogramming cells. Although the miRNA profiles represent clear difference from iPS cells to ESCs, the ability to reprogram cells throws the existing understanding of stemness into disarray.

Since transcriptional and post-transcriptional regulatory network could be determined in silico [60, 61], miR-145 controls Nanog, Sox2 and Oct4 in ESCs as well as in tumourigenesis [62, 63]. Murine miR-134, miR-296 and miR-470 target the protein-coding sequence of Nanog, Oct4 and Sox2 in various combinations, leading to transcriptional and morphological changes in differentiating murine ESCs [32, 64]. Although it has been shown that Nanog, Sox2 and Oct4 transcriptional regulators collaborate to form regulatory circuitry consisting of autoregulatory and feed-forward loops, which contribute to pluripotency and self-renewal of human ESCs [65], the miRNA findings discussed above suggest that reprogramming of differentiated adult cells by transduction may be a “phantom of the operon.” Later to further distinguish between transformation and reprogramming, the DGCR6 and Dicer gene were knocked out in ESCs.

DGCR8 OR DICER-DEFICIENT ESCs

DCGR8 cooperates with the RNase III enzyme Drosha in the biogenesis of pre-miRNA from pri-miRNA transcribed from non-coding regions of the human genome. As shown in Fig. (), ESCs in which the DCGR8 gene had been knockout (KO) are arrested in the G1 phase of the cell cycle, but these DGCR8-null ESCs remain alive [66]. It has been suggested that resident miRNAs are present independently in the cytoplasm, and these miRNAs may be able to self-proliferate the division of ESCs. Expression of markers of primitive ectodermal-FGF5, and epithelial-KRT18 was delayed in embryoid bodies (EBs), and DGCR8-KO EBs did not form a cyst. Therefore, DGCR8-null ESCs have a defect in their differentiation ability. The pluripotent markers of Oct4, Rex1, Sox2 and Nanog were not expressed [66]. It has been suggested that the DGCR8 gene is essential for the self-renewal silencing that normally occurs with the induction of differentiation. Surprisingly, it has been reported that by aberration of another RNase III enzyme Dicer ESCs are viable and retain their colony-forming ability and the ability to express Bmp4, HINF4A, GATA1, Brachyury and HPRT ESC-specific genes [67]. The Dicer-null ESCs show a prolonged G1 phase of the cell cycle and differentiation defects [68]. These results strongly support two elements of the RNA wave model: 1) the presence of resident miRNAs that might play an important role the self-renewal function of ESCs, and 2) the existence of genomic miRNAs, which are necessary for ESC differentiation (Fig. ).

By the high through-put analysis of the miRNA, ESC-specific miRNAs have been cloned and isolated. When human ESCs-specific miRNAs grouped by ‘AAAGUGCU’ seed sequence were shown, miR-17-5p, miR-18a, miR-18b, miR-20a, miR-20b, miR-93, miR-100a, miR-106b, miR-302a, miR-302b, miR-302c, miR-302d, miR-373, miR-519b, miR-519c, miR-519d, miR-520b and miR-520e including the same seed region were upregulated. These prominent upregulated miRNAs were clustered in human Chr 4, 13, 19 and X [69]. Comparison between human and mouse ESCs revealed that human miR-302b, 302c and miR-302d on Chr 4 are related to miR-302 on murine Chr 3. Another cluster of miRNAs, miR-371, miR-372, miR-373*, and miR-373 on human Chr 19 is the murine homologues of miR-290, miR-291a, miR-291b, miR-292, miR-293, miR-295 on Chr 7 [70, 71]. These miRNAs belong to a miRNA gene family common to human and murine (miRBase 13.0; see Fig. ). The miR-302 and miR-371 genes are regulated by KLF4, OCT4 and SOX2 transcriptional factors and miR-145 suppresses these three transcriptional factors in human ESCs [63, 70]. In the case of murine, miR-470 and miR-134 target the coding region of Oct4 and Sox2 mRNAs, respectively. Further, miR-470, miR-296 and miR-134 target to 3’UTR of Nanog mRNA [64]. According to miRBase, Nanog and SOX2 are targeted by 12 and 26 kinds of multi-miRNAs in human ESCs, and 17 and 26 kinds of multi-miRNAs in mouse, respectively (miRBase 13.0). In these miRNAs, miR-21 of both human and murine has binding sites in the 3’UTR of the mRNA that codes for SOX2, and probably OCT4 [72]. Like polycomb [73], the RE1-silencing transcriptional factor (REST) interacts with cis-elements in promoter of the miR-21 gene, REST represses expression of miR-21. Not only REST, but also the ESC-related transcriptional factors OCT4/SOX2/Nanog/TCF3 bind to the promoter of ESC-specific miRNAs [74]. For example, expression of pri-let-7g is dependent on OCT4 and let-7g is abundant in ESCs [75]. Although Lin28 inhibits expression of let-7g, ectopical expression of Lin28 promotes the reprogramming of human somatic cells into iPS cells [76] and is associated with human multiple cancer types [77]. Further, miR-302 reprogrammed human skin cancer cells into iPS cells [78]. Although quite recently it has been reported that down-regulation of the p53 guardian by siRNAs of p53 can increase the emergence rate of iPS cells produced by skin cells [79-83], we well know that the dominantly inherited Li-Fraumeni multi-carcinogenesis syndrome (MIM151623) as a case of mutations of p53 in human [84, 85]. Since the p53 was regulated by miR-34a and miR-125b [86, 87], these data suggest that production and proliferation of iPS cells may closely be controlled by miRNAs, and that ESCs are maintained not just by ESC-specific transcriptional factors but also by ESC-specific miRNAs as multiple factors. However, it is still not clear how the miRNA genes switch stemness on or off for a change and tuning of cell phenotypes.

One possible answer is an epigenetic regulation mechanism such as that involving miR-290. MiR-290-295 cluster in murine and are indirectly involved in epigenetic gene regulation. This cluster repressed DNA methyltransferase (DNMT), retinoblastoma like-2 (RBL2) and CpG methylation of the promoter by suppression of the DNMT regulated expression of OCT4 [88, 89]. Further, Lin28 regulates the production of the histone, H2a [63], suggesting that miRNAs may be indirectly involved in epigenetic gene regulation in ESCs. But there is a further question: why was only Lin28 finally suppressed via regulation of histone?

THE KEY FOR STEMNESS MAY BE FOUND IN THE LATENCY OF HIV-1

DNA and histone methylation is one of the epigenetic modifications of chromatin, which may be related to ESC programming [88, 89]. In ESCs, programming of cells results in loss of competency for the abbreviated G1 phase, such as independence of the release of RB protein from E2F transcriptional factor [90]. Since complete silencing of the typical cell cycle gene expression is necessary for the programming of ESCs, it is of great interest that miRNAs could be related to the epigenome. The heterochromatin domains are highly methylated and are also show histone modifications characteristic of silenced chromatin, such as trimethylated histone H3 lysine 9 (H3K9me+) [91]. Further, Dicer is essential for the formation of heterochromatin [92], suggesting that miRNAs control the methylation mechanisms of chromatin [93, 94]. In order to change from self-renewal to differentiation, self-renewal set of miRNAs derived from Alu-miRNAs must be switched off and differentiation set of miRNAs from Alu-miRNAs should be switched on. Recently, human miRNA-offset RNAs (moRNA) have been discovered [95]. Although chordate moRNAs are a relative small fraction of the miRNA family, moRNAs are produced from pri-miRNAs probably in Alu or L1 and could target to miRNAs and mRNAs, suggesting that in the case of human, TEs may play an important role in the epigenetic silencing of miRNA and/or mRNA expression. As in the above-described MADE1, human siRNA candidates from MADE1 have been revealed by computer analysis. Therefore, these data suggest that human moRNAs or siRNAs in TEs sought to be related to epigenetic silencing. Thus, the resident miRNAs, moRNAs, siRNAs and unusually small RNAs (usRNAs), which is much smaller than canonical miRNA [96] and the so-called nanoRNAs (nnRNAs) source from TEs are speculated as precursors of transcripts from retrovirus in ncDNA regions (Fig. ), suggesting that the function of HIV-1 proviruses may be similar to that of TEs in epigenetic regulation.

The typical model of complete and epigenetic silencing in TEs is like to be the latency of HIV-1 infection [97]. The integrated provirus corresponding to MADE1 could produce miRNAs from the 3’UTR (Fig. ). The targeting the 5’UTR by miRNAs can lead to histone 3 protein H3K9 methylation (H3K9me+) (Fig. ) and/or DNA methylation. Subsequently, epigenetic LTR-miRNA (Alu-miRNA) switch-off is occurred following histone deacetylation by histone deacetylase (HDAC) and the profile of the resident miRNAs changes in the cells. This is how the latent state is induced and the latent state is in control of the TEs. In the case of ESCs, programming should be initiated from the TE silencing like the latency.

On the other hand, started up reprogramming in iPS cells, it is resembles the reactivation state of HIV-1. Activation of LTR-like Alu promoter and/or environmental stimulation such as that with IL-2 and TNF-α induce transcriptional factor HIV-1 Tat corresponding to OCT4, SOX2 and Nanog in iPS cells (Fig. ). High concentration of transcriptional factor in the cells inhibits the activity of normal expression of genomic miRNA genes via suppression of Dicer activities, like Tat as an activator of RAS [37]. The H3K9me+ is reversed [98] and H3K4 plus H3K36me+ as well as CpG DNA methylation induce silencing of the LTR-miRNA (Alu-miRNA) (Fig. ). Therefore, transduction of OCT4, SOX2 and Nanog might induce an aberrant TE-expression state in differentiated skin cells and iPS cells may be occasionally turned into the tumourigenic cells by an aberrant miRNA gene-expression in TEs because Dicer loss promotes tumourigenesis [99]. To maintain the self-renewal state, the Alu-miR371 family is needed to switch on of self-renewal and to control TEs. In transcriptional factors dominant state, the Alu-miR371 cannot work at all. Consequently, it is clear that reprogramming could never be achieved only by transformation such as that in iPS cells because even ESCs decrease their pluripotency by epigenetic alteration of the expression of TEs in murine cells [100]. The human genome in ESCs always represents latent states between miRNAs and TEs. In turn, ESC genome robustness may be based on a circuit of fragile miRNAs in TEs, in turn based upon the capacious memory of retrotransposons, which may be built on the RNA wave including the niches’ mobile miRNA information.

MOBILE miRNA CODE

Observations of the expression of tens of thousands of miRNAs-TEs suggest that the classical promoter-transcription-protein view is contradicted because the miRNAs in TEs identified have been studied in relative depth and the functions have gradually been investigated in genomic in vivo and in silico studies. In light of the low number of validated miRNAs, bioinformatics may play a key role albeit in small ncRNAs [101], not obeying simple Watson-Crick DNA base pairing. A subset of the Watson-Crick model in the seed sequences is currently used to explain the incomplete target pairing of miRNA. A rather broader set of additional rules may be useful to explain how some small RNAs to develop bioinformatics approaches. The seed sequences of miRNAs target the 3’UTR as well as the 5’UTR and protein coding regions (CDS) based on computer analysis [98]. The miRNA genes could target CDS and/or 3’UTR of OCT4, Sox2 and Nanog [34]. Further, a miRNA, such as miR-24 regulates expression of E2F2 by recognizing out of the seed rule [102]. The resident miR-24 gene with miR-155 and miR-17 is specifically involved in the differentiation of granulocytic, monocytic megakaryocytic, and pro-B lymphoid as somatic stem cells [103]. In zebrafish, miR-214 without canonical seed pairing can effectively target a mRNA for silencing [104]. These data suggest that tuning by miRNAs is probably a high-speed event. Therefore, except for complete and partial seed base-pairing, we should add a neo-mechanism for the regulation by miRNA.

Recently, Fabian et al. [105] have reported that miRNA-mediated translational repression precedes target mRNA deadenylation. However, their model still depends on the seed sequences of miRNA-the target ones of mRNA. Further, there is no specificity of miRISC to bind to a specific miRNA. This is different from aminoacyl-tRNA synthetase, which recognizes a specific amino acid and the tRNAs corresponding to that amino acid. This yields the hypothesis that circular shaped mRNA produced by the cap-binding factor eIF4E plus eIF4G and poly(A)-binding protein PABP could make miRISC into a linear strand of mRNA. Since it is well known that episomal HIV-1 DNA forms the ring, if miRNA is also circular and miRISC including helicase activities links a miRNA ring, as shown in Fig. (), the miRNA rings may enter into the mRNA strand, resulting in probable inhibition of mRNA circularization [106].

Quite recently, it has been reported that HIV-1 replication and reactivation in patients with chronic HIV-1 infection are inhibited by shRNA (shN366) including HIV-1 miRNA MIRN367 [107]. This report suggests that TE-miRNA autoregulating system is involved in the latency model in humans. Therefore, for example of the miRNA gene ring, we represent here the miRNA rings in HIV-1 viral RNA [37], additionally in differentiation and in self-renewal specific sets into ESC (Fig. ). In Fig. (), according to electrostatic potential of each atom of the RNA bases with methyl residue instead of ribose 5-phosphate by the fragment molecular orbit methods (FMO) based on quantum theory [108], at first total electrostatic potential of A (approximate - 0.05 milli-Hartree), U (- 0.03), G (0.258) and C (- 0,002) base were calculated. Although in the original code, hypothesized long ago, the archaic code (G, C) would expand toward the complete variants (C, G, A, U) in the early evolution of mRNA [109, 110], only G had a positive potential. The calculated potencies were represented as a radar shape (Fig. ). We focused on the neighbouring bases and checked G multiplication in sequences of each miRNA (see black dots in Fig. ). In the seed sequences (see shaded diamonds in Fig. ), there was no specific relation between the G multiplication and miRNA specific sets in ESCs. But except for MIRH1 and MIR29b-1, black dots were not involved in the seed sequences of other miRNAs. Further, we found unbalanced and balanced electron potencies of the G base cluster in some miRNAs. In such miRNAs, electrons in the circular miRNAs move from one base to a neighboring base according to electrostatic potential of each atom of the RNA bases with FMO, whereas electrons may spin round. Depended on the direction of rotating electron, the direction of the magnetic field would go up or down. The superposing of the magnetic direction of electrons according to quantum theory may be the code of RNA (Fig. ), like quantum computing [19]. The pulse waves of the electrons in the organ hydrogen (H) atoms of water in magnetic resonance image (MRI), could, in part, be derived from a source of the magnetic waves produced by miRNA superposing for silencing gene expression, such as aberrant presence of miRNAs in tumor [17]. By measure miRNA electron spin, it may be available for investigation of the relationship between stemness and physicochemical characters of tuning by miRNAs and furthermore therapy of RNA diseases.

CONCLUSIONS

The RNA wave consists of four alternative concepts for crosstalk among cells: 1) mobile miRNAs induce transcriptional and post-transcriptional silencing processes; 2) the induced silencing information expands to intra-cellular, inter-cellular and intra-organ, inter-organ communication by free mobile miRNAs; 3) the mobile miRNAs can self-proliferate; and 4) cells contain both resident and genomic miRNAs. Thus, the non-selfish miRNAs are regarded as the basic unit of the evolution to move horizontally in genomic space. The miRNAs would horizontally transpose in TEs by feeding of miRNAs-containing foods, viruses, etc. The retroposable miRNAs might create phenotype by changing and tuning of the gene expression in the genomic space, in which its genotype can be vertically inherited under pressures of the environment corresponding to Darwininsm.

To clarify the capacity of ESCs to renew themselves and to differentiate into various cell lineages, miRNAs and TEs as RNA genes have an important role. The highlight of the importance of miRNAs in TEs has been shown in controlling ESC growth and differentiation by miRNAs. Although the significance of miRNA function has not yet been evaluated, retroposable miRNAs in TEs are expressed in ESCs, suggesting that TEs may switch on or off the alteration of ESCs from self-renewal, differentiation and reprogramming including the epigenetic regulation. This paper will aid in our understanding of the relationship between oncoviral miRNA genes and TEs. Understanding the seed theory for miRNA targets and further breaking the RNA code under the RNA wave will shed light on the quality of stemness in ESCs.

This article also contains supplementary material (video animation) and it can be viewed online along with the article at publisher’s website.

Table 1.

The General Abbreviations and Glossary for This Paper

Abbreviation	Technical Terms	Explanation
ESC	Embryonic stem cell	Undifferentiated, pluripotent cells from a mammalian embryo. A key cell of retroelements and microRNAs interplay in this paper.
miRNA (miR)	MicroRNA	Mobile and small genetic RNA elements. Small RNA elements are encoded from non-coding DNA region including retroelements as genes. Exosomes and viruses contain microRNAs therefore that can be mobile and control transcription and translation in host.
RE or TE	Retroelement or transposable retroelement	Mobile genetic RNA element. RNA elements are reversetranscribed by reverse transcriptase that makes a cDNA copies. The DNA copies integrated into prototype genome by integrase that are composed of about half of human genome. The DNA elements in genome contain pol gene that encodes reverse transcriptase and integrase and can therefore be copied by its own.
(-)	Retrovirus	An RNA virus. The genetic composition and copying process of viruses are the same as that of retroelements. HIV-1 is a typical retrovirus. A particle of virus contains about 10 kilobases of two single-stranded (+) RNA.
Chr	Chromosome	The coiled coil DNA and histone proteins that is believed to be as the basic conformation. The number of chromosome in a human somatic cell is twice 22 plus X and X (female) or X and Y (male). The DNA of all chromosomes is composed of a genome.
ncRNA	Non-coding RNA	About 98% DNA of 3.2 gigabases of genome is protein non-coding region that therefore has been believed to be nonfunctional region as junk. The region encodes non-coding RNA containing retroelement and microRNA.
DNMT	DNA methyltransferase	The sequence 5'CG3' CpG dinucleotides within genomic DNA are target of a DNA methyltransferase. DNA methylation means conversion of cytosine into 5-methyl-cytosine.
CDS	Protein coding regions	1-1.5% of DNA sequences in genome encode the amino acid sequences of a polypeptide via messenger RNA intermediates. The nucleotide triplet specifies an amino acid or translation stop signal or/and translation start signal.