Reduced subgenomic RNA expression is a molecular indicator of asymptomatic SARS-CoV-2 infection

Demultiplexing and barcode-specific adaptive sampling for nanopore direct RNA sequencing

Nanopore direct RNA sequencing (dRNA-seq) enables unique insights into RNA biology. However, applications are currently limited by the lack of accurate and cost-effective sample multiplexing. Here we introduce WarpDemuX, an ultra-fast and highly accurate adapter-barcoding and demultiplexing approach for dRNA-seq with SQK-RNA002 and SQK-RNA004 chemistries. WarpDemuX enhances speed and accuracy by fast processing of the raw nanopore signal, use of a light-weight machine-learning algorithm and design of optimized barcode sets. We demonstrate its utility by performing rapid phenotypic profiling of different SARS-CoV-2 viruses through multiplexed sequencing of longitudinal samples on a single flowcell, identifying systematic differences in transcript abundance and poly(A) tail lengths during infection. Additionally, integrating WarpDemuX into sequencing control software enables real-time enrichment of target molecules through barcode-specific adaptive sampling, which we demonstrate by enriching low abundance viral RNA. In summary, WarpDemuX represents a broadly applicable, high-performance, economical multiplexing solution for dRNA-seq, facilitating advanced (epi-) transcriptomic research.

Vaccination Shapes Within-Host SARS-CoV-2 Diversity of Omicron BA.2.2 Breakthrough Infection

BACKGROUND

Low-frequency intrahost single-nucleotide variants of SARS-CoV-2 have been recognized as predictive indicators of selection. However, the impact of vaccination on the intrahost evolution of SARS-CoV-2 remains uncertain at present.

METHODS

We investigated the genetic variation of SARS-CoV-2 in individuals who were unvaccinated, partially vaccinated, or fully vaccinated during Shanghai's Omicron BA.2.2 wave. We substantiated the connection between particular amino acid substitutions and immune-mediated selection through a pseudovirus neutralization assay or by cross-verification with the human leukocyte antigen-associated T-cell epitopes.

RESULTS

In contrast to those with immunologic naivety or partial vaccination, participants who were fully vaccinated had intrahost variant spectra characterized by reduced diversity. Nevertheless, the distribution of mutations in the fully vaccinated group was enriched in the spike protein. The distribution of intrahost single-nucleotide variants in individuals who were immunocompetent did not demonstrate notable signs of positive selection, in contrast to the observed adaptation in 2 participants who were immunocompromised who had an extended period of viral shedding.

CONCLUSIONS

In SARS-CoV-2 infections, vaccine-induced immunity was associated with decreased diversity of within-host variant spectra, with milder inflammatory pathophysiology. The enrichment of mutations in the spike protein gene indicates selection pressure exerted by vaccination on the evolution of SARS-CoV-2.

Vertical transmission of SARS-CoV-2 delta-variant in a preterm infant

BACKGROUND

As SARS-CoV-2 continues to be relevant and cause illnesses, the effect of emerging virus variants on perinatal health remains to be elucidated. It was demonstrated that vertical transmission of SARS-CoV-2 is a relatively rare event in the original SARS-CoV-2 strain. However, very few reports describe vertical transmission related to the delta-variant.

CASE PRESENTATION

We report a case of a preterm male neonate born to a mother with positive SARS-CoV-2 and mild respiratory complications. The neonate was born by cesarean section due to fetal distress. The rupture of the amniotic membrane was at delivery. The neonate had expected prematurity-related complications. His nasopharyngeal swabs for RT-PCR were positive from birth till three weeks of age. RT-ddPCR of the Placenta showed a high load of the SARS-CoV-2 virus with subgenomic viral RNA. RNAscope technique demonstrated both the positive strand of the S gene and the orf1ab negative strand. Detection of subgenomic RNA and the orf1ab negative strand indicats active viral replication in the placenta.

CONCLUSIONS

Our report demonstrates active viral replication of the SARS-CoV-2 delta-variant in the placenta associated with vertical transmission in a preterm infant.

Quantitative and qualitative subgenomic RNA profiles of SARS-CoV-2 in respiratory samples: A comparison between Omicron BA.2 and non-VOC-D614G

The SARS-CoV-2 Omicron variants are notorious for their transmissibility, but little is known about their subgenomic RNA (sgRNA) expression. This study applied RNA-seq to delineate the quantitative and qualitative profiles of canonical sgRNA of 118 respiratory samples collected from patients infected with Omicron BA.2 and compared with 338 patients infected with non-variant of concern (non-VOC)-D614G. A unique characteristic profile depicted by the relative abundance of 9 canonical sgRNAs was reproduced by both BA.2 and non-VOC-D614G regardless of host gender, age and presence of pneumonia. Remarkably, such profile was lost in samples with low viral load, suggesting a potential application of sgRNA pattern to indicate viral activity of individual patient at a specific time point. A characteristic qualitative profile of canonical sgRNAs was also reproduced by both BA.2 and non-VOC-D614G. The presence of a full set of canonical sgRNAs carried a coherent correlation with crude viral load (AUC ​= ​0.91, 95% CI 0.88-0.94), and sgRNA ORF7b was identified to be the best surrogate marker allowing feasible routine application in characterizing the infection status of individual patient. Further potentials in using sgRNA as a target for vaccine and antiviral development are worth pursuing.

Generation and Functional Analysis of Defective Viral Genomes during SARS-CoV-2 Infection

Defective viral genomes (DVGs) have been identified in many RNA viruses as a major factor influencing antiviral immune response and viral pathogenesis. However, the generation and function of DVGs in SARS-CoV-2 infection are less known. In this study, we elucidated DVG generation in SARS-CoV-2 and its relationship with host antiviral immune response. We observed DVGs ubiquitously from transcriptome sequencing (RNA-seq) data sets of in vitro infections and autopsy lung tissues of COVID-19 patients. Four genomic hot spots were identified for DVG recombination, and RNA secondary structures were suggested to mediate DVG formation. Functionally, bulk and single-cell RNA-seq analysis indicated the interferon (IFN) stimulation of SARS-CoV-2 DVGs. We further applied our criteria to the next-generation sequencing (NGS) data set from a published cohort study and observed a significantly higher amount and frequency of DVG in symptomatic patients than those in asymptomatic patients. Finally, we observed exceptionally diverse DVG populations in one immunosuppressive patient up to 140 days after the first positive test of COVID-19, suggesting for the first time an association between DVGs and persistent viral infections in SARS-CoV-2. Together, our findings strongly suggest a critical role of DVGs in modulating host IFN responses and symptom development, calling for further inquiry into the mechanisms of DVG generation and into how DVGs modulate host responses and infection outcome during SARS-CoV-2 infection. IMPORTANCE Defective viral genomes (DVGs) are generated ubiquitously in many RNA viruses, including SARS-CoV-2. Their interference activity to full-length viruses and IFN stimulation provide the potential for them to be used in novel antiviral therapies and vaccine development. SARS-CoV-2 DVGs are generated through the recombination of two discontinuous genomic fragments by viral polymerase complex, and this recombination is also one of the major mechanisms for the emergence of new coronaviruses. Focusing on the generation and function of SARS-CoV-2 DVGs, these studies identify new hot spots for nonhomologous recombination and strongly suggest that the secondary structures within viral genomes mediate the recombination. Furthermore, these studies provide the first evidence for IFN stimulation activity of de novo DVGs during natural SARS-CoV-2 infection. These findings set up the foundation for further mechanism studies of SARS-CoV-2 recombination and provide evidence to harness the immunostimulatory potential of DVGs in the development of a vaccine and antivirals for SARS-CoV-2.

sgRNAs: A SARS-CoV-2 emerging issue

Like for other coronaviruses, SARS-CoV-2 gene expression strategy is based on the synthesis of a nested set of subgenomic mRNA species (sgRNAs). These sgRNA are synthesized using a "discontinuous transcription" mechanism that relies on template switching at Transcription Regulatory Sequences (TRS). Both canonical (c-sgRNA) and non-canonical (nc-sgRNA, less numerous) subgenomic RNA species can be produced. Currently, sgRNAs are investigated on the basis of sequence data obtained through next generation sequencing (NGS), and bioinformatic tools are crucial for their identification, characterization and quantification. To date, few software have been developed to this aim, whose reliability and applicability to all the available NGS platforms need to be established, to build confidence on the information resulting from such tools. In fact, these information may be crucial for the in depth elucidation of viral expression strategy, particularly in respect of the significance of nc-sgRNAs, and for the possible use of sgRNAs as potential markers of virus replicative activity in infected patients.

Evaluation of SARS-CoV-2 ORF7a Deletions from COVID-19-Positive Individuals and Its Impact on Virus Spread in Cell Culture

The spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), causing the COVID-19 outbreak, posed a primary concern of public health worldwide. The most common changes in SARS-CoV-2 are single nucleotide substitutions, also reported insertions and deletions. This work investigates the presence of SARS-CoV-2 ORF7a deletions identified in COVID-19-positive individuals. Sequencing of SARS-CoV-2 complete genomes showed three different ORF7a size deletions (190-nt, 339-nt and 365-nt). Deletions were confirmed through Sanger sequencing. The ORF7a∆190 was detected in a group of five relatives with mild symptoms of COVID-19, and the ORF7a∆339 and ORF7a∆365 in a couple of co-workers. These deletions did not affect subgenomic RNAs (sgRNA) production downstream of ORF7a. Still, fragments associated with sgRNA of genes upstream of ORF7a showed a decrease in size when corresponding to samples with deletions. In silico analysis suggests that the deletions impair protein proper function; however, isolated viruses with partial deletion of ORF7a can replicate in culture cells similarly to wild-type viruses at 24 hpi, but with less infectious particles after 48 hpi. These findings on deleted ORF7a accessory protein gene, contribute to understanding SARS-CoV-2 phenotypes such as replication, immune evasion and evolutionary fitness as well insights into the role of SARS-CoV-2_ORF7a in the mechanism of virus-host interactions.

Generation and functional analysis of defective viral genomes during SARS-CoV-2 infection

Defective viral genomes (DVGs) have been identified in many RNA viruses as a major factor influencing antiviral immune response and viral pathogenesis. However, the generation and function of DVGs in SARS-CoV-2 infection are less known. In this study, we elucidated DVG generation in SARS-CoV-2 and its relationship with host antiviral immune response. We observed DVGs ubiquitously from RNA-seq datasets of in vitro infections and autopsy lung tissues of COVID-19 patients. Four genomic hotspots were identified for DVG recombination and RNA secondary structures were suggested to mediate DVG formation. Functionally, bulk and single cell RNA-seq analysis indicated the IFN stimulation of SARS-CoV-2 DVGs. We further applied our criteria to the NGS dataset from a published cohort study and observed significantly higher DVG amount and frequency in symptomatic patients than that in asymptomatic patients. Finally, we observed unusually high DVG frequency in one immunosuppressive patient up to 140 days after admitted to hospital due to COVID-19, first-time suggesting an association between DVGs and persistent viral infections in SARS-CoV-2. Together, our findings strongly suggest a critical role of DVGs in modulating host IFN responses and symptom development, calling for further inquiry into the mechanisms of DVG generation and how DVGs modulate host responses and infection outcome during SARS-CoV-2 infection.

Importance

Defective viral genomes (DVGs) are ubiquitously generated in many RNA viruses, including SARS-CoV-2. Their interference activity to full-length viruses and IFN stimulation provide them the potential for novel antiviral therapies and vaccine development. SARS-CoV-2 DVGs are generated through the recombination of two discontinuous genomic fragments by viral polymerase complex and the recombination is also one of the major mechanisms for the emergence of new coronaviruses. Focusing on the generation and function of SARS-CoV-2 DVGs, these studies identify new hotspots for non-homologous recombination and strongly suggest that the secondary structures within viral genomes mediate the recombination. Furthermore, these studies provide the first evidence for IFN stimulation activity of de novo DVGs during natural SARS-CoV-2 infection. These findings set up the foundation for further mechanism studies of SARS-CoV-2 recombination and provide the evidence to harness DVGs’ immunostimulatory potential in the development of vaccine and antivirals for SARS-CoV-2.

Detection of subgenomic mRNA from endemic human coronavirus OC43 and NL63 compared to viral genomic loads, single virus detection and clinical manifestations in children with respiratory tract infections

BACKGROUND

The importance of endemic human coronavirus (HCoV) in children has been insufficiently elucidated upon. Our aims were to develop subgenomic (sg) mRNA tests for HCoV species OC43 and NL63, and to evaluate the relationships to HCoV genomic loads, single HCoV detections and clinical manifestations.

METHODS

We have used an 11-yearlong cohort study of children admitted with respiratory tract infection (RTI) and hospital controls. Nasopharyngeal aspirates were analyzed for HCoV subtypes OC43 and NL63 with in-house diagnostic PCR. Positive samples were tested with newly developed real-time PCRs targeting sg mRNA coding for the nucleocapsid protein.

RESULTS

OC43 sg mRNA was detected in 86% (105/122) of available OC43-positive samples in the RTI group, and in 63% (12/19) of control samples. NL63 sg mRNA was detected in 72% (71/98) and 71% (12/17) of available NL63-positive patient and control samples, respectively. In RTI samples, sg mRNA detection was strongly associated with a Ct value <32 in both diagnostic PCR tests (OC43: OR = 54, 95% CI [6.8-428]; NL63: OR = 42, 95% CI [9.0-198]) and single NL63 detections (OR = 6.9, 95% CI [1.5-32]). Comparing RTI and controls, only OC43 was associated with RTI when adjusted for age (aOR = 3.2, 95% CI [1.1-9.4]).

CONCLUSION

We found strong associations between OC43 and NL63 sg mRNA and high viral genomic loads. sg mRNA for OC43 was associated with RTI. The association between sg mRNA and clinical manifestations needs further evaluation.

SARS-CoV-2 Within-Host and in vitro Genomic Variability and Sub-Genomic RNA Levels Indicate Differences in Viral Expression Between Clinical Cohorts and in vitro Culture

Background

Low frequency intrahost single nucleotide variants (iSNVs) of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) have been increasingly recognised as predictive indicators of positive selection. Particularly as growing numbers of SARS-CoV-2 variants of interest (VOI) and concern (VOC) emerge. However, the dynamics of subgenomic RNA (sgRNA) expression and its impact on genomic diversity and infection outcome remain poorly understood. This study aims to investigate and quantify iSNVs and sgRNA expression in single and longitudinally sampled cohorts over the course of mild and severe SARS-CoV-2 infection, benchmarked against an in vitro infection model.

Methods

Two clinical cohorts of SARS-CoV-2 positive cases in New South Wales, Australia collected between March 2020 and August 2021 were sequenced. Longitudinal samples from cases hospitalised due to SARS-CoV-2 infection (severe) (n = 16) were analysed and compared with cases that presented with SARS-CoV-2 symptoms but were not hospitalised (mild) (n = 23). SARS-CoV-2 genomic diversity profiles were also examined from daily sampling of culture experiments for three SARS-CoV-2 variants (Lineage A, B.1.351, and B.1.617.2) cultured in VeroE6 C1008 cells (n = 33).

Results

Intrahost single nucleotide variants were detected in 83% (19/23) of the mild cohort cases and 100% (16/16) of the severe cohort cases. SNP profiles remained relatively fixed over time, with an average of 1.66 SNPs gained or lost, and an average of 4.2 and 5.9 low frequency variants per patient were detected in severe and mild infection, respectively. sgRNA was detected in 100% (25/25) of the mild genomes and 92% (24/26) of the severe genomes. Total sgRNA expressed across all genes in the mild cohort was significantly higher than that of the severe cohort. Significantly higher expression levels were detected in the spike and the nucleocapsid genes. There was significantly less sgRNA detected in the culture dilutions than the clinical cohorts.

Discussion and Conclusion

The positions and frequencies of iSNVs in the severe and mild infection cohorts were dynamic overtime, highlighting the importance of continual monitoring, particularly during community outbreaks where multiple SARS-CoV-2 variants may co-circulate. sgRNA levels can vary across patients and the overall level of sgRNA reads compared to genomic RNA can be less than 1%. The relative contribution of sgRNA to the severity of illness warrants further investigation given the level of variation between genomes. Further monitoring of sgRNAs will improve the understanding of SARS-CoV-2 evolution and the effectiveness of therapeutic and public health containment measures during the pandemic.

Virus-like Particles: Measures and Biological Functions

Virus-like particles resemble infectious virus particles in size, shape, and molecular composition; however, they fail to productively infect host cells. Historically, the presence of virus-like particles has been inferred from total particle counts by microscopy, and infectious particle counts or plaque-forming-units (PFUs) by plaque assay; the resulting ratio of particles-to-PFUs is often greater than one, easily 10 or 100, indicating that most particles are non-infectious. Despite their inability to hijack cells for their reproduction, virus-like particles and the defective genomes they carry can exhibit a broad range of behaviors: interference with normal virus growth during co-infections, cell killing, and activation or inhibition of innate immune signaling. In addition, some virus-like particles become productive as their multiplicities of infection increase, a sign of cooperation between particles. Here, we review established and emerging methods to count virus-like particles and characterize their biological functions. We take a critical look at evidence for defective interfering virus genomes in natural and clinical isolates, and we review their potential as antiviral therapeutics. In short, we highlight an urgent need to better understand how virus-like genomes and particles interact with intact functional viruses during co-infection of their hosts, and their impacts on the transmission, severity, and persistence of virus-associated diseases.

Predicting infectivity: comparing four PCR-based assays to detect culturable SARS-CoV-2 in clinical samples

With the COVID-19 pandemic caused by SARS-CoV-2 now in its second year, there remains an urgent need for diagnostic testing that can identify infected individuals, particularly those who harbor infectious virus. Various RT-PCR strategies have been proposed to identify specific viral RNA species that may predict the presence of infectious virus, including detection of transcriptional intermediates (e.g., subgenomic RNA [sgRNA]) and replicative intermediates (e.g., negative-strand RNA species). Using a novel primer/probe set for detection of subgenomic (sg)E transcripts, we successfully identified 100% of specimens containing culturable SARS-CoV-2 from a set of 126 clinical samples (total sgE CT values ranging from 12.3 to 37.5). This assay showed superior performance compared to a previously published sgRNA assay and to a negative-strand RNA assay, both of which failed to detect target RNA in a subset of samples from which we isolated live virus. In addition, total levels of viral RNA (genome, negative-strand, and sgE) detected with the WHO/Charité primer-probe set correlated closely with levels of infectious virus. Specifically, infectious virus was not detected in samples with a CT above 31.0. Clinical samples with higher levels of viral RNA also displayed cytopathic effect (CPE) more quickly than those with lower levels of viral RNA. Finally, we found that the infectivity of SARS-CoV-2 samples is significantly dependent on the cell type used for viral isolation, as Vero E6 cells expressing TMRPSS2 extended the analytical sensitivity of isolation by more than 3 CT compared to parental Vero E6 cells and resulted in faster isolation. Our work shows that using a total viral RNA Ct cutoff of > 31 or specifically testing for sgRNA can serve as an effective rule-out test for the presence of culturable virus.

Nanopore-Based Direct RNA-Sequencing Reveals a High-Resolution Transcriptional Landscape of Porcine Reproductive and Respiratory Syndrome Virus

The TRS-mediated discontinuous transcription process is a hallmark of Arteriviruses. Precise assessment of the intricate subgenomic RNA (sg mRNA) populations is required to understand the kinetics of viral transcription. It is difficult to reconstruct and comprehensively quantify splicing events using short-read sequencing, making the identification of transcription-regulatory sequences (TRS) particularly problematic. Here, we applied long-read direct RNA sequencing to characterize the recombined RNA molecules produced in porcine alveolar macrophages during early passage infection of porcine reproductive and respiratory syndrome virus (PRRSV). Based on sequencing two PRRSV isolates, namely XM-2020 and GD, we revealed a high-resolution and diverse transcriptional landscape in PRRSV. The data revealed intriguing differences in subgenomic recombination types between the two PRRSVs while also demonstrating TRS-independent heterogeneous subpopulation not previously observed in Arteriviruses. We find that TRS usage is a regulated process and share the common preferred TRS in both strains. This study also identified a substantial number of TRS-mediated transcript variants, including alternative-sg mRNAs encoding the same annotated ORF, as well as putative sg mRNAs encoded nested internal ORFs, implying that the genetic information encoded in PRRSV may be more intensively expressed. Epigenetic modifications have emerged as an essential regulatory layer in gene expression. Here, we gained a deeper understanding of m5C modification in poly(A) RNA, elucidating a potential link between methylation and transcriptional regulation. Collectively, our findings provided meaningful insights for redefining the transcriptome complexity of PRRSV. This will assist in filling the research gaps and developing strategies for better control of the PRRS.