Pangenome graph construction from genome alignments with Minigraph-Cactus

Pangenome references address biases of reference genomes by storing a representative set of diverse haplotypes and their alignment, usually as a graph. Alternate alleles determined by variant callers can be used to construct pangenome graphs, but advances in long-read sequencing are leading to widely available, high-quality phased assemblies. Constructing a pangenome graph directly from assemblies, as opposed to variant calls, leverages the graph's ability to represent variation at different scales. Here we present the Minigraph-Cactus pangenome pipeline, which creates pangenomes directly from whole-genome alignments, and demonstrate its ability to scale to 90 human haplotypes from the Human Pangenome Reference Consortium. The method builds graphs containing all forms of genetic variation while allowing use of current mapping and genotyping tools. We measure the effect of the quality and completeness of reference genomes used for analysis within the pangenomes and show that using the CHM13 reference from the Telomere-to-Telomere Consortium improves the accuracy of our methods. We also demonstrate construction of a Drosophila melanogaster pangenome.

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.

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).

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.

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.

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.

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.

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.

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.

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.

Failure to Detect Mutations in U2AF1 due to Changes in the GRCh38 Reference Sequence

The U2AF1 gene is a core part of mRNA splicing machinery and frequently contains somatic mutations that contribute to oncogenesis in myelodysplastic syndrome, acute myeloid leukemia, and other cancers. A change introduced in the GRCh38 version of the human reference build prevents detection of mutations in this gene, and others, by variant calling pipelines. This study describes the problem in detail and shows that a modified GRCh38 reference build with unchanged coordinates can be used to ameliorate the issue.

Focal disruption of DNA methylation dynamics at enhancers in IDH-mutant AML cells

Recurrent mutations in IDH1 or IDH2 in acute myeloid leukemia (AML) are associated with increased DNA methylation, but the genome-wide patterns of this hypermethylation phenotype have not been comprehensively studied in AML samples. We analyzed whole-genome bisulfite sequencing data from 15 primary AML samples with IDH1 or IDH2 mutations, which identified ~4000 focal regions that were uniquely hypermethylated in IDH^mut samples vs. normal CD34+ cells and other AMLs. These regions had modest hypermethylation in AMLs with biallelic TET2 mutations, and levels of 5-hydroxymethylation that were diminished in IDH and TET-mutant samples, indicating that this hypermethylation results from inhibition of TET-mediated demethylation. Focal hypermethylation in IDH^mut AMLs occurred at regions with low methylation in CD34+ cells, implying that DNA methylation and demethylation are active at these loci. AML samples containing IDH and DNMT3A^R882 mutations were significantly less hypermethylated, suggesting that IDH^mut-associated hypermethylation is mediated by DNMT3A. IDH^mut-specific hypermethylation was highly enriched for enhancers that form direct interactions with genes involved in normal hematopoiesis and AML, including MYC and ETV6. These results suggest that focal hypermethylation in IDH-mutant AML occurs by altering the balance between DNA methylation and demethylation, and that disruption of these pathways at enhancers may contribute to AML pathogenesis.

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.).

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.

Rapid and extraction-free detection of SARS-CoV-2 from saliva with colorimetric LAMP

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 (PPE), instrumentation, and labor. Here we present an approach to overcome these challenges with the development of 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 our optimizations of the LAMP reaction and saliva pretreatment protocols that enabled rapid and sensitive detection of < 102 viral genomes per reaction in contrived saliva controls. Moreover, our saliva pretreatment protocol enabled sensitive viral detection by conventional quantitative reverse transcription polymerase chain reaction (qRT-PCR) without RNA extraction. We validated the high performance of these assays on clinical samples and demonstrate a promising approach to overcome the current bottlenecks limiting widespread testing.

The clonal evolution of metastatic colorectal cancer

Tumor heterogeneity and evolution drive treatment resistance in metastatic colorectal cancer (mCRC). Patient-derived xenografts (PDXs) can model mCRC biology; however, their ability to accurately mimic human tumor heterogeneity is unclear. Current genomic studies in mCRC have limited scope and lack matched PDXs. Therefore, the landscape of tumor heterogeneity and its impact on the evolution of metastasis and PDXs remain undefined. We performed whole-genome, deep exome, and targeted validation sequencing of multiple primary regions, matched distant metastases, and PDXs from 11 patients with mCRC. We observed intricate clonal heterogeneity and evolution affecting metastasis dissemination and PDX clonal selection. Metastasis formation followed both monoclonal and polyclonal seeding models. In four cases, metastasis-seeding clones were not identified in any primary region, consistent with a metastasis-seeding-metastasis model. PDXs underrepresented the subclonal heterogeneity of parental tumors. These suggest that single sample tumor sequencing and current PDX models may be insufficient to guide precision medicine.

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-3. Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.

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.

No evidence that G6PD deficiency affects the efficacy or safety of daunorubicin in acute lymphoblastic leukemia induction therapy

BACKGROUND/OBJECTIVES

Anthracyclines are used in induction therapy of pediatric acute lymphoblastic leukemia (ALL) and are known to generate oxidative stress; whether this translates into enhanced antileukemic activity or hemolytic effects in patients with glucose-6-phosphate dehydrogenase (G6PD) deficiency is unknown.

DESIGN/METHODS

Among 726 pediatric patients with newly diagnosed ALL treated at St. Jude Children's Research Hospital, 22 had deficient G6PD activity. We compared the prevalence of positive minimal residual disease (MRD) ≥1% at Day 15/Day 19 of induction or ≥0.01% at Day 42/Day 46 (end of induction) and the number of red blood cell (RBC) transfusions after daunorubicin in induction between patients with or without G6PD deficiency, adjusting for ALL risk group, treatment protocol, age, and gender.

RESULTS

There was no difference in Day 15/19 (P = 1) or end of induction MRD (P = 0.76) nor in the number of RBC transfusions (P = 0.73); the lack of association with MRD was confirmed in a dataset of 1192 newly diagnosed male patients enrolled in a Children's Oncology Group trial (P = 0.78).

CONCLUSION

We found no evidence that G6PD deficiency affects daunorubicin activity during induction treatment for ALL.

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.).

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.

Recurrent WNT pathway alterations are frequent in relapsed small cell lung cancer

Nearly all patients with small cell lung cancer (SCLC) eventually relapse with chemoresistant disease. The molecular mechanisms driving chemoresistance in SCLC remain un-characterized. Here, we describe whole-exome sequencing of paired SCLC tumor samples procured at diagnosis and relapse from 12 patients, and unpaired relapse samples from 18 additional patients. Multiple somatic copy number alterations, including gains in ABCC1 and deletions in MYCL, MSH2, and MSH6, are identifiable in relapsed samples. Relapse samples also exhibit recurrent mutations and loss of heterozygosity in regulators of WNT signaling, including CHD8 and APC. Analysis of RNA-sequencing data shows enrichment for an ASCL1-low expression subtype and WNT activation in relapse samples. Activation of WNT signaling in chemosensitive human SCLC cell lines through APC knockdown induces chemoresistance. Additionally, in vitro-derived chemoresistant cell lines demonstrate increased WNT activity. Overall, our results suggest WNT signaling activation as a mechanism of chemoresistance in relapsed SCLC.

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.).

High-resolution comparative analysis of great ape genomes

Genetic studies of human evolution require high-quality contiguous ape genome assemblies that are not guided by the human reference. We coupled long-read sequence assembly and full-length complementary DNA sequencing with a multiplatform scaffolding approach to produce ab initio chimpanzee and orangutan genome assemblies. By comparing these with two long-read de novo human genome assemblies and a gorilla genome assembly, we characterized lineage-specific and shared great ape genetic variation ranging from single- to mega-base pair-sized variants. We identified ~17,000 fixed human-specific structural variants identifying genic and putative regulatory changes that have emerged in humans since divergence from nonhuman apes. Interestingly, these variants are enriched near genes that are down-regulated in human compared to chimpanzee cerebral organoids, particularly in cells analogous to radial glial neural progenitors.

Resistance-promoting effects of ependymoma treatment revealed through genomic analysis of multiple recurrences in a single patient

As in other brain tumors, multiple recurrences after complete resection and irradiation of supratentorial ependymoma are common and frequently result in patient death. This standard-of-care treatment was established in the pregenomic era without the ability to evaluate the effect that mutagenic therapies may exert on tumor evolution and in promoting resistance, recurrence, and death. We seized a rare opportunity to characterize treatment effects and the evolution of a single patient's ependymoma across four recurrences after different therapies. A combination of high-depth whole-genome and exome-based DNA sequencing of germline and tumor specimens, RNA sequencing of tumor specimens, and advanced computational analyses were used. Treatment with radiation and chemotherapies resulted in a substantial increase in mutational burden and diversification of the tumor subclonal architecture without eradication of the founding clone. Notable somatic alterations included a MEN1 driver, several epigenetic modifiers, and therapy-induced mutations that impacted multiple other cancer-relevant pathways and altered the neoantigen landscape. These genomic data provided new mechanistic insights into the genesis of ependymoma and pathways of resistance. They also revealed that radiation and chemotherapy were significant forces in shaping the increased subclonal complexity of each tumor recurrence while also failing to eradicate the founding clone. This raises the question of whether standard-of-care treatments have similar consequences in other patients with ependymoma and other types of brain tumors. If so, the perspective obtained by real-time genomic characterization of a tumor may be essential for making effective patient-specific and adaptive clinical decisions.

A multiple myeloma-specific capture sequencing platform discovers novel translocations and frequent, risk-associated point mutations in IGLL5

Multiple myeloma (MM) is a disease of copy number variants (CNVs), chromosomal translocations, and single-nucleotide variants (SNVs). To enable integrative studies across these diverse mutation types, we developed a capture-based sequencing platform to detect their occurrence in 465 genes altered in MM and used it to sequence 95 primary tumor-normal pairs to a mean depth of 104×. We detected cases of hyperdiploidy (23%), deletions of 1p (8%), 6q (21%), 8p (17%), 14q (16%), 16q (22%), and 17p (4%), and amplification of 1q (19%). We also detected IGH and MYC translocations near expected frequencies and non-silent SNVs in NRAS (24%), KRAS (21%), FAM46C (17%), TP53 (9%), DIS3 (9%), and BRAF (3%). We discovered frequent mutations in IGLL5 (18%) that were mutually exclusive of RAS mutations and associated with increased risk of disease progression (p = 0.03), suggesting that IGLL5 may be a stratifying biomarker. We identified novel IGLL5/IGH translocations in two samples. We subjected 15 of the pairs to ultra-deep sequencing (1259×) and found that although depth correlated with number of mutations detected (p = 0.001), depth past ~300× added little. The platform provides cost-effective genomic analysis for research and may be useful in individualizing treatment decisions in clinical settings.

Subclones dominate at MDS progression following allogeneic hematopoietic cell transplant

Allogeneic hematopoietic cell transplantation (alloHCT) is a potentially curative treatment for myelodysplastic syndromes (MDS), but patients who relapse after transplant have poor outcomes. In order to understand the contribution of tumor clonal evolution to disease progression,we applied exome and error-corrected targeted sequencing coupled with copy number analysis to comprehensively define changes in the clonal architecture of MDS in response to therapy using 51 serially acquired tumor samples from 9 patients who progressed after an alloHCT. We show that small subclones before alloHCT can drive progression after alloHCT. Notably, at least one subclone expanded or emerged at progression in all patients. Newly acquired structural variants (SVs) were present in an emergent/expanding subclone in 8 of 9 patients at progression, implicating the acquisition of SVs as important late subclonal progression events. In addition, pretransplant therapy with azacitidine likely influenced the mutation spectrum and evolution of emergent subclones after alloHCT. Although subclone evolution is common, founding clone mutations are always present at progression and could be detected in the bone marrow as early as 30 and/or 100 days after alloHCT in 6 of 8 (75%) patients, often prior to clinical progression. In conclusion, MDS progression after alloHCT is characterized by subclonal expansion and evolution, which can be influenced by pretransplant therapy.

An "off-the-shelf" fratricide-resistant CAR-T for the treatment of T cell hematologic malignancies

T cell malignancies represent a group of hematologic cancers with high rates of relapse and mortality in patients for whom no effective targeted therapies exist. The shared expression of target antigens between chimeric antigen receptor (CAR) T cells and malignant T cells has limited the development of CAR-T because of unintended CAR-T fratricide and an inability to harvest sufficient autologous T cells. Here we describe a fratricide resistant ‘off-the-shelf’ CAR-T (or UCART7) that targets CD7+ T cell malignancies and, through CRISPR/Cas9 gene editing, lacks both CD7 and T cell receptor alpha chain (TRAC) expression. UCART7 demonstrates efficacy against human T cell acute lymphoblastic leukemia (T-ALL) cell lines and primary T-ALL in vitro and in vivo without the induction of xenogeneic GvHD. Fratricide resistant, allo-tolerant ‘off-the-shelf’ CAR-T represents a strategy for treatment of relapsed and refractory T-ALL and non-Hodgkin’s T cell lymphoma without a requirement for autologous T cells.

Cellular stressors contribute to the expansion of hematopoietic clones of varying leukemic potential

Hematopoietic clones harboring specific mutations may expand over time. However, it remains unclear how different cellular stressors influence this expansion. Here we characterize clonal hematopoiesis after two different cellular stressors: cytotoxic therapy and hematopoietic transplantation. Cytotoxic therapy results in the expansion of clones carrying mutations in DNA damage response genes, including TP53 and PPM1D. Analyses of sorted populations show that these clones are typically multilineage and myeloid-biased. Following autologous transplantation, most clones persist with stable chimerism. However, DNMT3A mutant clones often expand, while PPM1D mutant clones often decrease in size. To assess the leukemic potential of these expanded clones, we genotyped 134 t-AML/t-MDS samples. Mutations in non-TP53 DNA damage response genes are infrequent in t-AML/t-MDS despite several being commonly identified after cytotoxic therapy. These data suggest that different hematopoietic stressors promote the expansion of distinct long-lived clones, carrying specific mutations, whose leukemic potential depends partially on the mutations they harbor.

Cellular stressors can impact clonal hematopoiesis. Here, the authors explore the impact of cytotoxic therapy and hematopoietic transplantation on clonal expansion, suggesting different stressors can promote expansion of distinct long-lived clones.

Clinical outcomes and differential effects of PI3K pathway mutation in obese versus non-obese patients with cervical cancer

The purpose of this study was to evaluate the effect of obesity and obesity-associated factors on the outcomes of patients with cervical cancer. Outcomes were evaluated in 591 patients with FIGO Ib to IV cervical cancer treated uniformly with definitive radiation. Patients were stratified into 3 groups based upon pretreatment Body Mass Index (BMI): A ≤ 18.5; B 18.6 – 34.9; and C ≥ 35. The 5-year freedom from failure rates were 58, 59, and 73% for BMI groups A, B, and C (p = 0.01). Overall survival rates were 50, 59, and 68%, respectively (p = 0.02). High expression of phosphorylated AKT (pAKT) was associated with poor outcomes only in non-obese patients. Obese patients with PI3K pathway mutant tumors had a trend toward favorable outcomes, while a similar effect was not observed in non-obese patients. Compared to similar tumors from non-obese hosts, PIK3CA and PTEN mutant tumors from obese patients failed to express high levels of phosphorylated AKT and its downstream targets. These results show that patients with obesity at the time of diagnosis of cervical cancer exhibit improved outcomes after radiation. PI3K/AKT pathway mutations are common in obese patients, but are not associated with activation of AKT signaling.

Truncating Prolactin Receptor Mutations Promote Tumor Growth in Murine Estrogen Receptor-Alpha Mammary Carcinomas

Estrogen receptor alpha-positive (ERα+) luminal tumors are the most frequent subtype of breast cancer. Stat1(-/-) mice develop mammary tumors that closely recapitulate the biological characteristics of this cancer subtype. To identify transforming events that contribute to tumorigenesis, we performed whole genome sequencing of Stat1(-/-) primary mammary tumors and matched normal tissues. This investigation identified somatic truncating mutations affecting the prolactin receptor (PRLR) in all tumor and no normal samples. Targeted sequencing confirmed the presence of these mutations in precancerous lesions, indicating that this is an early event in tumorigenesis. Functional evaluation of these heterozygous mutations in Stat1(-/-) mouse embryonic fibroblasts showed that co-expression of truncated and wild-type PRLR led to aberrant STAT3 and STAT5 activation downstream of the receptor, cellular transformation in vitro, and tumor formation in vivo. In conclusion, truncating mutations of PRLR promote tumor growth in a model of human ERα+ breast cancer and warrant further investigation.

Comprehensive discovery of noncoding RNAs in acute myeloid leukemia cell transcriptomes

To detect diverse and novel RNA species comprehensively, we compared deep small RNA and RNA sequencing (RNA-seq) methods applied to a primary acute myeloid leukemia (AML) sample. We were able to discover previously unannotated small RNAs using deep sequencing of a library method using broader insert size selection. We analyzed the long noncoding RNA (lncRNA) landscape in AML by comparing deep sequencing from multiple RNA-seq library construction methods for the sample that we studied and then integrating RNA-seq data from 179 AML cases. This identified lncRNAs that are completely novel, differentially expressed, and associated with specific AML subtypes. Our study revealed the complexity of the noncoding RNA transcriptome through a combined strategy of strand-specific small RNA and total RNA-seq. This dataset will serve as an invaluable resource for future RNA-based analyses.

Brief Report: The Role of Rare Protein-Coding Variants in Anti-Tumor Necrosis Factor Treatment Response in Rheumatoid Arthritis

OBJECTIVE

In many rheumatoid arthritis (RA) patients, disease is controlled with anti-tumor necrosis factor (anti-TNF) biologic therapies. However, in a significant number of patients, the disease fails to respond to anti-TNF therapy. We undertook the present study to examine the hypothesis that rare and low-frequency genetic variants might influence response to anti-TNF treatment.

METHODS

We sequenced the coding region of 750 genes in 1,094 RA patients of European ancestry who were treated with anti-TNF. After quality control, 690 genes were included in the analysis. We applied single-variant association and gene-based association tests to identify variants associated with anti-TNF treatment response. In addition, given the key mechanistic role of TNF, we performed gene set analyses of 27 TNF pathway genes.

RESULTS

We identified 14,420 functional variants, of which 6,934 were predicted as nonsynonymous 2,136 of which were further predicted to be "damaging." Despite the fact that the study was well powered, no single variant or gene showed study-wide significant association with change in the outcome measures disease activity or European League Against Rheumatism response. Intriguingly, we observed 3 genes, of 27 with nominal signals of association (P < 0.05), that were involved in the TNF signaling pathway. However, when we performed a rigorous gene set enrichment analysis based on association P value ranking, we observed no evidence of enrichment of association at genes involved in the TNF pathway (P_enrichment  = 0.15, based on phenotype permutations).

CONCLUSION

Our findings suggest that rare and low-frequency protein-coding variants in TNF signaling pathway genes or other genes do not contribute substantially to anti-TNF treatment response in patients with RA.

Abstract 2444: Recurrent mutations and clinical outcomes in patients with follicular lymphoma

Follicular lymphoma (FL), the most common indolent non-Hodgkin lymphoma (NHL), is associated with a highly heterogeneous clinical course. The genetic landscape of FL is less understood than other B-cell NHLs such as chronic lymphocytic leukemia and diffuse large B-cell lymphoma. With ongoing clinical trials using therapies targeting B cell receptor signaling (e.g., ibrutinib) in FL, recurrent mutations could impact the interpretation of these trial outcomes and the identification of specific patient subpopulations that have increased or decreased clinical benefit.

To identify novel recurrent mutations in FL, we designed a large custom capture lymphoma-related gene panel (WUSM-LP) and queried 113 predominantly treatment-naïve FL samples. To design the WUSM-LP reagent (NimbleGen), we targeted the coding regions of 1716 genes identified in exome sequencing data from a cohort of 24 patients (28 samples) with FL combined with genes recurrently mutated in B-cell non-Hodgkin lymphoma in the literature (Krysiak K et al., Blood, 2016). Previously, we applied the WUSM-LP to the discovery cohort and an additional 81, primarily treatment-naïve, FL samples. Using mutations identified in these 105 individuals, we identified 39 significantly mutated genes (SMGs). These included well-characterized genes in FL (e.g., KMT2D, EZH2, TNFRSF14), and identified novel mutations (BTK, HVCN1, and others) thereby expanding the repertoire of mutations affecting genes in the BCR signaling pathway, SWI/SNF complex and histone family. Including novel BTK mutations, we identified 23 (22%) individuals harboring mutations that likely up-regulate signaling downstream of BTK, potentially affecting response rates to BTK-targeting therapies at the point of FL diagnosis.

We also focused on the subset FL patients (N=59) that were sequenced within 1 year of front-line therapy in order to identify clinical associations between recurrent mutations and clinical outcomes, comparing mutation status of SMGs to progression-free survival (PFS) rates. Improved PFS was observed in patients harboring HVCN1 mutations (p&lt;0.05) in the treated cohort. CREBBP mutations were associated with significantly reduced PFS (p&lt;0.05).

Since our initial analysis, we have completed sequencing of an additional 174 FL samples with clinical data available. Confirming our earlier observations with an independent cohort, we observed 9 TNFAIP3 truncating, 18 CARD11 coiled-coil domain, 4 BCL10 truncating and 8 BTK mutations in these samples. These mutations affect 18.4% of FL patients in this cohort, are likely to up-regulate NFkappaB signaling and may impact ibrutinib sensitivity. With the additional power of 8 HVCN1 and 110 CREBBP mutations we have identified, we will evaluate the previously identified associations. The analyses of recurrent mutations and clinical correlations in this combined cohort of nearly 300 FL patients are ongoing and will be presented.

Citation Format: Kilannin Krysiak, Cody Ramirez, Felicia Gomez, Christopher A. Miller, Robert S. Fulton, Friederike Kreisel, Amanda F. Cashen, Nancy L. Bartlett, Ana Ruano, Eric D. Hsi, Malachi Griffith, Obi L. Griffith, Todd A. Fehniger. Recurrent mutations and clinical outcomes in patients with follicular lymphoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 2444. doi:10.1158/1538-7445.AM2017-2444

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.

ATRX in Diffuse Gliomas With its Mosaic/Heterogeneous Expression in a Subset

This study aims (1) to evaluate ATRX expression in different grades and subtypes of gliomas and correlate with other hallmark genetic alterations, (2) to identify and characterize mosaic/heterogeneous staining in gliomas in terms of mutation status. One hundred seventy six cases of glioma were assessed for ATRX immunohistochemistry and subdivided into positive, negative and mosaic/heterogeneous staining patterns. Five cases with heterogeneous staining were further subjected to next generation sequencing. Higher frequency of ATRX immune-negativity was detected in grade II/III astrocytic, oligoastrocytic tumors and secondary glioblastomas (GBMs), while infrequent in primary GBMs and rare in oligodendrogliomas. Loss of expression was significantly associated with IDH1 and/or TP53 mutation, while mutually exclusive with 1p/19q codeletion. Mosaic/heterogeneous staining was detected exclusively in GBMs (21.2%). Two different types of mosaic staining were identified (1) Admixture of positive and negative nuclei or intermixed mosaic and (2) Separate fragments with positive and negative/intermixed mosaic staining. ATRX mutation was identified in 2/5 (40%) cases with mosaic staining while one case showed DAXX mutation. All these cases were characterized by distinctly separate immune-negative and positive/intermixed foci. Hence, it is suggested that cases with heterogeneous staining (especially those with distinctly negative fragments) should be subjected to mutation analysis.

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, DNMT3A^R882H, 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 DNMT3A^R882 mutations. Non-leukemic hematopoietic cells with DNMT3A^R882H 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 DNMT3A^R882 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 DNMT3A^R882, while DNMT3A-dependent CpG island hypermethylation is a consequence of AML progression.

Correction: Tumor Evolution in Two Patients with Basal-like Breast Cancer: A Retrospective Genomics Study of Multiple Metastases

[This corrects the article DOI: 10.1371/journal.pmed.1002174.].

Tumor Evolution in Two Patients with Basal-like Breast Cancer: A Retrospective Genomics Study of Multiple Metastases

BACKGROUND

Metastasis is the main cause of cancer patient deaths and remains a poorly characterized process. It is still unclear when in tumor progression the ability to metastasize arises and whether this ability is inherent to the primary tumor or is acquired well after primary tumor formation. Next-generation sequencing and analytical methods to define clonal heterogeneity provide a means for identifying genetic events and the temporal relationships between these events in the primary and metastatic tumors within an individual.

METHODS AND FINDINGS

We performed DNA whole genome and mRNA sequencing on two primary tumors, each with either four or five distinct tissue site-specific metastases, from two individuals with triple-negative/basal-like breast cancers. As evidenced by their case histories, each patient had an aggressive disease course with abbreviated survival. In each patient, the overall gene expression signatures, DNA copy number patterns, and somatic mutation patterns were highly similar across each primary tumor and its associated metastases. Almost every mutation found in the primary was found in a metastasis (for the two patients, 52/54 and 75/75). Many of these mutations were found in every tumor (11/54 and 65/75, respectively). In addition, each metastasis had fewer metastatic-specific events and shared at least 50% of its somatic mutation repertoire with the primary tumor, and all samples from each patient grouped together by gene expression clustering analysis. TP53 was the only mutated gene in common between both patients and was present in every tumor in this study. Strikingly, each metastasis resulted from multiclonal seeding instead of from a single cell of origin, and few of the new mutations, present only in the metastases, were expressed in mRNAs. Because of the clinical differences between these two patients and the small sample size of our study, the generalizability of these findings will need to be further examined in larger cohorts of patients.

CONCLUSIONS

Our findings suggest that multiclonal seeding may be common amongst basal-like breast cancers. In these two patients, mutations and DNA copy number changes in the primary tumors appear to have had a biologic impact on metastatic potential, whereas mutations arising in the metastases were much more likely to be passengers.

Deregulation of DUX4 and ERG in acute lymphoblastic leukemia

Chromosomal rearrangements deregulating hematopoietic transcription factors are common in acute lymphoblastic leukemia (ALL). Here we show that deregulation of the homeobox transcription factor gene DUX4 and the ETS transcription factor gene ERG is a hallmark of a subtype of B-progenitor ALL that comprises up to 7% of B-ALL. DUX4 rearrangement and overexpression was present in all cases and was accompanied by transcriptional deregulation of ERG, expression of a novel ERG isoform, ERGalt, and frequent ERG deletion. ERGalt uses a non-canonical first exon whose transcription was initiated by DUX4 binding. ERGalt retains the DNA-binding and transactivation domains of ERG, but it inhibits wild-type ERG transcriptional activity and is transforming. These results illustrate a unique paradigm of transcription factor deregulation in leukemia in which DUX4 deregulation results in loss of function of ERG, either by deletion or induced expression of an isoform that is a dominant-negative inhibitor of wild-type ERG function.

TP53 and Decitabine in Acute Myeloid Leukemia and Myelodysplastic Syndromes

BACKGROUND

The molecular determinants of clinical responses to decitabine therapy in patients with acute myeloid leukemia (AML) or myelodysplastic syndromes (MDS) are unclear.

METHODS

We enrolled 84 adult patients with AML or MDS in a single-institution trial of decitabine to identify somatic mutations and their relationships to clinical responses. Decitabine was administered at a dose of 20 mg per square meter of body-surface area per day for 10 consecutive days in monthly cycles. We performed enhanced exome or gene-panel sequencing in 67 of these patients and serial sequencing at multiple time points to evaluate patterns of mutation clearance in 54 patients. An extension cohort included 32 additional patients who received decitabine in different protocols.

RESULTS

Of the 116 patients, 53 (46%) had bone marrow blast clearance (<5% blasts). Response rates were higher among patients with an unfavorable-risk cytogenetic profile than among patients with an intermediate-risk or favorable-risk cytogenetic profile (29 of 43 patients [67%] vs. 24 of 71 patients [34%], P<0.001) and among patients with TP53 mutations than among patients with wild-type TP53 (21 of 21 [100%] vs. 32 of 78 [41%], P<0.001). Previous studies have consistently shown that patients with an unfavorable-risk cytogenetic profile and TP53 mutations who receive conventional chemotherapy have poor outcomes. However, in this study of 10-day courses of decitabine, neither of these risk factors was associated with a lower rate of overall survival than the rate of survival among study patients with intermediate-risk cytogenetic profiles.

CONCLUSIONS

Patients with AML and MDS who had cytogenetic abnormalities associated with unfavorable risk, TP53 mutations, or both had favorable clinical responses and robust (but incomplete) mutation clearance after receiving serial 10-day courses of decitabine. Although these responses were not durable, they resulted in rates of overall survival that were similar to those among patients with AML who had an intermediate-risk cytogenetic profile and who also received serial 10-day courses of decitabine. (Funded by the National Cancer Institute and others; ClinicalTrials.gov number, NCT01687400 .).