Enhanced methods for unbiased deep sequencing of Lassa and Ebola RNA viruses from clinical and biological samples

Zanamivir and Baloxavir Combination to cure Persistent Influenza and Coronavirus Infections after Hematopoietic Stem Cell Transplant

OBJECTIVES

. Immunocompromised patients may experience prolonged shedding of influenza virus potentially leading to severe infections. Alternatives to monotherapy with neuraminidase inhibitors should be evaluated to entirely suppress viral replication and prevent drug-resistant mutations.

METHODS

. We investigated the clinical and virological evolution in a case of persistent influenza A and human coronavirus OC43 (HCoV-OC43) coinfection in a hematopoietic stem cell transplant recipient after different therapeutic strategies.

RESULTS

. Successive oseltamivir and zanamivir monotherapies failed to control both infections, with positive results persisting for over 110 days each. This led to the emergence of highly resistant oseltamivir strains due to neuraminidase mutations (E119V and R292K) followed by a deletion (del245-248), while maintaining sensitivity to zanamivir. The intra-host viral diversity data showed that the treatments impacted viral diversity of influenza virus, but not of HCoV-OC43. Considering the patient's underlying condition and the impact of prolonged viral shedding on pulmonary function, eradicating the influenza virus was necessary. A 10-day regimen combining zanamivir and baloxavir-marboxil effectively controlled influenza virus replication and was associated with the clearance of HCoV-OC43, finally resulting in comprehensive respiratory recovery.

CONCLUSION

These observations underscore the importance of further investigating combination treatments as the primary approach to achieve influenza eradication in immunocompromised patients.

Distinct Escherichia coli transcriptional profiles in the guts of recurrent UTI sufferers revealed by pangenome hybrid selection

Low-abundance members of microbial communities are difficult to study in their native habitats. This includes Escherichia coli, a minor, but common inhabitant of the gastrointestinal tract and opportunistic pathogen, including of the urinary tract, where it is the primary pathogen. While multi-omic analyses have detailed critical interactions between uropathogenic Escherichia coli (UPEC) and the bladder that mediate UTI outcome, comparatively little is known about UPEC in its pre-infection reservoir, partly due to its low abundance there (<1% relative abundance). To accurately and sensitively explore the genomes and transcriptomes of diverse E. coli in gastrointestinal communities, we developed E. coli PanSelect which uses a set of probes designed to specifically recognize and capture E. coli's broad pangenome from sequencing libraries. We demonstrated the ability of E. coli PanSelect to enrich, by orders of magnitude, sequencing data from diverse E. coli using a mock community and a set of human stool samples collected as part of a cohort study investigating drivers of recurrent urinary tract infections (rUTI). Comparisons of genomes and transcriptomes between E. coli residing in the gastrointestinal tracts of women with and without a history of rUTI suggest that rUTI gut E. coli are responding to increased levels of oxygen and nitrate, suggestive of mucosal inflammation, which may have implications for recurrent disease. E. coli PanSelect is well suited for investigations of native in vivo biology of E. coli in other environments where it is at low relative abundance, and the framework described here has broad applicability to other highly diverse, low abundance organisms.

A proof of concept for a targeted enrichment approach to the simultaneous detection and characterization of rickettsial pathogens from clinical specimens

Infection with either Rickettsia prowazekii or Orientia tsutsugamushi is common, yet diagnostic capabilities are limited due to the short window for positive identification. Until now, although targeted enrichment had been applied to increase sensitivity of sequencing-based detection for various microorganisms, it had not been applied to sequencing of R. prowazekii in clinical samples. Additionally, hybridization-based targeted enrichment strategies had only scarcely been applied to qPCR of any pathogens in clinical samples. Therefore, we tested a targeted enrichment technique as a proof of concept and found that it dramatically reduced the limits of detection of these organisms by both qPCR and high throughput sequencing. The enrichment methodology was first tested in contrived clinical samples with known spiked-in concentrations of R. prowazekii and O. tsutsugamushi DNA. This method was also evaluated using clinical samples, resulting in the simultaneous identification and characterization of O. tsutsugamushi directly from clinical specimens taken from sepsis patients. We demonstrated that the targeted enrichment technique is helpful by lowering the limit of detection, not only when applied to sequencing, but also when applied to qPCR, suggesting the technique could be applied more broadly to include other assays and/or microbes for which there are limited diagnostic or detection modalities.

Hybrid-Capture Target Enrichment in Human Pathogens: Identification, Evolution, Biosurveillance, and Genomic Epidemiology

High-throughput sequencing (HTS) has revolutionised the field of pathogen genomics, enabling the direct recovery of pathogen genomes from clinical and environmental samples. However, pathogen nucleic acids are often overwhelmed by those of the host, requiring deep metagenomic sequencing to recover sufficient sequences for downstream analyses (e.g., identification and genome characterisation). To circumvent this, hybrid-capture target enrichment (HC) is able to enrich pathogen nucleic acids across multiple scales of divergences and taxa, depending on the panel used. In this review, we outline the applications of HC in human pathogens-bacteria, fungi, parasites and viruses-including identification, genomic epidemiology, antimicrobial resistance genotyping, and evolution. Importantly, we explored the applicability of HC to clinical metagenomics, which ultimately requires more work before it is a reliable and accurate tool for clinical diagnosis. Relatedly, the utility of HC was exemplified by COVID-19, which was used as a case study to illustrate the maturity of HC for recovering pathogen sequences. As we unravel the origins of COVID-19, zoonoses remain more relevant than ever. Therefore, the role of HC in biosurveillance studies is also highlighted in this review, which is critical in preparing us for the next pandemic. We also found that while HC is a popular tool to study viruses, it remains underutilised in parasites and fungi and, to a lesser extent, bacteria. Finally, weevaluated the future of HC with respect to bait design in the eukaryotic groups and the prospect of combining HC with long-read HTS.

An amplicon-based nanopore sequencing workflow for rapid tracking of avian influenza outbreaks, France, 2020-2022

During the recent avian influenza epizootics that occurred in France in 2020/21 and 2021/22, the virus was so contagiousness that it was impossible to control its spread between farms. The preventive slaughter of millions of birds consequently was the only solution available. In an effort to better understand the spread of avian influenza viruses (AIVs) in a rapid and innovative manner, we established an amplicon-based MinION sequencing workflow for the rapid genetic typing of circulating AIV strains. An amplicon-based MinION sequencing workflow based on a set of PCR primers targeting primarily the hemagglutinin gene but also the entire influenza virus genome was developed. Thirty field samples from H5 HPAIV outbreaks in France, including environmental samples, were sequenced using the MinION MK1C. A real-time alignment of the sequences with MinKNOW software allowed the sequencing run to be stopped as soon as enough data were generated. The consensus sequences were then generated and a phylogenetic analysis was conducted to establish links between the outbreaks. The whole sequence of the hemagglutinin gene was obtained for the 30 clinical samples of H5Nx HPAIV belonging to clade 2.3.4.4b. The consensus sequences comparison and the phylogenetic analysis demonstrated links between some outbreaks. While several studies have shown the advantages of MinION for avian influenza virus sequencing, this workflow has been applied exclusively to clinical field samples, without any amplification step on cell cultures or embryonated eggs. As this type of testing pipeline requires only a short amount of time to link outbreaks or demonstrate a new introduction, it could be applied to the real-time management of viral epizootics.

Natural history of Ebola virus disease in rhesus monkeys shows viral variant emergence dynamics and tissue-specific host responses

Ebola virus (EBOV) causes Ebola virus disease (EVD), marked by severe hemorrhagic fever; however, the mechanisms underlying the disease remain unclear. To assess the molecular basis of EVD across time, we performed RNA sequencing on 17 tissues from a natural history study of 21 rhesus monkeys, developing new methods to characterize host-pathogen dynamics. We identified alterations in host gene expression with previously unknown tissue-specific changes, including downregulation of genes related to tissue connectivity. EBOV was widely disseminated throughout the body; using a new, broadly applicable deconvolution method, we found that viral load correlated with increased monocyte presence. Patterns of viral variation between tissues differentiated primary infections from compartmentalized infections, and several variants impacted viral fitness in a EBOV/Kikwit minigenome system, suggesting that functionally significant variants can emerge during early infection. This comprehensive portrait of host-pathogen dynamics in EVD illuminates new features of pathogenesis and establishes resources to study other emerging pathogens.

TMPRSS2 is a functional receptor for human coronavirus HKU1

Four endemic seasonal human coronaviruses causing common colds circulate worldwide: HKU1, 229E, NL63 and OC43 (ref. 1). After binding to cellular receptors, coronavirus spike proteins are primed for fusion by transmembrane serine protease 2 (TMPRSS2) or endosomal cathepsins2-9. NL63 uses angiotensin-converting enzyme 2 as a receptor10, whereas 229E uses human aminopeptidase-N11. HKU1 and OC43 spikes bind cells through 9-O-acetylated sialic acid, but their protein receptors remain unknown12. Here we show that TMPRSS2 is a functional receptor for HKU1. TMPRSS2 triggers HKU1 spike-mediated cell-cell fusion and pseudovirus infection. Catalytically inactive TMPRSS2 mutants do not cleave HKU1 spike but allow pseudovirus infection. Furthermore, TMPRSS2 binds with high affinity to the HKU1 receptor binding domain (Kd 334 and 137 nM for HKU1A and HKU1B genotypes) but not to SARS-CoV-2. Conserved amino acids in the HKU1 receptor binding domain are essential for binding to TMPRSS2 and pseudovirus infection. Newly designed anti-TMPRSS2 nanobodies potently inhibit HKU1 spike attachment to TMPRSS2, fusion and pseudovirus infection. The nanobodies also reduce infection of primary human bronchial cells by an authentic HKU1 virus. Our findings illustrate the various evolution strategies of coronaviruses, which use TMPRSS2 to either directly bind to target cells or prime their spike for membrane fusion and entry.

Pan-viral ORFs discovery using Massively Parallel Ribosome Profiling

Unveiling the complete proteome of viruses is crucial to our understanding of the viral life cycle and interaction with the host. We developed Massively Parallel Ribosome Profiling (MPRP) to experimentally determine open reading frames (ORFs) in 20,170 designed oligonucleotides across 679 human-associated viral genomes. We identified 5,381 ORFs, including 4,208 non-canonical ORFs, and show successful detection of both annotated coding sequences (CDSs) and reported non-canonical ORFs. By examining immunopeptidome datasets of infected cells, we found class I human leukocyte antigen (HLA-I) peptides originating from non-canonical ORFs identified through MPRP. By inspecting ribosome occupancies on the 5'UTR and CDS regions of annotated viral genes, we identified hundreds of upstream ORFs (uORFs) that negatively regulate the synthesis of canonical viral proteins. The unprecedented source of viral ORFs across a wide range of viral families, including highly pathogenic viruses, expands the repertoire of vaccine targets and exposes new cis-regulatory sequences in viral genomes.

Survey of white-footed mice in Connecticut, USA reveals low SARS-CoV-2 seroprevalence and infection with divergent betacoronaviruses

Diverse mammalian species display susceptibility to and infection with SARS-CoV-2. Potential SARS-CoV-2 spillback into rodents is understudied despite their host role for numerous zoonoses and human proximity. We assessed exposure and infection among white-footed mice (Peromyscus leucopus) in Connecticut, USA. We observed 1% (6/540) wild-type neutralizing antibody seroprevalence among 2020-2022 residential mice with no cross-neutralization of variants. We detected no SARS-CoV-2 infections via RT-qPCR, but identified non-SARS-CoV-2 betacoronavirus infections via pan-coronavirus PCR among 1% (5/468) of residential mice. Sequencing revealed two divergent betacoronaviruses, preliminarily named Peromyscus coronavirus-1 and -2. Both belong to the Betacoronavirus 1 species and are ~90% identical to the closest known relative, Porcine hemagglutinating encephalomyelitis virus. Low SARS-CoV-2 seroprevalence suggests white-footed mice may not be sufficiently susceptible or exposed to SARS-CoV-2 to present a long-term human health risk. However, the discovery of divergent, non-SARS-CoV-2 betacoronaviruses expands the diversity of known rodent coronaviruses and further investigation is required to understand their transmission extent.

Advancing Precision Medicine in Myocarditis: Current Status and Future Perspectives in Endomyocardial Biopsy-Based Diagnostics and Therapeutic Approaches

The diagnosis and specific and causal treatment of myocarditis and inflammatory cardiomyopathy remain a major clinical challenge. Despite the rapid development of new imaging techniques, endomyocardial biopsies remain the gold standard for accurate diagnosis of inflammatory myocardial disease. With the introduction and continued development of immunohistochemical inflammation diagnostics in combination with viral nucleic acid testing, myocarditis diagnostics have improved significantly since their introduction. Together with new technologies such as miRNA and gene expression profiling, quantification of specific immune cell markers, and determination of viral activity, diagnostic accuracy and patient prognosis will continue to improve in the future. In this review, we summarize the current knowledge on the pathogenesis and diagnosis of myocarditis and inflammatory cardiomyopathies and highlight future perspectives for more in-depth and specialized biopsy diagnostics and precision, personalized medicine approaches.

Single-cell profiling of lncRNA expression during Ebola virus infection in rhesus macaques

Long non-coding RNAs (lncRNAs) are involved in numerous biological processes and are pivotal mediators of the immune response, yet little is known about their properties at the single-cell level. Here, we generate a multi-tissue bulk RNAseq dataset from Ebola virus (EBOV) infected and not-infected rhesus macaques and identified 3979 novel lncRNAs. To profile lncRNA expression dynamics in immune circulating single-cells during EBOV infection, we design a metric, Upsilon, to estimate cell-type specificity. Our analysis reveals that lncRNAs are expressed in fewer cells than protein-coding genes, but they are not expressed at lower levels nor are they more cell-type specific when expressed in the same number of cells. In addition, we observe that lncRNAs exhibit similar changes in expression patterns to those of protein-coding genes during EBOV infection, and are often co-expressed with known immune regulators. A few lncRNAs change expression specifically upon EBOV entry in the cell. This study sheds light on the differential features of lncRNAs and protein-coding genes and paves the way for future single-cell lncRNA studies.

Towards a Rapid-Turnaround Low-Depth Unbiased Metagenomics Sequencing Workflow on the Illumina Platforms

Unbiased metagenomic sequencing is conceptually well-suited for first-line diagnosis as all known and unknown infectious entities can be detected, but costs, turnaround time and human background reads in complex biofluids, such as plasma, hinder widespread deployment. Separate preparations of DNA and RNA also increases costs. In this study, we developed a rapid unbiased metagenomics next-generation sequencing (mNGS) workflow with a human background depletion method (HostEL) and a combined DNA/RNA library preparation kit (AmpRE) to address this issue. We enriched and detected bacterial and fungal standards spiked in plasma at physiological levels with low-depth sequencing (<1 million reads) for analytical validation. Clinical validation also showed 93% of plasma samples agreed with the clinical diagnostic test results when the diagnostic qPCR had a Ct < 33. The effect of different sequencing times was evaluated with the 19 h iSeq 100 paired end run, a more clinically palatable simulated iSeq 100 truncated run and the rapid 7 h MiniSeq platform. Our results demonstrate the ability to detect both DNA and RNA pathogens with low-depth sequencing and that iSeq 100 and MiniSeq platforms are compatible with unbiased low-depth metagenomics identification with the HostEL and AmpRE workflow.

Phylogeographic reconstruction of the emergence and spread of Powassan virus in the northeastern United States

Powassan virus is an emerging tick-borne virus of concern for public health, but very little is known about its transmission patterns and ecology. Here, we expanded the genomic dataset by sequencing 279 Powassan viruses isolated from Ixodes scapularis ticks from the northeastern United States. Our phylogeographic reconstructions revealed that Powassan virus lineage II was likely introduced or emerged from a relict population in the Northeast between 1940 and 1975. Sequences strongly clustered by sampling location, suggesting a highly focal geographical distribution. Our analyses further indicated that Powassan virus lineage II emerged in the northeastern United States mostly following a south-to-north pattern, with a weighted lineage dispersal velocity of ~3 km/y. Since the emergence in the Northeast, we found an overall increase in the effective population size of Powassan virus lineage II, but with growth stagnating during recent years. The cascading effect of population expansion of white-tailed deer and I. scapularis populations likely facilitated the emergence of Powassan virus in the northeastern United States.

High-depth sequencing characterization of viral dynamics across tissues in fatal COVID-19 reveals compartmentalized infection

SARS-CoV-2 distribution and circulation dynamics are not well understood due to challenges in assessing genomic data from tissue samples. We develop experimental and computational workflows for high-depth viral sequencing and high-resolution genomic analyses from formalin-fixed, paraffin-embedded tissues and apply them to 120 specimens from six subjects with fatal COVID-19. To varying degrees, viral RNA is present in extrapulmonary tissues from all subjects. The majority of the 180 viral variants identified within subjects are unique to individual tissue samples. We find more high-frequency (>10%) minor variants in subjects with a longer disease course, with one subject harboring ten such variants, exclusively in extrapulmonary tissues. One tissue-specific high-frequency variant was a nonsynonymous mutation in the furin-cleavage site of the spike protein. Our findings suggest adaptation and/or compartmentalized infection, illuminating the basis of extrapulmonary COVID-19 symptoms and potential for viral reservoirs, and have broad utility for investigating human pathogens.

Genetic characterization and pathogenicity of a Eurasian avian-like H1N1 swine influenza reassortant virus

BACKGROUND

Swine influenza viruses (SIV), considered the "mixing vessels" of influenza viruses, posed a significant threat to global health systems and are dangerous pathogens. Eurasian avian-like H1N1(EA-H1N1) viruses have become predominant in swine populations in China since 2016.

METHODS

Lung tissue samples were obtained from pregnant sows with miscarriage and respiratory disease in Heilongjiang province, and pathogens were detected by Next-generation sequencing (NGS) and PCR. The nucleic acid of isolates was extracted to detect SIV by RT-PCR. Then, SIV-positive samples were inoculated into embryonated chicken eggs. After successive generations, the isolates were identified by RT-PCR, IFA, WB and TEM. The genetic evolution and pathogenicity to mice of A/swine/Heilongjiang/GN/2020 were analyzed.

RESULTS

The major pathogens were influenza virus (31%), Simbu orthobunyavirus (15%) and Jingmen tick virus (8%) by NGS, while the pathogen that can cause miscarriage and respiratory disease was influenza virus. The SIV(A/swine/Heilongjiang/GN/2020) with hemagglutination activity was isolated from lung samples and was successfully identified by RT-PCR, IFA, WB and TEM. Homology and phylogenetic analysis showed that A/swine/Heilongjiang/GN/2020 is most closely related to A/swine/Henan/SN/10/2018 and belonged to EA-H1N1. Pathogenicity in mice showed that the EA-H1N1 could cause lethal or exhibit extrapulmonary virus spread and cause severe damage to respiratory tracts effectively proliferating in lung and trachea.

CONCLUSION

A/swine/Heilongjiang/GN/2020 (EA-H1N1) virus was isolated from pregnant sows with miscarriage and respiratory disease in Heilongjiang province, China. Clinical signs associated with influenza infection were observed during 14 days with A/swine/Heilongjiang/GN/2020 infected mice. These data suggest that A/swine/Heilongjiang/GN/2020 (EA-H1N1) had high pathogenicity and could be systemic spread in mice.

Next Generation Sequencing Approaches to Characterize the Respiratory Tract Virome

The COVID-19 pandemic and heightened perception of the risk of emerging viral infections have boosted the efforts to better understand the virome or complete repertoire of viruses in health and disease, with a focus on infectious respiratory diseases. Next-generation sequencing (NGS) is widely used to study microorganisms, allowing the elucidation of bacteria and viruses inhabiting different body systems and identifying new pathogens. However, NGS studies suffer from a lack of standardization, in particular, due to various methodological approaches and no single format for processing the results. Here, we review the main methodological approaches and key stages for studies of the human virome, with an emphasis on virome changes during acute respiratory viral infection, with applications for clinical diagnostics and epidemiologic analyses.

Improvements in metagenomic virus detection by simple pretreatment methods

Early detection of pathogens at the point of care helps reduce the threats to human and animal health from emerging pathogens. Initially, the disease-causing agent will be unknown and needs to be identified; this often requires specific laboratory facilities. Here we describe the development of an unbiased detection assay for RNA and DNA viruses using metagenomic Nanopore sequencing and simple methods that can be transferred into a field setting. Human clinical samples containing the RNA virus SARS-CoV-2 or the DNA viruses human papillomavirus (HPV) and molluscum contagiosum virus (MCV) were used as a test of concept. Firstly, the virus detection potential was optimized by investigating different pretreatments for reducing non-viral nucleic acid components. DNase I pretreatment followed by filtration increased the proportion of SARS-CoV-2 sequenced reads > 500-fold compared with no pretreatments. This was sufficient to achieve virus detection with high confidence and allowed variant identification. Next, we tested individual SARS-CoV-2 samples with various viral loads (measured as CT-values determined by RT-qPCR). Lastly, we tested the assay on clinical samples containing the DNA virus HPV and co-infection with MCV to show the assay's detection potential for DNA viruses. This protocol is fast (same day results). We hope to apply this method in other settings for point of care detection of virus pathogens, thus eliminating the need for transport of infectious samples, cold storage and a specialized laboratory.

Non-target RNA depletion strategy to improve sensitivity of next-generation sequencing for the detection of RNA viruses in poultry

PCR-based assays have become the benchmark for detecting pathogens of poultry and other livestock; however, these techniques are limited in their ability to detect multiple infecting agents, provide limited genetic information on the pathogen, and, for RNA viruses, must be reviewed frequently to assure high sensitivity and specificity. In contrast, untargeted, high-throughput sequencing can rapidly detect all infecting agents in a sample while providing genomic sequence information to allow more in-depth characterization of viruses. Although next-generation sequencing (NGS) offers many advantages, one of its primary limitations is low sensitivity to pathogens given the abundance of host and other non-target sequences in sequencing libraries. We explored methods for improving the sensitivity of NGS to detect respiratory and enteric viruses in poultry from RNA extracts of swab samples. We employed commercial and custom-designed negative enrichment strategies to selectively deplete the most abundant rRNA reads from the host and non-target bacteria; host RNA was diminished from up to 40% of total reads to as low as 3%, and the total number of reads assigned to abundant bacterial classes were reduced greatly. Our treatment resulted in up to a 700-fold increase in the number of viral reads, detection of a greater number of viral agents, and higher average genome coverage for pathogens. Depletion assays added only 2 h to the NGS library preparation workflow. Custom depletion probe design offered significant cost savings (US$7-12 per sample) compared to commercial kits (US$30-50 per sample).

Multiplexed CRISPR-based microfluidic platform for clinical testing of respiratory viruses and identification of SARS-CoV-2 variants

The coronavirus disease 2019 (COVID-19) pandemic has demonstrated a clear need for high-throughput, multiplexed and sensitive assays for detecting severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and other respiratory viruses and their emerging variants. Here, we present a cost-effective virus and variant detection platform, called microfluidic Combinatorial Arrayed Reactions for Multiplexed Evaluation of Nucleic acids (mCARMEN), which combines CRISPR-based diagnostics and microfluidics with a streamlined workflow for clinical use. We developed the mCARMEN respiratory virus panel to test for up to 21 viruses, including SARS-CoV-2, other coronaviruses and both influenza strains, and demonstrated its diagnostic-grade performance on 525 patient specimens in an academic setting and 166 specimens in a clinical setting. We further developed an mCARMEN panel to enable the identification of 6 SARS-CoV-2 variant lineages, including Delta and Omicron, and evaluated it on 2,088 patient specimens with near-perfect concordance to sequencing-based variant classification. Lastly, we implemented a combined Cas13 and Cas12 approach that enables quantitative measurement of SARS-CoV-2 and influenza A viral copies in samples. The mCARMEN platform enables high-throughput surveillance of multiple viruses and variants simultaneously, enabling rapid detection of SARS-CoV-2 variants.

Application of nanopore sequencing in diagnosis of secondary infections in patients with severe COVID-19

To explore the application value of nanopore sequencing technique in the diagnosis and treatment of secondary infections in patients with severe coronavirus disease 2019 (COVID-19). A total of 77 clinical specimens from 3 patients with severe COVID-19 were collected. After heat inactivation, all samples were subjected to total nucleic acid extraction based on magnetic bead enrichment. The extracted DNA was used for DNA library construction, then nanopore real-time sequencing detection was performed. The sequencing data were subjected to Centrifuge software database species matching and R program differential analysis to obtain potential pathogen identification. Nanopore sequencing results were compared with respiratory pathogen qPCR panel screening and conventional microbiological testing results to verify the effectiveness of nanopore sequencing detection. Nanopore sequencing results showed that positive pathogen were obtained in 44 specimens (57.1%). The potential pathogens identified by nanopore sequencing included , , and , et al. , , were also detected in clinical microbiological culture-based detection; was detected in respiratory pathogen screening qPCR panel; was only detected by the nanopore sequencing technique. Comprehensive considerations with the clinical symptoms, the patient was treated with antibiotics against , and the infection was controlled. Nanopore sequencing may assist the diagnosis and treatment of severe COVID-19 patients through rapid identification of potential pathogens.

Synthetic DNA spike-ins (SDSIs) enable sample tracking and detection of inter-sample contamination in SARS-CoV-2 sequencing workflows

The global spread and continued evolution of SARS-CoV-2 has driven an unprecedented surge in viral genomic surveillance. Amplicon-based sequencing methods provide a sensitive, low-cost and rapid approach but suffer a high potential for contamination, which can undermine laboratory processes and results. This challenge will increase with the expanding global production of sequences across a variety of laboratories for epidemiological and clinical interpretation, as well as for genomic surveillance of emerging diseases in future outbreaks. We present SDSI + AmpSeq, an approach that uses 96 synthetic DNA spike-ins (SDSIs) to track samples and detect inter-sample contamination throughout the sequencing workflow. We apply SDSIs to the ARTIC Consortium's amplicon design, demonstrate their utility and efficiency in a real-time investigation of a suspected hospital cluster of SARS-CoV-2 cases and validate them across 6,676 diagnostic samples at multiple laboratories. We establish that SDSI + AmpSeq provides increased confidence in genomic data by detecting and correcting for relatively common, yet previously unobserved modes of error, including spillover and sample swaps, without impacting genome recovery.

Multicenter assessment of shotgun metagenomics for pathogen detection

BACKGROUND

Shotgun metagenomics has been used clinically for diagnosing infectious diseases. However, most technical assessments have been limited to individual sets of reference standards, experimental workflows, and laboratories.

METHODS

A reference panel and performance metrics were designed and used to examine the performance of shotgun metagenomics at 17 laboratories in a coordinated collaborative study. We comprehensively assessed the reliability, key performance determinants, reproducibility, and quantitative potential.

FINDINGS

Assay performance varied significantly across sites and microbial classes, with a read depth of 20 millions as a generally cost-efficient assay setting. Results of mapped reads by shotgun metagenomics could indicate relative and intra-site (but not absolute or inter-site) microbial abundance.

INTERPRETATION

Assay performance was significantly impacted by the microbial type, the host context, and read depth, which emphasizes the importance of these factors when designing reference reagents and benchmarking studies. Across sites, workflows and platforms, false positive reporting and considerable site/library effects were common challenges to the assay's accuracy and quantifiability. Our study also suggested that laboratory-developed shotgun metagenomics tests for pathogen detection should aim to detect microbes at 500 CFU/mL (or copies/mL) in a clinically relevant host context (10^5 human cells/mL) within a 24h turn-around time, and with an efficient read depth of 20M.

FUNDING

This work was supported by National Science and Technology Major Project of China (2018ZX10102001).

Intrathecal inflammatory responses in the absence of SARS-CoV-2 nucleic acid in the CSF of COVID-19 hospitalized patients

OBJECTIVE

Little is known about CSF profiles in patients with acute COVID-19 infection and neurological symptoms. Here, CSF was tested for SARS-CoV-2 RNA and inflammatory cytokines and chemokines and compared to controls and patients with known neurotropic pathogens.

METHODS

CSF from twenty-seven consecutive patients with COVID-19 and neurological symptoms was assayed for SARS-CoV-2 RNA using quantitative reverse transcription PCR (RT-qPCR) and unbiased metagenomic sequencing. Assays for blood brain barrier (BBB) breakdown (CSF:serum albumin ratio (Q-Alb)), and proinflammatory cytokines and chemokines (IL-6, IL-8, IL-15, IL-16, monocyte chemoattractant protein -1 (MCP-1) and monocyte inhibitory protein - 1β (MIP-1β)) were performed in 23 patients and compared to CSF from patients with HIV-1 (16 virally suppressed, 5 unsuppressed), West Nile virus (WNV) (n = 4) and 16 healthy controls (HC).

RESULTS

Median CSF cell count for COVID-19 patients was 1 white blood cell/μL; two patients were infected with a second pathogen (Neisseria, Cryptococcus neoformans). No CSF samples had detectable SARS-CoV-2 RNA by either detection method. In patients with COVID-19 only, CSF IL-6, IL-8, IL-15, and MIP-1β levels were higher than HC and suppressed HIV (corrected-p < 0.05). MCP-1 and MIP-1β levels were higher, while IL-6, IL-8, IL-15 were similar in COVID-19 compared to WNV patients. Q-Alb correlated with all proinflammatory markers, with IL-6, IL-8, and MIP-1β (r ≥ 0.6, p < 0.01) demonstrating the strongest associations.

CONCLUSIONS

Lack of SARS-CoV-2 RNA in CSF is consistent with pre-existing literature. Evidence of intrathecal proinflammatory markers in a subset of COVID-19 patients with BBB breakdown despite minimal CSF pleocytosis is atypical for neurotropic pathogens.

A novel SARS-CoV-2 related coronavirus in bats from Cambodia

Knowledge of the origin and reservoir of the coronavirus responsible for the ongoing COVID-19 pandemic is still fragmentary. To date, the closest relatives to SARS-CoV-2 have been detected in Rhinolophus bats sampled in the Yunnan province, China. Here we describe the identification of SARS-CoV-2 related coronaviruses in two Rhinolophus shameli bats sampled in Cambodia in 2010. Metagenomic sequencing identifies nearly identical viruses sharing 92.6% nucleotide identity with SARS-CoV-2. Most genomic regions are closely related to SARS-CoV-2, with the exception of a region of the spike, which is not compatible with human ACE2-mediated entry. The discovery of these viruses in a bat species not found in China indicates that SARS-CoV-2 related viruses have a much wider geographic distribution than previously reported, and suggests that Southeast Asia represents a key area to consider for future surveillance for coronaviruses.

In this study, Delaune et al., isolate and characterise a SARS-CoV-2-related coronavirus from two bats sampled in Cambodia. Their findings suggest that the geographic distribution of SARS-CoV-2-related viruses is wider than previously reported.

A novel SARS-CoV-2 related coronavirus in bats from Cambodia

Knowledge of the origin and reservoir of the coronavirus responsible for the ongoing COVID-19 pandemic is still fragmentary. To date, the closest relatives to SARS-CoV-2 have been detected in Rhinolophus bats sampled in the Yunnan province, China. Here we describe the identification of SARS-CoV-2 related coronaviruses in two Rhinolophus shameli bats sampled in Cambodia in 2010. Metagenomic sequencing identifies nearly identical viruses sharing 92.6% nucleotide identity with SARS-CoV-2. Most genomic regions are closely related to SARS-CoV-2, with the exception of a region of the spike, which is not compatible with human ACE2-mediated entry. The discovery of these viruses in a bat species not found in China indicates that SARS-CoV-2 related viruses have a much wider geographic distribution than previously reported, and suggests that Southeast Asia represents a key area to consider for future surveillance for coronaviruses.

Genomic surveillance of enterovirus associated with aseptic meningitis cases in southern Spain, 2015-2018

New circulating Enterovirus (EV) strains often emerge through recombination. Upsurges of recombinant non-polio enteroviruses (NPEVs) associated with neurologic manifestations such as EVA71 or Echovirus 30 (E30) are a growing public health concern in Europe. Only a few complete genomes of EVs circulating in Spain are available in public databases, making it difficult to address the emergence of recombinant EVs, understand their evolutionary relatedness and the possible implication in human disease. We have used metagenomic (untargeted) NGS to generate full-length EV genomes from CSF samples of EV-positive aseptic meningitis cases in Southern Spain between 2015 and 2018. Our analyses reveal the co-circulation of multiple Enterovirus B (EV-B) types (E6, E11, E13 and E30), including a novel E13 recombinant form. We observed a genetic turnover where emergent lineages (C1 for E6 and I [tentatively proposed in this study] for E30) replaced previous lineages circulating in Spain, some concomitant with outbreaks in other parts of Europe. Metagenomic sequencing provides an effective approach for the analysis of EV genomes directly from PCR-positive CSF samples. The detection of a novel, disease-associated, recombinant form emphasizes the importance of genomic surveillance to monitor spread and evolution of EVs.

Molecular archeology of human viruses

The evolution of human-virus associations is usually reconstructed from contemporary patterns of genomic diversity. An intriguing, though still rarely implemented, alternative is to search for the genetic material of viruses in archeological and medical archive specimens to document evolution as it happened. In this chapter, we present lessons from ancient DNA research and incorporate insights from virology to explore the potential range of applications and likely limitations of archeovirological approaches. We also highlight the numerous questions archeovirology will hopefully allow us to tackle in the near future, and the main expected roadblocks to these avenues of research.

Utilizing the VirIdAl Pipeline to Search for Viruses in the Metagenomic Data of Bat Samples

According to various estimates, only a small percentage of existing viruses have been discovered, naturally much less being represented in the genomic databases. High-throughput sequencing technologies develop rapidly, empowering large-scale screening of various biological samples for the presence of pathogen-associated nucleotide sequences, but many organisms are yet to be attributed specific loci for identification. This problem particularly impedes viral screening, due to vast heterogeneity in viral genomes. In this paper, we present a new bioinformatic pipeline, VirIdAl, for detecting and identifying viral pathogens in sequencing data. We also demonstrate the utility of the new software by applying it to viral screening of the feces of bats collected in the Moscow region, which revealed a significant variety of viruses associated with bats, insects, plants, and protozoa. The presence of alpha and beta coronavirus reads, including the MERS-like bat virus, deserves a special mention, as it once again indicates that bats are indeed reservoirs for many viral pathogens. In addition, it was shown that alignment-based methods were unable to identify the taxon for a large proportion of reads, and we additionally applied other approaches, showing that they can further reveal the presence of viral agents in sequencing data. However, the incompleteness of viral databases remains a significant problem in the studies of viral diversity, and therefore necessitates the use of combined approaches, including those based on machine learning methods.

Enhanced Virus Detection and Metagenomic Sequencing in Patients with Meningitis and Encephalitis

Meningitis and encephalitis are leading causes of central nervous system (CNS) disease and often result in severe neurological compromise or death. Traditional diagnostic workflows largely rely on pathogen-specific tests, sometimes over days to weeks, whereas metagenomic next-generation sequencing (mNGS) profiles all nucleic acid in a sample. In this single-center, prospective study, 68 hospitalized patients with known (n = 44) or suspected (n = 24) CNS infections underwent mNGS from RNA and DNA to identify potential pathogens and also targeted sequencing of viruses using hybrid capture. Using a computational metagenomic classification pipeline based on KrakenUniq and BLAST, we detected pathogen nucleic acid in cerebrospinal fluid (CSF) from 22 subjects, 3 of whom had no clinical diagnosis by routine workup. Among subjects diagnosed with infection by serology and/or peripheral samples, we demonstrated the utility of mNGS to detect pathogen nucleic acid in CSF, importantly for the Ixodes scapularis tick-borne pathogens Powassan virus, Borrelia burgdorferi, and Anaplasma phagocytophilum. We also evaluated two methods to enhance the detection of viral nucleic acid, hybrid capture and methylated DNA depletion. Hybrid capture nearly universally increased viral read recovery. Although results for methylated DNA depletion were mixed, it allowed the detection of varicella-zoster virus DNA in two samples that were negative by standard mNGS. Overall, mNGS is a promising approach that can test for multiple pathogens simultaneously, with efficacy similar to that of pathogen-specific tests, and can uncover geographically relevant infectious CNS disease, such as tick-borne infections in New England. With further laboratory and computational enhancements, mNGS may become a mainstay of workup for encephalitis and meningitis. IMPORTANCE Meningitis and encephalitis are leading global causes of central nervous system (CNS) disability and mortality. Current diagnostic workflows remain inefficient, requiring costly pathogen-specific assays and sometimes invasive surgical procedures. Despite intensive diagnostic efforts, 40 to 60% of people with meningitis or encephalitis have no clear cause of CNS disease identified. As diagnostic uncertainty often leads to costly inappropriate therapies, the need for novel pathogen detection methods is paramount. Metagenomic next-generation sequencing (mNGS) offers the unique opportunity to circumvent these challenges using unbiased laboratory and computational methods. Here, we performed comprehensive mNGS from 68 prospectively enrolled patients with known (n = 44) or suspected (n = 24) CNS viral infection from a single center in New England and evaluated enhanced methods to improve the detection of CNS pathogens, including those not traditionally identified in the CNS by nucleic acid detection. Overall, our work helps elucidate how mNGS can become integrated into the diagnostic toolkit for CNS infections.

A single-cell and spatial atlas of autopsy tissues reveals pathology and cellular targets of SARS-CoV-2

The SARS-CoV-2 pandemic has caused over 1 million deaths globally, mostly due to acute lung injury and acute respiratory distress syndrome, or direct complications resulting in multiple-organ failures. Little is known about the host tissue immune and cellular responses associated with COVID-19 infection, symptoms, and lethality. To address this, we collected tissues from 11 organs during the clinical autopsy of 17 individuals who succumbed to COVID-19, resulting in a tissue bank of approximately 420 specimens. We generated comprehensive cellular maps capturing COVID-19 biology related to patients' demise through single-cell and single-nucleus RNA-Seq of lung, kidney, liver and heart tissues, and further contextualized our findings through spatial RNA profiling of distinct lung regions. We developed a computational framework that incorporates removal of ambient RNA and automated cell type annotation to facilitate comparison with other healthy and diseased tissue atlases. In the lung, we uncovered significantly altered transcriptional programs within the epithelial, immune, and stromal compartments and cell intrinsic changes in multiple cell types relative to lung tissue from healthy controls. We observed evidence of: alveolar type 2 (AT2) differentiation replacing depleted alveolar type 1 (AT1) lung epithelial cells, as previously seen in fibrosis; a concomitant increase in myofibroblasts reflective of defective tissue repair; and, putative TP63+ intrapulmonary basal-like progenitor (IPBLP) cells, similar to cells identified in H1N1 influenza, that may serve as an emergency cellular reserve for severely damaged alveoli. Together, these findings suggest the activation and failure of multiple avenues for regeneration of the epithelium in these terminal lungs. SARS-CoV-2 RNA reads were enriched in lung mononuclear phagocytic cells and endothelial cells, and these cells expressed distinct host response transcriptional programs. We corroborated the compositional and transcriptional changes in lung tissue through spatial analysis of RNA profiles in situ and distinguished unique tissue host responses between regions with and without viral RNA, and in COVID-19 donor tissues relative to healthy lung. Finally, we analyzed genetic regions implicated in COVID-19 GWAS with transcriptomic data to implicate specific cell types and genes associated with disease severity. Overall, our COVID-19 cell atlas is a foundational dataset to better understand the biological impact of SARS-CoV-2 infection across the human body and empowers the identification of new therapeutic interventions and prevention strategies.

A streamlined clinical metagenomic sequencing protocol for rapid pathogen identification

Metagenomic next-generation sequencing (mNGS) holds promise as a diagnostic tool for unbiased pathogen identification and precision medicine. However, its medical utility depends largely on assay simplicity and reproducibility. In the current study, we aimed to develop a streamlined Illumina and Oxford Nanopore-based DNA/RNA library preparation protocol and rapid data analysis pipeline. The Illumina sequencing-based mNGS method was first developed and evaluated using a set of samples with known aetiology. Its sensitivity for RNA viruses (influenza A, H1N1) was < 6.4 × 10² EID50/mL, and a good correlation between viral loads and mapped reads was observed. Then, the rapid turnaround time of Nanopore sequencing was tested by sequencing influenza A virus and adenoviruses. Furthermore, 11 respiratory swabs or sputum samples pre-tested for a panel of pathogens were analysed, and the pathogens identified by Illumina sequencing showed 81.8% concordance with qPCR results. Additional sequencing of cerebrospinal fluid (CSF) samples from HIV-1-positive patients with meningitis/encephalitis detected HIV-1 RNA and Toxoplasma gondii sequences. In conclusion, we have developed a simplified protocol that realizes efficient metagenomic sequencing of a variety of clinical samples and pathogen identification in a clinically meaningful time frame.

High Throughput Sequencing for the Detection and Characterization of RNA Viruses

This review aims to assess and recommend approaches for targeted and agnostic High Throughput Sequencing of RNA viruses in a variety of sample matrices. HTS also referred to as deep sequencing, next generation sequencing and third generation sequencing; has much to offer to the field of environmental virology as its increased sequencing depth circumvents issues with cloning environmental isolates for Sanger sequencing. That said however, it is important to consider the challenges and biases that method choice can impart to sequencing results. Here, methodology choices from RNA extraction, reverse transcription to library preparation are compared based on their impact on the detection or characterization of RNA viruses.

Genomic investigation of atypical malaria cases in Kanel, northern Senegal

BACKGROUND

The diagnosis of malaria cases in regions where the malaria burden has decreased significantly and prevalence is very low is more challenging, in part because of reduced clinical presumption of malaria. The appearance of a cluster of malaria cases with atypical symptoms in Mbounguiel, a village in northern Senegal where malaria transmission is low, in September 2018 exemplifies this scenario. The collaboration between the National Malaria Control Programme (NMCP) at the Senegal Ministry of Health and the Laboratory of Parasitology and Mycology at Cheikh Anta Diop University worked together to evaluate this cluster of malaria cases using molecular and serological tools.

METHODS

Malaria cases were diagnosed primarily by rapid diagnostic test (RDT), and confirmed by photo-induced electron transfer-polymerase chain reaction (PET-PCR). 24 single nucleotide polymorphisms (SNPs) barcoding was used for Plasmodium falciparum genotyping. Unbiased metagenomic sequencing and Luminex-based multi-pathogen antibody and antigen profiling were used to assess exposure to other pathogens.

RESULTS

Nine patients, of 15 suspected cases, were evaluated, and all nine samples were found to be positive for P. falciparum only. The 24 SNPs molecular barcode showed the predominance of polygenomic infections, with identifiable strains being different from one another. All patients tested positive for the P. falciparum antigens. No other pathogenic infection was detected by either the serological panel or metagenomic sequencing.

CONCLUSIONS

This work, undertaken locally within Senegal as a collaboration between the NMCP and a research laboratory at University of Cheikh Anta Diop (UCAD) revealed that a cluster of malaria cases were caused by different strains of P. falciparum. The public health response in real time demonstrates the value of local molecular and genomics capacity in affected countries for disease control and elimination.

Genomic surveillance of Nevada patients revealed prevalence of unique SARS-CoV-2 variants bearing mutations in the RdRp gene

Patients with signs of COVID-19 were tested through diagnostic RT-PCR for SARS-CoV-2 using RNA extracted from the nasopharyngeal/nasal swabs. To determine the variants of SARS-CoV-2 circulating in the state of Nevada, specimens from 200 COVID-19 patients were sequenced through our robust sequencing platform, which enabled sequencing of SARS-CoV-2 from specimens with even very low viral loads, without the need of culture-based amplification. High genome coverage allowed the identification of single and multi-nucleotide variants in SARS-CoV-2 in the community and their phylogenetic relationships with other variants present during the same period of the outbreak. We report the occurrence of a novel mutation at 323aa (314aa of orf1b) of nsp12 (RNA-dependent RNA polymerase) changed to phenylalanine (F) from proline (P), in the first reported isolate of SARS-CoV-2, Wuhan-Hu-1. This 323F variant was present at a very high frequency in Northern Nevada. Structural modeling determined this mutation in the interface domain, which is important for the association of accessory proteins required for the polymerase. In conclusion, we report the introduction of specific SARS-CoV-2 variants at very high frequency in distinct geographic locations, which is important for understanding the evolution and circulation of SARS-CoV-2 variants of public health importance, while it circulates in humans.

Phylogenetic analysis of SARS-CoV-2 in Boston highlights the impact of superspreading events

Analysis of 772 complete severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes from early in the Boston-area epidemic revealed numerous introductions of the virus, a small number of which led to most cases. The data revealed two superspreading events. One, in a skilled nursing facility, led to rapid transmission and significant mortality in this vulnerable population but little broader spread, whereas other introductions into the facility had little effect. The second, at an international business conference, produced sustained community transmission and was exported, resulting in extensive regional, national, and international spread. The two events also differed substantially in the genetic variation they generated, suggesting varying transmission dynamics in superspreading events. Our results show how genomic epidemiology can help to understand the link between individual clusters and wider community spread.

Molecular epidemiology of Plasmodium falciparum by multiplexed amplicon deep sequencing in Senegal

BACKGROUND

Molecular epidemiology can provide important information regarding the genetic diversity and transmission of Plasmodium falciparum, which can assist in designing and monitoring elimination efforts. However, malaria molecular epidemiology including understanding the genetic diversity of the parasite and performing molecular surveillance of transmission has been poorly documented in Senegal. Next Generation Sequencing (NGS) offers a practical, fast and high-throughput approach to understand malaria population genetics. This study aims to unravel the population structure of P. falciparum and to estimate the allelic diversity, multiplicity of infection (MOI), and evolutionary patterns of the malaria parasite using the NGS platform.

METHODS

Multiplex amplicon deep sequencing of merozoite surface protein 1 (PfMSP1) and merozoite surface protein 2 (PfMSP2) in fifty-three P. falciparum isolates from two epidemiologically different areas in the South and North of Senegal, was carried out.

RESULTS

A total of 76 Pfmsp1 and 116 Pfmsp2 clones were identified and 135 different alleles were found, 56 and 79 belonged to the pfmsp1 and pfmsp2 genes, respectively. K1 and IC3D7 allelic families were most predominant in both sites. The local haplotype diversity (Hd) and nucleotide diversity (π) were higher in the South than in the North for both genes. For pfmsp1, a high positive Tajima's D (TD) value was observed in the South (D = 2.0453) while negative TD value was recorded in the North (D = - 1.46045) and F-Statistic (Fst) was 0.19505. For pfmsp2, non-directional selection was found with a highly positive TD test in both areas and Fst was 0.02111. The mean MOI for both genes was 3.07 and 1.76 for the South and the North, respectively, with a statistically significant difference between areas (p = 0.001).

CONCLUSION

This study revealed a high genetic diversity of pfmsp1 and pfmsp2 genes and low genetic differentiation in P. falciparum population in Senegal. The MOI means were significantly different between the Southern and Northern areas. Findings also showed that multiplexed amplicon deep sequencing is a useful technique to investigate genetic diversity and molecular epidemiology of P. falciparum infections.

Powassan Virus Neuropathology and Genomic Diversity in Patients With Fatal Encephalitis

Background

Powassan virus (POWV) is an emerging cause of severe encephalitis; very little is known about human pathogenicity due to challenges in diagnosis and viral RNA recovery. We present 3 patients with fatal encephalitis due to POWV lineage II (deer tick virus).

Methods

We obtained 27 unique samples, including from brain biopsy and autopsy, and used metagenomic sequencing, quantitative reverse transcriptase polymerase chain reaction, and a newly developed CRISPR-based diagnostic assay to perform the first detailed characterization of POWV compartmentalization and genomics between and within human subjects.

Results

In all 3 patients, imaging and histopathology findings were notable for profound cerebellar involvement. All patients were initially diagnosed with POWV by metagenomic sequencing, and 2 of the 3 had negative clinical testing by serology. We detected POWV RNA in 13 clinical samples; levels were highest in the cerebellum, and there was very little involvement of peripheral tissue. We assembled complete POWV genomes from 8 samples, providing unique information about the strains of POWV lineage II (deer tick virus) that infect humans.

Conclusions

We demonstrate the utility of molecular assays for detecting POWV infection, including in seronegative patients, and nominate viral genomic features that may relate to human infection and neuropathogenicity. The cerebellum was identified as a key target POWV in fatal infection, by radiological and histopathological findings as well as molecular testing.

Phylogenetic analysis of SARS-CoV-2 in the Boston area highlights the role of recurrent importation and superspreading events

SARS-CoV-2 has caused a severe, ongoing outbreak of COVID-19 in Massachusetts with 111,070 confirmed cases and 8,433 deaths as of August 1, 2020. To investigate the introduction, spread, and epidemiology of COVID-19 in the Boston area, we sequenced and analyzed 772 complete SARS-CoV-2 genomes from the region, including nearly all confirmed cases within the first week of the epidemic and hundreds of cases from major outbreaks at a conference, a nursing facility, and among homeless shelter guests and staff. The data reveal over 80 introductions into the Boston area, predominantly from elsewhere in the United States and Europe. We studied two superspreading events covered by the data, events that led to very different outcomes because of the timing and populations involved. One produced rapid spread in a vulnerable population but little onward transmission, while the other was a major contributor to sustained community transmission, including outbreaks in homeless populations, and was exported to several other domestic and international sites. The same two events differed significantly in the number of new mutations seen, raising the possibility that SARS-CoV-2 superspreading might encompass disparate transmission dynamics. Our results highlight the failure of measures to prevent importation into MA early in the outbreak, underscore the role of superspreading in amplifying an outbreak in a major urban area, and lay a foundation for contact tracing informed by genetic data.

Genomic surveillance revealed prevalence of unique SARS-CoV- 2 variants bearing mutation in the RdRp gene among Nevada patients.. [PMID: 32869037 PMCID: PMC7457609 DOI: 10.1101/2020.08.21.20178863]

Patients with signs of COVID-19 were tested with CDC approved diagnostic RT-PCR for SARS-CoV-2 using RNA extracted from nasopharyngeal/nasal swabs. In order to determine the variants of SARS-CoV-2 circulating in the state of Nevada, 200 patient specimens from COVID-19 patients were sequenced through our robust protocol for sequencing SARS-CoV-2 genomes. Our protocol enabled sequencing of SARS-CoV-2 genome directly from the specimens, with even very low viral loads, without the need of culture-based amplification. This allowed the identification of specific nucleotide variants including those coding for D614G and clades defining mutations. These sequences were further analyzed for determining SARS-CoV-2 variants circulating in the state of Nevada and their phylogenetic relationships with other variants present in the united states and the world during the same period of the outbreak. Our study reports the occurrence of a novel variant in the nsp12 (RNA dependent RNA Polymerase) protein at residue 323 (314aa of orf1b) to Phenylalanine (F) from Proline (P), present in the original isolate of SARS-CoV-2 (Wuhan-Hu-1). This 323F variant is found at a very high frequency (46% of the tested specimen) in Northern Nevada. Functional significance of this unique and highly prevalent variant of SARS-CoV-2 with RdRp mutation is currently under investigation but structural modeling showed this 323aa residue in the interface domain of RdRp, which is required for association with accessory proteins. In conclusion, we report the introduction of specific SARS-CoV-2 variants at a very high frequency within a distinct geographic location, which is important for clinical and public health perspectives in understanding the evolution of SARS-CoV-2 while in circulation.

Deployable CRISPR-Cas13a diagnostic tools to detect and report Ebola and Lassa virus cases in real-time

Recent outbreaks of viral hemorrhagic fevers (VHFs), including Ebola virus disease (EVD) and Lassa fever (LF), highlight the urgent need for sensitive, deployable tests to diagnose these devastating human diseases. Here we develop CRISPR-Cas13a-based (SHERLOCK) diagnostics targeting Ebola virus (EBOV) and Lassa virus (LASV), with both fluorescent and lateral flow readouts. We demonstrate on laboratory and clinical samples the sensitivity of these assays and the capacity of the SHERLOCK platform to handle virus-specific diagnostic challenges. We perform safety testing to demonstrate the efficacy of our HUDSON protocol in heat-inactivating VHF viruses before SHERLOCK testing, eliminating the need for an extraction. We develop a user-friendly protocol and mobile application (HandLens) to report results, facilitating SHERLOCK’s use in endemic regions. Finally, we successfully deploy our tests in Sierra Leone and Nigeria in response to recent outbreaks.

Outbreaks of viral hemorrhagic fevers highlight the need for sensitive, field-deployable diagnostics. Here the authors present a CRISPR-based SHERLOCK platform with field protocol and mobile app for Ebola and Lassa fever outbreaks.

RNAseq Analysis of Rhizomania-Infected Sugar Beet Provides the First Genome Sequence of Beet Necrotic Yellow Vein Virus from the USA and Identifies a Novel Alphanecrovirus and Putative Satellite Viruses

“Rhizomania” of sugar beet is a soilborne disease complex comprised of beet necrotic yellow vein virus (BNYVV) and its plasmodiophorid vector, Polymyxa betae. Although BNYVV is considered the causal agent of rhizomania, additional viruses frequently accompany BNYVV in diseased roots. In an effort to better understand the virus cohort present in sugar beet roots exhibiting rhizomania disease symptoms, five independent RNA samples prepared from diseased beet seedlings reared in a greenhouse or from field-grown adult sugar beet plants and enriched for virus particles were subjected to RNAseq. In all but a healthy control sample, the technique was successful at identifying BNYVV and provided sequence reads of sufficient quantity and overlap to assemble > 98% of the published genome of the virus. Utilizing the derived consensus sequence of BNYVV, infectious RNA was produced from cDNA clones of RNAs 1 and 2. The approach also enabled the detection of beet soilborne mosaic virus (BSBMV), beet soilborne virus (BSBV), beet black scorch virus (BBSV), and beet virus Q (BVQ), with near-complete genome assembly afforded to BSBMV and BBSV. In one field sample, a novel virus sequence of 3682 nt was assembled with significant sequence similarity and open reading frame (ORF) organization to members within the subgenus Alphanecrovirus (genus Necrovirus; family Tombusviridae). Construction of a DNA clone based on this sequence led to the production of the novel RNA genome in vitro that was capable of inducing local lesion formation on leaves of Chenopodium quinoa. Additionally, two previously unreported satellite viruses were revealed in the study; one possessing weak similarity to satellite maize white line mosaic virus and a second possessing moderate similarity to satellite tobacco necrosis virus C. Taken together, the approach provides an efficient pipeline to characterize variation in the BNYVV genome and to document the presence of other viruses potentially associated with disease severity or the ability to overcome resistance genes used for sugar beet rhizomania disease management.

Massively multiplexed nucleic acid detection with Cas13

The great majority of globally circulating pathogens go undetected, undermining patient care and hindering outbreak preparedness and response. To enable routine surveillance and comprehensive diagnostic applications, there is a need for detection technologies that can scale to test many samples1-3 while simultaneously testing for many pathogens4-6. Here, we develop Combinatorial Arrayed Reactions for Multiplexed Evaluation of Nucleic acids (CARMEN), a platform for scalable, multiplexed pathogen detection. In the CARMEN platform, nanolitre droplets containing CRISPR-based nucleic acid detection reagents7 self-organize in a microwell array8 to pair with droplets of amplified samples, testing each sample against each CRISPR RNA (crRNA) in replicate. The combination of CARMEN and Cas13 detection (CARMEN-Cas13) enables robust testing of more than 4,500 crRNA-target pairs on a single array. Using CARMEN-Cas13, we developed a multiplexed assay that simultaneously differentiates all 169 human-associated viruses with at least 10 published genome sequences and rapidly incorporated an additional crRNA to detect the causative agent of the 2020 COVID-19 pandemic. CARMEN-Cas13 further enables comprehensive subtyping of influenza A strains and multiplexed identification of dozens of HIV drug-resistance mutations. The intrinsic multiplexing and throughput capabilities of CARMEN make it practical to scale, as miniaturization decreases reagent cost per test by more than 300-fold. Scalable, highly multiplexed CRISPR-based nucleic acid detection shifts diagnostic and surveillance efforts from targeted testing of high-priority samples to comprehensive testing of large sample sets, greatly benefiting patients and public health9-11.

Lassa Virus Genetics

In a pattern repeated across a range of ecological niches, arenaviruses have evolved a compact four-gene genome to orchestrate a complex life cycle in a narrow range of susceptible hosts. A number of mammalian arenaviruses cross-infect humans, often causing a life-threatening viral hemorrhagic fever. Among this group of geographically bound zoonoses, Lassa virus has evolved a unique niche that leads to significant and sustained human morbidity and mortality. As a biosafety level 4 pathogen, direct study of the pathogenesis of Lassa virus is limited by the sparse availability, high operating costs, and technical restrictions of the high-level biocontainment laboratories required for safe experimentation. In this chapter, we introduce the relationship between genome structure and the life cycle of Lassa virus and outline reverse genetic approaches used to probe and describe functional elements of the Lassa virus genome. We then review the tools used to obtain viral genomic sequences used for phylogeny and molecular diagnostics, before shifting to a population perspective to assess the contributions of phylogenetic analysis in understanding the evolution and ecology of Lassa virus in West Africa. We finally consider the future outlook and clinical applications for genetic study of Lassa virus.

Novel genome sequences of cell-fusing agent virus allow comparison of virus phylogeny with the genetic structure of Aedes aegypti populations

Flaviviruses encompass not only medically relevant arthropod-borne viruses (arboviruses) but also insect-specific flaviviruses (ISFs) that are presumably maintained primarily through vertical transmission in the insect host. Interestingly, ISFs are commonly found infecting important arbovirus vectors such as the mosquito Aedes aegypti. Cell-fusing agent virus (CFAV) was the first described ISF of mosquitoes more than four decades ago. Despite evidence for widespread CFAV infections in A.aegypti populations and for CFAV potential to interfere with arbovirus transmission, little is known about CFAV evolutionary history. Here, we generated six novel CFAV genome sequences by sequencing three new virus isolates and subjecting three mosquito samples to untargeted viral metagenomics. We used these new genome sequences together with published ones to perform a global phylogenetic analysis of CFAV genetic diversity. Although there was some degree of geographical clustering among CFAV sequences, there were also notable discrepancies between geography and phylogeny. In particular, CFAV sequences from Cambodia and Thailand diverged significantly, despite confirmation that A.aegypti populations from both locations are genetically close. The apparent phylogenetic discrepancy between CFAV and its A.aegypti host in Southeast Asia indicates that other factors than host population structure shape CFAV genetic diversity.

Amplicon deep sequencing of kelch13 in Plasmodium falciparum isolates from Senegal

BACKGROUND

In 2006, the Senegalese National Malaria Control Programme recommended artemisinin-based combination therapy (ACT) with artemether-lumefantrine as the first-line treatment for uncomplicated Plasmodium falciparum malaria. To date, multiple mutations associated with artemisinin delayed parasite clearance have been described in Southeast Asia in the Pfk13 gene, such as Y493H, R539T, I543T and C580Y. Even though ACT remains clinically and parasitologically efficacious in Senegal, the spread of resistance is possible as shown by the earlier emergence of resistance to chloroquine in Southeast Asia that subsequently spread to Africa. Therefore, surveillance of artemisinin resistance in malaria endemic regions is crucial and requires the implementation of sensitive tools, such as next-generation sequencing (NGS) which can detect novel mutations at low frequency.

METHODS

Here, an amplicon sequencing approach was used to identify mutations in the Pfk13 gene in eighty-one P. falciparum isolates collected from three different regions of Senegal.

RESULTS

In total, 10 SNPs around the propeller domain were identified; one synonymous SNP and nine non-synonymous SNPs, and two insertions. Three of these SNPs (T478T, A578S and V637I) were located in the propeller domain. A578S, is the most frequent mutation observed in Africa, but has not previously been reported in Senegal. A previous study has suggested that A578S could disrupt the function of the Pfk13 propeller region.

CONCLUSION

As the genetic basis of possible artemisinin resistance may be distinct in Africa and Southeast Asia, further studies are necessary to assess the new SNPs reported in this study.

Spatially distinct physiology of Bacteroides fragilis within the proximal colon of gnotobiotic mice

A complex microbiota inhabits various microenvironments of the gut, with some symbiotic bacteria having evolved traits to invade the epithelial mucus layer and reside deep within intestinal tissue of animals. Whether these distinct bacterial communities across gut biogeographies exhibit divergent behaviors remains largely unknown. Global transcriptomic analysis to investigate microbial physiology in specific mucosal niches has been hampered technically by overabundance of host RNA. Herein, we employed hybrid selection RNA sequencing (hsRNA-Seq) to enable detailed spatial transcriptomic profiling of a prominent human commensal as it colonizes the colonic lumen, mucus or epithelial tissue of mice. Compared to conventional RNA-Seq, hsRNA-Seq increased reads mapping to the Bacteroides fragilis genome by 48- and 154-fold in mucus and tissue, respectively, allowing for high fidelity comparisons across biogeographic sites. Near the epithelium, B. fragilis up-regulated numerous genes involved in protein synthesis, indicating that bacteria inhabiting the mucosal niche are metabolically active. Further, a specific sulfatase (BF3086) and glycosyl hydrolase (BF3134) were highly induced in mucus and tissue compared to bacteria in the lumen. In-frame deletion of these genes impaired in vitro growth on mucus as a carbon source, as well as mucosal colonization of mice. Mutants in either B. fragilis gene displayed a fitness defect in competing for colonization against bacterial challenge, revealing the importance of site-specific gene expression for robust host-microbial symbiosis. As a versatile tool, hsRNA-Seq can be deployed to explore the in vivo spatial physiology of numerous bacterial pathogens or commensals.

Microbiome analyses of blood and tissues suggest cancer diagnostic approach

Systematic characterization of the cancer microbiome provides the opportunity to develop techniques that exploit non-human, microorganism-derived molecules in the diagnosis of a major human disease. Following recent demonstrations that some types of cancer show substantial microbial contributions1-10, we re-examined whole-genome and whole-transcriptome sequencing studies in The Cancer Genome Atlas11 (TCGA) of 33 types of cancer from treatment-naive patients (a total of 18,116 samples) for microbial reads, and found unique microbial signatures in tissue and blood within and between most major types of cancer. These TCGA blood signatures remained predictive when applied to patients with stage Ia-IIc cancer and cancers lacking any genomic alterations currently measured on two commercial-grade cell-free tumour DNA platforms, despite the use of very stringent decontamination analyses that discarded up to 92.3% of total sequence data. In addition, we could discriminate among samples from healthy, cancer-free individuals (n = 69) and those from patients with multiple types of cancer (prostate, lung, and melanoma; 100 samples in total) solely using plasma-derived, cell-free microbial nucleic acids. This potential microbiome-based oncology diagnostic tool warrants further exploration.

Current Trends in Diagnostics of Viral Infections of Unknown Etiology

Viruses are evolving at an alarming rate, spreading and inconspicuously adapting to cutting-edge therapies. Therefore, the search for rapid, informative and reliable diagnostic methods is becoming urgent as ever. Conventional clinical tests (PCR, serology, etc.) are being continually optimized, yet provide very limited data. Could high throughput sequencing (HTS) become the future gold standard in molecular diagnostics of viral infections? Compared to conventional clinical tests, HTS is universal and more precise at profiling pathogens. Nevertheless, it has not yet been widely accepted as a diagnostic tool, owing primarily to its high cost and the complexity of sample preparation and data analysis. Those obstacles must be tackled to integrate HTS into daily clinical practice. For this, three objectives are to be achieved: (1) designing and assessing universal protocols for library preparation, (2) assembling purpose-specific pipelines, and (3) building computational infrastructure to suit the needs and financial abilities of modern healthcare centers. Data harvested with HTS could not only augment diagnostics and help to choose the correct therapy, but also facilitate research in epidemiology, genetics and virology. This information, in turn, could significantly aid clinicians in battling viral infections.

Combining genomics and epidemiology to track mumps virus transmission in the United States

Unusually large outbreaks of mumps across the United States in 2016 and 2017 raised questions about the extent of mumps circulation and the relationship between these and prior outbreaks. We paired epidemiological data from public health investigations with analysis of mumps virus whole genome sequences from 201 infected individuals, focusing on Massachusetts university communities. Our analysis suggests continuous, undetected circulation of mumps locally and nationally, including multiple independent introductions into Massachusetts and into individual communities. Despite the presence of these multiple mumps virus lineages, the genomic data show that one lineage has dominated in the US since at least 2006. Widespread transmission was surprising given high vaccination rates, but we found no genetic evidence that variants arising during this outbreak contributed to vaccine escape. Viral genomic data allowed us to reconstruct mumps transmission links not evident from epidemiological data or standard single-gene surveillance efforts and also revealed connections between apparently unrelated mumps outbreaks.

Evaluation of a portable nanopore-based sequencer for detection of viruses in water

The newly emerged nanopore sequencing technology such as MinION™ allows for real-time detection of long DNA/RNA fragments on a portable device, yet few have examined its performance for environmental viromes. Here we seeded one RNA virus bacteriophage MS2 and one DNA virus bacteriophage PhiX174 into 10 L well water at three levels ranging from 1 to 21,100 plaque-forming units (PFU)/mL. Two workflows were established to maximize the number of sequencing reads of RNA and DNA viruses using MinION™. With dead-end ultrafiltration, PEG precipitation, and random amplification, MinION™ was capable of detecting MS2 at 155 PFU/mL and PhiX174 at 1-2 PFU/mL. While the DNA workflow only detected PhiX174, the RNA workflow detected both MS2 and PhiX174. The virus concentration, or relative abundance of viral nucleic acids in total nucleic acids, is critical to the proportion of viral reads in sequencing results. Our findings also highlight the importance of including control samples in sequencing runs for environmental water samples with low virus abundance.