Novel microRNA families expanded in the human genome
© Du et al.; licensee BioMed Central Ltd. 2013
Received: 4 June 2012
Accepted: 31 January 2013
Published: 12 February 2013
Skip to main content
© Du et al.; licensee BioMed Central Ltd. 2013
Received: 4 June 2012
Accepted: 31 January 2013
Published: 12 February 2013
Most studies on the origin and evolution of microRNA in the human genome have been focused on its relationship with repetitive elements and segmental duplications. However, duplication events at a smaller scale (<1 kb) could also contribute to microRNA expansion, as demonstrated in this study.
Using comparative genome analysis and bioinformatics methods, we found nine novel expanded microRNA families enriched in short duplicated sequences in the human genome. Furthermore, novel genomic regions were found to contain microRNA paralogs for microRNA families previously analyzed to be related to segmental duplications. We found that for microRNA families expanded in the human genome, 14 families are specific to the primate lineage, and nine are non-specific, respectively. Two microRNA families (hsa-mir-1233 and hsa-mir-622) appear to be further expanded in the human genome, and were confirmed by fluorescence in situ hybridization. These novel microRNA families expanded in the human genome were mostly embedded in or close to proteins with conserved functions. Furthermore, besides the Alu element, L1 elements could also contribute to the origination of microRNA paralog families.
Together, we found that small duplication events could also contribute to microRNA expansion, which could provide us novel insights on the evolution of human genome structure and function.
Although only ~22 nucleotides long, microRNA are vital to the developmental process of animals and plants through post-transcriptional gene regulation. The efficiency of microRNA transcription and stability is regulated in a tissue-specific manner [1–4]. After processing from precursor microRNA molecules, single-stranded microRNA can bind to 3'-UTR of messenger RNA via its seed region, in turn affecting mRNA translation or stability [1, 5]. When aberrantly expressed, microRNA can lead to the progression of certain diseases [6, 7]. They also participate in pathogen-host interaction , and exert systemic effects through intercellular trafficking . Due to its significant roles in a wide variety of biological processes, genome evolution of microRNA structure and function have been studied extensively [1, 2, 10–18].
Both duplicated sequence fragments and repetitive elements in the genome could contribute to the expansion of microRNA families [10–18]. Species-specific microRNA expansion could have functional importance, and it has been demonstrated recently that an expanded microRNA cluster is fundamental to the maintenance of embryonic stem cell development in mice . However, most duplication events studied for microRNA evolution are segmental or tandem duplications [15–17]. Segmental duplications are tremendously important in elucidating the evolution of protein-coding genes, and their roles in providing novel genomic mechanisms to cope with selection pressure and environmental changes [19–24]. However, segmental duplications are defined as > 1 kb and > 90% identity, and computational methods are designed accordingly and limited in studying only these sequences. Thus, small-scale duplication events are not covered, but which could also play important roles in the evolution of small molecules, especially for microRNA. The relevance and importance of these small duplication events and their relationship with microRNA evolution has not been reported.
Here, we used a systematic approach, and found that small duplication events in the human genome could also contribute to microRNA expansion. In total, nine novel microRNA families were found to be expanded in the human genome, and additional new genomic regions were discovered to be related to the expansion of microRNA families reported previously. We found that novel microRNA family expanded in the human genome are located close to proteins with conserved function, and confirmed two of the microRNA expansion events by fluorescence in situ hybridization. These results could render us novel insights on the evolution of human genome structure and function.
To explore the relationship between short duplicated genomic sequences and the expansion of certain microRNA families in the human genome, we used repeat-masked reference genome for pairwise comparison (see Methods), since microRNA associated with repetitive elements have been well characterized [16–18]. We focused on detecting microRNA paralogs enriched in short duplication events (<1 kb) in the human genome, which were overlooked previously by studies on the role of segmental duplications (emerged recently, >1 kb and > 90% identity) to the origin and evolution of microRNA families [10–17].
In total, nine novel microRNA families were found to be enriched in short duplicated fragments (Additional file 1). We detected 26 microRNA families, previously found to be related to segmental duplications , in which new microRNA paralogs could also be enriched in short genomic fragments (Additional file 2).
Furthermore, our methods were able to detect 4 mitochondria-related microRNA families (Additional file 3), and 20 other microRNA families already deposited in miRBase (Additional file 4). However, three of the four mitochondria-related microRNA families, hsa-mir-1974, hsa-mir-1977, and hsa-mir-1978, were not further curated in miRBase, since they overlap with transfer-RNA (t-RNA) sequences. Despite the exclusion of these three microRNA in the miRBase, experimental evidence suggests that microRNA can be derived from t-RNA, and still demonstrate biological functions [25, 26].
We examined the location of these microRNA paralogs relative to protein-coding genes, the sequence conservation rate, and the flanking genomic sequences as well, to see if there exist certain patterns related to the expansion of these detected microRNA paralog families.
Direct analyses on microRNA paralog sequence features in all expanded microRNA families found that mature microRNA sequences and the seed regions of microRNA paralogs could be under the influence of different selective forces . For instance, in two microRNA families, hsa-mir-1233 and hsa-mir-1244 (Figure 3, Additional file 7), we can clearly see that most of the microRNA paralogs had nearly identical mature sequences, while others demonstrate dissimilar features. These sequence changes could affect the stability of the stem-loop structure of the microRNA paralogs. Further examination on the flanking genomic regions of these microRNA families found they are highly conserved (Figure 3, Additional file 7).
Furthermore, the analysis on repetitive elements in the flanking genomic regions may contribute to the understanding of the origin of these microRNA families. Repetitive elements, such as Alu elements, were discovered previously to potentially affect the expansion of certain microRNA clusters as well as the genes located in duplicated genomic segments [11, 30]. We found that Alu elements are also enriched in the flanking regions of microRNA paralogs (ANOVA) (Additional file 8). Surprisingly, the L1 elements have also been enriched in these sequences, which implied that they could also function in the evolutionary path of formation of certain microRNA paralogs. Furthermore, L1 elements had the highest percentage (21%), but not significantly different with AluS elements (t-test, P > 0.05).
Using information provided by short duplication events in the human genome, we discovered novel microRNA paralogs for different microRNA families, which could have fundamental importance for human biology and the study of complex phenotypic traits. However, since the number of microRNAs deposited in the public database are far less than predicted, new microRNA paralogs and clusters may still exist in the human genome. Furthermore, this method can be extended for other non-coding RNAs, which can potentially reveal interesting patterns of specific important duplications in the human genome.
The study of microRNA from a comparative genomics standpoint is not without some limitations. The majority of the computational prediction methods use either precursor or mature microRNA sequences, combined with the examination of predicted secondary-structure . The drawback of this approach is that only highly conserved potential microRNAs within or across species can be discovered. Experimental methods to date have only sampled limited tissues at few time-points, which will restrict the identification of novel microRNA . Even though high-throughput small RNA sequencing has led to the identification of numerous novel microRNAs, the number and speed of discovery are still limited . Our methods could to some degree help the identification and interpretation of novel microRNA paralogs in the human genome.
Alu elements have been proposed to be involved in the propagation of microRNA clusters and segmental duplications [11, 30]. We also found in this study that AluS is the most abundant subfamily (17%), followed by AluJ (8%) and AluY (4%) (Additional file 8), which indicates that microRNA expansion events discovered here could follow a similar evolutionary path as reported in . Furthermore, we provide evidence here that in addition to Alu elements, L1 repetitive elements are also associated with duplication of specific microRNA sequences. While numerous molecular mechanisms could be hypothesized for microRNA expansion, there remains a strong association that specific repetitive elements may have played important roles in the evolution of microRNA containing loci in humans. These duplicated paralogs could obtain partial and/or novel functions or exert dosage effects, acquire novel regulatory elements for tissue-specific expression, possess modifications in seed sequence resulting in novel microRNAs, and potentially beneficial or detrimental to certain individuals or a specific population . The conserved flanking sequences indicate that similar regulatory mechanisms could be involved in the tuning of these paralogs, but the changed sequences of the mature microRNA paralog could target different sets of genes and affect different gene clusters. However, detailed functional analysis is still needed to understand the outcome of the duplicated microRNA molecules .
We found that novel microRNA could potentially be generated through the processing of pseudogene transcripts. These expanded microRNA families in the human genome could provide novel insights on the evolution and regulation of the human genome and its interaction with the environment. Particularly, these expanded microRNA families could add another level of complexity into the regulatory network formed among messenger RNAs, transcribed pseudogenes, long noncoding RNAs, through the amplification of microRNA abundance .
Limited amounts of population resequencing data restricted our understanding of intra-species variation regarding one of the structural variations, the expansion of microRNA family, which has important practical implications (e.g. disease genetics). However, the rapid developments in sequencing technology will alleviate this problem in the near future, as already shown in the 1000 genomes project [35–37]. Further detailed analysis on the association of microRNA expansion with features unique to each ethnic group will potentially reveal their biological importance in human diversity.
Taken together, we found that small duplication events in the human genome may contribute to microRNA expansion, which could provide novel insights on the evolution of genome architecture, in addition to the development of human diseases.
The whole procedure for the analysis can be found in Additional file 9. The human genome and 25 other animal genome sequences with repetitive elements masked were downloaded from Ensembl (Additional file 10). We searched for genomic duplications by Megablast using default parameters . All duplicated fragments were kept without restriction on length and identity of the sequence alignments, which is different from the detection methods for segmental duplications (length >1 kb and identity >90%) .
Coordinates of human microRNAs were retrieved from miRBase (GRCh37) , and used to enrich for those duplicated genomic fragments overlapping with microRNAs using in-house Perl scripts. The extracted final coordinates were compared again to those of microRNAs deposited in miRBase, to find potential novel duplicated microRNAs, as well as their genomic locations with regard to exonic, intronic or untranslated regions of known genes.
To search for paralogs and orthologs of microRNAs duplicated in the human genome, we retrieved the regional genomic sequences with repetitive sequences masked surrounding a representative microRNA deposited in miRBase, 3,000 base pairs (kb) upstream and downstream. The sequences were used to compare to the human genome sequence using Megablast, and an iterative procedure was used to enrich for genomic sequences related to microRNA paralogs. Furthermore, the retrieved genomic sequences were used to detect orthologs in other species using Megablast, including nine primates (Pan troglodytes, Pongo pygmaeus, Gorilla gorilla, Macaca mulatta, Callithrix jacchus, Tarsius syrichta, Tupaia belangeri, Microcebus murinus, and Otolemur garnettii), two rodents (Rattus norvegicus and Mus musculus), and three domestic livestock species (Equus caballus, Bos taurus and Sus scrofa) (Additional file 10). cDNA evidence was also searched in NCBI by using the retrieved genomic sequences containing microRNAs detected to be expanded in the human genome.
To detect human microRNAs in miRBase overlapping with the repetitive elements, we retrieved the genomic positions of the repetitive elements in the human genome from the UCSC genome table browser (group: variation and repeats, track: RepeatMasker) , and compared results to the coordinates of human microRNAs in miRBase. To further explore the relationship of repetitive elements with the origin of duplicated microRNA paralogs, the coordinates of repeats distributed in the regions of microRNA paralogs (upstream and downstream 5 kb regions) were compared and examined. Fractions of repetitive elements to the selected genomic sequences were calculated by their nucleotide lengths, and enrichment test was performed using Chi-square test, by comparing to the fraction of repetitive elements in the whole human genome. Secondary structure of the microRNA precursor was predicted using RNAfold .
To validate the results obtained by computational prediction, we selected two probes for the two microRNA families (hsa-mir-1233 and hsa-mir-622), and labeled them with TAMRA (red) and FITC (green), respectively (Additional file 11). FISH was performed on a normal primary neonatal dermal fibroblast cell line of European origin (PCS-201-010 from ATCC® Primary Cell Solutions™) following standard procedures (Creative™ Biolabs) . Due to the short probe length (~1 kb), we optimized the hybridization condition several times, until consistent results were obtained.
This project was supported in part by National Research Initiative Competitive Grant no. 2008-35205-05309 and 2008-35205-18712 from the USDA National Institute of Food and Agriculture. We also appreciate the financial support provided by the State of Iowa and Hatch Funding.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.