Validation of our viral metagenomics sequencing and analysis pipeline for SARS-CoV-2 detection

In this article we’re sharing results from the validation study scientists at Diversigen and DNA Genotek undertook to ensure Diversigen’s viral metagenomics capabilities could be leveraged to help during the COVID-19 pandemic.

When I was a boy and I would see scary things in the news, my mother would say to me, “Look for the helpers. You will always find people who are helping.

 ― Fred Rogers (a.k.a. Mr Rogers)

Like we hear in the news so often, these are extraordinary (and scary) times. We’ve all had to adapt our work and home lives to stay safe and healthy during the COVID-19 global pandemic. Here at Diversigen we’re working hard to ensure we can continue to offer the level of service and support you expect from us, while keeping our coworkers and loved ones safe. We’ve also watched as scientists, health care professionals, and essential service workers around the world have switched their focus from their own research, work or life goals and are asking “how can I help?”  We’re doing the same.

Validation of Diversigen’s viral metagenomics sequencing and analysis pipeline for SARS-CoV-2 detection

Lisa Gamwell1 PhD,     Aaron Garoutte2 PhD,     Jean Macklaim3 PhD,     and Emily Hollister4PhD

1. Product Manager, DNA Genotek 2. Bioinformatician, Diversigen 3. Bioinformatician, DNA Genotek 4. Director, R&D,Computational Biology, Diversigen

Diversigen has years of experience processing, sequencing and analysing samples using our proprietary viral metagenomic pipeline. Depending on the genomes of the viruses of interest, collected samples, extracted RNA or DNA, cDNA or DNA libraries and sequence data (shallow or deep, RNA or DNA based) can all be used as inputs for this pipeline. In addition to nasal nasopharyngeal samples, there are reports that SARS-CoV-2, the virus that causes COVID-19, can be detected in samples that are often of interest to microbiome researchers, including saliva and feces

In response to the COVID-19 global pandemic, we updated the Diversigen viral database and validated our viral metagenomics pipeline to ensure they can return information about SARS-CoV-2. We determined that the pipeline can detect SARS-CoV-2 in confirmed positive samples. Using a mock community that includes 9 coronavirus strains, we determined that it can differentiate between coronavirus species. In addition, it can recover the SARS-CoV-2 genome at a depth and coverage that allows for genome assembly, which is crucial for strain and/or mutation rate tracking.

Methods and Results

To validate the viral metagenomics pipeline for the detection of SARS-CoV-2, we analyzed 5 publicly available, SARS-CoV-2 positive and 2 SARS-CoV-2 negative bronchoalveolar lavage fluid samples from The Wuhan Institute of Virology and Institute of Pathogen Biology, Chinese Academy of Medical Sciences and Peking Union Medical College.  These samples were sequenced using an RNA sequencing approach. Using our viral classification pipeline, we detected SARS-CoV-2 in all positive samples but did not detect SARS-CoV-2 in the negative samples (Table 1).

Table 1: Summary of viral classification analysis of SARS-CoV-2 positive and negative samples

Sample ID SRA run identifier(s) Expected SARS-CoV-2 Status SARS-CoV-2 detected DNA viruses detected Other RNA viruses detected
WIV02 SRR11092058, SRR11092063 ➕ ✅ ✅ ✅
WIV04 SRR11092057, SRR11092062 ➕ ✅ ✅ ✅
WIV05 SRR11092061 ➕ ✅ ✅ ✅
WIV06 SRR11092056, SRR11092060 ➕ ✅ ✅ ✅
WIV07 SRR11092064, SRR11092059 ➕ ✅ ✅ ✅
COPD18 SRR5677628 ➖ ➖ ✅ ✅
COPD25 SRR5677642 ➖ ➖ ✅ ➖

In addition to demonstrating that the SARS-CoV-2 virus could be detected in RNA sequence libraries from confirmed COVID-19 cases, we also compared recovery of the virus and its genome from sequence libraries of varying sequence depth. At depths ranging from 1-2 Gigabases (MiSeq) at the shallow end of the spectrum to 10-20 Gigabases (MGISEQ-2000RS) at the deep end, they found that they could recover 70 to 100% of the SARS-COV-2 genome with average depth of coverage ranging from 3 to >100x . The most complete recovery came from the deep sequencing samples (e.g., Table 1, SRR11092062), where we recovered a complete SARS-CoV-2 genome with average coverage of >100x. The viral metagenomics pipeline can recover and return reads sufficient for assembling high quality viral genomes. This level of coverage can facilitate strain tracking and/or mutation rate tracking, using tools like NextStrain.

 

Figure 1: Genome coverage of SARS-CoV-2 using Diversigen’s viral classification pipeline to analyse cDNA libraries sequenced at variable depths.

Plot of genome coverage of SARS-CoV-2 using Diversigen’s viral classification pipeline to analyse cDNA libraries sequenced at variable depths.

Although the occurrence of multiple coronaviruses in a single human biological specimen is likely to be rare, we created a coronavirus mock community (in silico) to evaluate the ability of our pipeline to differentiate among a variety of coronaviruses. This mock community consisted of 9 coronavirus strains, including one Alphacoronavirus, four Betacoronaviruses (including SARS-CoV-2), two Gammacoronaviruses, one Deltacoronavirus, and one unclassified Coronavirus (Table 2).

 

Table 2: Coronavirus Mock Community sources and descriptions

Isolate (NCBI accession)

Genome size (bp)

Alphacoronavirus; Feline coronavirus isolate XXN (MN165107)

29286

Betacoronavirus; Middle East respiratory syndrome-related coronavirus isolate llama-passaged-Qatar15 (MN507638)

30033

Betacoronavirus; Dromedary camel coronavirus HKU23 isolate DcCoV-HKU23/camel/Ethiopia/CAC1019/2015 (MN514962)

31021

Betacoronavirus; Middle East respiratory syndrome-related coronavirus isolate Hu/Riyadh-KSA-18013832/2018 (MN723544)

30123

Betacoronavirus; Severe acute respiratory syndrome coronavirus 2 (MT121215)

29945

Gammacoronavirus; Infectious bronchitis virus isolate GA/1359/199 (MN566147)

27631

Gammacoronavirus; Infectious bronchitis virus isolate GA/1476/2015 (MN599049)

27666

Deltacoronavirus; Porcine deltacoronavirus (MN025260)

25399

Unclassified Coronavirinae; Hypsugo bat coronavirus HKU25 isolate YD131305 (KX442564)

30498

 

Illumina reads (2×150 bp) providing 10x coverage of each genome were generated, combined, and profiled (Figure 2).

 

 

 

Figure 2: Expected and observed relative abundances of the Coronavirus Mock Community

 

 

 

Visualization of expected and observed relative abundances of the Coronavirus Mock Community

Each virus in the mock community was recovered using our analysis pipeline. It distinguished among species within the same genus (e.g., Betacoronavirus), and the majority of observed reads counts were similar to their respective expected read counts. Reads that map ambiguously (i.e., to multiple organisms in the database) are filtered during analysis, which can lead to lower than expected recovery in some cases.

 

Summary

Diversigen‘s viral metagenomic pipeline can not only detect SARS-CoV-2 in confirmed positive samples, it can also differentiate between coronavirus species. In addition it can be used in combination with deep sequencing to recover the SARS-CoV-2 genome at a depth and coverage that allows for genome assembly, which is crucial for tracking genetic evolution of the virus. At Diversigen, we hope we can help researchers understand how SARS-CoV-2 may be impacting their clinical trials or human research studies and support screening of donor-derived microbiome products. In addition we hope to help the research community contribute to the global efforts to track this virus. 

For more information on our viral metagenomic pipeline and how we can help during the COVID-19 pandemic, please contact info@diversigen.com.

Additional Resources

Data Sharing Considerations

As always, we keep all customer data confidential, but we strongly encourage all researchers with SARS-CoV2 genetic sequence data to follow the World Health Organization (WHO) code of conduct for open and timely sharing of pathogen genetic data.  

For some suggestions on where, how and why you should share SARS-CoV2 genetic sequences, please check out these links:

WHERE

The Global Initiative on Sharing all Influenza Data (GISAID),and the National Institute of Health’s genetic sequence databases, GenBank and Sequence Read Archive (SRA), host viral genomes from ongoing outbreaks, including SARS-Cov-2. 

NextStrain, from the research groups of Trevor Bedford (Fred Hutchinson Cancer Research Center) and Richard Nehar (University of Basel, Max Planck Institute for Developmental Biology), provides a real-time view into the evolution and spread of a range of viral pathogens of high public health importance. Their COVID-19 specific app, Next hCoV-19 App, uses SAR-Cov2 genetic sequence data from GISAID to visualize the outbreak.

HOW 

Instructions on how to upload data to the GISAID database can be found here: Submitting Data to EpiFlu™

Instructions on how to upload data to GenBank can be found here: GenBank: BankIt 

WHY

Interested in learning more about how genetic epidemiology is useful during an outbreak? 

Genomic Study Points to Natural Origin of COVID-19

Software and Genetic Sequencing Track the Coronavirus’s Path

 

NEED FUNDING? MORE INFO?

The National Institutes of Health has funding opportunities specific to COVID-19 and maintains updated information and guidance for people involved in NIH research 

Coronavirus Disease 2019 (COVID-19): Information for NIH Applicants and Recipients of NIH Funding

The FDA also maintains updated information and guidance documents related to Coronavirus Disease 2019 (COVID-19), including Guidance on Conducting Clinical Trials during COVID-19 Pandemic 

 

Do you have an important resource to share with our community?  Email us and we’ll post it: info@diversigen.com

The information presented in this blog is made available for general information purposes and may include content provided by third parties. We do not warrant the accuracy, completeness or usefulness of this information and disclaim all liability arising from reliance on it. The statements and opinions expressed in these materials, other than the content provided by us, are solely those of the person or entity making the statements.