In-depth characterization of intratumoral heterogeneity in refractory B-cell non-Hodgkin lymphoma through the lens of a Research Autopsy Program.

Intratumoral heterogeneity (ITH) provides the substrate for tumor evolution and treatment resistance, yet is remarkably understudied in lymphoma, due to the often limited amount of tissue that gets sampled during the routine diagnostic process, generally from a single nodal or extranodal site. Furthermore, the trajectory of how lymphoma, and especially non-Hodgkin lymphoma, spreads throughout the human body remains poorly understood. Here, we present a detailed characterization of ITH by applying whole-genome sequencing to spatially separated tumor samples harvested at the time of autopsy (n=24) and/or diagnosis (n=3) in three patients presenting with refractory B-cell non-Hodgkin lymphoma. Through deconvolution of bulk samples into clonal mixtures and inference of phylogenetic trees, we found evidence that polyclonal seeding underlies tumor dissemination in lymphoma. We identify mutation signatures associated with ancestral and descendant clones. In our series of patients with highly refractory lymphoma, the determinants of resistance were often harbored by founding clones, although there was also evidence of positive selection of driver mutations, likely under the influence of therapy. Lastly, we show that circulating tumor DNA is suitable for the detection of ancestral mutations but may miss a significant proportion of private mutations that can be detected in tissue. Our study clearly shows the existence of intricate patterns of regional and anatomical evolution that can only be disentangled through multi-regional tumor tissue profiling.

Prospective observational study and serosurvey of SARS-CoV-2 infection in asymptomatic healthcare workers at a Canadian tertiary care center.

Health care workers (HCWs) are at higher risk for SARS-CoV-2 infection and may play a role in transmitting the infection to vulnerable patients and members of the community. This is particularly worrisome in the context of asymptomatic infection. We performed a cross-sectional study looking at asymptomatic SARS-CoV-2 infection in HCWs. We screened asymptomatic HCWs for SARS-CoV-2 via PCR. Complementary viral genome sequencing was performed on positive swab specimens. A seroprevalence analysis was also performed using multiple assays. Asymptomatic health care worker cohorts had a combined swab positivity rate of 29/5776 (0.50%, 95%CI 0.32-0.75) relative to a comparative cohort of symptomatic HCWs, where 54/1597 (3.4%) tested positive for SARS-CoV-2 (ratio of symptomatic to asymptomatic 6.8:1). SARS-CoV-2 seroprevalence among 996 asymptomatic HCWs with no prior known exposure to SARS-CoV-2 was 1.4-3.4%, depending on assay. A novel in-house Coronavirus protein microarray showed differing SARS-CoV-2 protein reactivities and helped define likely true positives vs. suspected false positives. Our study demonstrates the utility of routine screening of asymptomatic HCWs, which may help to identify a significant proportion of infections.

Landscape of somatic single nucleotide variants and indels in colorectal cancer and impact on survival.

Colorectal cancer (CRC) is a biologically heterogeneous disease. To characterize its mutational profile, we conduct targeted sequencing of 205 genes for 2,105 CRC cases with survival data. Our data shows several findings in addition to enhancing the existing knowledge of CRC. We identify PRKCI, SPZ1, MUTYH, MAP2K4, FETUB, and TGFBR2 as additional genes significantly mutated in CRC. We find that among hypermutated tumors, an increased mutation burden is associated with improved CRC-specific survival (HR = 0.42, 95% CI: 0.21-0.82). Mutations in TP53 are associated with poorer CRC-specific survival, which is most pronounced in cases carrying TP53 mutations with predicted 0% transcriptional activity (HR = 1.53, 95% CI: 1.21-1.94). Furthermore, we observe differences in mutational frequency of several genes and pathways by tumor location, stage, and sex. Overall, this large study provides deep insights into somatic mutations in CRC, and their potential relationships with survival and tumor features.

GATA6 Expression Distinguishes Classical and Basal-like Subtypes in Advanced Pancreatic Cancer.

PURPOSE: To determine the impact of basal-like and classical subtypes in advanced pancreatic ductal adenocarcinoma (PDAC) and to explore GATA6 expression as a surrogate biomarker. EXPERIMENTAL DESIGN: Within the COMPASS trial, patients proceeding to chemotherapy for advanced PDAC undergo tumor biopsy for RNA-sequencing (RNA-seq). Overall response rate (ORR) and overall survival (OS) were stratified by subtypes and according to chemotherapy received. Correlation of GATA6 with the subtypes using gene expression profiling, in situ hybridization (ISH) was explored. RESULTS: Between December 2015 and May 2019, 195 patients (95%) had enough tissue for RNA-seq; 39 (20%) were classified as basal-like and 156 (80%) as classical. RECIST response data were available for 157 patients; 29 basal-like and 128 classical where the ORR was 10% versus 33%, respectively (P = 0.02). In patients with basal-like tumors treated with modified FOLFIRINOX (n = 22), the progression rate was 60% compared with 15% in classical PDAC (P = 0.0002). Median OS in the intention-to-treat population (n = 195) was 9.3 months for classical versus 5.9 months for basal-like PDAC (HR, 0.47; 95% confidence interval, 0.32-0.69; P = 0.0001). GATA6 expression by RNA-seq highly correlated with the classifier (P < 0.001) and ISH predicted the subtypes with sensitivity of 89% and specificity of 83%. In a multivariate analysis, GATA6 expression was prognostic (P = 0.02). In exploratory analyses, basal-like tumors, could be identified by keratin 5, were more hypoxic and enriched for a T-cell-inflamed gene expression signature. CONCLUSIONS: The basal-like subtype is chemoresistant and can be distinguished from classical PDAC by GATA6 expression.See related commentary by Collisson, p. 4715.

Author Correction: Transcription phenotypes of pancreatic cancer are driven by genomic events during tumor evolution.

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

Transcription phenotypes of pancreatic cancer are driven by genomic events during tumor evolution.

Pancreatic adenocarcinoma presents as a spectrum of a highly aggressive disease in patients. The basis of this disease heterogeneity has proved difficult to resolve due to poor tumor cellularity and extensive genomic instability. To address this, a dataset of whole genomes and transcriptomes was generated from purified epithelium of primary and metastatic tumors. Transcriptome analysis demonstrated that molecular subtypes are a product of a gene expression continuum driven by a mixture of intratumoral subpopulations, which was confirmed by single-cell analysis. Integrated whole-genome analysis uncovered that molecular subtypes are linked to specific copy number aberrations in genes such as mutant KRAS and GATA6. By mapping tumor genetic histories, tetraploidization emerged as a key mutational process behind these events. Taken together, these data support the premise that the constellation of genomic aberrations in the tumor gives rise to the molecular subtype, and that disease heterogeneity is due to ongoing genomic instability during progression.

Combined EsophaCap cytology and MUC2 immunohistochemistry for screening of intestinal metaplasia, dysplasia and carcinoma.

Purpose: The incidence of esophageal adenocarcinoma (EAC) has increased by 700% in Western countries over the last 30 years. Although clinical guidelines call for endoscopic surveillance for EAC among high-risk populations, fewer than 5% of new EAC patients are under surveillance at the time of diagnosis. We studied the accuracy of combined cytopathology and MUC2 immunohistochemistry (IHC) for screening of Intestinal Metaplasia (IM), dysplasia and EAC, using specimens collected from the EsophaCap swallowable encapsulated cytology sponge from Canada and United States. Patients and methods: By comparing the EsophaCap cytological diagnosis with concurrent endoscopic biopsies performed on the same patients in 28 cases, we first built up the cytology diagnostic categories and criteria. Based on these criteria, 136 cases were evaluated by both cytology and MUC2 IHC with blinded to patient biopsy diagnosis. Results: We first set up categories and criteria for cytological diagnosis of EscophaCap samples. Based on these, we divided our evaluated cytological samples into two groups: non-IM group and IM or dysplasia or adenocarcinoma group. Using the biopsy as our gold standard to screen IM, dysplasia and EAC by combined cytology and MUC2 IHC, the sensitivity and specificity were 68% and 91%, respectively, which is in the range of clinically useful cytological screening tests such as the cervical Pap smear. Conclusions: Combined EsophaCap cytology and MUC2 IHC could be a good screening test for IM and Beyond.

Detection of Circulating Tumor DNA in Plasma: A Potential Biomarker for Esophageal Adenocarcinoma.

BACKGROUND: Recent literature has demonstrated the potential of "liquid biopsy" and detection of circulating tumor (ct)DNA as a cancer biomarker. However, to date there is a lack of data specific to esophageal adenocarcinoma (EAC). This study was conducted to determine how detection and quantification of ctDNA changes with disease burden in patients with EAC and evaluate its potential as a biomarker in this population. METHODS: Blood samples were obtained from patients with stage I to IV EAC. Longitudinal blood samples were collected from a subset of patients. Imaging studies and pathology reports were reviewed to determine disease course. Tumor samples were sequenced to identify mutations. Mutations in plasma DNA were detected using custom, barcoded, patient-specific sequencing libraries. Mutations in plasma were quantified, and associations with disease stage and response to therapy were explored. RESULTS: Plasma samples from a final cohort of 38 patients were evaluated. Baseline plasma samples were ctDNA positive for 18 patients (47%) overall, with tumor allele frequencies ranging from 0.05% to 5.30%. Detection frequency of ctDNA and quantity of ctDNA increased with stage. Data from longitudinal samples indicate that ctDNA levels correlate with and precede evidence of response to therapy or recurrence. CONCLUSIONS: ctDNA can be detected in plasma of EAC patients and correlates with disease burden. Detection of ctDNA in early-stage EAC is challenging and may limit diagnostic applications. However, our data demonstrate the potential of ctDNA as a dynamic biomarker to monitor treatment response and disease recurrence in patients with EAC.

Regional perturbation of gene transcription is associated with intrachromosomal rearrangements and gene fusion transcripts in high grade ovarian cancer.

Genomic rearrangements are a hallmark of cancer biology and progression, allowing cells to rapidly transform through alterations in regulatory structures, changes in expression patterns, reprogramming of signaling pathways, and creation of novel transcripts via gene fusion events. Though functional gene fusions encoding oncogenic proteins are the most dramatic outcomes of genomic rearrangements, we investigated the relationship between rearrangements evidenced by fusion transcripts and local expression changes in cancer using transcriptome data alone. 9,953 gene fusion predictions from 418 primary serious ovarian cancer tumors were analyzed, identifying depletions of gene fusion breakpoints within coding regions of fused genes as well as an N-terminal enrichment of breakpoints within fused genes. We identified 48 genes with significant fusion-associated upregulation and furthermore demonstrate that significant regional overexpression of intact genes in patient transcriptomes occurs within 1 megabase of 78 novel gene fusions that function as central markers of these regions. We reveal that cancer transcriptomes select for gene fusions that preserve protein and protein domain coding potential. The association of gene fusion transcripts with neighboring gene overexpression supports rearrangements as mechanism through which cancer cells remodel their transcriptomes and identifies a new way to utilize gene fusions as indicators of regional expression changes in diseased cells with only transcriptomic data.

Integration of Genomic and Transcriptional Features in Pancreatic Cancer Reveals Increased Cell Cycle Progression in Metastases.

We integrated clinical, genomic, and transcriptomic data from 224 primaries and 95 metastases from 289 patients to characterize progression of pancreatic ductal adenocarcinoma (PDAC). Driver gene alterations and mutational and expression-based signatures were preserved, with truncations, inversions, and translocations most conserved. Cell cycle progression (CCP) increased with sequential inactivation of tumor suppressors, yet remained higher in metastases, perhaps driven by cell cycle regulatory gene variants. Half of the cases were hypoxic by expression markers, overlapping with molecular subtypes. Paired tumor heterogeneity showed cancer cell migration by Halstedian progression. Multiple PDACs arising synchronously and metachronously in the same pancreas were actually intra-parenchymal metastases, not independent primary tumors. Established clinical co-variates dominated survival analyses, although CCP and hypoxia may inform clinical practice.

Plasma circulating tumor DNA as a potential tool for disease monitoring in head and neck cancer.

BACKGROUND: Recommendations for perioperative therapy in head and neck cancer are not explicit and recurrence occurs frequently. Circulating tumor DNA is an emerging cancer biomarker, but has not been extensively explored for detection of recurrence in head and neck cancer. METHODS: Patients diagnosed with head and neck squamous cell carcinoma were recruited into the study protocol. Tumors were sequenced to identify patient-specific mutations. Mutations were then identified in plasma circulating tumor DNA from pre-treatment blood samples and longitudinally during standard follow-up. Circulating tumor DNA status during follow-up was correlated to disease recurrence. RESULTS: Samples were taken from eight patients. Tumor mutations were verified in seven patients. Baseline circulating tumor DNA was positive in six patients. Recurrence occurred in four patients, two of whom had detectable circulating tumor DNA prior to recurrence. CONCLUSION: Circulating tumor DNA is a potential tool for disease and recurrence monitoring following curative therapy in head and neck cancer, allowing for better prognostication, and/or modification of treatment strategies.

Prediction of acute myeloid leukaemia risk in healthy individuals.

The incidence of acute myeloid leukaemia (AML) increases with age and mortality exceeds 90% when diagnosed after age 65. Most cases arise without any detectable early symptoms and patients usually present with the acute complications of bone marrow failure1. The onset of such de novo AML cases is typically preceded by the accumulation of somatic mutations in preleukaemic haematopoietic stem and progenitor cells (HSPCs) that undergo clonal expansion2,3. However, recurrent AML mutations also accumulate in HSPCs during ageing of healthy individuals who do not develop AML, a phenomenon referred to as age-related clonal haematopoiesis (ARCH)4-8. Here we use deep sequencing to analyse genes that are recurrently mutated in AML to distinguish between individuals who have a high risk of developing AML and those with benign ARCH. We analysed peripheral blood cells from 95 individuals that were obtained on average 6.3 years before AML diagnosis (pre-AML group), together with 414 unselected age- and gender-matched individuals (control group). Pre-AML cases were distinct from controls and had more mutations per sample, higher variant allele frequencies, indicating greater clonal expansion, and showed enrichment of mutations in specific genes. Genetic parameters were used to derive a model that accurately predicted AML-free survival; this model was validated in an independent cohort of 29 pre-AML cases and 262 controls. Because AML is rare, we also developed an AML predictive model using a large electronic health record database that identified individuals at greater risk. Collectively our findings provide proof-of-concept that it is possible to discriminate ARCH from pre-AML many years before malignant transformation. This could in future enable earlier detection and monitoring, and may help to inform intervention.

Impact of Chemical-Induced Mutational Load Increase on Immune Checkpoint Therapy in Poorly Responsive Murine Tumors.

A recurring historic finding in cancer drug development is encouraging antitumor effects observed in tumor-bearing mice that fail to translate into the clinic. An intriguing exception to this pattern is immune checkpoint therapy, as the sustained tumor regressions observed in subsets of cancer patients are rare in mice. Reasoning that this may be due in part to relatively low mutational loads of mouse tumors, we mutagenized transplantable mouse tumor cell lines EMT-6/P, B16F1, RENCA, CT26, and MC38 in vitro with methylnitro-nitrosoguanidine (MNNG) or ethylmethane sulfonate (EMS) and tested their responsiveness to PD-L1 blockade. Exome sequencing confirmed an increase in somatic mutations by mutagen treatment, an effect mimicked in EMT-6 variants chronically exposed in vivo to cisplatin or cyclophosphamide. Certain mutagenized variants of B16F1, EMT-6/P, CT26, and MC38 (but not RENCA) were more immunogenic than their parents, yet anti-PD-L1 sensitization developed only in some EMT-6/P and B16F1 variants. Treatment response patterns corresponded with changes in immune cell infiltration and especially increases in CD8+ T cells. Chronically cisplatin-exposed EMT-6 variants were also more responsive to anti-PD-L1 therapy. Although tumor PD-L1 expression was upregulated in in vivo chemotherapy-exposed variants, PD-L1 expression levels were not consistently associated with anti-PD-L1 treatment activity across mutagenized or chemotherapy-exposed variants. In summary, mutagenized and more immunogenic mouse tumors were not universally sensitized to PD-L1 blockade. Chemically mutagenized variants may be useful to evaluate the impact of immunologically "hot" or "cold" tumors with a high mutational load, to which certain chemotherapy agents may contribute, on immunotherapy outcomes. Mol Cancer Ther; 17(4); 869-82. ©2018 AACR.

Genomics-Driven Precision Medicine for Advanced Pancreatic Cancer: Early Results from the COMPASS Trial.

Purpose: To perform real-time whole genome sequencing (WGS) and RNA sequencing (RNASeq) of advanced pancreatic ductal adenocarcinoma (PDAC) to identify predictive mutational and transcriptional features for better treatment selection.Experimental Design: Patients with advanced PDAC were prospectively recruited prior to first-line combination chemotherapy. Fresh tumor tissue was acquired by image-guided percutaneous core biopsy for WGS and RNASeq. Laser capture microdissection was performed for all cases. Primary endpoint was feasibility to report WGS results prior to first disease assessment CT scan at 8 weeks. The main secondary endpoint was discovery of patient subsets with predictive mutational and transcriptional signatures.Results: Sixty-three patients underwent a tumor biopsy between December 2015 and June 2017. WGS and RNASeq were successful in 62 (98%) and 60 (95%), respectively. Genomic results were reported at a median of 35 days (range, 19-52 days) from biopsy, meeting the primary feasibility endpoint. Objective responses to first-line chemotherapy were significantly better in patients with the classical PDAC RNA subtype compared with those with the basal-like subtype (P = 0.004). The best progression-free survival was observed in those with classical subtype treated with m-FOLFIRINOX. GATA6 expression in tumor measured by RNA in situ hybridization was found to be a robust surrogate biomarker for differentiating classical and basal-like PDAC subtypes. Potentially actionable genetic alterations were found in 30% of patients.Conclusions: Prospective genomic profiling of advanced PDAC is feasible, and our early data indicate that chemotherapy response differs among patients with different genomic/transcriptomic subtypes. Clin Cancer Res; 24(6); 1344-54. ©2017 AACR.

Simple multiplexed PCR-based barcoding of DNA for ultrasensitive mutation detection by next-generation sequencing.

Detection of extremely rare variant alleles within a complex mixture of DNA molecules is becoming increasingly relevant in many areas of clinical and basic research, such as the detection of circulating tumor DNA in the plasma of cancer patients. Barcoding of DNA template molecules early in next-generation sequencing (NGS) library construction provides a way to identify and bioinformatically remove polymerase errors that otherwise make detection of these rare variants very difficult. Several barcoding strategies have been reported, but all require long and complex library preparation protocols. Simple, multiplexed, PCR-based barcoding of DNA for sensitive mutation detection using sequencing (SiMSen-seq) was developed to generate targeted barcoded libraries with minimal DNA input, flexible target selection and a very simple, short (∼4 h) library construction protocol. The protocol comprises a three-cycle barcoding PCR step followed directly by adaptor PCR to generate the library and then bead purification before sequencing. Thus, SiMSen-seq allows detection of variant alleles at <0.1% frequency with easy customization of library content (from 1 to 40+ PCR amplicons) and a protocol that can be implemented in any molecular biology laboratory. Here, we provide a detailed protocol for assay development and describe software to process the barcoded sequence reads.

Interrogation of Functional Cell-Surface Markers Identifies CD151 Dependency in High-Grade Serous Ovarian Cancer.

The degree of genetic aberrations characteristic of high-grade serous ovarian cancer (HGSC) makes identification of the molecular features that drive tumor progression difficult. Here, we perform genome-wide RNAi screens and comprehensive expression analysis of cell-surface markers in a panel of HGSC cell lines to identify genes that are critical to their survival. We report that the tetraspanin CD151 contributes to survival of a subset of HGSC cell lines associated with a ZEB transcriptional program and supports the growth of HGSC tumors. Moreover, we show that high CD151 expression is prognostic of poor clinical outcome. This study reveals cell-surface vulnerabilities associated with HGSC, provides a framework for identifying therapeutic targets, and reports a role for CD151 in HGSC.

Simple, multiplexed, PCR-based barcoding of DNA enables sensitive mutation detection in liquid biopsies using sequencing.

Detection of cell-free DNA in liquid biopsies offers great potential for use in non-invasive prenatal testing and as a cancer biomarker. Fetal and tumor DNA fractions however can be extremely low in these samples and ultra-sensitive methods are required for their detection. Here, we report an extremely simple and fast method for introduction of barcodes into DNA libraries made from 5 ng of DNA. Barcoded adapter primers are designed with an oligonucleotide hairpin structure to protect the molecular barcodes during the first rounds of polymerase chain reaction (PCR) and prevent them from participating in mis-priming events. Our approach enables high-level multiplexing and next-generation sequencing library construction with flexible library content. We show that uniform libraries of 1-, 5-, 13- and 31-plex can be generated. Utilizing the barcodes to generate consensus reads for each original DNA molecule reduces background sequencing noise and allows detection of variant alleles below 0.1% frequency in clonal cell line DNA and in cell-free plasma DNA. Thus, our approach bridges the gap between the highly sensitive but specific capabilities of digital PCR, which only allows a limited number of variants to be analyzed, with the broad target capability of next-generation sequencing which traditionally lacks the sensitivity to detect rare variants.

Downregulation of histone H2A and H2B pathways is associated with anthracycline sensitivity in breast cancer.

BACKGROUND: Drug resistance in breast cancer is the major obstacle to effective treatment with chemotherapy. While upregulation of multidrug resistance genes is an important component of drug resistance mechanisms in vitro, their clinical relevance remains to be determined. Therefore, identifying pathways that could be targeted in the clinic to eliminate anthracycline-resistant breast cancer remains a major challenge. METHODS: We generated paired native and epirubicin-resistant MDA-MB-231, MCF7, SKBR3 and ZR-75-1 epirubicin-resistant breast cancer cell lines to identify pathways contributing to anthracycline resistance. Native cell lines were exposed to increasing concentrations of epirubicin until resistant cells were generated. To identify mechanisms driving epirubicin resistance, we used a complementary approach including gene expression analyses to identify molecular pathways involved in resistance, and small-molecule inhibitors to reverse resistance. In addition, we tested its clinical relevance in a BR9601 adjuvant clinical trial. RESULTS: Characterisation of epirubicin-resistant cells revealed that they were cross-resistant to doxorubicin and SN-38 and had alterations in apoptosis and cell-cycle profiles. Gene expression analysis identified deregulation of histone H2A and H2B genes in all four cell lines. Histone deacetylase small-molecule inhibitors reversed resistance and were cytotoxic for epirubicin-resistant cell lines, confirming that histone pathways are associated with epirubicin resistance. Gene expression of a novel 18-gene histone pathway module analysis of the BR9601 adjuvant clinical trial revealed that patients with low expression of the 18-gene histone module benefited from anthracycline treatment more than those with high expression (hazard ratio 0.35, 95 % confidence interval 0.13-0.96, p = 0.042). CONCLUSIONS: This study revealed a key pathway that contributes to anthracycline resistance and established model systems for investigating drug resistance in all four major breast cancer subtypes. As the histone modification can be targeted with small-molecule inhibitors, it represents a possible means of reversing clinical anthracycline resistance. TRIAL REGISTRATION: ClinicalTrials.gov identifier NCT00003012 . Registered on 1 November 1999.