[Study on transcriptome characteristics of respiratory syncytial virus bronchiolitis in children by RNA sequencing]

To explore the biological characteristics related to the pathogenesis and severity of respiratory syncytial virus (RSV) bronchiolitis by RNA sequencing of white blood cells in children with RSV bronchiolitis. This study is a case-control study. A total of 87 children diagnosed with bronchiolitis and RSV antigen positive and/or RSV nucleic acid positive in the pediatric respiratory department of the Second Affiliated Hospital of Wenzhou Medical University from October 2019 to April 2022 were selected as the case group. The case group was divided into three groups based on the condition: mild, moderate, and severe, and there were two groups according to the presence or absence of atopic symptoms: the atopic group and the non-atopic group, forty healthy children in the same period were selected as the control group. The whole blood leukocyte RNA of the children in the case group and the control group was extracted for RNA sequencing, and the data were analyzed to obtain differentially expressed genes (DEGs). Then, the immunobiological pathways and genes related to the pathogenesis, disease condition, and atopy were screened through Gene Ontology (GO) annotation, Kyoto Gene and Genome Encyclopedia (KEGG) annotation, and protein interaction network (PPI) construction methods. Construct the weighted gene co-expression network analysis (WGCNA) module to identify potential biological indicators related to disease severity.Compared with the control group, the case group had a total of 1 782 DEGs, including 1 586 upregulated genes and 196 downregulated genes. The GO pathway enrichment of DEGs is mainly enriched in molecular functions such as peroxidase activity and oxidoreductase activity. In the cytological components, it is mainly enriched in cytoplasmic vesicle lumen and secretory granule lumen. In biological processes, it is mainly enriched in processes such as neutrophil activation involved in immune responses, neutrophil degranulation, and neutrophil activation. KEGG analysis is mainly concentrated in the signal pathway of the viral protein interaction with cytokine and cytokine receptor. A PPI network was constructed to screen four genes at the core position, including CCL2, IL-10, MMP9 and JUN. The DEGs obtained by comparing different disease groups with the control group are mainly enriched in retrograde endocannabinoid signaling and cell apoptosis pathways. WGCNA analysis showed that the brown module related to oxygen saturation was most closely related to the disease, and its gene was mainly enriched in the RNA helicase retinoic acid inducible gene-I (RIG-I) like receptor signal pathway. There are 230 specific DEGs in the atopic group and 444 in the non-atopic group. KEGG enrichment analysis results show that both groups are enriched to NF-κB signaling pathway, the characteristic does not cause significant changes in immune response and transcriptome characteristics in children with RSV bronchiolitis. In conclusion, neutrophil activation, degranulation pathway and signal pathway of interaction between viral protein and cytokine and cytokine receptor are involved in the immune response of RSV bronchiolitis host. CCL2, IL-10, MMP9 and JUN genes may be associated with the pathogenesis. They might be potential biomarkers related to disease severity in RIG-I like receptors, cell apoptosis, and endogenous cannabinoid related signaling pathways.