miRNA Target Binding Sites Research Articles

Dissecting the Pathogenesis of Diabetic Retinopathy Based on the Biological ceRNA Network and Genome Variation Disturbance.

Background Diabetic retinopathy (DR) is the most important manifestation of diabetic microangiopathy. It is essential to explore the gene regulatory relationship and genomic variation disturbance of biological networks in DR progression. Methods In this study, we constructed a comprehensive lncRNA-mRNA ceRNA network of DR procession (CLMN) and explored its topological characteristics. Results Modular and functional analysis indicated that the organization of CLMN performed fundamental and specific functions in diabetes and DR pathology. The differential expression of hub ceRNA nodes and positive correlation reveals the highly connected ceRNA regulation and important roles in the regulating of DR pathology. A large proportion of SNPs in the TFBS, DHS, and enhancer regions of lncRNAs will affect lncRNA transcription and further cause expression variation. Some SNPs were found to disrupt the lncRNA functional elements such as miRNA target binding sites. These results indicate the complex nature of genotypic effects in the disturbing of CLMN and further contribute to gene expression variation and different disease phenotypes. Conclusion The identification of individual genomic variations and analysis of biological network disturbance by these genomic variations will help provide more personalized treatment plans and promote the development of precision medicine for DR.

MiRNA target-binding sites as regulators of genes involved in the lipid metabolism might explain the hypercholesterolaemia in FH patients

Background and Aims: Familial hypercholesterolemia (FH) is due to mutations in LDLR, APOB and PCSK9, however, about 50% of clinical FH patients do not have an identifiable genetic cause. MicroRNAs (small non-coding RNAs) are negative regulators of gene expression and creation of new biding sites can be the cause of hypercholesterolaemia in mutation-negative families. The present work aims to analyse miRNAs targets as regulators of genes involved in lipid metabolism in FH patients.

MiRNASNP-v3: a comprehensive database for SNPs and disease-related variations in miRNAs and miRNA targets.

MicroRNAs (miRNAs) related single-nucleotide variations (SNVs), including single-nucleotide polymorphisms (SNPs) and disease-related variations (DRVs) in miRNAs and miRNA-target binding sites, can affect miRNA functions and/or biogenesis, thus to impact on phenotypes. miRNASNP is a widely used database for miRNA-related SNPs and their effects. Here, we updated it to miRNASNP-v3 (http://bioinfo.life.hust.edu.cn/miRNASNP/) with tremendous number of SNVs and new features, especially the DRVs data. We analyzed the effects of 7 161 741 SNPs and 505 417 DRVs on 1897 pre-miRNAs (2630 mature miRNAs) and 3′UTRs of 18 152 genes. miRNASNP-v3 provides a one-stop resource for miRNA-related SNVs research with the following functions: (i) explore associations between miRNA-related SNPs/DRVs and diseases; (ii) browse the effects of SNPs/DRVs on miRNA-target binding; (iii) functional enrichment analysis of miRNA target gain/loss caused by SNPs/DRVs; (iv) investigate correlations between drug sensitivity and miRNA expression; (v) inquire expression profiles of miRNAs and their targets in cancers; (vi) browse the effects of SNPs/DRVs on pre-miRNA secondary structure changes; and (vii) predict the effects of user-defined variations on miRNA-target binding or pre-miRNA secondary structure. miRNASNP-v3 is a valuable and long-term supported resource in functional variation screening and miRNA function studies.

Evaluation of polymorphisms in microRNA-binding sites and pancreatic cancer risk in Chinese population.

As promising biomarkers and therapy targets, microRNAs (miRNAs) are involved in various physiological and tumorigenic processes. Genetic variants in miRNA‐binding sites can lead to dysfunction of miRNAs and contribute to disease. However, systematic investigation of the miRNA‐related single nucleotide polymorphisms (SNPs) for pancreatic cancer (PC) risk remains elusive. We performed integrative bioinformatics analyses to select 31 SNPs located in miRNA‐target binding sites using the miRNASNP v2.0, a solid database providing miRNA‐related SNPs for genetic research, and investigated their associations with risk of PC in two large case‐control studies totally including 1847 cases and 5713 controls. We observed that the SNP rs3802266 is significantly associated with increased risk of PC (odds ratio (OR) = 1.21, 95% confidence intervals (CI) = 1.11‐1.31, P = 1.29E‐05). Following luciferase reporter gene assays show that rs3802266‐G creates a stronger binding site for miR‐181a‐2‐3p in 3′ untranslated region (3′UTR) of the gene ZHX2. Expression quantitative trait loci (eQTL) analysis suggests that ZHX2 expression is lower in individuals carrying rs3802266‐G with increased PC risk. In conclusion, our findings highlight the involvement of miRNA‐binding SNPs in PC susceptibility and provide new clues for PC carcinogenesis.

Functional Polymorphism Located in the microRNA Binding Site of the Insulin Receptor (INSR) Gene Confers Risk for Type 2 Diabetes Mellitus in the Bangladeshi Population.

Bangladesh has the second largest number of adults with diabetes in South Asia. Compelling evidence suggest that miRNAs contribute to the etiology of Type 2 diabetes mellitus (T2DM) by regulating many aspects of glucose homeostasis. Hence, we hypothesized that genetic polymorphisms in the diabetes-related miRNA target-binding sites could be associated with the risk of T2DM in Bangladesh. The reference Single nucleotide polymorphism (SNP) data from the Insulin Receptor (INSR) gene were downloaded from the ENSEMBL genome browser release 88 and further analyzed in silico for identifying SNPs with deleterious effect and clinical relationships. Further, case-control study using the microRNA-binding site polymorphism rs1366600 (T > C) located at the 3' UTR of the INSR gene was carried out in 217 T2DM patients and 237 healthy controls from Bangladesh. Genotyping was performed using the real time PCR based allele discrimination method. The results showed that the minor allele 'C' is associated with increased risk of T2DM [Odds ratio (OR) 1.87; 95% confidence intervals (CI) 1.28-2.74; P = 0.0010]. When we dissected our analysis to include the dominant model (CC + TC genotype against the TT genotype), we found that the CC and TC genotypes were associated with increased risk of T2DM in Bangladeshi population (OR 2.01; 95% CI 1.31-3.07; P = 0.0012). However, in recessive model (CC vs TT + TC); the effect was not statistically significant (OR 2.23; 95% CI 0.66-7.51; P = 0.1848). Stratification of our data based on the gender of the cases and controls showed similar degree of risk association with respect to different genotypes and alleles. Our study showed that the miRNA binding site polymorphism rs1366600 located at the 3'-UTR region of the INSR gene is associated with increased risk of T2DM in Bangladeshi individuals.

Genome-wide analysis of microRNA targeting impacted by SNPs in cucumber genome

BackgroundmicroRNAs (miRNAs) are endogenous small RNAs that play important regulatory functions in plant development. Genetic variations in miRNAs sequences or their target-binding sites (microRNA-target interaction sites) can alter miRNA targets in animal and human. Whether these single nucleotide polymorphisms (SNPs) in plant are functional have not yet been determined.ResultsIn this study, we constructed leaf, root, and stem-derived small libraries of cucumber (Cucumis sativus) line 9930 (cultivated China-group cucumber) and C. sativus var. hardwickii (wild India group cucumber). A total of 22 conserved miRNA families, nine less-conserved miRNA families, and 49 cucumber-specific miRNAs were identified in both line 9930 and hardwickii. We employed cucumber resequencing data to perform a genome-wide scan for SNPs in cucumber miRNA-target interaction sites, including miRNA mature sequences and miRNA-target binding sites. As a result, we identified a total of 19 SNPs in mature miRNA sequences and 113 SNPs in miRNA-target binding sites with the potential to affect miRNA-target interactions. Furthermore, we experimentally confirmed that these SNPs produced 14 9930-unique targets mRNAs and 15 hardwickii-unique targets mRNA for cucumber miRNAs. This is the first experimental validation of SNPs in miRNA-target interaction sites affecting miRNA-target binding in plants.ConclusionsOur results indicate that SNPs can alter miRNA function and produce unique miRNA targets in cultivated and wild cucumbers. Therefore, miRNA-related SNPs may have played important in events that led to the agronomic differences between domestic and wild cucumber.

Tumor-promoting function of single nucleotide polymorphism rs1836724 (C3388T) alters multiple potential legitimate microRNA binding sites at the 3'-untranslated region of ErbB4 in breast cancer.

ErbB4 can act as either a tumor-suppressor gene or an oncogene in breast cancer. Multiple genetic factors including single nucleotide polymorphisms (SNPs) affect gene expression patterns. Multiple 3'-untranslated region (3'-UTR) SNPs reside within the target binding site of microRNAs, which can strengthen or weaken binding to target genes. The present study aimed to predict potential 3'‑UTR variants of ErbB4 that alter the target binding site of microRNAs (miRNAs) and to clarify the association of the potential variant with the risk of developing breast cancer. Insilico prediction was performed to identify potential functional SNPs within miRNA target binding sites in the 3'‑UTR of ErbB4. Thus, 146patients and controls were genotyped using restriction fragment length polymorphism-polymerase chain reaction. In addition to the Cochran-Armitage test for trend, allele and genotype frequency differences were determined to investigate the association between rs1836724 and the susceptibility to breast cancer. Bioinformatics analysis identified rs1836724 to be a polymorphism in the seed region of four miRNA binding sites (hsa-miR335-5p, hsa-miR-28-5p, has‑miR‑708‑5p and has‑miR‑665), which may participate in the development of breast cancer. Logistic regression data indicated that the Tallele of the polymorphism [OR (95% CI)=1.72 (1.056‑2.808), P=0.029] is associated with the risk of breast cancer. Using bioinformatics tools, a correlation was indicated between the presence of the T allele and a reduction in ErbB4 RNA silencing based on miRNA interaction. Furthermore, case subgroup data analysis revealed an association between the C/T genotype and an ER positive phenotype [OR (95% CI)=6.00 (1.082‑33.274), P=0.028] compared with the T/T genotype. ErbB4 and estrogen receptor 1 (ESR1) are regulated by identical miRNAs thus there may be a competition for binding sites. Due to this pattern, if the interaction between miRNAs with one gene is reduced, it may be consistent with the increase in interaction with another one. Therefore, more interaction with rs1836724 Cvariant within ErbB4 may be associated with higher expression of ESR1 (ER‑positive phenotype). miRNAs interact with ErbB4 mRNA more frequently when it carries C allele at the rs1836724 position compared with the T carriers. Therefore, the identical miRNA interacts with ESR1 less frequently when ErbB4 mRNA has a C allele. Therefore, ESR1 expression may be higher when ErbB4 mRNA has a C allele.

Detection of canonical A-to-G editing events at 3' UTRs and microRNA target sites in human lungs using next-generation sequencing.

RNA editing is a post-transcriptional modification of RNA. The majority of these changes result from adenosine deaminase acting on RNA (ADARs) catalyzing the conversion of adenosine residues to inosine in double-stranded RNAs (dsRNAs). Massively parallel sequencing has enabled the identification of RNA editing sites in human transcriptomes. In this study, we sequenced DNA and RNA from human lungs and identified RNA editing sites with high confidence via a computational pipeline utilizing stringent analysis thresholds. We identified a total of 3,447 editing sites that overlapped in three human lung samples, and with 50% of these sites having canonical A-to-G base changes. Approximately 27% of the edited sites overlapped with Alu repeats, and showed A-to-G clustering (>3 clusters in 100 bp). The majority of edited sites mapped to either 3' untranslated regions (UTRs) or introns close to splice sites; whereas, only few sites were in exons resulting in non-synonymous amino acid changes. Interestingly, we identified 652 A-to-G editing events in the 3' UTR of 205 target genes that mapped to 932 potential miRNA target binding sites. Several of these miRNA edited sites were validated in silico. Additionally, we validated several A-to-G edited sites by Sanger sequencing. Altogether, our study suggests a role for RNA editing in miRNA-mediated gene regulation and splicing in human lungs. In this study, we have generated a RNA editome of human lung tissue that can be compared with other RNA editomes across different lung tissues to delineate a role for RNA editing in normal and diseased states.

Variations in mismatch repair genes and colorectal cancer risk and clinical outcome.

DNA mismatch repair (MMR) deficiency is one of the best understood forms of genetic instability in colorectal cancer (CRC). CRC is routinely cured by 5-fluorouracil (5-FU)-based chemotherapy, with a prognostic effect and resistance to such therapy conferred by MMR status. In this study, we aimed to analyse the effect of genetic variants in classical coding regions or in less-explored predicted microRNA (miRNA)-binding sites in the 3' untranslated region (3'UTR) of MMR genes on the risk of CRC, prognosis and the efficacy of 5-FU therapy. Four single nucleotide polymorphisms (SNPs) in MMR genes were initially tested for susceptibility to CRC in a case-control study (1095 cases and 1469 healthy controls). Subsequently, the same SNPs were analysed for their role in survival on a subset of patients with complete follow-up. Two SNPs in MLH3 and MSH6 were associated with clinical outcome. Among cases with colon and sigmoideum cancer, carriers of the CC genotype of rs108621 in the 3'UTR of MLH3 showed a significantly increased survival compared to those with the CT + TT genotype (log-rank test, P = 0.05). Moreover, this polymorphism was also associated with an increased risk of relapse or metastasis in patients with heterozygous genotype (log-rank test, P = 0.03). Patients carrying the CC genotype for MSH6 rs1800935 (D180D) and not undergoing 5-FU-based chemotherapy showed a decreased number of recurrences (log-rank test, P = 0.03). No association with CRC risk was observed. We provide the first evidence that variations in potential miRNA target-binding sites in the 3'UTR of MMR genes may contribute to modulate CRC prognosis and predictivity of therapy.

A miR-151 binding site polymorphism in the 3′-untranslated region of the cyclin E1 gene associated with nasopharyngeal carcinoma

Genetic alterations in nasopharyngeal carcinoma (NPC) have been reported in previous works. However, it remains unclear whether polymorphisms within the miRNA-target binding sites are associated with individual NPC risk. In this study, new experimental and computational approaches were developed to assess the polymorphism frequency distribution within the miRNA sites in NPC patients, and to explore its association with NPC risk. We focused on 220 single-nucleotide polymorphisms (SNPs) in the 3′-untranslated regions (3′UTRs) of 32 genes carrying putative miRNA-binding sites by specialized algorithms. A total of 9 candidate genes were selected for further investigation, which were reportedly overexpressed in NPC, including EGFR, COX2, CCNE1, hTERT, MMP2, MMP9, NF-κB VEGF, and WNT3. SNPs in 3′UTRs were genotyped by direct polymerase chain reaction sequencing of the genomic DNA of 24 cases and 24 controls. Then, EGFR rs884225, CCNE1 rs3218073, and MMP2 rs7201 were screened with large samples. Based on the analysis of a series of 167 NPC cases and 171 controls from Guangdong Province, statistically significant associations were found between NPC risk and variant genotypes of CCNE1 rs3218073 for TC+TT (OR=1.585; 95% CI=1.023–2.458; P=0.046) and for T-allele (OR=1.464; 95% CI=1.012–2.118; P=0.042). In addition, a significant association among rs3218073 genotype TC (OR=1.959, P=0.043), T-allele (OR=2.123, P=0.006), and primary tumor (T3–T4) was retrieved. Genotype TC (OR=1.959, P=0.043) and T-allele (OR=2.123, P=0.006) of rs3218073 were correlated with increased risk of higher NPC stage (III to IV). In support of the postulation that the 3′UTR SNP directly affected miRNA-binding site, luciferase reporter assay indicated that CCNE1 was a direct target of miR-151, and the rs3218073 T>C change resulted in altered regulation of CCNE1 expression. By contrast, no statistically significant association with NPC risk was found for MMP2 rs7201 and EGFR rs884225 polymorphisms (P>0.05). In conclusion, our data demonstrate that CCNE1 rs3218073 polymorphism located at miRNA-151 binding site is associated with NPC susceptibility and is correlated with NPC stage. These results suggest that CCNE1 rs3218073 polymorphism can be exploited as a novel biomarker for future NPC diagnosis and prognosis.

Single-nucleotide polymorphisms inside microRNA target sites influence the susceptibility to type 2 diabetes

MicroRNAs (miRNAs) negatively regulate target gene expression by binding to 3'-untranslated region of target mRNA. Single-nucleotide polymorphisms (SNPs) that reside in the miRNA target sites can affect the bindings of miRNAs to mRNAs. Compelling evidence has shown that miRNAs contribute to the etiology of type 2 diabetes (T2DM). We hypothesized that SNPs in diabetes-related miRNA target-binding sites could be associated with the risk of T2DM. We selected 10 SNPs on miRNA-binding sites by using bioinformatics software. Genotypes of T2DM patients (n=1017) and normal controls (n=1059) were analyzed by TaqMan assay. The variant genotypes rs1366600CC and TC/CC in the insulin receptor (INSR) gene, rs2292899GA in the acyl-CoA synthetase 1 (ACSL1) gene and rs11724758AA in the fatty-acid-binding protein 2 (FABP2) gene were associated with T2DM (adjusted odds ratios (ORs) (95% confidence intervals)=2.03 (1.02-4.01), 1.28 (1.04-1.57), 1.22 (1.004-1.49) and 0.76 (0.58-0.997), respectively). The analysis stratified by age, gender, waist circumstance and living habits also revealed these genotypes' effect. Furthermore, crossover analysis indicated those who exposed to both environmental factor and putative risk genotypes did have the highest risk of T2DM. A cumulative effect of SNPs rs1366600, rs2292899 and rs11724758 was observed with individuals carrying 2, 3 and 4-6 risk alleles having a gradually increased risk of T2DM (OR=1.52, 1.81 and 2.28, P for trend <0.001). Our result suggested that INSR rs1366600, ACSL1 rs2292899 and FABP2 rs11724758 could influence the susceptibility to T2DM in Chinese Han population, most likely through their effects on the specific miRNA-binding sites and functional characterizations of three genes are needed.

Polymorphisms in miRNA-binding sites of nucleotide excision repair genes and colorectal cancer risk

Reduced DNA repair capacity and DNA damage accumulation may lead to cancer development. Regulation of and coordination between genes involved in DNA repair pathways is fundamental for maintaining genome stability, and post-transcriptional gene regulation by microRNAs (miRNAs) may therefore be of particular relevance. In this context, the presence of single nucleotide polymorphisms (SNPs) within the 3'untranslated regions of target DNA repair genes could alter the binding with specific miRNAs, modulating gene expression and ultimately affecting cancer susceptibility. In this study, we investigated the role of genetic variations in miRNA-binding sites of nucleotide excision repair (NER) genes in association with colorectal cancer (CRC) risk. From 28 NER genes, we screened among SNPs residing in their 3'untranslated regions and simultaneously located in miRNA-binding sites, with an in silico approach. Through the calculation of different binding free energy according to both alleles of identified SNPs, and with global binding free energies median providing a threshold, we selected nine NER gene variants. We tested those SNPs in 1098 colorectal cancer cases and 1469 healthy controls from the Czech Republic. Rs7356 in RPA2 and rs4596 in GTF2H1 were associated with colorectal cancer risk. After stratification for tumor location, the association of both SNPs was significant only for rectal cancer (rs7356: OR 1.52, 95% CI 1.02-2.26, P = 0.04 and rs4596: OR 0.69, 95% CI 0.50-0.94, P = 0.02; results not adjusted for multiple testing). Variation in miRNA target binding sites in the 3'untranslated region of NER genes may be important for modulating colorectal cancer risk, with a different relevance according to tumor location.

Implication of the miR-184 and miR-204 Competitive RNA Network in Control of Mouse Secondary Cataract

The high recurrence rate of secondary cataract (SC) is caused by the intrinsic differentiation activity of residual lens epithelial cells after extra-capsular lens removal. The objective of this study was to identify changes in the microRNA (miRNA) expression profile during mouse SC formation and to selectively manipulate miRNA expression for potential therapeutic intervention. To model SC, mouse cataract surgery was performed and temporal changes in the miRNA expression pattern were determined by microarray analysis. To study the potential SC counterregulative effect of miRNAs, a lens capsular bag in vitro model was used. Within the first 3 wks after cataract surgery, microarray analysis demonstrated SC-associated expression pattern changes of 55 miRNAs. Of the identified miRNAs, miR-184 and miR-204 were chosen for further investigations. Manipulation of miRNA expression by the miR-184 inhibitor (anti-miR-184) and the precursor miRNA for miR-204 (pre-miR-204) attenuated SC-associated expansion and migration of lens epithelial cells and signs of epithelial to mesenchymal transition such as α-smooth muscle actin expression. In addition, pre-miR-204 attenuated SC-associated expression of the transcription factor Meis homeobox 2 (MEIS2). Examination of miRNA target binding sites for miR-184 and miR-204 revealed an extensive range of predicted target mRNA sequences that were also a target to a complex network of other SC-associated miRNAs with possible opposing functions. The identification of the SC-specific miRNA expression pattern together with the observed in vitro attenuation of SC by anti-miR-184 and pre-miR-204 suggest that miR-184 and miR-204 play a significant role in the control of SC formation in mice that is most likely regulated by a complex competitive RNA network.

An integrated map of cattle candidate genes for mastitis: a step forward to new genetic markers

To facilitate the development of new genetic markers for mastitis resistance or susceptibility we used genome-wide comparative approach to review all known mastitis-associated loci. We assembled into a map 233 loci that were identified by six different study approaches (QTLs, association studies, expression experiments, AFLP-, miRNA- and epigenetic- studies). To integrate data from different sources and to identify overlapping regions we presented the results in the form of genetic map. The collected data represent genetic background for mastitis-related traits in cattle. Thirty most promising candidate genes (associated with mastitis in different study approaches or at least in two independent studies, or/and overlapping with QTL regions) were selected from database and in silico searched for genetic variability and putative miRNA target sites in 3’UTR. Thirty-one SNPs in the putative promoter/5' UTR (up to 2 kb upstream) of eight candidate genes were found. Bioinformatic analysis revealed that some promoter SNPs might potentially cause gain/loss of the putative transcription factors. Promoter SNPs were also present in CG dinucleotides and therefore possibly involved in gain/loss of CpG sites. In ten mastitis candidate genes we found 56 SNPs in exons of which 21 were non-synonymous substitutions. Additionally, 23 SNPs in intronic regions and 21 SNPs in 3’UTR were found. MiRNA target analysis revealed 89 putative target sites in 18 candidate genes; however, current SNPs were not identified in the miRNA target binding sites or miRNAs expressed in mammary gland. For SNPs with a putative regulatory role found in candidate genes, functional analyses and association studies are needed to facilitate identification of mastitis resistance or susceptibility alleles possibly involved in mastitis regulatory pathways.

Polymorphisms within micro-RNA-binding sites and risk of sporadic colorectal cancer

Recent evidence indicate that small non-coding RNA molecules, called micro-RNAs (miRNAs), can bind to the 3' untranslated regions (UTRs) of messenger RNAs and interfere with their translation, thereby regulating cell growth, differentiation, apoptosis and tumorigenesis. Genetic polymorphisms can reside on miRNA-binding sites. Thus, it is conceivable that the miRNA regulation may be affected by polymorphisms on the 3' UTRs. Since gene deregulation is one of the key mechanisms by which cells can progress to cancer, we hypothesize that common polymorphisms within miRNA-target binding sites could play a role in the individual risk of cancer. In the present study, we selected the 3' UTRs of 104 genes candidate for colorectal cancer (CRC) and we identified putative miRNA-binding sites by specialized algorithms (PicTar, DianaMicroT, miRBase, miRanda, TargetScan and microInspector). Fifty-seven single-nucleotide polymorphisms (SNPs) were identified in miRNA-binding sites. We evaluated the SNPs for their ability to affect the binding of the miRNA with its target, by assessing the variation of Gibbs free energy between the two alleles of each SNP. We found eight common polymorphisms that were further investigated by a case-control association studies. The study was carried out on a series of cases and controls from Czech Republic, a population with the highest worldwide incidence of CRC. We found statistically significant associations between risk of CRC and variant alleles of CD86 [odds ratio (OR) = 2.74; 95% confidence interval (CI) = 1.24-6.04, for the variant homozygotes] and INSR genes (OR = 1.94; 95% CI = 1.03-3.66, for the variant homozygotes). These results are the first reporting positive association between miRNA-binding SNPs sequences and cancer risk.