Diverse mutational landscapes in human lymphocytes

Somatic mutations in aging and disease

Time always leaves its mark, and our genome is no exception. Mutations in the genome of somatic cells were first hypothesized to be the cause of aging in the 1950s, shortly after the molecular structure of DNA had been described. Somatic mutation theories of aging are based on the fact that mutations in DNA as the ultimate template for all cellular functions are irreversible. However, it took until the 1990s to develop the methods to test if DNA mutations accumulate with age in different organs and tissues and estimate the severity of the problem. By now, numerous studies have documented the accumulation of somatic mutations with age in normal cells and tissues of mice, humans, and other animals, showing clock-like mutational signatures that provide information on the underlying causes of the mutations. In this review, we will first briefly discuss the recent advances in next-generation sequencing that now allow quantitative analysis of somatic mutations. Second, we will provide evidence that the mutation rate differs between cell types, with a focus on differences between germline and somatic mutation rate. Third, we will discuss somatic mutational signatures as measures of aging, environmental exposure, and activities of DNA repair processes. Fourth, we will explain the concept of clonally amplified somatic mutations, with a focus on clonal hematopoiesis. Fifth, we will briefly discuss somatic mutations in the transcriptome and in our other genome, i.e., the genome of mitochondria. We will end with a brief discussion of a possible causal contribution of somatic mutations to the aging process.

Mutational mechanisms in multiply relapsed pediatric acute lymphoblastic leukemia

Pediatric acute lymphoblastic leukemia (ALL) is marked by low mutational load at initial diagnosis, which increases at relapse. To determine which processes are active in (relapsed) ALL and how they behave during disease progression before and after therapy, we performed whole genome sequencing on 97 tumor samples of 29 multiply relapsed ALL patients. Mutational load increased upon relapse in 28 patients and upon every subsequent relapse in 22 patients. In addition to two clock-like mutational processes, we identified UV-like damage, APOBEC activity, reactive oxygen species, thiopurine-associated damage and an unknown therapy component as drivers of mutagenesis. Mutational processes often affected patients over longer time periods, but could also occur in isolated events, suggesting the requirement of additional triggers. Thiopurine exposure was the most prominent source of new mutations in relapse, affecting over half of the studied patients in first and/or later relapse and causing potential relapse-driving mutations in multiple patients. Our data demonstrate that multiple mutational processes frequently act in parallel as prominent secondary drivers with dynamic activity during ALL development and progression.

Genomic characterization of AML with aberrations of chromosome 7: a multinational cohort of 519 patients

BACKGROUND

Deletions and partial losses of chromosome 7 (chr7) are frequent in acute myeloid leukemia (AML) and are linked to dismal outcome. However, the genomic landscape and prognostic impact of concomitant genetic aberrations remain incompletely understood.

METHODS

To discover genetic lesions in adult AML patients with aberrations of chromosome 7 [abn(7)], 60 paired diagnostic/remission samples were investigated by whole-exome sequencing in the exploration cohort. Subsequently, a gene panel including 66 genes and a SNP backbone for copy-number variation detection was designed and applied to the remaining samples of the validation cohort. In total, 519 patients were investigated, of which 415 received intensive induction treatment, typically containing a combination of cytarabine and anthracyclines.

RESULTS

In the exploration cohort, the most frequently mutated gene was TP53 (33%), followed by epigenetic regulators (DNMT3A, KMT2C, IDH2) and signaling genes (NRAS, PTPN11). Thirty percent of 519 patients harbored ≥ 1 mutation in genes located in commonly deleted regions of chr7-most frequently affecting KMT2C (16%) and EZH2 (10%). KMT2C mutations were often subclonal and enriched in patients with del(7q), de novo or core-binding factor AML (45%). Cancer cell fraction analysis and reconstruction of mutation acquisition identified TP53 mutations as mainly disease-initiating events, while del(7q) or -7 appeared as subclonal events in one-third of cases. Multivariable analysis identified five genetic lesions with significant prognostic impact in intensively treated AML patients with abn(7). Mutations in TP53 and PTPN11 (11%) showed the strongest association with worse overall survival (OS, TP53: hazard ratio [HR], 2.53 [95% CI 1.66-3.86]; P < 0.001; PTPN11: HR, 2.24 [95% CI 1.56-3.22]; P < 0.001) and relapse-free survival (RFS, TP53: HR, 2.3 [95% CI 1.25-4.26]; P = 0.008; PTPN11: HR, 2.32 [95% CI 1.33-4.04]; P = 0.003). By contrast, IDH2-mutated patients (9%) displayed prolonged OS (HR, 0.51 [95% CI 0.30-0.88]; P = 0.0015) and durable responses (RFS: HR, 0.5 [95% CI 0.26-0.96]; P = 0.036).

CONCLUSION

This work unraveled formerly underestimated genetic lesions and provides a comprehensive overview of the spectrum of recurrent gene mutations and their clinical relevance in AML with abn(7). KMT2C mutations are among the most frequent gene mutations in this heterogeneous AML subgroup and warrant further functional investigation.

Transient Differentiation-State Plasticity Occurs during Acute Lymphoblastic Leukemia Initiation

Leukemia is characterized by oncogenic lesions that result in a block of differentiation, whereas phenotypic plasticity is retained. A better understanding of how these two phenomena arise during leukemogenesis in humans could help inform diagnosis and treatment strategies. Here, we leveraged the well-defined differentiation states during T-cell development to pinpoint the initiation of T-cell acute lymphoblastic leukemia (T-ALL), an aggressive form of childhood leukemia, and study the emergence of phenotypic plasticity. Single-cell whole genome sequencing of leukemic blasts was combined with multiparameter flow cytometry to couple cell identity and clonal lineages. Irrespective of genetic events, leukemia-initiating cells altered their phenotypes by differentiation and dedifferentiation. The construction of the phylogenies of individual leukemias using somatic mutations revealed that phenotypic diversity is reflected by the clonal structure of cancer. The analysis also indicated that the acquired phenotypes are heritable and stable. Together, these results demonstrate a transient period of plasticity during leukemia initiation, where phenotypic switches seem unidirectional. Significance: A method merging multicolor flow cytometry with single-cell whole genome sequencing to couple cell identity with clonal lineages uncovers differentiation-state plasticity in leukemia, reconciling blocked differentiation with phenotypic plasticity in cancer.

The biological and clinical impact of deletions before and after large chromosomal gains in multiple myeloma

ABSTRACT

Acquisition of a hyperdiploid (HY) karyotype or immunoglobulin heavy chain (IgH) translocations are considered key initiating events in multiple myeloma (MM). To explore if other genomic events can precede these events, we analyzed whole-genome sequencing data from 1173 MM samples. By integrating molecular time and structural variants within early chromosomal duplications, we indeed identified pregain deletions in 9.4% of patients with an HY karyotype without IgH translocations, challenging acquisition of an HY karyotype as the earliest somatic event. Remarkably, these deletions affected tumor suppressor genes (TSGs) and/or oncogenes in 2.4% of patients with an HY karyotype without IgH translocations, supporting their role in MM pathogenesis. Furthermore, our study points to postgain deletions as novel driver mechanisms in MM. Using multiomics approaches to investigate their biologic impact, we found associations with poor clinical outcome in newly diagnosed patients and profound effects on both the oncogene and TSG activity despite the diploid gene status. Overall, this study provides novel insights into the temporal dynamics of genomic alterations in MM.

Motive and opportunity: MYC rearrangements in high-grade B-cell lymphoma with MYC and BCL2 rearrangements (an LLMPP study)

ABSTRACT

Rearrangements that place the oncogenes MYC, BCL2, or BCL6 adjacent to superenhancers are common in mature B-cell lymphomas. Lymphomas with diffuse large B-cell lymphoma (DLBCL) or high-grade morphology with both MYC and BCL2 rearrangements are classified as high-grade B-cell lymphoma with MYC and BCL2 rearrangements ("double hit"; HGBCL-DH-BCL2) and are associated with aggressive disease and poor outcomes. Although it is established that MYC rearrangements involving immunoglobulin (IG) loci are associated with inferior outcomes relative to those involving other non-IG superenhancers, the frequency of and mechanisms driving IG vs non-IG MYC rearrangements have not been elucidated. Here, we used custom targeted capture and/or whole-genome sequencing to characterize oncogene rearrangements across 883 mature B-cell lymphomas including Burkitt lymphoma, follicular lymphoma, DLBCL, and HGBCL-DH-BCL2 tumors. We demonstrate that, although BCL2 rearrangement topology is consistent across entities, HGBCL-DH-BCL2 have distinct MYC rearrangement architecture relative to tumors with single MYC rearrangements or with both MYC and BCL6 rearrangements (HGBCL-DH-BCL6), including both a higher frequency of non-IG rearrangements and different architecture of MYC::IGH rearrangements. The distinct MYC rearrangement patterns in HGBCL-DH-BCL2 occur on the background of high levels of somatic hypermutation across MYC partner loci in HGBCL-DH-BCL2, creating more opportunity to form these rearrangements. Furthermore, because 1 IGH allele is already disrupted by the existing BCL2 rearrangement, the MYC rearrangement architecture in HGBCL-DH-BCL2 likely reflects selective pressure to preserve both BCL2 and B-cell receptor expression. These data provide new mechanistic explanations for the distinct patterns of MYC rearrangements observed across different lymphoma entities.

Complement C3d enables protective immunity capable of distinguishing spontaneously transformed from non-transformed cells

Immune-surveillance depends in part on the recognition of peptide variants by T cell antigen receptors. Given that both normal B cells and malignant B cells accumulate mutations we chose a murine model of multiple myeloma to test conditions to induce cell-mediated immunity targeting malignant plasma cell (PC) clones but sparing of normal PCs. Revealing a novel function for intracellular C3d, we discovered that C3d engaged T cell responses against malignant plasma cells in the bone marrow of mice that had developed multiple myeloma spontaneously. Our results show that C3d internalized by cells augments immune surveillance by several mechanisms. In one, C3d induces a master transcription regulator, E2f1, to increase the expression of long non-coding (lnc) RNAs, to generate peptides for MHC-I presentation and increase MHC-I expression. In another, C3d increases expression of RNAs encoding ribosomal proteins linked to processing of defective ribosomal products (DRiPs) that arise from non-canonical translation and known to promote immunosurveillance. Cancer cells are uniquely susceptible to increased expression and presentation of mutant peptides given the extent of protein misfolding and accumulation of somatic mutations. Accordingly, although C3d can be internalized by any cell, C3d preferentially targets malignant clones by evoking specific T cell mediated immunity (CMI) and sparing most non-transformed polyclonal B cells and plasma cells with lower mutation loads. Malignant plasma cell deletion was blocked by cyclosporin or by CD8 depletion confirming that endogenous T cells mediated malignant clone clearance. Besides the potential for therapeutic application our results highlight how intracellular C3d modifies cellular metabolism to augment immune surveillance.

One Sentence Summary

We show that intracellular soluble fragment 3d of complement (C3d) induces regression of spontaneous multiple myeloma in mice reducing tumor burden by 10 fold, after 8 weeks. C3d enables cell-mediated immunity to target multiple myeloma clones sparing non-transformed polyclonal B cells and plasma cells with lower mutation loads. We show that C3d increases the expression of ribosomal subunits associated with the translation of defective ribosomal products (DRiPs). C3d also decreases expression of protein arginine methyl transferase (PRMT) 5 which in turn relieves E2f1 repression increasing the expression of Lnc RNAs and derived peptides that evoke anti-tumor cellular immunity. The approach increases MHC-I expression by tumor cells and generates a CMI response that overcomes tumor immune-evasion strategies.

Significance

Tumors are immunogenic in part because of somatic mutations that originate novel peptides that once presented on MHC engage cell-mediated immunity (CMI). However, in spite of the higher mutation load most tumors evade immunity. We discovered that a component of the complement system (C3d) overcomes tumor immune evasion by augmenting expression of ribosomal proteins and lncRNAs linked to the presentation of novel peptides by tumor cells. C3d induced CMI targets cancer cells sparing non transformed cells uncovering a novel function for complement in immune surveillance.

Selective pressures of platinum compounds shape the evolution of therapy-related myeloid neoplasms

Therapy-related myeloid neoplasms (t-MN) arise as a complication of chemo- and/or radiotherapy. Although t-MN can occur both in adult and childhood cancer survivors, the mechanisms driving therapy-related leukemogenesis likely vary across different ages. Chemotherapy is thought to induce driver mutations in children, whereas in adults pre-existing mutant clones are selected by the exposure. However, selective pressures induced by chemotherapy early in life are less well studied. Here, we use single-cell whole genome sequencing and phylogenetic inference to show that the founding cell of t-MN in children starts expanding after cessation of platinum exposure. In patients with Li-Fraumeni syndrome, characterized by a germline TP53 mutation, we find that the t-MN already expands during treatment, suggesting that platinum-induced growth inhibition is TP53-dependent. Our results demonstrate that germline aberrations can interact with treatment exposures in inducing t-MN, which is important for the development of more targeted, patient-specific treatment regimens and follow-up.

Taming AID mutator activity in somatic hypermutation

Activation-induced cytidine deaminase (AID) initiates somatic hypermutation (SHM) by introducing base substitutions into antibody genes, a process enabling antibody affinity maturation in immune response. How a mutator is tamed to precisely and safely generate programmed DNA lesions in a physiological process remains unsettled, as its dysregulation drives lymphomagenesis. Recent research has revealed several hidden features of AID-initiated mutagenesis: preferential activity on flexible DNA substrates, restrained activity within chromatin loop domains, unique DNA repair factors to differentially decode AID-caused lesions, and diverse consequences of aberrant deamination. Here, we depict the multifaceted regulation of AID activity with a focus on emerging concepts/factors and discuss their implications for the design of base editors (BEs) that install somatic mutations to correct deleterious genomic variants.

Early-life tobacco exposure is causally implicated in aberrant RAG-mediated recombination in childhood acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) is the most common cancer in children, yet few environmental risk factors have been identified. We previously found an association between early-life tobacco smoke exposure and frequency of somatic deletions of 8 leukemia driver genes among childhood ALL patients in the California Childhood Leukemia Study. To expand analysis genome-wide and examine potential mechanisms, we conducted tumor whole-genome sequencing in 35 ALL patients, including 18 with high prenatal tobacco exposure and 17 with low exposure as determined by established epigenetic biomarkers. High tobacco exposure patients had significantly more structural variants (P < .001) and deletions (P = .001) genome-wide than low exposure patients. Investigation of off-target RAG recombination revealed that 41% of deletions in the high tobacco exposure patients were putatively RAG-mediated (full RAG motif identified at one or both breakpoints) compared with only 21% in the low exposure group (P = .001). In a multilevel model, deletions in high tobacco exposure patients were 2.44-fold (95% CI:1.13-5.38) more likely to be putatively RAG-mediated than deletions in low exposure patients. No point mutational signatures were associated with prenatal tobacco exposure. Our findings suggest that early-life tobacco smoke exposure may promote leukemogenesis by driving development of somatic deletions in pre-leukemic lymphocytes via off-target RAG recombination.

Clonal associations between lymphocyte subsets and functional states in rheumatoid arthritis synovium

Rheumatoid arthritis (RA) is an autoimmune disease involving antigen-specific T and B cells. Here, we perform single-cell RNA and repertoire sequencing on paired synovial tissue and blood samples from 12 seropositive RA patients. We identify clonally expanded CD4 + T cells, including CCL5+ cells and T peripheral helper (Tph) cells, which show a prominent transcriptomic signature of recent activation and effector function. CD8 + T cells show higher oligoclonality than CD4 + T cells, with the largest synovial clones enriched in GZMK+ cells. CD8 + T cells with possibly virus-reactive TCRs are distributed across transcriptomic clusters. In the B cell compartment, NR4A1+ activated B cells, and plasma cells are enriched in the synovium and demonstrate substantial clonal expansion. We identify synovial plasma cells that share BCRs with synovial ABC, memory, and activated B cells. Receptor-ligand analysis predicted IFNG and TNFRSF members as mediators of synovial Tph-B cell interactions. Together, these results reveal clonal relationships between functionally distinct lymphocyte populations that infiltrate the synovium of patients with RA.

Somatic mutations associate with clonal expansion of CD8+ T cells

Somatic mutations in T cells can cause cancer but also have implications for immunological diseases and cell therapies. The mutation spectrum in nonmalignant T cells is unclear. Here, we examined somatic mutations in CD4+ and CD8+ T cells from 90 patients with hematological and immunological disorders and used T cell receptor (TCR) and single-cell sequencing to link mutations with T cell expansions and phenotypes. CD8+ cells had a higher mutation burden than CD4+ cells. Notably, the biggest variant allele frequency (VAF) of non-synonymous variants was higher than synonymous variants in CD8+ T cells, indicating non-random occurrence. The non-synonymous VAF in CD8+ T cells strongly correlated with the TCR frequency, but not age. We identified mutations in pathways essential for T cell function and often affected lymphoid neoplasia. Single-cell sequencing revealed cytotoxic TEMRA phenotypes of mutated T cells. Our findings suggest that somatic mutations contribute to CD8+ T cell expansions without malignant transformation.

Analysis of somatic mutations in whole blood from 200,618 individuals identifies pervasive positive selection and novel drivers of clonal hematopoiesis

Human aging is marked by the emergence of a tapestry of clonal expansions in dividing tissues, particularly evident in blood as clonal hematopoiesis (CH). CH, linked to cancer risk and aging-related phenotypes, often stems from somatic mutations in a set of established genes. However, the majority of clones lack known drivers. Here we infer gene-level positive selection in whole blood exomes from 200,618 individuals in UK Biobank. We identify 17 additional genes, ZBTB33, ZNF318, ZNF234, SPRED2, SH2B3, SRCAP, SIK3, SRSF1, CHEK2, CCDC115, CCL22, BAX, YLPM1, MYD88, MTA2, MAGEC3 and IGLL5, under positive selection at a population level, and validate this selection pattern in 10,837 whole genomes from single-cell-derived hematopoietic colonies. Clones with mutations in these genes grow in frequency and size with age, comparable to classical CH drivers. They correlate with heightened risk of infection, death and hematological malignancy, highlighting the significance of these additional genes in the aging process.

Developmental origins shape the paediatric cancer genome

In the past two decades, technological advances have brought unprecedented insights into the paediatric cancer genome revealing characteristics distinct from those of adult cancer. Originating from developing tissues, paediatric cancers generally have low mutation burden and are driven by variants that disrupt the transcriptional activity, chromatin state, non-coding cis-regulatory regions and other biological functions. Within each tumour, there are multiple populations of cells with varying states, and the lineages of some can be tracked to their fetal origins. Genome-wide genetic screening has identified vulnerabilities associated with both the cell of origin and transcription deregulation in paediatric cancer, which have become a valuable resource for designing new therapeutic approaches including those for small molecules, immunotherapy and targeted protein degradation. In this Review, we present recent findings on these facets of paediatric cancer from a pan-cancer perspective and provide an outlook on future investigations.

Extrinsic and intrinsic drivers of natural killer cell clonality

Clonal expansion of antigen-specific lymphocytes is the fundamental mechanism enabling potent adaptive immune responses and the generation of immune memory. Accompanied by pronounced epigenetic remodeling, the massive proliferation of individual cells generates a critical mass of effectors for the control of acute infections, as well as a pool of memory cells protecting against future pathogen encounters. Classically associated with the adaptive immune system, recent work has demonstrated that innate immune memory to human cytomegalovirus (CMV) infection is stably maintained as large clonal expansions of natural killer (NK) cells, raising questions on the mechanisms for clonal selection and expansion in the absence of re-arranged antigen receptors. Here, we discuss clonal NK cell memory in the context of the mechanisms underlying clonal competition of adaptive lymphocytes and propose alternative selection mechanisms that might decide on the clonal success of their innate counterparts. We propose that the integration of external cues with cell-intrinsic sources of heterogeneity, such as variegated receptor expression, transcriptional states, and somatic variants, compose a bottleneck for clonal selection, contributing to the large size of memory NK cell clones.

Contrasting somatic mutation patterns in aging human neurons and oligodendrocytes

Characterizing somatic mutations in the brain is important for disentangling the complex mechanisms of aging, yet little is known about mutational patterns in different brain cell types. Here, we performed whole-genome sequencing (WGS) of 86 single oligodendrocytes, 20 mixed glia, and 56 single neurons from neurotypical individuals spanning 0.4-104 years of age and identified >92,000 somatic single-nucleotide variants (sSNVs) and small insertions/deletions (indels). Although both cell types accumulate somatic mutations linearly with age, oligodendrocytes accumulated sSNVs 81% faster than neurons and indels 28% slower than neurons. Correlation of mutations with single-nucleus RNA profiles and chromatin accessibility from the same brains revealed that oligodendrocyte mutations are enriched in inactive genomic regions and are distributed across the genome similarly to mutations in brain cancers. In contrast, neuronal mutations are enriched in open, transcriptionally active chromatin. These stark differences suggest an assortment of active mutagenic processes in oligodendrocytes and neurons.

DNA flexibility can shape the preferential hypermutation of antibody genes

Antibody-coding genes accumulate somatic mutations to achieve antibody affinity maturation. Genetic dissection using various mouse models has shown that intrinsic hypermutations occur preferentially and are predisposed in the DNA region encoding antigen-contacting residues. The molecular basis of nonrandom/preferential mutations is a long-sought question in the field. Here, we summarize recent findings on how single-strand (ss)DNA flexibility facilitates activation-induced cytidine deaminase (AID) activity and fine-tunes the mutation rates at a mesoscale within the antibody variable domain exon. We propose that antibody coding sequences are selected based on mutability during the evolution of adaptive immunity and that DNA mechanics play a noncoding role in the genome. The mechanics code may also determine other cellular DNA metabolism processes, which awaits future investigation.

Deciphering cell states and genealogies of human haematopoiesis

The human blood system is maintained through the differentiation and massive amplification of a limited number of long-lived haematopoietic stem cells (HSCs)1. Perturbations to this process underlie diverse diseases, but the clonal contributions to human haematopoiesis and how this changes with age remain incompletely understood. Although recent insights have emerged from barcoding studies in model systems2-5, simultaneous detection of cell states and phylogenies from natural barcodes in humans remains challenging. Here we introduce an improved, single-cell lineage-tracing system based on deep detection of naturally occurring mitochondrial DNA mutations with simultaneous readout of transcriptional states and chromatin accessibility. We use this system to define the clonal architecture of HSCs and map the physiological state and output of clones. We uncover functional heterogeneity in HSC clones, which is stable over months and manifests as both differences in total HSC output and biases towards the production of different mature cell types. We also find that the diversity of HSC clones decreases markedly with age, leading to an oligoclonal structure with multiple distinct clonal expansions. Our study thus provides a clonally resolved and cell-state-aware atlas of human haematopoiesis at single-cell resolution, showing an unappreciated functional diversity of human HSC clones and, more broadly, paving the way for refined studies of clonal dynamics across a range of tissues in human health and disease.

Impact of Rare Structural Variant Events in Newly Diagnosed Multiple Myeloma

PURPOSE

Whole-genome sequencing (WGS) of patients with newly diagnosed multiple myeloma (NDMM) has shown recurrent structural variant (SV) involvement in distinct regions of the genome (i.e., hotspots) and causing recurrent copy-number alterations. Together with canonical immunoglobulin translocations, these SVs are recognized as "recurrent SVs." More than half of SVs were not involved in recurrent events. The significance of these "rare SVs" has not been previously examined.

EXPERIMENTAL DESIGN

In this study, we utilize 752 WGS and 591 RNA sequencing data from patients with NDMM to determine the role of rare SVs in myeloma pathogenesis.

RESULTS

Ninety-four percent of patients harbored at least one rare SV event. Rare SVs showed an SV class-specific enrichment within genes and superenhancers associated with outlier gene expression. Furthermore, known myeloma driver genes recurrently impacted by point mutations were dysregulated by rare SVs.

CONCLUSIONS

Overall, we demonstrate the association of rare SVs with aberrant gene expression supporting a potential driver role in myeloma pathogenesis.

APOBEC Mutagenesis in Cancer Development and Susceptibility

APOBEC cytosine deaminases are prominent mutators in cancer, mediating mutations in over 50% of cancers. APOBEC mutagenesis has been linked to tumor heterogeneity, persistent cell evolution, and therapy responses. While emerging evidence supports the impact of APOBEC mutagenesis on cancer progression, the understanding of its contribution to cancer susceptibility and malignant transformation is limited. We examine the existing evidence for the role of APOBEC mutagenesis in carcinogenesis on the basis of the reported associations between germline polymorphisms in genes encoding APOBEC enzymes and cancer risk, insights into APOBEC activities from sequencing efforts of both malignant and non-malignant human tissues, and in vivo studies. We discuss key knowledge gaps and highlight possible ways to gain a deeper understanding of the contribution of APOBEC mutagenesis to cancer development.

The Diverse Roles of ETV6 Alterations in B-Lymphoblastic Leukemia and Other Hematopoietic Cancers

Genetic alterations of the repressive ETS family transcription factor gene ETV6 are recurrent in several categories of hematopoietic malignancy, including subsets of B-cell and T-cell acute lymphoblastic leukemias (B-ALL and T-ALL), myeloid neoplasms, and mature B-cell lymphomas. ETV6 is essential for adult hematopoietic stem cells (HSCs), contributes to specific functions of some mature immune cells, and plays a key role in thrombopoiesis as demonstrated by familial ETV6 mutations associated with thrombocytopenia and predisposition to hematopoietic cancers, particularly B-ALL. ETV6 appears to have a tumor suppressor role in several hematopoietic lineages, as demonstrated by recurrent somatic loss-of-function (LoF) and putative dominant-negative alterations in leukemias and lymphomas. ETV6 rearrangements contribute to recurrent fusion oncogenes such as the B-ALL-associated transcription factor (TF) fusions ETV6::RUNX1 and PAX5::ETV6, rare drivers such as ETV6::NCOA6, and a spectrum of tyrosine kinase gene fusions encoding hyperactive signaling proteins that self-associate via the ETV6 N-terminal pointed domain. Another subset of recurrent rearrangements involving the ETV6 gene locus appear to function primarily to drive overexpression of the partner gene. This review surveys what is known about the biochemical and genome regulatory properties of ETV6 as well as our current understanding of how alterations in these functions contribute to hematopoietic and nonhematopoietic cancers.

Clonal selection of hematopoietic stem cells after gene therapy for sickle cell disease

Gene therapy (GT) provides a potentially curative treatment option for patients with sickle cell disease (SCD); however, the occurrence of myeloid malignancies in GT clinical trials has prompted concern, with several postulated mechanisms. Here, we used whole-genome sequencing to track hematopoietic stem cells (HSCs) from six patients with SCD at pre- and post-GT time points to map the somatic mutation and clonal landscape of gene-modified and unmodified HSCs. Pre-GT, phylogenetic trees were highly polyclonal and mutation burdens per cell were elevated in some, but not all, patients. Post-GT, no clonal expansions were identified among gene-modified or unmodified cells; however, an increased frequency of potential driver mutations associated with myeloid neoplasms or clonal hematopoiesis (DNMT3A- and EZH2-mutated clones in particular) was observed in both genetically modified and unmodified cells, suggesting positive selection of mutant clones during GT. This work sheds light on HSC clonal dynamics and the mutational landscape after GT in SCD, highlighting the enhanced fitness of some HSCs harboring pre-existing driver mutations. Future studies should define the long-term fate of mutant clones, including any contribution to expansions associated with myeloid neoplasms.

Quantifying human genome parameters in aging

Healthy human longevity is a global goal of the world health system. Determining the causes and processes influencing human longevity is the primary fundamental goal facing the scientific community. Currently, the main efforts of the scientific community are aimed at identifying the qualitative characteristics of the genome that determine the trait. At the same time, when evaluating qualitative characteristics, there are many challenges that make it difficult to establish associations. Quantitative traits are burdened with such problems to a lesser extent, but they are largely overlooked in current genomic studies of aging and longevity. Although there is a wide repertoire of quantitative trait analyses based on genomic data, most opportunities are ignored by authors, which, along with the inaccessibility of published data, leads to the loss of this important information. This review focuses on describing quantitative traits important for understanding aging and necessary for analysis in further genomic studies, and recommends the inclusion of the described traits in the analysis. The review considers the relationship between quantitative characteristics of the mitochondrial genome and aging, longevity, and age-related neurodegenerative diseases, such as the frequency of extensive mitochondrial DNA (mtDNA) deletions, mtDNA half-life, the frequency of A>G replacements in the mtDNA heavy chain, the number of mtDNA copies; special attention is paid to the mtDNA methylation sign. A separate section of this review is devoted to the correlation of telomere length parameters with age, as well as the association of telomere length with the amount of mitochondrial DNA. In addition, we consider such a quantitative feature as the rate of accumulation of somatic mutations with aging in relation to the lifespan of living organisms. In general, it may be noted that there are quite serious reasons to suppose that various quantitative characteristics of the genome may be directly or indirectly associated with certain aspects of aging and longevity. At the same time, the available data are clearly insufficient for definitive conclusions and the determination of causal relationships.

Relapse Timing Is Associated With Distinct Evolutionary Dynamics in Diffuse Large B-Cell Lymphoma

PURPOSE

Diffuse large B-cell lymphoma (DLBCL) is cured in more than 60% of patients, but outcomes remain poor for patients experiencing disease progression or relapse (refractory or relapsed DLBCL [rrDLBCL]), particularly if these events occur early. Although previous studies examining cohorts of rrDLBCL have identified features that are enriched at relapse, few have directly compared serial biopsies to uncover biological and evolutionary dynamics driving rrDLBCL. Here, we sought to confirm the relationship between relapse timing and outcomes after second-line (immuno)chemotherapy and determine the evolutionary dynamics that underpin that relationship.

PATIENTS AND METHODS

Outcomes were examined in a population-based cohort of 221 patients with DLBCL who experienced progression/relapse after frontline treatment and were treated with second-line (immuno)chemotherapy with an intention-to-treat with autologous stem-cell transplantation (ASCT). Serial DLBCL biopsies from a partially overlapping cohort of 129 patients underwent molecular characterization, including whole-genome or whole-exome sequencing in 73 patients.

RESULTS

Outcomes to second-line therapy and ASCT are superior for late relapse (>2 years postdiagnosis) versus primary refractory (<9 months) or early relapse (9-24 months). Diagnostic and relapse biopsies were mostly concordant for cell-of-origin classification and genetics-based subgroup. Despite this concordance, the number of mutations exclusive to each biopsy increased with time since diagnosis, and late relapses shared few mutations with their diagnostic counterpart, demonstrating a branching evolution pattern. In patients with highly divergent tumors, many of the same genes acquired new mutations independently in each tumor, suggesting that the earliest mutations in a shared precursor cell constrain tumor evolution toward the same genetics-based subgroups at both diagnosis and relapse.

CONCLUSION

These results suggest that late relapses commonly represent genetically distinct and chemotherapy-naïve disease and have implications for optimal patient management.

Molecular mechanisms insulating proliferation from genotoxic stress in B lymphocytes

In mammals, B cells strictly segregate proliferation from somatic mutation as they develop within the bone marrow and then mature through germinal centers (GCs) in the periphery. Failure to do so risks autoimmunity and neoplastic transformation. Recent work has described how B cell progenitors transition between proliferation and mutation via cytokine signaling pathways, epigenetic chromatin regulation, and remodeling of 3D chromatin conformation. We propose a three-zone model of the GC that describes how proliferation and mutation are regulated. Using this model, we consider how recent mechanistic discoveries in B cell progenitors inform models of GC B cell function and reveal fundamental mechanisms underpinning humoral immunity, autoimmunity, and lymphomagenesis.

Convergent somatic evolution commences in utero in a germline ribosomopathy

Clonal tracking of cells using somatic mutations permits exploration of clonal dynamics in human disease. Here, we perform whole genome sequencing of 323 haematopoietic colonies from 10 individuals with the inherited ribosomopathy Shwachman-Diamond syndrome to reconstruct haematopoietic phylogenies. In ~30% of colonies, we identify mutually exclusive mutations in TP53, EIF6, RPL5, RPL22, PRPF8, plus chromosome 7 and 15 aberrations that increase SBDS and EFL1 gene dosage, respectively. Target gene mutations commence in utero, resulting in a profusion of clonal expansions, with only a few haematopoietic stem cell lineages (mean 8, range 1-24) contributing ~50% of haematopoietic colonies across 8 individuals (range 4-100% clonality) by young adulthood. Rapid clonal expansion during disease transformation is associated with biallelic TP53 mutations and increased mutation burden. Our study highlights how convergent somatic mutation of the p53-dependent nucleolar surveillance pathway offsets the deleterious effects of germline ribosomopathy but increases opportunity for TP53-mutated cancer evolution.

Diagnostic utility of whole genome sequencing in adults with B-other acute lymphoblastic leukemia

Genomic profiling during the diagnosis of B-cell precursor acute lymphoblastic leukemia (BCP-ALL) in adults is used to guide disease classification, risk stratification, and treatment decisions. Patients for whom diagnostic screening fails to identify disease-defining or risk-stratifying lesions are classified as having B-other ALL. We screened a cohort of 652 BCP-ALL cases enrolled in UKALL14 to identify and perform whole genome sequencing (WGS) of paired tumor-normal samples. For 52 patients with B-other, we compared the WGS findings with data from clinical and research cytogenetics. WGS identified a cancer-associated event in 51 of 52 patients, including an established subtype defining genetic alterations that were previously missed with standard-of-care (SoC) genetics in 5 of them. Of the 47 true B-other ALL, we identified a recurrent driver in 87% (41). A complex karyotype via cytogenetics emerges as a heterogeneous group, including distinct genetic alterations associated with either favorable (DUX4-r) or poor outcomes (MEF2D-r and IGK::BCL2). For a subset of 31 cases, we integrated the findings from RNA sequencing (RNA-seq) analysis to include fusion gene detection and classification based on gene expression. Compared with RNA-seq, WGS was sufficient to detect and resolve recurrent genetic subtypes; however, RNA-seq can provide orthogonal validation of findings. In conclusion, we demonstrated that WGS can identify clinically relevant genetic abnormalities missed with SoC testing as well as identify leukemia driver events in virtually all cases of B-other ALL.

Thermoresponsive Polypeptide Fused L-Asparaginase with Mitigated Immunogenicity and Enhanced Efficacy in Treating Hematologic Malignancies

L-Asparaginase (ASP) is well-known for its excellent efficacy in treating hematological malignancies. Unfortunately, the intrinsic shortcomings of ASP, namely high immunogenicity, severe toxicity, short half-life, and poor stability, restrict its clinical usage. Poly(ethylene glycol) conjugation (PEGylation) of ASP is an effective strategy to address these issues, but it is not ideal in clinical applications due to complex chemical synthesis procedures, reduced ASP activity after conjugation, and pre-existing anti-PEG antibodies in humans. Herein, the authors genetically engineered an elastin-like polypeptide (ELP)-fused ASP (ASP-ELP), a core-shell structured tetramer predicted by AlphaFold2, to overcome the limitations of ASP and PEG-ASP. Notably, the unique thermosensitivity of ASP-ELP enables the in situ formation of a sustained-release depot post-injection with zero-order release kinetics over a long time. The in vitro and in vivo studies reveal that ASP-ELP possesses increased activity retention, improved stability, extended half-life, mitigated immunogenicity, reduced toxicity, and enhanced efficacy compared to ASP and PEG-ASP. Indeed, ASP-ELP treatment in leukemia or lymphoma mouse models of cell line-derived xenograft (CDX) shows potent anti-cancer effects with significantly prolonged survival. The findings also indicate that artificial intelligence (AI)-assisted genetic engineering is instructive in designing protein-polypeptide conjugates and may pave the way to develop next-generation biologics to enhance cancer treatment.

Clonal origin of KMT2A wild-type lineage-switch leukemia following CAR-T cell and blinatumomab therapy

Children with acute lymphoblastic leukemia (ALL) undergoing anti-CD19 therapy occasionally develop acute myeloid leukemia (AML). The clonal origin of such lineage-switch leukemias1-4 remains unresolved. Here, we reconstructed the phylogeny of multiple leukemias in a girl who, following multiply relapsed ALL, received anti-CD19 cellular and antibody treatment and subsequently developed AML. Whole genome sequencing unambiguously revealed the AML derived from the initial ALL, with distinct driver mutations that were detectable before emergence. Extensive prior diversification and subsequent clonal selection underpins this fatal lineage switch. Genomic monitoring of primary leukemias and recurrences may predict therapy resistance, especially regarding anti-CD19 treatment.

The landscape of somatic mutations in lymphoblastoid cell lines

Somatic mutations have important biological ramifications while exerting substantial rate, type, and genomic location heterogeneity. Yet, their sporadic occurrence makes them difficult to study at scale and across individuals. Lymphoblastoid cell lines (LCLs), a model system for human population and functional genomics, harbor large numbers of somatic mutations and have been extensively genotyped. By comparing 1,662 LCLs, we report that the mutational landscape of the genome varies across individuals in terms of the number of mutations, their genomic locations, and their spectra; this variation may itself be modulated by somatic trans-acting mutations. Mutations attributed to the translesion DNA polymerase η follow two different modes of formation, with one mode accounting for the hypermutability of the inactive X chromosome. Nonetheless, the distribution of mutations along the inactive X chromosome appears to follow an epigenetic memory of the active form.

Cell-type specificity of the human mutation landscape with respect to DNA replication dynamics

The patterns of genomic mutations are associated with various genomic features, most notably late replication timing, yet it remains contested which mutation types and signatures relate to DNA replication dynamics and to what extent. Here, we perform high-resolution comparisons of mutational landscapes between lymphoblastoid cell lines, chronic lymphocytic leukemia tumors, and three colon adenocarcinoma cell lines, including two with mismatch repair deficiency. Using cell-type-matched replication timing profiles, we demonstrate that mutation rates exhibit heterogeneous replication timing associations among cell types. This cell-type heterogeneity extends to the underlying mutational pathways, as mutational signatures show inconsistent replication timing bias between cell types. Moreover, replicative strand asymmetries exhibit similar cell-type specificity, albeit with different relationships to replication timing than mutation rates. Overall, we reveal an underappreciated complexity and cell-type specificity of mutational pathways and their relationship to replication timing.

Ultraviolet radiation shapes dendritic cell leukaemia transformation in the skin

Tumours most often arise from progression of precursor clones within a single anatomical niche. In the bone marrow, clonal progenitors can undergo malignant transformation to acute leukaemia, or differentiate into immune cells that contribute to disease pathology in peripheral tissues1-4. Outside the marrow, these clones are potentially exposed to a variety of tissue-specific mutational processes, although the consequences of this are unclear. Here we investigate the development of blastic plasmacytoid dendritic cell neoplasm (BPDCN)-an unusual form of acute leukaemia that often presents with malignant cells isolated to the skin5. Using tumour phylogenomics and single-cell transcriptomics with genotyping, we find that BPDCN arises from clonal (premalignant) haematopoietic precursors in the bone marrow. We observe that BPDCN skin tumours first develop at sun-exposed anatomical sites and are distinguished by clonally expanded mutations induced by ultraviolet (UV) radiation. A reconstruction of tumour phylogenies reveals that UV damage can precede the acquisition of alterations associated with malignant transformation, implicating sun exposure of plasmacytoid dendritic cells or committed precursors during BPDCN pathogenesis. Functionally, we find that loss-of-function mutations in Tet2, the most common premalignant alteration in BPDCN, confer resistance to UV-induced cell death in plasmacytoid, but not conventional, dendritic cells, suggesting a context-dependent tumour-suppressive role for TET2. These findings demonstrate how tissue-specific environmental exposures at distant anatomical sites can shape the evolution of premalignant clones to disseminated cancer.

Molecular Evolution of Classic Hodgkin Lymphoma Revealed Through Whole-Genome Sequencing of Hodgkin and Reed Sternberg Cells

The rarity of malignant Hodgkin and Reed Sternberg (HRS) cells in classic Hodgkin lymphoma (cHL) limits the ability to study the genomics of cHL. To circumvent this, our group has previously optimized fluorescence-activated cell sorting to purify HRS cells. Using this approach, we now report the whole-genome sequencing landscape of HRS cells and reconstruct the chronology and likely etiology of pathogenic events leading to cHL. We identified alterations in driver genes not previously described in cHL, APOBEC mutational activity, and the presence of complex structural variants including chromothripsis. We found that high ploidy in cHL is often acquired through multiple, independent chromosomal gains events including whole-genome duplication. Evolutionary timing analyses revealed that structural variants enriched for RAG motifs, driver mutations in B2M, BCL7A, GNA13, and PTPN1, and the onset of AID-driven mutagenesis usually preceded large chromosomal gains. This study provides a temporal reconstruction of cHL pathogenesis.

SIGNIFICANCE

Previous studies in cHL were limited to coding sequences and therefore not able to comprehensively decipher the tumor complexity. Here, leveraging cHL whole-genome characterization, we identify driver events and reconstruct the tumor evolution, finding that structural variants, driver mutations, and AID mutagenesis precede chromosomal gains. This article is highlighted in the In This Issue feature, p. 171.

Relapse timing is associated with distinct evolutionary dynamics in DLBCL

Diffuse large B-cell lymphoma (DLBCL) is cured in over 60% of patients, but outcomes are poor for patients with relapsed or refractory disease (rrDLBCL). Here, we performed whole genome/exome sequencing (WGS/WES) on tumors from 73 serially-biopsied patients with rrDLBCL. Based on the observation that outcomes to salvage therapy/autologous stem cell transplantation are related to time-to-relapse, we stratified patients into groups according to relapse timing to explore the relationship to genetic divergence and sensitivity to salvage immunochemotherapy. The degree of mutational divergence increased with time between biopsies, yet tumor pairs were mostly concordant for cell-of-origin, oncogene rearrangement status and genetics-based subgroup. In patients with highly divergent tumors, several genes acquired exclusive mutations independently in each tumor, which, along with concordance of genetics-based subgroups, suggests that the earliest mutations in a shared precursor cell constrain tumor evolution. These results suggest that late relapses commonly represent genetically distinct and chemotherapy-naïve disease.

Lymphoma-associated mutations in autoreactive memory B cells of patients with Sjögren's syndrome

We recently demonstrated that normal memory B lymphocytes carry a substantial number of de novo mutations in the genome. Here, we performed exome-wide somatic mutation analyses of bona fide autoreactive rheumatoid factor (RF)-expressing memory B cells retrieved from patients with Sjӧgren's syndrome (SS). The amount and repertoire of the de novo exome mutations of RF B cells were found to be essentially different from those detected in healthy donor memory B cells. In contrast to the mutation spectra of normal B cells, which appeared random and non-selected, the mutations of the RF B cells were greater in number and enriched for mutations in genes also found mutated in B-cell non-Hodgkin lymphomas. During the study, one of the SS patients developed a diffuse large B-cell lymphoma (DLBCL) out of an RF clone that was identified 2 years earlier in an inflamed salivary gland biopsy. The successive oncogenic events in the RF precursor clone and the DLBCL were assessed. In conclusion, our findings of enhanced and selected genomic damage in growth-regulating genes in RF memory B cells of SS patients together with the documented transformation of an RF-precursor clone into DLBCL provide unique novel insight into the earliest stages of B-cell derailment and lymphomagenesis. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Detection of chromosomal alteration after infusion of gene-edited allogeneic CAR T cells

A chromosome 14 inversion was found in a patient who developed bone marrow aplasia following treatment with allogeneic chimeric antigen receptor (CAR) Tcells containing gene edits made with transcription activator-like effector nucleases (TALEN). TALEN editing sites were not involved at either breakpoint. Recombination signal sequences (RSSs) were found suggesting recombination-activating gene (RAG)-mediated activity. The inversion represented a dominant clone detected in the context of decreasing absolute CAR Tcell and overall lymphocyte counts. The inversion was not associated with clinical consequences and wasnot detected in the drug product administered to this patient or in any drug product used in this or other trials using the same manufacturing processes. Neither was the inversion detected in this patient at earlier time points or in any other patient enrolled in this or other trials treated with this or other product lots. This case illustrates that spontaneous, possibly RAG-mediated, recombination events unrelated to gene editing can occur in adoptive cell therapy studies, emphasizes the need for ruling out off-target gene editing sites, and illustrates that other processes, such as spontaneous V(D)J recombination, can lead to chromosomal alterations in infused cells independent of gene editing.

Dr.Nod: computational framework for discovery of regulatory non-coding drivers in tissue-matched distal regulatory elements

The discovery of cancer driver mutations is a fundamental goal in cancer research. While many cancer driver mutations have been discovered in the protein-coding genome, research into potential cancer drivers in the non-coding regions showed limited success so far. Here, we present a novel comprehensive framework Dr.Nod for detection of non-coding cis-regulatory candidate driver mutations that are associated with dysregulated gene expression using tissue-matched enhancer-gene annotations. Applying the framework to data from over 1500 tumours across eight tissues revealed a 4.4-fold enrichment of candidate driver mutations in regulatory regions of known cancer driver genes. An overarching conclusion that emerges is that the non-coding driver mutations contribute to cancer by significantly altering transcription factor binding sites, leading to upregulation of tissue-matched oncogenes and down-regulation of tumour-suppressor genes. Interestingly, more than half of the detected cancer-promoting non-coding regulatory driver mutations are over 20 kb distant from the cancer-associated genes they regulate. Our results show the importance of tissue-matched enhancer-gene maps, functional impact of mutations, and complex background mutagenesis model for the prediction of non-coding regulatory drivers. In conclusion, our study demonstrates that non-coding mutations in enhancers play a previously underappreciated role in cancer and dysregulation of clinically relevant target genes.

Contrasting patterns of somatic mutations in neurons and glia reveal differential predisposition to disease in the aging human brain

Characterizing the mechanisms of somatic mutations in the brain is important for understanding aging and disease, but little is known about the mutational patterns of different cell types. We performed whole-genome sequencing of 71 oligodendrocytes and 51 neurons from neurotypical individuals (0.4 to 104 years old) and identified >67,000 somatic single nucleotide variants (sSNVs) and small insertions and deletions (indels). While both cell types accumulate mutations with age, oligodendrocytes accumulate sSNVs 69% faster than neurons (27/year versus 16/year) whereas indels accumulate 42% slower (1.8/year versus 3.1/year). Correlation with single-cell RNA and chromatin accessibility from the same brains revealed that oligodendrocyte mutations are enriched in inactive genomic regions and are distributed similarly to mutations in brain cancers. In contrast, neuronal mutations are enriched in open, transcriptionally active chromatin. These patterns highlight differences in the mutagenic processes in glia and neurons and suggest cell type-specific, age-related contributions to neurodegeneration and oncogenesis.

Amoeboid migration in health and disease: Immune responses versus cancer dissemination

Cell migration is crucial for efficient immune responses and is aberrantly used by cancer cells during metastatic dissemination. Amoeboid migrating cells use myosin II-powered blebs to propel themselves, and change morphology and direction. Immune cells use amoeboid strategies to respond rapidly to infection or tissue damage, which require quick passage through several barriers, including blood, lymph and interstitial tissues, with complex and varied environments. Amoeboid migration is also used by metastatic cancer cells to aid their migration, dissemination and survival, whereby key mechanisms are hijacked from professionally motile immune cells. We explore important parallels observed between amoeboid immune and cancer cells. We also consider key distinctions that separate the lifespan, state and fate of these cell types as they migrate and/or fulfil their function. Finally, we reflect on unexplored areas of research that would enhance our understanding of how tumour cells use immune cell strategies during metastasis, and how to target these processes.

Impact of rare structural variant events in newly diagnosed multiple myeloma

Whole genome sequencing (WGS) of newly diagnosed multiple myeloma patients (NDMM) has shown recurrent structural variant (SV) involvement in distinct regions of the genome (i.e. hotspots) and causing recurrent copy number alterations. Together with canonical immunoglobulin translocations, these SVs are recognized as "recurrent SVs". More than half SVs were not involved in recurrent events. The significance of these "rare SVs" has not been previously examined. In this study, we utilize 752 WGS and 591 RNA-seq data from NDMM patients to determine the role of rare SVs in myeloma pathogenesis. 94% of patients harbored at least one rare SV event. Rare SVs showed an SV-class specific enrichment within genes and superenhancers associated with outlier gene expression. Furthermore, known myeloma driver genes recurrently impacted by point mutations were dysregulated by rare SVs. Overall, we demonstrate the association of rare SVs with aberrant gene expression supporting a driver role in myeloma pathogenesis.

SIGNIFICANCE

Characterization of multiple myeloma genome revealed that more than half structural variants are not involved in recurrent events. Here, we demonstrate that these rare SVs hold potential for myeloma pathogenesis through their gene expression impact. Rare SVs contribute to MM heterogeneity and have implications for development of individualized treatment.

Frequent somatic mosaicism in T lymphocyte subsets in individuals with and without multiple sclerosis

Background

Somatic variants are variations in an individual's genome acquired after the zygotic stadium and result from mitotic errors or not (fully) repaired DNA damage.

Objectives

To investigate whether somatic mosaicism in T lymphocyte subsets is enriched early in multiple sclerosis (MS).

Methods

We identified somatic variants with variant allele fractions ≥1% across the whole exome in CD4⁺ and CD8⁺ T lymphocytes of 21 treatment-naive MS patients with <5 years of disease duration and 16 partially age-matched healthy controls. We investigated the known somatic STAT3 variant p.Y640F in peripheral blood in a larger cohort of 446 MS patients and 259 controls.

Results

All subjects carried 1-142 variants in CD4⁺ or CD8⁺ T lymphocytes. Variants were more common, more abundant, and increased with age in CD8⁺ T lymphocytes. Somatic variants were common in the genes DNMT3A and especially STAT3. Overall, the presence or abundance of somatic variants, including the STAT3 p.Y640F variant, did not differ between MS patients and controls.

Conclusions

Somatic variation in T lymphocyte subsets is widespread in both control individuals and MS patients. Somatic mosaicism in T lymphocyte subsets is not enriched in early MS and thus unlikely to contribute to MS risk, but future research needs to address whether a subset of variants influences disease susceptibility.

Age-related somatic mutation burden in human tissues

The genome of multicellular organisms carries the hereditary information necessary for the development of all organs and tissues and to maintain function in adulthood. To ensure the genetic stability of the species, genomes are protected against changes in sequence information. However, genomes are not static. De novo mutations in germline cells are passed on to offspring and generate the variation needed in evolution. Moreover, postzygotic mutations occur in all somatic cells during development and aging. These somatic mutations remain limited to the individual, generating tissues that are genome mosaics. Insight into such mutations and their consequences has been limited due to their extremely low abundance, with most mutations unique for each cell. Recent advances in sequencing, including whole genome sequencing at the single-cell level, have now led to the first insights into somatic mutation burdens in human tissues. Here, we will first briefly describe the latest methodology for somatic mutation analysis, then review our current knowledge of somatic mutation burden in human tissues and, finally, briefly discuss the possible functional impact of somatic mutations on the aging process and age-related diseases, including cancer and diseases other than cancer.