Mapping functional to morphological variation reveals the basis of regional extracellular matrix subversion and nerve invasion in pancreatic cancer.

Intratumor morphological heterogeneity of pancreatic ductal adenocarcinoma (PDAC) predicts clinical outcomes but is only partially understood at the molecular level. To elucidate the gene expression programs underpinning intratumor morphological variation in PDAC, we investigated and deconvoluted at single cell level the molecular profiles of histologically distinct clusters of PDAC cells. We identified three major morphological and functional variants that co-exist in varying proportions in all PDACs, display limited genetic diversity, and are associated with a distinct organization of the extracellular matrix: a glandular variant with classical ductal features; a transitional variant displaying abortive ductal structures and mixed endodermal and myofibroblast-like gene expression; and a poorly differentiated variant lacking ductal features and basement membrane, and showing neuronal lineage priming. Ex vivo and in vitro evidence supports the occurrence of dynamic transitions among these variants in part influenced by extracellular matrix composition and stiffness and associated with local, specifically neural, invasion.

Tumoral and stromal hMENA isoforms impact tertiary lymphoid structure localization in lung cancer and predict immune checkpoint blockade response in patients with cancer.

BACKGROUND: Tertiary Lymphoid Structures (TLS) correlate with positive outcomes in patients with NSCLC and the efficacy of immune checkpoint blockade (ICB) in cancer. The actin regulatory protein hMENA undergoes tissue-specific splicing, producing the epithelial hMENA11a linked to favorable prognosis in early NSCLC, and the mesenchymal hMENAΔv6 found in invasive cancer cells and pro-tumoral cancer-associated fibroblasts (CAFs). This study investigates how hMENA isoforms in tumor cells and CAFs relate to TLS presence, localization and impact on patient outcomes and ICB response. METHODS: Methods involved RNA-SEQ on NSCLC cells with depleted hMENA isoforms. A retrospective observational study assessed tissues from surgically treated N0 patients with NSCLC, using immunohistochemistry for tumoral and stromal hMENA isoforms, fibronectin, and TLS presence. ICB-treated patient tumors were analyzed using Nanostring nCounter and GeoMx spatial transcriptomics. Multiparametric flow cytometry characterized B cells and tissue-resident memory T cells (TRM). Survival and ICB response were estimated in the cohort and validated using bioinformatics pipelines in different datasets. FINDINGS: Findings indicate that hMENA11a in NSCLC cells upregulates the TLS regulator LTßR, decreases fibronectin, and favors CXCL13 production by TRM. Conversely, hMENAΔv6 in CAFs inhibits LTßR-related NF-kB pathway, reduces CXCL13 secretion, and promotes fibronectin production. These patterns are validated in N0 NSCLC tumors, where hMENA11ahigh expression, CAF hMENAΔv6low, and stromal fibronectinlow are associated with intratumoral TLS, linked to memory B cells and predictive of longer survival. The hMENA isoform pattern, fibronectin, and LTßR expression broadly predict ICB response in tumors where TLS indicates an anti-tumor immune response. INTERPRETATION: This study uncovers hMENA alternative splicing as an unexplored contributor to TLS-related Tumor Immune Microenvironment (TIME) and a promising biomarker for clinical outcomes and likely ICB responsiveness in N0 patients with NSCLC. FUNDING: This work is supported by AIRC (IG 19822), ACC (RCR-2019-23669120), CAL.HUB.RIA Ministero Salute PNRR-POS T4, "Ricerca Corrente" granted by the Italian Ministry of Health.

Benchmarking and improving the performance of variant-calling pipelines with RecallME.

MOTIVATION: The steady increment of Whole Genome/Exome sequencing and the development of novel Next Generation Sequencing-based gene panels requires continuous testing and validation of variant calling (VC) pipelines and the detection of sequencing-related issues to be maintained up-to-date and feasible for the clinical settings. State of the art tools are reliable when used to compute standard performance metrics. However, the need for an automated software to discriminate between bioinformatic and sequencing issues and to optimize VC parameters remains unmet. RESULTS: The aim of the current work is to present RecallME, a bioinformatic suite that tracks down difficult-to-detect variants as insertions and deletions in highly repetitive regions, thus providing the maximum reachable recall for both single nucleotide variants and small insertion and deletions and to precisely guide the user in the pipeline optimization process. AVAILABILITY AND IMPLEMENTATION: Source code is freely available under MIT license at https://github.com/mazzalab-ieo/recallme. RecallME web application is available at https://translational-oncology-lab.shinyapps.io/recallme/. To use RecallME, users must obtain a license for ANNOVAR by themselves.

GASOLINE: detecting germline and somatic structural variants from long-reads data.

Long-read sequencing allows analyses of single nucleic-acid molecules and produces sequences in the order of tens to hundreds kilobases. Its application to whole-genome analyses allows identification of complex genomic structural-variants (SVs) with unprecedented resolution. SV identification, however, requires complex computational methods, based on either read-depth or intra- and inter-alignment signatures approaches, which are limited by size or type of SVs. Moreover, most currently available tools only detect germline variants, thus requiring separate computation of sample pairs for comparative analyses. To overcome these limits, we developed a novel tool (Germline And SOmatic structuraL varIants detectioN and gEnotyping; GASOLINE) that groups SV signatures using a sophisticated clustering procedure based on a modified reciprocal overlap criterion, and is designed to identify germline SVs, from single samples, and somatic SVs from paired test and control samples. GASOLINE is a collection of Perl, R and Fortran codes, it analyzes aligned data in BAM format and produces VCF files with statistically significant somatic SVs. Germline or somatic analysis of 30[Formula: see text] sequencing coverage experiments requires 4-5 h with 20 threads. GASOLINE outperformed currently available methods in the detection of both germline and somatic SVs in synthetic and real long-reads datasets. Notably, when applied on a pair of metastatic melanoma and matched-normal sample, GASOLINE identified five genuine somatic SVs that were missed using five different sequencing technologies and state-of-the art SV calling approaches. Thus, GASOLINE identifies germline and somatic SVs with unprecedented accuracy and resolution, outperforming currently available state-of-the-art WGS long-reads computational methods.

High-resolution Nanopore methylome-maps reveal random hyper-methylation at CpG-poor regions as driver of chemoresistance in leukemias.

Aberrant DNA methylation at CpG dinucleotides is a cancer hallmark that is associated with the emergence of resistance to anti cancer treatment, though molecular mechanisms and biological significance remain elusive. Genome scale methylation maps by currently used methods are based on chemical modification of DNA and are best suited for analyses of methylation at CpG rich regions (CpG islands). We report the first high coverage whole-genome map in cancer using the long read nanopore technology, which allows simultaneous DNA-sequence and -methylation analyses on native DNA. We analyzed clonal epigenomic/genomic evolution in Acute Myeloid Leukemias (AMLs) at diagnosis and relapse, after chemotherapy. Long read sequencing coupled to a novel computational method allowed definition of differential methylation at unprecedented resolution, and showed that the relapse methylome is characterized by hypermethylation at both CpG islands and sparse CpGs regions. Most differentially methylated genes, however, were not differentially expressed nor enriched for chemoresistance genes. A small fraction of under-expressed and hyper-methylated genes at sparse CpGs, in the gene body, was significantly enriched in transcription factors (TFs). Remarkably, these few TFs supported large gene-regulatory networks including 50% of all differentially expressed genes in the relapsed AMLs and highly-enriched in chemoresistance genes. Notably, hypermethylated regions at sparse CpGs were poorly conserved in the relapsed AMLs, under-represented at their genomic positions and showed higher methylation entropy, as compared to CpG islands. Analyses of available datasets confirmed TF binding to their target genes and conservation of the same gene-regulatory networks in large patient cohorts. Relapsed AMLs carried few patient specific structural variants and DNA mutations, apparently not involved in drug resistance. Thus, drug resistance in AMLs can be mainly ascribed to the selection of random epigenetic alterations at sparse CpGs of a few transcription factors, which then induce reprogramming of the relapsing phenotype, independently of clonal genomic evolution.

In vivo temporal resolution of acute promyelocytic leukemia progression reveals a role of Klf4 in suppressing early leukemic transformation.

Genome organization plays a pivotal role in transcription, but how transcription factors (TFs) rewire the structure of the genome to initiate and maintain the programs that lead to oncogenic transformation remains poorly understood. Acute promyelocytic leukemia (APL) is a fatal subtype of leukemia driven by a chromosomal translocation between the promyelocytic leukemia (PML) and retinoic acid receptor α (RARα) genes. We used primary hematopoietic stem and progenitor cells (HSPCs) and leukemic blasts that express the fusion protein PML-RARα as a paradigm to temporally dissect the dynamic changes in the epigenome, transcriptome, and genome architecture induced during oncogenic transformation. We found that PML-RARα initiates a continuum of topologic alterations, including switches from A to B compartments, transcriptional repression, loss of active histone marks, and gain of repressive histone marks. Our multiomics-integrated analysis identifies Klf4 as an early down-regulated gene in PML-RARα-driven leukemogenesis. Furthermore, we characterized the dynamic alterations in the Klf4 cis-regulatory network during APL progression and demonstrated that ectopic Klf4 overexpression can suppress self-renewal and reverse the differentiation block induced by PML-RARα. Our study provides a comprehensive in vivo temporal dissection of the epigenomic and topological reprogramming induced by an oncogenic TF and illustrates how topological architecture can be used to identify new drivers of malignant transformation.

Nanopore ReCappable sequencing maps SARS-CoV-2 5' capping sites and provides new insights into the structure of sgRNAs.

The SARS-CoV-2 virus has a complex transcriptome characterised by multiple, nested subgenomic RNAsused to express structural and accessory proteins. Long-read sequencing technologies such as nanopore direct RNA sequencing can recover full-length transcripts, greatly simplifying the assembly of structurally complex RNAs. However, these techniques do not detect the 5' cap, thus preventing reliable identification and quantification of full-length, coding transcript models. Here we used Nanopore ReCappable Sequencing (NRCeq), a new technique that can identify capped full-length RNAs, to assemble a complete annotation of SARS-CoV-2 sgRNAs and annotate the location of capping sites across the viral genome. We obtained robust estimates of sgRNA expression across cell lines and viral isolates and identified novel canonical and non-canonical sgRNAs, including one that uses a previously un-annotated leader-to-body junction site. The data generated in this work constitute a useful resource for the scientific community and provide important insights into the mechanisms that regulate the transcription of SARS-CoV-2 sgRNAs.

Lower probability and shorter duration of infections after COVID-19 vaccine correlate with anti-SARS-CoV-2 circulating IgGs.

The correlation between immune responses and protection from SARS-CoV-2 infections and its duration remains unclear. We performed a sanitary surveillance at the European Institute of Oncology (IEO) in Milan over a 17 months period. Pre-vaccination, in 1,493 participants, we scored 266 infections (17.8%) and 8 possible reinfections (3%). Post-vaccination, we identified 30 infections in 2,029 vaccinated individuals (1.5%). We report that the probability of infection post-vaccination is i) significantly lower compared to natural infection, ii) associated with a significantly shorter median duration of infection than that of first infection and reinfection, iii) anticorrelated with circulating antibody levels.

Fasting-Mimicking Diet Is Safe and Reshapes Metabolism and Antitumor Immunity in Patients with Cancer.

In tumor-bearing mice, cyclic fasting or fasting-mimicking diets (FMD) enhance the activity of antineoplastic treatments by modulating systemic metabolism and boosting antitumor immunity. Here we conducted a clinical trial to investigate the safety and biological effects of cyclic, five-day FMD in combination with standard antitumor therapies. In 101 patients, the FMD was safe, feasible, and resulted in a consistent decrease of blood glucose and growth factor concentration, thus recapitulating metabolic changes that mediate fasting/FMD anticancer effects in preclinical experiments. Integrated transcriptomic and deep-phenotyping analyses revealed that FMD profoundly reshapes anticancer immunity by inducing the contraction of peripheral blood immunosuppressive myeloid and regulatory T-cell compartments, paralleled by enhanced intratumor Th1/cytotoxic responses and an enrichment of IFNγ and other immune signatures associated with better clinical outcomes in patients with cancer. Our findings lay the foundations for phase II/III clinical trials aimed at investigating FMD antitumor efficacy in combination with standard antineoplastic treatments. SIGNIFICANCE: Cyclic FMD is well tolerated and causes remarkable systemic metabolic changes in patients with different tumor types and treated with concomitant antitumor therapies. In addition, the FMD reshapes systemic and intratumor immunity, finally activating several antitumor immune programs. Phase II/III clinical trials are needed to investigate FMD antitumor activity/efficacy.This article is highlighted in the In This Issue feature, p. 1.

TMBleR: a bioinformatic tool to optimize TMB estimation and predictive power.

MOTIVATION: Tumor mutational burden (TMB) has been proposed as a predictive biomarker for immunotherapy response in cancer patients, as it is thought to enrich for tumors with high neoantigen load. TMB assessed by whole-exome sequencing is considered the gold standard but remains confined to research settings. In the clinical setting, targeted gene panels sampling various genomic sizes along with diverse strategies to estimate TMB were proposed and no real standard has emerged yet. RESULTS: We provide the community with TMBleR, a tool to measure the clinical impact of various strategies of panel-based TMB measurement. AVAILABILITY AND IMPLEMENTATION: R package and docker container (GPL-3 Open Source license): https://acc-bioinfo.github.io/TMBleR/. Graphical-user interface website: https://bioserver.ieo.it/shiny/app/tmbler. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

VID22 counteracts G-quadruplex-induced genome instability.

Genome instability is a condition characterized by the accumulation of genetic alterations and is a hallmark of cancer cells. To uncover new genes and cellular pathways affecting endogenous DNA damage and genome integrity, we exploited a Synthetic Genetic Array (SGA)-based screen in yeast. Among the positive genes, we identified VID22, reported to be involved in DNA double-strand break repair. vid22Δ cells exhibit increased levels of endogenous DNA damage, chronic DNA damage response activation and accumulate DNA aberrations in sequences displaying high probabilities of forming G-quadruplexes (G4-DNA). If not resolved, these DNA secondary structures can block the progression of both DNA and RNA polymerases and correlate with chromosome fragile sites. Vid22 binds to and protects DNA at G4-containing regions both in vitro and in vivo. Loss of VID22 causes an increase in gross chromosomal rearrangement (GCR) events dependent on G-quadruplex forming sequences. Moreover, the absence of Vid22 causes defects in the correct maintenance of G4-DNA rich elements, such as telomeres and mtDNA, and hypersensitivity to the G4-stabilizing ligand TMPyP4. We thus propose that Vid22 is directly involved in genome integrity maintenance as a novel regulator of G4 metabolism.

Prospective Validation of the Italian Alliance Against Cancer Lung Panel in Patients With Advanced Non-Small-Cell Lung Cancer.

BACKGROUND: The deeper knowledge of non-small-cell lung cancer (NSCLC) biology and the discovery of driver molecular alterations have opened the era of precision medicine in lung oncology, thus significantly revolutionizing the diagnostic and therapeutic approach to NSCLC. In Italy, however, molecular assessment remains heterogeneous across the country, and numbers of patients accessing personalized treatments remain relatively low. Nationwide programs have demonstrated that the creation of consortia represent a successful strategy to increase the number of patients with a molecular classification. PATIENTS AND METHODS: The Alliance Against Cancer (ACC), a network of 25 Italian Research Institutes, has developed a targeted sequencing panel for the detection of genomic alterations in 182 genes in patients with a diagnosis of NSCLC (ACC lung panel). One thousand metastatic NSCLC patients will be enrolled onto a prospective trial designed to measure the sensitivity and specificity of the ACC lung panel as a tool for molecular screening compared to standard methods. RESULTS AND CONCLUSION: The ongoing trial is part of a nationwide strategy of ACC to develop infrastructures and improve competences to make the Italian research institutes independent for genomic profiling of cancer patients.

Author Correction: Differential epigenetic reprogramming in response to specific endocrine therapies promotes cholesterol biosynthesis and cellular invasion.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Acquired CYP19A1 amplification is an early specific mechanism of aromatase inhibitor resistance in ERα metastatic breast cancer.

Tumor evolution is shaped by many variables, potentially involving external selective pressures induced by therapies. After surgery, patients with estrogen receptor (ERα)-positive breast cancer are treated with adjuvant endocrine therapy, including selective estrogen receptor modulators (SERMs) and/or aromatase inhibitors (AIs). However, more than 20% of patients relapse within 10 years and eventually progress to incurable metastatic disease. Here we demonstrate that the choice of therapy has a fundamental influence on the genetic landscape of relapsed diseases. We found that 21.5% of AI-treated, relapsed patients had acquired CYP19A1 (encoding aromatase) amplification (CYP19A1amp). Relapsed patients also developed numerous mutations targeting key breast cancer-associated genes, including ESR1 and CYP19A1. Notably, CYP19A1amp cells also emerged in vitro, but only in AI-resistant models. CYP19A1 amplification caused increased aromatase activity and estrogen-independent ERα binding to target genes, resulting in CYP19A1amp cells showing decreased sensitivity to AI treatment. These data suggest that AI treatment itself selects for acquired CYP19A1amp and promotes local autocrine estrogen signaling in AI-resistant metastatic patients.

Differential epigenetic reprogramming in response to specific endocrine therapies promotes cholesterol biosynthesis and cellular invasion.

Endocrine therapies target the activation of the oestrogen receptor alpha (ERα) via distinct mechanisms, but it is not clear whether breast cancer cells can adapt to treatment using drug-specific mechanisms. Here we demonstrate that resistance emerges via drug-specific epigenetic reprogramming. Resistant cells display a spectrum of phenotypical changes with invasive phenotypes evolving in lines resistant to the aromatase inhibitor (AI). Orthogonal genomics analysis of reprogrammed regulatory regions identifies individual drug-induced epigenetic states involving large topologically associating domains (TADs) and the activation of super-enhancers. AI-resistant cells activate endogenous cholesterol biosynthesis (CB) through stable epigenetic activation in vitro and in vivo. Mechanistically, CB sparks the constitutive activation of oestrogen receptors alpha (ERα) in AI-resistant cells, partly via the biosynthesis of 27-hydroxycholesterol. By targeting CB using statins, ERα binding is reduced and cell invasion is prevented. Epigenomic-led stratification can predict resistance to AI in a subset of ERα-positive patients.

HDAC-regulated myomiRs control BAF60 variant exchange and direct the functional phenotype of fibro-adipogenic progenitors in dystrophic muscles.

Fibro-adipogenic progenitors (FAPs) are important components of the skeletal muscle regenerative environment. Whether FAPs support muscle regeneration or promote fibro-adipogenic degeneration is emerging as a key determinant in the pathogenesis of muscular diseases, including Duchenne muscular dystrophy (DMD). However, the molecular mechanism that controls FAP lineage commitment and activity is currently unknown. We show here that an HDAC-myomiR-BAF60 variant network regulates the fate of FAPs in dystrophic muscles of mdx mice. Combinatorial analysis of gene expression microarray, genome-wide chromatin remodeling by nuclease accessibility (NA) combined with next-generation sequencing (NA-seq), small RNA sequencing (RNA-seq), and microRNA (miR) high-throughput screening (HTS) against SWI/SNF BAF60 variants revealed that HDAC inhibitors (HDACis) derepress a "latent" myogenic program in FAPs from dystrophic muscles at early stages of disease. Specifically, HDAC inhibition induces two core components of the myogenic transcriptional machinery, MYOD and BAF60C, and up-regulates the myogenic miRs (myomiRs) (miR-1.2, miR-133, and miR-206), which target the alternative BAF60 variants BAF60A and BAF60B, ultimately directing promyogenic differentiation while suppressing the fibro-adipogenic phenotype. In contrast, FAPs from late stage dystrophic muscles are resistant to HDACi-induced chromatin remodeling at myogenic loci and fail to activate the promyogenic phenotype. These results reveal a previously unappreciated disease stage-specific bipotency of mesenchimal cells within the regenerative environment of dystrophic muscles. Resolution of such bipotency by epigenetic intervention with HDACis provides a molecular rationale for the in situ reprogramming of target cells to promote therapeutic regeneration of dystrophic muscles.

The DNA demethylating agent decitabine activates the TRAIL pathway and induces apoptosis in acute myeloid leukemia.

Although epigenetic drugs have been approved for use in selected malignancies, there is significant need for a better understanding of their mechanism of action. Here, we study the action of a clinically approved DNA-methyltransferase inhibitor - decitabine (DAC) - in acute myeloid leukemia (AML) cells. At low doses, DAC treatment induced apoptosis of NB4 Acute Promyelocytic Leukemia (APL) cells, which was associated with the activation of the extrinsic apoptotic pathway. Expression studies of the members of the Death Receptor family demonstrated that DAC induces the expression of TNF-related apoptosis-inducing ligand (TRAIL). Upregulation of TRAIL, upon DAC treatment, was associated with specific epigenetic modifications induced by DAC in the proximity of the TRAIL promoter, as demonstrated by DNA demethylation, increased DNaseI sensitivity and histone acetylation of a non-CpG island, CpG-rich region located 2kb upstream to the transcription start site. Luciferase assay experiments showed that this region behave as a DNA methylation sensitive transcriptional regulatory element. The CpG regulatory element was also found methylated in samples derived from APL patients. These findings have been confirmed in the non-APL, AML Kasumi cell line, suggesting that this regulatory mechanism may be extended to other AMLs. Our study suggests that DNA methylation is a regulatory mechanism relevant for silencing of the TRAIL apoptotic pathway in leukemic cells, and further elucidates the mechanism by which epigenetic drugs mediate their anti-leukemic effects.

Chromatin accessibility, p300, and histone acetylation define PML-RARα and AML1-ETO binding sites in acute myeloid leukemia.

Chromatin accessibility plays a key role in regulating cell type specific gene expression during hematopoiesis but has also been suggested to be aberrantly regulated during leukemogenesis. To understand the leukemogenic chromatin signature, we analyzed acute promyelocytic leukemia, a subtype of leukemia characterized by the expression of RARα-fusion proteins, such as PML-RARα. We used nuclease accessibility sequencing in cell lines as well as patient blasts to identify accessible DNA elements and identified > 100 000 accessible regions in each case. Using ChIP-seq, we identified H2A.Z as a histone modification generally associated with these accessible regions, whereas unsupervised clustering analysis of other chromatin features, including DNA methylation, H2A.Zac, H3ac, H3K9me3, H3K27me3, and the regulatory factor p300, distinguished 6 distinct clusters of accessible sites, each with a characteristic functional makeup. Of these, PML-RARα binding was found specifically at accessible chromatin regions characterized by p300 binding and hypoacetylated histones. Identifying regions with a similar epigenetic make up in t(8;21) acute myeloid leukemia (AML) cells, another subtype of AMLs, revealed that these regions are occupied by the oncofusion protein AML1-ETO. Together, our results suggest that oncofusion proteins localize to accessible regions and that chromatin accessibility together with p300 binding and histone acetylation characterize AML1-ETO and PML-RARα binding sites.