Transcriptional and epi-transcriptional dynamics of SARS-CoV-2 during cellular infection

Transcriptome-wide identification of 5-methylcytosine by deaminase and reader protein-assisted sequencing

5-Methylcytosine (m5C) is one of the posttranscriptional modifications in mRNA and is involved in the pathogenesis of various diseases. However, the capacity of existing assays for accurately and comprehensively transcriptome-wide m5C mapping still needs improvement. Here, we develop a detection method named DRAM (deaminase and reader protein assisted RNA methylation analysis), in which deaminases (APOBEC1 and TadA-8e) are fused with m5C reader proteins (ALYREF and YBX1) to identify the m5C sites through deamination events neighboring the methylation sites. This antibody-free and bisulfite-free approach provides transcriptome-wide editing regions which are highly overlapped with the publicly available bisulfite-sequencing (BS-seq) datasets and allows for a more stable and comprehensive identification of the m5C loci. In addition, DRAM system even supports ultralow input RNA (10 ng). We anticipate that the DRAM system could pave the way for uncovering further biological functions of m5C modifications.

Characterizing Host microRNA: Virus Interactions of Orthoavulavirus javaense

Post-transcriptional gene regulation mediated by microRNAs (miRNAs) relies on sequence complementarity between the miRNA seed site and the target gene transcript(s). This complementarity can completely inhibit or reduce translation into protein. We hypothesized that viruses employ sequence complementarity/similarity with host miRNAs to inhibit or increase the miRNA-mediated regulation of host gene expression specifically during viral infection(s). In this study, we focus on Orthoavulavirus javaense (OAVJ), the causative of Newcastle disease, a poultry disease with significant economic impact. A computational analysis of OAVJ genomes from low-virulence (lentogenic) versus virulent (velogenic) viruses was carried out to identify viral signature motifs that potentially either mimic or complement host miRNA seed sequences. Data show that OAVJ genomes harbor viral seed mimics (vSMs) or viral seed sponges (vSSs) and can mimic host miRNAs or inhibit their regulation of host genes, disrupting cellular pathways. Our analyses showed that velogens encode a statistically significant higher number of vSMs and a lower number of vSSs relative to lentogens. The number of vSMs or vSSs did not correlate with gene length. The analysis of the secondary structures flanking these vSMs and vSSs showed structural features common to miRNA precursors. The inhibition or upregulation of vSS-miR-27b-5p altered P gene expression in a sequence-dependent manner. These data demonstrate that viral transcripts can interact with host miRNAs to alter the outcomes of infection.

Adaptive truncation of the S gene in IBV during chicken embryo passaging plays a crucial role in its attenuation

Like all coronaviruses, infectious bronchitis virus, the causative agent of infectious bronchitis in chickens, exhibits a high mutation rate. Adaptive mutations that arise during the production of live attenuated vaccines against IBV often decrease virulence. The specific impact of these mutations on viral pathogenicity, however, has not been fully elucidated. In this study, we identified a mutation at the 3' end of the S gene in an IBV strain that was serially passaged in chicken embryos, and showed that this mutation resulted in a 9-aa truncation of the cytoplasmic tail (CT) of the S protein. This phenomenon of CT truncation has previously been observed in the production of attenuated vaccines against other coronaviruses such as the porcine epidemic diarrhea virus. We next discovered that the 9-aa truncation in the S protein CT resulted in the loss of the endoplasmic-reticulum-retention signal (KKSV). Rescue experiments with recombinant viruses confirmed that the deletion of the KKSV motif impaired the localization of the S protein to the endoplasmic-reticulum-Golgi intermediate compartment (ERGIC) and increased its expression on the cell surface. This significantly reduced the incorporation of the S protein into viral particles, impaired early subgenomic RNA and protein synthesis, and ultimately reduced viral invasion efficiency in CEK cells. In vivo experiments in chickens confirmed the reduced pathogenicity of the mutant IBV strains. Additionally, we showed that the adaptive mutation altered the TRS-B of ORF3 and impacted the transcriptional regulation of this gene. Our findings underscore the significance of this adaptive mutation in the attenuation of IBV infection and provide a novel strategy for the development of live attenuated IBV vaccines.

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.

Experimental and analytical pipeline for sub-genomic RNA landscape of coronavirus by Nanopore sequencer

Coronaviruses (CoVs), including severe acute respiratory syndrome coronavirus 2, can infect a variety of mammalian and avian hosts with significant medical and economic consequences. During the life cycle of CoV, a coordinated series of subgenomic RNAs, including canonical subgenomic messenger RNA and non-canonical defective viral genomes (DVGs), are generated with different biological implications. Studies that adopted the Nanopore sequencer (ONT) to investigate the landscape and dynamics of viral RNA subgenomic transcriptomes applied arbitrary bioinformatics parameters without justification or experimental validation. The current study used bovine coronavirus (BCoV), which can be performed under biosafety level 2 for library construction and experimental validation using traditional colony polymerase chain reaction and Sanger sequencing. Four different ONT protocols, including RNA direct and cDNA direct sequencing with or without exonuclease treatment, were used to generate RNA transcriptomic libraries from BCoV-infected cell lysates. Through rigorously examining the k-mer, gap size, segment size, and bin size, the optimal cutoffs for the bioinformatic pipeline were determined to remove the sequence noise while keeping the informative DVG reads. The sensitivity and specificity of identifying DVG reads using the proposed pipeline can reach 82.6% and 99.6% under the k-mer size cutoff of 15. Exonuclease treatment reduced the abundance of RNA transcripts; however, it was not necessary for future library preparation. Additional recovery of clipped BCoV nucleotide sequences with experimental validation expands the landscape of the CoV discontinuous RNA transcriptome, whose biological function requires future investigation. The results of this study provide the benchmarks for library construction and bioinformatic parameters for studying the discontinuous CoV RNA transcriptome.IMPORTANCEFunctional defective viral genomic RNA, containing all the cis-acting elements required for translation or replication, may play different roles in triggering cell innate immune signaling, interfering with the canonical subgenomic messenger RNA transcription/translation or assisting in establishing persistence infection. This study does not only provide benchmarks for library construction and bioinformatic parameters for studying the discontinuous coronavirus RNA transcriptome but also reveals the complexity of the bovine coronavirus transcriptome, whose functional assays will be critical in future studies.

Utilization of nanopore direct RNA sequencing to analyze viral RNA modifications

Modifications on viral RNAs (vRNAs), either genomic RNAs or RNA transcripts, have complex effects on the viral life cycle and cellular responses to viral infection. The advent of Oxford Nanopore Technologies Direct RNA Sequencing provides a new strategy for studying RNA modifications. To this end, multiple computational tools have been developed, but a systemic evaluation of their performance in mapping vRNA modifications is lacking. Here, 10 computational tools were tested using the Sindbis virus (SINV) RNAs isolated from infected mammalian (BHK-21) or mosquito (C6/36) cells, with in vitro-transcribed RNAs serving as modification-free control. Three single-mode approaches were shown to be inapplicable in the viral context, and three out of seven comparative methods required cutoff adjustments to reduce false-positive predictions. Utilizing optimized cutoffs, an integrated analysis of comparative tools suggested that the intersected predictions of Tombo_com and xPore were significantly enriched compared with the background. Consequently, a pipeline integrating Tombo_com and xPore was proposed for vRNA modification detection; the performance of which was supported by N6-methyladenosine prediction in severe acute respiratory syndrome coronavirus 2 RNAs using publicly available data. When applied to SINV RNAs, this pipeline revealed more intensive modifications in subgenomic RNAs than in genomic RNAs. Modified uridines were frequently identified, exhibiting substantive overlapping between vRNAs generated in different cell lines. On the other hand, the interpretation of other modifications remained unclear, underlining the limitations of the current computational tools despite their notable potential.IMPORTANCEComputational approaches utilizing Oxford Nanopore Technologies Direct RNA Sequencing data were almost exclusively designed to map eukaryotic epitranscriptomes. Therefore, extra caution must be exercised when using these tools to detect vRNA modifications, as in most cases, vRNA modification profiles should be regarded as unknown epitranscriptomes without prior knowledge. Here, we comprehensively evaluated the performance of 10 computational tools in detecting vRNA modification sites. All tested single-mode methods failed to differentiate native and in vitro-transcribed samples. Using optimized cutoff values, seven tested comparative tools generated very different predictions. An integrated analysis showed significant enrichment of Tombo_com and xPore predictions against the background. A pipeline for vRNA modification detection was proposed accordingly and applied to Sindbis virus RNAs. In conclusion, our study underscores the need for the careful application of computational tools to analyze viral epitranscriptomics. It also offers insights into alphaviral RNA modifications, although further validation is required.

SARS-CoV-2 Orphan Gene ORF10 Contributes to More Severe COVID-19 Disease

The orphan gene of SARS-CoV-2, ORF10, is the least studied gene in the virus responsible for the COVID-19 pandemic. Recent experimentation indicated ORF10 expression moderates innate immunity in vitro. However, whether ORF10 affects COVID-19 in humans remained unknown. We determine that the ORF10 sequence is identical to the Wuhan-Hu-1 ancestral haplotype in 95% of genomes across five variants of concern (VOC). Four ORF10 variants are associated with less virulent clinical outcomes in the human host: three of these affect ORF10 protein structure, one affects ORF10 RNA structural dynamics. RNA-Seq data from 2070 samples from diverse human cells and tissues reveals ORF10 accumulation is conditionally discordant from that of other SARS-CoV-2 transcripts. Expression of ORF10 in A549 and HEK293 cells perturbs immune-related gene expression networks, alters expression of the majority of mitochondrially-encoded genes of oxidative respiration, and leads to large shifts in levels of 14 newly-identified transcripts. We conclude ORF10 contributes to more severe COVID-19 clinical outcomes in the human host.

What do we know about the function of SARS-CoV-2 proteins?

The COVID-19 pandemic has highlighted the importance in the understanding of the biology of SARS-CoV-2. After more than two years since the first report of COVID-19, it remains crucial to continue studying how SARS-CoV-2 proteins interact with the host metabolism to cause COVID-19. In this review, we summarize the findings regarding the functions of the 16 non-structural, 6 accessory and 4 structural SARS-CoV-2 proteins. We place less emphasis on the spike protein, which has been the subject of several recent reviews. Furthermore, comprehensive reviews about COVID-19 therapeutic have been also published. Therefore, we do not delve into details on these topics; instead we direct the readers to those other reviews. To avoid confusions with what we know about proteins from other coronaviruses, we exclusively report findings that have been experimentally confirmed in SARS-CoV-2. We have identified host mechanisms that appear to be the primary targets of SARS-CoV-2 proteins, including gene expression and immune response pathways such as ribosome translation, JAK/STAT, RIG-1/MDA5 and NF-kβ pathways. Additionally, we emphasize the multiple functions exhibited by SARS-CoV-2 proteins, along with the limited information available for some of these proteins. Our aim with this review is to assist researchers and contribute to the ongoing comprehension of SARS-CoV-2's pathogenesis.

Unambiguous detection of SARS-CoV-2 subgenomic mRNAs with single-cell RNA sequencing

Single-cell RNA sequencing (scRNA-Seq) studies have provided critical insight into the pathogenesis of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of coronavirus disease 2019 (COVID-19). scRNA-Seq library preparation methods and data processing workflows are generally designed for the detection and quantification of eukaryotic host mRNAs and not viral RNAs. Here, we compare different scRNA-Seq library preparation methods for their ability to quantify and detect SARS-CoV-2 RNAs with a focus on subgenomic mRNAs (sgmRNAs). We show that compared to 10X Genomics Chromium Next GEM Single Cell 3' (10X 3') libraries or 10X Genomics Chromium Next GEM Single Cell V(D)J (10X 5') libraries sequenced with standard read configurations, 10X 5' libraries sequenced with an extended length read 1 (R1) that covers both cell barcode and transcript sequence (termed "10X 5' with extended R1") increase the number of unambiguous reads spanning leader-sgmRNA junction sites. We further present a data processing workflow, single-cell coronavirus sequencing (scCoVseq), which quantifies reads unambiguously assigned to viral sgmRNAs or viral genomic RNA (gRNA). We find that combining 10X 5' with extended R1 library preparation/sequencing and scCoVseq data processing maximizes the number of viral UMIs per cell quantified by scRNA-Seq. Corresponding sgmRNA expression levels are highly correlated with expression in matched bulk RNA-Seq data sets quantified with established tools for SARS-CoV-2 analysis. Using this scRNA-Seq approach, we find that SARS-CoV-2 gene expression is highly correlated across individual infected cells, which suggests that the proportion of viral sgmRNAs remains generally consistent throughout infection. Taken together, these results and corresponding data processing workflow enable robust quantification of coronavirus sgmRNA expression at single-cell resolution, thereby supporting high-resolution studies of viral RNA processes in individual cells. IMPORTANCE Single-cell RNA sequencing (scRNA-Seq) has emerged as a valuable tool to study host-virus interactions, especially for coronavirus disease 2019 (COVID-19). Here we compare the performance of different scRNA-Seq library preparation methods and sequencing strategies to detect SARS-CoV-2 RNAs and develop a data processing workflow to quantify unambiguous sequence reads derived from SARS-CoV-2 genomic RNA and subgenomic mRNAs. After establishing a workflow that maximizes the detection of SARS-CoV-2 subgenomic mRNAs, we explore patterns of SARS-CoV-2 gene expression across cells with variable levels of total viral RNA, assess host gene expression differences between infected and bystander cells, and identify non-canonical and lowly abundant SARS-CoV-2 RNAs. The sequencing and data processing strategies developed here can enhance studies of coronavirus RNA biology at single-cell resolution and thereby contribute to our understanding of viral pathogenesis.

Secure reversal of immune evasion from refractory NSCLC and highly contagious CoV-2 mutants by using 3D-engineered multifunctional biologics

There is an imperative choice to develop a secure feasible strategy to address evasion dynamics of refractory tumors and SARS-CoV-2-variants, while stem cell-based protocol may be more reliable as its unique ability for resetting multifunctional immunity to address progressive tumor and the constantly-evolving virus. In this study, spheroid-embryonoid stem cells from mature somatic cells were engineered as multifunctional biologics (3D-E/BSC) and inoculated in senile rhesus to identify secure potential against immune-evasion from viral-variants. Meanwhile, a cohort of eligible patients with stage IV NSCLC were approved for phase I clinical trials. Subsequently, long-lasting security and efficacy were validated by primate and clinical trials (p < 0.01) in that it could not only stimulate serological immunity, but also reset core immunity for hosts to address variant evasion after 3D-E/BSC withdrawal. Particularly, illustrated by single-cell evolving trajectory, 3D-E/BSC had securely reset senile thymus of aging hosts to remodel core immunity by rearranging naive rhythm to evolve TRGC2+/JCHAIN+NKT clusters to abolish tumoral and viral evasion dynamics with path-feedbacks of NSCLC and COVID-19 simultaneously activated, leading to continuous blockade of breakthrough infection of viral-mutants and long-term survival in one-third of terminal patients without adjuvant required. Our study may pioneer a practical multifunctional strategy to eliminate evasion of SARS-CoV-2 variants and refractory NSCLC so as for victims to restart a new life-equation.

Nanopore sequencing technology and its applications

Since the development of Sanger sequencing in 1977, sequencing technology has played a pivotal role in molecular biology research by enabling the interpretation of biological genetic codes. Today, nanopore sequencing is one of the leading third-generation sequencing technologies. With its long reads, portability, and low cost, nanopore sequencing is widely used in various scientific fields including epidemic prevention and control, disease diagnosis, and animal and plant breeding. Despite initial concerns about high error rates, continuous innovation in sequencing platforms and algorithm analysis technology has effectively addressed its accuracy. During the coronavirus disease (COVID-19) pandemic, nanopore sequencing played a critical role in detecting the severe acute respiratory syndrome coronavirus-2 virus genome and containing the pandemic. However, a lack of understanding of this technology may limit its popularization and application. Nanopore sequencing is poised to become the mainstream choice for preventing and controlling COVID-19 and future epidemics while creating value in other fields such as oncology and botany. This work introduces the contributions of nanopore sequencing during the COVID-19 pandemic to promote public understanding and its use in emerging outbreaks worldwide. We discuss its application in microbial detection, cancer genomes, and plant genomes and summarize strategies to improve its accuracy.

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.

High temporal resolution Nanopore sequencing dataset of SARS-CoV-2 and host cell RNAs

BACKGROUND

Recent studies have disclosed the genome, transcriptome, and epigenetic compositions of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and the effect of viral infection on gene expression of the host cells. It has been demonstrated that, besides the major canonical transcripts, the viral genome also codes for noncanonical RNA molecules. While the structural characterizations have revealed a detailed transcriptomic architecture of the virus, the kinetic studies provided poor and often misleading results on the dynamics of both the viral and host transcripts due to the low temporal resolution of the infection event and the low virus/cell ratio (multiplicity of infection [MOI] = 0.1) applied for the infection. It has never been tested whether the alteration in the host gene expressions is caused by aging of the cells or by the viral infection.

FINDINGS

In this study, we used Oxford Nanopore's direct cDNA and direct RNA sequencing methods for the generation of a high-coverage, high temporal resolution transcriptomic dataset of SARS-CoV-2 and of the primate host cells, using a high infection titer (MOI = 5). Sixteen sampling time points ranging from 1 to 96 hours with a varying time resolution and 3 biological replicates were used in the experiment. In addition, for each infected sample, corresponding noninfected samples were employed. The raw reads were mapped to the viral and to the host reference genomes, resulting in 49,661,499 mapped reads (54,62 Gbs). The genome of the viral isolate was also sequenced and phylogenetically classified.

CONCLUSIONS

This dataset can serve as a valuable resource for profiling the SARS-CoV-2 transcriptome dynamics, the virus-host interactions, and the RNA base modifications. Comparison of expression profiles of the host gene in the virally infected and in noninfected cells at different time points allows making a distinction between the effect of the aging of cells in culture and the viral infection. These data can provide useful information for potential novel gene annotations and can also be used for studying the currently available bioinformatics pipelines.

Systematic benchmarking of nanopore Q20+ kit in SARS-CoV-2 whole genome sequencing

Whole genome sequencing provides rapid insight into key information about the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), such as virus typing and key mutation site, and this information is important for precise prevention, control and tracing of coronavirus disease 2019 (COVID-19) outbreak in conjunction with the epidemiological information of the case. Nanopore sequencing is widely used around the world for its short sample-to-result time, simple experimental operation and long sequencing reads. However, because nanopore sequencing is a relatively new sequencing technology, many researchers still have doubts about its accuracy. The combination of the newly launched nanopore sequencing Q20+ kit (LSK112) and flow cell R10.4 is a qualitative improvement over the accuracy of the previous kits. In this study, we firstly used LSK112 kit with flow cell R10.4 to sequence the SARS-CoV-2 whole genome, and summarized the sequencing results of the combination of LSK112 kit and flow cell R10.4 for the 1200bp amplicons of SARS-CoV-2. We found that the proportion of sequences with an accuracy of more than 99% reached 30.1%, and the average sequence accuracy reached 98.34%, while the results of the original combination of LSK109 kit and flow cell R9.4.1 were 0.61% and 96.52%, respectively. The mutation site analysis showed that it was completely consistent with the final consensus sequence of next generation sequencing (NGS). The results showed that the combination of LSK112 kit and flow cell R10.4 allowed rapid whole-genome sequencing of SARS-CoV-2 without the need for verification of NGS.

Analyzing viral epitranscriptomes using nanopore direct RNA sequencing

RNA modifications are a common occurrence across all domains of life. Several chemical modifications, including N6-methyladenosine, have also been found in viral transcripts and viral RNA genomes. Some of the modifications increase the viral replication efficiency while also helping the virus to evade the host immune system. Nonetheless, there are numerous examples in which the host's RNA modification enzymes function as antiviral factors. Although established methods like MeRIP-seq and miCLIP can provide a transcriptome- wide overview of how viral RNA is modified, it is difficult to distinguish between the complex overlapping viral transcript isoforms using the short read-based techniques. Nanopore direct RNA sequencing (DRS) provides both long reads and direct signal readings, which may carry information about the modifications. Here, we describe a refined protocol for analyzing the RNA modifications in viral transcriptomes using nanopore technology.

The epitranscriptome of Vero cells infected with SARS-CoV-2 assessed by direct RNA sequencing reveals m6A pattern changes and DRACH motif biases in viral and cellular RNAs

The epitranscriptomics of the SARS-CoV-2 infected cell reveals its response to viral replication. Among various types of RNA nucleotide modifications, the m6A is the most common and is involved in several crucial processes of RNA intracellular location, maturation, half-life and translatability. This epitranscriptome contains a mixture of viral RNAs and cellular transcripts. In a previous study we presented the analysis of the SARS-CoV-2 RNA m6A methylation based on direct RNA sequencing and characterized DRACH motif mutations in different viral lineages. Here we present the analysis of the m6A transcript methylation of Vero cells (derived from African Green Monkeys) and Calu-3 cells (human) upon infection by SARS-CoV-2 using direct RNA sequencing data. Analysis of these data by nonparametric statistics and two computational methods (m6anet and EpiNano) show that m6A levels are higher in RNAs of infected cells. Functional enrichment analysis reveals increased m6A methylation of transcripts involved in translation, peptide and amine metabolism. This analysis allowed the identification of differentially methylated transcripts and m6A unique sites in the infected cell transcripts. Results here presented indicate that the cell response to viral infection not only changes the levels of mRNAs, as previously shown, but also its epitranscriptional pattern. Also, transcriptome-wide analysis shows strong nucleotide biases in DRACH motifs of cellular transcripts, both in Vero and Calu-3 cells, which use the signature GGACU whereas in viral RNAs the signature is GAACU. We hypothesize that the differences of DRACH motif biases, might force the convergent evolution of the viral genome resulting in better adaptation to target sequence preferences of writer, reader and eraser enzymes. To our knowledge, this is the first report on m6A epitranscriptome of the SARS-CoV-2 infected Vero cells by direct RNA sequencing, which is the sensu stricto RNA-seq.

Hydrophobic Alpha-Helical Short Peptides in Overlapping Reading Frames of the Coronavirus Genome

In this study, we show that the coronavirus (CoV) genome may encode many functional hydrophobic alpha-helical peptides (HAHPs) in overlapping reading frames of major coronaviral proteins throughout the entire viral genome. These HAHPs can theoretically be expressed from non-canonical sub-genomic (sg)RNAs that are synthesized in substantial amounts in infected cells. We selected and analyzed five and six HAHPs encoded in the S gene regions of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and Middle East respiratory syndrome coronavirus (MERS-CoV), respectively. Two and three HAHPs derived from SARS-CoV-2 and MERS-CoV, respectively, specifically interacted with both the SARS-CoV-2 and MERS-CoV S proteins and inhibited their membrane fusion activity. Furthermore, one of the SARS-CoV-2 HAHPs specifically inhibited viral RNA synthesis by accumulating at the site of viral RNA synthesis. Our data show that a group of HAHPs in the coronaviral genome potentially has a regulatory role in viral propagation.

A DNA-based non-infectious replicon system to study SARS-CoV-2 RNA synthesis

The coronavirus disease-2019 (COVID-19) pandemic caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has seriously affected public health around the world. In-depth studies on the pathogenic mechanisms of SARS-CoV-2 is urgently necessary for pandemic prevention. However, most laboratory studies on SARS-CoV-2 have to be carried out in bio-safety level 3 (BSL-3) laboratories, greatly restricting the progress of relevant experiments. In this study, we used a bacterial artificial chromosome (BAC) method to assemble a SARS-CoV-2 replication and transcription system in Vero E6 cells without virion envelope formation, thus avoiding the risk of coronavirus exposure. Furthermore, an improved real-time quantitative reverse transcription PCR (RT-qPCR) approach was used to distinguish the replication of full-length replicon RNAs and transcription of subgenomic RNAs (sgRNAs). Using the SARS-CoV-2 replicon, we demonstrated that the nucleocapsid (N) protein of SARS-CoV-2 facilitates the transcription of sgRNAs in the discontinuous synthesis process. Moreover, two high-frequency mutants of N protein, R203K and S194L, can obviously enhance the transcription level of the replicon, hinting that these mutations likely allow SARS-CoV-2 to spread and reproduce more quickly. In addition, remdesivir and chloroquine, two well-known drugs demonstrated to be effective against coronavirus in previous studies, also inhibited the transcription of our replicon, indicating the potential applications of this system in antiviral drug discovery. Overall, we developed a bio-safe and valuable replicon system of SARS-CoV-2 that is useful to study the mechanisms of viral RNA synthesis and has potential in novel antiviral drug screening.

Nanopore-Based Detection of Viral RNA Modifications

The chemical modification of ribonucleotides plays an integral role in the biology of diverse viruses and their eukaryotic host cells. Mapping the precise identity, location, and abundance of modified ribonucleotides remains a key goal of many studies aimed at characterizing the function and importance of a given modification. While mapping of specific RNA modifications through short-read sequencing approaches has powered a wealth of new discoveries in the past decade, this approach is limited by inherent biases and an absence of linkage information. Moreover, in viral contexts, the challenge is increased due to the compact nature of viral genomes giving rise to many overlapping transcript isoforms that cannot be adequately resolved using short-read sequencing approaches. The recent emergence of nanopore sequencing, specifically the ability to directly sequence native RNAs from virus-infected host cells, provides not just a new methodology for mapping modified ribonucleotides but also a new conceptual framework for what can be derived from the resulting sequencing data. In this minireview, we provide a detailed overview of how nanopore direct RNA sequencing works, the computational approaches applied to identify modified ribonucleotides, and the core concepts underlying both. We further highlight recent studies that have applied this approach to interrogating viral biology and finish by discussing key experimental considerations and how we predict that these methodologies will continue to evolve.

A Unique Robust Dual-Promoter-Driven and Dual-Reporter-Expressing SARS-CoV-2 Replicon: Construction and Characterization

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2, SARS2) remains a great global health threat and demands identification of more effective and SARS2-targeted antiviral drugs, even with successful development of anti-SARS2 vaccines. Viral replicons have proven to be a rapid, safe, and readily scalable platform for high-throughput screening, identification, and evaluation of antiviral drugs against positive-stranded RNA viruses. In the study, we report a unique robust HIV long terminal repeat (LTR)/T7 dual-promoter-driven and dual-reporter firefly luciferase (fLuc) and green fluorescent protein (GFP)-expressing SARS2 replicon. The genomic organization of the replicon was designed with quite a few features that were to ensure the replication fidelity of the replicon, to maximize the expression of the full-length replicon, and to offer the monitoring flexibility of the replicon replication. We showed the success of the construction of the replicon and expression of reporter genes fLuc and GFP and SARS structural N from the replicon DNA or the RNA that was in vitro transcribed from the replicon DNA. We also showed detection of the negative-stranded genomic RNA (gRNA) and subgenomic RNA (sgRNA) intermediates, a hallmark of replication of positive-stranded RNA viruses from the replicon. Lastly, we showed that expression of the reporter genes, N gene, gRNA, and sgRNA from the replicon was sensitive to inhibition by Remdesivir. Taken together, our results support use of the replicon for identification of anti-SARS2 drugs and development of new anti-SARS strategies targeted at the step of virus replication.

Long-Read RNA Sequencing Identifies Polyadenylation Elongation and Differential Transcript Usage of Host Transcripts During SARS-CoV-2 In Vitro Infection

Better methods to interrogate host-pathogen interactions during Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) infections are imperative to help understand and prevent this disease. Here we implemented RNA-sequencing (RNA-seq) using Oxford Nanopore Technologies (ONT) long-reads to measure differential host gene expression, transcript polyadenylation and isoform usage within various epithelial cell lines permissive and non-permissive for SARS-CoV-2 infection. SARS-CoV-2-infected and mock-infected Vero (African green monkey kidney epithelial cells), Calu-3 (human lung adenocarcinoma epithelial cells), Caco-2 (human colorectal adenocarcinoma epithelial cells) and A549 (human lung carcinoma epithelial cells) were analyzed over time (0, 2, 24, 48 hours). Differential polyadenylation was found to occur in both infected Calu-3 and Vero cells during a late time point (48 hpi), with Gene Ontology (GO) terms such as viral transcription and translation shown to be significantly enriched in Calu-3 data. Poly(A) tails showed increased lengths in the majority of the differentially polyadenylated transcripts in Calu-3 and Vero cell lines (up to ~101 nt in mean poly(A) length, padj = 0.029). Of these genes, ribosomal protein genes such as RPS4X and RPS6 also showed downregulation in expression levels, suggesting the importance of ribosomal protein genes during infection. Furthermore, differential transcript usage was identified in Caco-2, Calu-3 and Vero cells, including transcripts of genes such as GSDMB and KPNA2, which have previously been implicated in SARS-CoV-2 infections. Overall, these results highlight the potential role of differential polyadenylation and transcript usage in host immune response or viral manipulation of host mechanisms during infection, and therefore, showcase the value of long-read sequencing in identifying less-explored host responses to disease.

Control of animal virus replication by RNA adenosine methylation

Methylation at the N⁶-position of either adenosine (m⁶A) or 2'-O-methyladenosine (m⁶Am) represents two of the most abundant internal modifications of coding and non-coding RNAs, influencing their maturation, stability and function. Additionally, although less abundant and less well-studied, monomethylation at the N¹-position (m¹A) can have profound effects on RNA folding. It has been known for several decades that RNAs produced by both DNA and RNA viruses can be m⁶A/m⁶Am modified and the list continues to broaden through advances in detection technologies and identification of the relevant methyltransferases. Recent studies have uncovered varied mechanisms used by viruses to manipulate the m⁶A pathway in particular, either to enhance virus replication or to antagonize host antiviral defenses. As such, RNA modifications represent an important frontier of exploration in the broader realm of virus-host interactions, and this new knowledge already suggests exciting opportunities for therapeutic intervention. In this review we summarize the principal mechanisms by which m⁶A/m⁶Am can promote or hinder viral replication, describe how the pathway is actively manipulated by biomedically important viruses, and highlight some remaining gaps in understanding how adenosine methylation of RNA controls viral replication and pathogenesis.

Experimental Models of COVID-19

COVID-19 is the most consequential pandemic of the 21st century. Since the earliest stage of the 2019-2020 epidemic, animal models have been useful in understanding the etiopathogenesis of SARS-CoV-2 infection and rapid development of vaccines/drugs to prevent, treat or eradicate SARS-CoV-2 infection. Early SARS-CoV-1 research using immortalized in-vitro cell lines have aided in understanding different cells and receptors needed for SARS-CoV-2 infection and, due to their ability to be easily manipulated, continue to broaden our understanding of COVID-19 disease in in-vivo models. The scientific community determined animal models as the most useful models which could demonstrate viral infection, replication, transmission, and spectrum of illness as seen in human populations. Until now, there have not been well-described animal models of SARS-CoV-2 infection although transgenic mouse models (i.e. mice with humanized ACE2 receptors with humanized receptors) have been proposed. Additionally, there are only limited facilities (Biosafety level 3 laboratories) available to contribute research to aid in eventually exterminating SARS-CoV-2 infection around the world. This review summarizes the most successful animal models of SARS-CoV-2 infection including studies in Non-Human Primates (NHPs) which were found to be susceptible to infection and transmitted the virus similarly to humans (e.g., Rhesus macaques, Cynomolgus, and African Green Monkeys), and animal models that do not require Biosafety level 3 laboratories (e.g., Mouse Hepatitis Virus models of COVID-19, Ferret model, Syrian Hamster model). Balancing safety, mimicking human COVID-19 and robustness of the animal model, the Murine Hepatitis Virus-1 Murine model currently represents the most optimal model for SARS-CoV-2/COVID19 research. Exploring future animal models will aid researchers/scientists in discovering the mechanisms of SARS-CoV-2 infection and in identifying therapies to prevent or treat COVID-19.

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.

Isoform Age - Splice Isoform Profiling Using Long-Read Technologies

Alternative splicing (AS) of RNA is a key mechanism that results in the expression of multiple transcript isoforms from single genes and leads to an increase in the complexity of both the transcriptome and proteome. Regulation of AS is critical for the correct functioning of many biological pathways, while disruption of AS can be directly pathogenic in diseases such as cancer or cause risk for complex disorders. Current short-read sequencing technologies achieve high read depth but are limited in their ability to resolve complex isoforms. In this review we examine how long-read sequencing (LRS) technologies can address this challenge by covering the entire RNA sequence in a single read and thereby distinguish isoform changes that could impact RNA regulation or protein function. Coupling LRS with technologies such as single cell sequencing, targeted sequencing and spatial transcriptomics is producing a rapidly expanding suite of technological approaches to profile alternative splicing at the isoform level with unprecedented detail. In addition, integrating LRS with genotype now allows the impact of genetic variation on isoform expression to be determined. Recent results demonstrate the potential of these techniques to elucidate the landscape of splicing, including in tissues such as the brain where AS is particularly prevalent. Finally, we also discuss how AS can impact protein function, potentially leading to novel therapeutic targets for a range of diseases.