High-Dimensional Analysis Delineates Myeloid and Lymphoid Compartment Remodeling during Successful Immune-Checkpoint Cancer Therapy.

Although current immune-checkpoint therapy (ICT) mainly targets lymphoid cells, it is associated with a broader remodeling of the tumor micro-environment. Here, using complementary forms of high-dimensional profiling, we define differences across all hematopoietic cells from syngeneic mouse tumors during unrestrained tumor growth or effective ICT. Unbiased assessment of gene expression of tumor-infiltrating cells by single-cell RNA sequencing (scRNAseq) and longitudinal assessment of cellular protein expression by mass cytometry (CyTOF) revealed significant remodeling of both the lymphoid and myeloid intratumoral compartments. Surprisingly, we observed multiple subpopulations of monocytes/macrophages, distinguishable by the markers CD206, CX3CR1, CD1d, and iNOS, that change over time during ICT in a manner partially dependent on IFNγ. Our data support the hypothesis that this macrophage polarization/activation results from effects on circulatory monocytes and early macrophages entering tumors, rather than on pre-polarized mature intratumoral macrophages.

CpG Island Hypermethylation Mediated by DNMT3A Is a Consequence of AML Progression.

DNMT3A mutations occur in ∼25% of acute myeloid leukemia (AML) patients. The most common mutation, DNMT3AR882H, has dominant negative activity that reduces DNA methylation activity by ∼80% in vitro. To understand the contribution of DNMT3A-dependent methylation to leukemogenesis, we performed whole-genome bisulfite sequencing of primary leukemic and non-leukemic cells in patients with or without DNMT3AR882 mutations. Non-leukemic hematopoietic cells with DNMT3AR882H displayed focal methylation loss, suggesting that hypomethylation antedates AML. Although virtually all AMLs with wild-type DNMT3A displayed CpG island hypermethylation, this change was not associated with gene silencing and was essentially absent in AMLs with DNMT3AR882 mutations. Primary hematopoietic stem cells expanded with cytokines were hypermethylated in a DNMT3A-dependent manner, suggesting that hypermethylation may be a response to, rather than a cause of, cellular proliferation. Our findings suggest that hypomethylation is an initiating phenotype in AMLs with DNMT3AR882, while DNMT3A-dependent CpG island hypermethylation is a consequence of AML progression.

Sequencing the mouse Y chromosome reveals convergent gene acquisition and amplification on both sex chromosomes.

We sequenced the MSY (male-specific region of the Y chromosome) of the C57BL/6J strain of the laboratory mouse Mus musculus. In contrast to theories that Y chromosomes are heterochromatic and gene poor, the mouse MSY is 99.9% euchromatic and contains about 700 protein-coding genes. Only 2% of the MSY derives from the ancestral autosomes that gave rise to the mammalian sex chromosomes. Instead, all but 45 of the MSY's genes belong to three acquired, massively amplified gene families that have no homologs on primate MSYs but do have acquired, amplified homologs on the mouse X chromosome. The complete mouse MSY sequence brings to light dramatic forces in sex chromosome evolution: lineage-specific convergent acquisition and amplification of X-Y gene families, possibly fueled by antagonism between acquired X-Y homologs. The mouse MSY sequence presents opportunities for experimental studies of a sex-specific chromosome in its entirety, in a genetically tractable model organism.

A draft human pangenome reference.

Here the Human Pangenome Reference Consortium presents a first draft of the human pangenome reference. The pangenome contains 47 phased, diploid assemblies from a cohort of genetically diverse individuals1. These assemblies cover more than 99% of the expected sequence in each genome and are more than 99% accurate at the structural and base pair levels. Based on alignments of the assemblies, we generate a draft pangenome that captures known variants and haplotypes and reveals new alleles at structurally complex loci. We also add 119 million base pairs of euchromatic polymorphic sequences and 1,115 gene duplications relative to the existing reference GRCh38. Roughly 90 million of the additional base pairs are derived from structural variation. Using our draft pangenome to analyse short-read data reduced small variant discovery errors by 34% and increased the number of structural variants detected per haplotype by 104% compared with GRCh38-based workflows, which enabled the typing of the vast majority of structural variant alleles per sample.

The Human Pangenome Project: a global resource to map genomic diversity.

The human reference genome is the most widely used resource in human genetics and is due for a major update. Its current structure is a linear composite of merged haplotypes from more than 20 people, with a single individual comprising most of the sequence. It contains biases and errors within a framework that does not represent global human genomic variation. A high-quality reference with global representation of common variants, including single-nucleotide variants, structural variants and functional elements, is needed. The Human Pangenome Reference Consortium aims to create a more sophisticated and complete human reference genome with a graph-based, telomere-to-telomere representation of global genomic diversity. Here we leverage innovations in technology, study design and global partnerships with the goal of constructing the highest-possible quality human pangenome reference. Our goal is to improve data representation and streamline analyses to enable routine assembly of complete diploid genomes. With attention to ethical frameworks, the human pangenome reference will contain a more accurate and diverse representation of global genomic variation, improve gene-disease association studies across populations, expand the scope of genomics research to the most repetitive and polymorphic regions of the genome, and serve as the ultimate genetic resource for future biomedical research and precision medicine.

Semi-automated assembly of high-quality diploid human reference genomes.

The current human reference genome, GRCh38, represents over 20 years of effort to generate a high-quality assembly, which has benefitted society1,2. However, it still has many gaps and errors, and does not represent a biological genome as it is a blend of multiple individuals3,4. Recently, a high-quality telomere-to-telomere reference, CHM13, was generated with the latest long-read technologies, but it was derived from a hydatidiform mole cell line with a nearly homozygous genome5. To address these limitations, the Human Pangenome Reference Consortium formed with the goal of creating high-quality, cost-effective, diploid genome assemblies for a pangenome reference that represents human genetic diversity6. Here, in our first scientific report, we determined which combination of current genome sequencing and assembly approaches yield the most complete and accurate diploid genome assembly with minimal manual curation. Approaches that used highly accurate long reads and parent-child data with graph-based haplotype phasing during assembly outperformed those that did not. Developing a combination of the top-performing methods, we generated our first high-quality diploid reference assembly, containing only approximately four gaps per chromosome on average, with most chromosomes within ±1% of the length of CHM13. Nearly 48% of protein-coding genes have non-synonymous amino acid changes between haplotypes, and centromeric regions showed the highest diversity. Our findings serve as a foundation for assembling near-complete diploid human genomes at scale for a pangenome reference to capture global genetic variation from single nucleotides to structural rearrangements.

The origin and evolution of mutations in acute myeloid leukemia.

Most mutations in cancer genomes are thought to be acquired after the initiating event, which may cause genomic instability and drive clonal evolution. However, for acute myeloid leukemia (AML), normal karyotypes are common, and genomic instability is unusual. To better understand clonal evolution in AML, we sequenced the genomes of M3-AML samples with a known initiating event (PML-RARA) versus the genomes of normal karyotype M1-AML samples and the exomes of hematopoietic stem/progenitor cells (HSPCs) from healthy people. Collectively, the data suggest that most of the mutations found in AML genomes are actually random events that occurred in HSPCs before they acquired the initiating mutation; the mutational history of that cell is "captured" as the clone expands. In many cases, only one or two additional, cooperating mutations are needed to generate the malignant founding clone. Cells from the founding clone can acquire additional cooperating mutations, yielding subclones that can contribute to disease progression and/or relapse.

Author Correction: Exome sequencing of Finnish isolates enhances rare-variant association power.

An Amendment to this paper has been published and can be accessed via a link at the top of the paper.

Exome sequencing of Finnish isolates enhances rare-variant association power.

Exome-sequencing studies have generally been underpowered to identify deleterious alleles with a large effect on complex traits as such alleles are mostly rare. Because the population of northern and eastern Finland has expanded considerably and in isolation following a series of bottlenecks, individuals of these populations have numerous deleterious alleles at a relatively high frequency. Here, using exome sequencing of nearly 20,000 individuals from these regions, we investigate the role of rare coding variants in clinically relevant quantitative cardiometabolic traits. Exome-wide association studies for 64 quantitative traits identified 26 newly associated deleterious alleles. Of these 26 alleles, 19 are either unique to or more than 20 times more frequent in Finnish individuals than in other Europeans and show geographical clustering comparable to Mendelian disease mutations that are characteristic of the Finnish population. We estimate that sequencing studies of populations without this unique history would require hundreds of thousands to millions of participants to achieve comparable association power.

Genome Sequencing as an Alternative to Cytogenetic Analysis in Myeloid Cancers.

BACKGROUND: Genomic analysis is essential for risk stratification in patients with acute myeloid leukemia (AML) or myelodysplastic syndromes (MDS). Whole-genome sequencing is a potential replacement for conventional cytogenetic and sequencing approaches, but its accuracy, feasibility, and clinical utility have not been demonstrated. METHODS: We used a streamlined whole-genome sequencing approach to obtain genomic profiles for 263 patients with myeloid cancers, including 235 patients who had undergone successful cytogenetic analysis. We adapted sample preparation, sequencing, and analysis to detect mutations for risk stratification using existing European Leukemia Network (ELN) guidelines and to minimize turnaround time. We analyzed the performance of whole-genome sequencing by comparing our results with findings from cytogenetic analysis and targeted sequencing. RESULTS: Whole-genome sequencing detected all 40 recurrent translocations and 91 copy-number alterations that had been identified by cytogenetic analysis. In addition, we identified new clinically reportable genomic events in 40 of 235 patients (17.0%). Prospective sequencing of samples obtained from 117 consecutive patients was performed in a median of 5 days and provided new genetic information in 29 patients (24.8%), which changed the risk category for 19 patients (16.2%). Standard AML risk groups, as defined by sequencing results instead of cytogenetic analysis, correlated with clinical outcomes. Whole-genome sequencing was also used to stratify patients who had inconclusive results by cytogenetic analysis into risk groups in which clinical outcomes were measurably different. CONCLUSIONS: In our study, we found that whole-genome sequencing provided rapid and accurate genomic profiling in patients with AML or MDS. Such sequencing also provided a greater diagnostic yield than conventional cytogenetic analysis and more efficient risk stratification on the basis of standard risk categories. (Funded by the Siteman Cancer Research Fund and others.).

IL-1ß expression in bone marrow dendritic cells is induced by TLR2 agonists and regulates HSC function.

Hematopoietic stem/progenitor cells (HSPCs) reside in localized microenvironments, or niches, in the bone marrow that provide key signals regulating their activity. A fundamental property of hematopoiesis is the ability to respond to environmental cues such as inflammation. How these cues are transmitted to HSPCs within hematopoietic niches is not well established. Here, we show that perivascular bone marrow dendritic cells (DCs) express a high basal level of Toll-like receptor-1 (TLR1) and TLR2. Systemic treatment with a TLR1/2 agonist induces HSPC expansion and mobilization. It also induces marked alterations in the bone marrow microenvironment, including a decrease in osteoblast activity and sinusoidal endothelial cell numbers. TLR1/2 agonist treatment of mice in which Myd88 is deleted specifically in DCs using Zbtb46-Cre show that the TLR1/2-induced expansion of multipotent HPSCs, but not HSPC mobilization or alterations in the bone marrow microenvironment, is dependent on TLR1/2 signaling in DCs. Interleukin-1ß (IL-1ß) is constitutively expressed in both murine and human DCs and is further induced after TLR1/2 stimulation. Systemic TLR1/2 agonist treatment of Il1r1-/- mice show that TLR1/2-induced HSPC expansion is dependent on IL-1ß signaling. Single-cell RNA-sequencing of low-risk myelodysplastic syndrome bone marrow revealed that IL1B and TLR1 expression is increased in DCs. Collectively, these data suggest a model in which TLR1/2 stimulation of DCs induces secretion of IL-1ß and other inflammatory cytokines into the perivascular niche, which in turn, regulates multipotent HSPCs. Increased DC TLR1/2 signaling may contribute to altered HSPC function in myelodysplastic syndrome by increasing local IL-1ß expression.

Immunosuppression and outcomes in adult patients with de novo acute myeloid leukemia with normal karyotypes.

Acute myeloid leukemia (AML) patients rarely have long first remissions (LFRs; >5 y) after standard-of-care chemotherapy, unless classified as favorable risk at presentation. Identification of the mechanisms responsible for long vs. more typical, standard remissions may help to define prognostic determinants for chemotherapy responses. Using exome sequencing, RNA-sequencing, and functional immunologic studies, we characterized 28 normal karyotype (NK)-AML patients with >5 y first remissions after chemotherapy (LFRs) and compared them to a well-matched group of 31 NK-AML patients who relapsed within 2 y (standard first remissions [SFRs]). Our combined analyses indicated that genetic-risk profiling at presentation (as defined by European LeukemiaNet [ELN] 2017 criteria) was not sufficient to explain the outcomes of many SFR cases. Single-cell RNA-sequencing studies of 15 AML samples showed that SFR AML cells differentially expressed many genes associated with immune suppression. The bone marrow of SFR cases had significantly fewer CD4+ Th1 cells; these T cells expressed an exhaustion signature and were resistant to activation by T cell receptor stimulation in the presence of autologous AML cells. T cell activation could be restored by removing the AML cells or blocking the inhibitory major histocompatibility complex class II receptor, LAG3. Most LFR cases did not display these features, suggesting that their AML cells were not as immunosuppressive. These findings were confirmed and extended in an independent set of 50 AML cases representing all ELN 2017 risk groups. AML cell-mediated suppression of CD4+ T cell activation at presentation is strongly associated with unfavorable outcomes in AML patients treated with standard chemotherapy.

Remethylation of Dnmt3a-/- hematopoietic cells is associated with partial correction of gene dysregulation and reduced myeloid skewing.

Mutations in the DNA methyltransferase 3A (DNMT3A) gene are the most common cause of age-related clonal hematopoiesis (ARCH) in older individuals, and are among the most common initiating events for acute myeloid leukemia (AML). The most frequent DNMT3A mutation in AML patients (R882H) encodes a dominant-negative protein that reduces methyltransferase activity by ∼80% in cells with heterozygous mutations, causing a focal, canonical DNA hypomethylation phenotype; this phenotype is partially recapitulated in murine Dnmt3a-/- bone marrow cells. To determine whether the hypomethylation phenotype of Dnmt3a-/- hematopoietic cells is reversible, we developed an inducible transgene to restore expression of DNMT3A in transplanted bone marrow cells from Dnmt3a-/- mice. Partial remethylation was detected within 1 wk, but near-complete remethylation required 6 mo. Remethylation was accurate, dynamic, and highly ordered, suggesting that differentially methylated regions have unique properties that may be relevant for their functions. Importantly, 22 wk of DNMT3A addback partially corrected dysregulated gene expression, and mitigated the expansion of myeloid cells. These data show that restoring DNMT3A expression can alter the epigenetic "state" created by loss of Dnmt3a activity; this genetic proof-of-concept experiment suggests that this approach could be relevant for patients with ARCH or AML caused by loss-of-function DNMT3A mutations.

Immune Escape of Relapsed AML Cells after Allogeneic Transplantation.

BACKGROUND: As consolidation therapy for acute myeloid leukemia (AML), allogeneic hematopoietic stem-cell transplantation provides a benefit in part by means of an immune-mediated graft-versus-leukemia effect. We hypothesized that the immune-mediated selective pressure imposed by allogeneic transplantation may cause distinct patterns of tumor evolution in relapsed disease. METHODS: We performed enhanced exome sequencing on paired samples obtained at initial presentation with AML and at relapse from 15 patients who had a relapse after hematopoietic stem-cell transplantation (with transplants from an HLA-matched sibling, HLA-matched unrelated donor, or HLA-mismatched unrelated donor) and from 20 patients who had a relapse after chemotherapy. We performed RNA sequencing and flow cytometry on a subgroup of these samples and on additional samples for validation. RESULTS: On exome sequencing, the spectrum of gained and lost mutations observed with relapse after transplantation was similar to the spectrum observed with relapse after chemotherapy. Specifically, relapse after transplantation was not associated with the acquisition of previously unknown AML-specific mutations or structural variations in immune-related genes. In contrast, RNA sequencing of samples obtained at relapse after transplantation revealed dysregulation of pathways involved in adaptive and innate immunity, including down-regulation of major histocompatibility complex (MHC) class II genes ( HLA-DPA1, HLA-DPB1, HLA-DQB1, and HLA-DRB1) to levels that were 3 to 12 times lower than the levels seen in paired samples obtained at presentation. Flow cytometry and immunohistochemical analysis confirmed decreased expression of MHC class II at relapse in 17 of 34 patients who had a relapse after transplantation. Evidence suggested that interferon-γ treatment could rapidly reverse this phenotype in AML blasts in vitro. CONCLUSIONS: AML relapse after transplantation was not associated with the acquisition of relapse-specific mutations in immune-related genes. However, it was associated with dysregulation of pathways that may influence immune function, including down-regulation of MHC class II genes, which are involved in antigen presentation. These epigenetic changes may be reversible with appropriate therapy. (Funded by the National Cancer Institute and others.).

Mutation Clearance after Transplantation for Myelodysplastic Syndrome.

BACKGROUND: Allogeneic hematopoietic stem-cell transplantation is the only curative treatment for patients with myelodysplastic syndrome (MDS). The molecular predictors of disease progression after transplantation are unclear. METHODS: We sequenced bone marrow and skin samples from 90 adults with MDS who underwent allogeneic hematopoietic stem-cell transplantation after a myeloablative or reduced-intensity conditioning regimen. We detected mutations before transplantation using enhanced exome sequencing, and we evaluated mutation clearance by using error-corrected sequencing to genotype mutations in bone marrow samples obtained 30 days after transplantation. In this exploratory study, we evaluated the association of a mutation detected after transplantation with disease progression and survival. RESULTS: Sequencing identified at least one validated somatic mutation before transplantation in 86 of 90 patients (96%); 32 of these patients (37%) had at least one mutation with a maximum variant allele frequency of at least 0.5% (equivalent to 1 heterozygous mutant cell in 100 cells) 30 days after transplantation. Patients with disease progression had mutations with a higher maximum variant allele frequency at 30 days than those who did not (median maximum variant allele frequency, 0.9% vs. 0%; P<0.001). The presence of at least one mutation with a variant allele frequency of at least 0.5% at day 30 was associated with a higher risk of progression (53.1% vs. 13.0%; conditioning regimen-adjusted hazard ratio, 3.86; 95% confidence interval [CI], 1.96 to 7.62; P<0.001) and a lower 1-year rate of progression-free survival than the absence of such a mutation (31.3% vs. 59.3%; conditioning regimen-adjusted hazard ratio for progression or death, 2.22; 95% CI, 1.32 to 3.73; P=0.005). The rate of progression-free survival was lower among patients who had received a reduced-intensity conditioning regimen and had at least one persistent mutation with a variant allele frequency of at least 0.5% at day 30 than among patients with other combinations of conditioning regimen and mutation status (P≤0.001). Multivariate analysis confirmed that patients who had a mutation with a variant allele frequency of at least 0.5% detected at day 30 had a higher risk of progression (hazard ratio, 4.48; 95% CI, 2.21 to 9.08; P<0.001) and a lower 1-year rate of progression-free survival than those who did not (hazard ratio for progression or death, 2.39; 95% CI, 1.40 to 4.09; P=0.002). CONCLUSIONS: The risk of disease progression was higher among patients with MDS in whom persistent disease-associated mutations were detected in the bone marrow 30 days after transplantation than among those in whom these mutations were not detected. (Funded by the Leukemia and Lymphoma Society and others.).

Rapid and Extraction-Free Detection of SARS-CoV-2 from Saliva by Colorimetric Reverse-Transcription Loop-Mediated Isothermal Amplification.

BACKGROUND: Rapid, reliable, and widespread testing is required to curtail the ongoing COVID-19 pandemic. Current gold-standard nucleic acid tests are hampered by supply shortages in critical reagents including nasal swabs, RNA extraction kits, personal protective equipment, instrumentation, and labor. METHODS: To overcome these challenges, we developed a rapid colorimetric assay using reverse-transcription loop-mediated isothermal amplification (RT-LAMP) optimized on human saliva samples without an RNA purification step. We describe the optimization of saliva pretreatment protocols to enable analytically sensitive viral detection by RT-LAMP. We optimized the RT-LAMP reaction conditions and implemented high-throughput unbiased methods for assay interpretation. We tested whether saliva pretreatment could also enable viral detection by conventional reverse-transcription quantitative polymerase chain reaction (RT-qPCR). Finally, we validated these assays on clinical samples. RESULTS: The optimized saliva pretreatment protocol enabled analytically sensitive extraction-free detection of SARS-CoV-2 from saliva by colorimetric RT-LAMP or RT-qPCR. In simulated samples, the optimized RT-LAMP assay had a limit of detection of 59 (95% confidence interval: 44-104) particle copies per reaction. We highlighted the flexibility of LAMP assay implementation using 3 readouts: naked-eye colorimetry, spectrophotometry, and real-time fluorescence. In a set of 30 clinical saliva samples, colorimetric RT-LAMP and RT-qPCR assays performed directly on pretreated saliva samples without RNA extraction had accuracies greater than 90%. CONCLUSIONS: Rapid and extraction-free detection of SARS-CoV-2 from saliva by colorimetric RT-LAMP is a simple, sensitive, and cost-effective approach with broad potential to expand diagnostic testing for the virus causing COVID-19.

Role of TP53 mutations in the origin and evolution of therapy-related acute myeloid leukaemia.

Therapy-related acute myeloid leukaemia (t-AML) and therapy-related myelodysplastic syndrome (t-MDS) are well-recognized complications of cytotoxic chemotherapy and/or radiotherapy. There are several features that distinguish t-AML from de novo AML, including a higher incidence of TP53 mutations, abnormalities of chromosomes 5 or 7, complex cytogenetics and a reduced response to chemotherapy. However, it is not clear how prior exposure to cytotoxic therapy influences leukaemogenesis. In particular, the mechanism by which TP53 mutations are selectively enriched in t-AML/t-MDS is unknown. Here, by sequencing the genomes of 22 patients with t-AML, we show that the total number of somatic single-nucleotide variants and the percentage of chemotherapy-related transversions are similar in t-AML and de novo AML, indicating that previous chemotherapy does not induce genome-wide DNA damage. We identified four cases of t-AML/t-MDS in which the exact TP53 mutation found at diagnosis was also present at low frequencies (0.003-0.7%) in mobilized blood leukocytes or bone marrow 3-6 years before the development of t-AML/t-MDS, including two cases in which the relevant TP53 mutation was detected before any chemotherapy. Moreover, functional TP53 mutations were identified in small populations of peripheral blood cells of healthy chemotherapy-naive elderly individuals. Finally, in mouse bone marrow chimaeras containing both wild-type and Tp53(+/-) haematopoietic stem/progenitor cells (HSPCs), the Tp53(+/-) HSPCs preferentially expanded after exposure to chemotherapy. These data suggest that cytotoxic therapy does not directly induce TP53 mutations. Rather, they support a model in which rare HSPCs carrying age-related TP53 mutations are resistant to chemotherapy and expand preferentially after treatment. The early acquisition of TP53 mutations in the founding HSPC clone probably contributes to the frequent cytogenetic abnormalities and poor responses to chemotherapy that are typical of patients with t-AML/t-MDS.

Evaluation of GRCh38 and de novo haploid genome assemblies demonstrates the enduring quality of the reference assembly.

The human reference genome assembly plays a central role in nearly all aspects of today's basic and clinical research. GRCh38 is the first coordinate-changing assembly update since 2009; it reflects the resolution of roughly 1000 issues and encompasses modifications ranging from thousands of single base changes to megabase-scale path reorganizations, gap closures, and localization of previously orphaned sequences. We developed a new approach to sequence generation for targeted base updates and used data from new genome mapping technologies and single haplotype resources to identify and resolve larger assembly issues. For the first time, the reference assembly contains sequence-based representations for the centromeres. We also expanded the number of alternate loci to create a reference that provides a more robust representation of human population variation. We demonstrate that the updates render the reference an improved annotation substrate, alter read alignments in unchanged regions, and impact variant interpretation at clinically relevant loci. We additionally evaluated a collection of new de novo long-read haploid assemblies and conclude that although the new assemblies compare favorably to the reference with respect to continuity, error rate, and gene completeness, the reference still provides the best representation for complex genomic regions and coding sequences. We assert that the collected updates in GRCh38 make the newer assembly a more robust substrate for comprehensive analyses that will promote our understanding of human biology and advance our efforts to improve health.