HRV16 Infection Induces Changes in the Expression of Multiple piRNAs

Human rhinovirus (HRV) is one of the most important cold-causing pathogens in humans. Piwi-interacting RNAs (piRNAs) are a recently discovered class of small non-coding RNAs whose best-understood function is to repress mobile element (ME) activity in animal germline. However, the profile of human/host piRNA during HRV infection is largely unknown. Here we performed high-throughput sequencing of piRNAs from H1-HeLa cells infected with HRV16 at 12 h, 24 h, and 36 h. The results showed that 22,151,664, 24,362,486 and 22,726,546 piRNAs displayed differential expression after HRV16 infection for three time points. A significant differential expression of 21 piRNAs was found in all time points and further verified by RT-qPCR, including 7 known piRNAs and 14 newly found piRNAs. In addition, piRNA prediction was performed on Piano using the SVM algorithm and transposon information. It found that novel_pir78110, novel_pir78107, novel_pir78097, novel_pir78094 and novel_pir76584 are associated with the DNA/hobo of Drosophila, Ac of maize and Tam3 of snapdragon (hAT)-Charlie transposon. The novel_pir97924, novel_pir105705 and novel_pir105700 recognize long interspersed nuclear elements 1 (LINE-1). The novel_pir33182 and novel_pir46604 are related to the long terminal repeat (LTR)/(Endogenous Retrovirus1) ERV1 repetitive element. The novel_pir73855 is related to the LTR/ERVK repetitive element. Both novel_pir70108 and novel_pir70106 are associated with the LTR/ERVL-MaLR repetitive element. The novel_pir15900 is associated with the DNA/hAT-Tip100 repetitive element. Overall, our results indicated that rhinovirus infection could reduce the expression of some piRNAs to facilitate upregulation of LINE-1 transcription or retrotransposons' expression, which is helpful to further explore the mechanism of rhinovirus infection.