Transcriptomics analysis of Toxoplasma gondii-infected mouse macrophages reveals coding and noncoding signatures in the presence and absence of MyD88

Background Toxoplasma gondii is a globally distributed protozoan parasite that establishes life-long asymptomatic infection in humans, often emerging as a life-threatening opportunistic pathogen during immunodeficiency. As an intracellular microbe, Toxoplasma establishes an intimate relationship with its host cell from the outset of infection. Macrophages are targets of infection and they are important in early innate immunity and possibly parasite dissemination throughout the host. Here, we employ an RNA-sequencing approach to identify host and parasite transcriptional responses during infection of mouse bone marrow-derived macrophages (BMDM). We incorporated into our analysis infection with the high virulence Type I RH strain and the low virulence Type II strain PTG. Because the well-known TLR-MyD88 signaling axis is likely of less importance in humans, we examined transcriptional responses in both MyD88(+/+) and MyD88(-/-) BMDM. Long noncoding (lnc) RNA molecules are emerging as key regulators in infection and immunity, and were, therefore, included in our analysis. Results We found significantly more host genes were differentially expressed in response to the highly virulent RH strain rather than with the less virulent PTG strain (335 versus 74 protein coding genes for RH and PTG, respectively). Enriched in these protein coding genes were subsets associated with the immune response as well as cell adhesion and migration. We identified 249 and 83 non-coding RNAs as differentially expressed during infection with RH and PTG strains, respectively. Although the majority of these are of unknown function, one conserved lncRNA termed mir17hg encodes the mir17 microRNA gene cluster that has been implicated in down-regulating host cell apoptosis during T. gondii infection. Only a minimal number of transcripts were differentially expressed between MyD88 knockout and wild type cells. However, several immune genes were among the differences. While transcripts for parasite secretory proteins were amongst the most highly expressed T. gondii genes during infection, no differentially expressed parasite genes were identified when comparing infection in MyD88 knockout and wild type host BMDM. Conclusions The large dataset presented here lays the groundwork for continued studies on both the MyD88-independent immune response and the function of lncRNAs during Toxoplasma gondii infection.

Automated summary

The study identified differential gene expression in host cells infected with highly virulent RH strain versus less virulent PTG strain, with more immune response-related genes being upregulated in response to RH strain. Additionally, non-coding RNAs showed differential expression between infection with RH and PTG strains. The results provide insights into the immune response mechanisms and the role of long noncoding RNAs during Toxoplasma gondii infection.