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.

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.

The molecular hallmarks and clinical consequences of tumor hypoxia in prostate cancer

81

Background: Localised prostate cancers are classified into risk-groups using clinical measurements like grade and stage to inform treatment decisions. However, these groupings are imprecise: ~30% of intermediate-risk patients suffer relapse of their disease despite precision image-guided radiotherapy or radical prostatectomy. One reason for this variability in response to treatment is the underlying cellular and molecular heterogeneity of tumours. Prostate tumour cells exist within a microenvironment characterized by gradients of oxygen levels and prostate tumours with low levels of oxygen (hypoxia) have poor clinical outcomes. Methods: Hypoxia was measured using multiple mRNA-based signatures. We examined 548 patients with localised prostate cancer and statistically assessed the association of hypoxia with copy-number alterations (CNAs), single-nucleotide variants (SNVs), genomic rearrangements, focal genomic events ( i.e. kataegis, chromothripsis), telomere length, clinical indices ( i.e. grade, stage) and subclonal architecture. Results: Elevated hypoxia was associated with allelic loss of PTEN, higher rates of chromothripsis and intraductal and cribriform carcinoma (IDC-CA). To translate these findings into a biomarker for prostate cancer precision medicine, we integrated tumour microenvironmental data with genomic and pathological information to stratify patients into distinct prognostic groups. Patients with localized prostate cancer that have polyclonal tumours with elevated hypoxia, allelic loss of PTEN and IDC-CA were at the highest risk of rapid biochemical failure (P = 3.48 x10-3, Logrank test) and metastasis (P = 4.61 x 10-3, Logrank test), even after controlling for T-category, Gleason score and pre-treatment PSA. Conclusions: These data suggest that the aggressiveness of prostate cancers is driven by the interplay of the tumour microenvironment, tumour evolutionary trajectories and its genomic mutational profile.

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.

Abstract P2-10-04: Targeted sequencing in early breast cancer: Identification of novel candidate mutations predictive of anthracycline benefit

Background The use of chemotherapies such as anthracyclines and taxanes have improved overall and disease free survival in breast cancer. For all patients, anthracyclines can have significant toxicities including cardiotoxicity and leukemia. It is therefore essential to select the subset of patients who will receive the optimal overall benefit from anthracycline therapy and to identify molecular pathways driving resistance. To fully understand the impact of mutations in the context of current breast cancer therapy, requires a comprehensive mapping of key molecular events in the context of treatment. We sequenced 101 genes, that were prioritized based on not only gene frequency, but also taking into account the importance of amino acid substitution, type of mutation and network connectivity, in 692 primary tumours to both identify driver genes and pathway cassettes and to understand their clinical significance in response to anthracycline treatment.

Methods We performed targeted sequencing in patients from the BR9601 (n=374) and CCTG MA.5 (n=703) clinical trials. The BR9601 and MA.5 clinical trials examined the effectiveness of combination chemotherapy consisting of CMF (cyclophosphamide, methotrexate and 5-fluorouracil) with or without epirubicin. DNA was extracted, samples were sequenced using AmpliSeq Technology adapted to Illumina and somatic mutations were called using a novel mutation calling pipeline (ISOWN). A priori analyses were performed using distant recurrence free survival (DRFS) as the primary endpoint.

Results: In 692 successfully analysed samples 509 (73.6%) samples exhibited at least one single nucleotide mutation (range 0-54). 94/101 genes were mutated in at least one patient. Only variants in PIK3CA, TP53, CDH1, TLE6, MLL3 and USH2A were detected in 5% or more of samples. TSC22D1, RB1 and ZNF565 were associated with increased risk of distant relapse in multivariate analyses corrected for clinic-pathological variables. No single genes were predictive of anthracycline treatment compared to CMF in multivariate analyses corrected for clinic-pathological variables. Signaling cassettes/modules were designed based on the pathway database, Reactome. Within the signaling cassettes one module was predictive of anthracycline failure. Patients with one or more mutations in this module had an increased risk of distant relapse (HR 0.52, 95% CI 0.29-0.95, p=0.034) when treated with an anthracycline containing chemotherapy regimen compared to CMF (HR 1.34 95% CI 1.05-1.72, p=0.019).

Conclusions: We successfully performed a signaling pathway-based targeted sequencing analysis within predefined signaling modules. We identified a single signaling cassette linked to anthracycline resistance in early breast cancer. However, further work to validate this study in a separate clinical trial is warranted.

Citation Format: Spears M, Kalatskaya I, Trinh QM, Liao L, Chong TM, Crozier C, Dion D, Heisler L, Timms L, Stein LD, Pritchard KI, Levine MN, Shepherd L, Twelves CJ, Bartlett JMS. Targeted sequencing in early breast cancer: Identification of novel candidate mutations predictive of anthracycline benefit [abstract]. In: Proceedings of the 2017 San Antonio Breast Cancer Symposium; 2017 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2018;78(4 Suppl):Abstract nr P2-10-04.

Targeted sequencing in a phase III trial of luminal breast cancer: Identification of novel targets

505

Background: The International Cancer Genome Consortium and The Cancer Genome Atlas have had a global transformative impact on our understanding of cancer. These programs have mapped the genomic landscape of common and rare tumors setting the scene for a comprehensive change in the approach to cancer diagnosis and treatment. However, the task remains incomplete until these mutational events are linked to clinical outcomes in the context of current therapeutic intervention to drive future stratified medicine approaches. Methods: We performed targeted sequencing in patients from the Tamoxifen Exemestane Adjuvant Multicentre trial. DNA was extracted and a 101 gene panel analysed using a novel mutation calling pipeline. Both a priori and machine learning analyses were performed using distant recurrence free survival as the primary endpoint. Results: In 1,491 successfully analyzed samples 1,070 (71.76%) samples exhibited at least one single nucleotide mutation (range 0-94, 1.828+/-0.133, mean+/-s.e.). 98/101 genes were mutated in at least one patient. Only variants in PIK3CA, TP53, MLL3, CDH1 were detected in 5% or more of samples. Twenty genes were associated with increased risk of recurrence in multivariate analyses corrected for clinic-pathological variables, 50% of these genes were involved in transcriptional regulation or RNA/protein processing. In a multivariate analysis, two combined signalling modules were independently prognostic for residual risk following hormone therapy (HRvalidation 3.10 95%CI 1.78-5.40 and HRvalidation 2.70 95%CI 1.57-4.64). Conclusions: We successfully performed a signalling pathway-based targeted sequencing analysis within predefined signalling modules. In supervised and unsupervised analyses we identified multiple signalling cassettes linked to poor outcome in patients with ER+ve breast cancers treated with modern endocrine therapy in the context of a phase III clinical trial. These results identify novel candidates as targets to treat endocrine refractory breast cancers.

Integrated somatic subtypes of localized intermediate-risk prostate cancer

e560

Background: Approximately two thirds of intermediate risk prostate cancer patients are over- or under- treated because we cannot correctly prognose this risk group; therefore we require novel biomarkers to better direct patient therapies and avoid subjecting patients to side effects without benefit. One reason genomic biomarkers are not currently used in clinical settings is because they are notoriously difficult to validate in follow-up studies. Furthermore, the lack of clear prostate cancer subtypes prevents the development of subtype specific biomarkers as is standard practice in breast cancer. We aim to improve biomarker validation rates by defining prostate cancer subtypes that can be used to create subtype specific biomarkers. Methods: First, we assess large scale genomic patterns using whole genome sequencing and methylation data and create integrative subtypes for intermediate risk prostate cancer. Second, we assess associations between specific aberrations and subtypes, and determine whether certain types of molecular aberrations are more important background aberrations for subtype specific biomarker development. Finally, we assess biases in prognostic performance of the published Lalonde biomarker between groups associated with patient subtypes to show that subtype aware biomarkers are necessary. Results: We demonstrate that the Lalonde biomarker is biased by the cohorts’ proportion of TMPRSS2-ERG (T2E) aberrations illustrating the need to develop different biomarkers for patients with T2E and patients without T2E. Further, we suggest integrative subtypes can be used to select patients with similar genomic profiles for biomarker analysis to improve biomarker validation rates. Conclusions: This analysis provides direct guidance for future biomarker development and addresses an important barrier to clinical use of genomic biomarkers for prostate cancer.

A renewed model of pancreatic cancer evolution based on genomic rearrangement patterns

Pancreatic cancer, a highly aggressive tumour type with uniformly poor prognosis, exemplifies the classically held view of stepwise cancer development. The current model of tumorigenesis, based on analyses of precursor lesions, termed pancreatic intraepithelial neoplasm (PanINs) lesions, makes two predictions: first, that pancreatic cancer develops through a particular sequence of genetic alterations (KRAS, followed by CDKN2A, then TP53 and SMAD4); and second, that the evolutionary trajectory of pancreatic cancer progression is gradual because each alteration is acquired independently. A shortcoming of this model is that clonally expanded precursor lesions do not always belong to the tumour lineage, indicating that the evolutionary trajectory of the tumour lineage and precursor lesions can be divergent. This prevailing model of tumorigenesis has contributed to the clinical notion that pancreatic cancer evolves slowly and presents at a late stage. However, the propensity for this disease to rapidly metastasize and the inability to improve patient outcomes, despite efforts aimed at early detection, suggest that pancreatic cancer progression is not gradual. Here, using newly developed informatics tools, we tracked changes in DNA copy number and their associated rearrangements in tumour-enriched genomes and found that pancreatic cancer tumorigenesis is neither gradual nor follows the accepted mutation order. Two-thirds of tumours harbour complex rearrangement patterns associated with mitotic errors, consistent with punctuated equilibrium as the principal evolutionary trajectory. In a subset of cases, the consequence of such errors is the simultaneous, rather than sequential, knockout of canonical preneoplastic genetic drivers that are likely to set-off invasive cancer growth. These findings challenge the current progression model of pancreatic cancer and provide insights into the mutational processes that give rise to these aggressive tumours.

Distinct transcriptional profiles in ex vivo CD4+ and CD8+ T cells are established early in human immunodeficiency virus type 1 infection and are characterized by a chronic interferon response as well as extensive transcriptional changes in CD8+ T cells

Changes in T-cell function are a hallmark of human immunodeficiency virus type 1 (HIV-1) infection, but the pathogenic mechanisms leading to these changes are unclear. We examined the gene expression profiles in ex vivo human CD4+ and CD8+ T cells from untreated HIV-1-infected individuals at different clinical stages and rates of disease progression. Profiles of pure CD4+ and CD8+ T-cell subsets from HIV-1-infected nonprogressors with controlled viremia were indistinguishable from those of individuals not infected with HIV-1. Similarly, no gene clusters could distinguish T cells from individuals with early infection from those seen in chronic progressive HIV-1 infection, whereas differences were observed between uninfected individuals or nonprogressors versus early or chronic progressors. In early and chronic HIV-1 infection, three characteristic gene expression signatures were observed. (i) CD4+ and CD8+ T cells showed increased expression of interferon-stimulated genes (ISGs). However, some ISGs, including CXCL9, CXCL10, and CXCL11, and the interleukin-15 alpha receptor were not upregulated. (ii) CD4+ and CD8+ T cells showed a cluster similar to that observed in thymocytes. (iii) More genes were differentially regulated in CD8+ T cells than in CD4+ T cells, including a cluster of genes downregulated exclusively in CD8+ T cells. In conclusion, HIV-1 infection induces a persistent T-cell transcriptional profile, early in infection, characterized by a dramatic but potentially aberrant interferon response and a profile suggesting an active thymic output. These findings highlight the complexity of the host-virus relationship in HIV-1 infection.

Comparison of detection platforms and post-polymerase chain reaction DNA purification methods for use in conjunction with Cleavase fragment length polymorphism analysis

The removal of impurities and contaminants from PCR-amplified fragments is important for mutation detection methods which identify mutations based on shifts in electrophoretic mobility. This is particularly critical for assays and detection methods which use target DNA that is labeled prior to analysis and electrophoretic detection. We examined several procedures for purifying DNA amplified by the polymerase chain reaction (PCR) and their use in conjunction with a novel DNA scanning method, the Cleavase fragment length polymorphism (CFLP)* assay. In this study, a 480 bp DNA fragment, fluorescently labeled on the 5'-end of one strand, was amplified and subjected to various widely used purification procedures, including several commercially available clean-up kits. We demonstrate that visualization of the fluorescent label, as opposed to simple ethidium bromide staining, reveals the presence of considerable levels of labeled, truncated, amplification products. The various procedures were evaluated on the basis of their ability to remove these unwanted DNA fragments as well as on the degree to which they inhibited or promoted the CFLP reaction. Several procedures are recommended for use with CFLP analysis, including isopropanol precipitation, gel excision, and several commercially available spin columns. Concurrently, we evaluated (compared) a number of commonly used visualization platforms, including fluorescence imaging, chemiluminescence, and post-electrophoretic staining, for the ability to detect CFLP pattern changes. The advantages and disadvantages of different methods are discussed and amounts of DNA to be used for CFLP analysis on different detection platforms are recommended.

Comparison of detection platforms and post-polymerase chain reaction DNA purification methods for use in conjunction with Cleavase fragment length polymorphism analysis

The removal of impurities and contaminants from PCR-amplified fragments is important for mutation detection methods which identify mutations based on shifts in electrophoretic mobility. This is particularly critical for assays and detection methods which use target DNA that is labeled prior to analysis and electrophoretic detection. We examined several procedures for purifying DNA amplified by the polymerase chain reaction (PCR) and their use in conjunction with a novel DNA scanning method, the Cleavase fragment length polymorphism (CFLP)* assay. In this study, a 480 bp DNA fragment, fluorescently labeled on the 5'-end of one strand, was amplified and subjected to various widely used purification procedures, including several commercially available clean-up kits. We demonstrate that visualization of the fluorescent label, as opposed to simple ethidium bromide staining, reveals the presence of considerable levels of labeled, truncated, amplification products. The various procedures were evaluated on the basis of their ability to remove these unwanted DNA fragments as well as on the degree to which they inhibited or promoted the CFLP reaction. Several procedures are recommended for use with CFLP analysis, including isopropanol precipitation, gel excision, and several commercially available spin columns. Concurrently, we evaluated (compared) a number of commonly used visualization platforms, including fluorescence imaging, chemiluminescence, and post-electrophoretic staining, for the ability to detect CFLP pattern changes. The advantages and disadvantages of different methods are discussed and amounts of DNA to be used for CFLP analysis on different detection platforms are recommended.

Detection of p53 gene mutations: analysis by single-strand conformation polymorphism and Cleavase fragment length polymorphism

We have generated a collection of clones containing single point mutations within the exon 5-9 hot spot regions of the p53 gene by using polymerase chain reaction (PCR) to amplify select regions of the gene from characterized cell lines. These clones were then used to address the sensitivity of mutation detection using slab-gel single-strand conformation polymorphism (SSCP) and Cleavase fragment length polymorphism (CFLP) assay systems. Both methods exhibited high sensitivities for the detection of mutations in cloned p53 mutations in this study: 97% for CFLP and 94% for SSCP. In addition to resulting in higher sensitivity of mutation detection, CFLP has the capability to analyze longer fragments. In this study, CFLP identified five intronic mutations which were not investigated in the exon-specific SSCP assay. These results agree with those found elsewhere and demonstrate that CFLP scanning can have practical advantages when used for the identification of sequence alterations within the p53 gene.

Inverse Relationship between Polyamine Levels and the Degree of Phenotypic Alteration Induced by the Root-Inducing, Left-Hand Transferred DNA from Agrobacterium rhizogenes

Floral induction in plants is a paradigm for signal perception, transduction, and physiological response. The introduction of root-inducing, left-hand transferred DNA (Ri T-DNA) into the genomes of several plants results in modifications of flowering (D Tepfer [1984] Cell 47: 959-967), including a delay in flowering in tobacco (Nicotiana tabacum). Conjugated polyamines are markers for flowering in numerous species of plants. In tobacco their accumulation is correlated with the onset of flowering (F Cabanne et al. [1981] Physiol Plant 53: 399-404). Using tobacco, we have explored the possibility of a correlation between the expression of Ri TL-DNA and changes in polyamine metabolism. We made use of two levels of phenotypic change, designated T and T', that retard flowering by 5 to 10 and 15 to 20 days, respectively. We show that delay in flowering is correlated with a reduction in polyamine accumulation and with a delay in appearance of conjugated polyamines, and we propose that genes carried by the Ri TL-DNA intervene either directly in polyamine metabolism or that polyamine metabolism is closely linked to direct effects of Ri T-DNA expression.

Asymptomatic gonorrhea in the male

In a prospective study, 133 male contacts of 102 women with gonorrhea were examined. A diagnosis of gonorrhea was made in 63 patients, indicating an infection rate of 47% after exposure. Forty-three percent of the patients with gonorrhea were asymptomatic. The duration of the asymptomatic period ranged from 3 to 154 days. The epidemiologic importance of this hidden infectious reservoir is stressed.