Implications of noncoding regulatory functions in the development of insulinomas.

Insulinomas are rare neuroendocrine tumors arising from pancreatic ß cells, characterized by aberrant proliferation and altered insulin secretion, leading to glucose homeostasis failure. With the aim of uncovering the role of noncoding regulatory regions and their aberrations in the development of these tumors, we coupled epigenetic and transcriptome profiling with whole-genome sequencing. As a result, we unraveled somatic mutations associated with changes in regulatory functions. Critically, these regions impact insulin secretion, tumor development, and epigenetic modifying genes, including polycomb complex components. Chromatin remodeling is apparent in insulinoma-selective domains shared across patients, containing a specific set of regulatory sequences dominated by the SOX17 binding motif. Moreover, many of these regions are H3K27me3 repressed in ß cells, suggesting that tumoral transition involves derepression of polycomb-targeted domains. Our work provides a compendium of aberrant cis-regulatory elements affecting the function and fate of ß cells in their progression to insulinomas and a framework to identify coding and noncoding driver mutations.

Strand-resolved mutagenicity of DNA damage and repair.

DNA base damage is a major source of oncogenic mutations1. Such damage can produce strand-phased mutation patterns and multiallelic variation through the process of lesion segregation2. Here we exploited these properties to reveal how strand-asymmetric processes, such as replication and transcription, shape DNA damage and repair. Despite distinct mechanisms of leading and lagging strand replication3,4, we observe identical fidelity and damage tolerance for both strands. For small alkylation adducts of DNA, our results support a model in which the same translesion polymerase is recruited on-the-fly to both replication strands, starkly contrasting the strand asymmetric tolerance of bulky UV-induced adducts5. The accumulation of multiple distinct mutations at the site of persistent lesions provides the means to quantify the relative efficiency of repair processes genome wide and at single-base resolution. At multiple scales, we show DNA damage-induced mutations are largely shaped by the influence of DNA accessibility on repair efficiency, rather than gradients of DNA damage. Finally, we reveal specific genomic conditions that can actively drive oncogenic mutagenesis by corrupting the fidelity of nucleotide excision repair. These results provide insight into how strand-asymmetric mechanisms underlie the formation, tolerance and repair of DNA damage, thereby shaping cancer genome evolution.

Identification of Clonal Hematopoiesis Driver Mutations through In Silico Saturation Mutagenesis.

Clonal hematopoiesis (CH) is a phenomenon of clonal expansion of hematopoietic stem cells driven by somatic mutations affecting certain genes. Recently, CH has been linked to the development of hematologic malignancies, cardiovascular diseases, and other conditions. Although the most frequently mutated CH driver genes have been identified, a systematic landscape of the mutations capable of initiating this phenomenon is still lacking. Here, we trained machine-learning models for 12 of the most recurrent CH genes to identify their driver mutations. These models outperform expert-curated rules based on prior knowledge of the function of these genes. Moreover, their application to identify CH driver mutations across almost half a million donors of the UK Biobank reproduces known associations between CH driver mutations and age, and the prevalence of several diseases and conditions. We thus propose that these models support the accurate identification of CH across healthy individuals.

Genomic deletions explain the generation of alternative BRAF isoforms conferring resistance to MAPK inhibitors in melanoma.

Resistance to MAPK inhibitors (MAPKi), the main cause of relapse in BRAF-mutant melanoma, is associated with the production of alternative BRAF mRNA isoforms (altBRAFs) in up to 30% of patients receiving BRAF inhibitor monotherapy. These altBRAFs have been described as being generated by alternative pre-mRNA splicing, and splicing modulation has been proposed as a therapeutic strategy to overcome resistance. In contrast, we report that altBRAFs are generated through genomic deletions. Using different in vitro models of altBRAF-mediated melanoma resistance, we demonstrate the production of altBRAFs exclusively from the BRAF V600E allele, correlating with corresponding genomic deletions. Genomic deletions are also detected in tumor samples from melanoma and breast cancer patients expressing altBRAFs. Along with the identification of altBRAFs in BRAF wild-type and in MAPKi-naive melanoma samples, our results represent a major shift in our understanding of mechanisms leading to the generation of BRAF transcripts variants associated with resistance in melanoma.

Origins of Second Malignancies in Children and Mutational Footprint of Chemotherapy in Normal Tissues.

Pediatric cancers are rare diseases, and children without known germline predisposing conditions who develop a second malignancy during developmental ages are extremely rare. We present four such clinical cases and, through whole-genome and error-correcting ultra-deep duplex sequencing of tumor and normal samples, we explored the origin of the second malignancy in four children, uncovering different routes of development. The exposure to cytotoxic therapies was linked to the emergence of a secondary acute myeloid leukemia. A common somatic mutation acquired early during embryonic development was the driver of two solid malignancies in another child. In two cases, the two tumors developed from completely independent clones diverging during embryogenesis. Importantly, we demonstrate that platinum-based therapies contributed at least one order of magnitude more mutations per day of exposure than aging to normal tissues in these children. SIGNIFICANCE: Using whole-genome and error-correcting ultra-deep duplex sequencing, we uncover different origins for second neoplasms in four children. We also uncover the presence of platinum-related mutations across 10 normal tissues of exposed individuals, highlighting the impact that the use of cytotoxic therapies may have on cancer survivors. See related commentary by Pacyna and Nangalia, p. 900. This article is featured in Selected Articles from This Issue, p. 897.

Modulating the immune response to SARS-CoV-2 by different nanocarriers delivering an mRNA expressing trimeric RBD of the spike protein: COVARNA Consortium.

Vaccines based on mRNA technology have revolutionized the field. In fact, lipid nanoparticles (LNP) formulated with mRNA are the preferential vaccine platform used in the fight against SARS-CoV-2 infection, with wider application against other diseases. The high demand and property right protection of the most potent cationic/ionizable lipids used for LNP formulation of COVID-19 mRNA vaccines have promoted the design of alternative nanocarriers for nucleic acid delivery. In this study we have evaluated the immunogenicity and efficacy of different rationally designed lipid and polymeric-based nanoparticle prototypes against SARS-CoV-2 infection. An mRNA coding for a trimeric soluble form of the receptor binding domain (RBD) of the spike (S) protein from SARS-CoV-2 was encapsulated using different components to form nanoemulsions (NE), nanocapsules (NC) and lipid nanoparticles (LNP). The toxicity and biological activity of these prototypes were evaluated in cultured cells after transfection and in mice following homologous prime/boost immunization. Our findings reveal good levels of RBD protein expression with most of the formulations. In C57BL/6 mice immunized intramuscularly with two doses of formulated RBD-mRNA, the modified lipid nanoparticle (mLNP) and the classical lipid nanoparticle (LNP-1) were the most effective delivery nanocarriers at inducing binding and neutralizing antibodies against SARS-CoV-2. Both prototypes fully protected susceptible K18-hACE2 transgenic mice from morbidity and mortality following a SARS-CoV-2 challenge. These results highlight that modulation of mRNAs immunogenicity can be achieved by using alternative nanocarriers and support further assessment of mLNP and LNP-1 prototypes as delivery vehicles for mRNA vaccines.

Hotspot propensity across mutational processes.

The sparsity of mutations observed across tumours hinders our ability to study mutation rate variability at nucleotide resolution. To circumvent this, here we investigated the propensity of mutational processes to form mutational hotspots as a readout of their mutation rate variability at single base resolution. Mutational signatures 1 and 17 have the highest hotspot propensity (5-78 times higher than other processes). After accounting for trinucleotide mutational probabilities, sequence composition and mutational heterogeneity at 10 Kbp, most (94-95%) signature 17 hotspots remain unexplained, suggesting a significant role of local genomic features. For signature 1, the inclusion of genome-wide distribution of methylated CpG sites into models can explain most (80-100%) of the hotspot propensity. There is an increased hotspot propensity of signature 1 in normal tissues and de novo germline mutations. We demonstrate that hotspot propensity is a useful readout to assess the accuracy of mutation rate models at nucleotide resolution. This new approach and the findings derived from it open up new avenues for a range of somatic and germline studies investigating and modelling mutagenesis.

Cancer Evolution: A Multifaceted Affair.

Cancer cells adapt and survive through the acquisition and selection of molecular modifications. This process defines cancer evolution. Building on a theoretical framework based on heritable genetic changes has provided insights into the mechanisms supporting cancer evolution. However, cancer hallmarks also emerge via heritable nongenetic mechanisms, including epigenetic and chromatin topological changes, and interactions between tumor cells and the tumor microenvironment. Recent findings on tumor evolutionary mechanisms draw a multifaceted picture where heterogeneous forces interact and influence each other while shaping tumor progression. A comprehensive characterization of the cancer evolutionary toolkit is required to improve personalized medicine and biomarker discovery. SIGNIFICANCE: Tumor evolution is fueled by multiple enabling mechanisms. Importantly, genetic instability, epigenetic reprogramming, and interactions with the tumor microenvironment are neither alternative nor independent evolutionary mechanisms. As demonstrated by findings highlighted in this perspective, experimental and theoretical approaches must account for multiple evolutionary mechanisms and their interactions to ultimately understand, predict, and steer tumor evolution.

Mutational signature dynamics shaping the evolution of oesophageal adenocarcinoma.

A variety of mutational processes drive cancer development, but their dynamics across the entire disease spectrum from pre-cancerous to advanced neoplasia are poorly understood. We explore the mutagenic processes shaping oesophageal adenocarcinoma tumorigenesis in 997 instances comprising distinct stages of this malignancy, from Barrett Oesophagus to primary tumours and advanced metastatic disease. The mutational landscape is dominated by the C[T > C/G]T substitution enriched signatures SBS17a/b, which are linked with TP53 mutations, increased proliferation, genomic instability and disease progression. The APOBEC mutagenesis signature is a weak but persistent signal amplified in primary tumours. We also identify prevalent alterations in DNA damage repair pathways, with homologous recombination, base and nucleotide excision repair and translesion synthesis mutated in up to 50% of the cohort, and surprisingly uncoupled from transcriptional activity. Among these, the presence of base excision repair deficiencies show remarkably poor prognosis in the cohort. In this work, we provide insights on the mutational aetiology and changes enabling the transition from pre-neoplastic to advanced oesophageal adenocarcinoma.

Copy number footprints of platinum-based anticancer therapies.

Recently, distinct mutational footprints observed in metastatic tumors, secondary malignancies and normal human tissues have been demonstrated to be caused by the exposure to several chemotherapeutic drugs. These characteristic mutations originate from specific lesions caused by these chemicals to the DNA of exposed cells. However, it is unknown whether the exposure to these chemotherapies leads to a specific footprint of larger chromosomal aberrations. Here, we address this question exploiting whole genome sequencing data of metastatic tumors obtained from patients exposed to different chemotherapeutic drugs. As a result, we discovered a specific copy number footprint across tumors from patients previously exposed to platinum-based therapies. This footprint is characterized by a significant increase in the number of chromosomal fragments of copy number 1-4 and size smaller than 10 Mb in exposed tumors with respect to their unexposed counterparts (median 14-387% greater across tumor types). The number of chromosomal fragments characteristic of the platinum-associated CN footprint increases significantly with the activity of the well known platinum-related footprint of single nucleotide variants across exposed tumors.

Genomics improves risk stratification of adults with T-cell acute lymphoblastic leukemia enrolled in measurable residual disease-oriented trials.

Genetic information has been crucial to understand the pathogenesis of T-cell acute lymphoblastic leukemia (T-ALL) at diagnosis and at relapse, but still nowadays has a limited value in a clinical context. Few genetic markers are associated with the outcome of T-ALL patients, independently of measurable residual disease (MRD) status after therapy. In addition, the prognostic relevance of genetic features may be modulated by the specific treatment used. We analyzed the genetic profile of 145 T-ALL patients by targeted deep sequencing. Genomic information was integrated with the clinicalbiological and survival data of a subset of 116 adult patients enrolled in two consecutive MRD-oriented trials of the Spanish PETHEMA (Programa Español de Tratamientos en Hematología) group. Genetic analysis revealed a mutational profile defined by DNMT3A/ N/KRAS/ MSH2/ U2AF1 gene mutations that identified refractory/resistant patients. Mutations in the DMNT3A gene were also found in the non-leukemic cell fraction of patients with T-ALL, revealing a possible mutational-driven clonal hematopoiesis event to prime T-ALL in elderly. The prognostic impact of this adverse genetic profile was independent of MRD status on day +35 of induction therapy. The combined worse-outcome genetic signature and MRD on day +35 allowed risk stratification of T-ALL into standard or high-risk groups with significantly different 5- year overall survival (OS) of 52% (95% confidence interval: 37-67) and 17% (95% confidence interval: 1-33), respectively. These results confirm the relevance of the tumor genetic profile in predicting patient outcome in adult T-ALL and highlight the need for novel gene-targeted chemotherapeutic schedules to improve the OS of poor-prognosis T-ALL patients.

Assessment of Human SARS CoV-2-Specific T-Cell Responses Elicited In Vitro by New Computationally Designed mRNA Immunogens (COVARNA).

The COVID-19 pandemic has brought significant changes and advances in the field of vaccination, including the implementation and widespread use of encapsidated mRNA vaccines in general healthcare practice. Here, we present two new mRNAs expressing antigenic parts of the SARS-CoV-2 spike protein and provide data supporting their functionality. The first mRNA, called RBD-mRNA, encodes a trimeric form of the virus spike protein receptor binding domain (RBD). The other mRNA, termed T-mRNA, codes for the relevant HLA I and II spike epitopes. The two mRNAs (COVARNA mRNAs) were designed to be used for delivery to cells in combination, with the RBD-mRNA being the primary source of antigen and the T-mRNA working as an enhancer of immunogenicity by supporting CD4 and CD8 T-cell activation. This innovative approach substantially differs from other available mRNA vaccines, which are largely directed to antibody production by the entire spike protein. In this study, we first show that both mRNAs are functionally transfected into human antigen-presenting cells (APCs). We obtained peripheral blood mononuclear cell (PBMC) samples from three groups of voluntary donors differing in their immunity against SARS-CoV-2: non-infected (naïve), infected-recovered (convalescent), and vaccinated. Using an established method of co-culturing autologous human dendritic cells (hDCs) with T-cells, we detected proliferation and cytokine secretion, thus demonstrating the ability of the COVARNA mRNAs to activate T-cells in an antigen-specific way. Interestingly, important differences in the intensity of the response between the infected-recovered (convalescent) and vaccinated donors were observed, with the levels of T-cell proliferation and cytokine secretion (IFNγ, IL-2R, and IL-13) being higher in the vaccinated group. In summary, our data support the further study of these mRNAs as a combined approach for future use as a vaccine.

Identification of Clonal Hematopoiesis Driver Mutations through In Silico Saturation Mutagenesis.

Clonal hematopoiesis (CH) is a phenomenon of clonal expansion of hematopoietic stem cells driven by somatic mutations affecting certain genes. Recently, CH has been linked to the development of a number of hematologic malignancies, cardiovascular diseases and other conditions. Although the most frequently mutated CH driver genes have been identified, a systematic landscape of the mutations capable of initiating this phenomenon is still lacking. Here, we train high-quality machine-learning models for 12 of the most recurrent CH driver genes to identify their driver mutations. These models outperform an experimental base-editing approach and expert-curated rules based on prior knowledge of the function of these genes. Moreover, their application to identify CH driver mutations across almost half a million donors of the UK Biobank reproduces known associations between CH driver mutations and age, and the prevalence of several diseases and conditions. We thus propose that these models support the accurate identification of CH across healthy individuals.

High p16 expression and heterozygous RB1 loss are biomarkers for CDK4/6 inhibitor resistance in ER+ breast cancer.

CDK4/6 inhibitors combined with endocrine therapy have demonstrated higher antitumor activity than endocrine therapy alone for the treatment of advanced estrogen receptor-positive breast cancer. Some of these tumors are de novo resistant to CDK4/6 inhibitors and others develop acquired resistance. Here, we show that p16 overexpression is associated with reduced antitumor activity of CDK4/6 inhibitors in patient-derived xenografts (n = 37) and estrogen receptor-positive breast cancer cell lines, as well as reduced response of early and advanced breast cancer patients to CDK4/6 inhibitors (n = 89). We also identified heterozygous RB1 loss as biomarker of acquired resistance and poor clinical outcome. Combination of the CDK4/6 inhibitor ribociclib with the PI3K inhibitor alpelisib showed antitumor activity in estrogen receptor-positive non-basal-like breast cancer patient-derived xenografts, independently of PIK3CA, ESR1 or RB1 mutation, also in drug de-escalation experiments or omitting endocrine therapy. Our results offer insights into predicting primary/acquired resistance to CDK4/6 inhibitors and post-progression therapeutic strategies.

Detection of early seeding of Richter transformation in chronic lymphocytic leukemia.

Richter transformation (RT) is a paradigmatic evolution of chronic lymphocytic leukemia (CLL) into a very aggressive large B cell lymphoma conferring a dismal prognosis. The mechanisms driving RT remain largely unknown. We characterized the whole genome, epigenome and transcriptome, combined with single-cell DNA/RNA-sequencing analyses and functional experiments, of 19 cases of CLL developing RT. Studying 54 longitudinal samples covering up to 19 years of disease course, we uncovered minute subclones carrying genomic, immunogenetic and transcriptomic features of RT cells already at CLL diagnosis, which were dormant for up to 19 years before transformation. We also identified new driver alterations, discovered a new mutational signature (SBS-RT), recognized an oxidative phosphorylation (OXPHOS)high-B cell receptor (BCR)low-signaling transcriptional axis in RT and showed that OXPHOS inhibition reduces the proliferation of RT cells. These findings demonstrate the early seeding of subclones driving advanced stages of cancer evolution and uncover potential therapeutic targets for RT.