Co-expression gene module analysis in response to attenuated cercaria vaccine reveals a critical role for NK cells in protection against Schistosoma mansoni

Abstract

Background Despite decades of research, an effective schistosomiasis vaccine remains elusive. The radiation-attenuated (RA) cercarial vaccine remains the best model for eliciting high levels of protection. We have recently explored this model in mice to identify potentially protective pathways by examining gene expression patterns in peripheral blood mononuclear cells (PBMC). Methods Herein, we reanalyzed the transcriptomic data from PBMC obtained from vaccinated and infected C57BL/6 mice in three timepoints (Days 7 and 17 after infection or vaccination and Day 7 post-challenge). In addition, we generated new data on PBMC collected 35 days after infection. Deconvolution analysis was performed to estimate immune cell composition by CIBERSORTx. Gene co-expression networks and over-representation analysis (ORA) were performed using the CEMiTool package. Protein-protein interaction networks were constructed using STRING, and the hub proteins for each module were identified using Cytoscape. Results Co-expression network analysis identified a module (M2) associated with the infection process, grouping genes related to a Th2 immune response, and a second module (M6) associated with the vaccination process, displaying pathways related to a Th1 response, CD8 + T cells and NK cells. Within each module, five hub proteins were identified based on protein-protein interaction networks. The M2 infection module revealed Chil3, Il4, Cx3cr1, Emr1 and Ccl2 as hubs, while module M6, associated with vaccination, disclosed Prf1, Klrc1, IFN-γ, Ncr1 and Tbx21 as hub proteins. Conclusions Our data point to the potentiald role of NK cells that may contribute to the RA vaccine response through the production of IFN-γ orchestrated by the T-bet transcription factor (Tbx21). Graphical Abstract