Cluster analysis of splenocyte microRNAs in the pig reveals key signal regulators of immunomodulation in the host during acute and chronic Toxoplasma gondii infection
Authors: Hou Zhaofeng , Zhang Hui , Xu Kangzhi , Zhu Shifan , Wang Lele , Su Dingzeyang , Liu Jiantao , Su Shijie , Liu Dandan , Huang Siyang , Xu Jinjun , Pan Zhiming , Tao Jianping 
Published: February 2022

Background

Toxoplasma gondiiis an obligate intracellular protozoan parasite that can cause a geographically widespread zoonosis. Our previous splenocyte microRNA profile analyses of pig infected withT. gondiirevealed that the coordination of a large number of miRNAs regulates the host immune response during infection. However, the functions of other miRNAs involved in the immune regulation duringT. gondiiinfection are not yet known.

Methods

Clustering analysis was performed byK-means, self-organizing map (SOM), and hierarchical clustering to obtain miRNA groups with the similar expression patterns. Then, the target genes of the miRNA group in each subcluster were further analyzed for functional enrichment by Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Reactome pathway to recognize the key signaling molecules and the regulatory signatures of the innate and adaptive immune responses of the host duringT. gondiiinfection.

Results

A total of 252 miRNAs were successfully divided into 22 subclusters byK-means clustering (designated as K1–K22), 29 subclusters by SOM clustering (designated as SOM1–SOM29), and six subclusters by hierarchical clustering (designated as H1–H6) based on their dynamic expression levels in the different infection stages. A total of 634, 660, and 477 GO terms, 15, 26, and 14 KEGG pathways, and 16, 15, and 7 Reactome pathways were significantly enriched byK-means, SOM, and hierarchical clustering, respectively. Of note, up to 22 miRNAs mainly showing downregulated expression at 50 days post-infection (dpi) were grouped into one subcluster (namely subcluster H3-K17-SOM1) through the three algorithms. Functional analysis revealed that a large group of immunomodulatory signaling molecules were controlled by the different miRNA groups to regulate multiple immune processes, for instance, IL-1-mediated cellular response and Th1/Th2 cell differentiation partly depending on Notch signaling transduction for subclusters K1 and K2, innate immune response involved in neutrophil degranulation and TLR4 cascade signaling for subcluster K15, B cell activation for subclusters SOM17, SOM1, and SOM25, leukocyte migration, and chemokine activity for subcluster SOM9, cytokine–cytokine receptor interaction for subcluster H2, and interleukin production, chemotaxis of immune cells, chemokine signaling pathway, and C-type lectin receptor signaling pathway for subcluster H3-K17-SOM1.

Conclusions

Cluster analysis of splenocyte microRNAs in the pig revealed key regulatory properties of subcluster miRNA molecules and important features in the immune regulation induced by acute and chronicT. gondiiinfection. These results contribute new insight into the identification of physiological immune responses and maintenance of tolerance in pig spleen tissues duringT. gondiiinfection.

Graphical Abstract