The Canadian VirusSeq Data Portal & Duotang: open resources for SARS-CoV-2 viral sequences and genomic epidemiology.

The COVID-19 pandemic led to a large global effort to sequence SARS-CoV-2 genomes from patient samples to track viral evolution and inform public health response. Millions of SARS-CoV-2 genome sequences have been deposited in global public repositories. The Canadian COVID-19 Genomics Network (CanCOGeN - VirusSeq), a consortium tasked with coordinating expanded sequencing of SARS-CoV-2 genomes across Canada early in the pandemic, created the Canadian VirusSeq Data Portal, with associated data pipelines and procedures, to support these efforts. The goal of VirusSeq was to allow open access to Canadian SARS-CoV-2 genomic sequences and enhanced, standardized contextual data that were unavailable in other repositories and that meet FAIR standards (Findable, Accessible, Interoperable and Reusable). In addition, the Portal data submission pipeline contains data quality checking procedures and appropriate acknowledgement of data generators that encourages collaboration. From inception to execution, the portal was developed with a conscientious focus on strong data governance principles and practices. Extensive efforts ensured a commitment to Canadian privacy laws, data security standards, and organizational processes. This Portal has been coupled with other resources like Viral AI and was further leveraged by the Coronavirus Variants Rapid Response Network (CoVaRR-Net) to produce a suite of continually updated analytical tools and notebooks. Here we highlight this Portal, including its contextual data not available elsewhere, and the 'Duotang', a web platform that presents key genomic epidemiology and modeling analyses on circulating and emerging SARS-CoV-2 variants in Canada. Duotang presents dynamic changes in variant composition of SARS-CoV-2 in Canada and by province, estimates variant growth, and displays complementary interactive visualizations, with a text overview of the current situation. The VirusSeq Data Portal and Duotang resources, alongside additional analyses and resources computed from the Portal (COVID-MVP, CoVizu), are all open-source and freely available. Together, they provide an updated picture of SARS-CoV-2 evolution to spur scientific discussions, inform public discourse, and support communication with and within public health authorities. They also serve as a framework for other jurisdictions interested in open, collaborative sequence data sharing and analyses.

Detecting Somatic Mutations Without Matched Normal Samples Using Long Reads.

DNA sequencing of tumours to identify somatic mutations has become a critical tool to guide the type of treatment given to cancer patients. The gold standard for mutation calling is comparing sequencing data from the tumour to a matched normal sample to avoid mis-classifying inherited SNPs as mutations. This procedure works extremely well, but in certain situations only a tumour sample is available. While approaches have been developed to find mutations without a matched normal, they have limited accuracy or require specific types of input data (e.g. ultra-deep sequencing). Here we explore the application of single molecule long read sequencing to calling somatic mutations without matched normal samples. We develop a simple theoretical framework to show how haplotype phasing is an important source of information for determining whether a variant is a somatic mutation. We then use simulations to assess the range of experimental parameters (tumour purity, sequencing depth) where this approach is effective. These ideas are developed into a prototype somatic mutation caller, smrest, and its use is demonstrated on two highly mutated cancer cell lines. Finally, we argue that this approach has potential to measure clinically important biomarkers that are based on the genome-wide distribution of mutations: tumour mutation burden and mutation signatures.

Prevalence of Low-Frequency, Antiviral Resistance Variants in SARS-CoV-2 Isolates in Ontario, Canada, 2020-2023.

Importance: Nirmatrelvir-ritonavir is an oral antiviral medication that improves outcomes in SARS-CoV-2 infections. However, there is concern that antiviral resistance will develop and that these viruses could be selected for after treatment. Objective: To determine the prevalence of low-frequency SARS-CoV-2 variants in patient samples that could be selected for by nirmatrelvir-ritonavir. Design, Setting, and Participants: This retrospective cohort study was conducted at 4 laboratories that serve community hospitals, academic tertiary care centers, and COVID-19 assessment centers in Ontario, Canada. Participants included symptomatic or asymptomatic patients who tested positive for SARS-CoV-2 virus and submitted virus samples for diagnostic testing between March 2020 and January 2023. Exposure: SARS-CoV-2 infection. Main Outcomes and Measures: Samples with sufficient viral load underwent next-generation genome sequencing to identify low-frequency antiviral resistance variants that could not be identified through conventional sequencing. Results: This study included 78 866 clinical samples with next-generation whole-genome sequencing data for SARS-CoV-2. Low-frequency variants in the viral nsp5 gene were identified in 128 isolates (0.16%), and no single variant associated with antiviral resistance was predominate. Conclusions and Relevance: This cohort study of low-frequency variants resistant to nirmatrelvir-ritonavir found that these variants were very rare in samples from patients with SARS-CoV-2, suggesting that selection of these variants by nirmatrelvir-ritonavir following the initiation of treatment may also be rare. Surveillance efforts that involve sequencing of viral isolates should continue to monitor for novel resistance variants as nirmatrelvir-ritonavir is used more broadly.

Coronavirus disease 2019 (COVID-19) outbreak on an in-patient medical unit associated with unrecognized exposures in common areas-Epidemiological and whole-genome sequencing investigation.

OBJECTIVE: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) hospital outbreaks have been common and devastating during the coronavirus disease 2019 (COVID-19) pandemic. Understanding SARS-CoV-2 transmission in these environments is critical for preventing and managing outbreaks. DESIGN: Outbreak investigation through epidemiological mapping and whole-genome sequencing phylogeny. SETTING: Hospital in-patient medical unit outbreak in Toronto, Canada, from November 2020 to January 2021. PARTICIPANTS: The outbreak involved 8 patients and 10 staff and was associated with 3 patient deaths. RESULTS: Patients being cared for in geriatric chairs at the nursing station were at high risk for both acquiring and transmitting SARS-CoV-2 to other patients and staff. Furthermore, given the informal nature of these transmissions, they were not initially recognized, which led to further transmission and missing the opportunity for preventative COVID-19 therapies. CONCLUSIONS: During outbreak prevention and management, the risk of informal patient care settings, such as geriatric chairs, should be considered. During high-risk periods or during outbreaks, efforts should be made to care for patients in their rooms when possible.

A Novel Approach to Detect IDH Point Mutations in Gliomas Using Nanopore Sequencing: Test Validation for the Clinical Laboratory.

The use of standard next-generation sequencing technologies to detect key mutations in IDH genes for glioma diagnosis imposes several challenges, including high capital cost and turnaround delays associated with the need for batch testing. For both glioma testing and testing in other tumor types where highly specific mutation identification is required, the high-throughput nature of next-generation sequencing limits the feasibility of using it as a primary approach in clinical laboratories. We hypothesized that third-generation nanopore sequencing by Oxford Nanopore Technologies has the capability to overcome these limitations. This study aimed to develop and validate a nanopore-based IDH mutation detection assay for clinical practice using glioma formalin-fixed, paraffin-embedded (FFPE) tissue. Glioma FFPE (n = 66) samples with confirmed IDH gene mutational status were sequenced on the MinION device using an amplicon-based approach. All cases were concordant when compared with the reference results. Limit of blank and limit of detection for the variant allele fraction were 1.5% and 3.3%, respectively, at 500× read depth per gene. Total sequencing cost per sample was CAD$50 to CAD$134 with results being available in 9 to 15 hours. These findings demonstrate that nanopore-sequencing technology can be leveraged to develop low-cost, high-performance clinical sequencing-based assays with quick turnaround times to support the detection of targeted mutations in FFPE tumor tissue.

Emergence of a mutation in the nucleocapsid gene of SARS-CoV-2 interferes with PCR detection in Canada.

The emergence of Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV-2) was met with rapid development of robust molecular-based detection assays. Many SARS-CoV-2 molecular tests target multiple genetic regions of the virus to maximize detection and protect against diagnostic escape. Despite the relatively moderate mutational rate of SARS-CoV-2, numerous mutations with known negative impact on diagnostic assays have been identified. In early 2021, we identified four samples positive for SARS-CoV-2 with a nucleocapsid (N) gene drop out on Cepheid Xpert® Xpress SARS-CoV-2 assay. Sequencing revealed a single common mutation in the N gene C29200T. Spatiotemporal analysis showed that the mutation was found in at least six different Canadian provinces from May 2020 until May 2021. Phylogenetic analysis showed that this mutation arose multiple times in Canadian samples and is present in six different variants of interest and of concern. The Cepheid testing platform is commonly used in Canada including in remote regions. As such, the existence of N gene mutation dropouts required further investigation. While commercial SARS-CoV-2 molecular detection assays have contributed immensely to the response effort, many vendors are reluctant to make primer/probe sequences publicly available. Proprietary primer/probe sequences create diagnostic 'blind spots' for global SARS-CoV-2 sequence monitoring and limits the ability to detect and track the presence and prevalence of diagnostic escape mutations. We hope that our industry partners will seriously consider making primer/probe sequences available, so that diagnostic escape mutants can be identified promptly and responded to appropriately to maintain diagnostic accuracy.

Temporal Dynamics and Evolution of SARS-CoV-2 Demonstrate the Necessity of Ongoing Viral Genome Sequencing in Ontario, Canada.

Genome-wide variation in SARS-CoV-2 reveals evolution and transmission dynamics which are critical considerations for disease control and prevention decisions. Here, we review estimates of the genome-wide viral mutation rates, summarize current COVID-19 case load in the province of Ontario, Canada (5 January 2021), and analyze published SARS-CoV-2 genomes from Ontario (collected prior to 24 November 2020) to test for more infectious genetic variants or lineages. The reported mutation rate (∼10-6 nucleotide [nt]-1 cycle-1) for SARS-CoV-2 is typical for coronaviruses. Analysis of published SARS-CoV-2 genomes revealed that the G614 spike protein mutation has dominated infections in Ontario and that SARS-CoV-2 lineages present in Ontario have not differed significantly in their rate of spread. These results suggest that the SARS-CoV-2 population circulating in Ontario has not changed significantly to date. However, ongoing genome monitoring is essential for identification of new variants and lineages that may contribute to increased viral transmission.

Revealing nascent RNA processing dynamics with nano-COP.

During maturation, eukaryotic precursor RNAs undergo processing events including intron splicing, 3'-end cleavage, and polyadenylation. Here we describe nanopore analysis of co-transcriptional processing (nano-COP), a method for probing the timing and patterns of RNA processing. An extension of native elongating transcript sequencing, which quantifies transcription genome-wide through short-read sequencing of nascent RNA 3' ends, nano-COP uses long-read nascent RNA sequencing to observe global patterns of RNA processing. First, nascent RNA is stringently purified through a combination of 4-thiouridine metabolic labeling and cellular fractionation. In contrast to cDNA or short-read-based approaches relying on reverse transcription or amplification, the sample is sequenced directly through nanopores to reveal the native context of nascent RNA. nano-COP identifies both active transcription sites and splice isoforms of single RNA molecules during synthesis, providing insight into patterns of intron removal and the physical coupling between transcription and splicing. The nano-COP protocol yields data within 3 d.

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.

New Twists in Detecting mRNA Modification Dynamics.

Modified nucleotides in mRNA are an essential addition to the standard genetic code of four nucleotides in animals, plants, and their viruses. The emerging field of epitranscriptomics examines nucleotide modifications in mRNA and their impact on gene expression. The low abundance of nucleotide modifications and technical limitations, however, have hampered systematic analysis of their occurrence and functions. Selective chemical and immunological identification of modified nucleotides has revealed global candidate topology maps for many modifications in mRNA, but further technical advances to increase confidence will be necessary. Single-molecule sequencing introduced by Oxford Nanopore now promises to overcome such limitations, and we summarize current progress with a particular focus on the bioinformatic challenges of this novel sequencing technology.

Author Correction: Retrospective evaluation of whole exome and genome mutation calls in 746 cancer samples.

Correction to this paper has been published: https://doi.org/10.1038/s41467-020-20128-w.

Simultaneous profiling of chromatin accessibility and methylation on human cell lines with nanopore sequencing.

Probing epigenetic features on DNA has tremendous potential to advance our understanding of the phased epigenome. In this study, we use nanopore sequencing to evaluate CpG methylation and chromatin accessibility simultaneously on long strands of DNA by applying GpC methyltransferase to exogenously label open chromatin. We performed nanopore sequencing of nucleosome occupancy and methylome (nanoNOMe) on four human cell lines (GM12878, MCF-10A, MCF-7 and MDA-MB-231). The single-molecule resolution allows footprinting of protein and nucleosome binding, and determination of the combinatorial promoter epigenetic signature on individual molecules. Long-read sequencing makes it possible to robustly assign reads to haplotypes, allowing us to generate a fully phased human epigenome, consisting of chromosome-level allele-specific profiles of CpG methylation and chromatin accessibility. We further apply this to a breast cancer model to evaluate differential methylation and accessibility between cancerous and noncancerous cells.

Retrospective evaluation of whole exome and genome mutation calls in 746 cancer samples.

The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) curated consensus somatic mutation calls using whole exome sequencing (WES) and whole genome sequencing (WGS), respectively. Here, as part of the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium, which aggregated whole genome sequencing data from 2,658 cancers across 38 tumour types, we compare WES and WGS side-by-side from 746 TCGA samples, finding that ~80% of mutations overlap in covered exonic regions. We estimate that low variant allele fraction (VAF < 15%) and clonal heterogeneity contribute up to 68% of private WGS mutations and 71% of private WES mutations. We observe that ~30% of private WGS mutations trace to mutations identified by a single variant caller in WES consensus efforts. WGS captures both ~50% more variation in exonic regions and un-observed mutations in loci with variable GC-content. Together, our analysis highlights technological divergences between two reproducible somatic variant detection efforts.

GPU accelerated adaptive banded event alignment for rapid comparative nanopore signal analysis.

BACKGROUND: Nanopore sequencing enables portable, real-time sequencing applications, including point-of-care diagnostics and in-the-field genotyping. Achieving these outcomes requires efficient bioinformatic algorithms for the analysis of raw nanopore signal data. However, comparing raw nanopore signals to a biological reference sequence is a computationally complex task. The dynamic programming algorithm called Adaptive Banded Event Alignment (ABEA) is a crucial step in polishing sequencing data and identifying non-standard nucleotides, such as measuring DNA methylation. Here, we parallelise and optimise an implementation of the ABEA algorithm (termed f5c) to efficiently run on heterogeneous CPU-GPU architectures. RESULTS: By optimising memory, computations and load balancing between CPU and GPU, we demonstrate how f5c can perform ∼3-5 × faster than an optimised version of the original CPU-only implementation of ABEA in the Nanopolish software package. We also show that f5c enables DNA methylation detection on-the-fly using an embedded System on Chip (SoC) equipped with GPUs. CONCLUSIONS: Our work not only demonstrates that complex genomics analyses can be performed on lightweight computing systems, but also benefits High-Performance Computing (HPC). The associated source code for f5c along with GPU optimised ABEA is available at https://github.com/hasindu2008/f5c .

Candidate Cancer Driver Mutations in Distal Regulatory Elements and Long-Range Chromatin Interaction Networks.

A comprehensive catalog of cancer driver mutations is essential for understanding tumorigenesis and developing therapies. Exome-sequencing studies have mapped many protein-coding drivers, yet few non-coding drivers are known because genome-wide discovery is challenging. We developed a driver discovery method, ActiveDriverWGS, and analyzed 120,788 cis-regulatory modules (CRMs) across 1,844 whole tumor genomes from the ICGC-TCGA PCAWG project. We found 30 CRMs with enriched SNVs and indels (FDR < 0.05). These frequently mutated regulatory elements (FMREs) were ubiquitously active in human tissues, showed long-range chromatin interactions and mRNA abundance associations with target genes, and were enriched in motif-rewiring mutations and structural variants. Genomic deletion of one FMRE in human cells caused proliferative deficiencies and transcriptional deregulation of cancer genes CCNB1IP1, CDH1, and CDKN2B, validating observations in FMRE-mutated tumors. Pathway analysis revealed further sub-significant FMREs at cancer genes and processes, indicating an unexplored landscape of infrequent driver mutations in the non-coding genome.

A community effort to create standards for evaluating tumor subclonal reconstruction.

Tumor DNA sequencing data can be interpreted by computational methods that analyze genomic heterogeneity to infer evolutionary dynamics. A growing number of studies have used these approaches to link cancer evolution with clinical progression and response to therapy. Although the inference of tumor phylogenies is rapidly becoming standard practice in cancer genome analyses, standards for evaluating them are lacking. To address this need, we systematically assess methods for reconstructing tumor subclonality. First, we elucidate the main algorithmic problems in subclonal reconstruction and develop quantitative metrics for evaluating them. Then we simulate realistic tumor genomes that harbor all known clonal and subclonal mutation types and processes. Finally, we benchmark 580 tumor reconstructions, varying tumor read depth, tumor type and somatic variant detection. Our analysis provides a baseline for the establishment of gold-standard methods to analyze tumor heterogeneity.

Author Correction: Nanopore native RNA sequencing of a human poly(A) transcriptome.

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

Nanopore native RNA sequencing of a human poly(A) transcriptome.

High-throughput complementary DNA sequencing technologies have advanced our understanding of transcriptome complexity and regulation. However, these methods lose information contained in biological RNA because the copied reads are often short and modifications are not retained. We address these limitations using a native poly(A) RNA sequencing strategy developed by Oxford Nanopore Technologies. Our study generated 9.9 million aligned sequence reads for the human cell line GM12878, using thirty MinION flow cells at six institutions. These native RNA reads had a median length of 771 bases, and a maximum aligned length of over 21,000 bases. Mitochondrial poly(A) reads provided an internal measure of read-length quality. We combined these long nanopore reads with higher accuracy short-reads and annotated GM12878 promoter regions to identify 33,984 plausible RNA isoforms. We describe strategies for assessing 3' poly(A) tail length, base modifications and transcript haplotypes.