Direct RNA sequencing reveals m6A modifications and isoform changes in SARS-CoV-2-infected HEK cells

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection triggers complex host responses, including alterations in RNA transcription and modification. Understanding these changes is crucial for elucidating viral pathogenesis and identifying potential therapeutic targets. We used direct RNA sequencing to comprehensively profile the transcriptomic and epitranscriptomic landscapes of human HEK-AT cells infected with SARS-CoV-2 at 8 h post-infection, compared to mock controls. We analysed viral and host transcriptomes, focusing on gene and transcript expression, isoform usage and RNA m6A modifications. Viral RNA sequencing reads showed 3' end-biassed coverage indicative of subgenomic RNA synthesis, with high expression of N gene subgenomic RNA reads. Sixteen m6A modification sites were consistently identified in the viral genome, primarily within the ORF1ab and S genes. In the human transcriptome, we found 254 positions with significantly altered m6A modification rates, with 119 showing decreased modification and 135 showing increased modification in infected cells. Genes with decreased m6A modifications were enriched in the neurotrophin signalling pathway. Transcript-level analysis identified 19 upregulated and 12 downregulated transcripts. Notably, transcript discovery and quantification revealed a novel isoform of the HIST1H2BK gene, which was significantly more expressed in infected cells compared to mock controls. Isoform switching analysis revealed 24 significant switches involving 21 genes, implicating mitochondrial reprogramming and immune-related pathways. In conclusion, this study provides a detailed, direct RNA sequencing-based characterization of host-virus RNA interactions, revealing key insights into SARS-CoV-2 infection mechanisms and potential therapeutic targets.