Interactions of intergenic microRNAs with mRNAs of genes involved in carcinogenesis

miRNAs regulate gene expression by binding with mRNAs of many genes. Studying their effects on genes involved in oncogenesis is important in cancer diagnostics and therapeutics. The RNAHybrid 2.1 program was used to predict the strong miRNA binding sites (p < 0.0005) in target mRNAs. The program Finder 2.2 was created to verify 784 intergenic miRNAs (ig-miRNA) origin. Among 54 considered oncogenes and tumor suppressor genes, 47 genes are the best targets for ig-miRNAs. Accordingly, these genes are strongly regulated by 111 ig-miRNAs. Some miRNAs bind several mRNAs, and some mRNAs have several binding sites for miRNAs. Of the 54 mRNAs, 21.8%, 43.0%, and 35.2% of the miRNA binding sites are present in the 5'UTRs, CDSes, and 3'UTRs, respectively. The average density of the binding sites for miRNAs in the 5'UTR was 4.4 times and 4.1 times greater than in the CDS and the 3'UTR, respectively. Three types of interactions between miRNAs and mRNAs were identified, which differ according to the region of the miRNA bound to the mRNA: 1) binding occurs predominantly via the 3'-region of the miRNA; 2) binding occurs predominantly through the central region of the miRNA; and 3) binding occurs predominantly via the 5'-region of the miRNA. Several miRNAs effectively regulate only one gene, and this information could be useful in molecular medicine to modulate translation of the target mRNA. We recommend described new sites for validation by experimental investigation.


Background:
Increasing interest in the biological role of miRNA has resulted in a number of publications dedicated to these unique regulators of gene expression [1]. Correlations have been shown between changes in miRNA concentration and the development of certain cancers, such as lung [2], esophageal [3], stomach [4], colon [5], liver [6], prostate [7], and breast [8,9]. Cancer diagnostic methods have been developed based on changes in miRNA expression levels during the stages of oncogenesis [10][11][12][13][14]. Searching for binding sites for miRNAs in mRNAs is an important challenge. The conventional point of view is that miRNAs bind primarily to the 3'UTR of mRNA, while other regions of mRNA (specifically the 5'UTR and CDS) do not have a significant role in the regulation of translation by miRNAs [15]. However, there are publications demonstrating miRNA and mRNA interactions in the 5'UTR and CDS [16][17].  inhibits mRNA translation of the lin-41 gene in Zebrafish embryos by binding to a site in the 3'UTR. If this binding site for let-7 appears in both the protein-coding sequence and the 5'UTR, it is likely to also be functional in the 5'UTR [18]. In Drosophila, it has been shown that binding sites in the 5'UTR and CDS are similarly functional. MiR-2 represses translation initiation in the 5'UTR and stimulates mRNA deadenylation [19]. Another article presented data using a reported gene assay that confirmed the presence of binding sites for miR-181a in the coding sequence of ZNF mRNA in several cell lines [20]. Information about the number of miRNA binding sites and the features of miRNA binding could improve the current understanding of interactions between miRNAs and mRNAs. The number of genes regulated by an individual miRNA and the number of miRNAs affecting the expression of an individual mRNA remain unknown [21]. In this research, the binding sites of more 700 intergenic miRNAs within the 5'UTRs, CDSes, and 3'UTRs of mRNAs for 54 human oncogenes were investigated. These genes participate in the development of esophageal, stomach, colon, liver, breast, prostate, and ovarian cancers. The aims of the present work are as follows: a) to identify features of interactions between miRNAs and different regions of mRNAs; b) to establish differences in the ability of miRNAs to bind to different mRNAs; and c) to determine features of nucleotide interactions during the formation of complexes between miRNA and mRNA

Methodology:
Nucleotide sequences of mRNAs of 54 human genes (Homo sapiens Genome build 37.2.) were obtained from Genbank. This set comprises well-known oncogenes and tumor suppressor genes provided by literature review. Nucleotide sequences of ig-miRNAs were obtained from the miRBase database. The program Finder 2.2 was created and used to find genomic origin of ig-miRNAs. All miRNA:mRNA pairs were found by RNAHybrid 2.1 program that provides position of potential binding sites, the free energy value (ΔG) of binding and interaction schemes. As a comparative quantitative criterion of binding was inserted the ΔG/ΔG m value (%), where ΔG m equals the binding energy for a miRNA with a completely complementary nucleotide sequence. Searching for miRNA binding sites was performed along the total mRNA sequence. miRNA binding sites were selected for each target sequence based on the ΔG value and its standard deviation. The binding ability of ig-miRNAs, in-miRNAs, and eх-miRNAs to the 5'UTRs, CDSes, and 3'UTRs of the 54 mRNAs was estimated with a significance of p < 0.0005 (Student's criterion). The density of the binding sites in the 5'UTR, CDS, and 3'UTR per 1000 nucleotides was computed as the number of sites divided by the nucleotide length of the region and multiplied by 10 3 (s/l).

Results: Interactions between intergenic miRNAs and the mRNAs of genes involved in oncogenesis
The interactions between 784 intergenic miRNAs and the mRNAs of 54 protein-coding human genes were explored. It was revealed that 47 mRNAs are targets for these miRNAs, while the mRNAs of seven genes (ABCB1, MSH2, MSH3, MYC, PROM1, SNAI1, and TNFSF10) have no binding sites for ig-miRNAs based on established interaction criteria (see Materials and Methods). Only 111 miRNAs out of the 784 ig-miRNAs bind to these 47 mRNAs with a high interaction energy, at a total of 165 binding sites. Table 1 (see supplementary material) represents the characteristics of ig-miRNA binding sites that have an ΔG/ΔG m value is greater than or equal to 75% (p<0.0005). Each mRNA binds one or several ig-miRNAs. Some of these mRNAs are bound by six or more ig-miRNAs. For example, the mRNAs of the AXIN1, CCND1, CDH1, FLCN, MMP2, SMAD4, and SRC genes have from 6 to 13 sites for ig-miRNAs. This number is greater than the average number of ig-miRNA binding sites for the 47 genes, which is 3.3. Table 1 (see supplementary material) shows that each of the mRNAs has one binding site for a miRNA. Several miRNAs associate with several mRNAs. For example, miR-4472 binds to the mRNAs of seven genes, and miR-1279 has five target mRNAs. There is no relation between the length of the mRNA and the number of ig-miRNA binding sites because the correlation coefficient equals -0.061. The number of binding sites is not proportional to mRNA length, and consequently the established binding sites are not random. Shuffling the mRNA nucleotide sequences reduced the number of binding sites. Of the binding sites for ig-miRNAs, 21.8%, 43.0%, and 35.2% are located in the 5'UTRs, CDSes, and 3'UTRs of the 47 mRNAs, respectively. The density of the binding sites for miRNAs varies among mRNAs. Data from Additional Table 1 (see supplementary material) shows that the density of sites varies from 0.14 s/l (for the mRNA of BRCA1) to 3.09 s/l (for the mRNA of SRC), with an average density of 0.89 s/l for all 47 mRNAs. The mRNA of SRC has the highest binding site density, showing an important role for intergenic miRNAs in the regulation of the expression of this gene.
The mRNAs of the studied genes have different numbers of binding sites for ig-miRNAs in their 5'UTRs, CDSes, and 3'UTRs. The average binding site densities for ig-miRNAs in the 5'UTR, CDS, and 3'UTR of the 47 mRNAs are 2.60 s/l, 0.65 s/l, and 0.68 s/l. The average number of binding sites for miRNAs in the 5'UTR is 4.4 times more than in the CDS and 4.1 times more than in the 3'UTR. These data show that miRNAs can bind not only in the 3'UTR but also in the 5'UTR and CDS of target mRNAs. mRNAs of the GNAS, and PTEN genes interact with five ig-miRNAs at sites all located in their 5'UTRs. The mRNAs of CCND1, MTHFR, SMAD4, and SRC associate with ig-miRNAs preferably in their 3'UTRs Table 4 (see supplementary material). These data indicate the specificity of the interactions between mRNAs and miRNAs. The nucleotide sequences for binding sites between some miRNAs and mRNAs are shown in Table 2 Table 2 show that ig-miRNAs can bind in the 5'UTR, CDS, and 3'UTR with a high interaction energy.

Discussion:
The number of miRNAs identified in humans is constantly increasing, and currently more than 1500 miRNAs are known. Most studies use near 300 well-described miRNAs; our investigation provided for 784 miRNAs the best binding sites with 54 target genes. The identification of target genes for miRNAs depends on the quality of prediction data in silico. Several miRNAs are able to interact with multiple mRNAs [22]. This data suggests that the majority of human genes are regulated by miRNAs. It is known that many miRNAs appear

Conclusions:
Present research data suggest that the miRNAs described here may regulate gene expression more effectively than miRNAs studied in other experimental investigations. ADAM29, BAX, MSH2, PIK3CA genes have only one strong potential site for considered miRNAs. We recommend that described new sites are validated by experimental investigation. Several miRNAs affect only one gene, and this information could be useful in molecular medicine to modulate translation of the target mRNA.

Conflict of interest:
The authors declare that they have no conflict of interest.