Monitoring clonal burden as an alternative to blast count for myelodysplastic neoplasm treatment response.

Accurate assessment of therapy response in myelodysplastic neoplasm (MDS) has been challenging. Directly monitoring mutational disease burden may be useful, but is not currently included in MDS response criteria, and the correlation of mutational burden and traditional response endpoints is not completely understood. Here, we used genome-wide and targeted next-generation sequencing (NGS) to monitor clonal and subclonal molecular disease burden in 452 samples from 32 patients prospectively treated in a clinical trial. Molecular responses were compared with International Working Group (IWG) 2006-defined response assessments. We found that myeloblast percentage consistently underestimates MDS molecular disease burden and that mutational clearance patterns for marrow complete remission (mCR), which depends on myeloblast assessment, was not different than stable disease or bone marrow aplasia, underscoring a major limitation of using mCR. In contrast, achieving a complete remission (CR) was associated with the highest level of mutation clearance and lowest residual mutational burden in higher-risk MDS patients. A targeted gene panel approach was inferior to genome-wide sequencing in defining subclones and their molecular responses but may be adequate for monitoring molecular disease burden when a targeted gene is present in the founding clone. Our work supports incorporating serial NGS-based monitoring into prospective MDS clinical trials.

Differential chromatin accessibility and transcriptional dynamics define breast cancer subtypes and their lineages.

Breast cancer (BC) is defined by distinct molecular subtypes with different cells of origin. The transcriptional networks that characterize the subtype-specific tumor-normal lineages are not established. In this work, we applied bulk, single-cell and single-nucleus multi-omic techniques as well as spatial transcriptomics and multiplex imaging on 61 samples from 37 patients with BC to show characteristic links in gene expression and chromatin accessibility between BC subtypes and their putative cells of origin. Regulatory network analysis of transcription factors underscored the importance of BHLHE40 in luminal BC and luminal mature cells and KLF5 in basal-like tumors and luminal progenitor cells. Furthermore, we identify key genes defining the basal-like (SOX6 and KCNQ3) and luminal A/B (FAM155A and LRP1B) lineages. Exhausted CTLA4-expressing CD8+ T cells were enriched in basal-like BC, suggesting an altered means of immune dysfunction. These findings demonstrate analysis of paired transcription and chromatin accessibility at the single-cell level is a powerful tool for investigating cancer lineage and highlight transcriptional networks that define basal and luminal BC lineages.

Somatic mutation phasing and haplotype extension using linked-reads in multiple myeloma.

Somatic mutation phasing informs our understanding of cancer-related events, like driver mutations. We generated linked-read whole genome sequencing data for 23 samples across disease stages from 14 multiple myeloma (MM) patients and systematically assigned somatic mutations to haplotypes using linked-reads. Here, we report the reconstructed cancer haplotypes and phase blocks from several MM samples and show how phase block length can be extended by integrating samples from the same individual. We also uncover phasing information in genes frequently mutated in MM, including DIS3, HIST1H1E, KRAS, NRAS, and TP53, phasing 79.4% of 20,705 high-confidence somatic mutations. In some cases, this enabled us to interpret clonal evolution models at higher resolution using pairs of phased somatic mutations. For example, our analysis of one patient suggested that two NRAS hotspot mutations occurred on the same haplotype but were independent events in different subclones. Given sufficient tumor purity and data quality, our framework illustrates how haplotype-aware analysis of somatic mutations in cancer can be beneficial for some cancer cases.

Neoantigen landscape supports feasibility of personalized cancer vaccine for follicular lymphoma.

ABSTRACT: Personalized cancer vaccines designed to target neoantigens represent a promising new treatment paradigm in oncology. In contrast to classical idiotype vaccines, we hypothesized that "polyvalent" vaccines could be engineered for the personalized treatment of follicular lymphoma (FL) using neoantigen discovery by combined whole-exome sequencing (WES) and RNA sequencing (RNA-seq). Fifty-eight tumor samples from 57 patients with FL underwent WES and RNA-seq. Somatic and B-cell clonotype neoantigens were predicted and filtered to identify high-quality neoantigens. B-cell clonality was determined by the alignment of B-cell receptor (BCR) CDR3 regions from RNA-seq data, grouping at the protein level, and comparison with the BCR repertoire from healthy individuals using RNA-seq data. An average of 52 somatic mutations per patient (range, 2-172) were identified, and ≥2 (median, 15) high-quality neoantigens were predicted for 56 of 58 FL samples. The predicted neoantigen peptides were composed of missense mutations (77%), indels (9%), gene fusions (3%), and BCR sequences (11%). Building off of these preclinical analyses, we initiated a pilot clinical trial using personalized neoantigen vaccination combined with PD-1 blockade in patients with relapsed or refractory FL (#NCT03121677). Synthetic long peptide vaccines targeting predicted high-quality neoantigens were successfully synthesized for and administered to all 4 patients enrolled. Initial results demonstrate feasibility, safety, and potential immunologic and clinical responses. Our study suggests that a genomics-driven personalized cancer vaccine strategy is feasible for patients with FL, and this may overcome prior challenges in the field. This trial was registered at www.ClinicalTrials.gov as #NCT03121677.

Ultra-Deep Sequencing Reveals the Mutational Landscape of Classical Hodgkin Lymphoma.

The malignant Hodgkin and Reed Sternberg (HRS) cells of classical Hodgkin lymphoma (cHL) are scarce in affected lymph nodes, creating a challenge to detect driver somatic mutations. As an alternative to cell purification techniques, we hypothesized that ultra-deep exome sequencing would allow genomic study of HRS cells, thereby streamlining analysis and avoiding technical pitfalls. To test this, 31 cHL tumor/normal pairs were exome sequenced to approximately 1,000× median depth of coverage. An orthogonal error-corrected sequencing approach verified >95% of the discovered mutations. We identified mutations in genes novel to cHL including: CDH5 and PCDH7, novel stop gain mutations in IL4R, and a novel pattern of recurrent mutations in pathways regulating Hippo signaling. As a further application of our exome sequencing, we attempted to identify expressed somatic single-nucleotide variants (SNV) in single-nuclei RNA sequencing (snRNA-seq) data generated from a patient in our cohort. Our snRNA analysis identified a clear cluster of cells containing a somatic SNV identified in our deep exome data. This cluster has differentially expressed genes that are consistent with genes known to be dysregulated in HRS cells (e.g., PIM1 and PIM3). The cluster also contains cells with an expanded B-cell clonotype further supporting a malignant phenotype. This study provides proof-of-principle that ultra-deep exome sequencing can be utilized to identify recurrent mutations in HRS cells and demonstrates the feasibility of snRNA-seq in the context of cHL. These studies provide the foundation for the further analysis of genomic variants in large cohorts of patients with cHL. SIGNIFICANCE: Our data demonstrate the utility of ultra-deep exome sequencing in uncovering somatic variants in Hodgkin lymphoma, creating new opportunities to define the genes that are recurrently mutated in this disease. We also show for the first time the successful application of snRNA-seq in Hodgkin lymphoma and describe the expression profile of a putative cluster of HRS cells in a single patient.

Differential chromatin accessibility and transcriptional dynamics define breast cancer subtypes and their lineages.

Breast cancer is a heterogeneous disease, and treatment is guided by biomarker profiles representing distinct molecular subtypes. Breast cancer arises from the breast ductal epithelium, and experimental data suggests breast cancer subtypes have different cells of origin within that lineage. The precise cells of origin for each subtype and the transcriptional networks that characterize these tumor-normal lineages are not established. In this work, we applied bulk, single-cell (sc), and single-nucleus (sn) multi-omic techniques as well as spatial transcriptomics and multiplex imaging on 61 samples from 37 breast cancer patients to show characteristic links in gene expression and chromatin accessibility between breast cancer subtypes and their putative cells of origin. We applied the PAM50 subtyping algorithm in tandem with bulk RNA-seq and snRNA-seq to reliably subtype even low-purity tumor samples and confirm promoter accessibility using snATAC. Trajectory analysis of chromatin accessibility and differentially accessible motifs clearly connected progenitor populations with breast cancer subtypes supporting the cell of origin for basal-like and luminal A and B tumors. Regulatory network analysis of transcription factors underscored the importance of BHLHE40 in luminal breast cancer and luminal mature cells, and KLF5 in basal-like tumors and luminal progenitor cells. Furthermore, we identify key genes defining the basal-like ( PRKCA , SOX6 , RGS6 , KCNQ3 ) and luminal A/B ( FAM155A , LRP1B ) lineages, with expression in both precursor and cancer cells and further upregulation in tumors. Exhausted CTLA4-expressing CD8+ T cells were enriched in basal-like breast cancer, suggesting altered means of immune dysfunction among breast cancer subtypes. We used spatial transcriptomics and multiplex imaging to provide spatial detail for key markers of benign and malignant cell types and immune cell colocation. These findings demonstrate analysis of paired transcription and chromatin accessibility at the single cell level is a powerful tool for investigating breast cancer lineage development and highlight transcriptional networks that define basal and luminal breast cancer lineages.

Genomic landscape of TP53-mutated myeloid malignancies.

TP53-mutated myeloid malignancies are associated with complex cytogenetics and extensive structural variants, which complicates detailed genomic analysis by conventional clinical techniques. We performed whole-genome sequencing (WGS) of 42 acute myeloid leukemia (AML)/myelodysplastic syndromes (MDS) cases with paired normal tissue to better characterize the genomic landscape of TP53-mutated AML/MDS. WGS accurately determines TP53 allele status, a key prognostic factor, resulting in the reclassification of 12% of cases from monoallelic to multihit. Although aneuploidy and chromothripsis are shared with most TP53-mutated cancers, the specific chromosome abnormalities are distinct to each cancer type, suggesting a dependence on the tissue of origin. ETV6 expression is reduced in nearly all cases of TP53-mutated AML/MDS, either through gene deletion or presumed epigenetic silencing. Within the AML cohort, mutations of NF1 are highly enriched, with deletions of 1 copy of NF1 present in 45% of cases and biallelic mutations in 17%. Telomere content is increased in TP53-mutated AMLs compared with other AML subtypes, and abnormal telomeric sequences were detected in the interstitial regions of chromosomes. These data highlight the unique features of TP53-mutated myeloid malignancies, including the high frequency of chromothripsis and structural variation, the frequent involvement of unique genes (including NF1 and ETV6) as cooperating events, and evidence for altered telomere maintenance.

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.

Persistent Molecular Disease in Adult Patients With AML Evaluated With Whole-Exome and Targeted Error-Corrected DNA Sequencing.

PURPOSE: Persistent molecular disease (PMD) after induction chemotherapy predicts relapse in AML. In this study, we used whole-exome sequencing (WES) and targeted error-corrected sequencing to assess the frequency and mutational patterns of PMD in 30 patients with AML. MATERIALS AND METHODS: The study cohort included 30 patients with adult AML younger than 65 years who were uniformly treated with standard induction chemotherapy. Tumor/normal WES was performed for all patients at presentation. PMD analysis was evaluated in bone marrow samples obtained during clinicopathologic remission using repeat WES and analysis of patient-specific mutations and error-corrected sequencing of 40 recurrently mutated AML genes (MyeloSeq). RESULTS: WES for patient-specific mutations detected PMD in 63% of patients (19/30) using a minimum variant allele fraction (VAF) of 2.5%. In comparison, MyeloSeq identified persistent mutations above 0.1% VAF in 77% of patients (23/30). PMD was usually present at relatively high levels (>2.5% VAFs), such that WES and MyeloSeq agreed for 73% of patients despite differences in detection limits. Mutations in DNMT3A, ASXL1, and TET2 (ie, DTA mutations) were persistent in 16 of 17 patients, but WES also detected non-DTA mutations in 14 of these patients, which for some patients distinguished residual AML cells from clonal hematopoiesis. Surprisingly, MyeloSeq detected additional variants not identified at presentation in 73% of patients that were consistent with new clonal cell populations after chemotherapy. CONCLUSION: PMD and clonal hematopoiesis are both common in patients with AML in first remission. These findings demonstrate the importance of baseline testing for accurate interpretation of mutation-based tumor monitoring assays for patients with AML and highlight the need for clinical trials to determine whether these complex mutation patterns correlate with clinical outcomes in AML.

Utility of targeted gene sequencing to differentiate myeloid malignancies from other cytopenic conditions.

The National Heart, Lung, and Blood Institute-funded National MDS Natural History Study (NCT02775383) is a prospective cohort study enrolling patients with cytopenia with suspected myelodysplastic syndromes (MDS) to evaluate factors associated with disease. Here, we sequenced 53 genes in bone marrow samples harvested from 1298 patients diagnosed with myeloid malignancy, including MDS and non-MDS myeloid malignancy or alternative marrow conditions with cytopenia based on concordance between independent histopathologic reviews (local, centralized, and tertiary to adjudicate disagreements when needed). We developed a novel 2-stage diagnostic classifier based on mutational profiles in 18 of 53 sequenced genes that were sufficient to best predict a diagnosis of myeloid malignancy and among those with a predicted myeloid malignancy, predict whether they had MDS. The classifier achieved a positive predictive value (PPV) of 0.84 and negative predictive value (NPV) of 0.8 with an area under the receiver operating characteristic curve (AUROC) of 0.85 when classifying patients as having myeloid vs no myeloid malignancy based on variant allele frequencies (VAFs) in 17 genes and a PPV of 0.71 and NPV of 0.64 with an AUROC of 0.73 when classifying patients as having MDS vs non-MDS malignancy based on VAFs in 10 genes. We next assessed how this approach could complement histopathology to improve diagnostic accuracy. For 99 of 139 (71%) patients (PPV of 0.83 and NPV of 0.65) with local and centralized histopathologic disagreement in myeloid vs no myeloid malignancy, the classifier-predicted diagnosis agreed with the tertiary pathology review (considered the internal gold standard).

Genomic landscape of TP53 -mutated myeloid malignancies.

TP53 -mutated myeloid malignancies are most frequently associated with complex cytogenetics. The presence of complex and extensive structural variants complicates detailed genomic analysis by conventional clinical techniques. We performed whole genome sequencing of 42 AML/MDS cases with paired normal tissue to characterize the genomic landscape of TP53 -mutated myeloid malignancies. The vast majority of cases had multi-hit involvement at the TP53 genetic locus (94%), as well as aneuploidy and chromothripsis. Chromosomal patterns of aneuploidy differed significantly from TP53 -mutated cancers arising in other tissues. Recurrent structural variants affected regions that include ETV6 on chr12p, RUNX1 on chr21, and NF1 on chr17q. Most notably for ETV6 , transcript expression was low in cases of TP53 -mutated myeloid malignancies both with and without structural rearrangements involving chromosome 12p. Telomeric content is increased in TP53 -mutated AML/MDS compared other AML subtypes, and telomeric content was detected adjacent to interstitial regions of chromosomes. The genomic landscape of TP53 -mutated myeloid malignancies reveals recurrent structural variants affecting key hematopoietic transcription factors and telomeric repeats that are generally not detected by panel sequencing or conventional cytogenetic analyses. Key Points: WGS comprehensively determines TP53 mutation status, resulting in the reclassification of 12% of cases from mono-allelic to multi-hit Chromothripsis is more frequent than previously appreciated, with a preference for specific chromosomes ETV6 is deleted in 45% of cases, with evidence for epigenetic suppression in non-deleted cases NF1 is mutated in 48% of cases, with multi-hit mutations in 17% of these cases TP53 -mutated AML/MDS is associated with altered telomere content compared with other AMLs.

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.

Nine out of ten samples were mistakenly switched by The Orang-utan Genome Consortium.

The Sumatran orang-utan (Pongo abelii) reference genome was first published in 2011, in conjunction with ten re-sequenced genomes from unrelated wild-caught individuals. Together, these published data have been utilized in almost all great ape genomic studies, plus in much broader comparative genomic research. Here, we report that the original sequencing Consortium inadvertently switched nine of the ten samples and/or resulting re-sequenced genomes, erroneously attributing eight of these to the wrong source individuals. Among them is a genome from the recently identified Tapanuli (P. tapanuliensis) species: thus, this genome was sequenced and published a full six years prior to the species' description. Sex was wrongly assigned to five known individuals; the numbers in one sample identifier were swapped; and the identifier for another sample most closely resembles that of a sample from another individual entirely. These errors have been reproduced in countless subsequent manuscripts, with noted implications for studies reliant on data from known individuals.

Spatially restricted drivers and transitional cell populations cooperate with the microenvironment in untreated and chemo-resistant pancreatic cancer.

Pancreatic ductal adenocarcinoma is a lethal disease with limited treatment options and poor survival. We studied 83 spatial samples from 31 patients (11 treatment-naïve and 20 treated) using single-cell/nucleus RNA sequencing, bulk-proteogenomics, spatial transcriptomics and cellular imaging. Subpopulations of tumor cells exhibited signatures of proliferation, KRAS signaling, cell stress and epithelial-to-mesenchymal transition. Mapping mutations and copy number events distinguished tumor populations from normal and transitional cells, including acinar-to-ductal metaplasia and pancreatic intraepithelial neoplasia. Pathology-assisted deconvolution of spatial transcriptomic data identified tumor and transitional subpopulations with distinct histological features. We showed coordinated expression of TIGIT in exhausted and regulatory T cells and Nectin in tumor cells. Chemo-resistant samples contain a threefold enrichment of inflammatory cancer-associated fibroblasts that upregulate metallothioneins. Our study reveals a deeper understanding of the intricate substructure of pancreatic ductal adenocarcinoma tumors that could help improve therapy for patients with this disease.

Convergent Clonal Evolution of Signaling Gene Mutations Is a Hallmark of Myelodysplastic Syndrome Progression.

Progression from myelodysplastic syndromes (MDS) to secondary acute myeloid leukemia (AML) is associated with the acquisition and expansion of subclones. Our understanding of subclone evolution during progression, including the frequency and preferred order of gene mutation acquisition, remains incomplete. Sequencing of 43 paired MDS and secondary AML samples identified at least one signaling gene mutation in 44% of MDS and 60% of secondary AML samples, often below the level of standard sequencing detection. In addition, 19% of MDS and 47% of secondary AML patients harbored more than one signaling gene mutation, almost always in separate, coexisting subclones. Signaling gene mutations demonstrated diverse patterns of clonal evolution during disease progression, including acquisition, expansion, persistence, and loss of mutations, with multiple patterns often coexisting in the same patient. Multivariate analysis revealed that MDS patients who had a signaling gene mutation had a higher risk of AML progression, potentially providing a biomarker for progression. SIGNIFICANCE: Subclone expansion is a hallmark of progression from MDS to secondary AML. Subclonal signaling gene mutations are common at MDS (often at low levels), show complex and convergent patterns of clonal evolution, and are associated with future progression to secondary AML. See related article by Guess et al., p. 316 (33). See related commentary by Romine and van Galen, p. 270. This article is highlighted in the In This Issue feature, p. 265.

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.

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.

The complete sequence of a human genome.

Since its initial release in 2000, the human reference genome has covered only the euchromatic fraction of the genome, leaving important heterochromatic regions unfinished. Addressing the remaining 8% of the genome, the Telomere-to-Telomere (T2T) Consortium presents a complete 3.055 billion-base pair sequence of a human genome, T2T-CHM13, that includes gapless assemblies for all chromosomes except Y, corrects errors in the prior references, and introduces nearly 200 million base pairs of sequence containing 1956 gene predictions, 99 of which are predicted to be protein coding. The completed regions include all centromeric satellite arrays, recent segmental duplications, and the short arms of all five acrocentric chromosomes, unlocking these complex regions of the genome to variational and functional studies.