Targeted genomic sequencing with probe capture for discovery and surveillance of coronaviruses in bats

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 variant of concern lineages disrupt genomic sequencing of SARS-CoV-2 clinical specimens

Mutations in emerging SARS-CoV-2 lineages can interfere with the laboratory methods used to generate high-quality genome sequences for COVID-19 surveillance. Here, we identify 46 mutations in current variant of concern lineages affecting the widely used laboratory protocols for SARS-CoV-2 genomic sequencing by Freed et al. and the ARTIC network. We provide laboratory data showing how three of these mutations disrupted sequencing of P.1 lineage specimens during a recent outbreak in British Columbia, Canada, and we also demonstrate how we modified the Freed et al. protocol to restore performance.

The Contrasting Role of Nasopharyngeal Angiotensin Converting Enzyme 2 (ACE2) Expression in SARS-CoV-2 Infection: A Cross-Sectional Study of People Tested for COVID-19 in British Columbia

BackgroundAngiotensin converting enzyme 2 (ACE2) serves as the host receptor for SARS-CoV-2, with a critical role in viral infection. We aim to understand population level variation of nasopharyngeal ACE2 expression in people tested for COVID-19 and the relationship between ACE2 expression and SARS-CoV-2 viral RNA load, while adjusting for expression of the complementary protease, Transmembrane serine protease 2 (TMPRSS2), soluble ACE2, age, and biological sex. MethodsA cross-sectional study of n=424 participants aged 1-104 years referred for COVID-19 testing was performed in British Columbia, Canada. Participants who tested negative or positive for COVID-19 were matched by age and biological sex. Viral and host gene expression was measured by quantitative reverse-transcriptase polymerase chain reaction. Bivariate analysis and multiple linear regression were performed to understand the role of nasopharyngeal ACE2 expression in SARS-CoV-2 infection. The ACE2 gene was targeted to measure expression of transmembrane and soluble transcripts. FindingsAnalysis shows no association between age and nasopharyngeal ACE2 expression in those who tested negative for COVID-19 (P=0{middle dot}092). Mean expression of transmembrane (P=1{middle dot}2e-4), soluble ACE2 (P<0{middle dot}0001) and TMPRSS2 (P<0{middle dot}0001) differed between COVID-19-negative and -positive groups. In bivariate analysis of COVID-19-positive participants, expression of transmembrane ACE2 positively correlated with SARS-CoV-2 RNA viral load (P<0{middle dot}0001), expression of soluble ACE2 negatively correlated (P<0{middle dot}0001), and no correlation was found with TMPRSS2 (P=0{middle dot}694). Multivariable analysis showed that the greatest viral RNA loads were observed in participants with high transmembrane ACE2 expression (B=0{middle dot}886, 95%CI:[0{middle dot}596 to 1{middle dot}18]), while expression of soluble ACE2 may protect against high viral RNA load in the upper respiratory tract (B= -0{middle dot}0990, 95%CI:[-0{middle dot}176 to -0{middle dot}0224]). InterpretationNasopharyngeal ACE2 expression plays a dual, contrasting role in SARS-CoV-2 infection of the upper respiratory tract. Transmembrane ACE2 positively correlates, while soluble ACE2 negatively correlates with viral RNA load after adjusting for age, biological sex and expression of TMPRSS2. FundingThis project (COV-55) was funded by Genome British Columbia as part of their COVID-19 rapid response initiative.