Direct RNA sequencing of Respiratory Syncytial Virus infected human cells generates a detailed overview of RSV polycistronic mRNA and transcript abundance

To characterize species of viral mRNA transcripts generated during respiratory syncytial virus (RSV) infection, human fibroblast-like MRC5 lung cells were infected with subgroup A RSV for 6, 16 and 24 hours. Total RNA was harvested and polyadenylated mRNA was enriched and sequenced by direct RNA sequencing on an Oxford nanopore device. This yielded over 150,000 direct mRNA transcript reads which were mapped to the viral genome and analysed to determine relative mRNA levels of viral genes using our in-house ORF-centric pipeline. We were also able to examine frequencies with which polycistronic readthrough mRNAs were generated and to assess the length of the polyadenylated tails for each group of transcripts. We show that there is a general but non-linear decline in gene transcript abundance across the viral genome, as predicted by the model of RSV gene transcription. However, the decline in transcript abundance is not consistent. We show that the polyadenylate tails generated by the viral polymerase are similar in length to those generated by the host cells polyadenylation machinery and broadly declined in length for most transcripts as infection progressed. Finally, we observed that the steady state abundance of transcripts with very short polyadenylate tails is much less for N, SH and G transcripts compared to NS1, NS2, P, M, F and M2 which may reflect differences in mRNA stability and/or translation rates.