Background
Toxoplasma gondii is an obligate intracellular protozoan parasite that can cause a geographically widespread zoonosis. Our previous splenocyte microRNA profile analyses of pig infected with T. gondii revealed 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 during T. gondii infection are not yet known.
Methods
Clustering analysis was performed by K-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 during T. gondii infection.
Results
A total of 252 miRNAs were successfully divided into 22 subclusters by K-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 by K-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 chronic T. gondii infection. These results contribute new insight into the identification of physiological immune responses and maintenance of tolerance in pig spleen tissues during T. gondii infection.
Graphical Abstract