HRV16 Infection Induces Changes in the Expression of Multiple piRNAs.

Abstract

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 12h, 24h, and 36h. The results showed that22,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 knownpiRNAs 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 theexpression of some piRNAs to facilitate upregulation of LINE-1 transcription or retrotransposons' expression, which is helpful to further explore the mechanism of rhinovirus infection.