The genomic landscape of 2,023 colorectal cancers.

Colorectal carcinoma (CRC) is a common cause of mortality1, but a comprehensive description of its genomic landscape is lacking2-9. Here we perform whole-genome sequencing of 2,023 CRC samples from participants in the UK 100,000 Genomes Project, thereby providing a highly detailed somatic mutational landscape of this cancer. Integrated analyses identify more than 250 putative CRC driver genes, many not previously implicated in CRC or other cancers, including several recurrent changes outside the coding genome. We extend the molecular pathways involved in CRC development, define four new common subgroups of microsatellite-stable CRC based on genomic features and show that these groups have independent prognostic associations. We also characterize several rare molecular CRC subgroups, some with potential clinical relevance, including cancers with both microsatellite and chromosomal instability. We demonstrate a spectrum of mutational profiles across the colorectum, which reflect aetiological differences. These include the role of Escherichia colipks+ colibactin in rectal cancers10 and the importance of the SBS93 signature11-13, which suggests that diet or smoking is a risk factor. Immune-escape driver mutations14 are near-ubiquitous in hypermutant tumours and occur in about half of microsatellite-stable CRCs, often in the form of HLA copy number changes. Many driver mutations are actionable, including those associated with rare subgroups (for example, BRCA1 and IDH1), highlighting the role of whole-genome sequencing in optimizing patient care.

Heterologous mRNA/MVA delivering trimeric-RBD as effective vaccination regimen against SARS-CoV-2: COVARNA Consortium.

Despite the high efficiency of current SARS-CoV-2 mRNA vaccines in reducing COVID-19 morbidity and mortality, waning immunity and the emergence of resistant variants underscore the need for novel vaccination strategies. This study explores a heterologous mRNA/Modified Vaccinia virus Ankara (MVA) prime/boost regimen employing a trimeric form of the receptor binding domain (RBD) of the SARS-CoV-2 spike (S) protein compared to a homologous MVA/MVA regimen. In C57BL/6 mice, the RBD was delivered during priming via an mRNA vector encapsulated in nanoemulsions (NE) or lipid nanoparticles (LNP), followed by a booster with a replication-deficient MVA-based recombinant virus (MVA-RBD). This heterologous mRNA/MVA regimen elicited strong anti-RBD binding and neutralizing antibodies (BAbs and NAbs) against both the ancestral SARS-CoV-2 strain and different variants of concern (VoCs). Additionally, this protocol induced robust and polyfunctional RBD-specific CD4 and CD8 T cell responses, particularly in animals primed with mLNP-RBD. In K18-hACE2 transgenic mice, the LNP-RBD/MVA combination provided complete protection from morbidity and mortality following a live SARS-CoV-2 challenge compared with the partial protection observed with mNE-RBD/MVA or MVA/MVA regimens. Although the mNE-RBD/MVA regimen only protects half of the animals, it was able to induce antibodies with Fc-mediated effector functions besides NAbs. Moreover, viral replication and viral load in the respiratory tract were markedly reduced and decreased pro-inflammatory cytokine levels were observed. These results support the efficacy of heterologous mRNA/MVA vaccine combinations over homologous MVA/MVA regimen, using alternative nanocarriers that circumvent intellectual property restrictions of current mRNA vaccine formulations.

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.

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.

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.

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.

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.

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.