Molecular Determinants and Heterogeneity Underlying Host Response to EV-A71 Infection at Single-cell Resolution

The pathogenic human enterovirus EV-A71 has raised serious public health concerns. A hallmark of EV-A71 infection is the distortion of host transcriptomes in favour of viral replication. While high-throughput approaches have been exploited to dissect these gene dysregulations, they do not fully capture molecular perturbations at the single-cell level and in a physiologically relevant context. In this study, we applied a single-cell RNA sequencing approach on infected differentiated enterocyte cells (C2BBe1), which model the gastrointestinal epithelium targeted initially by EV-A71. Our single-cell analysis of EV-A71-infected culture provided several lines of illuminating observations: 1) This systems approach demonstrated extensive cell-to-cell variation in a single culture upon viral infection and delineated transcriptomic differences between the EV-A71-infected and bystander cells. 2) By analysing expression profiles of known EV-A71 receptors and entry facilitation factors, we found that ANXA2 was closely correlated in expression with the viral RNA in the infected population, supporting its role in EV-A71 entry in the enteric cells. 3) We further catalogued dysregulated lncRNAs elicited by EV-A71 infection and demonstrated the functional implication of lncRNA CYTOR in promoting EV-A71 replication. Viewed together, our single-cell transcriptomic analysis illustrated at the single-cell resolution the heterogeneity of host susceptibility to EV-A71 and revealed the involvement of lncRNAs in host antiviral response.