Public health emergencies like SARS, MERS, and COVID-19 have prioritized surveillance of zoonotic coronaviruses, resulting in extensive genomic characterization of coronavirus diversity in bats. Sequencing viral genomes directly from animal specimens remains a laboratory challenge, however, and most bat coronaviruses have been characterized solely by PCR amplification of small regions from the best-conserved gene. This has resulted in limited phylogenetic resolution and left viral genetic factors relevant to threat assessment undescribed. In this study, we evaluated whether a technique called hybridization probe capture can achieve more extensive genome recovery from surveillance specimens. Using a custom panel of 20,000 probes, we captured and sequenced coronavirus genomic material in 21 swab specimens collected from bats in the Democratic Republic of the Congo. For 15 of these specimens, probe capture recovered more genome sequence than had been previously generated with standard amplicon sequencing protocols, providing a median 6.1-fold improvement (ranging up to 69.1-fold). Probe capture data also identified five novel alpha- and betacoronaviruses in these specimens, and their full genomes were recovered with additional deep sequencing. Based on these experiences, we discuss how probe capture could be effectively operationalized alongside other sequencing technologies for high-throughput, genomics-based discovery and surveillance of bat coronaviruses.

Mutations in emerging severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lineages can interfere with laboratory methods used to generate viral genome sequences for public health surveillance. We identified 20 mutations that are widespread in variant of concern lineages and affect widely used sequencing protocols by the ARTIC network and Freed et al. Three of these mutations disrupted sequencing of P.1 lineage specimens during a recent outbreak in British Columbia, Canada. We provide laboratory validation of protocol modifications that restored sequencing performance. The study findings indicate that genomic sequencing protocols require immediate updating to address emerging mutations. This work also suggests that routine monitoring and protocol updates will be necessary as SARS-CoV-2 continues to evolve. The bioinformatic and laboratory approaches used here provide guidance for this kind of assay maintenance.

Background: Sequencing viruses in many specimens is hindered by excessive background material from hosts, microbiota, and environmental organisms. Consequently, enrichment of target genomic material is necessary for practical high-throughput viral genome sequencing. Hybridization probes are widely used for enrichment in many fields, but their application to viral sequencing faces a major obstacle: it is difficult to design panels of probe oligo sequences that broadly target many viral taxa due to their rapid evolution, extensive diversity, and genetic hypervariability. To address this challenge, we created ProbeTools, a package of bioinformatic tools for generating effective viral capture panels, and for assessing coverage of target sequences by probe panel designs in silico. In this study, we validated ProbeTools by designing a panel of 3600 probes for subtyping the hypervariable haemagglutinin (HA) and neuraminidase (NA) genome segments of avian-origin influenza A viruses (AIVs). Using in silico assessment of AIV reference sequences and in vitro capture on egg-cultured viral isolates, we demonstrated effective performance by our custom AIV panel and ProbeTools’ suitability for challenging viral probe design applications. Results: Based on ProbeTool’s in silico analysis, our panel provided broadly inclusive coverage of 14,772 HA and 11,967 NA reference sequences. For each reference sequence, we calculated the percentage of nucleotide positions covered by our panel in silico; 90% of HA and NA references sequences had at least 90.8 and 95.1% of their nucleotide positions covered respectively. We also observed effective in vitro capture on a representative collection of 23 egg-cultured AIVs that included isolates from wild birds, poultry, and humans and representatives from all HA and NA subtypes. Forty-two of forty-six HA and NA segments had over 98.3% of their nucleotide positions significantly enriched by our custom panel. These in vitro results were further used to validate ProbeTools’ in silico coverage assessment algorithm; 89.2% of in silico predictions were concordant with in vitro results. Conclusions: ProbeTools generated an effective panel for subtyping AIVs that can be deployed for genomic surveillance, outbreak prevention, and pandemic preparedness. Effective probe design against hypervariable AIV targets also validated ProbeTools’ design and coverage assessment algorithms, demonstrating their suitability for other challenging viral capture applications.

Diverse influenza A viruses (IAVs) circulate in wild birds, including highly pathogenic strains that infect poultry and humans. Consequently, surveillance of IAVs in wild birds is a cornerstone of agricultural biosecurity and pandemic preparedness. Surveillance is traditionally done by testing wild birds directly, but obtaining these specimens is labor intensive, detection rates can be low, and sampling is often biased toward certain avian species. As a result, local incursions of dangerous IAVs are rarely detected before outbreaks begin. Testing environmental specimens from wild bird habitats has been proposed as an alternative surveillance strategy. These specimens are thought to contain diverse IAVs deposited by a broad range of avian hosts, including species that are not typically sampled by surveillance programs. To enable this surveillance strategy, we developed a targeted genomic sequencing method for characterizing IAVs in these challenging environmental specimens. It combines custom hybridization probes, unique molecular index-based library construction, and purpose-built bioinformatic tools, allowing IAV genomic material to be enriched and analyzed with single-fragment resolution. We demonstrated our method on 90 sediment specimens from wetlands around Vancouver, Canada. We recovered 2,312 IAV genome fragments originating from all eight IAV genome segments. Eleven hemagglutinin subtypes and nine neuraminidase subtypes were detected, including H5, the current global surveillance priority. Our results demonstrate that targeted genomic sequencing of environmental specimens from wild bird habitats could become a valuable complement to avian influenza surveillance programs. IMPORTANCE: In this study, we developed genome sequencing tools for characterizing avian influenza viruses in sediment from wild bird habitats. These tools enable an environment-based approach to avian influenza surveillance. This could improve early detection of dangerous strains in local wild birds, allowing poultry producers to better protect their flocks and prevent human exposures to potential pandemic threats. Furthermore, we purposefully developed these methods to contend with viral genomic material that is diluted, fragmented, incomplete, and derived from multiple strains and hosts. These challenges are common to many environmental specimens, making these methods broadly applicable for genomic pathogen surveillance in diverse contexts.

Highly pathogenic avian influenza (HPAI) is a viral disease that causes significant rates of morbidity and mortality in domestic poultry and wild birds, with occasional spillover into mammals, including humans. Beginning in November 2021, Canada experienced its longest and largest outbreak of HPAI in history. A portion of this outbreak (H5N1, clade 2.3.4.4b) occurred in western Canada, specifically in British Columbia (B.C.) and the Yukon, between April 12 and September 11, 2022, which was classified as the “first wave” in this region. Wild birds and mammals identified through passive surveillance and suspect domestic poultry flocks were screened for avian influenza virus (AIV), typed H5 by qPCR, and positive cases were whole genome sequenced. Descriptive epidemiological and phylodynamic analyses were performed to: (1) understand the taxonomic and geographic extent of wild species involved; and (2) examine the origins and probable transmission networks of HPAI viruses introduced into B.C./Yukon by comparing local viruses with those circulating elsewhere in North America. This outbreak included 21 species of wild birds, 2 species of wild mammals, 4 commercial, and 12 domestic small flock infected premises. Canada geese (Branta canadensis) and bald eagles (Haliaeetus leucocephalus) were the most common wild species detected. We demonstrate that north‐south avian migration via the Pacific Flyway is the probable route of multiple incursions into this region. Phylogenetic analysis of the hemagglutinin (HA) segment revealed that the B.C./Yukon viruses detected formed five distinct genetic clusters which were maintained across the whole genome. Although, the genome segments were predominantly Eurasian in origin, NP and PB2 segments from all samples, as well as NS and PB1 segments from Cluster 3, had North American origins. Overall, we demonstrate the utility of genomic epidemiology to inform HPAI transmission dynamics across Western Canada and discuss potential knowledge gaps that exist in passive surveillance strategies for HPAI.