Single-cell genomics details the maturation block in BCP-ALL and identifies therapeutic vulnerabilities in DUX4-r cases.

ABSTRACT: B-cell progenitor acute lymphoblastic leukemia (BCP-ALL) is the most common childhood malignancy and is driven by multiple genetic alterations that cause maturation arrest and accumulation of abnormal progenitor B cells. Current treatment protocols with chemotherapy have led to favorable outcomes but are associated with significant toxicity and risk of side effects, highlighting the necessity for highly effective, less toxic, targeted drugs, even in subtypes with a favorable outcome. Here, we used multimodal single-cell sequencing to delineate the transcriptional, epigenetic, and immunophenotypic characteristics of 23 childhood BCP-ALLs belonging to the BCR::ABL1+, ETV6::RUNX1+, high hyperdiploid, and recently discovered DUX4-rearranged (DUX4-r) subtypes. Projection of the ALL cells along the normal hematopoietic differentiation axis revealed a diversity in the maturation pattern between the different BCP-ALL subtypes. Although the BCR::ABL1+, ETV6::RUNX1+, and high hyperdiploidy cells mainly showed similarities to normal pro-B cells, DUX4-r ALL cells also displayed transcriptional signatures resembling mature B cells. Focusing on the DUX4-r subtype, we found that the blast population displayed not only multilineage priming toward nonhematopoietic cells, myeloid, and T-cell lineages, but also an activation of phosphatidylinositol 3-kinase (PI3K)/AKT signaling that sensitized the cells to PI3K inhibition in vivo. Given the multilineage priming of DUX4-r blasts with aberrant expression of myeloid marker CD371 (CLL-1), we generated chimeric antigen receptor T cells, which effectively eliminated DUX4-r ALL cells in vivo. These results provide a detailed characterization of BCP-ALL at the single-cell level and reveal therapeutic vulnerabilities in the DUX4-r subtype, with implications for the understanding of ALL biology and new therapeutic strategies.

SMSSVD: SubMatrix Selection Singular Value Decomposition.

Motivation: High throughput biomedical measurements normally capture multiple overlaid biologically relevant signals and often also signals representing different types of technical artefacts like e.g. batch effects. Signal identification and decomposition are accordingly main objectives in statistical biomedical modeling and data analysis. Existing methods, aimed at signal reconstruction and deconvolution, in general, are either supervised, contain parameters that need to be estimated or present other types of ad hoc features. We here introduce SubMatrix Selection Singular Value Decomposition (SMSSVD), a parameter-free unsupervised signal decomposition and dimension reduction method, designed to reduce noise, adaptively for each low-rank-signal in a given data matrix, and represent the signals in the data in a way that enable unbiased exploratory analysis and reconstruction of multiple overlaid signals, including identifying groups of variables that drive different signals. Results: The SMSSVD method produces a denoised signal decomposition from a given data matrix. It also guarantees orthogonality between signal components in a straightforward manner and it is designed to make automation possible. We illustrate SMSSVD by applying it to several real and synthetic datasets and compare its performance to golden standard methods like PCA (Principal Component Analysis) and SPC (Sparse Principal Components, using Lasso constraints). The SMSSVD is computationally efficient and despite being a parameter-free method, in general, outperforms existing statistical learning methods. Availability and implementation: A Julia implementation of SMSSVD is openly available on GitHub (https://github.com/rasmushenningsson/SubMatrixSelectionSVD.jl). Supplementary information: Supplementary data are available at Bioinformatics online.

Increasing Clinical Severity during a Dengue Virus Type 3 Cuban Epidemic: Deep Sequencing of Evolving Viral Populations.

UNLABELLED: During the dengue virus type 3 (DENV-3) epidemic that occurred in Havana in 2001 to 2002, severe disease was associated with the infection sequence DENV-1 followed by DENV-3 (DENV-1/DENV-3), while the sequence DENV-2/DENV-3 was associated with mild/asymptomatic infections. To determine the role of the virus in the increasing severity demonstrated during the epidemic, serum samples collected at different time points were studied. A total of 22 full-length sequences were obtained using a deep-sequencing approach. Bayesian phylogenetic analysis of consensus sequences revealed that two DENV-3 lineages were circulating in Havana at that time, both grouped within genotype III. The predominant lineage is closely related to Peruvian and Ecuadorian strains, while the minor lineage is related to Venezuelan strains. According to consensus sequences, relatively few nonsynonymous mutations were observed; only one was fixed during the epidemic at position 4380 in the NS2B gene. Intrahost genetic analysis indicated that a significant minor population was selected and became predominant toward the end of the epidemic. In conclusion, greater variability was detected during the epidemic's progression in terms of significant minority variants, particularly in the nonstructural genes. An increasing trend of genetic diversity toward the end of the epidemic was observed only for synonymous variant allele rates, with higher variability in secondary cases. Remarkably, significant intrahost genetic variation was demonstrated within the same patient during the course of secondary infection with DENV-1/DENV-3, including changes in the structural proteins premembrane (PrM) and envelope (E). Therefore, the dynamic of evolving viral populations in the context of heterotypic antibodies could be related to the increasing clinical severity observed during the epidemic. IMPORTANCE: Based on the evidence that DENV fitness is context dependent, our research has focused on the study of viral factors associated with intraepidemic increasing severity in a unique epidemiological setting. Here, we investigated the intrahost genetic diversity in acute human samples collected at different time points during the DENV-3 epidemic that occurred in Cuba in 2001 to 2002 using a deep-sequencing approach. We concluded that greater variability in significant minor populations occurred as the epidemic progressed, particularly in the nonstructural genes, with higher variability observed in secondary infection cases. Remarkably, for the first time significant intrahost genetic variation was demonstrated within the same patient during the course of secondary infection with DENV-1/DENV-3, including changes in structural proteins. These findings indicate that high-resolution approaches are needed to unravel molecular mechanisms involved in dengue pathogenesis.

Group Selection and Contribution of Minority Variants during Virus Adaptation Determines Virus Fitness and Phenotype.

Understanding how a pathogen colonizes and adapts to a new host environment is a primary aim in studying emerging infectious diseases. Adaptive mutations arise among the thousands of variants generated during RNA virus infection, and identifying these variants will shed light onto how changes in tropism and species jumps can occur. Here, we adapted Coxsackie virus B3 to a highly permissive and less permissive environment. Using deep sequencing and bioinformatics, we identified a multi-step adaptive process to adaptation involving residues in the receptor footprints that correlated with receptor availability and with increase in virus fitness in an environment-specific manner. We show that adaptation occurs by selection of a dominant mutation followed by group selection of minority variants that together, confer the fitness increase observed in the population, rather than selection of a single dominant genotype.

The complement receptor C3AR constitutes a novel therapeutic target in NPM1-mutated AML.

Mutated nucleophosmin 1 (NPM1) is the most common genetic alteration in acute myeloid leukemia (AML), found in ∼30% of cases. Although mutations in this gene are considered favorable according to current risk stratification guidelines, a large fraction of patients will experience relapse, demonstrating the urgent need for new treatment options. Therefore, we aimed to identify cell surface proteins specifically expressed on NPM1-mutated AML cells, allowing for potential targeting with antibody-based therapies. Herein, we report on an arrayed flow cytometry-based screen directed to 362 cell surface markers. In comparing the cell surface expression on NPM1-mutated AML cells with primitive (CD34+ CD38-) normal bone marrow cells, we identified the complement receptor C3AR as being specifically expressed in NPM1-mutated AML. By flow cytometry and single-cell RNA sequencing, we further show that normal hematopoietic stem and progenitor cells lack detectable C3AR gene and protein expression, making it particularly suitable as a target for antibody therapy. We also demonstrate that C3AR in combination with GPR56 distinguishes the leukemic stem cells (LSCs) in NPM1-mutated AML from the normal hematopoietic stem cells, defining the LSC population, as shown by transplantation into immunodeficient mice. Mechanistically, the stimulation of C3AR-expressing cells with C3a, the ligand of C3AR, leads to the activation of ERK1/2 and increased survival of AML cells, suggesting that this is an important signaling axis in this subtype of AML. Finally, we show that antibodies directed against C3AR efficiently elicit natural killer cell-mediated killing of primary AML cells ex vivo, highlighting C3AR as a candidate therapeutic target in NPM1-mutated AML.

Clonal competition within complex evolutionary hierarchies shapes AML over time.

Clonal heterogeneity and evolution has major implications for disease progression and relapse in acute myeloid leukemia (AML). To model clonal dynamics in vivo, we serially transplanted 23 AML cases to immunodeficient mice and followed clonal composition for up to 15 months by whole-exome sequencing of 84 xenografts across two generations. We demonstrate vast changes in clonality that both progress and reverse over time, and define five patterns of clonal dynamics: Monoclonal, Stable, Loss, Expansion and Burst. We also show that subclonal expansion in vivo correlates with a more adverse prognosis. Furthermore, clonal expansion enabled detection of very rare clones with AML driver mutations that were undetectable by sequencing at diagnosis, demonstrating that the vast majority of AML cases harbor multiple clones already at diagnosis. Finally, the rise and fall of related clones enabled deconstruction of the complex evolutionary hierarchies of the clones that compete to shape AML over time.

Attenuation of RNA viruses by redirecting their evolution in sequence space.

RNA viruses pose serious threats to human health. Their success relies on their capacity to generate genetic variability and, consequently, on their adaptive potential. We describe a strategy to attenuate RNA viruses by altering their evolutionary potential. We rationally altered the genomes of Coxsackie B3 and influenza A viruses to redirect their evolutionary trajectories towards detrimental regions in sequence space. Specifically, viral genomes were engineered to harbour more serine and leucine codons with nonsense mutation targets: codons that could generate Stop mutations after a single nucleotide substitution. Indeed, these viruses generated more Stop mutations both in vitro and in vivo, accompanied by significant losses in viral fitness. In vivo, the viruses were attenuated, generated high levels of neutralizing antibodies and protected against lethal challenge. Our study demonstrates that cornering viruses in 'risky' areas of sequence space may be implemented as a broad-spectrum vaccine strategy against RNA viruses.

Whole-Genome Sequencing Analysis from the Chikungunya Virus Caribbean Outbreak Reveals Novel Evolutionary Genomic Elements.

BACKGROUND: Chikungunya virus (CHIKV), an alphavirus and member of the Togaviridae family, is capable of causing severe febrile disease in humans. In December of 2013 the Asian Lineage of CHIKV spread from the Old World to the Americas, spreading rapidly throughout the New World. Given this new emergence in naïve populations we studied the viral genetic diversity present in infected individuals to understand how CHIKV may have evolved during this continuing outbreak. METHODOLOGY/PRINCIPLE FINDINGS: We used deep-sequencing technologies coupled with well-established bioinformatics pipelines to characterize the minority variants and diversity present in CHIKV infected individuals from Guadeloupe and Martinique, two islands in the center of the epidemic. We observed changes in the consensus sequence as well as a diverse range of minority variants present at various levels in the population. Furthermore, we found that overall diversity was dramatically reduced after single passages in cell lines. Finally, we constructed an infectious clone from this outbreak and identified a novel 3' untranslated region (UTR) structure, not previously found in nature, that led to increased replication in insect cells. CONCLUSIONS/SIGNIFICANCE: Here we preformed an intrahost quasispecies analysis of the new CHIKV outbreak in the Caribbean. We identified novel variants present in infected individuals, as well as a new 3'UTR structure, suggesting that CHIKV has rapidly evolved in a short period of time once it entered this naïve population. These studies highlight the need to continue viral diversity surveillance over time as this epidemic evolves in order to understand the evolutionary potential of CHIKV.

Identification of ETV6-RUNX1-like and DUX4-rearranged subtypes in paediatric B-cell precursor acute lymphoblastic leukaemia.

Fusion genes are potent driver mutations in cancer. In this study, we delineate the fusion gene landscape in a consecutive series of 195 paediatric B-cell precursor acute lymphoblastic leukaemia (BCP ALL). Using RNA sequencing, we find in-frame fusion genes in 127 (65%) cases, including 27 novel fusions. We describe a subtype characterized by recurrent IGH-DUX4 or ERG-DUX4 fusions, representing 4% of cases, leading to overexpression of DUX4 and frequently co-occurring with intragenic ERG deletions. Furthermore, we identify a subtype characterized by an ETV6-RUNX1-like gene-expression profile and coexisting ETV6 and IKZF1 alterations. Thus, this study provides a detailed overview of fusion genes in paediatric BCP ALL and adds new pathogenetic insights, which may improve risk stratification and provide therapeutic options for this disease.