Next Article in Journal
The Virus-Induced Transcription Factor SHE1 Interacts with and Regulates Expression of the Inhibitor of Virus Replication (IVR) in N Gene Tobacco
Previous Article in Journal
A Polydnavirus Protein Tyrosine Phosphatase Negatively Regulates the Host Phenoloxidase Pathway
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Parent tRNA Modification Status Determines the Induction of Functional tRNA-Derived RNA by Respiratory Syncytial Virus Infection

1
Department of Pediatrics, University of Texas Medical Branch, Galveston, TX 77555, USA
2
Department of Biological and Chemical Science, New York Institute of Technology, New York, NY 10023, USA
3
Institute of Translational Science, University of Texas Medical Branch, Galveston, TX 77555, USA
4
Institute for Human Infections Immunity, University of Texas Medical Branch, Galveston, TX 77555, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Viruses 2023, 15(1), 57; https://doi.org/10.3390/v15010057
Submission received: 30 October 2022 / Revised: 16 December 2022 / Accepted: 21 December 2022 / Published: 24 December 2022
(This article belongs to the Section Viral Immunology, Vaccines, and Antivirals)

Abstract

:
tRNA-derived RNA fragments (tRFs) are a recently discovered family of small noncoding RNAs (sncRNAs). We previously reported that respiratory syncytial virus (RSV) infection induces functional tRFs, which are derived from a limited subset of parent tRNAs, in airway epithelial cells. Such induction is also observed in nasopharyngeal wash samples from RSV patients and correlates to RSV genome copies, suggesting a clinical significance of tRFs in RSV infection. This work also investigates whether the modification of parent tRNAs is changed by RSV to induce tRFs, using one of the most inducible tRFs as a model. We discovered that RSV infection changed the methylation modification of adenine at position 57 in tRNA glutamic acid, with a codon of CTC (tRNA-GluCTC), and the change is essential for its cleavage. AlkB homolog 1, a previously reported tRNA demethylase, appears to remove methyladenine from tRNA-GluCTC, prompting the subsequent production of tRFs from the 5′-end of tRNA-GluCTC, a regulator of RSV replication. This study demonstrates for the first time the importance of post-transcriptional modification of tRNAs in tRF biogenesis following RSV infection, providing critical insights for antiviral strategy development.

1. Introduction

Respiratory syncytial virus (RSV) is a negative-sense single-stranded non-segmented RNA virus. As an orthopneumovirus, it belongs to the Pneumoviridae family and is the most significant cause of lower respiratory tract infections (LRTIs) in the pediatric population [1,2]. Other than palivizumab, which prevents RSV infection in high-risk infants, there is no vaccine or virus-specific pharmaceutic treatment for RSV [3], demonstrating the need to expand our current knowledge on viral replication mechanisms.
Noncoding RNAs (ncRNAs) comprise most of the human genome (as compared with protein-coding RNAs) [4,5], and their functions in physiological and pathological settings are increasingly recognized. tRNA-derived RNA fragments (tRFs) are a recently discovered category of ncRNAs whose importance in viral infectious diseases is becoming more evident [6,7,8,9,10]. We have recently shown that RSV infection induces a significant increase in tRFs. These induced tRFs are not byproducts of tRNA degradation but functional molecules because (1) the induction is virus-specific, (2) only a selective subset of parent tRNAs are cleaved in a tRNA structure-dependent manner, and (3) some are functionally important in viral replication and regulate genes using mechanisms different from miRNAs [8,11]. We also found that ribonuclease angiogenin (ANG) controls most cleaved tRNAs in RSV while ElaC ribonuclease Z 2 (ELAC2) is essential for the cleavage of tRNA-GlnCTG [8,11,12]. However, the cause of tRNAs being prone to cleavage is not fully understood.
tRF derived from the 5′-end of tRNA-encoding GluCTC (tRF5-GluCTC) is heavily impacted by RSV. It has also shown its importance in other diseases, including cancer and neurodegenerative disease [13,14,15], and is also inducible by SARS-CoV-2 infection. Therefore, we used tRF5-GluCTC as a start to study the biogenesis mechanism of tRF5s. Furthermore, since tRNAs are highly modified molecules, and the modification status of RNAs plays a role in RNA stability [16,17,18], we investigated whether RSV infection ignites the cleavage of parent tRNA-GluCTC by changing its modification status. Herein, our cutting-edge de novo mass spectrometric (MS)-based sequencing methods, known as MS-Seq, enabled us to directly detect a decrease in adenine (A) methylation at position 57 in tRNA-GluCTC. AlkB homolog 1 (ALKBH1), a reported tRNA demethylase with a demethylation preference for methyladenine present on a T loop structure of tRNAs [19], played a critical role in tRNA-GluCTC cleavage following the RSV infection, demonstrating that ALKBH1-mediated A demethylation is critical in the induction of tRF5-GluCTC. We also found a novel modification change on two pseudoruidine residues by RSV, awaiting the identification of associated modification enzymes. Since tRF5-GluCTC promotes RSV replication [8,13], the importance of tRNA modification to tRF biogenesis provides new insights for viral replication control.

2. Materials and Methods

2.1. Nasopharyngeal Secretions (NPS) from Patients

Samples of NPS were collected by gentle suction of the nasopharynx from hospitalized children who were enrolled in a study on the pathogenesis of RSV bronchiolitis under protocol #03-117 (approved by the Institutional Review Board of UTMB). The Luminex NxTAG™ Respiratory Pathogen Panel (RPP, Luminex Molecular Diagnostics, Austin, TX, USA) was used to detect RSV and other viruses in NPS samples. Only RSV-positive NPS (n = 9) were used as positive samples and pathogen-negative NPS (n = 8) as controls. The information on the patient’s gender and age is shown in Table 1.

2.2. Cells and Virus

A549 (human alveolar type II-like epithelial) and HEp-2 (human epithelial type 2) were from the ATCC (Manassas, VA) and maintained as previously described [20]. RSV long strain was propagated in HEp-2 cells at 37 °C, and the virus purification was done by sucrose gradient as described previously [20,21]. RSV titer was determined by immunostaining in HEp-2 cells using polyclonal biotin-conjugated goat anti-RSV antibody (Bio-Rad, Hercules, CA, USA) and streptavidin peroxidase polymer, as previously described [20,21].

2.3. RNA Preparation

To prepare RNAs from patient samples, RNAs from 300 µL of NPS, after adding cel-miR-39 as reference controls for sample preparation, were harvested using mirVanaTM Protein and RNA Isolation System from Invitrogen (Waltham, MA, USA, Cat# AM1556), following the manufacturer’s protocol. To harvest RNAs from cell cultures, we used Trizol from Invitrogen (Cat# 15596-026) to extract RNAs as we previously described in [20,22,23].

2.4. qRT-PCR

Trizol-extracted RNAs or RNAs from patient NPS were used to examine viral replication by qRT-PCR, as previously described for RSV or its family member, human metapeumovirus [20,22,23]. In brief, to quantify viral genomic copies of RSV, we used primers to assess the genomic part between the N and P regions. First-strand cDNA was transcribed with a P-specific primer: 5′-CGTCTCAGCCAATCCCTGGTGATTATGAGTAATTAAAAAATGGGACAAG-3′. The underlined letters indicate the chloramphenicol resistance (Cmr) tag sequence, which was designed as part of the assay to prevent self-priming exhibited by viral RNA [24]. QPCRs were performed using the forward primer 5′-CGTCTCAGCCAATCCCTGG-3′ and the reverse primer 5′-GCTTCATTACCCATGAAAAGAATATC-3′.
To quantify NPS tRF expression, we used modified qRT-PCR that can differentiate tRF signals from its parent tRNA, as previously described in [9]. In brief, RNAs were treated with T4 PNK and then ligated to a 3′-RNA linker using T4 RNA ligase from Thermo Fisher Scientific (Waltham, MA, USA). The product was used as a template for RT with a linker-specific reverse primer using TaqMan Reverse Transcription Reagents from Thermo Fisher Scientific. The cDNA was subjected to SYBR Green qPCR using iTaqTM Universal SYBR Green Supermix kit from Bio-Rad (Hercules, CA, USA) and primers specific to the 5′-end of tRFs and RNA linker. U6 and externally added cel-miR-39 were used for normalization.

2.5. tRNA-GluCTC Sample Preparation

Ten T-150 flasks of confluent A549 cells were infected with RSV at a multiplicity of infection (MOI) of 1. An equivalent amount of 30% sucrose solution was added to uninfected A549, which was used as a control. At 15 h post-infection (p.i.), total cellular RNAs were harvested by Trizol. The RNA pellets were then suspended in 10 mL of 5 × SSC buffer (diluted from 20 × SSC buffer, Corning Inc., Life Sciences, Acton, MA, USA). To purify tRNA-GluCTC, biotinylated DNA oligos (which are complementary to tRF5-GluCTC) were added to the RNA suspension. The mixture was then heated to 95 °C for 3 min and subsequently incubated at 65 °C for 15 min, followed by another incubation at room temperature for 3 h. The DNA-tRNA complex was then pulled down by adding streptavidin agarose beads (Thermo Scientific, Waltham, MA, USA), which were pre-washed twice in 5 × SSC buffer, followed by rotation at 4 °C overnight. The beads were washed once with 1 × SSC buffer and three times with 0.1 × SSC buffer before DNase I treatment to degrade biotinylated-DNA oligos. The released RNAs from the DNA-RNA complex were then extracted using Trizol LS reagents, followed by 15% denaturing polyacrylamide gel. The tRNA bands were cut and followed by PAGE gel elution using ZR small-RNATM PAGE Recovery Kit (Zymo Research, Orange, CA, USA). The purified tRNAs were then under acid hydrolysis degradation and lyophilization as previously described for RNA 2D-HELS MS sequencing [25,26,27].

2.6. tRNA-GluCTC Sequencing

tRNAs are enriched with chemical modifications, which are dynamic and sensitive to multiple post-transcriptional regulatory mechanisms. MS, in particular combined liquid chromatography and tandem MS spectrometry (LC-MS/MS), is recognized as the ‘gold standard’ for RNA modification analysis. However, it involves the full digestion of RNAs to component nucleotides and cannot provide full location or modification information. Recent works by Zhang et al. have developed a series of novel MS-based sequencing methods, referred to as MassSpec-Seq, that can simultaneously sequence RNA directly and quantify different RNA modifications with single-nucleotide stoichiometric precision through a complete set of MS ladders obtained from random RNA cuts and sequence mapping [25,26,28]. This method was further improved by an MS ladder complementation sequencing approach (MLC-Seq) [29]. By using MLC-Seq on tRNA-GluCTC obtained from human mock-infected or RSV-infected A549 cells, we tracked the location and stoichiometric changes of tRNA-GluCTC nucleotides caused by RSV infection. In brief, samples were analyzed using previously published protocols with minor modifications [29], and data were collected using an Orbitrap Exploris 240 mass spectrometer (ThermoFisher Scientific, Bremen, Germany) coupled to a Vanquish Horizon UHPLC. These experiments were performed using a DNAPac reversed phase (RP) column (2.1 mm × 50 mm, ThermoFisher Scientific) with 2% HFIP and 0.1% DIPEA as eluent A and methanol, 0.075% HFIP, and 0.0375% DIPEA as eluent B. A gradient of 20% to 80% B over 6.7 min was used for the analysis of intact RNA samples, while acid-degraded samples used a gradient from 15% to 35% B over 20 min.

2.7. NGS Sequencing of tRNA-GluCTC Sample

Pulled-down RNAs were delivered to the Genomics Core of UTMB for small RNA sequencing. RNAs were treated with T4 polynucleotide kinase (T4PNK) to convert the 3’ cyclic phosphate group into a hydroxyl group, as described in [12,15], before small RNA libraries were created using the New England Biolabs small RNA library protocol (New England Biolabs, Ipswich, MA, USA). Library construction used a two-step ligation process to create templates compatible with Illumina-based next-generation sequence (NGS) analysis. Where appropriate, RNA samples were quantified using a Qubit™ fluorometric assay (Thermo Fisher Scientific). RNA quality was assessed using a 2100 Bioanalyzer and the RNA 6000 Pico LabChip ® (Agilent Technologies, Santa Clara, CA, USA). Library creation was conducted by sequentially adding a 3′ adapter sequence followed by a 5′ adapter sequence. A cDNA copy was then synthesized using ProtoScript® II reverse transcriptase (New England Biolabs) and a primer complementary to a segment of the 3′ adapter. The template population was amplified in 15 cycles (94 °C for 30 s; 62 °C for 30 s; 70 °C for 30 s). The libraries were not size-selected, and all NGS libraries were indexed. The final concentration of all NGS libraries was determined using a Qubit™ fluorometric assay, and the DNA fragment size of each library was assessed using a DNA 1000 high-sensitivity chip and a 2100 Bioanalyzer (Agilent Technologies). Sequence analysis (2  ×  50 bases) was performed on an Illumina Hi-Seq 1500 using the TruSeq SBS v3 kit. Sequence counts of 10 or more after adaptor sequence removal were used for further classification, and small RNAs were mapped using Novoalign software as previously described [8].

2.8. Northern Blot (NB)

Northern hybridization for tRF/tRNA was performed as we previously described [8,30,31]. To summarize the process: RNAs, pulled down by antisense DNA oligos, were separated in 15% denaturing polyacrylamide gel with 7 mol/L urea and then transferred to a positively charged nylon membrane (Amersham Biosciences, NJ, USA). The membrane was hybridized with a 32P-labeled DNA probe inversely complementary to the tRF of interest in ULTRAhyb-Oligo solution (Life Technologies, NY, USA), followed by washing according to the manufacturer’s instructions.

2.9. RNA Interference

Transfections of siRNA into A549 cells were carried out at a final concentration of 100 nM, targeting ALKBH1 or scrambled negative control (Sigma, Woodland, TX, USA), using A549 transfection reagent Lipofectamine 2000 (ThermoFisher Scientific, Waltham, MA, USA) according to the manufacturer’s recommendations. After 24 h, A549 cells were mock infected or infected with RSV for 15 h at an MOI of 1.

2.10. Statistical Analysis

Experimental results were analyzed using Graphpad Prism 5 software. Patient group comparison was done by non-parametric statistics methods. Specifically, an unpaired two-tailed Mann-Whitney U test was used for the comparison of two independent groups. Single, double, and triple asterisks represent a p-value of 0.05, 0.01, and 0.001, respectively. Means ± standard errors (SE) are shown. Correlation analyses were performed using Spearman’s rank correlation test, and Spearman’s rank correlation coefficient (RS) was used to determine correlations. A p-value of less than 0.05 was considered significant. For the statistical analysis to determine the impact of ALKBH1 on RSV replication and tRF induction in cells, the statistical significance was determined using analysis of variance (ANOVA). A p-value of less than 0.05 was considered significant. Means ± SE are shown.

3. Results

3.1. tRF Expression in Nasopharyngeal Swab (NPS) Samples

The induction of functional tRF5s from a limited set of parent tRNAs by RSV has been demonstrated in airway epithelial cells [8,11]. The present study compared the tRF abundance in NPS collected from RSV-infected or virus-negative children to investigate whether such an induction also occurs in the course of naturally acquired infections. Table 1 and Table 2 list the patients’ information and qRT-PCR primers to measure tRFs in NPSs. As shown in Figure 1A–C, the expression of tRF5-CysGCA, tRF5-GlyCCC, and Gly5-GluCTC in NPSs was significantly higher in RSV-positive NPS samples than in virus-negative ones, consistent with RSV-changed tRFs in infected cells [8]. Such an increase was not observed for tRF5-HisGTG. We also found that tRF5-GluCTC expression correlated to the genome copies of RSV (Figure 1D), which is consistent with the reported function of tRF5-GluCTC in promoting RSV replication [8,13]. These preliminary data suggest a possible detrimental role of tRF in RSV infection.

3.2. Modifications of tRNA-GluCTC

tRF induction by RSV is a tightly regulated process, as the infection only leads to the cleavage of one or very few isoforms among tRNA isoacceptors (tRNAs that have different anticodons but still carry the same amino acid) and isodecoders (tRNA genes that share the same anticodon but have different body sequences). tRNAs are enriched with chemical modifications, and their modification status plays a vital role in determining their interaction with ribonuclease and associated cleavage [32,33]. Therefore, the selective cleavage may result from the modification changes of parent tRNAs by RSV. Since tRF5-GluCTC is one of the most abundantly induced tRFs by RSV and also functionally important to RSV replication [8,13] and showed a correlation with viral genome copies in clinical samples (Figure 1D), the modifications of tRF5-GluCTC’s parent tRNA were investigated before and after RSV infection. To achieve a high resolution on modification identification, the parent tRNA of tRF5-GluCTC was enriched as shown in Figure 2. tRNA-GluCTC was incubated with biotinylated antisense oligo (Figure 2A) and subsequently pulled down by streptavidin beads, which was followed by gel separation and purification (Figure 2B). Figure 2C shows that tRNA-GluCTC from control and infected cells was successfully enriched, as the bands of purified samples can be detected by the probe against tRNA-GluCTC (middle panel) but not by the probe against tRNA-GlyGCC (right panel).
The use of MLC-Seq on tRNA-GluCTC revealed a change in the methylation level of adenine at position 57 of the tRNA-GluCTC. About half of the tRNA-GluCTC stayed methylated in A57 for mock-infected cells (Figure 3B); this decreased to 11.5% after the infection (Figure 3C). The detailed measurement of methylation levels at position A57 in patient samples is provided in (Supplementary Excel). Additionally, a novel modification was discovered at positions 19 and 20 of tRNA-GluCTC. Both sites in mock-infected cells were identified to be 100% peudourantine (D). After infection, for the mass of 33% of D19 and 71% of D20 decreased by 1 Da. This was presumably caused by 3,4-dihydrocytidine (C′) at position 16, as described in Zhang et al. [29] (Figure 3D). Some methylation modifications were also observed on C and U. However, they were not changed by RSV infection.

3.3. Impact of Methylation Change on RSV Replication

The demethylation of A57 by RSV may play an important role in tRF5-GluCTC induction, as tRNA cleavage has been shown to be sensitive to this modification [34]. The change is thought to impact tRNA structure, stability, and function, which in turn can have a wide range of effects at the cellular level [35]. Of these modifications, m1A is common and has a major impact on tRNA stability [36]. AlkB Homolog 1 (ALKBH1) has been reported to be a tRNA demethylase by removing N(1)-methyladenine from various tRNAs, with a preference for N(1)-methyladenine at position 57/58 (m1A) present on a stem-loop structure of tRNAs. Therefore, ALKBH1-specific siRNAs were used first to suppress its expression and investigate its role in tRF5-GluCTC induction. In scrambled siRNA-treated cells, tRF5-GluCTC was detected in RSV-infected cells but not in mock-infected cells, indicating that the RSV infection induced tRF5-GluCTC (Figure 4A, first two lanes). This induction was suppressed by ALKBH1 knockdown (Figure 4A, last two lanes), suggesting the importance of ALKBH1 in tRNA-GluCTC cleavage. The effectiveness of ALKBH1-specific siRNA in silencing ALKBH1 expression was confirmed by Western blot using the antibody specifically against ALKBH1 (Cat#: ab126596, Abcam, Boston, MA) (Figure 4B,C). Suppression of ALKBH1 expression also led to decreased viral proteins (Figure 4B,C) and virus replication (Figure 4D). To investigate whether RSV uses ALKBH1 to control replication via controlling ALKBH1 expression, we measured ALKBH1 expression in the context of RSV infection in both cytosol and nuclear compartments. As shown in Figure 4E, the expression of ALKBH1 was not changed by RSV infection, suggesting that RSV-altered ALKBH1 activity affects the replication. This will be investigated in the future. These results are consistent with our previous findings showing the promotion of RSV replication by tRF5-GluCTC.

4. Discussion

tRFs are recently discovered sncRNAs and have quickly been found to be relevant to many diseases, including neurodegenerative diseases, cancers, and infectious diseases [7,9,15,37,38,39,40,41,42]. Since our first publication, which demonstrated that the most RSV-impacted sncRNAs in airway epithelial cells belong to tRFs [8], our team has investigated the functions of induced tRFs [8,11,12], identified the ribonucleases responsible for the tRNA cleavage [8,11,12], determined the platform through which tRFs carry out their gene regulatory functions [12], and outlined the flow work of tRF target identification [13]. The present study demonstrates the possible relevance of tRFs in RSV infection by showing tRF induction in NPS samples of infected patients and/or the direct correlation between expression levels of tRF5-GluCTC-CTC and RSV genome copies (Figure 1). In addition to tRF5-GluCTC, we also checked the expression of tRF5-GlyGCC-3, a tRF that is strongly induced by RSV in cells (~30 folds) [8]. We found that its expression levels also correlated with RSV replication, further supporting the clinical significance of tRFs in RSV infections (Supplementary Figure S1). This isoform of tRF5-GlyGCC-3 cannot be induced by SARS-CoV-2 and other viruses [9,43], suggesting virus-specific induction. In this study, no correlation between tRF5-CysGCA and RSV genome copies was found (Figure 1F), consistent with our previous finding showing that tRF5-CysGCA is not involved in regulating RSV replication. There was a trend showing the correlation between the expression of tRF5-GlyCCC and RSV genome (p = 0.08, Figure 1E), which may require more samples for confirmation. Nevertheless, the expression-viral genome correlation using limited sample numbers highlighted the possibility of using specific tRFs as biomarkers for viral infectious diseases.
In this study, we tried different methods to enrich tRNA-GluCTC as discussed in Drino et al. [44] and found the protocol summarized in Figure 2B reached a better result compared to the alternative Method 2 listed in Supplementary Figure S2. In brief, Method 2 used ion exchange chromatography followed by NaCl gradient elution as an essential step to enrich total tRNAs [44]. This step was found to be unnecessary and detrimental, however, as many tRNAs were wasted when 500 mM or greater NaCl was used to elute the larger RNAs (Supplementary Figure S2, right panels of B and C). In addition, the elongated processing may lead to potential RNA degradation. The method shown in Figure 2B and Method 1 in Supplementary Figure S2 were both found to enrich the desired tRNA; the method shown in Figure 2B was chosen because of its shorter exposure of RNAs to the annealing temperature of 60 °C.
As mentioned, tRNAs are enriched with chemical modifications. Such complexity, along with the existence of isoacceptor and isodecoder tRNAs, makes it hard to identify impacted tRNA modifications in response to environmental changes. Next-generation sequencing (NGS)-based RNA sequencing (RNA-Seq) has advanced the discovery of novel functional RNAs. However, it is limited in identifying modifications because the sequencing substrates are the reverse transcribed cDNA rather than the RNA itself. RNA-seq in combination with a pull-down using a modification-specific antibody is also limited in its ability to identify multiple tRNA modifications simultaneously. Nanopore-based RNA-Seq has been used recently for the sequencing of tRNAs [45], but it suffers from a high error rate and cannot reach a single-nucleotide resolution for RNA modifications. LC-MS/MS, as mentioned, cannot locate the modifications because it is restrained to the ribonucleoside level, losing information regarding the location and co-occurrence of modified nucleotides [46,47,48]. The MLC-Seq methods under continuing development by Shenglong Zhang’s lab, which combine MS analysis of RNA ladders obtained from randomly cut RNAs with complementary NSG mapping, are capable of sequencing RNA directly and simultaneously quantifying different RNA modifications with single-nucleotide stoichiometric precision [25,26,28,29]. In addition to identifying RSV-demethylated A57, this approach revealed a novel modification, i.e., a 1 Da reduction of D19 and D20; the associated modification mechanisms, regulatory enzymes, and the role of 1 Da reduction by RSV infections would need to be clarified in the future.
Overall, this study revealed a new mechanism controlling RSV replication. The importance of post-transcriptional modification of tRNAs in tRF biogenesis and the identification of the associated enzyme responsible for the modification provide critical insight that can lead to better strategies for antiviral developments to counter RSV infection.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v15010057/s1, Figure S1: RSV infection affects the expression of tRF5-GlyGCC-3.; Figure S2: Additional tRNA purification methods.; Supplementary Excel.

Author Contributions

Conceptualization, E.-J.C. and X.B.; methodology, E.-J.C., W.W., K.Z., X.Y., J.D., D.I., D.L.B. and K.S.T.; software, X.Y. and S.Z.; formal analysis, E.-J.C. and X.Y.; investigation, X.B., R.P.G. and S.Z.; writing—original draft preparation, E.-J.C.; writing—review and editing, X.B.; All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by grants from the US National Institutes of Health (NIH) R01 AI116812, R21 AI166543, AG069226, and TARRC Investigator-Initiated Research Award to X.B, the grant support from NIH R21 HG009576, R56 HG011099, and RM1 HG011563 to S.Z, and AI62885 to RPG.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and children who were enrolled in a study on the pathogenesis of RSV bronchiolitis under protocol #03-117 (approved by the Institutional Review Board of UTMB), for studies involving tRF measurement and RSV genome copies.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study and written informed consent has been obtained from the patient(s) to publish this paper” if applicable.

Data Availability Statement

The data presented in this study are openly available in Viruses.

Acknowledgments

We also thank Benjamin Johnson for editing the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Staebler, S.; Blake, S. Respiratory Syncytial Virus Disease: Immunoprophylaxis Policy Review and Public Health Concerns in Preterm and Young Infants. Policy Polit. Nurs. Pract. 2021, 22, 41–50. [Google Scholar] [CrossRef] [PubMed]
  2. Domachowske, J.B.; Anderson, E.J.; Goldstein, M. The Future of Respiratory Syncytial Virus Disease Prevention and Treatment. Infect. Dis. Ther. 2021, 10, 47–60. [Google Scholar] [CrossRef] [PubMed]
  3. Goldstein, M.; Harding, B.; Fayard, E. Guidance for palivizumab prophylaxis and implications for compliance. Pediatr. Pulmonol. 2021, 56, 3575–3576. [Google Scholar] [CrossRef] [PubMed]
  4. Baptista, B.; Riscado, M.; Queiroz, J.A.; Pichon, C.; Sousa, F. Non-coding RNAs: Emerging from the discovery to therapeutic applications. Biochem. Pharmacol. 2021, 189, 114469. [Google Scholar] [CrossRef] [PubMed]
  5. Hombach, S.; Kretz, M. Non-coding RNAs: Classification, Biology and Functioning. Adv. Exp. Med. Biol. 2016, 937, 3–17. [Google Scholar]
  6. Jackowiak, P.; Hojka-Osinska, A.; Philips, A.; Zmienko, A.; Budzko, L.; Maillard, P.; Budkowska, A.; Figlerowicz, M. Small RNA fragments derived from multiple RNA classes—The missing element of multi-omics characteristics of the hepatitis C virus cell culture model. BMC Genom. 2017, 18, 502. [Google Scholar] [CrossRef]
  7. Ruggero, K.; Guffanti, A.; Corradin, A.; Sharma, V.K.; De Bellis, G.; Corti, G.; Grassi, A.; Zanovello, P.; Bronte, V.; Ciminale, V.; et al. Small noncoding RNAs in cells transformed by human T-cell leukemia virus type 1: A role for a tRNA fragment as a primer for reverse transcriptase. J. Virol. 2014, 88, 3612–3622. [Google Scholar] [CrossRef] [Green Version]
  8. Wang, Q.; Lee, I.; Ren, J.; Ajay, S.S.; Lee, Y.S.; Bao, X. Identification and functional characterization of tRNA-derived RNA fragments (tRFs) in respiratory syncytial virus infection. Mol. Ther. 2013, 21, 368–379. [Google Scholar] [CrossRef] [Green Version]
  9. Wu, W.; Choi, E.J.; Wang, B.; Zhang, K.; Adam, A.; Huang, G.; Tunkle, L.; Huang, P.; Goru, R.; Imirowicz, I.; et al. Changes of Small Non-coding RNAs by Severe Acute Respiratory Syndrome Coronavirus 2 Infection. Front. Mol. Biosci. 2022, 9, 821137. [Google Scholar] [CrossRef]
  10. Eng, M.W.; Clemons, A.; Hill, C.; Engel, R.; Severson, D.W.; Behura, S.K. Multifaceted functional implications of an endogenously expressed tRNA fragment in the vector mosquito Aedes aegypti. PLoS Negl. Trop. Dis. 2018, 12, e0006186. [Google Scholar] [CrossRef] [Green Version]
  11. Zhou, J.; Liu, S.; Chen, Y.; Fu, Y.; Silver, A.J.; Hill, M.S.; Lee, I.; Lee, Y.S.; Bao, X. Identification of two novel functional tRNA-derived fragments induced in response to respiratory syncytial virus infection. J. Gen. Virol. 2017, 98, 1600–1610. [Google Scholar] [CrossRef] [PubMed]
  12. Choi, E.J.; Wu, W.; Zhang, K.; Lee, I.; Kim, I.-H.; Lee, Y.S.; Bao, X. ELAC2, an Enzyme for tRNA Maturation, Plays a Role in the Cleavage of a Mature tRNA to Produce a tRNA-Derived RNA Fragment During Respiratory Syncytial Virus Infection. Front. Mol. Biosci. 2020, 7, 609732. [Google Scholar] [CrossRef] [PubMed]
  13. Deng, J.; Ptashkin, R.N.; Chen, Y.; Cheng, Z.; Liu, G.; Phan, T.; Deng, X.; Zhou, J.; Lee, I.; Lee, Y.S.; et al. Respiratory Syncytial Virus Utilizes a tRNA Fragment to Suppress Antiviral Responses Through a Novel Targeting Mechanism. Mol. Ther. 2015, 23, 1622–1629. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Zhou, K.; Diebel, K.W.; Holy, J.; Skildum, A.; Odean, E.; Hicks, D.A.; Schotl, B.; Abrahante, J.E.; Spillman, M.A.; Bemis, L.T. A tRNA fragment, tRF5-Glu, regulates BCAR3 expression and proliferation in ovarian cancer cells. Oncotarget 2017, 8, 95377–95391. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Wu, W.; Lee, I.; Spratt, H.; Fang, X.; Bao, X. tRNA-Derived Fragments in Alzheimer’s Disease: Implications for New Disease Biomarkers and Neuropathological Mechanisms. J. Alzheimers Dis. 2021, 79, 793–806. [Google Scholar] [CrossRef]
  16. Ranjan, N.; Rodnina, M.V. tRNA wobble modifications and protein homeostasis. Translation 2016, 4, e1143076. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Jonkhout, N.; Tran, J.; Smith, M.A.; Schonrock, N.; Mattick, J.S.; Novoa, E.M. The RNA modification landscape in human disease. RNA 2017, 23, 1754–1769. [Google Scholar] [CrossRef] [Green Version]
  18. Tuorto, F.; Liebers, R.; Musch, T.; Schaefer, M.; Hofmann, S.; Kellner, S.; Frye, M.; Helm, M.; Stoecklin, G.; Lyko, F. RNA cytosine methylation by Dnmt2 and NSun2 promotes tRNA stability and protein synthesis. Nat. Struct. Mol. Biol. 2012, 19, 900–905. [Google Scholar] [CrossRef]
  19. Liu, F.; Clark, W.; Luo, G.; Wang, X.; Fu, Y.; Wei, J.; Wang, X.; Hao, Z.; Dai, Q.; Zheng, G.; et al. ALKBH1-Mediated tRNA Demethylation Regulates Translation. Cell 2016, 167, 816–828.e816. [Google Scholar] [CrossRef] [Green Version]
  20. Choi, E.J.; Ren, Y.; Chen, Y.; Chen, Y.; Liu, S.; Wu, W.; Ren, J.; Wang, P.; Garofalo, R.P.; Zhou, J.; et al. Exchange Proteins Directly Activated by cAMP and Their Roles in Respiratory Syncytial Virus Infection. J. Virol. 2018, 92, e01200-18. [Google Scholar] [CrossRef] [Green Version]
  21. Bao, X.; Indukuri, H.; Liu, T.; Liao, S.L.; Tian, B.; Brasier, A.R.; Garofalo, R.P.; Casola, A. IKKepsilon modulates RSV-induced NF-kappaB-dependent gene transcription. Virology 2010, 408, 224–231. [Google Scholar] [CrossRef] [PubMed]
  22. Ren, J.; Wang, Q.; Kolli, D.; Prusak, D.J.; Tseng, C.-T.K.; Chen, Z.J.; Li, K.; Wood, T.G.; Bao, X. Human metapneumovirus M2-2 protein inhibits innate cellular signaling by targeting MAVS. J. Virol. 2012, 86, 13049–13061. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Choi, E.J.; Wu, W.; Cong, X.; Zhang, K.; Luo, J.; Ye, S.; Wang, P.; Suresh, A.; Ullah, U.; Zhou, J.; et al. Broad Impact of Exchange Protein Directly Activated by cAMP 2 (EPAC2) on Respiratory Viral Infections. Viruses 2021, 13, 1179. [Google Scholar] [CrossRef] [PubMed]
  24. Bannister, R.; Rodrigues, D.; Murray, E.J.; Laxton, C.; Westby, M.; Bright, H. Use of a highly sensitive strand-specific quantitative PCR to identify abortive replication in the mouse model of respiratory syncytial virus disease. Virol. J. 2010, 7, 250. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Zhang, N.; Shi, S.; Yoo, B.; Yuan, X.; Li, W.; Zhang, S. 2D-HELS MS Seq: A General LC-MS-Based Method for Direct and de novo Sequencing of RNA Mixtures with Different Nucleotide Modifications. J. Vis. Exp. 2020, 161, e61281. [Google Scholar] [CrossRef] [PubMed]
  26. Zhang, N.; Shi, S.; Wang, X.; Ni, W.; Yuan, X.; Duan, J.; Jia, T.Z.; Yoo, B.; Ziegler, A.; Russo, J.; et al. Direct sequencing of tRNA by 2D-HELS-AA MS Seq reveals its different isoforms and dynamic base modifications. ACS Chem. Biol. 2020, 15, 1464–1472. [Google Scholar] [CrossRef] [PubMed]
  27. Zhang, N.; Shi, S.; Yuan, X.; Ni, W.; Wang, X.; Yoo, B.; Jia, T.Z.; Li, W.; Zhang, S. A General LC-MS-Based Method for Direct and De Novo Sequencing of RNA Mixtures Containing both Canonical and Modified Nucleotides. Methods Mol. Biol. 2021, 2298, 261–277. [Google Scholar] [PubMed]
  28. Zhang, N.; Shi, S.; Jia, T.Z.; Ziegler, A.; Yoo, B.; Yuan, X.; Li, W.; Zhang, S. A general LC-MS-based RNA sequencing method for direct analysis of multiple-base modifications in RNA mixtures. Nucleic Acids Res. 2019, 47, e125. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Zhang, S.; Yuan, X.; Su, Y.; Zhang, X.; Turkel, S.; Shi, S.; Wang, X.; Choi, E.-J.; Wu, W.; Liu, H.; et al. MLC-Seq: De novo Sequencing of Full-Length tRNAs and Quantitative Mapping of Multiple RNA Modifications b. Res. Sq. 2021. (preprint and under revision). [Google Scholar] [CrossRef]
  30. Liu, S.; Chen, Y.; Ren, Y.; Zhou, J.; Ren, J.; Lee, I.; Bao, X. A tRNA-derived RNA Fragment Plays an Important Role in the Mechanism of Arsenite -induced Cellular Responses. Sci. Rep. 2018, 8, 16838. [Google Scholar] [CrossRef] [Green Version]
  31. Deng, J.; Ptashkin, R.N.; Wang, Q.; Liu, G.; Zhang, G.; Lee, I.; Lee, Y.S.; Bao, X. Human metapneumovirus infection induces significant changes in small noncoding RNA expression in airway epithelial cells. Mol. Ther. Nucleic Acids 2014, 3, e163. [Google Scholar] [CrossRef] [PubMed]
  32. Okamoto, M.; Fujiwara, M.; Hori, M.; Okada, K.; Yazama, F.; Konishi, H.; Xiao, Y.; Qi, G.; Shimamoto, F.; Ota, T.; et al. tRNA modifying enzymes, NSUN2 and METTL1, determine sensitivity to 5-fluorouracil in HeLa cells. PLoS Genet. 2014, 10, e1004639. [Google Scholar] [CrossRef] [PubMed]
  33. Chen, Z.; Qi, M.; Shen, B.; Luo, G.; Wu, Y.; Li, J.; Lu, Z.; Zheng, Z.; Dai, Q.; Wang, H. Transfer RNA demethylase ALKBH3 promotes cancer progression via induction of tRNA-derived small RNAs. Nucleic Acids Res. 2019, 47, 2533–2545. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Rashad, S.; Niizuma, K.; Tominaga, T. tRNA cleavage: A new insight. Neural Regen. Res. 2020, 15, 47–52. [Google Scholar] [CrossRef] [PubMed]
  35. Oerum, S.; Degut, C.; Barraud, P.; Tisne, C. m1A Post-Transcriptional Modification in tRNAs. Biomolecules 2017, 7, 20. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Zhang, C.; Jia, G. Reversible RNA Modification N(1)-methyladenosine (m(1)A) in mRNA and tRNA. Genom. Proteom. Bioinform. 2018, 16, 155–161. [Google Scholar] [CrossRef]
  37. Prehn, J.H.M.; Jirstrom, E. Angiogenin and tRNA fragments in Parkinson’s disease and neurodegeneration. Acta Pharmacol. Sin. 2020, 41, 442–446. [Google Scholar] [CrossRef]
  38. Olvedy, M.; Scaravilli, M.; Hoogstrate, Y.; Visakorpi, T.; Jenster, G.; Martens-Uzunova, E. A comprehensive repertoire of tRNA-derived fragments in prostate cancer. Oncotarget 2016, 7, 24766. [Google Scholar] [CrossRef] [Green Version]
  39. Goodarzi, H.; Liu, X.; Nguyen, H.C.; Zhang, S.; Fish, L.; Tavazoie, S.F. Endogenous tRNA-Derived Fragments Suppress Breast Cancer Progression via YBX1 Displacement. Cell 2015, 161, 790–802. [Google Scholar] [CrossRef] [Green Version]
  40. Choi, E.J.; Ren, J.; Zhang, K.; Zhang, K.; Wu, W.; Lee, Y.S.; Lee, I.; Bao, X. The Importance of AGO 1 and 4 in Post-Transcriptional Gene Regulatory Function of tRF5-GluCTC, an Respiratory Syncytial Virus-Induced tRNA-Derived RNA Fragment. Int. J. Mol. Sci. 2020, 21, 8766. [Google Scholar] [CrossRef] [PubMed]
  41. Wu, W.; Choi, E.J.; Lee, I.; Lee, Y.S.; Bao, X. Non-Coding RNAs and Their Role in Respiratory Syncytial Virus (RSV) and Human Metapneumovirus (hMPV) Infections. Viruses 2020, 12, 345. [Google Scholar] [CrossRef] [PubMed]
  42. Magee, R.; Londin, E.; Rigoutsos, I. TRNA-derived fragments as sex-dependent circulating candidate biomarkers for Parkinson’s disease. Park. Relat. Disord. 2019, 65, 203–209. [Google Scholar] [CrossRef] [PubMed]
  43. Selitsky, S.R.; Baran-Gale, J.; Honda, M.; Yamane, D.; Masaki, T.; Fannin, E.E.; Guerra, B.; Shirasaki, T.; Shimakami, T.; Kaneko, S.; et al. Small tRNA-derived RNAs are increased and more abundant than microRNAs in chronic hepatitis B and C. Sci. Rep. 2015, 5, 7675. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Drino, A.; Oberbauer, V.; Troger, C.; Janisiw, E.; Anrather, D.; Hartl, M.; Kaiser, S.; Kllner, S.; Schafer, M.R. Production and purification of endogenously modified tRNA-derived small RNAs. RNA Biol. 2020, 17, 1104–1115. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Thomas, N.; Poodari, V.; Jain, M.; Olsen, H.; Akeson, M.; Abu-Shumays, R. Direct Nanopore Sequencing of Individual Full Length tRNA Strands. ACS Nano 2021, 15, 16642–16653. [Google Scholar] [CrossRef]
  46. Su, D.; Chan, C.T.; Gu, C. Quantitative analysis of ribonucleoside modifications in tRNA by HPLC-coupled mass spectrometry. Nat. Protoc. 2014, 9, 828–841. [Google Scholar] [CrossRef]
  47. Lauman, R.; Garcia, B.A. Unraveling the RNA modification code with mass spectrometry. Mol. Omics. 2020, 16, 305–315. [Google Scholar] [CrossRef]
  48. Wetzel, C.; Limbach, P.A. Mass spectrometry of modified RNAs: Recent developments. Analyst 2016, 141, 16–23. [Google Scholar] [CrossRef]
Figure 1. Changes in the expression of tRF5s in NPS samples by RSV. qRT-PCR was performed to detect (A) tRF5-GluCTC, (B) tRF5-GlyCCC, and (C) tRF5-CysGCA in the NPS from RSV and control (CN) patients. An unpaired two-tailed Mann-Whitney U test was used to compare two independent groups. The three asterisks represent a p-value of 0.001. For correlation analyses for tRF5-GluCTC (D), tRF5-GlyCCC (E), or tRF5-CysGCA (F) with NPS viral genome copies, we performed Spearman’s rank correlation test. Spearman’s rank correlation coefficient (RS) was used to determine correlations. A p-value of less than 0.05 was considered significant.
Figure 1. Changes in the expression of tRF5s in NPS samples by RSV. qRT-PCR was performed to detect (A) tRF5-GluCTC, (B) tRF5-GlyCCC, and (C) tRF5-CysGCA in the NPS from RSV and control (CN) patients. An unpaired two-tailed Mann-Whitney U test was used to compare two independent groups. The three asterisks represent a p-value of 0.001. For correlation analyses for tRF5-GluCTC (D), tRF5-GlyCCC (E), or tRF5-CysGCA (F) with NPS viral genome copies, we performed Spearman’s rank correlation test. Spearman’s rank correlation coefficient (RS) was used to determine correlations. A p-value of less than 0.05 was considered significant.
Viruses 15 00057 g001
Figure 2. The schematic summaries on tRNA-GluCTC preparation. (A) Sequencing information on tRF5-GluCTC and its parent tRNA, the antisense probe used to detect tRNA by Northern blot, and biotinylated antisense oligo to pull down tRNA-GluCTC. (B) The workflow to purify tRNA-GluCTC. (C) The input and purified tRNA-GluCTC from mock- and RSV-infected cells. Left panel: SYBRTM Green II RNA staining of total RNAs for mock- and RSV-infected cells (left two lanes) and enriched tRNAs (right two lanes). Middle panel: input RNAs and enriched tRNAs were detected by antisense oligo against tRNA-GluCTC. Right panel: input RNAs and enriched tRNAs were detected by antisense against tRNA-GlyGCC.
Figure 2. The schematic summaries on tRNA-GluCTC preparation. (A) Sequencing information on tRF5-GluCTC and its parent tRNA, the antisense probe used to detect tRNA by Northern blot, and biotinylated antisense oligo to pull down tRNA-GluCTC. (B) The workflow to purify tRNA-GluCTC. (C) The input and purified tRNA-GluCTC from mock- and RSV-infected cells. Left panel: SYBRTM Green II RNA staining of total RNAs for mock- and RSV-infected cells (left two lanes) and enriched tRNAs (right two lanes). Middle panel: input RNAs and enriched tRNAs were detected by antisense oligo against tRNA-GluCTC. Right panel: input RNAs and enriched tRNAs were detected by antisense against tRNA-GlyGCC.
Viruses 15 00057 g002
Figure 3. Identified tRNA modification changes by RSV. (A) Sequencing workflow. More detailed information can be found in Materials and Methods. (B,C) Stoichiometric changes of m1A at position 57 methylated nucleotides by RSV infection. The percentage of m1A at position 57 was reduced from 51% in the mock-infected cells (B) to 11% in RSV-infected cells (C). (D) The putative change of 1 Da on the pseudouritidine D by RSV.
Figure 3. Identified tRNA modification changes by RSV. (A) Sequencing workflow. More detailed information can be found in Materials and Methods. (B,C) Stoichiometric changes of m1A at position 57 methylated nucleotides by RSV infection. The percentage of m1A at position 57 was reduced from 51% in the mock-infected cells (B) to 11% in RSV-infected cells (C). (D) The putative change of 1 Da on the pseudouritidine D by RSV.
Viruses 15 00057 g003
Figure 4. The impact of ALKBH1 on the RSV infection. (A) ALKBH1 gene knockdown prevents the induction of tRF5-GluCTC by RSV. The input of total RNAs in mock- and RSV-infected cells were checked by EB staining and induced tRF5-GluCTC was determined by Northern blot. (B,C) ALKBH1 gene knockdown decreased the expression of RSV proteins. The ALKBH1 gene knockdown was confirmed by Western blot using an anti-ALKBH1 antibody. The membrane was stripped and reprobed for RSV proteins and β-actin (B). Densitometric analysis of the band intensity of ALKBH1 was performed using the histogram function of Adobe Photoshop. Results are shown after corrections to β-actin. Results are representative of three independent experiments (C). (D) The impact of ALKBH1 on RSV virus titration was also determined. The double asterisk represents p values of 0.01. Data are shown as means ± SE. (E) The nuclear and cytosolic expression of ALKBH1 along RSV infection. The membranes were stripped and reprobed for lamin B as a loading control for nuclear fraction, while β-actin was used as a loading control for cytosol fraction.
Figure 4. The impact of ALKBH1 on the RSV infection. (A) ALKBH1 gene knockdown prevents the induction of tRF5-GluCTC by RSV. The input of total RNAs in mock- and RSV-infected cells were checked by EB staining and induced tRF5-GluCTC was determined by Northern blot. (B,C) ALKBH1 gene knockdown decreased the expression of RSV proteins. The ALKBH1 gene knockdown was confirmed by Western blot using an anti-ALKBH1 antibody. The membrane was stripped and reprobed for RSV proteins and β-actin (B). Densitometric analysis of the band intensity of ALKBH1 was performed using the histogram function of Adobe Photoshop. Results are shown after corrections to β-actin. Results are representative of three independent experiments (C). (D) The impact of ALKBH1 on RSV virus titration was also determined. The double asterisk represents p values of 0.01. Data are shown as means ± SE. (E) The nuclear and cytosolic expression of ALKBH1 along RSV infection. The membranes were stripped and reprobed for lamin B as a loading control for nuclear fraction, while β-actin was used as a loading control for cytosol fraction.
Viruses 15 00057 g004
Table 1. Patient gender, age, and race information.
Table 1. Patient gender, age, and race information.
Control (CN)RSV Patient
No. of patients89
Gender (Male:Female)5:35:4
Mean age (months (range))7.68 (0.75~14)7.88 (1.25~16)
Race (African American:Caucasian:Asian)2:5:13:6:0
Table 2. Primer information to detect tRF5s using qRT-PCR.
Table 2. Primer information to detect tRF5s using qRT-PCR.
tRF5–GlyCCCSequenceGCAUUGGUGGUUCAGUGGUAGAAUUCUCGCC
Forward primerGCATGGGTGGTTCAGTG
Reverse primerCGTCGGACTGTAGAACTCTCAAAGC
tRF5–GluCTCSequenceUCCCUGGUGGUCUAGUGGUUAGGAUUCGGCGCU
Forward primerTCCCTGGTGGTCTAGTG
Reverse primerCGTCGGACTGTAGAACTCTCAAAGC
tRF5–CysGCASequenceGGGUAUAGCUCAGUGGUAGAGCAUUUGACUGC
Forward primerAGTGGTAGAGCATTTGACTGC
Reverse primerCGTCGGACTGTAGAACTCTCAAAGC
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Choi, E.-J.; Wu, W.; Zhang, K.; Yuan, X.; Deng, J.; Ismail, D.; Buck, D.L.; Thomason, K.S.; Garofalo, R.P.; Zhang, S.; et al. Parent tRNA Modification Status Determines the Induction of Functional tRNA-Derived RNA by Respiratory Syncytial Virus Infection. Viruses 2023, 15, 57. https://doi.org/10.3390/v15010057

AMA Style

Choi E-J, Wu W, Zhang K, Yuan X, Deng J, Ismail D, Buck DL, Thomason KS, Garofalo RP, Zhang S, et al. Parent tRNA Modification Status Determines the Induction of Functional tRNA-Derived RNA by Respiratory Syncytial Virus Infection. Viruses. 2023; 15(1):57. https://doi.org/10.3390/v15010057

Chicago/Turabian Style

Choi, Eun-Jin, Wenzhe Wu, Ke Zhang, Xiaohong Yuan, Junfang Deng, Deena Ismail, Darby L. Buck, Kerrie S. Thomason, Roberto P. Garofalo, Shenglong Zhang, and et al. 2023. "Parent tRNA Modification Status Determines the Induction of Functional tRNA-Derived RNA by Respiratory Syncytial Virus Infection" Viruses 15, no. 1: 57. https://doi.org/10.3390/v15010057

APA Style

Choi, E. -J., Wu, W., Zhang, K., Yuan, X., Deng, J., Ismail, D., Buck, D. L., Thomason, K. S., Garofalo, R. P., Zhang, S., & Bao, X. (2023). Parent tRNA Modification Status Determines the Induction of Functional tRNA-Derived RNA by Respiratory Syncytial Virus Infection. Viruses, 15(1), 57. https://doi.org/10.3390/v15010057

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop