Table of Contents Author Guidelines Submit a Manuscript
International Journal of Genomics
Volume 2018 (2018), Article ID 2456076, 4 pages
https://doi.org/10.1155/2018/2456076
Research Article

Clustering Pattern and Functional Effect of SNPs in Human miRNA Seed Regions

1Bioinformatics Section, School of Basic Medical Sciences, Southern Medical University, Guangzhou, Guangdong 510515, China
2Department of Medical Genetics, School of Basic Medical Sciences, Southern Medical University, Guangzhou, Guangdong 510515, China
3Guangdong Key Laboratory of Biochip Technology, Southern Medical University, Guangzhou, Guangdong 510515, China
4Guangdong Province Key Laboratory of Psychiatric Disorders, Southern Medical University, Guangzhou, Guangdong 510515, China

Correspondence should be addressed to Jian Zhang; moc.liamxof@ums.gnahznaij

Received 5 December 2017; Accepted 7 February 2018; Published 6 March 2018

Academic Editor: Giulia Piaggio

Copyright © 2018 Sha He et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

miRNAs are a class of noncoding RNAs important in posttranscriptional repressors and involved in the regulation of almost every biological process by base paring with target genes through sequence in their seed regions. Genetic variations in the seed regions have vital effects on gene expression, phenotypic variation, and disease susceptibility in humans. The distribution pattern of genetic variation in miRNA seed regions might be related to miRNA function and is worth paying more attention to. We here employed computational analyses to explore the clustering pattern and functional effect of SNPs in human miRNA seed regions. A total of 1879 SNPs were mapped to 1226 human miRNA seed regions. We found that miRNAs with SNPs in their seed region are significantly enriched in miRNA clusters. We also found that SNPs in clustered miRNA seed regions have a lower functional effect than have SNPs in nonclustered miRNA seed regions. Additionally, we found that clustered miRNAs with SNPs in seed regions are involved in more pathways. Overall, our results demonstrate that SNPs in clustered miRNA seed regions can take part in more intricate and complex gene-regulating networks with lower functional cost by functional complementarity. Moreover, our results also broaden current knowledge on the genetic variation in human miRNA seed regions.

1. Introduction

miRNAs are small noncoding RNAs of 20–22 nucleotides (nt) in length in their mature form, processed from a longer hairpin structure, that act as posttranscriptional gene regulators by either mRNA degradation or translational repression [1]. It is estimated that miRNAs regulate more than 30~60% of all protein-coding genes, thus building complex regulatory networks that participate in the control of most biological processes and are related to shaping phenotypic variability and disease development [2, 3]. miRNA-mediated gene regulation depends on perfect matching between the seven nucleotides of its seed region (nucleotides 2 through 8 from the 5 end of the mature miRNAs) and the target sequence usually located at the 3 untranslated regions (UTRs) of the regulated genes [4, 5]. With the rapid advance of genome sequencing technologies in recent years, many genetic variants have been identified in miRNA seed regions [69]. These genetic variants may disrupt the interactions between miRNAs and their targets or create new targets, rewiring the miRNA regulatory networks and causing diseases. Moreover, previous studies have proven that genetic variations in the seed regions have profound and broad effects on gene expression, phenotypic variation, and disease susceptibility in humans [1013]. However, there is still little information on the distribution patterns and allele-dependent functional effect of the genetic variations in human miRNA seed regions.

Wang et al. [14] reported that the clustered miRNAs tend to be evolutionarily conserved and miRNAs in the same clusters tend to be coexpressed and regulate overlapping sets of target genes. Therefore, we hypothesized that clustered miRNAs tend to be more tolerant of genetic variations in the seed region due to functional complementary than are those nonclustered miRNAs. In order to test this hypothesis, we investigated the clustering patterns and functional effect of SNPs in human miRNA seed regions in this study. Our results demonstrate that SNPs in clustered miRNA seed regions can take part in more intricate gene-regulating networks with lower functional cost by functional complementarity.

2. Materials and Methods

2.1. Genome-Wide Identification of SNPs in Human miRNA Seed Regions

Genomic coordinates of human mature miRNAs were downloaded and extracted from miRBase version 21 (http://www.mirbase.org/). This release of miRBase comprises 1881 precursor sequences and 2813 mature miRNAs (2588 unique mature miRNA sequences). Genomic locations of miRNA seeds were determined from genomic locations of the 2nd and 8th bases of mature miRNAs. The SNP data for human assemblies GRCh38/hg38 were downloaded from NCBI dbSNP database (version 147, https://www.ncbi.nlm.nih.gov/SNP/). For each miRNA, we collected all SNPs in the seed regions from dbSNP using in-house Perl scripts.

2.2. Computational Predictions of Allele-Dependent miRNA Target and Percent Overlap

The miRNA sequences were downloaded from the miRBase (version 21). To determine the target gene spectrum for the reference and derived allele of the SNPs in the miRNA seed regions, we predicted targets on the human genome assembly (NCBI36/hg18, March 2006) using the online target prediction program, namely, TargetRank (http://hollywood.mit.edu/targetrank) [15]. The reference targets for a SNP were predicted using TargetRank to analyse the miRNA seed sequence carrying the reference allele of the SNP. On the other hand, the derived targets of a SNP were predicted when the derived allele of the SNP was in the miRNA seed sequence. We defined the common targets (overlapping genes) as predicted targets for both the reference and derived alleles of a SNP. Percent overlap between the reference and derived targets for a SNP in miRNA seeds was determined using cosine similarity [16], which is calculated by the total number of overlapping genes divided by the square root of the product of the number of targets of reference and derived alleles of a SNP. Taking the square root of the number of predicted targets reduces the influence of miRNAs with abnormally large numbers of targets and simultaneously normalizes the result, generating a score between 0 and 1.

2.3. Clustering Analysis of miRNAs with SNPs in Their Seed Region

The clustering information of human miRNAs was obtained from Wang et al. [14] to investigate the clustering patterns of the miRNAs with SNPs in their seed region. Specifically, clustering of miRNA genomic locations is determined if two neighboring miRNAs are located within 10 kb and are in the same strand. Based on this criterion, among all the 1881 precursor sequences and 2813 mature miRNAs annotated in humans, 352 miRNA genes including 634 mature miRNAs were grouped into 99 distinct clusters [14]. The significance of the difference in percentage of clustered miRNAs between miRNAs with or without SNPs in seed regions was calculated using the chi-square test. Additionally, a two-tailed Student’s t-test was used for comparisons between the percent overlap of the clustered and nonclustered miRNAs with SNPs in seed regions. value < 0.05 was considered statistically significantly.

2.4. Function and Pathway Analysis

The function annotation and pathway enrichment of the clustered and nonclustered miRNAs with SNPs in seed regions were performed using miRNA Enrichment Analysis and Annotation (miEAA) tool database (http://www.ccb.uni-saarland.de/mieaa_tool/) [17]. miEAA is a web-based system, which offers miRNA set enrichment analysis similar to gene set enrichment analysis (GSEA).

3. Results and Discussion

3.1. miRNAs with SNPs in Their Seed Region Are Significantly Enriched in Clusters

In total, we identified 1879 SNPs in 1226 (43.6%) human miRNA seed regions after mapping genetic variation onto human miRNA seed regions based on the genomic coordinates of SNPs in dbSNP human Build 147 and miRNAs in miRBase release 21 (Table S1). We found that most of the SNPs (1833 SNPs, 97.5%, Table S1) in miRNA seed regions were rare variants (defined as SNPs with minor allele frequency (MAF) < 5%). Recently, Torruella-Loran et al. [18] studied miRNA genetic variation in human populations and found that the seed regions tend to be depleted of high-frequency variants, which is consistent with our finding. We also found that there are 1587, 749, 340, 102, 31, and 4 miRNAs, which carry zero, one, two, three, four, and five SNPs, respectively, in their seed region (Figure S1, Table S1). This indicates that miRNA seed regions might be not so tolerant of genetic variants since most miRNAs have few or rare SNPs in their seed regions.

Further, we investigated the clustering patterns of the miRNAs with SNPs in their seed regions. An interesting observation is that the miRNAs with SNPs in their seed region are significantly enriched in clusters (Figure 1, Table S2). For the 1226 human miRNAs with SNPs in their seed region, 314 (25.6%) of them are located in miRNA clusters, whereas among the 1587 human miRNAs without SNPs in their seed region, only 320 (20.2%) of them are located in miRNA clusters (, χ2 test) (Table S2). miRNAs from the same cluster have the tendency to regulate the same sets of target genes and cooperatively repress expression levels of such genes [14]. Therefore, clustered miRNAs tend to be more tolerant of genetic variations in the seed region due to functional complementarity than are nonclustered miRNAs.

Figure 1: Clustering patterns of the miRNAs with SNPs in their seed region. Percentage of clustered miRNAs in miRNAs with or without variations in their seed region. miRNAs with SNPs in their seed region were significantly enriched in clusters. The y-axis is the percentage of miRNAs that are located in clusters in that group. .
3.2. SNPs in Clustered miRNA Seed Regions Have a Lower Functional Effect

We speculated that the functional effect of genetic variation in clustered miRNA seed regions may be relatively low compared with that in nonclustered miRNA seed regions. The functional effect or “cost” of SNPs in miRNA seed regions involves the loss of regulatory control over previously targeted mRNAs and/or the acquisition of novel regulatory control over previously untargeted mRNAs [16].

In order to test this hypothesis, we calculated the percent overlap (cosine similarity) of predicted targets for the reference and derived allele of each SNP by TargetRank. A lower percent overlap indicates large overall differences between the targets of reference and derived allele. The average percent overlap between reference target sets and derived target sets for all the 1879 SNPs is only 15.8% (Table S1). By considering the importance of miRNAs in regulating gene expression, this result suggests that as few as one nucleotide substitution within the seed region of miRNAs will cause a significant functional effect.

Then, we compared the percent overlap of SNPs in clustered miRNA seed regions with that of SNPs in nonclustered miRNA seed regions. We found that the percent overlap of SNPs in clustered miRNA seed regions is much higher than that of SNPs in nonclustered miRNA seed regions (Student’s t-test; ) (Figure 2). Therefore, it indicates that SNPs in clustered miRNA seed regions will produce a lower functional effect, which may be a result of functional complementarity and relaxed selection or adaptive evolution.

Figure 2: Difference in the functional effect of SNPs in clustered and nonclustered miRNA seed regions. .
3.3. Clustered miRNAs with SNPs in Seed Regions Are Involved in More Pathways

Additionally, we applied a gene set enrichment analysis to understand on which pathways the clustered and nonclustered miRNAs with SNPs in seed regions participate using miEAA [17]. Table 1 gives most enriched pathway information with value < 0.01. Clustered miRNAs with SNPs in seed regions are involved in more pathways, such as Alzheimer disease amyloid secretase pathway, and TCR and IL signaling pathways, whereas nonclustered miRNAs with SNPs in seed regions are only enriched in TCR and Id signaling pathways (). From the point of evolution, a gene can allow more variations and obtain new function during evolution through gene duplication [19]. We guess that compared with nonclustered miRNAs, clustered miRNAs with SNPs in seed regions can take part in more intricate and complex regulating networks by obtaining SNPs in seed regions and functional complementarity of the members in the same cluster.

Table 1: Pathway enrichment analysis of clustered or nonclustered miRNAs with SNPs in seed regions by using miEAA.

4. Conclusions

In this study, we took advantage of microRNA gene location and genetic variability obtained from miRBase 21 and dbSNP database to systematically identify all substitutions located in human miRNA seed regions and explore the clustering pattern and functional effect of these SNPs. In total, we have identified 1879 SNPs in 1226 human miRNA seed regions. We found that miRNAs with SNPs in their seed region are significantly enriched in clusters. We also found that functional cost of genetic variations in clustered miRNA seed regions was relatively low due to functional complementarity compared with that of variations in nonclustered miRNA seed regions. Additionally, we found that clustered miRNAs with SNPs in seed regions are involved in more pathways. Taken together, our results broaden current knowledge on the genetic variation in human miRNA seed regions and demonstrate that SNPs in clustered miRNA seed regions can take part in more intricate and complex networks with lower functional cost by functional complementarity.

Disclosure

The funders had no further role in the study design; in the collection, analysis, and interpretation of data; in the writing of the report; or in the decision to submit the paper for publication.

Conflicts of Interest

All authors have reported no financial interests or potential conflicts of interest regarding the publication of this article.

Acknowledgments

Funding for this study was provided by the National Science Foundation of China (81601175 and 31601032), the Natural Science Foundation of Guangdong Province (2016A030310401), the Guangdong Medical Research Foundation (A2016539), and the Scientific Research Staring Foundation of Southern Medical University (PY2015N002 and PY2015N006).

Supplementary Materials

Supplementary 1. Figure S1: SNPs in human miRNA seed regions. Number of miRNAs that carry zero, one, two, three, four, and five SNPs, respectively, in their seed region.

Supplementary 2. Table S1: genome-wide identification of SNPs in human miRNA seed regions.

Supplementary 3. Table S2: chi-square test of clustering patterns of miRNAs with SNPs in their seed region.

References

  1. J. Krol, I. Loedige, and W. Filipowicz, “The widespread regulation of microRNA biogenesis, function and decay,” Nature Reviews Genetics, vol. 11, no. 9, pp. 597–610, 2010. View at Publisher · View at Google Scholar · View at Scopus
  2. H. Siomi and M. C. Siomi, “Posttranscriptional regulation of microRNA biogenesis in animals,” Molecular Cell, vol. 38, no. 3, pp. 323–332, 2010. View at Publisher · View at Google Scholar · View at Scopus
  3. W. Filipowicz, S. N. Bhattacharyya, and N. Sonenberg, “Mechanisms of post-transcriptional regulation by microRNAs: are the answers in sight?” Nature Reviews Genetics, vol. 9, no. 2, pp. 102–114, 2008. View at Publisher · View at Google Scholar · View at Scopus
  4. B. P. Lewis, C. B. Burge, and D. P. Bartel, “Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets,” Cell, vol. 120, no. 1, pp. 15–20, 2005. View at Publisher · View at Google Scholar · View at Scopus
  5. D. P. Bartel, “MicroRNAs: genomics, biogenesis, mechanism, and function,” Cell, vol. 116, no. 2, pp. 281–297, 2004. View at Publisher · View at Google Scholar · View at Scopus
  6. A. Bhattacharya, J. D. Ziebarth, and Y. Cui, “PolymiRTS database 3.0: linking polymorphisms in microRNAs and their target sites with human diseases and biological pathways,” Nucleic Acids Research, vol. 42, Database issue, pp. D86–D91, 2013. View at Publisher · View at Google Scholar · View at Scopus
  7. J. Gong, C. Liu, W. Liu et al., “An update of miRNASNP database for better SNP selection by GWAS data, miRNA expression and online tools,” Database, vol. 2015, article bav029, 2015. View at Publisher · View at Google Scholar · View at Scopus
  8. M. Zorc, J. Obsteter, P. Dovc, and T. Kunej, “Genetic variability of microRNA genes in 15 animal species,” Journal of Genomics, vol. 3, pp. 51–56, 2015. View at Publisher · View at Google Scholar
  9. D. Bhartiya, S. V. Laddha, A. Mukhopadhyay, and V. Scaria, “miRvar: a comprehensive database for genomic variations in microRNAs,” Human Mutation, vol. 32, no. 6, pp. E2226–E2245, 2011. View at Publisher · View at Google Scholar · View at Scopus
  10. C. Borel and S. E. Antonarakis, “Functional genetic variation of human miRNAs and phenotypic consequences,” Mammalian Genome, vol. 19, no. 7-8, pp. 503–509, 2008. View at Publisher · View at Google Scholar · View at Scopus
  11. B. W. Iliff, S. A. Riazuddin, and J. D. Gottsch, “A single-base substitution in the seed region of miR-184 causes EDICT syndrome,” Investigative Ophthalmology & Visual Science, vol. 53, no. 1, pp. 348–353, 2012. View at Publisher · View at Google Scholar · View at Scopus
  12. Á. Mencía, S. Modamio-Høybjør, N. Redshaw et al., “Mutations in the seed region of human miR-96 are responsible for nonsyndromic progressive hearing loss,” Nature Genetics, vol. 41, no. 5, pp. 609–613, 2009. View at Publisher · View at Google Scholar · View at Scopus
  13. A. Bhattacharya and Y. Cui, “SomamiR 2.0: a database of cancer somatic mutations altering microRNA-ceRNA interactions,” Nucleic Acids Research, vol. 44, no. D1, pp. D1005–D1010, 2016. View at Publisher · View at Google Scholar · View at Scopus
  14. Y. Wang, J. Luo, H. Zhang, and J. Lu, “microRNAs in the same clusters evolve to coordinately regulate functionally related genes,” Molecular Biology and Evolution, vol. 33, no. 9, pp. 2232–2247, 2016. View at Publisher · View at Google Scholar · View at Scopus
  15. C. B. Nielsen, N. Shomron, R. Sandberg, E. Hornstein, J. Kitzman, and C. B. Burge, “Determinants of targeting by endogenous and exogenous microRNAs and siRNAs,” RNA, vol. 13, no. 11, pp. 1894–1910, 2007. View at Publisher · View at Google Scholar · View at Scopus
  16. C. G. Hill, N. Jabbari, L. V. Matyunina, and J. F. McDonald, “Functional and evolutionary significance of human MicroRNA seed region mutations,” PLoS One, vol. 9, no. 12, article e115241, 2014. View at Publisher · View at Google Scholar · View at Scopus
  17. C. Backes, Q. T. Khaleeq, E. Meese, and A. Keller, “miEAA: microRNA enrichment analysis and annotation,” Nucleic Acids Research, vol. 44, no. W1, pp. W110–W116, 2016. View at Publisher · View at Google Scholar · View at Scopus
  18. I. Torruella-Loran, H. Laayouni, B. Dobon et al., “MicroRNA genetic variation: from population analysis to functional implications of three allele variants associated with cancer,” Human Mutation, vol. 37, no. 10, pp. 1060–1073, 2016. View at Publisher · View at Google Scholar · View at Scopus
  19. C. Pegueroles, S. Laurie, and M. Mar Albà, “Accelerated evolution after gene duplication: a time-dependent process affecting just one copy,” Molecular Biology and Evolution, vol. 30, no. 8, pp. 1830–1842, 2013. View at Publisher · View at Google Scholar · View at Scopus