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.

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.

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 3600 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.

175)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. Here, based on analysis of the whole-exome sequencing of over 10,000 cancer/normal sample pairs deposited within the TCGA repository, we identified and characterized over 10,000 somatic mutations in miRNA genes and showed that some of the 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. Our 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 cancerdriver miRNA genes. The results published here are based upon data generated by the TCGA Research Network: https://www.cancer.gov/tcga. This work was supported by research grants from the Polish National Science Centre [2016/22/A/NZ2/00184 (to P.K.) and 2015/17/N/NZ3/03629 (to M.O.U-T.)] Authors contributions MUTparticipated in conceiving the study, wrote all the scripts, performed most of the computational and statistical analyses, drafted the manuscript, prepared figures, tables, and supplementary materials; PGMparticipated in conceiving the study, discussed the study on all steps of analyses, participated in manuscript preparation, performed some analyses (3D structures); PNdiscussed the analyses on all steps of the study, performed some analyses, participated in manuscript preparation; EKperformed the sequencing experiments, critically read and corrected manuscript; SSperformed the sequencing experiments, critically read and corrected manuscript; MGprovided samples for the hsa-miR-1324 and hsa-miR-142 sequencing validation experiments, advised on hematological neoplasms and head and neck cancer, critically read and corrected manuscript; PKreceived financing, participated in conceiving the study, supervised and coordinated the study, drafted the manuscript (with MUT).

Homologous recombination (HR) faithfully restores DNA double-strand breaks. Defects in this HR repair pathway are associated with cancer predisposition. In genetic engineering, HR has been used extensively to study gene function and it represents an ideal method of gene therapy for single gene disorders. Here, we present a novel assay to measure HR in living cells. The HR substrate consisted of a non-fluorescent 3’ truncated form of the eGFP gene and was integrated into the AAVS1 locus, known as a safe harbor. The donor DNA template comprised a 5’ truncated eGFP copy and was delivered via AAV particles. HR mediated repair restored full-length eGFP coding sequence, resulting in eGFP+ cells. The utility of our assay in quantifying HR events was validated by exploring the impact of the overexpression of HR promoters and the siRNA-mediated silencing of genes known to play a role in DNA repair on the frequency of HR. We conclude that this novel assay represents a useful tool to further investigate the mechanisms that control HR and test continually emerging tools for HR-mediated genome editing.