The utility of cancer whole genome and transcriptome sequencing (cWGTS) in oncology is increasingly recognized. However, implementation of cWGTS is challenged by the need to deliver results within ...clinically relevant timeframes, concerns about assay sensitivity, reporting and prioritization of findings. In a prospective research study we develop a workflow that reports comprehensive cWGTS results in 9 days. Comparison of cWGTS to diagnostic panel assays demonstrates the potential of cWGTS to capture all clinically reported mutations with comparable sensitivity in a single workflow. Benchmarking identifies a minimum of 80× as optimal depth for clinical WGS sequencing. Integration of germline, somatic DNA and RNA-seq data enable data-driven variant prioritization and reporting, with oncogenic findings reported in 54% more patients than standard of care. These results establish key technical considerations for the implementation of cWGTS as an integrated test in clinical oncology.
Clinical recommendations for Acute Myeloid Leukemia (AML) classification and risk-stratification remain heavily reliant on cytogenetic findings at diagnosis, which are present in <50% of patients. ...Using comprehensive molecular profiling data from 3,653 patients we characterize and validate 16 molecular classes describing 100% of AML patients. Each class represents diverse biological AML subgroups, and is associated with distinct clinical presentation, likelihood of response to induction chemotherapy, risk of relapse and death over time. Secondary AML-2, emerges as the second largest class (24%), associates with high-risk disease, poor prognosis irrespective of flow Minimal Residual Disease (MRD) negativity, and derives significant benefit from transplantation. Guided by class membership we derive a 3-tier risk-stratification score that re-stratifies 26% of patients as compared to standard of care. This results in a unified framework for disease classification and risk-stratification in AML that relies on information from cytogenetics and 32 genes. Last, we develop an open-access patient-tailored clinical decision support tool.
•UBA1 mutations were identified in 1% of patients with MDS, enriched in patients with few myeloid mutations or established disease classification.•Inflammatory clinical presentation and vacuoles were ...observed in 83% and 71%, respectively, of patients with pathogenic UBA1 mutations.
Display omitted
Mutations in UBA1, which are disease-defining for VEXAS (vacuoles, E1 enzyme, X-linked, autoinflammatory, somatic) syndrome, have been reported in patients diagnosed with myelodysplastic syndromes (MDSs). Here, we define the prevalence and clinical associations of UBA1 mutations in a representative cohort of patients with MDS. Digital droplet polymerase chain reaction profiling of a selected cohort of 375 male patients lacking MDS disease-defining mutations or established World Health Organization (WHO) disease classification identified 28 patients (7%) with UBA1 p.M41T/V/L mutations. Using targeted sequencing of UBA1 in a representative MDS cohort (n = 2027), we identified an additional 27 variants in 26 patients (1%), which we classified as likely/pathogenic (n = 12) and of unknown significance (n = 15). Among the total 40 patients with likely/pathogenic variants (2%), all were male and 63% were classified by WHO 2016 criteria as MDS with multilineage dysplasia or MDS with single-lineage dysplasia. Patients had a median of additional myeloid gene mutation, often in TET2 (n = 12), DNMT3A (n = 10), ASXL1 (n = 3), or SF3B1 (n = 3). Retrospective clinical review, if possible, showed that 83% (28/34) UBA1-mutant cases had VEXAS syndrome–associated diagnoses or inflammatory clinical presentation. The prevalence of UBA1 mutations in patients with MDS argues for systematic screening for UBA1 in the management of MDS.
The goal of therapy for patients with essential thrombocythemia (ET) and polycythemia vera (PV) is to reduce thrombotic events by normalizing blood counts. Hydroxyurea (HU) and interferon-α (IFN-α) ...are the most frequently used cytoreductive options for patients with ET and PV at high risk for vascular complications. Myeloproliferative Disorders Research Consortium 112 was an investigator-initiated, phase 3 trial comparing HU to pegylated IFN-α (PEG) in treatment-naïve, high-risk patients with ET/PV. The primary endpoint was complete response (CR) rate at 12 months. A total of 168 patients were treated for a median of 81.0 weeks. CR for HU was 37% and 35% for PEG (P = .80) at 12 months. At 24 to 36 months, CR was 20% to 17% for HU and 29% to 33% for PEG. PEG led to a greater reduction in JAK2V617F at 24 months, but histopathologic responses were more frequent with HU. Thrombotic events and disease progression were infrequent in both arms, whereas grade 3/4 adverse events were more frequent with PEG (46% vs 28%). At 12 months of treatment, there was no significant difference in CR rates between HU and PEG. This study indicates that PEG and HU are both effective treatments for PV and ET. With longer treatment, PEG was more effective in normalizing blood counts and reducing driver mutation burden, whereas HU produced more histopathologic responses. Despite these differences, both agents did not differ in limiting thrombotic events and disease progression in high-risk patients with ET/PV. This trial was registered at www.clinicaltrials.gov as #NCT01259856.
•Rates of thrombosis and progression were low in patients with ET/PV treated with either HU or IFN in this randomized study.•PEG was more effective in normalizing counts and reducing JAK2V617F VAF in PV whereas HU induced more HPRs in ET.
Display omitted
Recent genomic research efforts in multiple myeloma have revealed clinically relevant molecular subgroups beyond conventional cytogenetic classifications. Implementing these advances in clinical ...trial design and in routine patient care requires a new generation of molecular diagnostic tools. Here, we present a custom capture next-generation sequencing (NGS) panel designed to identify rearrangements involving the IGH locus, arm level, and focal copy number aberrations, as well as frequently mutated genes in multiple myeloma in a single assay. We sequenced 154 patients with plasma cell disorders and performed a head-to-head comparison with the results from conventional clinical assays, i.e., fluorescent in situ hybridization (FISH) and single-nucleotide polymorphism (SNP) microarray. Our custom capture NGS panel had high sensitivity (>99%) and specificity (>99%) for detection of IGH translocations and relevant chromosomal gains and losses in multiple myeloma. In addition, the assay was able to capture novel genomic markers associated with poor outcome such as bi-allelic events involving TP53. In summary, we show that a multiple myeloma designed custom capture NGS panel can detect IGH translocations and CNAs with very high concordance in relation to FISH and SNP microarrays and importantly captures the most relevant and recurrent somatic mutations in multiple myeloma rendering this approach highly suitable for clinical application in the modern era.
BACKGROUND: Risk stratification and therapeutic decision-making for myelodysplastic syndromes (MDS) are based on the International Prognostic Scoring System–Revised (IPSS-R), which considers ...hematologic parameters and cytogenetic abnormalities. Somatic gene mutations are not yet used in the risk stratification of patients with MDS. METHODS: To develop a clinical-molecular prognostic model (IPSS-Molecular IPSS-M), pretreatment diagnostic or peridiagnostic samples from 2957 patients with MDS were profiled for mutations in 152 genes. Clinical and molecular variables were evaluated for associations with leukemia-free survival, leukemic transformation, and overall survival. Feature selection was applied to determine the set of independent IPSS-M prognostic variables. The relative weights of the selected variables were estimated using a robust Cox multivariable model adjusted for confounders. The IPSS-M was validated in an external cohort of 754 Japanese patients with MDS. RESULTS: We mapped at least one oncogenic genomic alteration in 94% of patients with MDS. Multivariable analysis identified TP53multihit, FLT3 mutations, and MLLPTD as top genetic predictors of adverse outcomes. Conversely, SF3B1 mutations were associated with favorable outcomes, but this was modulated by patterns of comutation. Using hematologic parameters, cytogenetic abnormalities, and somatic mutations of 31 genes, the IPSS-M resulted in a unique risk score for individual patients. We further derived six IPSS-M risk categories with prognostic differences. Compared with the IPSS-R, the IPSS-M improved prognostic discrimination across all clinical end points and restratified 46% of patients. The IPSS-M was applicable in primary and secondary/therapy-related MDS. To simplify clinical use of the IPSS-M, we developed an open-access Web calculator that accounts for missing values. CONCLUSIONS: Combining genomic profiling with hematologic and cytogenetic parameters, the IPSS-M improves the risk stratification of patients with MDS and represents a valuable tool for clinical decision-making. (Funded by Celgene Corporation through the MDS Foundation, the Josie Robertson Investigators Program, the Edward P. Evans Foundation, the Projects of National Relevance of the Italian Ministry of University and Research, Associazione Italiana per la Ricerca sul Cancro, the Japan Agency for Medical Research and Development, Cancer Research UK, the Austrian Science Fund, the MEXT Japanese Ministry of Education, Culture, Sports, Science and Technology Program for Promoting Research on the Supercomputer Fugaku, the Japan Society for the Promotion of Science, the Taiwan Department of Health, and Celgene Corporation through the MDS Foundation.)
Multiple large-scale genomic profiling efforts have been undertaken in osteosarcoma to define the genomic drivers of tumorigenesis, therapeutic response, and disease recurrence. The spatial and ...temporal intratumor heterogeneity could also play a role in promoting tumor growth and treatment resistance. We conducted longitudinal whole-genome sequencing of 37 tumor samples from 8 patients with relapsed or refractory osteosarcoma. Each patient had at least one sample from a primary site and a metastatic or relapse site. Subclonal copy-number alterations were identified in all patients except one. In 5 patients, subclones from the primary tumor emerged and dominated at subsequent relapses. MYC gain/amplification was enriched in the treatment-resistant clones in 6 of 7 patients with multiple clones. Amplifications in other potential driver genes, such as CCNE1, RAD21, VEGFA, and IGF1R, were also observed in the resistant copy-number clones. A chromosomal duplication timing analysis revealed that complex genomic rearrangements typically occurred prior to diagnosis, supporting a macroevolutionary model of evolution, where a large number of genomic aberrations are acquired over a short period of time followed by clonal selection, as opposed to ongoing evolution. A mutational signature analysis of recurrent tumors revealed that homologous repair deficiency (HRD)-related SBS3 increases at each time point in patients with recurrent disease, suggesting that HRD continues to be an active mutagenic process after diagnosis. Overall, by examining the clonal relationships between temporally and spatially separated samples from patients with relapsed/refractory osteosarcoma, this study sheds light on the intratumor heterogeneity and potential drivers of treatment resistance in this disease.
The chemoresistant population in recurrent osteosarcoma is subclonal at diagnosis, emerges at the time of primary resection due to selective pressure from neoadjuvant chemotherapy, and is characterized by unique oncogenic amplifications.
SF3B1K700E is the most frequent mutation in myelodysplastic syndrome (MDS), but the mechanisms by which it drives MDS pathogenesis remain unclear. We derived a panel of 18 genetically matched ...SF3B1K700E- and SF3B1WT-induced pluripotent stem cell (iPSC) lines from patients with MDS with ring sideroblasts (MDS-RS) harboring isolated SF3B1K700E mutations and performed RNA and ATAC sequencing in purified CD34+/CD45+ hematopoietic stem/progenitor cells (HSPCs) derived from them. We developed a novel computational framework integrating splicing with transcript usage and gene expression analyses and derived a SF3B1K700E splicing signature consisting of 59 splicing events linked to 34 genes, which associates with the SF3B1 mutational status of primary MDS patient cells. The chromatin landscape of SF3B1K700E HSPCs showed increased priming toward the megakaryocyte- erythroid lineage. Transcription factor motifs enriched in chromatin regions more accessible in SF3B1K700E cells included, unexpectedly, motifs of the TEA domain (TEAD) transcription factor family. TEAD expression and transcriptional activity were upregulated in SF3B1-mutant iPSC-HSPCs, in support of a Hippo pathway-independent role of TEAD as a potential novel transcriptional regulator of SF3B1K700E cells. This study provides a comprehensive characterization of the transcriptional and chromatin landscape of SF3B1K700E HSPCs and nominates novel mis-spliced genes and transcriptional programs with putative roles in MDS-RS disease biology.
•Genetically matched MDS-RS and normal patient-specific iPSC-HSPCs are used to derive a mutant SF3B1 splicing signature.•Integrated transcriptomics and chromatin accessibility nominate TEAD as a putative novel transcriptional regulator of SF3B1K700E cells.
Display omitted
▪
The somatic hotspot mutation SF3B1K700E is characteristically found in myelodysplastic syndrome with ring sideroblasts (MDS-RS) and frequently occurs as an isolated mutation. However, our ...understanding of how this mutation drives MDS pathogenesis remains limited. To explore the downstream consequences of the SF3B1K700E mutation and its role in disease pathogenesis, we generated a panel of isogenic SF3B1K700E and SF3B1WT induced pluripotent stem cell (iPSC) lines from 3 MDS-RS patients with isolated SF3B1K700E mutation (3 SF3B1K700E and 3 SF3B1WT lines per patient, total 18). Upon hematopoietic differentiation, SF3B1K700E cells exhibited lower growth and colony-forming ability, compared to SF3B1WT cells, recapitulating hallmark phenotypes of MDS cells.
To investigate the effects of the SF3B1K700E mutation on the transcriptome and chromatin landscape, we performed RNA- and ATAC- sequencing in purified CD34+/CD45+ hematopoietic stem/progenitor cells (HSPCs) derived from the panel of the 18 isogenic SF3B1K700E and SF3B1WT iPSC lines. Principal component analysis (PCA) and hierarchical clustering based on gene expression grouped the iPSC lines primarily by genotype (SF3B1K700E vs SF3B1WT) and secondarily by genetic background. To assess the impact of the SF3B1K700E mutation at the exon, transcript and gene level, we developed an analytical framework integrating differential splicing with differential transcript usage and differential gene expression analyses. We thus discovered 59 splicing events linked to 34 genes (most statistically significant events that also mapped to differentially used transcripts and differentially expressed genes). This SF3B1K700E splicing signature includes genes previously reported as mis-spliced in SF3B1K700E cells (e.g BRD9, ABCB7), as well as novel genes. We tested this signature against a published dataset of primary MDS patient samples (Pellagatti et al.). PCA based on the inclusion level of the splicing events of our signature separated SF3B1-mutated MDS patients from patients without splicing factor mutations (SF-WT) or healthy individuals. Furthermore, it identified one patient erroneously annotated as SF-WT that clustered together with the SF3B1-mutated patients, who had a, previously overlooked, 6bp in-frame deletion spanning the K700E hotspot.
By comparing the chromatin accessibility profiles of SF3B1K700E and SF3B1WT iPSC-HSPCs to those defined in primary human cell types along the hematopoietic hierarchy (Corces et al.), we found that the chromatin landscape of SF3B1K700E HSPCs resembled more this of megakaryocyte-erythroid progenitor cells (MEPs) and erythroid cells, whereas that of SF3B1WT HSPCs resembled more granulocyte-monocyte progenitors (GMPs) and monocytes. This finding may underlie the more prominent involvement of the erythroid lineage in the pathology and clinical presentation of MDS-RS. To interrogate transcriptional programs in SF3B1K700E mutant cells, we performed transcription factor (TF) motif enrichment analysis. Motifs enriched in ATAC-Seq peaks more accessible in SF3B1K700E cells that were linked to genes upregulated in SF3B1K700E cells, included motifs of several TFs with known roles in hematopoiesis (GATA, ETS, STAT, AP-1). Unexpectedly, motifs of the TEAD TFs were also enriched. The TEAD family of TFs are best known as effectors of the Hippo signaling pathway, with important roles in various biological processes and malignancies, albeit no clear links to adult hematopoiesis or hematologic disease. TEAD2 and TEAD4 were upregulated in SF3B1-mutant, compared to the WT, iPSC-HSPCs and TEAD transcriptional activity, measured with a luciferase reporter construct, was higher in SF3B1K700E, compared to SF3B1WT iPSC-HSPCs. We did not find expression or activation of YAP or TAZ, which bind to DNA as a complex with TEAD upon Hippo pathway activation. These results support a Hippo-independent increase of TEAD expression and activity in SF3B1K700E cells.
In summary, we generated a panel of isogenic patient-derived iPSCs that allowed us to comprehensively characterize the transcriptome and chromatin landscape of SF3B1K700E HSPCs in an isogenic system, derive a splicing signature of SF3B1K700E and identify the TEAD TF as a new transcriptional regulator of SF3B1K700E mutant HSPCs.
G. Asimomitis, A.G. Deslauriers: shared 1 st authorship
E. Papaemmanuil, E.P. Papapetrou: shared senior authorship
Deslauriers: Novo Nordisk A/S: Current Employment. Hellström-Lindberg: Celgene: Research Funding. Papaemmanuil: Isabl Technologies: Divested equity in a private or publicly-traded company in the past 24 months; Kyowa Hakko Kirin Pharma: Consultancy.