sgDI-tector: defective interfering viral genome bioinformatics for detection of coronavirus subgenomic RNAs

SARS-CoV-2 subgenomic RNA: formation process and rapid molecular diagnostic methods

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) which caused coronavirus disease-2019 (COVID-19) is spreading worldwide and posing enormous losses to human health and socio-economic. Due to the limitations of medical and health conditions, it is still a huge challenge to develop appropriate discharge standards for patients with COVID-19 and to use medical resources in a timely and effective manner. Similar to other coronaviruses, SARS-CoV-2 has a very complex discontinuous transcription process to generate subgenomic RNA (sgRNA). Some studies support that sgRNA of SARS-CoV-2 can only exist when the virus is active and is an indicator of virus replication. The results of sgRNA detection in patients can be used to evaluate the condition of hospitalized patients, which is expected to save medical resources, especially personal protective equipment. There have been numerous investigations using different methods, especially molecular methods to detect sgRNA. Here, we introduce the process of SARS-CoV-2 sgRNA formation and the commonly used molecular diagnostic methods to bring a new idea for clinical detection in the future.

Accumulation Dynamics of Defective Genomes during Experimental Evolution of Two Betacoronaviruses

Virus-encoded replicases often generate aberrant RNA genomes, known as defective viral genomes (DVGs). When co-infected with a helper virus providing necessary proteins, DVGs can multiply and spread. While DVGs depend on the helper virus for propagation, they can in some cases disrupt infectious virus replication, impact immune responses, and affect viral persistence or evolution. Understanding the dynamics of DVGs alongside standard viral genomes during infection remains unclear. To address this, we conducted a long-term experimental evolution of two betacoronaviruses, the human coronavirus OC43 (HCoV-OC43) and the murine hepatitis virus (MHV), in cell culture at both high and low multiplicities of infection (MOI). We then performed RNA-seq at regular time intervals, reconstructed DVGs, and analyzed their accumulation dynamics. Our findings indicate that DVGs evolved to exhibit greater diversity and abundance, with deletions and insertions being the most common types. Notably, some high MOI deletions showed very limited temporary existence, while others became prevalent over time. We observed differences in DVG abundance between high and low MOI conditions in HCoV-OC43 samples. The size distribution of HCoV-OC43 genomes with deletions differed between high and low MOI passages. In low MOI lineages, short and long DVGs were the most common, with an additional cluster in high MOI lineages which became more prevalent along evolutionary time. MHV also showed variations in DVG size distribution at different MOI conditions, though they were less pronounced compared to HCoV-OC43, suggesting a more random distribution of DVG sizes. We identified hotspot regions for deletions that evolved at a high MOI, primarily within cistrons encoding structural and accessory proteins. In conclusion, our study illustrates the widespread formation of DVGs during betacoronavirus evolution, influenced by MOI and cell- and virus-specific factors.

Comparative analysis of bioinformatics tools to characterize SARS-CoV-2 subgenomic RNAs

During the replication of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), positive-sense genomic RNA and subgenomic RNAs (sgRNAs) are synthesized by a discontinuous process of transcription characterized by a template switch, regulated by transcription-regulating sequences (TRS). Although poorly known about makeup and dynamics of sgRNAs population and function of its constituents, next-generation sequencing approaches with the help of bioinformatics tools have made a significant contribution to expand the knowledge of sgRNAs in SARS-CoV-2. For this scope to date, Periscope, LeTRS, sgDI-tector, and CORONATATOR have been developed. However, limited number of studies are available to compare the performance of such tools. To this purpose, we compared Periscope, LeTRS, and sgDI-tector in the identification of canonical (c-) and noncanonical (nc-) sgRNA species in the data obtained with the Illumina ARTIC sequencing protocol applied to SARS-CoV-2-infected Caco-2 cells, sampled at different time points. The three software showed a high concordance rate in the identification and in the quantification of c-sgRNA, whereas more differences were observed in nc-sgRNA. Overall, LeTRS and sgDI-tector result to be adequate alternatives to Periscope to analyze Fastq data from sequencing platforms other than Nanopore.

DVGfinder: A Metasearch Tool for Identifying Defective Viral Genomes in RNA-Seq Data

The generation of different types of defective viral genomes (DVG) is an unavoidable consequence of the error-prone replication of RNA viruses. In recent years, a particular class of DVGs, those containing long deletions or genome rearrangements, has gain interest due to their potential therapeutic and biotechnological applications. Identifying such DVGs in high-throughput sequencing (HTS) data has become an interesting computational problem. Several algorithms have been proposed to accomplish this goal, though all incur false positives, a problem of practical interest if such DVGs have to be synthetized and tested in the laboratory. We present a metasearch tool, DVGfinder, that wraps the two most commonly used DVG search algorithms in a single workflow for the identification of the DVGs in HTS data. DVGfinder processes the results of ViReMa-a and DI-tector and uses a gradient boosting classifier machine learning algorithm to reduce the number of false-positive events. The program also generates output files in user-friendly HTML format, which can help users to explore the DVGs identified in the sample. We evaluated the performance of DVGfinder compared to the two search algorithms used separately and found that it slightly improves sensitivities for low-coverage synthetic HTS data and DI-tector precision for high-coverage samples. The metasearch program also showed higher sensitivity on a real sample for which a set of copy-backs were previously validated.