Evolution of Metastases in Space and Time under Immune Selection.

We examined how the immune microenvironment molds tumor evolution at different metastatic organs in a longitudinal dataset of colorectal cancer. Through multiplexed analyses, we showed that clonal evolution patterns during metastatic progression depend on the immune contexture at the metastatic site. Genetic evidence of neoantigen depletion was observed in the sites with high Immunoscore and spatial proximity between Ki67+ tumor cells and CD3+ cells. The immunoedited tumor clones were eliminated and did not recur, while progressing clones were immune privileged, despite the presence of tumor-infiltrating lymphocytes. Characterization of immune-privileged metastases revealed tumor-intrinsic and tumor-extrinsic mechanisms of escape. The lowest recurrence risk was associated with high Immunoscore, occurrence of immunoediting, and low tumor burden. We propose a parallel selection model of metastatic progression, where branched evolution could be traced back to immune-escaping clones. The findings could inform the understanding of cancer dissemination and the development of immunotherapeutics.

Molecular Profiling Reveals Biologically Discrete Subsets and Pathways of Progression in Diffuse Glioma.

Therapy development for adult diffuse glioma is hindered by incomplete knowledge of somatic glioma driving alterations and suboptimal disease classification. We defined the complete set of genes associated with 1,122 diffuse grade II-III-IV gliomas from The Cancer Genome Atlas and used molecular profiles to improve disease classification, identify molecular correlations, and provide insights into the progression from low- to high-grade disease. Whole-genome sequencing data analysis determined that ATRX but not TERT promoter mutations are associated with increased telomere length. Recent advances in glioma classification based on IDH mutation and 1p/19q co-deletion status were recapitulated through analysis of DNA methylation profiles, which identified clinically relevant molecular subsets. A subtype of IDH mutant glioma was associated with DNA demethylation and poor outcome; a group of IDH-wild-type diffuse glioma showed molecular similarity to pilocytic astrocytoma and relatively favorable survival. Understanding of cohesive disease groups may aid improved clinical outcomes.

Germline genetic contribution to the immune landscape of cancer.

Understanding the contribution of the host's genetic background to cancer immunity may lead to improved stratification for immunotherapy and to the identification of novel therapeutic targets. We investigated the effect of common and rare germline variants on 139 well-defined immune traits in ∼9000 cancer patients enrolled in TCGA. High heritability was observed for estimates of NK cell and T cell subset infiltration and for interferon signaling. Common variants of IFIH1, TMEM173 (STING1), and TMEM108 were associated with differential interferon signaling and variants mapping to RBL1 correlated with T cell subset abundance. Pathogenic or likely pathogenic variants in BRCA1 and in genes involved in telomere stabilization and Wnt-ß-catenin also acted as immune modulators. Our findings provide evidence for the impact of germline genetics on the composition and functional orientation of the tumor immune microenvironment. The curated datasets, variants, and genes identified provide a resource toward further understanding of tumor-immune interactions.

Reinventing gene expression connectivity through regulatory and spatial structural empowerment via principal node aggregation graph neural network.

The intricacies of the human genome, manifested as a complex network of genes, transcend conventional representations in text or numerical matrices. The intricate gene-to-gene relationships inherent in this complexity find a more suitable depiction in graph structures. In the pursuit of predicting gene expression, an endeavor shared by predecessors like the L1000 and Enformer methods, we introduce a novel spatial graph-neural network (GNN) approach. This innovative strategy incorporates graph features, encompassing both regulatory and structural elements. The regulatory elements include pair-wise gene correlation, biological pathways, protein-protein interaction networks, and transcription factor regulation. The spatial structural elements include chromosomal distance, histone modification and Hi-C inferred 3D genomic features. Principal Node Aggregation models, validated independently, emerge as frontrunners, demonstrating superior performance compared to traditional regression and other deep learning models. By embracing the spatial GNN paradigm, our method significantly advances the description of the intricate network of gene interactions, surpassing the performance, predictable scope, and initial requirements set by previous methods.

MOViDA: multiomics visible drug activity prediction with a biologically informed neural network model.

MOTIVATION: The process of drug development is inherently complex, marked by extended intervals from the inception of a pharmaceutical agent to its eventual launch in the market. Additionally, each phase in this process is associated with a significant failure rate, amplifying the inherent challenges of this task. Computational virtual screening powered by machine learning algorithms has emerged as a promising approach for predicting therapeutic efficacy. However, the complex relationships between the features learned by these algorithms can be challenging to decipher. RESULTS: We have engineered an artificial neural network model designed specifically for predicting drug sensitivity. This model utilizes a biologically informed visible neural network, thereby enhancing its interpretability. The trained model allows for an in-depth exploration of the biological pathways integral to prediction and the chemical attributes of drugs that impact sensitivity. Our model harnesses multiomics data derived from a different tumor tissue sources, as well as molecular descriptors that encapsulate the properties of drugs. We extended the model to predict drug synergy, resulting in favorable outcomes while retaining interpretability. Given the imbalanced nature of publicly available drug screening datasets, our model demonstrated superior performance to state-of-the-art visible machine learning algorithms. AVAILABILITY AND IMPLEMENTATION: MOViDA is implemented in Python using PyTorch library and freely available for download at https://github.com/Luigi-Ferraro/MOViDA. Training data, RIS score and drug features are archived on Zenodo https://doi.org/10.5281/zenodo.8180380.

Identification of therapeutic targets in osteoarthritis by combining heterogeneous transcriptional datasets, drug-induced expression profiles, and known drug-target interactions.

BACKGROUND: Osteoarthritis (OA) is a multifactorial, hypertrophic, and degenerative condition involving the whole joint and affecting a high percentage of middle-aged people. It is due to a combination of factors, although the pivotal mechanisms underlying the disease are still obscure. Moreover, current treatments are still poorly effective, and patients experience a painful and degenerative disease course. METHODS: We used an integrative approach that led us to extract a consensus signature from a meta-analysis of three different OA cohorts. We performed a network-based drug prioritization to detect the most relevant drugs targeting these genes and validated in vitro the most promising candidates. We also proposed a risk score based on a minimal set of genes to predict the OA clinical stage from RNA-Seq data. RESULTS: We derived a consensus signature of 44 genes that we validated on an independent dataset. Using network analysis, we identified Resveratrol, Tenoxicam, Benzbromarone, Pirinixic Acid, and Mesalazine as putative drugs of interest for therapeutics in OA for anti-inflammatory properties. We also derived a list of seven gene-targets validated with functional RT-qPCR assays, confirming the in silico predictions. Finally, we identified a predictive subset of genes composed of DNER, TNFSF11, THBS3, LOXL3, TSPAN2, DYSF, ASPN and HTRA1 to compute the patient's risk score. We validated this risk score on an independent dataset with a high AUC (0.875) and compared it with the same approach computed using the entire consensus signature (AUC 0.922). CONCLUSIONS: The consensus signature highlights crucial mechanisms for disease progression. Moreover, these genes were associated with several candidate drugs that could represent potential innovative therapeutics. Furthermore, the patient's risk scores can be used in clinical settings.

Mutational landscape of inflammatory breast cancer.

BACKGROUND: Inflammatory breast cancer (IBC) is the most pro-metastatic form of BC. Better understanding of its enigmatic pathophysiology is crucial. We report here the largest whole-exome sequencing (WES) study of clinical IBC samples. METHODS: We retrospectively applied WES to 54 untreated IBC primary tumor samples and matched normal DNA. The comparator samples were 102 stage-matched non-IBC samples from TCGA. We compared the somatic mutational profiles, spectra and signatures, copy number alterations (CNAs), HRD and heterogeneity scores, and frequencies of actionable genomic alterations (AGAs) between IBCs and non-IBCs. The comparisons were adjusted for the molecular subtypes. RESULTS: The number of somatic mutations, TMB, and mutational spectra were not different between IBCs and non-IBCs, and no gene was differentially mutated or showed differential frequency of CNAs. Among the COSMIC signatures, only the age-related signature was more frequent in non-IBCs than in IBCs. We also identified in IBCs two new mutational signatures not associated with any environmental exposure, one of them having been previously related to HIF pathway activation. Overall, the HRD score was not different between both groups, but was higher in TN IBCs than TN non-IBCs. IBCs were less frequently classified as heterogeneous according to heterogeneity H-index than non-IBCs (21% vs 33%), and clonal mutations were more frequent and subclonal mutations less frequent in IBCs. More than 50% of patients with IBC harbored at least one high-level of evidence (LOE) AGA (OncoKB LOE 1-2, ESCAT LOE I-II), similarly to patients with non-IBC. CONCLUSIONS: We provide the largest mutational landscape of IBC. Only a few subtle differences were identified with non-IBCs. The most clinically relevant one was the higher HRD score in TN IBCs than in TN non-IBCs, whereas the most intriguing one was the smaller intratumor heterogeneity of IBCs.

Whole-exome profiles of inflammatory breast cancer and pathological response to neoadjuvant chemotherapy.

BACKGROUND: Neoadjuvant chemotherapy (NACT) became a standard treatment strategy for patients with inflammatory breast cancer (IBC) because of high disease aggressiveness. However, given the heterogeneity of IBC, no molecular feature reliably predicts the response to chemotherapy. Whole-exome sequencing (WES) of clinical tumor samples provides an opportunity to identify genomic alterations associated with chemosensitivity. METHODS: We retrospectively applied WES to 44 untreated IBC primary tumor samples and matched normal DNA. The pathological response to NACT, assessed on operative specimen, distinguished the patients with versus without pathological complete response (pCR versus no-pCR respectively). We compared the mutational profiles, spectra and signatures, pathway mutations, copy number alterations (CNAs), HRD, and heterogeneity scores between pCR versus no-pCR patients. RESULTS: The TMB, HRD, and mutational spectra were not different between the complete (N = 13) versus non-complete (N = 31) responders. The two most frequently mutated genes were TP53 and PIK3CA. They were more frequently mutated in the complete responders, but the difference was not significant. Only two genes, NLRP3 and SLC9B1, were significantly more frequently mutated in the complete responders (23% vs. 0%). By contrast, several biological pathways involved in protein translation, PI3K pathway, and signal transduction showed significantly higher mutation frequency in the patients with pCR. We observed a higher abundance of COSMIC signature 7 (due to ultraviolet light exposure) in tumors from complete responders. The comparison of CNAs of the 3808 genes included in the GISTIC regions between both patients' groups identified 234 genes as differentially altered. The CIN signatures were not differentially represented between the complete versus non-complete responders. Based on the H-index, the patients with heterogeneous tumors displayed a lower pCR rate (11%) than those with less heterogeneous tumors (35%). CONCLUSIONS: This is the first study aiming at identifying correlations between the WES data of IBC samples and the achievement of pCR to NACT. Our results, obtained in this 44-sample series, suggest a few subtle genomic alterations associated with pathological response. Additional investigations are required in larger series.

Differential methylation of circulating free DNA assessed through cfMeDiP as a new tool for breast cancer diagnosis and detection of BRCA1/2 mutation.

BACKGROUND: Recent studies have highlighted the importance of the cell-free DNA (cfDNA) methylation profile in detecting breast cancer (BC) and its different subtypes. We investigated whether plasma cfDNA methylation, using cell-free Methylated DNA Immunoprecipitation and High-Throughput Sequencing (cfMeDIP-seq), may be informative in characterizing breast cancer in patients with BRCA1/2 germline mutations for early cancer detection and response to therapy. METHODS: We enrolled 23 BC patients with germline mutation of BRCA1 and BRCA2 genes, 19 healthy controls without BRCA1/2 mutation, and two healthy individuals who carried BRCA1/2 mutations. Blood samples were collected for all study subjects at the diagnosis, and plasma was isolated by centrifugation. Cell-free DNA was extracted from 1 mL of plasma, and cfMeDIP-seq was performed for each sample. Shallow whole genome sequencing was performed on the immuno-precipitated samples. Then, the differentially methylated 300-bp regions (DMRs) between 25 BRCA germline mutation carriers and 19 non-carriers were identified. DMRs were compared with tumor-specific regions from public datasets to perform an unbiased analysis. Finally, two statistical classifiers were trained based on the GLMnet and random forest model to evaluate if the identified DMRs could discriminate BRCA-positive from healthy samples. RESULTS: We identified 7,095 hypermethylated and 212 hypomethylated regions in 25 BRCA germline mutation carriers compared to 19 controls. These regions discriminate tumors from healthy samples with high accuracy and sensitivity. We show that the circulating tumor DNA of BRCA1/2 mutant breast cancers is characterized by the hypomethylation of genes involved in DNA repair and cell cycle. We uncovered the TFs associated with these DRMs and identified that proteins of the Erythroblast Transformation Specific (ETS) family are particularly active in the hypermethylated regions. Finally, we assessed that these regions could discriminate between BRCA positives from healthy samples with an AUC of 0.95, a sensitivity of 88%, and a specificity of 94.74%. CONCLUSIONS: Our study emphasizes the importance of tumor cell-derived DNA methylation in BC, reporting a different methylation profile between patients carrying mutations in BRCA1, BRCA2, and wild-type controls. Our minimally invasive approach could allow early cancer diagnosis, assessment of minimal residual disease, and monitoring of response to therapy.

Identification and bioinformatic characterization of a serum miRNA signature for early detection of laryngeal squamous cell carcinoma.

BACKGROUND: The growing understanding of cancer biology and the establishment of new treatment modalities has not yielded the expected results in terms of survival for Laryngeal Squamous Cell Cancer (LSCC). Early diagnosis, as well as prompt identification of patients with high risk of relapse would ensure greater chance of therapeutic success. However, this goal remains a challenge due to the absence of specific biomarkers for this neoplasm. METHODS: Serum samples from 45 LSCC patients and 23 healthy donors were collected for miRNA expression profiling by TaqMan Array analysis. Additional 20 patients and 42 healthy volunteers were included for the validation set, reaching an equal number of clinical samples for each group. The potential diagnostic ability of the such identified three-miRNA signature was confirmed by ROC analysis. Moreover, each miRNA was analyzed for the possible correlation with HNSCC patients' survival and TNM status by online databases Kaplan-Meier (KM) plotter and OncomiR. In silico analysis of common candidate targets and their network relevance to predict shared biological functions was finally performed by PANTHER and GeneMANIA software. RESULTS: We characterized serum miRNA profile of LSCC patients identifying a novel molecular signature, including miR-223, miR-93 and miR-532, as circulating marker endowed with high selectivity and specificity. The oncogenic effect and the prognostic significance of each miRNA was investigated by bioinformatic analysis, denoting significant correlation with OS. To analyse the molecular basis underlying the pro-tumorigenic role of the signature, we focused on the simultaneously regulated gene targets-IL6ST, GTDC1, MAP1B, CPEB3, PRKACB, NFIB, PURB, ATP2B1, ZNF148, PSD3, TBC1D15, PURA, KLF12-found by prediction tools and deepened for their functional role by pathway enrichment analysis. The results showed the involvement of 7 different biological processes, among which inflammation, proliferation, migration, apoptosis and angiogenesis. CONCLUSIONS: In conclusion, we have identified a possible miRNA signature for early LSCC diagnosis and we assumed that miR-93, miR-223 and miR-532 could orchestrate the regulation of multiple cancer-related processes. These findings encourage the possibility to deepen the molecular mechanisms underlying their oncogenic role, for the desirable development of novel therapeutic opportunities based on the use of short single-stranded oligonucleotides acting as non-coding RNA antagonists in cancer.

Epigenetic remodeling to improve the efficacy of immunotherapy in human glioblastoma: pre-clinical evidence for development of new immunotherapy approaches.

BACKGROUND: Glioblastoma multiforme (GBM) is a highly aggressive primary brain tumor, that is refractory to standard treatment and to immunotherapy with immune-checkpoint inhibitors (ICI). Noteworthy, melanoma brain metastases (MM-BM), that share the same niche as GBM, frequently respond to current ICI therapies. Epigenetic modifications regulate GBM cellular proliferation, invasion, and prognosis and may negatively regulate the cross-talk between malignant cells and immune cells in the tumor milieu, likely contributing to limit the efficacy of ICI therapy of GBM. Thus, manipulating the tumor epigenome can be considered a therapeutic opportunity in GBM. METHODS: Microarray transcriptional and methylation profiles, followed by gene set enrichment and IPA analyses, were performed to study the differences in the constitutive expression profiles of GBM vs MM-BM cells, compared to the extracranial MM cells and to investigate the modulatory effects of the DNA hypomethylating agent (DHA) guadecitabine among the different tumor cells. The prognostic relevance of DHA-modulated genes was tested by Cox analysis in a TCGA GBM patients' cohort. RESULTS: The most striking differences between GBM and MM-BM cells were found to be the enrichment of biological processes associated with tumor growth, invasion, and extravasation with the inhibition of MHC class II antigen processing/presentation in GBM cells. Treatment with guadecitabine reduced these biological differences, shaping GBM cells towards a more immunogenic phenotype. Indeed, in GBM cells, promoter hypomethylation by guadecitabine led to the up-regulation of genes mainly associated with activation, proliferation, and migration of T and B cells and with MHC class II antigen processing/presentation. Among DHA-modulated genes in GBM, 7.6% showed a significant prognostic relevance. Moreover, a large set of immune-related upstream-regulators (URs) were commonly modulated by DHA in GBM, MM-BM, and MM cells: DHA-activated URs enriched for biological processes mainly involved in the regulation of cytokines and chemokines production, inflammatory response, and in Type I/II/III IFN-mediated signaling; conversely, DHA-inhibited URs were involved in metabolic and proliferative pathways. CONCLUSIONS: Epigenetic remodeling by guadecitabine represents a promising strategy to increase the efficacy of cancer immunotherapy of GBM, supporting the rationale to develop new epigenetic-based immunotherapeutic approaches for the treatment of this still highly deadly disease.

A community challenge to predict clinical outcomes after immune checkpoint blockade in non-small cell lung cancer.

BACKGROUND: Predictive biomarkers of immune checkpoint inhibitor (ICI) efficacy are currently lacking for non-small cell lung cancer (NSCLC). Here, we describe the results from the Anti-PD-1 Response Prediction DREAM Challenge, a crowdsourced initiative that enabled the assessment of predictive models by using data from two randomized controlled clinical trials (RCTs) of ICIs in first-line metastatic NSCLC. METHODS: Participants developed and trained models using public resources. These were evaluated with data from the CheckMate 026 trial (NCT02041533), according to the model-to-data paradigm to maintain patient confidentiality. The generalizability of the models with the best predictive performance was assessed using data from the CheckMate 227 trial (NCT02477826). Both trials were phase III RCTs with a chemotherapy control arm, which supported the differentiation between predictive and prognostic models. Isolated model containers were evaluated using a bespoke strategy that considered the challenges of handling transcriptome data from clinical trials. RESULTS: A total of 59 teams participated, with 417 models submitted. Multiple predictive models, as opposed to a prognostic model, were generated for predicting overall survival, progression-free survival, and progressive disease status with ICIs. Variables within the models submitted by participants included tumor mutational burden (TMB), programmed death ligand 1 (PD-L1) expression, and gene-expression-based signatures. The best-performing models showed improved predictive power over reference variables, including TMB or PD-L1. CONCLUSIONS: This DREAM Challenge is the first successful attempt to use protected phase III clinical data for a crowdsourced effort towards generating predictive models for ICI clinical outcomes and could serve as a blueprint for similar efforts in other tumor types and disease states, setting a benchmark for future studies aiming to identify biomarkers predictive of ICI efficacy. TRIAL REGISTRATION: CheckMate 026; NCT02041533, registered January 22, 2014. CheckMate 227; NCT02477826, registered June 23, 2015.

A metabolic function of FGFR3-TACC3 gene fusions in cancer.

Chromosomal translocations that generate in-frame oncogenic gene fusions are notable examples of the success of targeted cancer therapies. We have previously described gene fusions of FGFR3-TACC3 (F3-T3) in 3% of human glioblastoma cases. Subsequent studies have reported similar frequencies of F3-T3 in many other cancers, indicating that F3-T3 is a commonly occuring fusion across all tumour types. F3-T3 fusions are potent oncogenes that confer sensitivity to FGFR inhibitors, but the downstream oncogenic signalling pathways remain unknown. Here we show that human tumours with F3-T3 fusions cluster within transcriptional subgroups that are characterized by the activation of mitochondrial functions. F3-T3 activates oxidative phosphorylation and mitochondrial biogenesis and induces sensitivity to inhibitors of oxidative metabolism. Phosphorylation of the phosphopeptide PIN4 is an intermediate step in the signalling pathway of the activation of mitochondrial metabolism. The F3-T3-PIN4 axis triggers the biogenesis of peroxisomes and the synthesis of new proteins. The anabolic response converges on the PGC1α coactivator through the production of intracellular reactive oxygen species, which enables mitochondrial respiration and tumour growth. These data illustrate the oncogenic circuit engaged by F3-T3 and show that F3-T3-positive tumours rely on mitochondrial respiration, highlighting this pathway as a therapeutic opportunity for the treatment of tumours with F3-T3 fusions. We also provide insights into the genetic alterations that initiate the chain of metabolic responses that drive mitochondrial metabolism in cancer.

A review of COVID-19 biomarkers and drug targets: resources and tools.

The stratification of patients at risk of progression of COVID-19 and their molecular characterization is of extreme importance to optimize treatment and to identify therapeutic options. The bioinformatics community has responded to the outbreak emergency with a set of tools and resource to identify biomarkers and drug targets that we review here. Starting from a consolidated corpus of 27 570 papers, we adopt latent Dirichlet analysis to extract relevant topics and select those associated with computational methods for biomarker identification and drug repurposing. The selected topics span from machine learning and artificial intelligence for disease characterization to vaccine development and to therapeutic target identification. Although the way to go for the ultimate defeat of the pandemic is still long, the amount of knowledge, data and tools generated so far constitutes an unprecedented example of global cooperation to this threat.

Network-based identification of key master regulators associated with an immune-silent cancer phenotype.

A cancer immune phenotype characterized by an active T-helper 1 (Th1)/cytotoxic response is associated with responsiveness to immunotherapy and favorable prognosis across different tumors. However, in some cancers, such an intratumoral immune activation does not confer protection from progression or relapse. Defining mechanisms associated with immune evasion is imperative to refine stratification algorithms, to guide treatment decisions and to identify candidates for immune-targeted therapy. Molecular alterations governing mechanisms for immune exclusion are still largely unknown. The availability of large genomic datasets offers an opportunity to ascertain key determinants of differential intratumoral immune response. We follow a network-based protocol to identify transcription regulators (TRs) associated with poor immunologic antitumor activity. We use a consensus of four different pipelines consisting of two state-of-the-art gene regulatory network inference techniques, regularized gradient boosting machines and ARACNE to determine TR regulons, and three separate enrichment techniques, including fast gene set enrichment analysis, gene set variation analysis and virtual inference of protein activity by enriched regulon analysis to identify the most important TRs affecting immunologic antitumor activity. These TRs, referred to as master regulators (MRs), are unique to immune-silent and immune-active tumors, respectively. We validated the MRs coherently associated with the immune-silent phenotype across cancers in The Cancer Genome Atlas and a series of additional datasets in the Prediction of Clinical Outcomes from Genomic Profiles repository. A downstream analysis of MRs specific to the immune-silent phenotype resulted in the identification of several enriched candidate pathways, including NOTCH1, TGF-$\beta $, Interleukin-1 and TNF-$\alpha $ signaling pathways. TGFB1I1 emerged as one of the main negative immune modulators preventing the favorable effects of a Th1/cytotoxic response.

Pan-cancer onco-signatures reveal a novel mitochondrial subtype of luminal breast cancer with specific regulators.

BACKGROUND: Somatic alterations in cancer cause dysregulation of signaling pathways that control cell-cycle progression, apoptosis, and cell growth. The effect of individual alterations in these pathways differs between individual tumors and tumor types. Recognizing driver events is a complex task requiring integrating multiple molecular data, including genomics, epigenomics, and functional genomics. A common hypothesis is that these driver events share similar effects on the hallmarks of cancer. The availability of large-scale multi-omics studies allows for inferring these common effects from data. Once these effects are known, one can then deconvolve in every individual patient whether a given genomics alteration is a driver event. METHODS: Here, we develop a novel data-driven approach to identify shared oncogenic expression signatures among tumors. We aim to identify gene onco-signature for classifying tumor patients in homogeneous subclasses with distinct prognoses and specific genomic alterations. We derive expression pan-cancer onco-signatures from TCGA gene expression data using a discovery set of 9107 primary pan-tumor samples together with respective matched mutational data and a list of known cancer-related genes from COSMIC database. RESULTS: We use the derived ono-signatures to state their prognostic significance and apply them to the TCGA breast cancer dataset as proof of principle of our approach. We uncover a "mitochondrial" sub-group of Luminal patients characterized by its biological features and regulated by specific genetic modulators. Collectively, our results demonstrate the effectiveness of onco-signatures-based methodologies, and they also contribute to a comprehensive understanding of the metabolic heterogeneity of Luminal tumors. CONCLUSIONS: These findings provide novel genomics evidence for developing personalized breast cancer patient treatments. The onco-signature approach, demonstrated here on breast cancer, is general and can be applied to other cancer types.

Pancancer network analysis reveals key master regulators for cancer invasiveness.

BACKGROUND: Tumor invasiveness reflects numerous biological changes, including tumorigenesis, progression, and metastasis. To decipher the role of transcriptional regulators (TR) involved in tumor invasiveness, we performed a systematic network-based pan-cancer assessment of master regulators of cancer invasiveness. MATERIALS AND METHODS: We stratified patients in The Cancer Genome Atlas (TCGA) into invasiveness high (INV-H) and low (INV-L) groups using consensus clustering based on an established robust 24-gene signature to determine the prognostic association of invasiveness with overall survival (OS) across 32 different cancers. We devise a network-based protocol to identify TRs as master regulators (MRs) unique to INV-H and INV-L phenotypes. We validated the activity of MRs coherently associated with INV-H phenotype and worse OS across cancers in TCGA on a series of additional datasets in the Prediction of Clinical Outcomes from the Genomic Profiles (PRECOG) repository. RESULTS: Based on the 24-gene signature, we defined the invasiveness score for each patient sample and stratified patients into INV-H and INV-L clusters. We observed that invasiveness was associated with worse survival outcomes in almost all cancers and had a significant association with OS in ten out of 32 cancers. Our network-based framework identified common invasiveness-associated MRs specific to INV-H and INV-L groups across the ten prognostic cancers, including COL1A1, which is also part of the 24-gene signature, thus acting as a positive control. Downstream pathway analysis of MRs specific to INV-H phenotype resulted in the identification of several enriched pathways, including Epithelial into Mesenchymal Transition, TGF-ß signaling pathway, regulation of Toll-like receptors, cytokines, and inflammatory response, and selective expression of chemokine receptors during T-cell polarization. Most of these pathways have connotations of inflammatory immune response and feasibility for metastasis. CONCLUSION: Our pan-cancer study provides a comprehensive master regulator analysis of tumor invasiveness and can suggest more precise therapeutic strategies by targeting the identified MRs and downstream enriched pathways for patients across multiple cancers.

The epitranscriptome of high-grade gliomas: a promising therapeutic target with implications from the tumor microenvironment to endogenous retroviruses.

Glioblastoma (GBM) comprises 45.6% of all primary malignant brain cancers and is one of the most common and aggressive intracranial tumors in adults. Intratumoral heterogeneity with a wide range of proteomic, genetic, and epigenetic dysregulation contributes to treatment resistance and poor prognosis, thus demanding novel therapeutic approaches. To date, numerous clinical trials have been developed to target the proteome and epigenome of high-grade gliomas with promising results. However, studying RNA modifications, or RNA epitranscriptomics, is a new frontier within neuro-oncology. RNA epitranscriptomics was discovered in the 1970s, but in the last decade, the extent of modification of mRNA and various non-coding RNAs has emerged and been implicated in transposable element activation and many other oncogenic processes within the tumor microenvironment. This review provides background information and discusses the therapeutic potential of agents modulating epitranscriptomics in high-grade gliomas. A particular emphasis will be placed on how combination therapies that include immune agents targeting hERV-mediated viral mimicry could improve the treatment of GBM.

Immunologic constant of rejection as a predictive biomarker of immune checkpoint inhibitors efficacy in non-small cell lung cancer.

BACKGROUND: Anti-PD1/PDL1 immune checkpoint inhibitors (ICI) transformed the prognosis of patients with advanced non-small cell lung cancer (NSCLC). However, the response rate remains disappointing and toxicity may be life-threatening, making urgent identification of biomarkers predictive for efficacy. Immunologic Constant of Rejection signature (ICR) is a 20-gene expression signature of cytotoxic immune response with prognostic value in some solid cancers. Our objective was to assess its predictive value for benefit from anti-PD1/PDL1 in patients with advanced NSCLC. METHODS: We retrospectively profiled 44 primary tumors derived from NSCLC patients treated with ICI as single-agent in at least the second-line metastatic setting. Transcriptomic analysis was performed using the nCounter® analysis system and the PanCancer Immune Profiling Panel. We then pooled our data with clinico-biological data from four public gene expression data sets, leading to a total of 162 NSCLC patients treated with single-agent anti-PD1/PDL1. ICR was applied to all samples and correlation was searched between ICR classes and the Durable Clinical Benefit (DCB), defined as stable disease or objective response according to RECIST 1.1 for a minimum of 6 months after the start of ICI. RESULTS: The DCB rate was 29%; 22% of samples were classified as ICR1, 30% ICR2, 22% ICR3, and 26% ICR4. These classes were not associated with the clinico-pathological variables, but showed enrichment from ICR1 to ICR4 in quantitative/qualitative markers of immune response. ICR2-4 class was associated with a 5.65-fold DCB rate when compared with ICR1 class. In multivariate analysis, ICR classification remained associated with DCB, independently from PDL1 expression and other predictive immune signatures. By contrast, it was not associated with disease-free survival in 556 NSCLC TCGA patients untreated with ICI. CONCLUSION: The 20-gene ICR signature was independently associated with benefit from anti-PD1/PDL1 ICI in patients with advanced NSCLC. Validation in larger retrospective and prospective series is warranted.

Genetic predisposition to cancer across people of different ancestries in Qatar: a population-based, cohort study.

BACKGROUND: Disparities in the genetic risk of cancer among various ancestry groups and populations remain poorly defined. This challenge is even more acute for Middle Eastern populations, where the paucity of genomic data could affect the clinical potential of cancer genetic risk profiling. We used data from the phase 1 cohort of the Qatar Genome Programme to investigate genetic variation in cancer-susceptibility genes in the Qatari population. METHODS: The Qatar Genome Programme generated high-coverage genome sequencing on DNA samples collected from 6142 native Qataris, stratified into six distinct ancestry groups: general Arab, Persian, Arabian Peninsula, Admixture Arab, African, and South Asian. In this population-based, cohort study, we evaluated the performance of polygenic risk scores for the most common cancers in Qatar (breast, prostate, and colorectal cancers). Polygenic risk scores were trained in The Cancer Genome Atlas (TCGA) dataset, and their distributions were subsequently applied to the six different genetic ancestry groups of the Qatari population. Rare deleterious variants within 1218 cancer susceptibility genes were analysed, and their clinical pathogenicity was assessed by ClinVar and the CharGer computational tools. FINDINGS: The cohort included in this study was recruited by the Qatar Biobank between Dec 11, 2012, and June 9, 2016. The initial dataset comprised 6218 cohort participants, and whole genome sequencing quality control filtering led to a final dataset of 6142 samples. Polygenic risk score analyses of the most common cancers in Qatar showed significant differences between the six ancestry groups (p<0·0001). Qataris with Arabian Peninsula ancestry showed the lowest polygenic risk score mean for colorectal cancer (-0·41), and those of African ancestry showed the highest average for prostate cancer (0·85). Cancer-gene rare variant analysis identified 76 Qataris (1·2% of 6142 individuals in the Qatar Genome Programme cohort) carrying ClinVar pathogenic or likely pathogenic variants in clinically actionable cancer genes. Variant analysis using CharGer identified 195 individuals carriers (3·17% of the cohort). Breast cancer pathogenic variants were over-represented in Qataris of Persian origin (22 [56·4%] of 39 BRCA1/BRCA2 variant carriers) and completely absent in those of Arabian Peninsula origin. INTERPRETATION: We observed a high degree of heterogeneity for cancer predisposition genes and polygenic risk scores across ancestries in this population from Qatar. Stratification systems could be considered for the implementation of national cancer preventive medicine programmes. FUNDING: Qatar Foundation.