Passenger Mutations in More Than 2,500 Cancer Genomes: Overall Molecular Functional Impact and Consequences.

The dichotomous model of "drivers" and "passengers" in cancer posits that only a few mutations in a tumor strongly affect its progression, with the remaining ones being inconsequential. Here, we leveraged the comprehensive variant dataset from the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) project to demonstrate that-in addition to the dichotomy of high- and low-impact variants-there is a third group of medium-impact putative passengers. Moreover, we also found that molecular impact correlates with subclonal architecture (i.e., early versus late mutations), and different signatures encode for mutations with divergent impact. Furthermore, we adapted an additive-effects model from complex-trait studies to show that the aggregated effect of putative passengers, including undetected weak drivers, provides significant additional power (∼12% additive variance) for predicting cancerous phenotypes, beyond PCAWG-identified driver mutations. Finally, this framework allowed us to estimate the frequency of potential weak-driver mutations in PCAWG samples lacking any well-characterized driver alterations.

Shared Immunogenic Poly-Epitope Frameshift Mutations in Microsatellite Unstable Tumors.

Microsatellite instability-high (MSI-H) tumors are characterized by high tumor mutation burden and responsiveness to checkpoint blockade. We identified tumor-specific frameshifts encoding multiple epitopes that originated from indel mutations shared among patients with MSI-H endometrial, colorectal, and stomach cancers. Epitopes derived from these shared frameshifts have high population occurrence rates, wide presence in many tumor subclones, and are predicted to bind to the most frequent MHC alleles in MSI-H patient cohorts. Neoantigens arising from these mutations are distinctly unlike self and viral antigens, signifying novel groups of potentially highly immunogenic tumor antigens. We further confirmed the immunogenicity of frameshift peptides in T cell stimulation experiments using blood mononuclear cells isolated from both healthy donors and MSI-H cancer patients. Our study uncovers the widespread occurrence and strong immunogenicity of tumor-specific antigens derived from shared frameshift mutations in MSI-H cancer and Lynch syndrome patients, suitable for the design of common "off-the-shelf" cancer vaccines.

A Compendium of Mutational Signatures of Environmental Agents.

Whole-genome-sequencing (WGS) of human tumors has revealed distinct mutation patterns that hint at the causative origins of cancer. We examined mutational signatures in 324 WGS human-induced pluripotent stem cells exposed to 79 known or suspected environmental carcinogens. Forty-one yielded characteristic substitution mutational signatures. Some were similar to signatures found in human tumors. Additionally, six agents produced double-substitution signatures and eight produced indel signatures. Investigating mutation asymmetries across genome topography revealed fully functional mismatch and transcription-coupled repair pathways. DNA damage induced by environmental mutagens can be resolved by disparate repair and/or replicative pathways, resulting in an assortment of signature outcomes even for a single agent. This compendium of experimentally induced mutational signatures permits further exploration of roles of environmental agents in cancer etiology and underscores how human stem cell DNA is directly vulnerable to environmental agents. VIDEO ABSTRACT.

Characterizing Mutational Signatures in Human Cancer Cell Lines Reveals Episodic APOBEC Mutagenesis.

Multiple signatures of somatic mutations have been identified in cancer genomes. Exome sequences of 1,001 human cancer cell lines and 577 xenografts revealed most common mutational signatures, indicating past activity of the underlying processes, usually in appropriate cancer types. To investigate ongoing patterns of mutational-signature generation, cell lines were cultured for extended periods and subsequently DNA sequenced. Signatures of discontinued exposures, including tobacco smoke and ultraviolet light, were not generated in vitro. Signatures of normal and defective DNA repair and replication continued to be generated at roughly stable mutation rates. Signatures of APOBEC cytidine deaminase DNA-editing exhibited substantial fluctuations in mutation rate over time with episodic bursts of mutations. The initiating factors for the bursts are unclear, although retrotransposon mobilization may contribute. The examined cell lines constitute a resource of live experimental models of mutational processes, which potentially retain patterns of activity and regulation operative in primary human cancers.

Coupled analysis of transcriptome and BCR mutations reveals role of OXPHOS in affinity maturation.

Antigen-activated B cells diversify variable regions of B cell antigen receptors by somatic hypermutation in germinal centers (GCs). The positive selection of GC B cells that acquire high-affinity mutations enables antibody affinity maturation. In spite of considerable progress, the genomic states underlying this process remain to be elucidated. Single-cell RNA sequencing and topic modeling revealed increased expression of the oxidative phosphorylation (OXPHOS) module in GC B cells undergoing mitoses. Coupled analysis of somatic hypermutation in immunoglobulin heavy chain (Igh) variable gene regions showed that GC B cells acquiring higher-affinity mutations had further elevated expression of OXPHOS genes. Deletion of mitochondrial Cox10 in GC B cells resulted in reduced cell division and impaired positive selection. Correspondingly, augmentation of OXPHOS activity with oltipraz promoted affinity maturation. We propose that elevated OXPHOS activity promotes B cell clonal expansion and also positive selection by tuning cell division times.

Genomic Patterns of De Novo Mutation in Simplex Autism.

To further our understanding of the genetic etiology of autism, we generated and analyzed genome sequence data from 516 idiopathic autism families (2,064 individuals). This resource includes >59 million single-nucleotide variants (SNVs) and 9,212 private copy number variants (CNVs), of which 133,992 and 88 are de novo mutations (DNMs), respectively. We estimate a mutation rate of ∼1.5 × 10-8 SNVs per site per generation with a significantly higher mutation rate in repetitive DNA. Comparing probands and unaffected siblings, we observe several DNM trends. Probands carry more gene-disruptive CNVs and SNVs, resulting in severe missense mutations and mapping to predicted fetal brain promoters and embryonic stem cell enhancers. These differences become more pronounced for autism genes (p = 1.8 × 10-3, OR = 2.2). Patients are more likely to carry multiple coding and noncoding DNMs in different genes, which are enriched for expression in striatal neurons (p = 3 × 10-3), suggesting a path forward for genetically characterizing more complex cases of autism.

Bedside Back to Bench: Building Bridges between Basic and Clinical Genomic Research.

Genome sequencing has revolutionized the diagnosis of genetic diseases. Close collaborations between basic scientists and clinical genomicists are now needed to link genetic variants with disease causation. To facilitate such collaborations, we recommend prioritizing clinically relevant genes for functional studies, developing reference variant-phenotype databases, adopting phenotype description standards, and promoting data sharing.

Ubiquitination-Deficient Mutations in Human Piwi Cause Male Infertility by Impairing Histone-to-Protamine Exchange during Spermiogenesis.

Genetic studies have elucidated critical roles of Piwi proteins in germline development in animals, but whether Piwi is an actual disease gene in human infertility remains unknown. We report germline mutations in human Piwi (Hiwi) in patients with azoospermia that prevent its ubiquitination and degradation. By modeling such mutations in Piwi (Miwi) knockin mice, we demonstrate that the genetic defects are directly responsible for male infertility. Mechanistically, we show that MIWI binds the histone ubiquitin ligase RNF8 in a Piwi-interacting RNA (piRNA)-independent manner, and MIWI stabilization sequesters RNF8 in the cytoplasm of late spermatids. The resulting aberrant sperm show histone retention, abnormal morphology, and severely compromised activity, which can be functionally rescued via blocking RNF8-MIWI interaction in spermatids with an RNF8-N peptide. Collectively, our findings identify Piwi as a factor in human infertility and reveal its role in regulating the histone-to-protamine exchange during spermiogenesis.

Mutational Strand Asymmetries in Cancer Genomes Reveal Mechanisms of DNA Damage and Repair.

Mutational processes constantly shape the somatic genome, leading to immunity, aging, cancer, and other diseases. When cancer is the outcome, we are afforded a glimpse into these processes by the clonal expansion of the malignant cell. Here, we characterize a less explored layer of the mutational landscape of cancer: mutational asymmetries between the two DNA strands. Analyzing whole-genome sequences of 590 tumors from 14 different cancer types, we reveal widespread asymmetries across mutagenic processes, with transcriptional ("T-class") asymmetry dominating UV-, smoking-, and liver-cancer-associated mutations and replicative ("R-class") asymmetry dominating POLE-, APOBEC-, and MSI-associated mutations. We report a striking phenomenon of transcription-coupled damage (TCD) on the non-transcribed DNA strand and provide evidence that APOBEC mutagenesis occurs on the lagging-strand template during DNA replication. As more genomes are sequenced, studying and classifying their asymmetries will illuminate the underlying biological mechanisms of DNA damage and repair.

Biochemical Basis for Dominant Inheritance, Variable Penetrance, and Maternal Effects in RBP4 Congenital Eye Disease.

Gestational vitamin A (retinol) deficiency poses a risk for ocular birth defects and blindness. We identified missense mutations in RBP4, encoding serum retinol binding protein, in three families with eye malformations of differing severity, including bilateral anophthalmia. The mutant phenotypes exhibit dominant inheritance, but incomplete penetrance. Maternal transmission significantly increases the probability of phenotypic expression. RBP normally delivers retinol from hepatic stores to peripheral tissues, including the placenta and fetal eye. The disease mutations greatly reduce retinol binding to RBP, yet paradoxically increase the affinity of RBP for its cell surface receptor, STRA6. By occupying STRA6 nonproductively, the dominant-negative proteins disrupt vitamin A delivery from wild-type proteins within the fetus, but also, in the case of maternal transmission, at the placenta. These findings establish a previously uncharacterized mode of maternal inheritance, distinct from imprinting and oocyte-derived mRNA, and define a group of hereditary disorders plausibly modulated by dietary vitamin A.

Computationally restoring the potency of a clinical antibody against Omicron.

The COVID-19 pandemic underscored the promise of monoclonal antibody-based prophylactic and therapeutic drugs1-3 and revealed how quickly viral escape can curtail effective options4,5. When the SARS-CoV-2 Omicron variant emerged in 2021, many antibody drug products lost potency, including Evusheld and its constituent, cilgavimab4-6. Cilgavimab, like its progenitor COV2-2130, is a class 3 antibody that is compatible with other antibodies in combination4 and is challenging to replace with existing approaches. Rapidly modifying such high-value antibodies to restore efficacy against emerging variants is a compelling mitigation strategy. We sought to redesign and renew the efficacy of COV2-2130 against Omicron BA.1 and BA.1.1 strains while maintaining efficacy against the dominant Delta variant. Here we show that our computationally redesigned antibody, 2130-1-0114-112, achieves this objective, simultaneously increases neutralization potency against Delta and subsequent variants of concern, and provides protection in vivo against the strains tested: WA1/2020, BA.1.1 and BA.5. Deep mutational scanning of tens of thousands of pseudovirus variants reveals that 2130-1-0114-112 improves broad potency without increasing escape liabilities. Our results suggest that computational approaches can optimize an antibody to target multiple escape variants, while simultaneously enriching potency. Our computational approach does not require experimental iterations or pre-existing binding data, thus enabling rapid response strategies to address escape variants or lessen escape vulnerabilities.

EXO5-DNA structure and BLM interactions direct DNA resection critical for ATR-dependent replication restart.

Stalled DNA replication fork restart after stress as orchestrated by ATR kinase, BLM helicase, and structure-specific nucleases enables replication, cell survival, and genome stability. Here we unveil human exonuclease V (EXO5) as an ATR-regulated DNA structure-specific nuclease and BLM partner for replication fork restart. We find that elevated EXO5 in tumors correlates with increased mutation loads and poor patient survival, suggesting that EXO5 upregulation has oncogenic potential. Structural, mechanistic, and mutational analyses of EXO5 and EXO5-DNA complexes reveal a single-stranded DNA binding channel with an adjacent ATR phosphorylation motif (T88Q89) that regulates EXO5 nuclease activity and BLM binding identified by mass spectrometric analysis. EXO5 phospho-mimetic mutant rescues the restart defect from EXO5 depletion that decreases fork progression, DNA damage repair, and cell survival. EXO5 depletion furthermore rescues survival of FANCA-deficient cells and indicates EXO5 functions epistatically with SMARCAL1 and BLM. Thus, an EXO5 axis connects ATR and BLM in directing replication fork restart.

RASGRP1 deficiency causes immunodeficiency with impaired cytoskeletal dynamics.

RASGRP1 is an important guanine nucleotide exchange factor and activator of the RAS-MAPK pathway following T cell antigen receptor (TCR) signaling. The consequences of RASGRP1 mutations in humans are unknown. In a patient with recurrent bacterial and viral infections, born to healthy consanguineous parents, we used homozygosity mapping and exome sequencing to identify a biallelic stop-gain variant in RASGRP1. This variant segregated perfectly with the disease and has not been reported in genetic databases. RASGRP1 deficiency was associated in T cells and B cells with decreased phosphorylation of the extracellular-signal-regulated serine kinase ERK, which was restored following expression of wild-type RASGRP1. RASGRP1 deficiency also resulted in defective proliferation, activation and motility of T cells and B cells. RASGRP1-deficient natural killer (NK) cells exhibited impaired cytotoxicity with defective granule convergence and actin accumulation. Interaction proteomics identified the dynein light chain DYNLL1 as interacting with RASGRP1, which links RASGRP1 to cytoskeletal dynamics. RASGRP1-deficient cells showed decreased activation of the GTPase RhoA. Treatment with lenalidomide increased RhoA activity and reversed the migration and activation defects of RASGRP1-deficient lymphocytes.

Mouse Watch: A Cautionary Tale.

In this issue of Immunity, Chisolm et al. (2019) issue a "three-alarm fire" warning to the immunology research community of unexpectedly widespread genetic variation in widely used congenic mouse strains and provide a simple method to identify such a variation through a re-analysis of existing RNA-seq and ChIP-seq datasets.

Endogenous retrotransposition activates oncogenic pathways in hepatocellular carcinoma.

LINE-1 (L1) retrotransposons are mobile genetic elements comprising ~17% of the human genome. New L1 insertions can profoundly alter gene function and cause disease, though their significance in cancer remains unclear. Here, we applied enhanced retrotransposon capture sequencing (RC-seq) to 19 hepatocellular carcinoma (HCC) genomes and elucidated two archetypal L1-mediated mechanisms enabling tumorigenesis. In the first example, 4/19 (21.1%) donors presented germline retrotransposition events in the tumor suppressor mutated in colorectal cancers (MCC). MCC expression was ablated in each case, enabling oncogenic ß-catenin/Wnt signaling. In the second example, suppression of tumorigenicity 18 (ST18) was activated by a tumor-specific L1 insertion. Experimental assays confirmed that the L1 interrupted a negative feedback loop by blocking ST18 repression of its enhancer. ST18 was also frequently amplified in HCC nodules from Mdr2(-/-) mice, supporting its assignment as a candidate liver oncogene. These proof-of-principle results substantiate L1-mediated retrotransposition as an important etiological factor in HCC.

Structures of Drosophila cryptochrome and mouse cryptochrome1 provide insight into circadian function.

Drosophila cryptochrome (dCRY) is a FAD-dependent circadian photoreceptor, whereas mammalian cryptochromes (CRY1/2) are integral clock components that repress mCLOCK/mBMAL1-dependent transcription. We report crystal structures of full-length dCRY, a dCRY loop deletion construct, and the photolyase homology region of mouse CRY1 (mCRY1). Our dCRY structures depict Phe534 of the regulatory tail in the same location as the photolesion in DNA-repairing photolyases and reveal that the sulfur loop and tail residue Cys523 plays key roles in the dCRY photoreaction. Our mCRY1 structure visualizes previously characterized mutations, an NLS, and MAPK and AMPK phosphorylation sites. We show that the FAD and antenna chromophore-binding regions, a predicted coiled-coil helix, the C-terminal lid, and charged surfaces are involved in FAD-independent mPER2 and FBXL3 binding and mCLOCK/mBMAL1 transcriptional repression. The structure of a mammalian cryptochrome1 protein may catalyze the development of CRY chemical probes and the design of therapeutic metabolic modulators.

Spatial and temporal mapping of de novo mutations in schizophrenia to a fetal prefrontal cortical network.

Genes disrupted in schizophrenia may be revealed by de novo mutations in affected persons from otherwise healthy families. Furthermore, during normal brain development, genes are expressed in patterns specific to developmental stage and neuroanatomical structure. We identified de novo mutations in persons with schizophrenia and then mapped the responsible genes onto transcriptome profiles of normal human brain tissues from age 13 weeks gestation to adulthood. In the dorsolateral and ventrolateral prefrontal cortex during fetal development, genes harboring damaging de novo mutations in schizophrenia formed a network significantly enriched for transcriptional coexpression and protein interaction. The 50 genes in the network function in neuronal migration, synaptic transmission, signaling, transcriptional regulation, and transport. These results suggest that disruptions of fetal prefrontal cortical neurogenesis are critical to the pathophysiology of schizophrenia. These results also support the feasibility of integrating genomic and transcriptome analyses to map critical neurodevelopmental processes in time and space in the brain.

Signatures of copy number alterations in human cancer.

Gains and losses of DNA are prevalent in cancer and emerge as a consequence of inter-related processes of replication stress, mitotic errors, spindle multipolarity and breakage-fusion-bridge cycles, among others, which may lead to chromosomal instability and aneuploidy1,2. These copy number alterations contribute to cancer initiation, progression and therapeutic resistance3-5. Here we present a conceptual framework to examine the patterns of copy number alterations in human cancer that is widely applicable to diverse data types, including whole-genome sequencing, whole-exome sequencing, reduced representation bisulfite sequencing, single-cell DNA sequencing and SNP6 microarray data. Deploying this framework to 9,873 cancers representing 33 human cancer types from The Cancer Genome Atlas6 revealed a set of 21 copy number signatures that explain the copy number patterns of 97% of samples. Seventeen copy number signatures were attributed to biological phenomena of whole-genome doubling, aneuploidy, loss of heterozygosity, homologous recombination deficiency, chromothripsis and haploidization. The aetiologies of four copy number signatures remain unexplained. Some cancer types harbour amplicon signatures associated with extrachromosomal DNA, disease-specific survival and proto-oncogene gains such as MDM2. In contrast to base-scale mutational signatures, no copy number signature was associated with many known exogenous cancer risk factors. Our results synthesize the global landscape of copy number alterations in human cancer by revealing a diversity of mutational processes that give rise to these alterations.

Genomics to select treatment for patients with metastatic breast cancer.

Cancer progression is driven in part by genomic alterations1. The genomic characterization of cancers has shown interpatient heterogeneity regarding driver alterations2, leading to the concept that generation of genomic profiling in patients with cancer could allow the selection of effective therapies3,4. Although DNA sequencing has been implemented in practice, it remains unclear how to use its results. A total of 1,462 patients with HER2-non-overexpressing metastatic breast cancer were enroled to receive genomic profiling in the SAFIR02-BREAST trial. Two hundred and thirty-eight of these patients were randomized in two trials (nos. NCT02299999 and NCT03386162) comparing the efficacy of maintenance treatment5 with a targeted therapy matched to genomic alteration. Targeted therapies matched to genomics improves progression-free survival when genomic alterations are classified as level I/II according to the ESMO Scale for Clinical Actionability of Molecular Targets (ESCAT)6 (adjusted hazards ratio (HR): 0.41, 90% confidence interval (CI): 0.27-0.61, P < 0.001), but not when alterations are unselected using ESCAT (adjusted HR: 0.77, 95% CI: 0.56-1.06, P = 0.109). No improvement in progression-free survival was observed in the targeted therapies arm (unadjusted HR: 1.15, 95% CI: 0.76-1.75) for patients presenting with ESCAT alteration beyond level I/II. Patients with germline BRCA1/2 mutations (n = 49) derived high benefit from olaparib (gBRCA1: HR = 0.36, 90% CI: 0.14-0.89; gBRCA2: HR = 0.37, 90% CI: 0.17-0.78). This trial provides evidence that the treatment decision led by genomics should be driven by a framework of target actionability in patients with metastatic breast cancer.

Candidate Cancer Driver Mutations in Distal Regulatory Elements and Long-Range Chromatin Interaction Networks.

A comprehensive catalog of cancer driver mutations is essential for understanding tumorigenesis and developing therapies. Exome-sequencing studies have mapped many protein-coding drivers, yet few non-coding drivers are known because genome-wide discovery is challenging. We developed a driver discovery method, ActiveDriverWGS, and analyzed 120,788 cis-regulatory modules (CRMs) across 1,844 whole tumor genomes from the ICGC-TCGA PCAWG project. We found 30 CRMs with enriched SNVs and indels (FDR < 0.05). These frequently mutated regulatory elements (FMREs) were ubiquitously active in human tissues, showed long-range chromatin interactions and mRNA abundance associations with target genes, and were enriched in motif-rewiring mutations and structural variants. Genomic deletion of one FMRE in human cells caused proliferative deficiencies and transcriptional deregulation of cancer genes CCNB1IP1, CDH1, and CDKN2B, validating observations in FMRE-mutated tumors. Pathway analysis revealed further sub-significant FMREs at cancer genes and processes, indicating an unexplored landscape of infrequent driver mutations in the non-coding genome.