Fig. 4

Co-expression analysis of genes differentially expressed between patient groups 3 and 4. Coseq was used to generate clusters of genes with positively correlated transcription profiles. 1493 DEGs formed 14 clusters; five clusters with dynamics that clearly distinguish the two patient groups are shown here. a, b Concern transcripts that are abundant in patient group 4 and positively correlated with PvMSP-7 expression. c–e Concern transcripts that were abundant in patient group 3, and were negatively correlated with PvMSP-7l. Boxplots (left) represent expression values for all genes in a cluster (Log2-transformed, normalized read counts) from an individual patient, shaded by patient group; the black line connecting boxplots is the mean expression value of all genes in the cluster. a PvMSP-7K, 7I, 7H and 7C (only 7H is shown) were co-expressed with 91 other genes including various invasion-related genes such as ETRAMP, SEA1, and RON4. b A fifth MSP-7 gene, PvMSP-7L, was co-expressed with 125 other genes including PMV, MOP, and RON5. c A gametocyte stage-specific marker, gamete antigen 27/25, was absent in patient group 4 but positively correlated with TRAG36 and PIR genes. d A cohort of 168 genes that are abundant in patient group 3 but absent in group 4. e Another stage-specific marker, liver stage antigen 3, is abundant in patient group 3 and co-expressed with 204 other genes. Pearson’s correlation coefficient (r) for abundance values of MSP-7H and selected genes are shown