Lengthening of 3' Untranslated Regions of mRNAs by Alternative Polyadenylation Is Associated With Tumor Progression and Poor Prognosis of Clear Cell Renal Cell Carcinoma.

Alternative polyadenylation (APA) is emerging as a major posttranscriptional mechanism for gene regulation in cancer. A prevailing hypothesis is that shortening of the 3' untranslated region (3'UTR) increases oncoprotein expression because of the loss of miRNA-binding sites (MBSs). We showed that the longer 3'UTR is associated with a more advanced tumor stage in patients with clear cell renal cell carcinoma (ccRCC). More surprisingly, 3'UTR shortening is correlated with better overall survival in patients with ccRCC. Furthermore, we identified a mechanism by which longer transcripts lead to increased oncogenic protein and decreased tumor-suppressive protein expression compared to the shorter transcripts. In our model, shortening of 3'UTRs by APA may increase the mRNA stability of the majority of the potential tumor-suppressor genes due to the loss of MBSs and AU-rich elements (AREs). Unlike potential tumor-suppressor genes, the potential oncogenes display much lower MBS and ARE density and globally much higher m6A density in distal 3'UTRs. As a result, 3'UTRs shortening decreases the mRNA stability of potential oncogenes and enhances the mRNA stability of potential tumor-suppressor genes. Our findings highlight the cancer-specific pattern of APA regulation and extend our understanding of the mechanism of APA-mediated 3'UTR length changes in cancer biology.

Nuclear Pore Glycoprotein 62 Genetic Variant rs9523 is Associated with Clinical Outcomes of Epidermal Growth Factor Receptor Tyrosine Kinase Inhibitors in Lung Adenocarcinoma Patients.

INTRODUCTION: Epidermal growth factor receptor (EGFR) tyrosine kinase inhibitors (TKIs) have represented the prototype of targeted therapy in NSCLC. Patients with EGFR-mutant lung adenocarcinoma extract an extraordinary clinical benefit from EGFR-TKIs. However, the extent and duration of these responses are heterogeneous, suggesting the existence of genetic modifiers affecting an individual's response to TKIs. We investigated whether genetic variants in miRNA binding sites are associated with the clinical outcome of EGFR-TKIs in lung adenocarcinoma patients. METHODS: One hundred SNPs at miRNA binding sites in cancer-related genes were selected for the analysis using the crosslinking, ligation and sequencing of hybrids (CLASH) and CancerGenes database. qRT-PCR and luciferase assays were conducted to evaluate the functional relevance of the SNPs. RESULTS: NUP62 rs9523A>G were significantly associated with worse response to EGFR-TKIs, overall survival (OS), and progression-free survival (PFS). The other three SNPs (DVL2 rs2074216G>A, ARF1 rs11541557G>T, and UHRF1 rs2261988C>A) were significantly associated with worse OS and PFS. The rs9523A>G was significantly associated with decreased NUP62 expression in tumor tissues. In addition, a significantly decreased luciferase activity was noted in NUP62 rs9523 G allele compared to A allele. CONCLUSION: Genetic variants in miRNA binding sites, especially NUP62 rs9523A>G, may be useful in predicting the clinical outcomes of EGFR-mutant lung adenocarcinoma patients treated with EGFR-TKIs.

miRNA Combinatorics and its Role in Cell State Control-A Probabilistic Approach.

A hallmark of cancer evolution is that the tumor may change its cell identity and improve its survival and fitness. Drastic change in microRNA (miRNA) composition and quantities accompany such dynamic processes. Cancer samples are composed of cells' mixtures of varying stages of cancerous progress. Therefore, cell-specific molecular profiling represents cellular averaging. In this study, we consider the degree to which altering miRNAs composition shifts cell behavior. We used COMICS, an iterative framework that simulates the stochastic events of miRNA-mRNA pairing, using a probabilistic approach. COMICS simulates the likelihood that cells change their transcriptome following many iterations (100 k). Results of COMICS from the human cell line (HeLa) confirmed that most genes are resistant to miRNA regulation. However, COMICS results suggest that the composition of the abundant miRNAs dictates the nature of the cells (across three cell lines) regardless of its actual mRNA steady-state. In silico perturbations of cell lines (i.e., by overexpressing miRNAs) allowed to classify genes according to their sensitivity and resilience to any combination of miRNA perturbations. Our results expose an overlooked quantitative dimension for a set of genes and miRNA regulation in living cells. The immediate implication is that even relatively modest overexpression of specific miRNAs may shift cell identity and impact cancer evolution.

Sfold Tools for MicroRNA Target Prediction.

Computational prediction of miRNA binding sites on target mRNAs facilitates experimental investigation of miRNA functions. In this chapter, we describe STarMir and STarMirDB, two application modules of the Sfold RNA package. STarMir is a Web server for performing miRNA binding site predictions for mRNA and target sequences submitted by users. STarMirDB is a database of precomputed transcriptome-scale predictions. Both STarMir and STarMirDB provide comprehensive sequence, thermodynamic, and target structure features, a logistic probability as a measure of confidence for each predicted site, and a publication-quality diagram of the predicted miRNA-target hybrid. In addition, STarMir now offers a new quantitative score to address combined regulatory effects of multiple seed and seedless sites. This score provides a quantitative measure of the overall regulatory effects of both seed and seedless sites on the target. STarMir and STarMirDB are freely available to all through the Sfold Web application server at http://sfold.wadsworth.org .

Comparative genomics reveals origin of MIR159A-MIR159B paralogy, and complexities of PTGS interaction between miR159 and target GA-MYBs in Brassicaceae.

Whole-genome and segmental duplications coupled with sequence and functional diversification are responsible for gene family expansion, and morphological and adaptive diversity. Although broad contours of such processes are understood, detailed investigations on regulatory elements, such as miRNA-transcription factor modules, especially in non-model crop plants with complex genomes, are few. The present study was performed to understand evolutionary history of MIR159 family, and changes in the miRNA-binding site (MBS) of the targets MYB33, MYB65, and MYB101 that may affect post-transcriptional gene silencing. We established orthology and paralogy between members of MIR159 family by reconstructing the phylogeny based on 240 precursor sequences sampled across green plants. An unambiguous paralogous relationship between MIR159A and MIR159B was observed only in Brassicaceae which prompted us to analyze the origin of this paralogy. Comparative micro-synteny of ca. 100 kb genomic segments surrounding MIR159A, MIR159B, and MIR159C loci across 15 genomes of Brassicaceae revealed segmental duplication that occurred in the common ancestor of Brassicaceae to be responsible for origin of MIR159A-MIR159B paralogy; extensive gene loss and rearrangements were also encountered. The impact of polyploidy was revealed when the three sub-genomes-least fractionated (LF), moderately fractionated (MF1), and most fractionated (MF2) sub-genomes of Brassica and Camelina sativa-were analyzed. Extensive gene loss was observed among sub-genomes of Brassica, whereas those in Camelina were largely conserved. Analysis of the target MYBs revealed the complete loss of MYB33 homologs in a Brassica lineage-specific manner. Our findings suggest that mature miR159a/b /c are capable of targeting MYB65 across Brassicaceae, MYB33 in all species except Brassica, and MYB101 only in Arabidopsis thaliana. Comparative analysis of the mature miRNA sequence and the miRNA-binding site (MBS) in MYB33, MYB65, and MYB101 showed the complexity of regulatory network that is dependent on strict sequence complementarity potentially leading to regulatory diversity.

HDL-cholesterol concentration in pregnant Chinese Han women of late second trimester associated with genetic variants in CETP, ABCA1, APOC3, and GALNT2.

OBJECTIVE: To investigate whether HDL-C level in pregnant Chinese Han women of late second trimester correlated with loci in high-density lipoprotein-cholesterol (HDL-C)-related genes found in genome-wide association studies (GWAS). METHODS: Seven single-nucleotide polymorphisms (rs3764261 in CETP, rs1532085 in LIPC, rs7241918 in LIPG, rs1883025 in ABCA1, rs4225 in APOC3, rs1059611 in LPL, and rs16851339 in GALNT2) were genotyped using the Sequenom MassArray system for 1,884 pregnant women. RESULTS: The following polymorphisms were statistically associated with HDL-C level after adjusting for age, gestational week, pre-pregnancy BMI and state of GDM or HOMAIR: (i) rs3764261 (b = -0.055 mmol/L, 95% CI -0.101 to -0.008, p = 0.021), (ii) rs1883025 (b = -0.054 mmol/L, 95% CI -0.097 to -0.012, p = 0.013), (iii) rs4225 (b = -0.071 mmol/L, 95% CI -0.116 to -0.027, p = 1.79E-3) and (iv) rs16851339 (b = -0.064 mmol/L, 95% CI -0.120 to -0.008, p = 0.025). The more risk alleles the pregnant women have, the lower the plasma HDL-C levels of the subjects are. CONCLUSIONS: Several risk alleles found to be related to HDL-C in GWAS are also associated with HDL-C levels in pregnant Chinese Han women and these risk loci contribute additively to low HDL-C levels.

Dietary factors and microRNA-binding site polymorphisms in the IL13 gene: risk and prognosis analysis of colorectal cancer.

Long-term dietary intake influences the structure and activity of microorganisms residing in the human gut. The immune response and gut microbiota have a mutual influence on the risk of colorectal cancer (CRC). This study examines the association of gut microbiota-related dietary factors and polymorphisms in the microRNA-binding site of the interleukin 13 gene (IL13) with the risk and prognosis of CRC. Three polymorphisms (rs847, rs848, and rs1295685) were selected for genotyping in a case-control study (513 cases, 572 controls), and 386 CRC patients were followed up. Two dietary factors closely related with gut microbiota (allium vegetables, overnight meal) were significantly associated with CRC development. Although the three SNPs showed no statistically significant associations with the risk and prognosis of CRC, a significant antagonistic interaction was found between rs848 (G-T) and allium vegetable intake (ORi (odds ratio of interaction), 0.92; 95% CI (confidence interval): 0.86, 0.99; P = 0.03); moreover, significant combined and synergistic interactions were observed for all three SNPs and overnight meal intake. This is the first report of significant combined and interactive effects between dietary factors and polymorphisms in the microRNA binding site of IL13 in CRC and may provide direct guidance on intake of allium vegetable and overnight meals for individuals with specific genetic variants of IL13 to modify their susceptibility to CRC.

Genetic variants in microRNAs and their binding sites within gene 3'UTRs associate with susceptibility to age-related macular degeneration.

Age-related macular degeneration (AMD), the leading cause of blindness in the elderly, is a complex disease that results from multiple genetic and environmental factors. MicroRNAs (miRNAs) are small noncoding RNAs that post-transcriptionally regulate target mRNAs and are frequently implicated in human diseases. Here, we investigated the association of genetic variants in miRNAs and miRNA-binding sites within gene 3'-untranslated regions (3'UTRs) with AMD using data from the largest AMD genome-wide association study. First, we identified three variants in miRNAs significantly associated with AMD. These include rs2168518:G>A in the miR-4513 seed sequence, rs41292412:C>T in pre-miR-122/miR-3591, and rs4351242:C>T in the terminal-loop of pre-miR-3135b. We demonstrated that these variants reduce expression levels of the mature miRNAs in vitro and pointed the target genes that may mediate downstream effects of these miRNAs in AMD. Second, we identified 54 variants (in 31 genes) in miRNA-binding sites associated with AMD. Based on stringent prioritization criteria, we highlighted the variants that are more likely to have an impact on the miRNA-target interactions. Further, we selected rs4151672:C>T within the CFB 3'UTR and experimentally showed that while miR-210-5p downregulates expression of CFB, the variant decreases miR-210-5p-mediated repression of CFB. Together, our findings support the notion that miRNAs may play a role in AMD.

rs15869 at miRNA binding site in BRCA2 is associated with breast cancer susceptibility.

BRCA1 and BRCA2 mutations confer an increased lifetime risk of breast cancer; however, the associations of microRNA (miRNA) binding site single nucleotide polymorphisms (SNPs) in 3' untranslated region (3'-UTR) of BRCA1 and BRCA2 with breast cancer (BC) risk were rarely reported. In this case-control study (498 BC patients and 498 matched controls), three SNPs (rs8176318, rs12516 and rs15869) were selected in the 3'-UTR of BRCA1 and BRCA2 genes, which were within miRNA-binding seed regions and might have potential function to regulate the expression of BRCA1/BRCA2. Unconditional logistic regression model was used to analyze the association between three SNPs and BC risk with adjustment of reproductive factors, and Student's t test was performed to assess relative expression of BRCA2 in human breast cancer cell lines. Multifactor dimensionality reduction method was applied to calculate gene-reproductive factors interactions. A novel finding showed that AC [odds ratio (OR) 1.524; 95% confidence interval (CI) 1.141-2.035] genotype of rs15869 in BRCA2 could increase the risk of BC and recombinant plasmid-pGenesil-1-miR-627 could negatively regulate the expression of BRCA2 in MCF-7 and MDA-MB-231 cells. Gene-reproductive factors interactions analysis revealed that rs15869 together with age at menarche and number of pregnancy could increase the risk of BC by 2.39-fold and TT genotype (OR 0.316; 95% CI 0.130-0.767) of rs8176318 had a significant association with progesterone receptor status in BC patients. Our findings suggest that the miRNA-binding SNPs in BRCA1/BRCA2 and their interaction with reproductive factors might contribute to BC risk, and miR-627 might down-regulate BRCA2 expression in MCF-7 and MDA-MB-231 cells.

MiRNA-binding site functional polymorphisms in DNA repair genes RAD51, RAD52, and XRCC2 and breast cancer risk in Chinese population.

RAD51, RAD52, and XRCC2 are all involved in DNA homologous recombinational repair, and there are interactions among those genes. Polymorphisms in 3'-UTR of DNA repair genes may change DNA repair capacity by regulating gene expression. However, potential regulatory variants affecting their expression remain largely unexplored. Five miRNA-binding site SNPs (rs7180135 and rs45549040 in RAD51, rs1051669 and rs7963551 in RAD52 and rs3218550 in XRCC2) selected by bioinformatics method were genotyped in 498 breast cancer (BC) patients and 498 matched controls in Chinese population. Association between SNPs and BC risk was analyzed by adjusted odds ratios (ORs) and 95 % confidence intervals (CIs) in unconditional logistic regression model. Quantitative real-time (qRT) PCR and Western Blot assays were used to calculate the relative expression of RAD52 in recombinant plasmid-pGenesil-1-let-7b group and let-7b-inhibitor group. Gene-reproductive factors interactions were evaluated by multifactor dimensionality reduction (MDR) method. We found that individuals with AC (OR 0.684, 95%CI 0.492-0.951) and CC (OR 0.317, 95%CI 0.200-0.503) genotypes of rs7963551 had a significantly lower risk of breast cancer and qRT-PCR and Western Blot revealed that let-7b might downregulate the expression of RAD52 in MCF-7 and SKBR-3 cells. A significant interaction between the number of pregnancy (≥2) and rs7963551 (Ars7963551) was found to increase breast cancer risk by 2.63-fold (OR 2.63; 95%CI 2.03-3.42). In summary, the miRNA-binding SNPs in DNA repair genes RAD51, RAD52, and XRCC2 and their interaction with reproductive factors might play important roles in the development of BC, and let-7b might downregulate RAD52 expression in MCF-7 and SKBR-3 cells.

STarMir Tools for Prediction of microRNA Binding Sites.

MicroRNAs (miRNAs) are a class of endogenous short noncoding RNAs that regulate gene expression by targeting messenger RNAs (mRNAs), which results in translational repression and/or mRNA degradation. As regulatory molecules, miRNAs are involved in many mammalian biological processes and also in the manifestation of certain human diseases. As miRNAs play central role in the regulation of gene expression, understanding miRNA-binding patterns is essential to gain an insight of miRNA mediated gene regulation and also holds promise for therapeutic applications. Computational prediction of miRNA binding sites on target mRNAs facilitates experimental investigation of miRNA functions. This chapter provides protocols for using the STarMir web server for improved predictions of miRNA binding sites on a target mRNA. As an application module of the Sfold RNA package, the current version of STarMir is an implementation of logistic prediction models developed with high-throughput miRNA binding data from cross-linking immunoprecipitation (CLIP) studies. The models incorporated comprehensive thermodynamic, structural, and sequence features, and were found to make improved predictions of both seed and seedless sites, in comparison to the established algorithms (Liu et al., Nucleic Acids Res 41:e138, 2013). Their broad applicability was indicated by their good performance in cross-species validation. STarMir is freely available at http://sfold.wadsworth.org/starmir.html .

Genetic variants within microRNA-binding site of RAD51B are associated with risk of cervical cancer in Chinese women.

RAD51B plays a central role in homologous recombinational repair (HRR) of DNA double-strand breaks (DSBs), which is important to prevent genomic instability, a hallmark of cancer. Recent studies suggested that common genetic variants of RAD51B may contribute to cancer susceptibility. In this study, we aimed to investigate whether potentially functional variants within miRNA-binding sites of RAD51B are associated with risk of cervical cancer. A total of 1486 cervical cancer patients and 1536 cancer-free controls were enrolled, and two genetic variants, rs963917 (A > G) and rs963918 (T > C), were genotyped in all participants. Using multivariate logistic regression analyses, we found that G allele of rs963917 conferred lower risk of cervical cancer compared to A allele (adjusted OR = 0.89, 95% CI = 0.80-0.99, P = 0.039). Similarly, rs963918 allele C was associated with a decreased risk for cervical cancer compared with allele T (adjusted OR = 0.84, 95% CI = 0.74-0.94, P = 0.004). Haplotype analyses showed that haplotype GC was also correlated with lower risk (OR = 0.83, 95% CI = 0.73-0.95, P = 0.005) compared with the most common haplotype AT. In summary, our study suggested that miRNA-binding site genetic variants of RAD51B may modify the susceptibility to cervical cancer, which is important to identify individuals with differential risk for this malignancy and to improve the effectiveness of preventive intervention.

Genetic Variants in MicroRNAs and Their Binding Sites Are Associated with the Risk of Parkinson Disease.

MicroRNAs (miRNAs) are small noncoding RNAs that serve as key regulators of gene expression. They have been shown to be involved in a wide range of biological processes including neurodegenerative diseases. Genetic variants in miRNAs or miRNA-binding sites on their target genes could affect miRNA function and contribute to disease risk. Here, we investigated the association of miRNA-related genetic variants with Parkinson disease (PD) using data from the largest GWAS on PD. Of 243 miRNA variants, we identified rs897984:T>C in miR-4519 (P value = 1.3×10(-5) and OR = 0.93) and rs11651671:A>G in miR-548at-5p (P value = 1.1×10(-6) and OR = 1.09) to be associated with PD. We showed that the variant's mutant alleles change the secondary structure and decrease expression level of their related miRNAs. Subsequently, we highlighted target genes that might mediate the effects of miR-4519 and miR-548at-5p on PD. Among them, we experimentally showed that NSF is a direct target of miR-4519. Furthermore, among 48,844 miRNA-binding site variants, we found 32 variants (within 13 genes) that are associated with PD. Four of the host genes, CTSB, STX1B, IGSF9B, and HSD3B7, had not previously been reported to be associated with PD. We provide evidence supporting the potential impact of the identified miRNA-binding site variants on miRNA-mediated regulation of their host genes.