Abstract
MicroRNAs target specific mRNA(s) to silence its expression and thereby regulate various cellular processes. We have investigated miRNA gene counts in chromosomes for 20 different species and observed wide variation. Certain chromosomes have extremely high number of miRNA gene compared with others in all the species. For example, high number of miRNA gene in X chromosome and the least or absence of miRNA gene in Y chromosome was observed in all species. To search the criteria governing such variation of miRNA gene counts in chromosomes, we have selected three parameters- length, number of non-coding and coding genes in a chromosome. We have calculated Pearson's correlation coefficient of miRNA gene counts with length, number of non-coding and coding genes in a chromosome for all 20 species. Major number of species showed that number of miRNA gene was not correlated with chromosome length. Eighty five percent of species under study showed strong positive correlation coefficient (r ≥ 0.5) between the numbers of miRNA gene vs. non-coding gene in chromosomes as expected because miRNA is a sub-set of non-coding genes. 55% species under study showed strong positive correlation coefficient (r ≥ 0.5) between numbers of miRNA gene vs. coding gene. We hypothesize biogenesis of miRNA largely depends on coding genes, an evolutionary conserved process. Chromosomes having higher number of miRNA genes will be most likely playing regulatory roles in several cellular processes including different disorders. In humans, cancer and cardiovascular disease associated miRNAs are mostly intergenic and located in Chromosome 19, X, 14, and 1.
Introduction
MicroRNAs are the major key players to regulate the expression of coding genes (Lagos-Quintana et al., ). Genes are divided into two groups—coding genes that code proteins and non-coding genes that do not code functional proteins. More than 98% of human genome is constituted of non-coding DNA (Elgar and Vavouri, ). Recent data from the Encyclopedia of DNA Elements (ENCODE) project suggests that 80.4% of the DNA in human genome participates in “at least one bio-chemical RNA and/or chromatin associated event” (ENCODE Project Consortium, ) implicating the evidence of their regulatory functions. So, it has been established fact that non-coding genes can be transcribed to form non-coding RNAs that can play regulatory role over expression of coding genes. Such non-coding RNAs are found in all eukaryotes like fungi, plants, protozoans, and metazoans (Carrington and Ambros, ). miRNA is a sub-set of non-coding RNAs that are transcribed by RNA polymerase II to form primary miRNA (Pri-miRNA). Nuclear RNase III enzyme Drosha and its co-factor process pri-miRNA into 60 nt length precursor miRNA (pre-miRNA) (Chang and Mendell, ). The precursor miRNA has stem and loop structure and it is exported into cytoplasm by exportin 5 followed by cleavage with Dicer-TRBP complex to form mature miRNA (Zhang et al., 2007). This mature miRNA makes complex with Argonaute (Ago) proteins to form RNA induced silencing complex (RISC) (Nasser et al., ). The functional strand of miRNA stably associated with RISC and the other strand is called passenger strand. Now the mature miRNA guides the RISC complex to repress its targets, mainly at 3′-untranslated region (3′-UTR) of mRNA via short complementary sequence of 6–8 nucleotides (called “seed” sequence), inducing mRNA destabilization, degradation and/or inhibiting translation for protein synthesis (Kim, ; Nilsen, ; Filipowicz et al., ; Fabian et al., ). One miRNA can post-transcriptionally silence several 1000 of genes expressions and a single gene can be silenced by more than one miRNA too (Cui et al., ; Wouters et al., 2011). Tissue-specific expression of miRNA was also observed (Choudhury et al., ; Salvi et al., 2013).
miRNAs can be categorized into four sub-types such as intergenic (dme-mir-1, mmu-mir-7b, hsa-mir-7-2 etc.), intronic (dme-mir-2a-1, mmu-mir-199b, hsa-let-7d etc.), exonic (dme-let-7, mmu-mir-9-2, hsa-let-7a-2 etc.), and others (dme-mir-4949, mmu-mir-331, hsa-mir-632 etc.) depending on their genomic origin (Ying et al., 2010; Chien et al., ). miRNA precursors located in intronic and exonic region of protein coding gene are called intronic and exonic miRNA respectively. The miRNA precursors located in between two consecutive protein coding genes are termed as intergenic miRNA. The others category includes 3′UTR, 5′UTR and combinations of any two from intron, exon, 3′UTR and 5′UTR (Griffiths-Jones et al., ; Godnic et al., ). Various miRNAs are produced as per requirement of the cell but the regulation or detail mechanism(s) of biogenesis of miRNAs are yet to be elucidated. miRNA genes can have transcription start site (TSS) like the coding genes and transcription factors that are used for transcription of mRNA are also used for miRNA transcription (Aguda et al., ; Pichiorri et al., ; Wang et al., 2010a,b). Again, miRNA could also be produced in cell as a by-product of post-transcriptional processing of coding genes using splicing machinery and lariat de-branching enzymes bypassing the conventional nuclear miRNA biogenesis pathway by Drosha cleavage as mentioned above (Okamura et al., ; Flynt et al., ). These miRNAs are termed as “Mirtrons” (Berezikov et al., ; Chan and Slack, ; Ruby et al., 2007; Westholm and Lai, 2011; Havens et al., ).Therefore, it is not very clear whether biogenesis of miRNA is a separate, parallel process as that of coding gene expression or linked with it.
Numerous studies to date have established the role of miRNAs in diverse cellular processes like stem cell differentiation, heart development (Chen et al., ; Zhao et al., 2007; Ivey et al., ), insulin secretion (Plaisance et al., 2006), apoptosis (Lukiw and Pogue, ; Tarasov et al., 2007), aging (Kumamoto et al., ; Maes et al., ), immunity (Chen et al., ; Rodriguez et al., 2007), cell proliferation, metabolism (Hwang and Mendell, ; Taganov et al., 2007; Stefani and Slack, 2008) etc. Recent studies also revealed that certain miRNAs are strongly associated with various diseases such as diabetes (Ciccacci et al., ; van de Bunt et al., 2013), cancer (Calin and Croce, ; Blenkiron and Miska, ; Wuchty et al., 2012; Benetatos et al., ), cardiovascular disease (Zhao et al., 2005; Sayed et al., 2007; Zhao et al., 2007), neurodegenerative diseases (Hébert et al., ; Wang et al., 2008a,b) and hence miRNAs can be treated as potential biomarker or diagnostic tool (Keller et al., ; Fassina et al., ; Jones et al., ; Weiland et al., 2012; Cuk et al., ; Endo et al., ).
There is as such no report available regarding detailed chromosome-specific localization of miRNA genes or precursors in different species. Here we have extensively studied the miRNA gene counts in different chromosomes in various species including human and also collected the miRNAs associated with two diseases-cancer and cardiovascular in human only. We have calculated the correlation coefficient between miRNA gene counts with chromosome length, number of coding and non-coding genes in a chromosome and proposed the evolutionary conserved predominant way of biogenesis of miRNA.
Results
Variation of miRNA gene counts in the chromosomes of different species of metazoa
We have extensively searched miRBase (Griffiths-Jones et al., ) to locate the miRNA precursors in different chromosomes and details of karyotype data are collected from Ensembl genome browser (Flicek et al., ) in 20 different species including human under metazoan. The detail classification of 20 species under study is shown in Figure 1. We have plotted number of miRNA precursors vs. chromosome number in all 20 species as shown in Figures 2A–E. It is evident from those figures that numbers of miRNA genes or precursors are widely varied throughout the all chromosomes of a particular species under study. Instead, there are certain chromosomes having higher counts of miRNA precursors ranging from C. elegans to H. sapiens. The number of miRNA precursors varies from 0 to 161 in a particular chromosome. For example, the highest number, i.e., 161 miRNA precursors are observed in Chromosome 1 whereas only 2 miRNA precursors are observed in Y chromosome in human. Chromosome 6 of R. norvegicus has the highest 67 miRNA precursors out of all the chromosomes but chromosome 16 has only 2 miRNA precursors in the same species. Importantly, the lowest or no miRNA precursors are observed in Y chromosomes, on the contrary high number of miRNAs is observed in X chromosomes in all the species. This observation suggests that presence of high number of miRNA genes in X chromosome and the least or absence of miRNA gene(s) in Y chromosome is an evolutionary conserved phenomenon. Different species have different number of chromosomes and the miRNA gene counts in chromosomes are non-uniform as evident from Figures 2A–E. Furthermore, we have selected four chromosomes of each species having top most number of miRNA genes as shown in Table 1 and calculated percentage of miRNA genes in the said four chromosomes. For example, Chromosome 1, 2, 4, and 3 are the top most miRNA genes containing chromosomes of G. gallus in decreasing order and as high as 39.66% miRNA genes come from these four chromosomes of this species; rest 30 chromosomes contribute ~60% miRNA genes.
Figure 1
Figure 2
Table 1
| Sl no. | Species | Total no. of miRNA (A) | Ch no. (highest no. of miRNA) (B) | Ch no. (2nd highest no. of miRNA) (C) | Ch no. (3rd highest no. of miRNA) (D) | Ch no. (4th highest no. of miRNA) (E) | % miRNA contributed by top four highest miRNA gene containing Ch no. [(B+C+D+E)/A]×100 (%) |
|---|---|---|---|---|---|---|---|
| 1 | S. mansoni | 19 | Ch W(5) | Ch 1(5) | Ch 5(3) | Ch 6(2) | 78.94 |
| 2 | D. melanogaster | 236 | Ch 2L(53) | Ch 2R(52) | Ch 3R(52) | Ch X(50) | 87.71 |
| 3 | C. elegans | 223 | Ch X(67) | Ch II(41) | Ch IV(40) | Ch I(27) | 78.47 |
| 4 | C. intestinalis | 272 | Ch 7(66) | Ch 1(43) | Ch 3(32) | Ch 5(20) | 59.19 |
| 5 | A. carolinensis | 156 | Ch 2(47) | Ch 1(35) | Ch 4(22) | Ch 3(18) | 78.20 |
| 6 | G. gallus | 706 | Ch 1(90) | Ch 2(73) | Ch 4(60) | Ch 3(57) | 39.66 |
| 7 | T. guttata | 225 | Ch 2(25) | Ch 1(24) | Ch 3(18) | Ch Z(16) | 36.88 |
| 8 | D. rerio | 347 | Ch 4(69) | Ch 5(25) | Ch 10(21) | Ch 23(20) | 38.90 |
| 9 | O. latipes | 153 | Ch 7(18) | Ch 4(18) | Ch 14(14) | Ch 17(12) | 40.52 |
| 10 | T. nigroviridis | 113 | Ch 2(10) | Ch 1(10) | Ch 9(9) | Ch 13(9) | 33.62 |
| 11 | S. scrofa | 266 | Ch X(38) | Ch 2(23) | Ch 1(20) | Ch 12(19) | 37.59 |
| 12 | O. anatinus | 106 | Ch X1(38) | Ch 6(9) | Ch 7(9) | Ch 1(8) | 60.37 |
| 13 | M. domestica | 441 | Ch X(143) | Ch 2(63) | Ch 1(61) | Ch 4(40) | 69.61 |
| 14 | M. musculus | 1183 | Ch 2(146) | Ch 12(94) | Ch X(90) | Ch 11(87) | 35.24 |
| 15 | R. norvegicus | 437 | Ch 6(67) | Ch X(58) | Ch 1(49) | Ch 10(33) | 47.36 |
| 16 | B. taurus | 801 | Ch 21(77) | Ch X(62) | Ch 19(54) | Ch 5(52) | 30.58 |
| 17 | M. mulatta | 558 | Ch 7(79) | Ch X(73) | Ch 19(58) | Ch 1(47) | 46.05 |
| 18 | P. troglodytes | 622 | Ch X(72) | Ch 19(63) | Ch 14(59) | Ch 1(51) | 39.38 |
| 19 | G. gorilla | 327 | Ch 14(44) | Ch X(43) | Ch 1(28) | Ch 5(24) | 42.51 |
| 20 | H. sapiens | 1876 | Ch 1(161) | Ch19(142) | Ch X(118) | Ch 2(113) | 28.46 |
Percentage of miRNA gene counts in top four highest number of miRNA gene containing chromosomes.
Correlation of miRNA precursors with chromosomal length, number of coding/non-coding genes present in respective chromosome within a species
Different species have different number of chromosomes of variable length containing variable number of coding and non-coding genes. We have calculated the correlation coefficients of miRNA precursor counts with length of chromosome, number of coding/non-coding genes in all the 20 species as shown in Table 2. Majority of species (55%) under study showed that miRNA gene counts were not correlated with the chromosome length. So, longer chromosomes are not having higher number of miRNA precursors always. On the contrary, there are certain chromosomes in all species where higher number of miRNA genes located even if the length is shorter than rest of chromosomes and hence this is an evolutionary conserved event. For example, in M. musculus chromosome 1 (195.47 Mbps) is longer than chromosome 2 (182.11 Mbps) but the latter has the highest 146 miRNA genes whereas the earlier has 63 miRNA genes only. Then, we looked into the number of coding genes in different chromosomes. It has been observed that shorter chromosome can have more coding genes than longer one and maintained through evolution. The correlation coefficient (r-values, red square) of miRNA precursor counts with chromosomal length is just following the correlation coefficient (r-values, green square) of coding gene counts vs. chromosomal length in different species as shown in Figure 3. This data implicates that chromosomal localization of coding genes is someway related with that of miRNA genes in all the species under study.
Table 2
| Sl no. | Species | Correlation coefficient of number of miRNA precursors in respective chromosomes with | ||
|---|---|---|---|---|
| Chromosome length | Coding genes | Non-coding genes | ||
| 1 | S. mansoni | 0.76 | 0.77 | 0.62 |
| 2 | D. melanogaster | 0.85 | 0.85 | 0.75 |
| 3 | C. elegans | 0.09 | −0.37 | 0.19 |
| 4 | C. intestinalis | 0.54 | 0.50 | 0.81 |
| 5 | A. carolinensis | 0.84 | 0.94 | 0.92 |
| 6 | G. gallus | 0.97 | 0.97 | 0.98 |
| 7 | T. guttata | 0.86 | 0.79 | 0.88 |
| 8 | D. rerio | 0.31 | 0.50 | 0.92 |
| 9 | O. latipes | 0.22 | 0.43 | 0.67 |
| 10 | T. nigroviridis | 0.50 | 0.47 | 0.63 |
| 11 | S. scrofa | 0.26 | 0.41 | 0.64 |
| 12 | O. anatinus | 0.43 | 0.46 | 0.78 |
| 13 | M. domestica | −0.26 | −0.15 | −0.13 |
| 14 | M. musculus | 0.61 | 0.70 | 0.84 |
| 15 | R. norvegicus | 0.46 | 0.45 | 0.80 |
| 16 | B. taurus | 0.29 | 0.45 | 0.78 |
| 17 | M. mulatta | 0.27 | 0.55 | 0.66 |
| 18 | P. troglodytes | 0.26 | 0.55 | 0.50 |
| 19 | G. gorilla | 0.31 | 0.42 | 0.46 |
| 20 | H. sapiens | 0.53 | 0.91 | 0.78 |
Correlation coefficient of number of miRNA precursors vs. chromosome length, number of coding and non-coding genes in the respective chromosome.
Figure 3
Further, we calculated correlation coefficient of coding and non-coding genes with miRNA precursors present in a chromosome of all the species as shown in Table 2. The positive correlation has been noticed between miRNA vs. coding gene counts in all species under study except C. elegans. 55% (75% species if r ≥ 0.45) of species under consideration shows strong correlation coefficient (r ≥ 0.5) between number of miRNA genes and coding genes in a chromosome as evident from Table 2. Eighty five percent of species under consideration shows correlation coefficient (r ≥ 0.5) for miRNA vs. non-coding genes as expected because miRNA is sub-set of non-coding genes. We took four species showing r-value less than 0.5 and four species showing r-value greater than 0.5 and calculated percentage of each miRNA subtypes out of total miRNA precursors within the respective species as shown in Table 3. We have observed that r-value is greater than 0.5 when the sum of intronic and exonic miRNA precursors is greater than intergenic miRNA precursors in several species like H. sapiens, T. guttata, A. carolinensis, and D. melanogaster. On the contrary, r-value less than 0.5 is observed when number of intergenic miRNA precursors exceeds the total of intronic and exonic miRNA counts in species like O. latipes, S. scrofa, O. anatinus, and G. gorilla as shown in Table 3. These findings imply that intronic and exonic miRNA genes are closely related with coding genes in a chromosome and that too evolutionary conserved.
Table 3
| Species (r-value from Table 2) | Total miRNA precursors (under study) | Number of miRNA sub-types (%) | |||
|---|---|---|---|---|---|
| Intergenic % | Intronic % | Exonic % | Others % | ||
| H. sapiens (0.87) | 1876 | 445 (23.72) | 1059 (56.44) | 110 (10.38) | 262 (13.96) |
| T. guttata (0.78) | 242 | 56 (23.14) | 34 (14.04) | 101 (41.73) | 51 (21.07) |
| A. carolinensis (0.93) | 261 | 88 (33.71) | 64 (24.52) | 71 (27.20) | 38 (14.55) |
| D. melanogaster (0.85) | 236 | 75 (31.77) | 115 (48.72) | 21 (8.89) | 25 (10.59) |
| O. latipes (0.42) | 164 | 92 (56.09) | 21 (12.80) | 41 (25.00) | 10 (6.09) |
| S. scrofa (0.38) | 227 | 109 (48.01) | 53 (23.34) | 39 (17.18) | 26 (11.45) |
| O. anatinus (0.35) | 396 | 251 (63.38) | 106 (26.76) | 23 (5.80) | 16 (4.04) |
| G. gorilla (0.41) | 317 | 166 (52.36) | 89 (28.07) | 29 (9.14) | 33 (10.41) |
Number of miRNA sub-types in different species.
Diversity of miRNA sub-types in human chromosomes
We have observed that miRNA genes are highly concentrated into certain chromosomes in all the 20 species including human. We collected total 1876 miRNA precursors from miRBase (Griffiths-Jones et al., ) and categorized them into four groups according to their genomic location such as intergenic, intronic, exonic and others. Chromosome wise variation of sub-types of miRNA gene counts are given in Figure 4. The four top most intergenic miRNA genes containing chromosomes are 19 (14.61%), X (11.91%), 14 (9.44%), and 1 (5.39%) whereas top most intronic miRNA genes containing chromosomes are 1 (10.10%), 2 (7.93%), 11 (5.76%), and 3 (5.38%). The exonic miRNA counts are very less than intergenic and intronic counts in a chromosome. We calculated the correlation of miRNA sub-types with coding genes. Overall correlation coefficients of coding genes vs. intergenic, intronic, exonic and others miRNAs are given by 0.47, 0.91, 0.25, and 0.74 respectively, indicating that intronic and others (3′UTR/5′UTR or combination) miRNA precursors are strongly correlated with coding genes in human.
Figure 4
miRNAs associated with cancer and cardiovascular diseases in human chromosomes
We collected data of miRNA- disease association from HMDD v1.0 (January, 2012) which includes 617 miRNAs (Lu et al., ). We have categorized the miRNAs that are associated with cancer and cardiovascular disease separately and searched their location in chromosomes from miRBase database (Griffiths-Jones et al., ). The variation of counts of miRNA gene associated to the diseases in different human chromosomes is shown in Figure 5. We have observed that the counts of disease related miRNA gene widely varied into all chromosomes as expected but concentrated in certain chromosomes. Chromosomes 19, X, 14, 1 are the major cancer-associated miRNA gene containing chromosomes whereas chromosome 19, 14, 1, X are the major cardiovascular disease-associated miRNA gene containing chromosomes in decreasing order. Therefore, all the four chromosomes 19, 14, 1, and X are the major source of these two disease associated miRNA genes.
Figure 5
Furthermore, we have arranged the miRNAs associated with cancer and cardiovascular disease according to their genomic origin like intronic or intergenic using miRBase database resources and plotted number of miRNA sub-types associated with cancer (Figure 6A) or cardiovascular disease (Figure 6B) vs. chromosome number. Bar diagram of Figures 6A,B shows that disease-associated intronic or intergenic miRNA gene counts are different in different chromosomes. For example, maximum intergenic miRNAs associated with both cancer and cardiovascular diseases are from Chromosome 19 followed by chromosome 14 and X as shown in Figures 6A,B. On the contrary, the highest number of intronic miRNAs associated with cancer and cardiovascular diseases are from chromosome X and 9 respectively.
Figure 6
We calculated the percentages of miRNA sub-types (intergenic and intronic) associated with cancer and cardiovascular disease out of total miRNA precursors of that particular sub-type (Table 4). We see that 45.84% and 24.26% intergenic miRNAs are associated with cancer and cardiovascular disease respectively. In a nutshell, higher proportions of inetrgenic miRNAs are associated with these two diseases compared with intronic sub-type and their major sources are chromosome 19, 14, and X.
Table 4
| Sub-type | Total counts of miRNA in human | Associated with cancer | Associated with cardiovascular disease | ||
|---|---|---|---|---|---|
| Total miRNAs | % | Total miRNAs | % | ||
| Intronic | 1059 | 180 | 16.99 | 60 | 5.66 |
| Intergenic | 445 | 204 | 45.84 | 108 | 24.26 |
Percentage of prevalence of miRNA sub-types associated with diseases.
Discussion
We have investigated variation of miRNA gene counts in all the chromosomes of 20 different species. To investigate the criteria/logic behind such variation of miRNA gene counts in chromosomes, we have selected three parameters such as chromosomal length, number of coding and non-coding genes present in a chromosome and evaluated correlation coefficient of miRNA gene counts with length of chromosome, number of coding and non-coding genes. We have observed majority (55%) of species are not following the rule—“higher length chromosome contains higher number of miRNA genes.” So, length of chromosome is not only the criteria for chromosomal variation of miRNA gene counts. Surprisingly, there are certain chromosomes observed in all the species where the numbers of miRNA gene are extremely high. In human, Chromosome 1, 19, X, and 2 have high number of miRNA genes and constitute about 29% of total miRNA genes. Previously it was shown that miRNAs were non-randomly distributed in human genome and certain chromosomes had significantly higher number of miRNAs than other chromosomes (Calin et al., ). Similar study showed distribution of cancer-associated miRNAs in mouse genome (Sevignani et al., 2007). But, we have studied miRNA gene counts in chromosomes in 20 different species and observed certain chromosomes are having higher miRNA genes through evolution. For example, X chromosomes have high number of miRNA genes whereas Y chromosomes have least or no miRNA gene(s) in all the species. The reason behind such wide variation of miRNA gene counts in chromosomes is not clearly understood but this has come up as an evolutionary conserved phenomenon. It is reported that miRNAs are observed to be localized in chromosomes as clusters and co-expressed all the members from a single polycistronic transcript (Guo et al., ). The clustered miRNAs are expressed simultaneously and are involved in a particular pathway or a particular kind of cellular function. The number of miRNA genes in a cluster may be as high as 40 and such clustered pattern is maintained through evolution (Guo et al., ). The distribution of such clusters is again uneven and concentrated in to certain chromosomes in various species (Chan et al., ). For example, chromosome X and 14 has 40 and 31 number of clusters of miRNA in human (Chan et al., ). This data is supporting our observation that chromosome X and chromosome 14 have higher number of miRNA genes. Possibly, due to this clustering of miRNA we are getting high concentration of miRNA genes into certain chromosomes in various species.
Further, we evaluated Pearson's correlation coefficient between number of miRNA genes with number of coding and non-coding genes in a particular chromosome. About 85% species under study shows strong positive correlation (r ≥ 0.5) for non-coding genes. One obvious reason is that miRNA is a sub-set of non-coding RNA and a single primary transcript of non-coding RNA may produce more than one miRNA (Saini et al., 2007). We have observed strong positive correlation coefficient (r ≥ 0.5) between number of miRNA genes and coding genes in a chromosome as evident from Table 2. The relation between expressions of intronic, exonic and others miRNA with protein coding genes is very close. During the transcription of coding genes the introns left are used as source of intronic miRNAs. Intronic miRNAs or mirtrons can be spliced and debranched through Drosha-independent pathway to produce mature miRNA as observed in Drosophila (Flynt et al., ), nematodes (Ruby et al., 2007), avians (Glazov et al., ), mammals (Berezikov et al., ; Babiarz et al., ), and also plants (Zhu et al., 2008). Lots of reports showed co-expression of intronic miRNAs with corresponding host genes (protein coding genes) in various species, implicating that intronic miRNA and corresponding host protein coding genes are co-transcribed and both groups are under same regulatory elements such as promoter etc. (Lagos-Quintana et al., ; Lau et al., ; Rodriguez et al., 2004; Baskerville and Bartel, ; Chien et al., ). Therefore, array of mirtrons and host protein coding genes co-localize in chromosomes and we should get strong positive correlation between intronic miRNA-coding gene pair. We do get the correlation coefficient value r = 0.91 between intronic miRNA-mRNA pair in human. Earlier report showed strong positive (r > 0.5) Pearson's correlation coefficient of miRNA-mRNA pair in human using expression data (Wang and Li, 2009). There are reports where intronic miRNAs can also be transcribed using their own promoter regions and TSS (Corcoran et al., ; Wang et al., 2010a,b) although major portion of intronic miRNAs are co-transcribed with host protein coding genes. Similarly, exonic and others (includes 3′UTR, 5′ UTR and combinations of any two from intron/exon/3′UTR/5′UTR) miRNA are produced as a by-product of transcription and post-transcription processing of coding genes. It is to be noted that we do not check the expression of miRNAs or coding/non-coding genes but we have calculated correlation coefficient between miRNA-coding gene pair based on their chromosome-specific location. We have observed strong positive correlation (r > 0.5) when summation of intronic and exonic miRNA exceeds intergenic miRNA while r < 0.5 when intergenic exceeds the total of intronic and exonic miRNA genes in a particular species as evident from Table 3. miRNAs under others category are integral part of protein coding genes and can be expressed using transcription machinery of coding genes. Therefore, intronic, exonic and others miRNA genes co-localize with host protein coding genes through evolution and so there must be relation between expression of these miRNA genes with the corresponding coding genes. We hypothesize that intronic, exonic and others miRNAs are produced as by-product of transcription and post-transcriptional processing of corresponding host protein coding genes and this is an evolutionary conserved phenomenon.
Transcription of intergenic miRNAs and its relation with coding gene expression is poorly understood. Here, we have observed positive correlation (r = 0.47) between intergenic miRNA vs. coding genes in humans. The result indicates that expression of intergenic miRNA is less or no way related with expression of protein coding genes with compared to intronic miRNAs (discussed earlier). Previous studies showed that intergenic miRNAs can have its own promoters for transcription, independent of coding genes (Saini et al., 2007; Chien et al., ). But the promoters and TSS has similarity with that of coding genes (Corcoran et al., ). Sometimes more than one pre-miRNA may be produced from the same pri-miRNA transcript, forms cluster miRNAs (Altuvia et al., ), but how exactly intergenic miRNA transcription is related with that of coding gene is not known clearly. However, RNA polymerase II transcribes miRNA like mRNA and the same transcription factors are involved in miRNA transcription as that of mRNA transcription (Wang et al., 2010a,b).
The chromosomes having high number of miRNA genes are more crucial for several cellular processes including different diseases. We have observed that the miRNAs associated with cancer and cardiovascular disorders are concentrated into certain human chromosomes. An experimental observation showed that miRNAs were in cancer-associated genomic regions (CAGR) or in fragile sites (FRA) in different human chromosomes and certain chromosomes were found rich with cancer-associated miRNAs (Calin et al., ; Rossi et al., 2008). A statistically significant association between chromosomal location of miRNAs and tumor susceptibility loci was reported in mouse model and such miRNAs were unevenly distributed in the mouse chromosomes (Sevignani et al., 2007; Rossi et al., 2008). These findings corroborate our data that miRNAs associated with cancer/cardiovascular disease are unevenly distributed in the human chromosomes and certain chromosomes have relatively higher disease-associated miRNA counts like chromosome 19, 14, X, and 1. In other-words, certain chromosomes are most likely controlling major cellular processes including the disease onset and progression in different species through miRNA. Further, we see that almost 46% and 24% of total intergenic miRNAs are associated with cancer and cardiovascular diseases respectively whereas about 17% and 6% of total intronic miRNAs are involved in cancer and cardiovascular diseases. Although the number of intronic miRNAs is greater than intergenic miRNAs in human, the latter group contributes larger for these two disease onset and progression. One possible reason behind this observation could be as follows. Since the clustered miRNAs are mostly involved in any particular cellular pathway including disease onset/progression etc. (Willimott and Wagner, 2012; Godnic et al., ) and their origin is mostly from intergenic portion (Chien et al., ), and so we get intergenic miRNAs sub-type is predominantly associated with cancer and cardiovascular diseases. Since the intronic miRNAs are mostly by-product of transcription/post-transcription processing of corresponding host protein coding genes, one can look into detailed expression profile of the coding genes which hosts those disease-associated intronic/exonic/others miRNAs for better understanding of onset and progression of such diseases. Furthermore, regulation of biogenesis of miRNAs especially intergenic miRNAs will help us to understand some of the molecular mechanisms of pathogenesis of these two diseases.
Materials and methods
Data collection
Details of chromosomal information
We collected the details of each chromosome like length of the chromosome, number of coding genes and non-coding genes in that chromosome for 20 different species (Figure 1) in metazoans from the karyotype data using Ensembl genome browser (release 74- December 2013), jointly maintained by European Bioinformatics Institute (EBI) and Welcome Trust Sanger Institute (WTSI), UK (Flicek et al., ). We used the following Genome assemblies for each species in this work- ASM23792v1 [S. mansoni], BDGP5 (GCA_000001215.1) [D. melanogaster], WBcel235 (GCA_000002985.3) [C. elegans], KH (GCA_000224145.1) [C. intestinalis], AnoCar2.0 (GCA_000090745.1) [A. carolinensis], Galgal4 (GCA_000002315.2) [G. gallus], taeGut3.2.4 [T. guttata], Zv9 (GCA_000002035.2) [D. rerio], MEDAKA1 [O. latipes], TETRAODON8 [T. nigroviridis], Sscrofa10.2 (GCA_000003025.4) [S. scrofa], OANA5 (GCF_000002275.2) [O. anatinus], BROADO5 (GCF_000002295.2) [M. domestica], GRCm38 (GCA_000001635.4) [M. musculus], Rnor_5.0 (GCA_000001895.3) [R. norvegicus], UMD3.1 (GCA_000003055.3) [B. taurus], MMUL_1 [M. mulatta], CHIMP2.1.4 (GCA_000001515.4) [P. troglodytes], gorGor3.1 (GCA_000151905.1) [G. gorilla], and GRCh37 (GCA_000001405.14) [H. sapiens].
Resource for miRNAs
miRBase (Griffiths-Jones et al., ) is a vastly used online registry of miRNAs of a wide range of species focusing on miRNA nomenclature, annotation and target prediction. We collected information about the chromosomal location of individual miRNA of a particular species from miRBase (Release 20: June 2013). We also retrieved the details of each miRNAs and categorized them into four groups according to their genomic location such as intergenic, intronic, exonic and others. We created tables with columns of chromosome number, chromosome length, number of coding/non-coding genes and number of miRNAs in respective chromosomes for each species. We downloaded the miRNAs associated with cancer and cardiovascular diseases from HMDD v1.0—a human miRNA disease database (Lu et al., ).
Statistical analysis
We calculated Pearson product-moment correlation coefficient between column of miRNAs with the column of chromosome length, coding genes and non-coding genes for each species using Microsoft Office Excel 2003 software.
Statements
Acknowledgments
Atanu Ghorai is thankful to ICMR, (New Delhi) for his senior research fellowship.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
AgudaB. D.KimY.Piper-HunterM. G.FriedmanA.MarshC. B. (2008). MicroRNA regulation of a cancer network: consequences of the feedback loops involving miR-17-92, E2F, and Myc. Proc. Natl. Acad. Sci. U.S.A. 105, 19678–19683. 10.1073/pnas.0811166106
2
AltuviaY.LandgrafP.LithwickG.ElefantN.PfefferS.AravinA.et al. (2005). Clustering and conservation patterns of human microRNAs. Nucleic Acids Res. 33, 2697–2706. 10.1093/nar/gki567
3
BabiarzJ. E.RubyJ. G.WangY.BartelD. P.BlellochR. (2008). Mouse ES cells express endogenous shRNAs, siRNAs, and other Microprocessor-independent, Dicer-dependent small RNAs. Genes Dev. 22, 2773–2785. 10.1101/gad.1705308
4
BaskervilleS.BartelD. P. (2005). Microarray profiling of microRNAs reveals frequent coexpression with neighboring miRNAs and host genes. RNA11, 241–247. 10.1261/rna.7240905
5
BenetatosL.HatzimichaelE.LondinE.VartholomatosG.LoherP.RigoutsosI.et al. (2013). The microRNAs within the DLK1-DIO3 genomic region: involvement in disease pathogenesis. Cell. Mol. Life Sci. 70, 795–814. 10.1007/s00018-012-1080-8
6
BerezikovE.ChungW. J.WillsJ.CuppenE.LaiE. C. (2007). Mammalian Mirtron Genes. Molecular Cell. 28, 328–336. 10.1016/j.molcel.2007.09.028
7
BlenkironC.MiskaE. A. (2007). miRNAs in cancer: approaches, aetiology, diagnostics and therapy. Hum. Mol. Genet. 16, R106–R113. 10.1093/hmg/ddm056
8
CalinG. A.CroceC. M. (2006). MicroRNA signatures in human cancers. Nat. Rev. Cancer6, 857–866. 10.1038/nrc1997
9
CalinG. A.SevignaniC.DumitruC. D.HyslopT.NochE.YendamuriS.et al. (2004). Human microRNA genes are frequently located at fragile sites and genomic regions involved in cancers. Proc. Natl. Acad. Sci. U.S.A. 101, 2999–3004. 10.1073/pnas.0307323101
10
CarringtonJ. C.AmbrosV. (2003). Role of microRNAs in plant and animal development. Science301, 336–338. 10.1126/science.1085242
11
ChanS. P.SlackF. J. (2007). And now introducing mammalian mirtrons. Dev. Cell13, 605–607. 10.1016/j.devcel.2007.10.010
12
ChanW.HoM.LiS.TsaiK.LaiC.HsuC.et al. (2012). MetaMirClust: discovery of miRNA cluster patterns using a data-mining approach. Genomics100, 141–148. 10.1016/j.ygeno.2012.06.007
13
ChangT. C.MendellJ. T. (2007). microRNAs in vertebrate physiology and human disease. Annu. Rev. Genomics Hum. Genet. 8, 215–239. 10.1146/annurev.genom.8.080706.092351
14
ChenC. Z.LiL.LodishH. F.BartelD. P. (2004). Micrornas modulate hematopoietic lineage differentiation. Science303, 83–86. 10.1126/science.1091903
15
ChenJ. F.MandelE. M.ThomsonJ. M.WuQ.CallisT. E.HammondS. M.et al. (2006). The role of microrna-1 and microrna-133 in skeletal muscle proliferation and differentiation. Nat. Genet. 38, 228–233. 10.1038/ng1725
16
ChienC. H.SunY. M.ChangW. C.Chiang-HsiehP. Y.LeeT. Y.TsaiW. C.et al. (2011). Identifying transcriptional start sites of human microRNAs based on high-throughput sequencing data. Nucleic Acids Res. 39, 9345–9356. 10.1093/nar/gkr604
17
ChoudhuryN. R.de Lima AlvesF.de Andrés-AguayoL.GrafT.CáceresJ. F.RappsilberJ.et al. (2013). Tissue-specific control of brain-enriched miR-7 biogenesis. Genes Dev. 27, 24–38. 10.1101/gad.199190.112
18
CiccacciC.Di FuscoD.CacciottiL.MorgantiR.D'AmatoC.GrecoC.et al. (2013). MicroRNA genetic variations: association with type 2 diabetes. Acta Diabetol. 50, 867–872. 10.1007/s00592-013-0469-7
19
CorcoranD. L.PanditK. V.GordonB.BhattacharjeeA.KaminskiN.BenosP. V. (2009). Features of mammalian microRNA promoters emerge from polymerase II chromatin immunoprecipitation data. PLoS ONE4:e5279. 10.1371/journal.pone.0005279
20
CuiQ.YuZ.PurisimaE. O.WangE. (2006). Principles of microRNA regulation of a human cellular signaling network. Mol. Syst. Biol. 2, 46. 10.1038/msb4100089
21
CukK.ZucknickM.HeilJ.MadhavanD.SchottS.TurchinovichA.et al. (2013). Circulating microRNAs in plasma as early detection markers for breast cancer. Int. J. Cancer. 132, 1602–1612. 10.1002/ijc.27799
22
ElgarG.VavouriT. (2008). Tuning in to the signals: noncoding sequence conservation in vertebrate genomes. Trends Genet. 24, 344–352. 10.1016/j.tig.2008.04.005
23
ENCODE Project Consortium, BernsteinB. E.BirneyE.DunhamI.GreenE. D.GunterC.SnyderM. (2012). An integrated encyclopedia of DNA elements in the human genome. Nature489, 57–74. 10.1038/nature11247
24
EndoK.NaitoY.JiX.NakanishiM.NoguchiT.GotoY.et al. (2013). MicroRNA 210 as a biomarker for congestive heart failure. Biol. Pharm. Bull. 36, 48–54. 10.1248/bpb.b12-00578
25
FabianM. R.SonenbergN.FilipowiczW. (2010). Regulation of mRNA translation and stability by microRNAs. Annu. Rev. Biochem. 79, 351–379. 10.1146/annurev-biochem-060308-103103
26
FassinaA.MarinoF.SiriM.ZambelloR.VenturaL.FassanM.et al. (2012). The miR-17-92 microRNA cluster: a novel diagnostic tool in large B-cell malignancies. Lab Invest. 92, 1574–1582. 10.1038/labinvest.2012.129
27
FilipowiczW.BhattacharyyaS. N.SonenbergN. (2008). Mechanisms of post-transcriptional regulation by microRNAs: are the answers in sight?Nat. Rev. Genet. 9, 102–114. 10.1038/nrg2290
28
FlicekP.AmodeM. R.BarrellD.BealK.BrentS.Carvalho-SilvaD.et al. (2012). Ensembl 2012. Nucleic Acids Res. 40, D84–D90. 10.1093/nar/gkr991
29
FlyntA. S.GreimannJ. C.ChungW. J.LimaC. D.LaiE. C. (2010). MicroRNA biogenesis via splicing and exosome-mediated trimming in Drosophila. Mol Cell. 38, 900–907. 10.1016/j.molcel.2010.06.014
30
GlazovE. A.CotteeP. A.BarrisW. C.MooreR. J.DalrympleB. P.TizardM. L. (2008). A microRNA catalog of the developing chicken embryo identified by a deep sequencing approach. Genome Res. 18, 957–964. 10.1101/gr.074740.107
31
GodnicI.ZorcM.SkokD. J.CalinG. A.HorvatS.DovcP.et al. (2013). Genome-wide and species-Wide In Silico screening for intragenic microRNAs in human, mouse and chicken. PLoS ONE8:e65165. 10.1371/journal.pone.0065165
32
Griffiths-JonesS.SainiH. K.van DongenS.EnrightA. J. (2008). miRBase: tools for microRNA genomics. Nucleic Acids Res. 36, D154–D158. 10.1093/nar/gkm952
33
GuoL.ZhaoY.ZhangH.YangS.ChenF. (2014). Integrated evolutionary analysis of human miRNA gene clusters and families implicates evolutionary relationships. Gene534, 24–32. 10.1016/j.gene.2013.10.037
34
HavensM. A.ReichA. A.DuelliD. M.HastingsM. L. (2012). Biogenesis of mammalian microRNAs by a non-canonical processing pathway. Nucleic Acids Res. 40, 4626–4640. 10.1093/nar/gks026
35
HébertS. S.HorréK.NicolaïL.PapadopoulouA. S.MandemakersW.SilahtarogluA. N.et al. (2008). Loss of microRNA cluster mir-29a/b-1 in sporadic Alzheimer's disease correlates with increased BACE1/β-secretase expression. Proc. Natl. Acad. Sci. U.S.A. 105, 6415–6420. 10.1073/pnas.0710263105
36
HwangH. W.MendellJ. T. (2007). MicroRNAs in cell proliferation, cell death, and tumorigenesis. Br. J. Cancer96, R40–R44. 10.1038/sj.bjc.6603023
37
IveyK. N.MuthA.ArnoldJ.KingF. W.YehR. F.FishJ. E.et al. (2008). Microrna regulation of cell lineages in mouse and human embryonic stem cells. Cell Stem Cell2, 219–229. 10.1016/j.stem.2008.01.016
38
JonesC. I.ZabolotskayaM. V.KingA. J.StewartH. J.HorneG. A.ChevassutT. J.et al. (2012). Identification of circulating microRNAs as diagnostic biomarkers for use in multiple myeloma. Br. J. Cancer107, 1987–1996. 10.1038/bjc.2012.525
39
KellerA.LeidingerP.GislefossR.HaugenA.LangsethH.StaehlerP.et al. (2011). Stable serum miRNA profiles as potential tool for non-invasive lung cancer diagnosis. RNA Biol. 8, 506–516. 10.4161/rna.8.3.14994
40
KimV. N. (2005). MicroRNA biogenesis: coordinated cropping and dicing. Nat. Rev. Mol. Cell. Biol. 6, 376–385. 10.1038/nrm1644
41
KumamotoK.SpillareE. A.FujitaK.HorikawaI.YamashitaT.AppellaE.et al. (2008). Nutlin-3a activates p53 to both down-regulate inhibitor of growth 2 and up-regulate mir-34a, mir-34b, and mir-34c expression, and induce senescence. Cancer Res. 68, 3193–3203. 10.1158/0008-5472.CAN-07-2780
42
Lagos-QuintanaM.RauhutR.LendeckelW.TuschlT. (2001). Identification of novel genes coding for small expressed RNAs. Science294, 853–858. 10.1126/science.1064921
43
LauN. C.LimL. P.WeinsteinE. G.BartelD. P. (2001). An abundant class of tiny RNAs with probable regulatory roles in Caenorhabditis elegans. Science294, 858–862. 10.1126/science.1065062
44
LuM.ZhangQ.DengM.MiaoJ.GuoY.GaoW.et al. (2008). An analysis of human microRNA and disease associations. PLoS ONE3:e3420. 10.1371/journal.pone.0003420
45
LukiwW. J.PogueA. I. (2007). Induction of specific micro rna (mirna) species by ros-generating metal sulfates in primary human brain cells. J. Inorg. Biochem. 101, 1265–1269. 10.1016/j.jinorgbio.2007.06.004
46
MaesO. C.AnJ.SarojiniH.WangE. (2008). Murine micrornas implicated in liver functions and aging process. Mech. Ageing Dev. 129, 534–541. 10.1016/j.mad.2008.05.004
47
NasserM. W.DattaJ.NuovoG.KutayH.MotiwalaT.MajumderS.et al. (2008). Down-regulation of micro-RNA-1 (miR-1) in Lung Cancer. Suppression of tumorigenic property of lung cancer cells and their sensitization to doxorubicin-induced apoptosis by miR-1. J. Biol. Chem. 283, 33394–33405. 10.1074/jbc.M804788200
48
NilsenT. W. (2007). Mechanisms of microRNA-mediated gene regulation in animal cells. Trends Genet. 23, 243–249. 10.1016/j.tig.2007.02.011
49
OkamuraK.HagenJ. W.DuanH.TylerD. M.LaiE. C. (2007). The mirtron pathway generates microRNA-class regulatory RNAs in drosophila. Cell130, 89–100. 10.1016/j.cell.2007.06.028
50
PichiorriF.SuhS. S.RocciA.De LucaL.TaccioliC.SanthanamR.et al. (2010). Downregulation of p53-inducible microRNAs 192, 194, and 215 impairs the p53/MDM2 autoregulatory loop in multiple myeloma development. Cancer Cell18, 367–381. 10.1016/j.ccr.2010.09.005
51
PlaisanceV.AbderrahmaniA.Perret–MenoudV.JacqueminP.LemaigreF.RegazziR. (2006). Microrna-9 controls the expression of Granuphilin/Slp4 and the secretory response of insulinproducing cells. J. Biol. Chem. 281, 26932–26942. 10.1074/jbc.M601225200
52
RodriguezA.Griffiths-JonesS.AshurstJ. L.BradleyA. (2004). Identification of mammalian microRNA host genes and transcription units. Genome Res. 14, 1902–1910. 10.1101/gr.2722704
53
RodriguezA.VigoritoE.ClareS.WarrenM. V.CouttetP.SoondD. R.et al. (2007). Requirement of Bic/microrna-155 for normal immune function. Science316, 608–611. 10.1126/science.1139253
54
RossiS.SevignaniC.NnadiS. C.SiracusaL. D.CalinG. A. (2008). Cancer-associated genomic regions (CAGRs) and noncoding RNAs: bioinformatics and therapeutic implications. Mamm Genome19, 526–540. 10.1007/s00335-008-9119-8
55
RubyJ. G.JanC. H.BartelD. P. (2007). Intronic microRNA precursors that bypass Drosha processing. Nature448, 83–86. 10.1038/nature05983
56
SainiH. K.Griffiths-JonesS.EnrightA. J. (2007). Genomic analysis of human microRNA transcripts. Proc. Natl. Acad. Sci. U.S.A. 104, 17719–17724. 10.1073/pnas.0703890104
57
SalviA.AbeniE.PortolaniN.BarlatiS.De PetroG. (2013). Human hepatocellular carcinoma cell-specific miRNAs reveal the differential expression of miR-24 and miR-27a in cirrhotic/non-cirrhotic HCC. Int. J. Oncol. 42, 391–402. 10.3892/ijo.2012.1716
58
SaweraM.GorodkinJ.CireraS.FredholmM. (2005). Mapping and expression studies of the mir17-92 cluster on pig Chromosome 11. Mamm. Genome16, 594–598. 10.1007/s00335-005-0013-3
59
SayedD.HongC.ChenI. Y.LypowyJ.AbdellatifM. (2007). MicroRNAs play an essential role in the development of cardiac hypertrophy. Circ. Res. 100, 416–424. 10.1161/01.RES.0000257913.42552.23
60
SevignaniC.CalinG. A.NnadiS. C.ShimizuM.DavuluriR. V.HyslopT.et al. (2007). MicroRNA genes are frequently located near mouse cancer susceptibility loci. Proc. Natl. Acad. Sci. U.S.A. 104, 8017–8022. 10.1073/pnas.0702177104
61
StefaniG.SlackF. J. (2008). Small non-coding RNAs in animal development. Nat. Rev. Mol. Cell. Biol. 9, 219–230. 10.1038/nrm2347
62
TaganovK. D.BoldinM. P.BaltimoreD. (2007). MicroRNAs and immunity: tiny players in a big field. Immunity26, 133–137. 10.1016/j.immuni.2007.02.005
63
TarasovV.JungP.VerdoodtB.LodyginD.EpanchintsevA.MenssenA.et al. (2007). Differential regulation of micrornas by p53 revealed by massively parallel sequencing: mir-34a is a p53 target that induces apoptosis and G1-arrest. Cell Cycle6, 1586–1593. 10.4161/cc.6.13.4436
64
van de BuntM.GaultonK. J.PartsL.MoranI.JohnsonP. R.LindgrenC. M.et al. (2013). The miRNA profile of human pancreatic islets and beta-cells and relationship to type 2 diabetes pathogenesis. PLoS ONE8:e55272. 10.1371/journal.pone.0055272
65
WangG.van der WaltJ. M.MayhewG.LiY. J.ZüchnerS.ScottW. K.et al. (2008a). Variation in the miRNA-433 binding site of FGF20 confers risk for Parkinson disease by overexpression of alpha-synuclein. Am. J. Hum. Genet. 82, 283–289. 10.1016/j.ajhg.2007.09.021
66
WangG.WangY.ShenC.HuangY.HuangK.HuangT. H. M.et al. (2010a). RNA polymerase II binding patterns reveal genomic regions involved in microRNA gene regulation. PLoS ONE5:e13798. 10.1371/journal.pone.0013798
67
WangJ.LuM.QiuC.CuiQ. (2010b). TransmiR: a transcription factor-microRNA regulation database. Nucleic Acids Res. 38, D119–D122. 10.1093/nar/gkp803
68
WangW. X.RajeevB. W.StrombergA. J.RenN.TangG.HuangQ.et al. (2008b). The expression of microRNA mir-107 decreases early in Alzheimer's disease and may accelerate disease progression through regulation of beta-site amyloid precursor protein-cleaving enzyme 1. J. Neurosci. 28, 1213–1223. 10.1523/JNEUROSCI.5065-07.2008
69
WangY. P.LiK. B. (2009). Correlation of expression profiles between microRNAs and mRNA targets using NCI-60 data. BMC Genomics10:218. 10.1186/1471-2164-10-218
70
WeilandM.GaoX. H.ZhouL.MiQ. S. (2012). Small RNAs have a large impact: circulating microRNAs as biomarkers for human diseases. RNA Biol. 9, 850–859. 10.4161/rna.20378
71
WestholmJ. O.LaiE. C. (2011). Mirtrons: microRNA biogenesis via splicing. Biochimie. 93, 1897–1904. 10.1016/j.biochi.2011.06.017
72
WillimottS.WagnerS. D. (2012). Stromal cells and CD40 ligand (CD154) alter the miRNome and induce miRNA clusters including, miR-125b/miR-99a/let-7c and miR-17-92 in chronic lymphocytic leukaemia. Leukemia26, 1113–1116. 10.1038/leu.2011.299
73
WoutersM. D.van GentD. C.HoeijmakersJ. H.PothofJ. (2011). MicroRNAs, the DNA damage response and cancer. Mutat. Res. 717, 54–66. 10.1016/j.mrfmmm.2011.03.012
74
WuchtyS.ArjonaD.BozdagS.BauerP. O. (2012). Involvement of microRNA families in cancer. Nucleic Acids Res. 40, 8219–8226. 10.1093/nar/gks627
75
YingS. Y.ChangC. P.LinS. L. (2010). Intron-mediated RNA interference, intronic microRNAs, and applications. Methods Mol. Biol. 629, 205–237. 10.1007/978-1-60761-657-3_14
76
ZhangB.WangQ.PanX. (2007). MicroRNAs and their regulatory roles in animals and plants. J. Cell Physiol. 210, 279–289. 10.1002/jcp.20869
77
ZhaoY.RansomJ. F.LiA.VedanthamV.von DrehleM.MuthA. N.et al. (2007). Dysregulation of cardiogenesis, cardiac conduction, and cell cycle in mice lacking mirna-1–2. Cell129, 303–317. 10.1016/j.cell.2007.03.030
78
ZhaoY.SamalE.SrivastavaD. (2005). Serum response factor regulates a muscle-specific microRNA that targets Hand2 during cardiogenesis. Nature436, 214–220. 10.1038/nature03817
79
ZhuQ. H.SpriggsA.MatthewL.FanL.KennedyG.GublerF.et al. (2008). A diverse set of microRNAs and microRNA-like small RNAs in developing rice grains. Genome Res. 18, 1456–1465. 10.1101/gr.075572.107
Summary
Keywords
non-coding gene, miRNA, chromosome, intronic, intergenic, correlation coefficient, coding genes, disease-association
Citation
Ghorai A and Ghosh U (2014) miRNA gene counts in chromosomes vary widely in a species and biogenesis of miRNA largely depends on transcription or post-transcriptional processing of coding genes. Front. Genet. 5:100. doi: 10.3389/fgene.2014.00100
Received
08 October 2013
Accepted
07 April 2014
Published
29 April 2014
Volume
5 - 2014
Edited by
Ryan Michael O'Connell, University of Utah, USA
Reviewed by
Vladimir Benes, European Molecular Biology Laboratory, Germany; Alessio Paone, Sapienza University of Rome, Italy
Copyright
© 2014 Ghorai and Ghosh.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Utpal Ghosh, Department of Biochemistry and Biophysics, University of Kalyani, Kalyani 741235, India e-mail: utpal8ghosh@gmail.com; utpal8ghosh@klyuniv.ac.in
This article was submitted to Non-Coding RNA, a section of the journal Frontiers in Genetics.
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.