MicroRNA Genetic Variation: From Population Analysis to Functional Implications of Three Allele Variants Associated with Cancer

Signatures of genetic variation in human microRNAs point to processes of positive selection and population-specific disease risks

The occurrence of natural variation in human microRNAs has been the focus of numerous studies during the last 20 years. Most of them have been focused on the role of specific mutations in disease, while a minor proportion seek to analyse microRNA diversity in the genomes of human populations. We analyse the latest human microRNA annotations in the light of the most updated catalogue of genetic variation provided by the 1000 Genomes Project. By means of the in silico analysis of microRNA genetic variation we show that the level of evolutionary constraint of these sequences is governed by the interplay of different factors, like their evolutionary age or genomic location. The role of mutations in the shaping of microRNA-driven regulatory interactions is emphasized with the acknowledgement that, while the whole microRNA sequence is highly conserved, the seed region shows a pattern of higher genetic diversity that appears to be caused by the dramatic frequency shifts of a fraction of human microRNAs. We highlight the participation of these microRNAs in population-specific processes by identifying that not only the seed, but also the loop, are particularly differentiated regions among human populations. The quantitative computational comparison of signatures of population differentiation showed that candidate microRNAs with the largest differences are enriched in variants implicated in gene expression levels (eQTLs), selective sweeps and pathological processes. We explore the implication of these evolutionary-driven microRNAs and their SNPs in human diseases, such as different types of cancer, and discuss their role in population-specific disease risk.

MiRLog and dbmiR: prioritization and functional annotation tools to study human microRNA sequence variants

The recent identification of noncoding variants with pathogenic effects suggests that these variations could underlie a significant number of undiagnosed cases. Several computational methods have been developed to predict the functional impact of noncoding variants, but they exhibit only partial concordance and are not integrated with functional annotation resources, making the interpretation of these variants still challenging. MicroRNAs (miRNAs) are small noncoding RNA molecules that act as fine regulators of gene expression and play crucial functions in several biological processes, such as cell proliferation and differentiation. An increasing number of studies demonstrate a significant impact of miRNA single nucleotide variants (SNVs) both in Mendelian diseases and complex traits. To predict the functional effect of miRNA SNVs, we implemented a new meta‐predictor, MiRLog, and we integrated it into a comprehensive database, dbmiR, which includes a precompiled list of all possible miRNA allelic SNVs, providing their biological annotations at nucleotide and miRNA levels. MiRLog and dbmiR were used to explore the genetic variability of miRNAs in 15,708 human genomes included in the gnomAD project, finding several ultra‐rare SNVs with a potentially deleterious effect on miRNA biogenesis and function representing putative contributors to human phenotypes.

The Association between Five Genetic Variants in MicroRNAs (rs2910164, rs11614913, rs3746444, rs11134527, and rs531564) and Cervical Cancer Risk: A Meta-Analysis

The objective of this study was to conduct a meta-analysis to systematically summarize and investigate the association of miRNA-124 rs531564, miRNA-218 rs11134527, miRNA-146a rs2910164, miRNA-196a2 rs11614913, and miRNA-499 rs3746444 polymorphisms with cervical cancer. A systematic review was performed to identify relevant studies using Embase and PubMed databases. A chi-square-based Q-test combined with the inconsistency index (I 2) was used to check the heterogeneity between studies. A total of six case-control studies on rs2910164 and rs11614913, 4 studies on rs3746444 and rs11134527, and three studies on rs531564 were included. No evidence of association was found between miR-146a rs2910164, miR-196a2 rs11614913, miRNA-499 rs3746444, and miR-218 rs11134527 polymorphisms and cervical cancer risk in all the genetic models. The miR-124 rs531564 polymorphism was associated with a statistically increased risk of cervical cancer in a homozygote model (CC vs. GG: OR = 2.87, 95% CI: 1.40-5.91, P H = 0.887), dominant model (GC/CC vs. GG: OR = 1.38, 95% CI: 1.07-1.80, P H = 0.409), and recessive model (CC vs. GC/GG: OR = 2.26, 95% CI: 1.58-3.23, P H = 0.979). However, this finding should be interpreted with caution for limited samples and heterogeneity. Large-scale and well-designed studies are needed to validate our result.

Association between MIR499A rs3746444 polymorphism and breast cancer susceptibility: a meta-analysis

Numerous studies have investigated the association of MIR499A rs3746444 polymorphism with breast cancer susceptibility, but the results have been inconsistent. In this work, we performed a meta-analysis to obtain a more reliable estimate of the association between the polymorphism and susceptibility to breast cancer. A comprehensive literature search was conducted on PubMed, Scopus, Web of Science (WoS), China National Knowledge Infrastructure (CNKI), VIP and Wanfang databases up to January 2020. A total of 14 studies involving 6,797 cases and 8,534 controls were included for analysis under five genetic models: homozygous (GG vs. AA), heterozygous (AG vs. AA), dominant (AG + GG vs. AA), recessive (GG vs. AA + AG) and allele (G vs. A). A statistically significant association was observed between the polymorphism and an increased breast cancer susceptibility under all genetic models (homozygous, OR = 1.33, 95% CI = 1.03-1.71, P = 0.03; heterozygous, OR = 1.08, 95% CI = 1.00-1.16, P = 0.04; dominant, OR = 1.15, 95% CI = 1.02-1.30; P = 0.03; recessive, OR = 1.35, 95% CI = 1.06-1.72, P = 0.01; allele, OR = 1.12, 95% CI = 1.00-1.26, P = 0.04). Subgroup analysis based on ethnicity suggested that significant association was present only among Asians, but not Caucasians. In conclusion, MIR499A rs3746444 polymorphism was significantly associated with breast cancer susceptibility among Asians, suggesting its potential use as a genetic risk marker in this population.

Polymorphisms in miRNA genes play roles in the initiation and development of cervical cancer

MicroRNA deregulation is crucial for cancer development. Studies showed that polymorphisms in miRNA genes could affect miRNA expression, which might be associated with cancer development. In the current study, we investigated the association of seven single nucleotide polymorphisms (SNPs) in seven miRNA genes with the initiation and development of cervical cancer in a Chinese Han population. The SNPs of 358 cervical intraepithelial neoplasia (CIN) patients, 547 cervical cancer patients and 567 healthy individuals were genotyped using TaqMan assays. Moreover, we evaluated the association of the seven SNPs with the different stages of cervical cancer. Our results showed that rs4636297 in miR-126 was associated with susceptibility to CIN and cervical cancer (P=0.019 and 0.019, respectively) and that the T allele was associated with a higher risk of CIN (OR=1.334, 95% CI: 1.049-1.698) and cervical cancer (OR=1.296, 95% CI: 1.044-1.609). Similarly, rs11614913 in miR-125a was associated with CIN and cervical cancer (P=0.025 and 0.015, respectively), and the T allele might be the protective factor for CIN (OR=0.807, 95% CI: 0.669-0.974) and cervical cancer (OR=0.814, 95% CI: 0.689-0.961). Our results indicated that rs4636297 in miR-126 and rs11614913 in miR-196a2 play an important role only in the initiation of cervical cancer not in the development of CIN to cervical cancer.

rs12416605:C>T in MIR938 associates with gastric cancer through affecting the regulation of the CXCL12 chemokine gene

Background

MicroRNAs are small regulatory RNAs with important roles in carcinogenesis. Genetic variants in these regulatory molecules may contribute to disease. We aim to identify allelic variants in microRNAs as susceptibility factors to gastric cancer using association studies and functional approaches.

Methods

Twenty‐one single nucleotide variants potentially functional, because of their location in either the seed, mature or precursor region of 22 microRNAs, were selected for association studies. Genetic association with gastric cancer in 365 cases and 1,284 matched controls (European Prospective Investigation into Cancer and Nutrition Cohort) was analysed using logistic regression. MicroRNA overexpression, transcriptome analysis, and target gene validation experiments were performed for functional studies.

Results

rs3746444:T>C, in the seed of MIR499A and mature MIR499B, associated with the cardia adenocarcinoma location; rs12416605:C>T, in the seed of MIR938, associated with the diffuse subtype; and rs2114358:T>C, in the precursor MIR1206, associated with the noncardia phenotype. In all cases, the association was inverse, indicating a protective affect against gastric cancer of the three minor allelic variants. MIR499 rs3746444:T>C and MIR1206 rs2114358:T>C are reported to affect the expression of these miRNAs, but the effect of MIR938 rs12416605:C>T is unknown yet. Functional approaches showed that the expression of MIR938 is affected by rs12416605:C>T and revealed that MIR938 could regulate a subset of cancer‐related genes in an allele‐specific fashion. Furthermore, we demonstrated that CXCL12, a chemokine participating in gastric cancer metastasis, is specifically regulated by only one of the rs12416605:C>T alleles.

Conclusion

rs12416605 appears to be involved in gastric cancer by affecting the regulatory function of MIR938 on genes related to this cancer type, particularly on CXCL12 posttranscriptional regulation.

Mapping of microRNAs related to cervical cancer in Latin American human genomic variants

Background: MicroRNAs are related to human cancers, including cervical cancer (CC) caused by HPV. In 2018, approximately 56.075 cases and 28.252 deaths from this cancer were registered in Latin America and the Caribbean according to GLOBOCAN reports. The main molecular mechanism of HPV in CC is related to integration of viral DNA into the hosts' genome. However, the different variants in the human genome can result in different integration mechanisms, specifically involving microRNAs (miRNAs). Methods: The miRNAs associated with CC were obtained from literature, the miRNA sequences and four human genome variants (HGV) from Latin American populations were obtained from miRBase and 1000 Genomes Browser, respectively. HPV integration sites near cell cycle regulatory genes were identified. miRNAs were mapped on HGV. miRSNPs were identified in the miRNA sequences located at HPV integration sites on the Latin American HGV.  Results: Two hundred seventy-two miRNAs associated with CC were identified in 139 reports from different geographic locations. By mapping with Blast-Like Alignment Tool (BLAT), 2028 binding sites were identified from these miRNAs on the human genome (version GRCh38/hg38); 42 miRNAs were located on unique integration sites; and miR-5095, miR-548c-5p and miR-548d-5p were involved with multiple genes related to the cell cycle. Thirty-seven miRNAs were mapped on the Latin American HGV (PUR, MXL, CLM and PEL), but only miR-11-3p, miR-31-3p, miR-107, miR-133a-3p, miR-133a-5p, miR-133b, miR-215-5p, miR-491-3p, miR-548d-5p and miR-944 were conserved. Conclusions: Ten miRNAs were conserved in the four HGV. In the remaining 27 miRNAs, substitutions, deletions or insertions were observed. These variation patterns can imply differentiated mechanisms towards each genomic variant in human populations because of specific genomic patterns and geographic features. These findings may help in determining susceptibility for CC development. Further identification of cellular genes and signalling pathways involved in CC progression could lead new therapeutic strategies based on miRNAs.

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

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.

Genetic predisposition to lung cancer: comprehensive literature integration, meta-analysis, and multiple evidence assessment of candidate-gene association studies

More than 1000 candidate-gene association studies on genetic susceptibility to lung cancer have been published over the last two decades but with few consensuses for the likely culprits. We conducted a comprehensive review, meta-analysis and evidence strength evaluation of published candidate-gene association studies in lung cancer up to November 1, 2015. The epidemiological credibility of cumulative evidence was assessed using the Venice criteria. A total of 1018 publications with 2910 genetic variants in 754 different genes or chromosomal loci were eligible for inclusion. Main meta-analyses were performed on 246 variants in 138 different genes. Twenty-two variants from 21 genes (APEX1 rs1130409 and rs1760944, ATM rs664677, AXIN2 rs2240308, CHRNA3 rs6495309, CHRNA5 rs16969968, CLPTM1L rs402710, CXCR2 rs1126579, CYP1A1 rs4646903, CYP2E1 rs6413432, ERCC1 rs11615, ERCC2 rs13181, FGFR4 rs351855, HYKK rs931794, MIR146A rs2910164, MIR196A2 rs11614913, OGG1 rs1052133, PON1 rs662, REV3L rs462779, SOD2 rs4880, TERT rs2736098, and TP53 rs1042522) showed significant associations with lung cancer susceptibility with strong cumulative epidemiological evidence. No significant associations with lung cancer risk were found for other 150 variants in 98 genes; however, seven variants demonstrated strong cumulative evidence. Our findings provided the most updated summary of genetic risk effects on lung cancer and would help inform future research direction.