Prediction and pan-cancer analysis of mammalian transcripts involved in target directed miRNA degradation.

It is currently unknown how many RNA transcripts are able to induce degradation of microRNAs (miRNA) via the mechanism known as target-directed miRNA degradation (TDMD). We developed TDMDfinder, a computational pipeline that identifies 'high confidence' TDMD interactions in the Human and Mouse transcriptomes by combining sequence alignment and feature selection approaches. Our predictions suggested that TDMD is widespread, with potentially every miRNA controlled by endogenous targets. We experimentally tested 37 TDMDfinder predictions, of which 17 showed TDMD effects as measured by RT-qPCR and small RNA sequencing, linking the miR-17, miR-19, miR-30, miR-221, miR-26 and miR-23 families to novel endogenous TDMDs. In some cases, TDMD was found to affect different members of the same miRNA family selectively. Features like complementarity to the miRNA 3' region, bulge size and hybridization energy appeared to be the main factors determining sensitivity. Computational analyses performed using the multiomic TCGA platform substantiated the involvement of many TDMD transcripts in human cancer and highlighted 36 highly significant interactions, suggesting TDMD as a new potential oncogenic mechanism. In conclusion, TDMDfinder provides the first inventory of bona fide human and mouse TDMDs. Available as a free webtool, TDMDfinder allows users to search for any TDMD interaction of interest by customizing its selection criteria.

A novel microRNA signature for the detection of melanoma by liquid biopsy.

BACKGROUND: Melanoma is the deadliest form of skin cancer and metastatic disease is associated with a significant survival rate drop. There is an urgent need for consistent tumor biomarkers to scale precision medicine and reduce cancer mortality. Here, we aimed to identify a melanoma-specific circulating microRNA signature and assess its value as a diagnostic tool. METHODS: The study consisted of a discovery phase and two validation phases. Circulating plasma extracellular vesicles (pEV) associated microRNA profiles were obtained from a discovery cohort of metastatic melanoma patients and normal subjects as controls. A pEV-microRNA signature was obtained using a LASSO penalized logistic regression model. The pEV-microRNA signature was subsequently validated both in a publicly available dataset and in an independent internal cohort. RESULTS: We identified and validated in three independent cohorts a panel of melanoma-specific circulating microRNAs that showed high accuracy in differentiating melanoma patients from healthy subjects with an area under the curve (AUC) of 1.00, 0.94 and 0.75 respectively. Investigation of the function of the pEV-microRNA signature evidenced their possible immune suppressive role in melanoma patients. CONCLUSIONS: We demonstrate that a blood test based on circulating microRNAs can non-invasively detect melanoma, offering a novel diagnostic tool for improving standard care. Moreover, we revealed an immune suppressive role for melanoma pEV-microRNAs.

Deep learning predicts short non-coding RNA functions from only raw sequence data.

Small non-coding RNAs (ncRNAs) are short non-coding sequences involved in gene regulation in many biological processes and diseases. The lack of a complete comprehension of their biological functionality, especially in a genome-wide scenario, has demanded new computational approaches to annotate their roles. It is widely known that secondary structure is determinant to know RNA function and machine learning based approaches have been successfully proven to predict RNA function from secondary structure information. Here we show that RNA function can be predicted with good accuracy from a lightweight representation of sequence information without the necessity of computing secondary structure features which is computationally expensive. This finding appears to go against the dogma of secondary structure being a key determinant of function in RNA. Compared to recent secondary structure based methods, the proposed solution is more robust to sequence boundary noise and reduces drastically the computational cost allowing for large data volume annotations. Scripts and datasets to reproduce the results of experiments proposed in this study are available at: https://github.com/bioinformatics-sannio/ncrna-deep.

Structural Model for Recruitment of RIT1 to the LZTR1 E3 Ligase: Evidences from an Integrated Computational Approach.

Leucine-zipper transcription regulator 1 (LZTR1) is a highly mutated tumor suppressor gene, involved in the pathogenesis of several cancer types and developmental disorders. In proteasomal degradation, it acts as an adaptor protein responsible for the recognition and recruitment of substrates to be ubiquitinated in Cullin3-RING ligase E3 (CRL3) machinery. LZTR1 belongs to the BTB-Kelch family, a multi-domain protein where the Kelch propeller plays as the substrate recognition region and for which no experimental structure has been solved. Recently, large effort mutational analyses pointed to the role of disease-associated LZTR1 mutations in the RAS/MAPK signaling pathway and RIT1, a small Ras-related GTPase protein, has been identified by mass spectroscopy to interact with LZTR1. Hence, a better understanding of native structure, molecular mechanism, and substrate specificity would help clarifying the role of LZTR1 in pathological diseases, thus promoting advancement in the development of novel therapeutic strategies. Here, we address the interaction model between adaptor LZTR1 and substrate RIT1 by applying an integrated computational approach, including molecular modeling and docking techniques. We observe that the interaction model LZTR1-RIT1 is stabilized by an electrostatic bond network established between the two protein surfaces, which is reminiscent of homologous ubiquitin ligases complexes. Then, running MD simulations, we characterize differential conformational dynamics of the multi-domain LZTR1, offering interesting implications on the mechanistic role of specific point mutations. We identify G248R and R283Q as damaging mutations involved in the recognition process of the substrate RIT1 and R412C as a possible allosteric mutation from the Kelch to the C-term BTB-domain. Our findings provide important structural insights on targeting CRL3s for drug discovery.

MiR-449a antagonizes EMT through IL-6-mediated trans-signaling in laryngeal squamous cancer.

MicroRNAs (miRNAs) are involved in post-transcriptional gene expression regulation and in mechanisms of cancer growth and metastases. In this light, miRNAs could be promising therapeutic targets and biomarkers in clinical practice. Therefore, we investigated if specific miRNAs and their target genes contribute to laryngeal squamous cell carcinoma (LSCC) development. We found a significant decrease of miR-449a in LSCC patients with nodal metastases (63.3%) compared with patients without nodal involvement (44%). The AmpliSeq Transcriptome of HNO-210 miR-449a-transfected cell lines allowed the identification of IL6-R as a potential target. Moreover, the downregulation of IL6-R and the phosphorylation reduction of the downstream signaling effectors, suggested the inhibition of the IL-6 trans-signaling pathway. These biochemical effects were paralleled by a significant inhibition of invasion and migration in vitro and in vivo, supporting an involvement of epithelial-mesenchymal transition. These findings indicate that miR-449a contributes to suppress the metastasization of LSCC by the IL-6 trans-signaling block and affects sensitivity to external stimuli that mimic pro-inflammatory conditions.

Guadecitabine plus ipilimumab in unresectable melanoma: five-year follow-up and integrated multi-omic analysis in the phase 1b NIBIT-M4 trial.

Association with hypomethylating agents is a promising strategy to improve the efficacy of immune checkpoint inhibitors-based therapy. The NIBIT-M4 was a phase Ib, dose-escalation trial in patients with advanced melanoma of the hypomethylating agent guadecitabine combined with the anti-CTLA-4 antibody ipilimumab that followed a traditional 3 + 3 design (NCT02608437). Patients received guadecitabine 30, 45 or 60 mg/m2/day subcutaneously on days 1 to 5 every 3 weeks starting on week 0 for a total of four cycles, and ipilimumab 3 mg/kg intravenously starting on day 1 of week 1 every 3 weeks for a total of four cycles. Primary outcomes of safety, tolerability, and maximum tolerated dose of treatment were previously reported. Here we report the 5-year clinical outcome for the secondary endpoints of overall survival, progression free survival, and duration of response, and an exploratory integrated multi-omics analysis on pre- and on-treatment tumor biopsies. With a minimum follow-up of 45 months, the 5-year overall survival rate was 28.9% and the median duration of response was 20.6 months. Re-expression of immuno-modulatory endogenous retroviruses and of other repetitive elements, and a mechanistic signature of guadecitabine are associated with response. Integration of a genetic immunoediting index with an adaptive immunity signature stratifies patients/lesions into four distinct subsets and discriminates 5-year overall survival and progression free survival. These results suggest that coupling genetic immunoediting with activation of adaptive immunity is a relevant requisite for achieving long term clinical benefit by epigenetic immunomodulation in advanced melanoma patients.