SINGLE CELL DISSECTION OF DEVELOPMENTAL ORIGINS AND TRANSCRIPTIONAL HETEROGENEITY IN B-CELL ACUTE LYMPHOBLASTIC LEUKEMIA.

Sequencing of bulk tumor populations has improved genetic classification and risk assessment of B-ALL, but does not directly examine intratumor heterogeneity or infer leukemia cellular origins. We profiled 89 B-ALL samples by single-cell RNA-seq (scRNA-seq) and compared them to a reference map of normal human B-cell development established using both functional and molecular assays. Intra-sample heterogeneity was driven by cell cycle, metabolism, differentiation, and inflammation transcriptional programs. By inference of B lineage developmental state composition, nearly all samples possessed a high abundance of pro-B cells, with variation between samples mainly driven by sub-populations. However, ZNF384- r and DUX4- r B-ALL showed composition enrichment of hematopoietic stem cells, BCR::ABL1 and KMT2A -r ALL of Early Lymphoid progenitors, MEF2D -r and TCF3::PBX1 of Pre-B cells. Enrichment of Early Lymphoid progenitors correlated with high-risk clinical features. Understanding variation in transcriptional programs and developmental states of B-ALL by scRNA-seq refines existing clinical and genomic classifications and improves prediction of treatment outcome.

Single-cell chromatin accessibility profiling of acute myeloid leukemia reveals heterogeneous lineage composition upon therapy-resistance.

Acute myeloid leukemia (AML) is a heterogeneous disease characterized by high rate of therapy resistance. Since the cell of origin can impact response to therapy, it is crucial to understand the lineage composition of AML cells at time of therapy resistance. Here we leverage single-cell chromatin accessibility profiling of 22 AML bone marrow aspirates from eight patients at time of therapy resistance and following subsequent therapy to characterize their lineage landscape. Our findings reveal a complex lineage architecture of therapy-resistant AML cells that are primed for stem and progenitor lineages and spanning quiescent, activated and late stem cell/progenitor states. Remarkably, therapy-resistant AML cells are also composed of cells primed for differentiated myeloid, erythroid and even lymphoid lineages. The heterogeneous lineage composition persists following subsequent therapy, with early progenitor-driven features marking unfavorable prognosis in The Cancer Genome Atlas AML cohort. Pseudotime analysis further confirms the vast degree of heterogeneity driven by the dynamic changes in chromatin accessibility. Our findings suggest that therapy-resistant AML cells are characterized not only by stem and progenitor states, but also by a continuum of differentiated cellular lineages. The heterogeneity in lineages likely contributes to their therapy resistance by harboring different degrees of lineage-specific susceptibilities to therapy.

Transcriptomic classes of BCR-ABL1 lymphoblastic leukemia.

In BCR-ABL1 lymphoblastic leukemia, treatment heterogeneity to tyrosine kinase inhibitors (TKIs), especially in the absence of kinase domain mutations in BCR-ABL1, is poorly understood. Through deep molecular profiling, we uncovered three transcriptomic subtypes of BCR-ABL1 lymphoblastic leukemia, each representing a maturation arrest at a stage of B-cell progenitor differentiation. An earlier arrest was associated with lineage promiscuity, treatment refractoriness and poor patient outcomes. A later arrest was associated with lineage fidelity, durable leukemia remissions and improved patient outcomes. Each maturation arrest was marked by specific genomic events that control different transition points in B-cell development. Interestingly, these events were absent in BCR-ABL1+ preleukemic stem cells isolated from patients regardless of subtype, which supports that transcriptomic phenotypes are determined downstream of the leukemia-initialing event. Overall, our data indicate that treatment response and TKI efficacy are unexpected outcomes of the differentiation stage at which this leukemia transforms.

Distinct Assemblies of Heterodimeric Cytokine Receptors Govern Stemness Programs in Leukemia.

Leukemia stem cells (LSC) possess distinct self-renewal and arrested differentiation properties that are responsible for disease emergence, therapy failure, and recurrence in acute myeloid leukemia (AML). Despite AML displaying extensive biological and clinical heterogeneity, LSC with high interleukin-3 receptor (IL3R) levels are a constant yet puzzling feature, as this receptor lacks tyrosine kinase activity. Here, we show that the heterodimeric IL3Rα/ßc receptor assembles into hexamers and dodecamers through a unique interface in the 3D structure, where high IL3Rα/ßc ratios bias hexamer formation. Importantly, receptor stoichiometry is clinically relevant as it varies across the individual cells in the AML hierarchy, in which high IL3Rα/ßc ratios in LSCs drive hexamer-mediated stemness programs and poor patient survival, while low ratios mediate differentiation. Our study establishes a new paradigm in which alternative cytokine receptor stoichiometries differentially regulate cell fate, a signaling mechanism that may be generalizable to other transformed cellular hierarchies and of potential therapeutic significance. SIGNIFICANCE: Stemness is a hallmark of many cancers and is largely responsible for disease emergence, progression, and relapse. Our finding that clinically significant stemness programs in AML are directly regulated by different stoichiometries of cytokine receptors represents a hitherto unexplained mechanism underlying cell-fate decisions in cancer stem cell hierarchies. This article is highlighted in the In This Issue feature, p. 1749.

Simultaneous inhibition of Sirtuin 3 and cholesterol homeostasis targets acute myeloid leukemia stem cells by perturbing fatty acid ß-oxidation and inducing lipotoxicity.

Outcomes for patients with acute myeloid leukemia (AML) remain poor due to the inability of current therapeutic regimens to fully eradicate disease-initiating leukemia stem cells (LSC). Previous studies have demonstrated that oxidative phosphorylation (OXPHOS) is an essential process that is targetable in LSC. Sirtuin 3 (SIRT3), a mitochondrial deacetylase with a multi-faceted role in metabolic regulation, has been shown to regulate OXPHOS in cancer models; however, it has not yet been studied in the context of LSC. Thus, we sought to identify if SIRT3 is important for LSC function. Using RNAi and a SIRT3 inhibitor (YC8-02), we demonstrate that SIRT3 is a critical target for the survival of primary human LSC but is not essential for normal human hematopoietic stem and progenitor cell function. In order to elucidate the molecular mechanisms by which SIRT3 is essential in LSC we combined transcriptomic, proteomic, and lipidomic approaches, showing that SIRT3 is important for LSC function through the regulation of fatty acid oxidation (FAO) which is required to support OXPHOS and ATP production in human LSC. Further, we discovered two approaches to further sensitize LSC to SIRT3 inhibition. First, we found that LSC tolerate the toxic effects of fatty acid accumulation induced by SIRT3 inhibition by upregulating cholesterol esterification. Disruption of cholesterol homeostasis sensitizes LSC to YC8-02 and potentiates LSC death. Second, SIRT3 inhibition sensitizes LSC to the BCL-2 inhibitor venetoclax. Together, these findings establish SIRT3 as a regulator of lipid metabolism and potential therapeutic target in primitive AML cells.

KDM6 demethylases integrate DNA repair gene regulation and loss of KDM6A sensitizes human acute myeloid leukemia to PARP and BCL2 inhibition.

Acute myeloid leukemia (AML) is a heterogeneous, aggressive malignancy with dismal prognosis and with limited availability of targeted therapies. Epigenetic deregulation contributes to AML pathogenesis. KDM6 proteins are histone-3-lysine-27-demethylases that play context-dependent roles in AML. We inform that KDM6-demethylase function critically regulates DNA-damage-repair-(DDR) gene expression in AML. Mechanistically, KDM6 expression is regulated by genotoxic stress, with deficiency of KDM6A-(UTX) and KDM6B-(JMJD3) impairing DDR transcriptional activation and compromising repair potential. Acquired KDM6A loss-of-function mutations are implicated in chemoresistance, although a significant percentage of relapsed-AML has upregulated KDM6A. Olaparib treatment reduced engraftment of KDM6A-mutant-AML-patient-derived xenografts, highlighting synthetic lethality using Poly-(ADP-ribose)-polymerase-(PARP)-inhibition. Crucially, a higher KDM6A expression is correlated with venetoclax tolerance. Loss of KDM6A increased mitochondrial activity, BCL2 expression, and sensitized AML cells to venetoclax. Additionally, BCL2A1 associates with venetoclax resistance, and KDM6A loss was accompanied with a downregulated BCL2A1. Corroborating these results, dual targeting of PARP and BCL2 was superior to PARP or BCL2 inhibitor monotherapy in inducing AML apoptosis, and primary AML cells carrying KDM6A-domain mutations were even more sensitive to the combination. Together, our study illustrates a mechanistic rationale in support of a novel combination therapy for AML based on subtype-heterogeneity, and establishes KDM6A as a molecular regulator for determining therapeutic efficacy.

Precise single-cell transcriptomic mapping of normal and leukemic cell states reveals unconventional lineage priming in acute myeloid leukemia.

Initial classification of acute leukemia involves the assignment of blasts to cell states within the hematopoietic hierarchy based on morphological and immunophenotypic features. Yet, these traditional classification approaches lack precision, especially at the level of immature blasts. Single-cell RNA-sequencing (scRNA-seq) enables precise determination of cell state using thousands of markers, thus providing an opportunity to re-examine present-day classification schemes of acute leukemia. Here, we developed a detailed reference map of human bone marrow hematopoiesis from 263,519 single-cell transcriptomes spanning 55 cellular states. Cell state annotations were benchmarked against purified cell populations, and in-depth characterization of gene expression programs underlying hematopoietic differentiation was undertaken. Projection of single-cell transcriptomes from 175 samples spanning acute myeloid leukemia (AML), mixed phenotype acute leukemia (MPAL), and acute erythroid leukemia (AEL) revealed 11 subtypes involving distinct stages of hematopoietic differentiation. These included AML subtypes with notable lymphoid or erythroid lineage priming, challenging traditional diagnostic boundaries between AML, MPAL, and AEL. Quantification of lineage priming in bulk patient cohorts revealed specific genetic alterations associated with this unconventional lineage priming. Integration of transcriptional and genetic information at the single-cell level revealed how genetic subclones can induce lineage restriction, differentiation blocks, or expansion of mature myeloid cells. Furthermore, we demonstrate that distinct cellular hierarchies can co-exist within individual patients, providing insight into AML evolution in response to varying selection pressures. Together, precise mapping of hematopoietic cell states can serve as a foundation for refining disease classification in acute leukemia and understanding response or resistance to emerging therapies.

Identification and characterization of in vitro expanded hematopoietic stem cells.

Hematopoietic stem cells (HSCs) cultured outside the body are the fundamental component of a wide range of cellular and gene therapies. Recent efforts have achieved > 200-fold expansion of functional HSCs, but their molecular characterization has not been possible since the majority of cells are non-HSCs and single cell-initiated cultures have substantial clone-to-clone variability. Using the Fgd5 reporter mouse in combination with the EPCR surface marker, we report exclusive identification of HSCs from non-HSCs in expansion cultures. By directly linking single-clone functional transplantation data with single-clone gene expression profiling, we show that the molecular profile of expanded HSCs is similar to proliferating fetal HSCs and reveals a gene expression signature, including Esam, Prdm16, Fstl1, and Palld, that can identify functional HSCs from multiple cellular states. This "repopulation signature" (RepopSig) also enriches for HSCs in human datasets. Together, these findings demonstrate the power of integrating functional and molecular datasets to better derive meaningful gene signatures and opens the opportunity for a wide range of functional screening and molecular experiments previously not possible due to limited HSC numbers.

A cellular hierarchy framework for understanding heterogeneity and predicting drug response in acute myeloid leukemia.

The treatment landscape of acute myeloid leukemia (AML) is evolving, with promising therapies entering clinical translation, yet patient responses remain heterogeneous, and biomarkers for tailoring treatment are lacking. To understand how disease heterogeneity links with therapy response, we determined the leukemia cell hierarchy makeup from bulk transcriptomes of more than 1,000 patients through deconvolution using single-cell reference profiles of leukemia stem, progenitor and mature cell types. Leukemia hierarchy composition was associated with functional, genomic and clinical properties and converged into four overall classes, spanning Primitive, Mature, GMP and Intermediate. Critically, variation in hierarchy composition along the Primitive versus GMP or Primitive versus Mature axes were associated with response to chemotherapy or drug sensitivity profiles of targeted therapies, respectively. A seven-gene biomarker derived from the Primitive versus Mature axis was associated with response to 105 investigational drugs. Cellular hierarchy composition constitutes a novel framework for understanding disease biology and advancing precision medicine in AML.

Identification of the global miR-130a targetome reveals a role for TBL1XR1 in hematopoietic stem cell self-renewal and t(8;21) AML.

Gene expression profiling and proteome analysis of normal and malignant hematopoietic stem cells (HSCs) point to shared core stemness properties. However, discordance between mRNA and protein signatures highlights an important role for post-transcriptional regulation by microRNAs (miRNAs) in governing this critical nexus. Here, we identify miR-130a as a regulator of HSC self-renewal and differentiation. Enforced expression of miR-130a impairs B lymphoid differentiation and expands long-term HSCs. Integration of protein mass spectrometry and chimeric AGO2 crosslinking and immunoprecipitation (CLIP) identifies TBL1XR1 as a primary miR-130a target, whose loss of function phenocopies miR-130a overexpression. Moreover, we report that miR-130a is highly expressed in t(8;21) acute myeloid leukemia (AML), where it is critical for maintaining the oncogenic molecular program mediated by the AML1-ETO complex. Our study establishes that identification of the comprehensive miRNA targetome within primary cells enables discovery of genes and molecular networks underpinning stemness properties of normal and leukemic cells.

TFEB-mediated endolysosomal activity controls human hematopoietic stem cell fate.

It is critical to understand how human quiescent long-term hematopoietic stem cells (LT-HSCs) sense demand from daily and stress-mediated cues and then transition into bioenergetically active progeny to differentiate and meet these cellular needs. However, the demand-adapted regulatory circuits of these early steps of hematopoiesis are largely unknown. Here we show that lysosomes, sophisticated nutrient-sensing and signaling centers, are regulated dichotomously by transcription factor EB (TFEB) and MYC to balance catabolic and anabolic processes required for activating LT-HSCs and guiding their lineage fate. TFEB-mediated induction of the endolysosomal pathway causes membrane receptor degradation, limiting LT-HSC metabolic and mitogenic activation, promoting quiescence and self-renewal, and governing erythroid-myeloid commitment. In contrast, MYC engages biosynthetic processes while repressing lysosomal catabolism, driving LT-HSC activation. Our study identifies TFEB-mediated control of lysosomal activity as a central regulatory hub for proper and coordinated stem cell fate determination.

Enhancer Hijacking Drives Oncogenic BCL11B Expression in Lineage-Ambiguous Stem Cell Leukemia.

Lineage-ambiguous leukemias are high-risk malignancies of poorly understood genetic basis. Here, we describe a distinct subgroup of acute leukemia with expression of myeloid, T lymphoid, and stem cell markers driven by aberrant allele-specific deregulation of BCL11B, a master transcription factor responsible for thymic T-lineage commitment and specification. Mechanistically, this deregulation was driven by chromosomal rearrangements that juxtapose BCL11B to superenhancers active in hematopoietic progenitors, or focal amplifications that generate a superenhancer from a noncoding element distal to BCL11B. Chromatin conformation analyses demonstrated long-range interactions of rearranged enhancers with the expressed BCL11B allele and association of BCL11B with activated hematopoietic progenitor cell cis-regulatory elements, suggesting BCL11B is aberrantly co-opted into a gene regulatory network that drives transformation by maintaining a progenitor state. These data support a role for ectopic BCL11B expression in primitive hematopoietic cells mediated by enhancer hijacking as an oncogenic driver of human lineage-ambiguous leukemia. SIGNIFICANCE: Lineage-ambiguous leukemias pose significant diagnostic and therapeutic challenges due to a poorly understood molecular and cellular basis. We identify oncogenic deregulation of BCL11B driven by diverse structural alterations, including de novo superenhancer generation, as the driving feature of a subset of lineage-ambiguous leukemias that transcend current diagnostic boundaries.This article is highlighted in the In This Issue feature, p. 2659.

A latent subset of human hematopoietic stem cells resists regenerative stress to preserve stemness.

Continuous supply of immune cells throughout life relies on the delicate balance in the hematopoietic stem cell (HSC) pool between long-term maintenance and meeting the demands of both normal blood production and unexpected stress conditions. Here we identified distinct subsets of human long-term (LT)-HSCs that responded differently to regeneration-mediated stress: an immune checkpoint ligand CD112lo subset that exhibited a transient engraftment restraint (termed latency) before contributing to hematopoietic reconstitution and a primed CD112hi subset that responded rapidly. This functional heterogeneity and CD112 expression are regulated by INKA1 through direct interaction with PAK4 and SIRT1, inducing epigenetic changes and defining an alternative state of LT-HSC quiescence that serves to preserve self-renewal and regenerative capacity upon regeneration-mediated stress. Collectively, our data uncovered the molecular intricacies underlying HSC heterogeneity and self-renewal regulation and point to latency as an orchestrated physiological response that balances blood cell demands with preserving a stem cell reservoir.

Sphingosine-1-phosphate receptor 3 potentiates inflammatory programs in normal and leukemia stem cells to promote differentiation.

Acute myeloid leukemia (AML) is a caricature of normal hematopoiesis, driven from leukemia stem cells (LSC) that share some hematopoietic stem cell (HSC) programs including responsiveness to inflammatory signaling. Although inflammation dysregulates mature myeloid cells and influences stemness programs and lineage determination in HSC by activating stress myelopoiesis, such roles in LSC are poorly understood. Here, we show that S1PR3, a receptor for the bioactive lipid sphingosine-1-phosphate, is a central regulator which drives myeloid differentiation and activates inflammatory programs in both HSC and LSC. S1PR3-mediated inflammatory signatures varied in a continuum from primitive to mature myeloid states across AML patient cohorts, each with distinct phenotypic and clinical properties. S1PR3 was high in LSC and blasts of mature myeloid samples with linkages to chemosensitivity, while S1PR3 activation in primitive samples promoted LSC differentiation leading to eradication. Our studies open new avenues for therapeutic target identification specific for each AML subset.

Integration of intra-sample contextual error modeling for improved detection of somatic mutations from deep sequencing.

Sensitive mutation detection by next-generation sequencing is critical for early cancer detection, monitoring minimal/measurable residual disease (MRD), and guiding precision oncology. Nevertheless, because of artifacts introduced during library preparation and sequencing, the detection of low-frequency variants at high specificity is problematic. Here, we present Espresso, an error suppression method that considers local sequence features to accurately detect single-nucleotide variants (SNVs). Compared to other advanced error suppression techniques, Espresso consistently demonstrated lower numbers of false-positive mutation calls and greater sensitivity. We demonstrated Espresso's superior performance in detecting MRD in the peripheral blood of patients with acute myeloid leukemia (AML) throughout their treatment course. Furthermore, we showed that accurate mutation calling in a small number of informative genomic loci might provide a cost-efficient strategy for pragmatic risk prediction of AML development in healthy individuals. More broadly, we aim for Espresso to aid with accurate mutation detection in many other research and clinical settings.

Newsletter Fall 2020: Clinician Investigator Trainee Association of Canada (CITAC).

Message from the CITAC president To say that 2020 has been an unprecedented year is an understatement. The coronavirus disease 2019 (COVID-19) global pandemic and the major societal awakening on racial equity and justice have led us to reflect on our direction, goals and mission. Thanks to our talented and dedicated executive team, we were able to pivot our efforts and adapt to the changing landscape of research and advocacy. In April, we provided our members with a list of resources to help facilitate a smooth transition to working from home. In June, we published Clinician Investigator Trainee Association of Canada's (CITAC) press release on our role in combating anti-Black discrimination and racial injustice and have outlined specific advocacy efforts that we will be committing to over the next years (the full statement can be found on our website, https://www.citac-accfc.org).

Fundamental trends within falling match rates: Insights from the past decade of Canadian residency matching data.

BACKGROUND: The number of unmatched Canadian Medical Graduates (CMGs) has risen dramatically over the last decade. To identify long-term solutions to this problem, an understanding of the factors contributing to these rising unmatched rates is critical. METHODS: Using match and electives data from 2009-2019, we employed machine learning algorithms to identify three clusters of disciplines with distinct trends in match and electives behaviours. We assessed the relationships between unmatched rates, competitiveness, rates of parallel planning, and program selection practices at a discipline level. RESULTS: Across Canada, growth in CMGs has outpaced growth in residency seats, narrowing the seat-to-applicant ratio. Yet not all disciplines have been affected equally: a subset of surgical disciplines experienced a consistent decline in residency seats over time. Applicants to these disciplines are also at disproportionate risk of becoming unmatched, and this is associated with lower rates of parallel planning as quantified through clinical electives and match applications. This, in turn, is associated with the program selection practices of these disciplines. CONCLUSION: Long-term solutions to the unmatched CMG crisis require more nuance than indiscriminately increasing residency seats and should consider cluster specific match ratios as well as regulations around clinical electives and program selection practices.

CONTEXTE: Le nombre de diplômés canadiens en médecine (DCM) non jumelés a augmenté considérablement au cours des dix dernières années. Afin de trouver des solutions à long terme à ce problème, il est primordial de comprendre les facteurs qui contribuent à cette hausse. MÉTHODES: À l'aide des données de jumelages et de stages à option de 2009 à 2019, nous nous sommes servis d'algorithmes d'apprentissage automatique afin d'identifier trois groupes de disciplines démontrant des tendances distinctes en ce qui a trait aux jumelages et au choix des stages à option. Nous avons évalué les relations entre le taux de diplômés non jumelés, la compétitivité, les taux de planification parallèle et les pratiques de sélection des programmes pour chacune de ses disciplines. RÉSULTATS: Partout au Canada, la croissance des DCM a dépassé la croissance du nombre de postes de résidence, réduisant ainsi le ratio postes-candidats. Cependant, les disciplines n'ont pas toutes été touchées de la même manière: un sous-ensemble de disciplines en chirurgie a connu, au fil du temps, un déclin continu en ce qui a trait aux postes de résidence offerts. Les candidats de ces disciplines sont aussi exposés à un risque démesuré de ne pas être jumelés et ceci est lié à la réduction des taux de planification parallèle tels que quantifiés par les stages à option cliniques et les demandes de jumelage. Ceci est, par conséquent, lié aux pratiques de sélection des programmes de ces disciplines. CONCLUSION: Les solutions à long terme de la crise touchant les DCM non jumelés requièrent plus de subtilités que le simple fait d'augmenter sans distinction le nombre de postes de résidence. Elles devraient également prendre en compte les ratios de jumelage propres aux groupes de disciplines ainsi que les règlements concernant les stages à option et les pratiques de sélection des programmes.

Sphingolipid Modulation Activates Proteostasis Programs to Govern Human Hematopoietic Stem Cell Self-Renewal.

Cellular stress responses serve as crucial decision points balancing persistence or culling of hematopoietic stem cells (HSCs) for lifelong blood production. Although strong stressors cull HSCs, the linkage between stress programs and self-renewal properties that underlie human HSC maintenance remains unknown, particularly at quiescence exit when HSCs must also dynamically shift metabolic state. Here, we demonstrate distinct wiring of the sphingolipidome across the human hematopoietic hierarchy and find that genetic or pharmacologic modulation of the sphingolipid enzyme DEGS1 regulates lineage differentiation. Inhibition of DEGS1 in hematopoietic stem and progenitor cells during the transition from quiescence to cellular activation with N-(4-hydroxyphenyl) retinamide activates coordinated stress pathways that coalesce on endoplasmic reticulum stress and autophagy programs to maintain immunophenotypic and functional HSCs. Thus, our work identifies a linkage between sphingolipid metabolism, proteostatic quality control systems, and HSC self-renewal and provides therapeutic targets for improving HSC-based cellular therapeutics.

Somatic Mitochondrial DNA Mutations in Diffuse Large B-Cell Lymphoma.

Diffuse Large B-Cell Lymphoma (DLBCL) is an aggressive hematological cancer for which mitochondrial metabolism may play an important role. Mitochondrial DNA (mtDNA) encodes crucial mitochondrial proteins, yet the relationship between mtDNA and DLBCL remains unclear. We analyzed the functional consequences and mutational spectra of mtDNA somatic mutations and private constitutional variants in 40 DLBCL tumour-normal pairs. While private constitutional variants occurred frequently in the D-Loop, somatic mutations were randomly distributed across the mitochondrial genome. Heteroplasmic constitutional variants showed a trend towards loss of heteroplasmy in the corresponding tumour regardless of whether the reference or variant allele was being lost, suggesting that these variants are selectively neutral. The mtDNA mutational spectrum showed minimal support for ROS damage and revealed strand asymmetry with increased C > T and A > G transitions on the heavy strand, consistent with a replication-associated mode of mutagenesis. These heavy strand transitions carried higher proportions of amino acid changes - which were also more pathogenic - than equivalent substitutions on the light strand. Taken together, endogenous replication-associated events underlie mtDNA mutagenesis in DLBCL and preferentially generate functionally consequential mutations. Yet mtDNA somatic mutations remain selectively neutral, suggesting that mtDNA-encoded mitochondrial functions may not play an important role in DLBCL.