Differences in molecular evolutionary rates among microRNAs in the human and chimpanzee genomes.

BACKGROUND: The rise of the primate lineage is accompanied by an outstanding emergence of microRNAs, small non-coding RNAs with a prominent role in gene regulation. In spite of their biological importance little is known about the way in which natural selection has influenced microRNAs in the human lineage. To study the recent evolutionary history of human microRNAs and to analyze the signatures of natural selection in genomic regions harbouring microRNAs we have investigated the nucleotide substitution rates of 1,872 human microRNAs in the human and chimpanzee lineages. RESULTS: We produced a depurated set of microRNA alignments of human, chimpanzee and orang-utan orthologs combining BLAT and liftOver and selected 1,214 microRNA precursors presenting optimal secondary structures. We classified microRNAs in categories depending on their genomic organization, duplication status and conservation along evolution. We compared substitution rates of the aligned microRNAs between human and chimpanzee using Tajima's Relative Rate Test taking orang-utan as out-group and found several microRNAs with particularly high substitution rates in either the human or chimpanzee branches. We fitted different models of natural selection on these orthologous microRNA alignments and compared them using a likelihood ratio test that uses ancestral repeats and microRNA flanking regions as neutral sequences. We found that although a large fraction of human microRNAs is highly conserved among the three species studied, significant differences in rates of molecular evolution exist among microRNA categories. Particularly, primate-specific microRNAs, which are enriched in isolated and single copy microRNAs, more than doubled substitution rates of those belonging to older, non primate-specific microRNA families. CONCLUSIONS: Our results corroborate the remarkable conservation of microRNAs, a proxy of their functional relevance, and indicate that a subset of human microRNAs undergo nucleotide substitutions at higher rates, which may be suggestive of the action of positive selection.

An ancestral miR-1304 allele present in Neanderthals regulates genes involved in enamel formation and could explain dental differences with modern humans.

Genetic changes in regulatory elements are likely to result in phenotypic effects that might explain population-specific as well as species-specific traits. MicroRNAs (miRNAs) are posttranscriptional repressors involved in the control of almost every biological process. These small noncoding RNAs are present in various phylogenetic groups, and a large number of them remain highly conserved at the sequence level. MicroRNA-mediated regulation depends on perfect matching between the seven nucleotides of its seed region and the target sequence usually located at the 3' untranslated region of the regulated gene. Hence, even single changes in seed regions are predicted to be deleterious as they may affect miRNA target specificity. In accordance to this, purifying selection has strongly acted on these regions. Comparison between the genomes of present-day humans from various populations, Neanderthal, and other nonhuman primates showed an miRNA, miR-1304, that carries a polymorphism on its seed region. The ancestral allele is found in Neanderthal, nonhuman primates, at low frequency (~5%) in modern Asian populations and rarely in Africans. Using miRNA target site prediction algorithms, we found that the derived allele increases the number of putative target genes for the derived miRNA more than ten-fold, indicating an important functional evolution for miR-1304. Analysis of the predicted targets for derived miR-1304 indicates an association with behavior and nervous system development and function. Two of the predicted target genes for the ancestral miR-1304 allele are important genes for teeth formation, enamelin, and amelotin. MicroRNA overexpression experiments using a luciferase-based assay showed that the ancestral version of miR-1304 reduces the enamelin- and amelotin-associated reporter gene expression by 50%, whereas the derived miR-1304 does not have any effect. Deletion of the corresponding target sites for miR-1304 in these dental genes avoided their repression, which further supports their regulation by the ancestral miR-1304. Morphological studies described several differences in the dentition of Neanderthals and present-day humans like slower dentition timing and thicker enamel for present-day humans. The observed miR-1304-mediated regulation of enamelin and amelotin could at least partially underlie these differences between the two Homo species as well as other still-unraveled phenotypic differences among modern human populations.

Y-chromosome diversity in Native Mexicans reveals continental transition of genetic structure in the Americas.

The genetic characterization of Native Mexicans is important to understand multiethnic based features influencing the medical genetics of present Mexican populations, as well as to the reconstruct the peopling of the Americas. We describe the Y-chromosome genetic diversity of 197 Native Mexicans from 11 populations and 1,044 individuals from 44 Native American populations after combining with publicly available data. We found extensive heterogeneity among Native Mexican populations and ample segregation of Q-M242* (46%) and Q-M3 (54%) haplogroups within Mexico. The northernmost sampled populations falling outside Mesoamerica (Pima and Tarahumara) showed a clear differentiation with respect to the other populations, which is in agreement with previous results from mtDNA lineages. However, our results point toward a complex genetic makeup of Native Mexicans whose maternal and paternal lineages reveal different narratives of their population history, with sex-biased continental contributions and different admixture proportions. At a continental scale, we found that Arctic populations and the northernmost groups from North America cluster together, but we did not find a clear differentiation within Mesoamerica and the rest of the continent, which coupled with the fact that the majority of individuals from Central and South American samples are restricted to the Q-M3 branch, supports the notion that most Native Americans from Mesoamerica southwards are descendants from a single wave of migration. This observation is compatible with the idea that present day Mexico might have constituted an area of transition in the diversification of paternal lineages during the colonization of the Americas.

RNA editing independently occurs at three mir-376a-1 sites and may compromise the stability of the microRNA hairpin.

RNA editing is being recognized as an important post-transcriptional mechanism that may have crucial roles in introducing genetic variation and phenotypic diversity. Despite microRNA editing recurrence, defining its biological relevance is still under extended debate. To better understand microRNA editing function and regulation we performed an exhaustive characterization of the A-to-I site-specific patterns in mir-376a-1, a mammalian microRNA which RNA editing is involved in the regulation of development and in disease. Thorough an integrative approach based on high-throughput small RNA sequencing, Sanger sequencing and computer simulations we explored mir-376a-1 editing in samples from various individuals and primate species including human placenta and macaque, gorilla, chimpanzee and human brain cortex. We observed that mir-376a-1 editing is a common phenomenon in the mature and primary microRNA molecules and it is more frequently detected in brain than in placenta. Primary mir-376a-1 is edited at three positions, -1, +4 and +44. Editing frequency estimations and in silico simulations indicated that editing was not equally recurrent along the three mir-376a-1 sites, nevertheless no epistatic interactions among them were observed. Particularly, the +4 site, located in the seed region of the mature miR-376a-5p, reached the highest editing frequency in all samples. Secondary structure predictions revealed that the +4 position was the one that conferred the highest stability to the mir-376a-1 hairpin. We suggest that molecular stability might partially explain the editing recurrence observed in certain microRNAs and that editing events conferring new functional regulatory roles in particular tissues and species could have been conserved along evolution, as it might be the case of mir-376a-1 in primate brain cortex.