Assessing Bias and Reproducibility of Viral Metagenomics Methods for the Combined Detection of Faecal RNA and DNA Viruses

Whole transcriptome amplification (WTA2) and sequence-independent single primer amplification (SISPA) are two widely used methods for combined metagenomic sequencing of RNA and DNA viruses. However, information on the reproducibility and bias of these methods on diverse viruses in faecal samples is currently lacking. A mock community (MC) of diverse viruses was developed and used to spike faecal samples at different concentrations. Virus-like particles (VLPs) were extracted, nucleic acid isolated, reverse-transcribed, and PCR amplified using either WTA2 or SISPA and sequenced for metagenomic analysis. A bioinformatics pipeline measured the recovery of MC viruses in replicates of faecal samples from three human donors, analysing the consistency of viral abundance measures and taxonomy. Viruses had different recovery levels with VLP extraction introducing variability between replicates, while WTA2 and SISPA produced comparable results. In comparing WTA2- and SISPA-generated libraries, WTA2 gave more uniform coverage depth profiles and improved assembly quality and virus identification. SISPA produced more consistent abundance, with a 50% difference between replicates occurring in similar to 20% and similar to 10% of sequences for WTA2 and SISPA, respectively. In conclusion, a bioinformatics pipeline has been developed to assess the methodological variability and bias of WTA2 and SISPA, demonstrating higher sensitivity with WTA2 and higher consistency with SISPA.