We demonstrated that Hsa-miR-567 expression is significantly downregulated in poor prognosis breast cancer, compared to better prognosis breast cancer, having a role in the control of cell proliferation and migration by regulating KPNA4 gene.

Modern high-throughput genomic technologies represent a comprehensive hallmark of molecular changes in pan-cancer studies. Although different cancer gene signatures have been revealed, the mechanism of tumourigenesis has yet to be completely understood. Pathways and networks are important tools to explain the role of genes in functional genomic studies. However, few methods consider the functional non-equal roles of genes in pathways and the complex gene-gene interactions in a network.

Cancer progression involves the gradual loss of a differentiated phenotype and acquisition of progenitor and stem-cell-like features. Here, we provide novel stemness indices for assessing the degree of oncogenic dedifferentiation. We used an innovative one-class logistic regression (OCLR) machine-learning algorithm to extract transcriptomic and epigenetic feature sets derived from non-transformed pluripotent stem cells and their differentiated progeny. Using OCLR, we were able to identify previously undiscovered biological mechanisms associated with the dedifferentiated oncogenic state. Analyses of the tumor microenvironment revealed unanticipated correlation of cancer stemness with immune checkpoint expression and infiltrating immune cells. We found that the dedifferentiated oncogenic phenotype was generally most prominent in metastatic tumors. Application of our stemness indices to single-cell data revealed patterns of intra-tumor molecular heterogeneity. Finally, the indices allowed for the identification of novel targets and possible targeted therapies aimed at tumor differentiation.

Like other cancer diseases, prostate cancer (PC) is caused by the accumulation of genetic alterations in the cells that drives malignant growth. These alterations are revealed by gene profiling and copy number alteration (CNA) analysis. Moreover, recent evidence suggests that also microRNAs have an important role in PC development. Despite efforts to profile PC, the alterations (gene, CNA, and miRNA) and biological processes that correlate with disease development and progression remain partially elusive. Many gene signatures proposed as diagnostic or prognostic tools in cancer poorly overlap. The identification of co-expressed genes, that are functionally related, can identify a core network of genes associated with PC with a better reproducibility. By combining different approaches, including the integration of mRNA expression profiles, CNAs, and miRNA expression levels, we identified a gene signature of four genes overlapping with other published gene signatures and able to distinguish, in silico, high Gleason-scored PC from normal human tissue, which was further enriched to 19 genes by gene co-expression analysis. From the analysis of miRNAs possibly regulating this network, we found that hsa-miR-153 was highly connected to the genes in the network. Our results identify a four-gene signature with diagnostic and prognostic value in PC and suggest an interesting gene network that could play a key regulatory role in PC development and progression. Furthermore, hsa-miR-153, controlling this network, could be a potential biomarker for theranostics in high Gleason-scored PC.

MicroRNAs (miRNAs) are small noncoding RNAs that have emerged as new potential epigenetic biomarkers. Here, we evaluate the efficacy of six circulating miRNA previously described in the literature as biomarkers for the diagnosis of temporal lobe epilepsy (TLE) and/or as predictive biomarkers to antiepileptic drug response. We measured the differences in serum miRNA levels by quantitative reverse transcriptase polymerase chain reaction (qRT-PCR) assays in a cohort of 27 patients (14 women and 13 men; mean ± SD age: 43.65 ± 17.07) with TLE compared to 20 healthy controls (HC) matched for sex, age and ethnicity (11 women and 9 men; mean ± SD age: 47.5 ± 9.1). Additionally, patients were classified according to whether they had drug-responsive (n = 17) or drug-resistant (n = 10) TLE. We have investigated any correlations between miRNAs and several electroclinical parameters. Three miRNAs (miR-142, miR-146a, miR-223) were significantly upregulated in patients (expressed as average expression ± SD). In detail, miR-142 expression was 0.40 ± 0.29 versus 0.16 ± 0.10 in TLE patients compared to HC (t-test, p < 0.01), miR-146a expression was 0.15 ± 0.11 versus 0.07 ± 0.04 (t-test, p < 0.05), and miR-223 expression was 6.21 ± 3.65 versus 1.23 ± 0.84 (t-test, p < 0.001). Moreover, results obtained from a logistic regression model showed the good performance of miR-142 and miR-223 in distinguishing drug-sensitive vs. drug-resistant TLE. The results of this pilot study give evidence that miRNAs are suitable targets in TLE and offer the rationale for further confirmation studies in larger epilepsy cohorts.

Cancer driver gene alterations influence cancer development, occurring in oncogenes, tumor suppressors, and dual role genes. Discovering dual role cancer genes is difficult because of their elusive context-dependent behavior. We define oncogenic mediators as genes controlling biological processes. With them, we classify cancer driver genes, unveiling their roles in cancer mechanisms. To this end, we present Moonlight, a tool that incorporates multiple -omics data to identify critical cancer driver genes. With Moonlight, we analyze 8000+ tumor samples from 18 cancer types, discovering 3310 oncogenic mediators, 151 having dual roles. By incorporating additional data (amplification, mutation, DNA methylation, chromatin accessibility), we reveal 1000+ cancer driver genes, corroborating known molecular mechanisms. Additionally, we confirm critical cancer driver genes by analysing cell-line datasets. We discover inactivation of tumor suppressors in intron regions and that tissue type and subtype indicate dual role status. These findings help explain tumor heterogeneity and could guide therapeutic decisions.

Previous studies reported that Angiotensin converting enzyme 2 (ACE2) is the main cell receptor of SARS-CoV and SARS-CoV-2. It plays a key role in the access of the virus into the cell to produce the final infection. In the present study we investigated in silico the basic mechanism of ACE2 in the lung and provided evidences for new potentially effective drugs for Covid-19. Specifically, we used the gene expression profiles from public datasets including The Cancer Genome Atlas, Gene Expression Omnibus and Genotype-Tissue Expression, Gene Ontology and pathway enrichment analysis to investigate the main functions of ACE2-correlated genes. We constructed a protein-protein interaction network containing the genes co-expressed with ACE2. Finally, we focused on the genes in the network that are already associated with known drugs and evaluated their role for a potential treatment of Covid-19. Our results demonstrate that the genes correlated with ACE2 are mainly enriched in the sterol biosynthetic process, Aryldialkylphosphatase activity, adenosylhomocysteinase activity, trialkylsulfonium hydrolase activity, acetate-CoA and CoA ligase activity. We identified a network of 193 genes, 222 interactions and 36 potential drugs that could have a crucial role. Among possible interesting drugs for Covid-19 treatment, we found Nimesulide, Fluticasone Propionate, Thiabendazole, Photofrin, Didanosine and Flutamide.

We undertook a comprehensive proteogenomic characterization of 95 prospectively collected endometrial carcinomas, comprising 83 endometrioid and 12 serous tumors. This analysis revealed possible new consequences of perturbations to the p53 and Wnt/β-catenin pathways, identified a potential role for circRNAs in the epithelial-mesenchymal transition, and provided new information about proteomic markers of clinical and genomic tumor subgroups, including relationships to known druggable pathways. An extensive genome-wide acetylation survey yielded insights into regulatory mechanisms linking Wnt signaling and histone acetylation. We also characterized aspects of the tumor immune landscape, including immunogenic alterations, neoantigens, common cancer/testis antigens, and the immune microenvironment, all of which can inform immunotherapy decisions. Collectively, our multi-omic analyses provide a valuable resource for researchers and clinicians, identify new molecular associations of potential mechanistic significance in the development of endometrial cancers, and suggest novel approaches for identifying potential therapeutic targets.

Prostate cancer (PC) is a male common neoplasm and is the second leading cause of cancer death in American men. PC is traditionally diagnosed by the evaluation of prostate secreted antigen (PSA) in the blood. Due to the high levels of false positives, digital rectal examination and transrectal ultrasound guided biopsy are necessary in uncertain cases with elevated PSA levels. Nevertheless, the high mortality rate suggests that new PC biomarkers are urgently needed to help clinical diagnosis. In a previous study, we have identified a network of genes, altered in high Gleason Score (GS) PC (GS ≥ 7), being regulated by miR-153. Until now, no publication has explained the mechanism of action of miR-153 in PC. By in vitro studies, we found that the overexpression of miR-153 in high GS cell lines is required to control cell proliferation, migration and invasion rates, targeting Kruppel-like factor 5 (KLF5). Moreover, miR-153 could be secreted by exosomes and microvesicles in the microenvironment and, once entered into the surrounding tissue, could influence cellular growth. Being upregulated in high GS human PC, miR-153 could be proposed as a circulating biomarker for PC diagnosis.

Dysregulation of microRNAs (miRNAs) is involved in the initiation and progression of several human cancers, including breast cancer (BC), as strong evidence has been found that miRNAs can act as oncogenes or tumor suppressor genes. This review presents the state of the art on the role of miRNAs in the diagnosis, prognosis, and therapy of BC. Based on the results obtained in the last decade, some miRNAs are emerging as biomarkers of BC for diagnosis (i.e., miR-9, miR-10b, and miR-17-5p), prognosis (i.e., miR-148a and miR-335), and prediction of therapeutic outcomes (i.e., miR-30c, miR-187, and miR-339-5p) and have important roles in the control of BC hallmark functions such as invasion, metastasis, proliferation, resting death, apoptosis, and genomic instability. Other miRNAs are of interest as new, easily accessible, affordable, non-invasive tools for the personalized management of patients with BC because they are circulating in body fluids (e.g., miR-155 and miR-210). In particular, circulating multiple miRNA profiles are showing better diagnostic and prognostic performance as well as better sensitivity than individual miRNAs in BC. New miRNA-based drugs are also promising therapy for BC (e.g., miR-9, miR-21, miR34a, miR145, and miR150), and other miRNAs are showing a fundamental role in modulation of the response to other non-miRNA treatments, being able to increase their efficacy (e.g., miR-21, miR34a, miR195, miR200c, and miR203 in combination with chemotherapy).

Genes that encode proteins playing a role in more than one biological process are frequently dependent on their tissue context, and human diseases result from the altered interplay of tissue- and cell-specific processes. In this work, we performed a computational approach that identifies tissue-specific co-expression networks by integrating miRNAs, long-non-coding RNAs, and mRNAs in more than eight thousands of human samples from thirty normal tissue types. Our analysis (1) shows that long-non coding RNAs and miRNAs have a high specificity, (2) confirms several known tissue-specific RNAs, and (3) identifies new tissue-specific co-expressed RNAs that are currently still not described in the literature. Some of these RNAs interact with known tissue-specific RNAs or are crucial in key cancer functions, suggesting that they are implicated in tissue specification or cell differentiation.

An important challenge in cancer biology is to understand the complex aspects of the disease. It is increasingly evident that genes are not isolated from each other and the comprehension of how different genes are related to each other could explain biological mechanisms causing diseases. Biological pathways are important tools to reveal gene interaction and reduce the large number of genes to be studied by partitioning it into smaller paths. Furthermore, recent scientific evidence has proven that a combination of pathways, instead than a single element of the pathway or a single pathway, could be responsible for pathological changes in a cell.

We report a comprehensive proteogenomics analysis, including whole-genome sequencing, RNA sequencing, and proteomics and phosphoproteomics profiling, of 218 tumors across 7 histological types of childhood brain cancer: low-grade glioma (n = 93), ependymoma (32), high-grade glioma (25), medulloblastoma (22), ganglioglioma (18), craniopharyngioma (16), and atypical teratoid rhabdoid tumor (12). Proteomics data identify common biological themes that span histological boundaries, suggesting that treatments used for one histological type may be applied effectively to other tumors sharing similar proteomics features. Immune landscape characterization reveals diverse tumor microenvironments across and within diagnoses. Proteomics data further reveal functional effects of somatic mutations and copy number variations (CNVs) not evident in transcriptomics data. Kinase-substrate association and co-expression network analysis identify important biological mechanisms of tumorigenesis. This is the first large-scale proteogenomics analysis across traditional histological boundaries to uncover foundational pediatric brain tumor biology and inform rational treatment selection.

The Cancer Genome Atlas (TCGA) has catalyzed systematic characterization of diverse genomic alterations underlying human cancers. At this historic junction marking the completion of genomic characterization of over 11,000 tumors from 33 cancer types, we present our current understanding of the molecular processes governing oncogenesis. We illustrate our insights into cancer through synthesis of the findings of the TCGA PanCancer Atlas project on three facets of oncogenesis: (1) somatic driver mutations, germline pathogenic variants, and their interactions in the tumor; (2) the influence of the tumor genome and epigenome on transcriptome and proteome; and (3) the relationship between tumor and the microenvironment, including implications for drugs targeting driver events and immunotherapies. These results will anchor future characterization of rare and common tumor types, primary and relapsed tumors, and cancers across ancestry groups and will guide the deployment of clinical genomic sequencing.

We performed an extensive immunogenomic analysis of more than 10,000 tumors comprising 33 diverse cancer types by utilizing data compiled by TCGA. Across cancer types, we identified six immune subtypes-wound healing, IFN-γ dominant, inflammatory, lymphocyte depleted, immunologically quiet, and TGF-β dominant-characterized by differences in macrophage or lymphocyte signatures, Th1:Th2 cell ratio, extent of intratumoral heterogeneity, aneuploidy, extent of neoantigen load, overall cell proliferation, expression of immunomodulatory genes, and prognosis. Specific driver mutations correlated with lower (CTNNB1, NRAS, or IDH1) or higher (BRAF, TP53, or CASP8) leukocyte levels across all cancers. Multiple control modalities of the intracellular and extracellular networks (transcription, microRNAs, copy number, and epigenetic processes) were involved in tumor-immune cell interactions, both across and within immune subtypes. Our immunogenomics pipeline to characterize these heterogeneous tumors and the resulting data are intended to serve as a resource for future targeted studies to further advance the field.

Triple-negative breast cancer (TNBC) is a collection of biologically diverse cancers characterized by distinct transcriptional patterns, biology, and immune composition. TNBCs subtypes include two basal-like (BL1, BL2), a mesenchymal (M) and a luminal androgen receptor (LAR) subtype. Through a comprehensive analysis of mutation, copy number, transcriptomic, epigenetic, proteomic, and phospho-proteomic patterns we describe the genomic landscape of TNBC subtypes. Mesenchymal subtype tumors display high mutation loads, genomic instability, absence of immune cells, low PD-L1 expression, decreased global DNA methylation, and transcriptional repression of antigen presentation genes. We demonstrate that major histocompatibility complex I (MHC-I) is transcriptionally suppressed by H3K27me3 modifications by the polycomb repressor complex 2 (PRC2). Pharmacological inhibition of PRC2 subunits EZH2 or EED restores MHC-I expression and enhances chemotherapy efficacy in murine tumor models, providing a rationale for using PRC2 inhibitors in PD-L1 negative mesenchymal tumors. Subtype-specific differences in immune cell composition and differential genetic/pharmacological vulnerabilities suggest additional treatment strategies for TNBC.

Biotechnological advances in sequencing have led to an explosion of publicly available data via large international consortia such as The Cancer Genome Atlas (TCGA), The Encyclopedia of DNA Elements (ENCODE), and The NIH Roadmap Epigenomics Mapping Consortium (Roadmap). These projects have provided unprecedented opportunities to interrogate the epigenome of cultured cancer cell lines as well as normal and tumor tissues with high genomic resolution. The Bioconductor project offers more than 1,000 open-source software and statistical packages to analyze high-throughput genomic data. However, most packages are designed for specific data types (e.g. expression, epigenetics, genomics) and there is no one comprehensive tool that provides a complete integrative analysis of the resources and data provided by all three public projects. A need to create an integration of these different analyses was recently proposed. In this workflow, we provide a series of biologically focused integrative analyses of different molecular data. We describe how to download, process and prepare TCGA data and by harnessing several key Bioconductor packages, we describe how to extract biologically meaningful genomic and epigenomic data. Using Roadmap and ENCODE data, we provide a work plan to identify biologically relevant functional epigenomic elements associated with cancer. To illustrate our workflow, we analyzed two types of brain tumors: low-grade glioma (LGG) versus high-grade glioma (glioblastoma multiform or GBM). This workflow introduces the following Bioconductor packages: AnnotationHub, ChIPSeeker, ComplexHeatmap, pathview, ELMER, GAIA, MINET, RTCGAToolbox,  TCGAbiolinks.

Identifying molecular cancer drivers is critical for precision oncology. Multiple advanced algorithms to identify drivers now exist, but systematic attempts to combine and optimize them on large datasets are few. We report a PanCancer and PanSoftware analysis spanning 9,423 tumor exomes (comprising all 33 of The Cancer Genome Atlas projects) and using 26 computational tools to catalog driver genes and mutations. We identify 299 driver genes with implications regarding their anatomical sites and cancer/cell types. Sequence- and structure-based analyses identified >3,400 putative missense driver mutations supported by multiple lines of evidence. Experimental validation confirmed 60%-85% of predicted mutations as likely drivers. We found that >300 MSI tumors are associated with high PD-1/PD-L1, and 57% of tumors analyzed harbor putative clinically actionable events. Our study represents the most comprehensive discovery of cancer genes and mutations to date and will serve as a blueprint for future biological and clinical endeavors.

We present a proteogenomic study of 108 human papilloma virus (HPV)-negative head and neck squamous cell carcinomas (HNSCCs). Proteomic analysis systematically catalogs HNSCC-associated proteins and phosphosites, prioritizes copy number drivers, and highlights an oncogenic role for RNA processing genes. Proteomic investigation of mutual exclusivity between FAT1 truncating mutations and 11q13.3 amplifications reveals dysregulated actin dynamics as a common functional consequence. Phosphoproteomics characterizes two modes of EGFR activation, suggesting a new strategy to stratify HNSCCs based on EGFR ligand abundance for effective treatment with inhibitory EGFR monoclonal antibodies. Widespread deletion of immune modulatory genes accounts for low immune infiltration in immune-cold tumors, whereas concordant upregulation of multiple immune checkpoint proteins may underlie resistance to anti-programmed cell death protein 1 monotherapy in immune-hot tumors. Multi-omic analysis identifies three molecular subtypes with high potential for treatment with CDK inhibitors, anti-EGFR antibody therapy, and immunotherapy, respectively. Altogether, proteogenomics provides a systematic framework to inform HNSCC biology and treatment.

Glioblastoma (GBM) is the most aggressive nervous system cancer. Understanding its molecular pathogenesis is crucial to improving diagnosis and treatment. Integrated analysis of genomic, proteomic, post-translational modification and metabolomic data on 99 treatment-naive GBMs provides insights to GBM biology. We identify key phosphorylation events (e.g., phosphorylated PTPN11 and PLCG1) as potential switches mediating oncogenic pathway activation, as well as potential targets for EGFR-, TP53-, and RB1-altered tumors. Immune subtypes with distinct immune cell types are discovered using bulk omics methodologies, validated by snRNA-seq, and correlated with specific expression and histone acetylation patterns. Histone H2B acetylation in classical-like and immune-low GBM is driven largely by BRDs, CREBBP, and EP300. Integrated metabolomic and proteomic data identify specific lipid distributions across subtypes and distinct global metabolic changes in IDH-mutated tumors. This work highlights biological relationships that could contribute to stratification of GBM patients for more effective treatment.