Search

It is a well-known and intensively studied phenomenon that the levels of many miRNAs are differentiated in cancer. miRNA biogenesis and functional expression are complex processes orchestrated by many proteins cumulatively called miRNA biogenesis proteins. To characterize cancer somatic mutations in the miRNA biogenesis genes and investigate their potential impact on the levels of miRNAs, we analyzed whole-exome sequencing datasets of over 10,000 cancer/normal sample pairs deposited within the TCGA repository. We identified and characterized over 3,600 somatic mutations in 29 miRNA biogenesis genes and showed that some of the genes are overmutated in specific cancers and/or have recurrent hotspot mutations (e.g., SMAD4 in PAAD, COAD, and READ; DICER1 in UCEC; PRKRA in OV; and LIN28B in SKCM). We identified a list of miRNAs whose level is affected by particular types of mutations in either SMAD4, SMAD2, or DICER1 and showed that hotspot mutations in the RNase domains in DICER1 not only decrease the level of 5p-miRNAs but also increase the level of 3p-miRNAs, including many well-known cancer-related miRNAs. We also showed an association of the mutations with patient survival. Eventually, we created an atlas/compendium of miRNA biogenesis alterations providing a useful resource for different aspects of biomedical research.

Spinocerebellar ataxia type 3 (SCA3) is a progressive neurodegenerative disorder caused by a CAG repeat expansion in the ATXN3 gene encoding the ataxin-3 protein. Despite extensive research the exact pathogenic mechanisms of SCA3 are still not understood in depth. In the present study, to gain insight into the toxicity induced by the expanded CAG repeats in SCA3, we comprehensively investigated repeat-associated non-ATG (RAN) translation in various cellular models expressing translated or non-canonically translated ATXN3 sequences with an increasing number of CAG repeats. We demonstrate that two SCA3 RAN proteins, polyglutamine (polyQ) and polyalanine (polyA), are found only in the case of CAG repeats of pathogenic length. Despite having distinct cellular localization, RAN polyQ and RAN polyA proteins are very often coexpressed in the same cell, impairing nuclear integrity and inducing apoptosis. We provide for the first time mechanistic insights into SCA3 RAN translation indicating that ATXN3 sequences surrounding the repeat region have an impact on SCA3 RAN translation initiation and efficiency. We revealed that RAN translation of polyQ proteins starts at non-cognate codons upstream of the CAG repeats, whereas RAN polyA proteins are likely translated within repeats. Furthermore, integrated stress response activation enhances SCA3 RAN translation. Our findings suggest that the ATXN3 sequence context plays an important role in triggering SCA3 RAN translation and that SCA3 RAN proteins may cause cellular toxicity.

Spinocerebellar ataxia type 3 (SCA3) is a progressive neurodegenerative disorder caused by a CAG repeat expansion in theATXN3gene encoding the ataxin-3 protein. Despite extensive research the exact pathogenic mechanisms of SCA3 are still not understood in depth. In the present study, to gain insight into the toxicity induced by the expanded CAG repeats in SCA3, we comprehensively investigated repeat-associated non-ATG (RAN) translation in various cellular models expressing translated or non-canonically translatedATXN3sequences with an increasing number of CAG repeats. We demonstrate that two SCA3 RAN proteins, polyglutamine (polyQ) and polyalanine (polyA), are found only in the case of CAG repeats of pathogenic length. Despite having distinct cellular localization, RAN polyQ and RAN polyA proteins are very often coexpressed in the same cell, impairing nuclear integrity and inducing apoptosis. We provide for the first time mechanistic insights into SCA3 RAN translation indicating thatATXN3sequences surrounding the repeat region have an impact on SCA3 RAN translation initiation and efficiency. We revealed that RAN translation of polyQ proteins starts at non-cognate codons upstream of the CAG repeats, whereas RAN polyA proteins are likely translated within repeats. Furthermore, integrated stress response activation enhances SCA3 RAN translation. We suggest that RAN translation in SCA3 is a common event substantially contributing to SCA3 pathogenesis and that theATXN3sequence context plays an important role in triggering this unconventional translation.

A growing body of evidence indicates that miRNAs may either drive or suppress oncogenesis. However, little is known about somatic mutations in miRNA genes. To determine the frequency and potential consequences of miRNA gene mutations, we analyzed whole exome sequencing datasets of 569 lung adenocarcinoma (LUAD) and 597 lung squamous cell carcinoma (LUSC) samples generated in The Cancer Genome Atlas (TCGA) project. Altogether, we identified 1091 somatic sequence variants affecting 522 different miRNA genes and showed that half of all cancers had at least one such somatic variant/mutation. These sequence variants occurred in most crucial parts of miRNA precursors, including mature miRNA and seed sequences. Due to our findings, we hypothesize that seed mutations may affect miRNA:target interactions, drastically changing the pool of predicted targets. Mutations may also affect miRNA biogenesis by changing the structure of miRNA precursors, DROSHA and DICER cleavage sites, and regulatory sequence/structure motifs. We identified 10 significantly overmutated hotspot miRNA genes, including the miR-379 gene in LUAD enriched in mutations in the mature miRNA and regulatory sequences. The occurrence of mutations in the hotspot miRNA genes was also shown experimentally. We present a comprehensive analysis of somatic variants in miRNA genes and show that some of these genes are mutational hotspots, suggesting their potential role in cancer.

A growing body of evidence indicates that miRNAs may either drive or suppress oncogenesis. However, little is known about somatic mutations in miRNA genes. To determine the frequency and potential consequences of miRNA gene mutations, we analyzed whole exome sequencing datasets of ∼500 lung adenocarcinoma (LUAD) and ∼500 lung squamous cell carcinoma (LUSC) samples generated in the TCGA. Altogether, we identified >1000 mutations affecting ∼500 different miRNA genes and showed that half of all cancers had at least one such mutation. Mutations occurred in most crucial parts of miRNA precursors, including mature miRNA and seed sequences. We showed that seed mutations strongly affected miRNA:target interactions, drastically changing the pool of predicted targets. Mutations may also affect miRNA biogenesis by changing the structure of miRNA precursors, DROSHA and DICER cleavage sites, and regulatory sequence/structure motifs. We identified 10 significantly overmutated hotspot miRNA genes, including themiR-379gene in LUAD enriched in mutations in the mature miRNA and regulatory sequences. The occurrence of mutations in the hotspot miRNA genes was also shown experimentally. We present a comprehensive analysis of somatic mutations in miRNA genes and show that some of these genes are mutational hotspots, suggesting their potential role in cancer.

Background miRNAs are considered important players in oncogenesis, serving either as oncomiRs or suppressormiRs. Although the accumulation of somatic alterations is an intrinsic aspect of cancer development and many important cancer-driving mutations have been identified in protein-coding genes, the area of functional somatic mutations in miRNA genes is heavily understudied. Methods Here, based on the analysis of large genomic datasets, mostly the whole-exome sequencing of over 10,000 cancer/normal sample pairs deposited within the TCGA repository, we undertook an analysis of somatic mutations in miRNA genes. Findings We identified and characterized over 10,000 somatic mutations and showed that some of the miRNA genes are overmutated in Pan-Cancer and/or specific cancers. Nonrandom occurrence of the identified mutations was confirmed by a strong association of overmutated miRNA genes with KEGG pathways, most of which were related to specific cancer types or cancer-related processes. Additionally, we showed that mutations in some of the overmutated genes correlate with miRNA expression, cancer staging, and patient survival. Interpretation Our study is the first comprehensive Pan-Cancer study of cancer somatic mutations in miRNA genes. It may help to understand the consequences of mutations in miRNA genes and the identification of miRNA functional mutations. The results may also be the first step (form the basis and provide the resources) in the development of computational and/or statistical approaches/tools dedicated to the identification of cancer-driver miRNA genes. Funding This work was supported by research grants from the Polish National Science Centre 2016/22/A/NZ2/00184 and 2015/17/N/NZ3/03629.