Early Cancer Detection in Li-Fraumeni Syndrome with Cell-Free DNA.

People with Li-Fraumeni syndrome (LFS) harbor a germline pathogenic variant in the TP53 tumor suppressor gene, face a near 100% lifetime risk of cancer, and routinely undergo intensive surveillance protocols. Liquid biopsy has become an attractive tool for a range of clinical applications, including early cancer detection. Here, we provide a proof-of-principle for a multimodal liquid biopsy assay that integrates a targeted gene panel, shallow whole-genome, and cell-free methylated DNA immunoprecipitation sequencing for the early detection of cancer in a longitudinal cohort of 89 LFS patients. Multimodal analysis increased our detection rate in patients with an active cancer diagnosis over uni-modal analysis and was able to detect cancer-associated signal(s) in carriers prior to diagnosis with conventional screening (positive predictive value = 67.6%, negative predictive value = 96.5%). Although adoption of liquid biopsy into current surveillance will require further clinical validation, this study provides a framework for individuals with LFS. SIGNIFICANCE: By utilizing an integrated cell-free DNA approach, liquid biopsy shows earlier detection of cancer in patients with LFS compared with current clinical surveillance methods such as imaging. Liquid biopsy provides improved accessibility and sensitivity, complementing current clinical surveillance methods to provide better care for these patients. See related commentary by Latham et al., p. 23. This article is featured in Selected Articles from This Issue, p. 5.

Integrated, Longitudinal Analysis of Cell-free DNA in Uveal Melanoma.

Uveal melanomas are rare tumors arising from melanocytes that reside in the eye. Despite surgical or radiation treatment, approximately 50% of patients with uveal melanoma will progress to metastatic disease, most often to the liver. Cell-free DNA (cfDNA) sequencing is a promising technology due to the minimally invasive sample collection and ability to infer multiple aspects of tumor response. We analyzed 46 serial cfDNA samples from 11 patients with uveal melanoma over a 1-year period following enucleation or brachytherapy (n = ∼4/patient) using targeted panel, shallow whole genome, and cell-free methylated DNA immunoprecipitation sequencing. We found detection of relapse was highly variable using independent analyses (P = 0.06-0.46), whereas a logistic regression model integrating all cfDNA profiles significantly improved relapse detection (P = 0.02), with greatest power derived from fragmentomic profiles. This work provides support for the use of integrated analyses to improve the sensitivity of circulating tumor DNA detection using multi-modal cfDNA sequencing. Significance: Here, we demonstrate integrated, longitudinal cfDNA sequencing using multi-omic approaches is more effective than unimodal analysis. This approach supports the use of frequent blood testing using comprehensive genomic, fragmentomic, and epigenomic techniques.

Detection of the novel SARS-CoV-2 European lineage B.1.177 in Ontario, Canada.

Background: Travel-related dissemination of SARS-CoV-2 continues to contribute to the global pandemic. A novel SARS-CoV-2 lineage (B.1.177) reportedly arose in Spain in the summer of 2020, with subsequent spread across Europe linked to travel by infected individuals. Surveillance and monitoring through the use of whole genome sequencing (WGS) offers insights into the global and local movement of pathogens such as SARS-CoV-2 and can detect introductions of novel variants. Methods: We analysed the genomes of SARS-CoV-2 sequenced for surveillance purposes from specimens received by Public Health Ontario (Sept 6 - Oct 10, 2020), collected from individuals in eastern Ontario, which comprised the study sample. Taxonomic lineages were identified using pangolin (v2.08) and phylogenetic analysis incorporated publicly available genomes covering the same time period as the study sample. Epidemiological data collected from laboratory requisitions and standard reportable disease case investigation was integrated into the analysis. Results: Genomic surveillance identified a COVID-19 case with SARS-CoV-2 lineage B.1.177 from an individual in eastern Ontario in late September, 2020. The individual had recently returned from Europe. Genomic analysis with publicly available data indicate the most closely related genomes to this specimen were from Southern Europe. Genomic surveillance did not identify further cases with this lineage. Conclusions: Genomic surveillance allowed for early detection of a novel SARS-CoV-2 lineage in Ontario which was deemed to be travel related. This type of genomic-based surveillance is a key tool to measure the effectiveness of public health measures such as mandatory self-isolation for returned travellers, aimed at preventing onward transmission of newly introduced lineages of SARS-CoV-2.

Epigenomic, genomic, and transcriptomic landscape of schwannomatosis.

Schwannomatosis (SWNTS) is a genetic cancer predisposition syndrome that manifests as multiple and often painful neuronal tumors called schwannomas (SWNs). While germline mutations in SMARCB1 or LZTR1, plus somatic mutations in NF2 and loss of heterozygosity in chromosome 22q have been identified in a subset of patients, little is known about the epigenomic and genomic alterations that drive SWNTS-related SWNs (SWNTS-SWNs) in a majority of the cases. We performed multiplatform genomic analysis and established the molecular signature of SWNTS-SWNs. We show that SWNTS-SWNs harbor distinct genomic features relative to the histologically identical non-syndromic sporadic SWNs (NS-SWNS). We demonstrate the existence of four distinct DNA methylation subgroups of SWNTS-SWNs that are associated with specific transcriptional programs and tumor location. We show several novel recurrent non-22q deletions and structural rearrangements. We detected the SH3PXD2A-HTRA1 gene fusion in SWNTS-SWNs, with predominance in LZTR1-mutant tumors. In addition, we identified specific genetic, epigenetic, and actionable transcriptional programs associated with painful SWNTS-SWNs including PIGF, VEGF, MEK, and MTOR pathways, which may be harnessed for management of this syndrome.

Childhood cerebellar tumours mirror conserved fetal transcriptional programs.

Study of the origin and development of cerebellar tumours has been hampered by the complexity and heterogeneity of cerebellar cells that change over the course of development. Here we use single-cell transcriptomics to study more than 60,000 cells from the developing mouse cerebellum and show that different molecular subgroups of childhood cerebellar tumours mirror the transcription of cells from distinct, temporally restricted cerebellar lineages. The Sonic Hedgehog medulloblastoma subgroup transcriptionally mirrors the granule cell hierarchy as expected, while group 3 medulloblastoma resembles Nestin+ stem cells, group 4 medulloblastoma resembles unipolar brush cells, and PFA/PFB ependymoma and cerebellar pilocytic astrocytoma resemble the prenatal gliogenic progenitor cells. Furthermore, single-cell transcriptomics of human childhood cerebellar tumours demonstrates that many bulk tumours contain a mixed population of cells with divergent differentiation. Our data highlight cerebellar tumours as a disorder of early brain development and provide a proximate explanation for the peak incidence of cerebellar tumours in early childhood.

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.

Complex rearrangements and oncogene amplifications revealed by long-read DNA and RNA sequencing of a breast cancer cell line.

The SK-BR-3 cell line is one of the most important models for HER2+ breast cancers, which affect one in five breast cancer patients. SK-BR-3 is known to be highly rearranged, although much of the variation is in complex and repetitive regions that may be underreported. Addressing this, we sequenced SK-BR-3 using long-read single molecule sequencing from Pacific Biosciences and develop one of the most detailed maps of structural variations (SVs) in a cancer genome available, with nearly 20,000 variants present, most of which were missed by short-read sequencing. Surrounding the important ERBB2 oncogene (also known as HER2), we discover a complex sequence of nested duplications and translocations, suggesting a punctuated progression. Full-length transcriptome sequencing further revealed several novel gene fusions within the nested genomic variants. Combining long-read genome and transcriptome sequencing enables an in-depth analysis of how SVs disrupt the genome and sheds new light on the complex mechanisms involved in cancer genome evolution.

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.

De Novo Genome and Transcriptome Assembly of the Canadian Beaver (Castor canadensis).

The Canadian beaver (Castor canadensis) is the largest indigenous rodent in North America. We report a draft annotated assembly of the beaver genome, the first for a large rodent and the first mammalian genome assembled directly from uncorrected and moderate coverage (< 30 ×) long reads generated by single-molecule sequencing. The genome size is 2.7 Gb estimated by k-mer analysis. We assembled the beaver genome using the new Canu assembler optimized for noisy reads. The resulting assembly was refined using Pilon supported by short reads (80 ×) and checked for accuracy by congruency against an independent short read assembly. We scaffolded the assembly using the exon-gene models derived from 9805 full-length open reading frames (FL-ORFs) constructed from the beaver leukocyte and muscle transcriptomes. The final assembly comprised 22,515 contigs with an N50 of 278,680 bp and an N50-scaffold of 317,558 bp. Maximum contig and scaffold lengths were 3.3 and 4.2 Mb, respectively, with a combined scaffold length representing 92% of the estimated genome size. The completeness and accuracy of the scaffold assembly was demonstrated by the precise exon placement for 91.1% of the 9805 assembled FL-ORFs and 83.1% of the BUSCO (Benchmarking Universal Single-Copy Orthologs) gene set used to assess the quality of genome assemblies. Well-represented were genes involved in dentition and enamel deposition, defining characteristics of rodents with which the beaver is well-endowed. The study provides insights for genome assembly and an important genomics resource for Castoridae and rodent evolutionary biology.