Abstract
Respiratory syncytial virus (RSV) or measles virus (MeV) infection modifies host responses through small non-coding RNA (sncRNA) expression. We show that RSV or MeV infection of neuronal cells induces sncRNAs including various microRNAs and transfer RNA fragments (tRFs). We show that these tRFs originate from select tRNAs (GCC and CAC for glycine, CTT and AAC for Valine, and CCC and TTT for Lysine). Some of the tRNAs are rarely used by RSV or MeV as indicated by relative synonymous codon usage indices suggesting selective cleavage of the tRNAs occurs in infected neuronal cells. The data implies that differentially expressed sncRNAs may regulate host gene expression via multiple mechanisms in neuronal cells.
Introduction
Small non-coding RNAs (sncRNAs) are <200 nucleotides and include microRNAs (miRs), PIWI-interacting RNAs (piRs), and transfer RNA fragments (tRNA)-derived RNA fragments (tRFs) (Petrov et al., 2017). The expression of sncRNAs differs in their biogenesis and functional activity (Bartel, 2009; Bakre et al., 2019). Most sncRNAs induced by viral infection involve regulatory mechanisms affecting virus replication, persistence, or latency (Grassmann and Jeang, 2008; Roizman et al., 2011), and immune evasion (Cullen, 2013; Bernier and Sagan, 2018). RSV typically infects ciliated respiratory epithelial cells, while MeV may infect dendritic cells and subsequently infect the epithelia for virus release from the respiratory tract (Lay et al., 2016; Lin and Richardson, 2016; Singh et al., 2016; Taniguchi et al., 2019). These viruses cause a spectrum of diseases and may infect a spectrum of people, i.e., the very young to the elderly (Griffin, 2018; Mazur et al., 2018). RSV is a major cause of serious lower respiratory tract infections in childhood (Mirra et al., 2018; Paes, 2018), but also affects the elderly and immune-compromised (Haber, 2018). MeV causes acute viral infection with respiratory involvement leading to symptoms of rash and more serious complications, and in some cases to mortality (Bohmwald et al., 2018; Leung et al., 2018). Both RSV and MeV are paramyxoviruses with non-segmented negative-sense single-stranded RNA genomes (Russell et al., 2018). Unfortunately, there are no approved RSV vaccines and anti-viral treatments are limited. Although there is an effective MeV vaccine, incidences of MeV infection are increasing worldwide due to vaccine complacency.
Respiratory syncytial virus and MeV can infect neuronal cells, and MeV may cause persistent infection (Rota et al., 2017; Griffin et al., 2018). Compromised immunity may facilitate neuronal infection (Rima and Duprex, 2005; Omar et al., 2017; Ganesan et al., 2018; Ferren et al., 2019). MeV infection of neuronal cells may lead to subacute sclerosing panencephalitis (SSPE), and measles inclusion body encephalitis (MIBE) in immunocompromised individuals. In MIBE, MeV replicates in the central nervous system usually because of an inadequate immune response (McQuaid et al., 1998; Hardie et al., 2013; Rota et al., 2017). RSV infection in mice has been associated with exaggerated neurogenic inflammation in the airways, and studies have shown that the neuropeptide, substance P (SP) is a mediator of neurogenic inflammation (King et al., 2001; Tan et al., 2008).
Expression of the RSV G protein is associated with increased pulmonary expression of SP, and lung neurons express CX3CR1 (Tripp et al., 2000, 2002, 2003; Haynes et al., 2003). CX3CR1 is the receptor for fractalkine (aka neurotactin) and CX3CR1 has been shown to bind RSV G protein (Tripp et al., 2001; Harcourt et al., 2004; Choi et al., 2012; Bakre et al., 2017; Bergeron et al., 2021).
Studies in mice suggest that RSV may cause lasting infection (Dakhama et al., 1997; Li et al., 2006). RSV has been shown to infect mouse primary cortical neuronal cells as shown by co-localization of the RSV N protein and neuronal markers (Li et al., 2006). These findings support how RSV may cause neurological irregularities in patients (Eisenhut, 2006; Li et al., 2006; Morichi et al., 2017). Several studies have described an association between lower respiratory tract RSV infection in infancy and the subsequent development of persistent wheezing in children (Zhou N. et al., 2017; Zhong et al., 2018). In mice, RSV induces long-term airway disease characterized by chronic airway inflammation and airway hyperreactivity (Long et al., 2016; Zhou N. et al., 2017), and RSV antigens have been detected in the lungs >100 days (Schwarze et al., 2004). However, RSV persistence does not appear in the bronchial epithelium (the primary site of viral replication) but in deeper lung structures (Long et al., 2016; Zhou N. et al., 2017).
Respiratory syncytial virus and MeV infection modulate sncRNA expression inducing mRNA degradation or translation inhibition and have a role in determining the level of protein expression of host cells (Munday et al., 2012). It has been shown that tRNAs are cleaved during cellular stress, and that cleavage results in tRFs that contain the 5′ end (5′tRFs) or the 3′ end (3′tRFs) of the parent tRNA molecule (Sofos et al., 2015; Advani and Ivanov, 2019). Although the specific activity of tRFs is not well-understood, tRFs can behave as small interfering RNA leading to the degradation of transcripts (Ivanov et al., 2011) which can regulate ribosomal loading and protein chain elongation (Sobala and Hutvagner, 2013).
Accumulating data have shown that both coding and non-coding transcriptomes are modified during RSV infection. We have shown that RSV infection of human A549 lung cells modifies sncRNA expression and that the RSV G protein and NS1/NS2 proteins modulate miR expression (Bakre et al., 2012, 2015, 2017). Other studies have shown that the pattern of miR expression is modified following RSV infection of nasal epithelial cells (Inchley et al., 2015; Hasegawa et al., 2018), and PBMCs from children (Wang et al., 2017). Deregulated miR expression occurs in MeV-infected human neuroblastoma cells and PBMCs (Inchley et al., 2015; Yis et al., 2015; Naaman et al., 2017; Hasegawa et al., 2018), as well as in PBMCs from RSV-infected children (Wang et al., 2017). It is known that MeV persistence in neuroblastoma cells is assisted by the downregulation of CDK6, a component of cell cycle progression regulated by miRNA-124 (Riddell et al., 2007; Naaman et al., 2017; Griffin et al., 2018).
In this study, we examined sncRNA responses in neuronal SHSY5Y (SHS) cells infected with RSV or MeV and show that infection with either of these viruses modifies miRs and tRF expression in a temporal and virus-specific manner. Anti-microbial peptides (AMPs) are generated in the response to pathogens, and arginine and glycine are especially abundant in AMPs (Mishra and Wang, 2012). In this study, we show that RSV or MeV infection modifies the extent of transfer RNA fragments (tRFs), and the finding of decreased full-length tRNAs after infection indicates decreased translation rates possibly indicative of a lessened anti-viral state. Together, these findings strengthen the notion that sncRNAs have a key role in regulating the host cell response to infection and viral replication.
Materials and Methods
Cells and Viruses
Human SHSY5Y (SHS) neuroblastoma cells (ATCC CRL 2266) were grown in DMEM (Sigma, St. Louis, MO, United States) containing 10% FBS (Atlanta Biologicals, Atlanta, GA, United States). SHS cells were maintained at 37°C, 5% CO2. RSV strain A2 (ATCC CCL 81) was propagated in Vero cells (Oshansky et al., 2009; Dahiya and Atreya, 2014). Briefly, 80% confluent Vero cells were infected (MOI = 1.0) for 1 h at 37°C, 5% CO2. After infection, 2% FBS was added to the DMEM, and the cells were incubated for 3 days until syncytia were evident. On the day of harvest, the cells were removed and centrifuged at 500 × g for 15 min at 4°C. RSV strain A2 was used for the infection in SHS cells and rgRSV224 GFP (Hallak et al., 2000, 2007; Techaarpornkul et al., 2002) (a kind gift from Dr. Michael Teng, University of South Florida) for microscopy studies. The WT strain Dublin-3267 of MeV was grown in Vero cells expressing human CD150 (SLAM), a high-affinity receptor for MeV (Tatsuo et al., 2000).
RNA Isolation and Library Construction
Total RNA from mock-treated or virus-infected SHS cells was isolated using RNAzol RT (MRC, Cincinnati, OH, United States) per the manufacturer’s instructions and divided into small and large RNA fractions and quantified using a Qubit fluorimeter-based RNA assay (Invitrogen Life Tech, Carlsbad, CA, United States). Illumina next-generation sequencing (NGS) libraries were constructed using a TruSeq® small RNA library (Illumina, San Diego, CA, United States). Briefly, small RNA from mock-treated SHS cells (n = 2/group), or RSV- or MeV-infected SHS cells (n = 3/group) was isolated at 12 or 24 hpi and ligated with 3′ and 5′ RNA adaptors using T4 RNA ligase 2 (Epicenter Biotechnologies, Madison, WI, United States) at 28°C for 60 min. Ligated RNA was reverse transcribed using proprietary reverse primers and the first-strand cDNA was amplified with a combination of a proprietary forward primer and reverse primers containing unique index bar codes per the manufacturer’s instructions. Amplicons were run on a Tapestation 2200 (Agilent, Santa Clara, CA, United States) high-sensitivity DNA chip to visualize products. Indexed samples were pooled per the manufacturer’s protocol and electrophoresed on a 6% TBE-urea PAGE gel. Bands corresponding to 147 and 157 bp corresponding to adaptor-ligated mature miRNA and other non-coding RNAs were cut out and cDNAs were precipitated. Libraries were validated on a Tapestation 2200 (Agilent) using a high-sensitivity DNA chip. Libraries were denatured for 5 min using 5 M NaOH and loaded onto a MiSeq (Agilent) for sequencing by synthesis. Reads were trimmed to remove adaptor sequences and then analyzed for miRNA differential expression using Chimira (Vitsios and Enright, 2015) and OASIS 2.0 (Rahman et al., 2018a). Read counts were normalized to a total number of reads per sample and then differential expression was calculated using DESeq2 (Love et al., 2014). All reads used for analysis had quality scores (Q_ > 30. This represents less than or equal to 1 error per 1,000 nts corresponding to a 99.99% accuracy.
miR Analysis
The host genes and pathways predicted to be regulated by miRs were determined by DIANA miRPath analysis (Vlachos et al., 2015). miRs were examined using the DIANA-miRPath server and targets extracted from the Tarbase database (Karagkouni et al., 2018). Enriched KEGG pathway analysis was used to determine intersections (Kanehisa et al., 2016, 2017). We applied a p-value threshold of ≤0.05 along with a false-discovery rate correction.
Virus Quantitation
QRT-PCR of RSV M, or MeV M, and N gene primer-probes were examined for each replicate to confirm infection of SHS cells. Sequences of the primers/probes used were RSV M forward 5′- ACACCCAAGGGACCTTCACTAAGA, RSV M reverse 5′-GCCTTGATTTCACAGGGTGTGGTT-3′, RSV M probe/56-FAM/AAG TGC AGT GCT AGC ACA AAT GCC CA/3BHQ_1. For MeV, primer sequences were forward primer MeVEdNPF-5′-CAAGACCCTGAGGGATTCAA-3′ and reverse primer MeVEdNPR-5′-CTCTCGCATCACTTGCTCTG-3. For both viruses, 100 ng of RNA was reverse transcribed using Lunascript RT Supermix (New England Biolabs, Ipswitch, MA, United States) using the following conditions: primer annealing at 25°C for 2 min, reverse transcription at 55°C for 30 min followed by heat-inactivation of reverse transcriptase at 95°C for 1 min. Undiluted cDNA was used for PCR. PCR reactions contained cDNA, 0.4 μM of forward and reverse oligomers, 0.8 μM of the probe, 2X Luna Universal probe master mix (New England Biolabs) and were amplified as follows: initial denaturation at 95°C for 60 s, 40 cycles of denaturation at 95°C for 10 s, and combined annealing and extension at 60°C for 30 s. For MeV NP amplification, PCR reactions contained cDNA, 10 μM of forward and reverse primers, Q5 2X High-Fidelity master mix (New England Biolabs), and was amplified using the conditions: initial denaturation at 98°C for 30 s, 25 cycles of 98°C for 10 s, annealing at 65°C for 30 s and extension at 72°C for 2 min, and a final extension at 72°C for 2 min. MeV NP amplicons were resolved on a 1% agarose gel in low osmolality buffer at 5 V/cm and imaged using a Protein Simple imager (Invitrogen).
Results
sncRNAs and Virus-Infected SHS Cells
Previous studies from our laboratory have shown that RSV infection of A549 cells modifies miR expression (Bakre et al., 2012, 2015, 2017) however, information regarding miR expression in RSV- or MeV-infected neuronal SHS cells is inadequate (Naaman et al., 2017; Cakmak Genc et al., 2018; Haralambieva et al., 2018). Neuronal cells are not the primary cell type infected by RSV, but they can be infected and may act as a reservoir of infection. Mechanisms that contribute to this are poorly understood. Thus, we examined SHS cells that were RSV- or MeV-infected (MOI = 1.0), or mock-treated and sncRNA expression determined at 12 and 24 hpi by Illumina NGS. We confirmed RSV infection of SHS neuronal cells using RSV-GFP. SHS cells were infected with RSV-GFP at a MOI = 0.1 (Supplementary Figure 1). At 48 h and 72 hpi, virus replication was evident using combined dark field and fluorescence microscopy (UV Eclipse Ti2000-U microscope, Nikon).
It is important to note that production of infectious viral particles is restricted in the neuronal environment and viral spread is mediated by fusogenic events leading to syncytia formation that are common to both MeV and RSV infection (Sato et al., 2018). Total reads expressed in a text-based format (fastq) were used and the corresponding reads passing quality scores (Q > 30) are listed in Table 1. Raw fastq reads were mapped to the human genome, and miRs using Chimira (Vitsios and Enright, 2015) and Oasis 2.0 (Rahman et al., 2018a) were normalized to the total read counts and used to calculate differential expression using DESeq 2.0 (Love et al., 2014) relative to mock-treated cells. We used principal component analysis (PCA) to determine sample clustering based on similarity for miR expression profiling. PCA is a method of taking combinations of measured variables and lowering the dimensionality of the data to find the principal source of variance. As the human genome encodes > 2,300 miRs and we examined RSV and MeV and two time points, the dataset had > 9,200 variables. PCA identified two major component dimensions—PC1 consists of the time component in addition to virus and miRNA expression, and PC2 is the infection state. PC1 and PC2 are principal components and organized by the percent of variability they explain. Being derived from combinations of variables, PCs do not have labels/units. Samples 1–3 (12 h MeV), 4–6 (12 h RSV), 7–9 (24 h MeV), and 10–12 (24 h RSV) clustered together and away from mock infected cells (samples 13–14) suggesting that the sncRNA response to MeV or RSV infection is distinct from mock infected cells. In addition, clustering of 12 and 24 h MeV and RSV datasets suggests that infections elicit temporal responses. The PCA results indicate that the time post infection is a major component driving differential miRNA expression as shown on PC1 (x-axis with 54% variance) compared to pathogen (PC2 y-axis 13% variance). This is evident in the clustering among mock-treated vs. virus-infected cells (Figure 1). Similarly, miR expression profiles in RSV- or MeV-infected SHS cells at 12 and 24 hpi revealed distinct clustering. Most RSV samples clustered separately from MeV at each time-point (Figure 1).
TABLE 1
| Sample | Original # of reads | % Trimmed reads | % Uniquely mapped reads |
| 12 hpi–RSV-1 | 5,953 | 98% | 5% |
| 12 hpi–RSV-2 | 18,599 | 92% | 42% |
| 12 hpi–RSV-3 | 7,546 | 94% | 23% |
| 12 hpi–MeV-1 | 101,182 | 78% | 50% |
| 12 hpi–MeV-2 | 27,880 | 74% | 41% |
| 12 hpi–MeV-3 | 12,088 | 41% | 11% |
| 24 hpi–RSV-1 | 46,761 | 72% | 10% |
| 24 hpi–RSV-2 | 239,242 | 74% | 9% |
| 24 hpi–RSV-3 | 74,687 | 74% | 12% |
| 24 hpi–MeV-1 | 27,320 | 66% | 7% |
| 24 hpi–MeV-2 | 26,797 | 62% | 6% |
| 24 hpi–MeV-3 | 92,746 | 66% | 9% |
| SHS-1 | 57,494 | 79% | 51% |
| SHS-2 | 77,356 | 79% | 44% |
Summary of read numbers.
Numbers indicate raw read counts, percent of reads with adaptor trimming [59], and those mapping to unique loci on the genome.
FIGURE 1
The data were filtered to contain only read numbers > 5 (Table 1), and those having p-values ≤0.05, and fold-change ≤1.5 or ≥1.5 were used to help identify downregulated vs. upregulated miRs. RSV or MeV infection was confirmed by qPCR. The chief deregulated miRs expressed in SHS cells following RSV- or MeV-infection are shown in Table 2. Compared to mock-infection, expression of miR-10b-5p, miR-181c-5p, and miR-125b-5p were strikingly downregulated in SHS cells at 12 or 24 hpi following RSV or MeV infection. In contrast, expression of p-miR-330, p-miR-338, p-miR-317, p-miR-106, p-miR-113, p-miR-256, p-miR-232, p-miR-103 and miR-182-5p, miR-4532, miR-1273g-3p, miR-483-3p, miR-3141, miR-1246 miR-3196, miR-4516, and miR-4488 were substantially upregulated in SHS cells.
TABLE 2
| Mock | 12 h | 24 h | |||||
| Mature | log2 fold-change | padj | Mock | RSV | MeV | RSV | MeV |
| miR-10b-5p | −2 | 0.00 | 1,416, 1,636 | 94, 144, 78 | 369, 262, 128 | 165, 248, 158 | 166, 173, 243 |
| miR-181c-5p | −2 | 0.03 | 25, 26 | 12, 3, 6 | 6, 10, 12 | 9, 3, 9 | 10, 5, 6 |
| miR-125b-5p | −2 | 0.02 | 50, 34 | 12, 12, 10 | 9, 9, 27 | 12, 10, 10 | 7, 21, 12 |
| miR-375 | 2 | 0.01 | 56, 49 | 82, 225, 273 | 137, 165, 155 | 168, 132, 193 | 128, 163, 136 |
| p-miR-247 | 2 | 0.02 | 128, 147 | 222, 336, 458 | 131, 298, 530 | 560, 688, 694 | 416, 446, 387 |
| miR-21-5p | 2 | 0.03 | 11, 11 | 211, 19, 29 | 43, 23, 15 | 32, 65, 27 | 57, 26, 62 |
| miR-877-5p | 2 | 0.02 | 1, 3 | 12, 28, 45 | 5, 10, 19 | 5, 6, 8 | 14, 21, 4 |
| miR-7704 | 3 | 0.03 | 1, 1 | 0, 1, 6 | 1, 2, 4 | 25, 16, 17 | 14, 10, 18 |
| miR-100-5p | 3 | 0.05 | 0, 1 | 35, 1, 3 | 3, 2, 4 | 2, 1, 4 | 10, 10, 3 |
| p-miR-330 | 3 | 0.01 | 1, 2 | 23, 6, 10 | 3, 4, 15 | 17, 13, 6 | 17, 52, 16 |
| miR-182-5p | 3 | 0.00 | 3, 3 | 12, 22, 16 | 34, 30, 23 | 18, 37, 24 | 20, 47, 39 |
| miR-4532 | 4 | 0.02 | 0, 0 | 0, 0, 0 | 0, 0, 0 | 11, 29, 10 | 7, 26, 10 |
| miR-1273g-3p | 5 | 0.00 | 0, 0 | 0, 1, 3 | 0, 1, 0 | 15, 20, 18 | 24, 21, 22 |
| miR-483-3p | 4 | 0.00 | 0, 0 | 12, 11, 6 | 3, 10, 12 | 3, 4, 9 | 14, 0, 3 |
| miR-3141 | 5 | 0.05 | 0, 0 | 0, 0, 0 | 0, 0, 0 | 29, 32, 14 | 61, 47, 17 |
| p-miR-338 | 5 | 0.01 | 0, 0 | 0, 0, 0 | 0, 0, 0 | 23, 32, 10 | 24, 47, 11 |
| miR-1246 | 5 | 0.00 | 0, 1 | 0, 3, 3 | 0, 1, 0 | 32, 42, 22 | 44, 110, 32 |
| p-miR-317 | 6 | 0.00 | 0, 0 | 0, 0, 0 | 0, 0, 0 | 42, 44, 36 | 47, 68, 41 |
| miR-3196 | 6 | 0.02 | 0, 0 | 0, 0, 0 | 0, 0, 0 | 29, 52, 14 | 54, 68, 24 |
| miR-4516 | 6 | 0.01 | 0, 0 | 0, 0, 3 | 0, 0, 0 | 20, 20, 6 | 41, 26, 11 |
| p-miR-106 | 6 | 0.01 | 0, 0 | 0, 0, 0 | 0, 0, 0 | 22, 18, 12 | 27, 52, 7 |
| miR-4488 | 6 | 0.00 | 0, 0 | 0, 0, 0 | 0, 0, 0 | 25, 22, 11 | 34, 63, 15 |
| p-miR-113 | 6 | 0.00 | 0, 0 | 0, 0, 0 | 0, 0, 0 | 37, 61, 35 | 44, 68, 50 |
| p-miR-256 | 6 | 0.00 | 0, 0 | 12, 0, 3 | 0, 1, 0 | 12, 11, 10 | 10, 10, 3 |
| p-miR-232 | 6 | 0.00 | 0, 0 | 0, 0, 0 | 0, 0, 0 | 8, 16, 6 | 20, 26, 9 |
| p-miR-103 | 7 | 0.00 | 0, 0 | 0, 1, 0 | 0, 0, 0 | 22, 29, 21 | 17, 5, 34 |
Normalized read counts between mock, RSV- and MeV-infected samples as analyzed by OASIS 2.0 (Rahman et al., 2018a).
Negative log2 fold-change values indicate downregulation.
“p,” predicted miR; RSV, respiratory syncytial virus; MeV, measles virus; padj, Adjusted p-value.
Studies suggest that modifications involving Adenosine to Inosine (A- > I) or Cytosine to Uracil (C- > U) within miRs, N6 methylation, and oxidation affect miR stability, maturation and target specificity (Alarcon et al., 2015; Wang et al., 2015; Correia De Sousa et al., 2019). Post-transcriptional modifications to miRNAs have been shown to alter their binding and function (Pirouz et al., 2019; Briand et al., 2020; Cheray et al., 2020). We therefore examined epitranscriptional modifications of the deregulated miRs using Chimira (Figure 2; Vitsios and Enright, 2015). Most modifications were located on the 3p arms of the deregulated miRs. RSV infection in SHS cells decreased the abundance of 3p modifications at 12 hpi (Figure 2 and Supplementary Table 1).
FIGURE 2
RSV or MeV Infection Induces tRNA Fragments
Cellular stress triggers the degradation of tRNA to tRFs (Shen et al., 2018). tRFs may encompass the 5′, 3′, or internal anti-codon loop and are classified as 5′-, 3-, or class I tRFs (Shen et al., 2018). Functions of tRFs are poorly understood, but tRFs can perform as siRNAs degrading target transcripts (Maute et al., 2013), and regulate ribosomal loading and protein chain elongation (Sobala and Hutvagner, 2013). The biosynthesis of tRFs involves degradation of pre-tRNA molecules via the nuclear TRAMP pathway (Anderson, 2005; Maraia and Lamichhane, 2011) or cytosolic degradation of mature tRNAs via the rapid tRNA decay pathway. Analysis of the small RNA sequencing dataset in this study showed a differing abundance of tRFs following RSV or MeV infection. Analysis of tRFs showed that among the most abundant were 5′tRFs derived from tRNAs for glycine and valine (Figure 3), suggesting that tRNAs for these two amino acids were targeted during RSV and MeV infection. Levels of 5′tRFs for glycine and valine were substantially increased at 24 hpi following RSV or MeV infection. Frequencies of glycine and valine 5′tRFs increased with time, i.e., from 12 to 24 hpi in RSV- and MeV-infected SHS cells. Analysis of the anti-codon sequence for these tRFs showed that GCC, AAC, and CAC anticodons are enriched in 5′tRFs in SHS cells following RSV or MeV infection (Figure 4). A previous report showed that 5′tRF-Glu-CTC produced during RSV infection of A549 cells interacted with apolipoprotein E receptor 2 (APOER2) mRNA to inhibit translation of the gene, favoring RSV replication (Deng et al., 2015).
FIGURE 3
FIGURE 4
We postulated that tRFs had a direct effect on translation and compared tRF profiles in RSV- or MeV-infected SHS cells with synonymous codon usage. Our analysis compared how frequently each codon in the RSV or MeV genome is utilized relative to the standard genetic code frequencies that are given as a codon adaptation index (CAI). CAI ranges from 0 to 1.0, with 0 representing least used, and 1.0 representing most common codons. A comparison of RSCU profiles for RSV and MeV genomes using Mega X (Kumar et al., 2018) showed the codon usage to be considerably different than standard human codon usage frequencies (Figures 5A,B). The assessment showed that of the four possible codons that encode glycine (GGU, GGC, GGA, GGG), three codons (GGC, GGU, GGG) are under-utilized for RSV genes, and valine, RSV uses the GUC codon less frequently. For MeV, two of the four codons for glycine (GGU, GGC) and GUA for valine are under-utilized. The consequence of these findings is unknown but suggests that the tRNAs used during RSV or MeV infection are selectively processed to produce these tRNAs potentially for two reasons—to reduce host translation (Kim et al., 2017, 2019; Mesitov et al., 2017), and/or to function as siRNAs and silence host genes (Emara et al., 2010; Wang et al., 2013; Kuscu et al., 2018; Jehn et al., 2020). For valine, and the increased abundance of valine 5′tRFs (AAC and CAC), the corresponding codon usage was not low suggesting that valine tRFs might have functions other than reducing translation rates. RSCU analysis of individual RSV and MeV genes showed that glycine and valine tRNAs were under-utilized following RSV or MeV infection in SHS cells (Table 3). Together these findings suggest that tRNA cleavage is linked with codon usage. The sncRNA results show that the expression profiles for small RNAs differ in kinetics of expression for RSV and MeV infection of SHS cells.
FIGURE 5
TABLE 3
| Codon | RSV_NS1 | RSV_NS2 | RSV_NP | RSV_P | RSV_M | RSV_SH | RSV_G | RSV_F | RSV_M2 | RSV_L |
| GUU (V) | 0.8 | 0 | 1.2 | 0.9 | 0.7 | 0 | 0 | 2.0 | 1.3 | 1.6 |
| GUG (V) | 1.2 | 0 | 1.0 | 2.2 | 1.4 | 0 | 0.5 | 0.4 | 0.8 | 0.6 |
| GGC (G) | 1.14 | 0 | 0.8 | 0 | 0.7 | 0 | 0.7 | 0.7 | 0 | 0.5 |
| GGG (G) | 0.57 | 2.0 | 0.8 | 1.7 | 0.7 | 4.0 | 0.7 | 0.3 | 0.4 | 0.5 |
| Codon | MeV_N | MeV_P | MeV_M | MeV_F | MeV_H | MeV_L | ||||
| GUU (V) | 2.0 | 0.9 | 0.6 | 0.9 | 0.6 | 1.9 | ||||
| GUG (V) | 0.2 | 2.0 | 2.3 | 0.65 | 1.6 | 0.7 | ||||
| GGC (G) | 0.7 | 0.6 | 1.1 | 0.8 | 0.5 | 1.0 | ||||
| GGG (G) | 0.5 | 1.3 | 0.4 | 1.4 | 1.8 | 0.8 | ||||
RSCU indices for selected codons for MeV and RSV genes.
RSCU scores for individual RSV and MeV genes were determined using MEGA X (Kumar et al., 2018).
RSV_NS1, RSV Non-structural protein 1; RSV_NS2, RSV Non-structural gene 2; RSV_NP, RSV Nucleoprotein; RSV_P, RSV Phosphoprotein; RSV_M, RSV Matrix protein; RSV_SH, RSV Small hydrophobic protein; RSV_G, RSV glycoprotein; RSV_F, RSV Fusion protein; RSV_M2, RSV Matrix ORF2 protein; RSV_L, RSV Large polymerase protein; MeV_N, Measles virus Nucleocapsid protein; MeV_P, Measles virus Phosphoprotein; MeV_M, Measles virus Matrix protein; MeV_F, Measles virus fusion glycoprotein; MeV_H, Measles virus hemagglutinin; MeV_L, Measles virus large polymerase.
Discussion
In this study, we examined differentially expressed sncRNAs following RSV or MeV infection of SHS cells at 12 and 24 hpi using NGS. NGS showed unique miR expression profiles following RSV or MeV infection (Figure 1 and Table 2). We identified that (1) miR-10b-5p was expressed in mock-treated cells but showed reduced expression in RSV- or MeV-infected cells, (2) p-miR-247 and miR-375 had robust expression in RSV- or MeV-infected cells, (3) miRs had no/low expression in mock-treated cells, but were increased by RSV or MeV infection, and (4) miRs were not expressed in mock-infected cells at 12 hpi but showed robust copy numbers at 24 hpi. These findings show that RSV or MeV infection induces distinct miR expression in SHS cells that is virus-type and time-point specific. While the sequencing data from our study has identified these deregulated miRNAs, the mechanisms contributing to their biogenesis remains unknown and is out of the scope of this manuscript.
Evidence suggests that epitranscriptomics (or post-transcriptional modification of cellular RNAs) especially 3′ end uridylation and adenylation of miRs can alter both miRNA stability and target repertoire (Vitsios and Enright, 2015). We explored the nature of epitranscriptional miR modifications following RSV or MeV infection. The findings suggest that RSV and MeV induce sncRNA modifications affecting target cell specificity. The results showed extensive 3′ modifications (mostly uridylation and adenylation) in the miRs following RSV or MeV infection that may alter host mRNA or miRNA transcript stability. Further studies are needed to elucidate the function of these epitranscriptional modifications on miRNA function.
Post-transcriptional regulatory mechanisms control translation of >60% of the human transcriptome (Friedman et al., 2009). The identification of cleaved tRNAs during RSV or MeV infection of SHS cells is novel. The tRFs identified contained exclusively of 5′tRFs during later stages of RSV or MeV infection, and of these tRNAs, glycine, and valine were the primary targets for 5′tRF formation. Other studies showed that RSV infection can induce angiogenin (ANG)-mediated cleavage of tRNAs to produce tRFs (Deng et al., 2015; Zhou J. et al., 2017). These studies identified a pro-viral role for 5′tRF-GlyCCC, 5′tRF-LysCTT, and 5′tRF-GluCTC by depression of APOER2 mRNA. We postulate that tRNA cleavage to tRFs may have a broader anti-viral role. Degradation of host tRNAs may reduce the anti-viral protein response, and recent findings from proteomic studies in RSV-infected A549 cells have shown that the host nuclear proteins are considerably reduced following RSV infection (Munday et al., 2012, 2015). Similar results have also been shown for MeV, an attribute linked to the N protein (Sato et al., 2007; Okonski and Samuel, 2013). It is important to investigate why RSV or MeV infection induces selective cleavage of glycine and valine tRNAs as glycine, and valine tRNAs are selectively under-utilized by most RSV and MeV genes (Table 3). Transfer RNA fragments operate in mucosal immunity (Chen and Shen, 2021) regulating retrotransposon expression (Schorn et al., 2017) and viral adaptation to the host (Lauring et al., 2012; Pavon-Eternod et al., 2013; Rahman et al., 2018b; Le Nouën et al., 2020; Nunes et al., 2020).
The degradation of tRNAs may be a host stress response to viral infection, however, the increases in tRNA for glycine and valine hint that MeV and RSV may potentially target these tRNAs. The mechanisms that cause are unknown and are difficult to resolve in the context of this study. Currently, we do not understand the mechanisms involved in sncRNA deregulation, nor the functional outcome, thus further studies are necessary. Translation of mRNAs linked to the anti-viral response is key to control of virus infection and clearance. The findings show that RSV or MeV infection induces several classes of sncRNAs that regulate post-transcriptional gene expression and emphasizes the need to better understand the host-virus interface.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The data presented in the study are deposited in Gene Expression Omnibus (GEO) with accession numbers GSE167724 and GSM5112024–GSM-5112037.
Author contributions
AB contributed to conception, execution, and data analysis, and wrote manuscript drafts and revisions. CD contributed to execution. SC and RT contributed to study design, funding, data interpretation, manuscript writing, and review. HA contributed to validating infection using dark field fluorescence microscopy. All authors have seen the draft of the manuscript.
Funding
The studies were supported by a grant awarded to SC and RT (IE140652) to study “Paramyxovirus miRNA biogenesis and implications for host and viral gene expression” and the Georgia Research Alliance Foundation for RT.
Acknowledgments
The authors would like to acknowledge the kind gift of rgRSV224GFP from Michael Teng at University of South Florida (USF).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.671852/full#supplementary-material
Supplementary Figure 1RSV infects SHS neuronal cells. SHS cells were infected with rgRSV-GFP at MOI = 0.1 for indicated time points. Viral replication was visualized by dark field fluorescence microscopy on a Nikon UV Eclipse TE2000-U microscope at 10X magnification.
Supplementary Table 1MiRNA modification data patterns in RSV and MeV infected neuronal cells. Table lists samples, microRNAs, type of modifications, modification site, pattern, modification position, internal modification type, internal modification position, repeat modifications and counts of modifications for all deregulated miRNAs identified in the study.
References
1
AdvaniV. M.IvanovP. (2019). Translational Control under Stress: Reshaping the Translatome.Bioessays41:e1900009.
2
AlarconC. R.LeeH.GoodarziH.HalbergN.TavazoieS. F. (2015). N6-methyladenosine marks primary microRNAs for processing.Nature519482–485. 10.1038/nature14281
3
AndersonJ. T. (2005). RNA turnover: unexpected consequences of being tailed.Curr. Biol.15R635–R638.
4
BakreA. A.HarcourtJ. L.HaynesL. M.AndersonL. J.TrippR. A. (2017). The Central Conserved Region (CCR) of Respiratory Syncytial Virus (RSV) G Protein Modulates Host miRNA Expression and Alters the Cellular Response to Infection.Vaccines5:16. 10.3390/vaccines5030016
5
BakreA. A.MalekiA.TrippR. A. (2019). MicroRNA and Nonsense Transcripts as Putative Viral Evasion Mechanisms.Front. Cell Infect. Microbiol.9:152. 10.3389/fcimb.2019.00152
6
BakreA.MitchellP.ColemanJ. K.JonesL. P.SaavedraG.TengM.et al (2012). Respiratory syncytial virus modifies microRNAs regulating host genes that affect virus replication.J. Gen. Virol.932346–2356. 10.1099/vir.0.044255-0
7
BakreA.WuW.HiscoxJ.SpannK.TengM. N.TrippR. A. (2015). Human respiratory syncytial virus non-structural protein NS1 modifies miR-24 expression via transforming growth factor-beta.J. Gen. Virol.963179–3191. 10.1099/jgv.0.000261
8
BartelD. P. (2009). MicroRNAs: Target Recognition and Regulatory Functions.Cell136215–233. 10.1016/j.cell.2009.01.002
9
BergeronH. C.MurrayJ.Nunez CastrejonA. M.DuboisR. M.TrippR. A. (2021). Respiratory Syncytial Virus (RSV) G Protein Vaccines With Central Conserved Domain Mutations Induce CX3C-CX3CR1 Blocking Antibodies.Viruses13:352. 10.3390/v13020352
10
BernierA.SaganS. M. (2018). The Diverse Roles of microRNAs at the Host(-)Virus Interface.Viruses10:440. 10.3390/v10080440
11
BohmwaldK.GalvezN. M. S.RiosM.KalergisA. M. (2018). Neurologic Alterations Due to Respiratory Virus Infections.Front. Cell Neurosci.12:386. 10.3389/fncel.2018.00386
12
BriandJ.SérandourA. A.NadaradjaneA.Bougras-CartronG.HeymannD.OryB.et al (2020). N6-Adenosine Methylation of miRNA-200b-3p Influences Its Functionality and Is a Theranostic Tool.Mol. Ther. Nucleic Acids2272–83. 10.1016/j.omtn.2020.08.010
13
Cakmak GencG.DursunA.Karakas CelikS.CalikM.KokturkF.PiskinI. E. (2018). IL28B, IL29 and micro-RNA 548 in subacute sclerosing panencephalitis as a rare disease.Gene67873–78. 10.1016/j.gene.2018.07.062
14
ChenY.ShenJ. (2021). Mucosal immunity and tRNA, tRF, and tiRNA.J. Mol. Med.9947–56. 10.1007/s00109-020-02008-4
15
CherayM.EtcheverryA.JacquesC.PacaudR.Bougras-CartronG.AubryM.et al (2020). Cytosine methylation of mature microRNAs inhibits their functions and is associated with poor prognosis in glioblastoma multiforme.Mol. Cancer19:36.
16
ChoiY.MasonC. S.JonesL. P.CrabtreeJ.JorqueraP. A.TrippR. A. (2012). Antibodies to the central conserved region of respiratory syncytial virus (RSV) G protein block RSV G protein CX3C-CX3CR1 binding and cross-neutralize RSV A and B strains.Viral Immunol.25193–203.
17
Correia De SousaM.GjorgjievaM.DolickaD.SobolewskiC.FotiM. (2019). Deciphering miRNAs’ Action through miRNA Editing.Int. J. Mol. Sci.20:6249. 10.3390/ijms20246249
18
CullenB. R. (2013). MicroRNAs as mediators of viral evasion of the immune system.Nat. Immunol.14205–210. 10.1038/ni.2537
19
DahiyaN.AtreyaC. D. (2014). MicroRNAs and Major Blood-borne Infectious Viral Diseases.Microrna2212–218. 10.2174/2211536602666131118224225
20
DakhamaA.VitalisT. Z.HegeleR. G. (1997). Persistence of respiratory syncytial virus (RSV) infection and development of RSV-specific IgG1 response in a guinea-pig model of acute bronchiolitis.Eur. Respir. J.1020–26. 10.1183/09031936.97.10010020
21
DengJ.PtashkinR. N.ChenY.ChengZ.LiuG.PhanT.et al (2015). Respiratory Syncytial Virus Utilizes a tRNA Fragment to Suppress Antiviral Responses Through a Novel Targeting Mechanism.Mol. Ther.231622–1629. 10.1038/mt.2015.124
22
EisenhutM. (2006). Extrapulmonary manifestations of severe respiratory syncytial virus infection–a systematic review.Crit. Care10:R107.
23
EmaraM. M.IvanovP.HickmanT.DawraN.TisdaleS.KedershaN.et al (2010). Angiogenin-induced tRNA-derived stress-induced RNAs promote stress-induced stress granule assembly.J. Biol. Chem.28510959–10968. 10.1074/jbc.m109.077560
24
FerrenM.HorvatB.MathieuC. (2019). Measles Encephalitis: Towards New Therapeutics.Viruses11:1017. 10.3390/v11111017
25
FriedmanR. C.FarhK. K.BurgeC. B.BartelD. P. (2009). Most mammalian mRNAs are conserved targets of microRNAs.Genome Res.1992–105. 10.1101/gr.082701.108
26
GanesanP.ChandwaniM. N.CreisherP. S.BohnL.O’donnellL. A. (2018). The neonatal anti-viral response fails to control measles virus spread in neurons despite interferon-gamma expression and a Th1-like cytokine profile.J. Neuroimmunol.31680–97. 10.1016/j.jneuroim.2017.12.018
27
GrassmannR.JeangK. T. (2008). The roles of microRNAs in mammalian virus infection.Biochim. Biophys. Acta1779706–711. 10.1016/j.bbagrm.2008.05.005
28
GriffinD. E. (2018). Measles Vaccine.Viral Immunol.3186–95.
29
GriffinD. E.LinW. W.NelsonA. N. (2018). Understanding the causes and consequences of measles virus persistence.F1000Res7:237. 10.12688/f1000research.12094.1
30
HaberN. (2018). Respiratory syncytial virus infection in elderly adults.Med. Mal. Infect.48377–382. 10.1016/j.medmal.2018.01.008
31
HallakL. K.KwilasS. A.PeeplesM. E. (2007). Interaction between respiratory syncytial virus and glycosaminoglycans, including heparan sulfate.Methods Mol. Biol.37915–34. 10.1385/1-59745-393-5:15
32
HallakL. K.SpillmannD.CollinsP. L.PeeplesM. E. (2000). Glycosaminoglycan sulfation requirements for respiratory syncytial virus infection.J. Virol.7410508–10513. 10.1128/jvi.74.22.10508-10513.2000
33
HaralambievaI. H.KennedyR. B.SimonW. L.GoergenK. M.GrillD. E.OvsyannikovaI. G.et al (2018). Differential miRNA expression in B cells is associated with inter-individual differences in humoral immune response to measles vaccination.PLoS One13:e0191812. 10.1371/journal.pone.0191812
34
HarcourtJ. L.KarronR. A.TrippR. A. (2004). Anti-G protein antibody responses to respiratory syncytial virus infection or vaccination are associated with inhibition of G protein CX3C-CX3CR1 binding and leukocyte chemotaxis.J. Infect. Dis.1901936–1940. 10.1086/425516
35
HardieD. R.AlbertynC.HeckmannJ. M.SmutsH. E. (2013). Molecular characterisation of virus in the brains of patients with measles inclusion body encephalitis (MIBE).Virol J.10:283.
36
HasegawaK.Perez-LosadaM.HoptayC. E.EpsteinS.MansbachJ. M.TeachS. J.et al (2018). RSV vs. rhinovirus bronchiolitis: difference in nasal airway microRNA profiles and NFkappaB signaling.Pediatr. Res.83606–614. 10.1038/pr.2017.309
37
HaynesL. M.JonesL. P.BarskeyA.AndersonL. J.TrippR. A. (2003). Enhanced disease and pulmonary eosinophilia associated with formalin-inactivated respiratory syncytial virus vaccination are linked to G glycoprotein CX3C-CX3CR1 interaction and expression of substance P.J. Virol.779831–9844. 10.1128/jvi.77.18.9831-9844.2003
38
InchleyC. S.SonerudT.FjaerliH. O.NakstadB. (2015). Nasal mucosal microRNA expression in children with respiratory syncytial virus infection.BMC Infect. Dis.15:150. 10.1186/s12879-015-0878-z
39
IvanovP.EmaraM. M.VillenJ.GygiS. P.AndersonP. (2011). Angiogenin-induced tRNA fragments inhibit translation initiation.Mol. Cell43613–623. 10.1016/j.molcel.2011.06.022
40
JehnJ.TremlJ.WulschS.OttumB.ErbV.HewelC.et al (2020). 5’ tRNA halves are highly expressed in the primate hippocampus and might sequence-specifically regulate gene expression.RNA26694–707. 10.1261/rna.073395.119
41
KanehisaM.FurumichiM.TanabeM.SatoY.MorishimaK. (2017). KEGG: new perspectives on genomes, pathways, diseases and drugs.Nucleic Acids Res.45D353–D361.
42
KanehisaM.SatoY.KawashimaM.FurumichiM.TanabeM. (2016). KEGG as a reference resource for gene and protein annotation.Nucleic Acids Res.44D457–D462.
43
KaragkouniD.ParaskevopoulouM. D.ChatzopoulosS.VlachosI. S.TastsoglouS.KanellosI.et al (2018). DIANA-TarBase v8: a decade-long collection of experimentally supported miRNA-gene interactions.Nucleic Acids Res.46D239–D245.
44
KimH. K.FuchsG.WangS.WeiW.ZhangY.ParkH.et al (2017). A transfer-RNA-derived small RNA regulates ribosome biogenesis.Nature55257–62. 10.1038/nature25005
45
KimH. K.XuJ.ChuK.ParkH.JangH.LiP.et al (2019). A tRNA-Derived Small RNA Regulates Ribosomal Protein S28 Protein Levels after Translation Initiation in Humans and Mice.Cell Rep.293816.e–3824.e.
46
KingK. A.HuC.RodriguezM. M.RomagueraR.JiangX.PiedimonteG. (2001). Exaggerated neurogenic inflammation and substance P receptor upregulation in RSV-infected weanling rats.Am. J. Respir. Cell Mol. Biol.24101–107. 10.1165/ajrcmb.24.2.4264
47
KumarS.StecherG.LiM.KnyazC.TamuraK. (2018). MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms.Mol. Biol. Evol.351547–1549. 10.1093/molbev/msy096
48
KuscuC.KumarP.KiranM.SuZ.MalikA.DuttaA. (2018). tRNA fragments (tRFs) guide Ago to regulate gene expression post-transcriptionally in a Dicer-independent manner.RNA241093–1105. 10.1261/rna.066126.118
49
LauringA. S.AcevedoA.CooperS. B.AndinoR. (2012). Codon usage determines the mutational robustness, evolutionary capacity, and virulence of an RNA virus.Cell Host Microbe12623–632. 10.1016/j.chom.2012.10.008
50
LayM. K.BuenoS. M.GalvezN.RiedelC. A.KalergisA. M. (2016). New insights on the viral and host factors contributing to the airway pathogenesis caused by the respiratory syncytial virus.Crit. Rev. Microbiol.42800–812.
51
Le NouënC.LuongoC. L.YangL.MuellerS.WimmerE.DinapoliJ. M.et al (2020). Optimization of the Codon Pair Usage of Human Respiratory Syncytial Virus Paradoxically Resulted in Reduced Viral Replication In Vivo and Reduced Immunogenicity.J. Virol.94e1296–e1219.
52
LeungA. K.HonK. L.LeongK. F.SergiC. M. (2018). Measles: a disease often forgotten but not gone.Hong Kong Med. J.24512–520.
53
LiX. Q.FuZ. F.AlvarezR.HendersonC.TrippR. A. (2006). Respiratory syncytial virus (RSV) infects neuronal cells and processes that innervate the lung by a process involving RSV G protein.J. Virol.80537–540. 10.1128/jvi.80.1.537-540.2006
54
LinL. T.RichardsonC. D. (2016). The Host Cell Receptors for Measles Virus and Their Interaction with the Viral Hemagglutinin (H) Protein.Viruses8:250. 10.3390/v8090250
55
LongX.XieJ.ZhaoK.LiW.TangW.ChenS.et al (2016). NK cells contribute to persistent airway inflammation and AHR during the later stage of RSV infection in mice.Med. Microbiol. Immunol.205459–470. 10.1007/s00430-016-0459-9
56
LoveM. I.HuberW.AndersS. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2.Genome Biol.15:550.
57
MaraiaR. J.LamichhaneT. N. (2011). 3’ processing of eukaryotic precursor tRNAs.Wiley Interdiscip. Rev. RNA2362–375. 10.1002/wrna.64
58
MauteR. L.SchneiderC.SumazinP.HolmesA.CalifanoA.BassoK.et al (2013). tRNA-derived microRNA modulates proliferation and the DNA damage response and is down-regulated in B cell lymphoma.Proc. Natl. Acad. Sci. U S A.1101404–1409. 10.1073/pnas.1206761110
59
MazurN. I.HigginsD.NunesM. C.MeleroJ. A.LangedijkA. C.HorsleyN.et al (2018). The respiratory syncytial virus vaccine landscape: lessons from the graveyard and promising candidates.Lancet Infect. Dis.18e295–e311.
60
McQuaidS.CosbyS. L.KoffiK.HondeM.KirkJ.LucasS. B. (1998). Distribution of measles virus in the central nervous system of HIV-seropositive children.Acta Neuropathol.96637–642. 10.1007/s004010050945
61
MesitovM. V.SoldatovR. A.ZaichenkoD. M.MalakhoS. G.KlementyevaT. S.SokolovskayaA. A.et al (2017). Differential processing of small RNAs during endoplasmic reticulum stress.Sci. Rep.7:46080.
62
MirraV.UllmannN.CherchiC.OnofriA.PagliettiM. G.CutreraR. (2018). Respiratory syncytial virus prophylaxis and the “special population”.Minerva Pediatr.70589–599.
63
MishraB.WangG. (2012). The Importance of Amino Acid Composition in Natural AMPs: An Evolutional, Structural, and Functional Perspective.Front. Immunol.3:221. 10.3389/fimmu.2012.00221
64
MorichiS.MorishitaN.IshidaY.OanaS.YamanakaG.KashiwagiY.et al (2017). Examination of neurological prognostic markers in patients with respiratory syncytial virus-associated encephalopathy.Int. J. Neurosci.12744–50. 10.3109/00207454.2016.1138951
65
MundayD. C.HowellG.BarrJ. N.HiscoxJ. A. (2015). Proteomic analysis of mitochondria in respiratory epithelial cells infected with human respiratory syncytial virus and functional implications for virus and cell biology.J. Pharm. Pharmacol.67300–318.
66
MundayD. C.SurteesR.EmmottE.DoveB. K.DigardP.BarrJ. N.et al (2012). Using SILAC and quantitative proteomics to investigate the interactions between viral and host proteomes.Proteomics12666–672. 10.1002/pmic.201100488
67
NaamanH.RallG.MatulloC.Veksler-LublinskyI.Shemer-AvniY.GopasJ. (2017). MiRNA-124 is a link between measles virus persistent infection and cell division of human neuroblastoma cells.PLoS One12:e0187077. 10.1371/journal.pone.0187077
68
NunesA.RibeiroD. R.MarquesM.SantosM. A. S.RibeiroD.SoaresA. R. (2020). Emerging Roles of tRNAs in RNA Virus Infections.Trends Biochem. Sci.45794–805. 10.1016/j.tibs.2020.05.007
69
OkonskiK. M.SamuelC. E. (2013). Stress granule formation induced by measles virus is protein kinase PKR dependent and impaired by RNA adenosine deaminase ADAR1.J. Virol.87756–766. 10.1128/jvi.02270-12
70
OmarS.ClarkeR.AbdullahH.BradyC.CorryJ.WinterH.et al (2017). Respiratory virus infection up-regulates TRPV1, TRPA1 and ASICS3 receptors on airway cells.PLoS One12:e0171681. 10.1371/journal.pone.0171681
71
OshanskyC. M.KrunkoskyT. M.BarberJ.JonesL. P.TrippR. A. (2009). Respiratory syncytial virus proteins modulate suppressors of cytokine signaling 1 and 3 and the type I interferon response to infection by a toll-like receptor pathway.Viral Immunol.22147–161. 10.1089/vim.2008.0098
72
PaesB. (2018). Respiratory Syncytial Virus in Otherwise Healthy Prematurely Born Infants: A Forgotten Majority.Am. J. Perinatol.35541–544. 10.1055/s-0038-1637762
73
Pavon-EternodM.DavidA.DittmarK.BerglundP.PanT.BenninkJ. R.et al (2013). Vaccinia and influenza A viruses select rather than adjust tRNAs to optimize translation.Nucleic Acids Res.411914–1921. 10.1093/nar/gks986
74
PetrovA. I.KayS. J. E.KalvariI.HoweK. L.GrayK. A.BrufordE. A.et al (2017). RNAcentral: a comprehensive database of non-coding RNA sequences.Nucleic Acids Res.45D128–D134.
75
PirouzM.EbrahimiA. G.GregoryR. I. (2019). Unraveling 3’-end RNA uridylation at nucleotide resolution.Methods15510–19. 10.1016/j.ymeth.2018.10.024
76
RahmanR.-U.GautamA.BethuneJ.SattarA.FiosinsM.MagruderD. S.et al (2018a). Oasis 2: improved online analysis of small RNA-seq data.BMC Bioinformatics19:54. 10.1186/s12859-018-2047-z
77
RahmanS. U.YaoX.LiX.ChenD.TaoS. (2018b). Analysis of codon usage bias of Crimean-Congo hemorrhagic fever virus and its adaptation to hosts.Infect. Genet. Evol.581–16. 10.1016/j.meegid.2017.11.027
78
RiddellM. A.MossW. J.HauerD.MonzeM.GriffinD. E. (2007). Slow clearance of measles virus RNA after acute infection.J. Clin. Virol.39312–317. 10.1016/j.jcv.2007.05.006
79
RimaB. K.DuprexW. P. (2005). Molecular mechanisms of measles virus persistence.Virus Res.111132–147. 10.1016/j.virusres.2005.04.005
80
RoizmanB.ZhouG.DuT. (2011). Checkpoints in productive and latent infections with herpes simplex virus 1: conceptualization of the issues.J. Neurovirol.17512–517. 10.1007/s13365-011-0058-x
81
RotaP. A.RotaJ. S.GoodsonJ. L. (2017). Subacute Sclerosing Panencephalitis.Clin. Infect. Dis.65233–234.
82
RussellC. J.SimoesE. A. F.HurwitzJ. L. (2018). Vaccines for the Paramyxoviruses and Pneumoviruses: Successes, Candidates, and Hurdles.Viral Immunol.31133–141. 10.1089/vim.2017.0137
83
SatoH.MasudaM.KanaiM.Tsukiyama-KoharaK.YonedaM.KaiC. (2007). Measles virus N protein inhibits host translation by binding to eIF3-p40.J. Virol.8111569–11576. 10.1128/jvi.00570-07
84
SatoY.WatanabeS.FukudaY.HashiguchiT.YanagiY.OhnoS. (2018). Cell-to-Cell Measles Virus Spread between Human Neurons Is Dependent on Hemagglutinin and Hyperfusogenic Fusion Protein.J. Virol.92e2166–e2117.
85
SchornA. J.GutbrodM. J.LeblancC.MartienssenR. (2017). LTR-Retrotransposon Control by tRNA-Derived Small RNAs.Cell17061.e–71.e.
86
SchwarzeJ.O’donnellD. R.RohwedderA.OpenshawP. J. (2004). Latency and persistence of respiratory syncytial virus despite T cell immunity.Am. J. Respir. Crit. Care Med.169801–805. 10.1164/rccm.200308-1203oc
87
ShenY.YuX.ZhuL.LiT.YanZ.GuoJ. (2018). Transfer RNA-derived fragments and tRNA halves: biogenesis, biological functions and their roles in diseases.J. Mol. Med.961167–1176. 10.1007/s00109-018-1693-y
88
SinghB. K.LiN.MarkA. C.MateoM.CattaneoR.SinnP. L. (2016). Cell-to-Cell Contact and Nectin-4 Govern Spread of Measles Virus from Primary Human Myeloid Cells to Primary Human Airway Epithelial Cells.J. Virol.906808–6817.
89
SobalaA.HutvagnerG. (2013). Small RNAs derived from the 5’ end of tRNA can inhibit protein translation in human cells.RNA Biol.10553–563. 10.4161/rna.24285
90
SofosN.XuK.DedicE.BrodersenD. E. (2015). Cut to the chase–Regulating translation through RNA cleavage.Biochimie11410–17. 10.1016/j.biochi.2015.01.009
91
TanY. R.YangT.LiuS. P.XiangY.QuF.LiuH. J.et al (2008). Pulmonary peptidergic innervation remodeling and development of airway hyperresponsiveness induced by RSV persistent infection.Peptides2947–56. 10.1016/j.peptides.2007.10.020
92
TaniguchiM.YanagiY.OhnoS. (2019). Both type I and type III interferons are required to restrict measles virus growth in lung epithelial cells.Arch. Virol.164439–446. 10.1007/s00705-018-4087-0
93
TatsuoH.OnoN.TanakaK.YanagiY. (2000). SLAM (CDw150) is a cellular receptor for measles virus.Nature406893–897. 10.1038/35022579
94
TechaarpornkulS.CollinsP. L.PeeplesM. E. (2002). Respiratory syncytial virus with the fusion protein as its only viral glycoprotein is less dependent on cellular glycosaminoglycans for attachment than complete virus.Virology294296–304. 10.1006/viro.2001.1340
95
TrippR. A.BarskeyA.GossL.AndersonL. J. (2002). Substance P receptor expression on lymphocytes is associated with the immune response to respiratory syncytial virus infection.J. Neuroimmunol.129141–153. 10.1016/s0165-5728(02)00169-8
96
TrippR. A.DakhamaA.JonesL. P.BarskeyA.GelfandE. W.AndersonL. J. (2003). The G glycoprotein of respiratory syncytial virus depresses respiratory rates through the CX3C motif and substance P.J. Virol.776580–6584. 10.1128/jvi.77.11.6580-6584.2003
97
TrippR. A.JonesL. P.HaynesL. M.ZhengH.MurphyP. M.AndersonL. J. (2001). CX3C chemokine mimicry by respiratory syncytial virus G glycoprotein.Nat. Immunol.2732–738. 10.1038/90675
98
TrippR. A.MooreD.WinterJ.AndersonL. J. (2000). Respiratory syncytial virus infection and G and/or SH protein expression contribute to substance P, which mediates inflammation and enhanced pulmonary disease in BALB/c mice.J. Virol.741614–1622. 10.1128/jvi.74.4.1614-1622.2000
99
VitsiosD. M.EnrightA. J. (2015). Chimira: analysis of small RNA sequencing data and microRNA modifications.Bioinformatics313365–3367. 10.1093/bioinformatics/btv380
100
VlachosI. S.ZagganasK.ParaskevopoulouM. D.GeorgakilasG.KaragkouniD.VergoulisT.et al (2015). DIANA-miRPath v3.0: deciphering microRNA function with experimental support.Nucleic Acids Res.43W460–W466.
101
WangJ. X.GaoJ.DingS. L.WangK.JiaoJ. Q.WangY.et al (2015). Oxidative Modification of miR-184 Enables It to Target Bcl-xL and Bcl-w.Mol. Cell5950–61. 10.1016/j.molcel.2015.05.003
102
WangQ.LeeI.RenJ.AjayS. S.LeeY. S.BaoX. (2013). Identification and functional characterization of tRNA-derived RNA fragments (tRFs) in respiratory syncytial virus infection.Mol. Ther.21368–379. 10.1038/mt.2012.237
103
WangS.LiuP.YangP.ZhengJ.ZhaoD. (2017). Peripheral blood microRNAs expression is associated with infant respiratory syncytial virus infection.Oncotarget896627–96635. 10.18632/oncotarget.19364
104
YisU.TufekciU. K.GencS.CarmanK. B.BayramE.TopcuY.et al (2015). Expression patterns of micro-RNAs 146a, 181a, and 155 in subacute sclerosing panencephalitis.J. Child Neurol.3069–74. 10.1177/0883073814531329
105
ZhongQ.FengH.LuQ.LiuX.ZhaoQ.DuY.et al (2018). Recurrent wheezing in neonatal pneumonia is associated with combined infection with Respiratory Syncytial Virus and Staphylococcus aureus or Klebsiella pneumoniae.Sci. Rep.8:995.
106
ZhouJ.LiuS.ChenY.FuY.SilverA. J.HillM. S.et al (2017). Identification of two novel functional tRNA-derived fragments induced in response to respiratory syncytial virus infection.J. Gen. Virol.981600–1610. 10.1099/jgv.0.000852
107
ZhouN.LiW.RenL.XieX.LiuE. (2017). An Interaction of LPS and RSV Infection in Augmenting the AHR and Airway Inflammation in Mice.Inflammation401643–1653. 10.1007/s10753-017-0604-7
Summary
Keywords
respiratory syncytial virus, measles virus, neuronal cells, microRNAs, piwi-associated RNAs, transfer RNAs
Citation
Bakre AA, Duffy C, Abdullah H, Cosby SL and Tripp RA (2021) Small Non-coding RNA Expression Following Respiratory Syncytial Virus or Measles Virus Infection of Neuronal Cells. Front. Microbiol. 12:671852. doi: 10.3389/fmicb.2021.671852
Received
24 February 2021
Accepted
03 August 2021
Published
03 September 2021
Volume
12 - 2021
Edited by
Jianrong Li, The Ohio State University, United States
Reviewed by
Qiyun Zhu, Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Sciences (CAS), China; Mark Peeples, The Ohio State University, United States
Updates
Copyright
© 2021 Bakre, Duffy, Abdullah, Cosby and Tripp.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Abhijeet A. Bakre, bakre@uga.eduRalph A. Tripp, ratripp@uga.edu
This article was submitted to Virology, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.