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.

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.

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.

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.

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.

CaPSID: a bioinformatics platform for computational pathogen sequence identification in human genomes and transcriptomes.

BACKGROUND: It is now well established that nearly 20% of human cancers are caused by infectious agents, and the list of human oncogenic pathogens will grow in the future for a variety of cancer types. Whole tumor transcriptome and genome sequencing by next-generation sequencing technologies presents an unparalleled opportunity for pathogen detection and discovery in human tissues but requires development of new genome-wide bioinformatics tools. RESULTS: Here we present CaPSID (Computational Pathogen Sequence IDentification), a comprehensive bioinformatics platform for identifying, querying and visualizing both exogenous and endogenous pathogen nucleotide sequences in tumor genomes and transcriptomes. CaPSID includes a scalable, high performance database for data storage and a web application that integrates the genome browser JBrowse. CaPSID also provides useful metrics for sequence analysis of pre-aligned BAM files, such as gene and genome coverage, and is optimized to run efficiently on multiprocessor computers with low memory usage. CONCLUSIONS: To demonstrate the usefulness and efficiency of CaPSID, we carried out a comprehensive analysis of both a simulated dataset and transcriptome samples from ovarian cancer. CaPSID correctly identified all of the human and pathogen sequences in the simulated dataset, while in the ovarian dataset CaPSID's predictions were successfully validated in vitro.

Essential gene profiles in breast, pancreatic, and ovarian cancer cells.

UNLABELLED: Genomic analyses are yielding a host of new information on the multiple genetic abnormalities associated with specific types of cancer. A comprehensive description of cancer-associated genetic abnormalities can improve our ability to classify tumors into clinically relevant subgroups and, on occasion, identify mutant genes that drive the cancer phenotype ("drivers"). More often, though, the functional significance of cancer-associated mutations is difficult to discern. Genome-wide pooled short hairpin RNA (shRNA) screens enable global identification of the genes essential for cancer cell survival and proliferation, providing a "functional genomic" map of human cancer to complement genomic studies. Using a lentiviral shRNA library targeting ~16,000 genes and a newly developed, dynamic scoring approach, we identified essential gene profiles in 72 breast, pancreatic, and ovarian cancer cell lines. Integrating our results with current and future genomic data should facilitate the systematic identification of drivers, unanticipated synthetic lethal relationships, and functional vulnerabilities of these tumor types. SIGNIFICANCE: This study presents a resource of genome-scale, pooled shRNA screens for 72 breast, pancreatic, and ovarian cancer cell lines that will serve as a functional complement to genomics data, facilitate construction of essential gene profiles, help uncover synthetic lethal relationships, and identify uncharacterized genetic vulnerabilities in these tumor types. SIGNIFICANCE: This study presents a resource of genome-scale, pooled shRNA screens for 72 breast, pancreatic, and ovarian cancer cell lines that will serve as a functional complement to genomics data, facilitate construction of essential gene profiles, help uncover synthetic lethal relationships, and identify uncharacterized genetic vulnerabilities in these tumor types.

Computational approaches to discovering noncoding RNA.

New developments are being brought to the field of molecular biology with the mounting evidence that RNA transcripts not translated into protein (noncoding RNAs, ncRNAs) hold a variety of biological functions. Computational discovery of ncRNAs is one of these developments, fueled not only by the urge to characterize these sequences but also by necessity to prioritize ones with the most relevant functions for experimental verification. The heterogeneity in size and mode of activity of ncRNAs is reflected in the corresponding diversity of computational methods for their study. Sequence and structural analysis, conservation across species, and relative position to other genomic elements are being used for ncRNA detection. In addition, the recent development of techniques that allow deep sequencing of cell transcripts either globally or from isolated ncRNA-related material is leading the field toward increased use of such high-throughput data. We expect that imminent breakthroughs will include the classification of newer types of ncRNA and new insights into miRNA and piRNA biology, eventually leading toward the completion of a catalog of all human ncRNAs.

Loss of Tankyrase-mediated destruction of 3BP2 is the underlying pathogenic mechanism of cherubism.

Cherubism is an autosomal-dominant syndrome characterized by inflammatory destructive bony lesions resulting in symmetrical deformities of the facial bones. Cherubism is caused by mutations in Sh3bp2, the gene that encodes the adaptor protein 3BP2. Most identified mutations in 3BP2 lie within the peptide sequence RSPPDG. A mouse model of cherubism develops hyperactive bone-remodeling osteoclasts and systemic inflammation characterized by expansion of the myelomonocytic lineage. The mechanism by which cherubism mutations alter 3BP2 function has remained obscure. Here we show that Tankyrase, a member of the poly(ADP-ribose)polymerase (PARP) family, regulates 3BP2 stability through ADP-ribosylation and subsequent ubiquitylation by the E3-ubiquitin ligase RNF146 in osteoclasts. Cherubism mutations uncouple 3BP2 from Tankyrase-mediated protein destruction, which results in its stabilization and subsequent hyperactivation of the SRC, SYK, and VAV signaling pathways.