In silico approach for the identification of tRNA-derived small non-coding RNAs in SARS-CoV infection

tsRNAs (tRNA-derived small non-coding RNAs), including tRNA halves (tiRNAs) and tRNA fragments (tRFs), have been implicated in some viral infections, such as respiratory viral infections. However, their involvement in SARS-CoV infection is completely unknown. A comprehensive analysis was performed to determine tsRNA populations in a mouse model of SARS-CoV-infected samples containing the wild-type and attenuated viruses. Data from the Gene Expression Omnibus (GEO) dataset at NCBI (accession ID GSE90624 ) was used for this study. A count matrix was generated for the tRNAs. Differentially expressed tRNAs, followed by tsRNAs derived from each significant tRNAs at different conditions and time points between the two groups WT(SARS-CoV-MA15-WT) vs Mock and ΔE (SARS-CoV-MA15-ΔE) vs Mock were identified. Notably, significantly differentially expressed tRNAs at 2dpi but not at 4dpi. The tsRNAs originating from differentially expressed tRNAs across all the samples belonging to each condition (WT, ΔE, and Mock) were identified. Intriguingly, tRFs (tRNA-derived RNA fragments) exhibited higher levels compared to tiRNAs (tRNA-derived stress-induced RNAs) across all samples associated with WT SARS-CoV strain compared to ΔE and mock-infected samples. This discrepancy suggests a non-random formation of tsRNAs, hinting at a possible involvement of tsRNAs in SARS-CoV viral infection.