Fig. 2

Incubation with tylosin alters the transcriptome of Toxoplasma gondii RH tachyzoites. a Heatmap for the Pearson correlation coefficient between samples. b Principal component analysis (PCA) score scatter plot of Toxoplasma samples based on FPKM values. c Unsupervised hierarchical clustering of R-seq data. FPKM value is normalized so that blue represents low expression and reddish color represents high expression. Columns were hierarchically clustered based on a complete linkage using Pearson correlation coefficient as the distance measure. d Volcano plot showing q values (–log10) versus RNA expression ratios between tylosin treatment group and vehicle group (0.1% DMSO). e Gene Ontology (GO) analysis of DE genes. The X-axis label denotes –log10(q value), whereas the Y-axis label represents the corresponding GO terms. In addition, the yellow dot indicates the number of DE genes. f The top 20 significantly enriched KEGG pathways of the DE genes. The X-axis label shows the rich ratio. The Y-axis label shows the KEGG pathway terms. The color of the dots represents log10 (q value), and the size of the dot represents the number of DE genes enriched in the KEGG pathway