A pan-cancer atlas of somatic mutations in miRNA biogenesis genes.

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.

Essential thrombocythaemia progression to the fibrotic phase is associated with a decrease in JAK2 and PDL1 levels.

It has been postulated that the changes in the molecular characteristics of the malignant clone(s) and the abnormal activation of JAK-STAT signaling are responsible for myeloproliferative neoplasm progression to more advanced disease phases and the immune escape of the malignant clone. The continuous JAK-STAT pathway activation leads to enhanced activity of the promoter of CD274 coding programmed death-1 receptor ligand (PD-L1), increased PD-L1 level, and the immune escape of MPN cells. The aim of study was to evaluate the PDL1 mRNA and JAK2 mRNA level in molecularly defined essential thrombocythaemia (ET) patients (pts) during disease progression to post-ET- myelofibrosis (post-ET-MF). The study group consisted of 162 ET pts, including 30 pts diagnosed with post-ET-MF. The JAK2V617F, CALR, and MPL mutations were found in 59.3%, 19.1%, and 1.2% of pts, respectively. No copy-number alternations of the JAK2, PDL1, and PDCDL1G2 (PDL2) genes were found. The level of PD-L1 was significantly higher in the JAK2V617F than in the JAK2WT, CALR mutation-positive, and triple-negative pts. The PD-L1 mRNA level was weakly correlated with both the JAK2V617F variant allele frequency (VAF), and with the JAK2V617F allele mRNA level. The total JAK2 level in post-ET-MF pts was lower than in ET pts, despite the lack of differences in the JAK2V617F VAF. In addition, the PD-L1 level was lower in post-ET-MF. A detailed analysis has shown that the decrease in JAK2 and PDL1 mRNA levels depended on the bone marrow fibrosis grade. The PDL1 expression showed no differences in relation to the genotype of the JAK2 haplotypeGGCC_46/1, hemoglobin concentration, hematocrit value, leukocyte, and platelet counts. The observed drop of the total JAK2 and PDL1 levels during the ET progression to the post-ET-MF may reflect the changes in the JAK2V617F positive clone proliferative potential and the PD-L1 level-related immunosuppressive effect. The above-mentioned hypothesis is supported by The Cancer Genome Atlas (TCGA) data, confirming a strong positive association between CD274 (encoding PD-L1), CXCR3 (encoding CXCR3), and CSF1 (encoding M-CSF) expression levels, and recently published results documenting a drop in the CXCR3 level and circulating M-CSF in patients with post-ET-MF.

Profile of Basal Cell Carcinoma Mutations and Copy Number Alterations - Focus on Gene-Associated Noncoding Variants.

Basal cell carcinoma (BCC) of the skin is the most common cancer in humans, characterized by the highest mutation rate among cancers, and is mostly driven by mutations in genes involved in the hedgehog pathway. To date, almost all BCC genetic studies have focused exclusively on protein-coding sequences; therefore, the impact of noncoding variants on the BCC genome is unrecognized. In this study, with the use of whole-exome sequencing of 27 tumor/normal pairs of BCC samples, we performed an analysis of somatic mutations in both protein-coding sequences and gene-associated noncoding regions, including 5'UTRs, 3'UTRs, and exon-adjacent intron sequences. Separately, in each region, we performed hotspot identification, mutation enrichment analysis, and cancer driver identification with OncodriveFML. Additionally, we performed a whole-genome copy number alteration analysis with GISTIC2. Of the >80,000 identified mutations, ~50% were localized in noncoding regions. The results of the analysis generally corroborated the previous findings regarding genes mutated in coding sequences, including PTCH1, TP53, and MYCN, but more importantly showed that mutations were also clustered in specific noncoding regions, including hotspots. Some of the genes specifically mutated in noncoding regions were identified as highly potent cancer drivers, of which BAD had a mutation hotspot in the 3'UTR, DHODH had a mutation hotspot in the Kozak sequence in the 5'UTR, and CHCHD2 frequently showed mutations in the 5'UTR. All of these genes are functionally implicated in cancer-related processes (e.g., apoptosis, mitochondrial metabolism, and de novo pyrimidine synthesis) or the pathogenesis of UV radiation-induced cancers. We also found that the identified BAD and CHCHD2 mutations frequently occur in melanoma but not in other cancers via The Cancer Genome Atlas analysis. Finally, we identified a frequent deletion of chr9q, encompassing PTCH1, and unreported frequent copy number gain of chr9p, encompassing the genes encoding the immune checkpoint ligands PD-L1 and PD-L2. In conclusion, this study is the first systematic analysis of coding and noncoding mutations in BCC and provides a strong basis for further analyses of the variants in BCC and cancer in general.

Somatic Mutations in miRNA Genes in Lung Cancer-Potential Functional Consequences of Non-Coding Sequence Variants.

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.