Expression profiles of host miRNAs and circRNAs and ceRNA network during Toxoplasma gondii lytic cycle

Abstract

Toxoplasma gondii is an opportunistic protozoan parasite that is highly prevalent in the human population and can lead to adverse health consequences in immunocompromised patients and pregnant women. Noncoding RNAs, such as microRNAs (miRNAs) and circular RNAs (circRNAs), play important regulatory roles in the pathogenesis of many infections. However, the differentially expressed (DE) miRNAs and circRNAs implicated in the host cell response during the lytic cycle of T. gondii are unknown. In this study, we profiled the expression of miRNAs and circRNAs in human foreskin fibroblasts (HFFs) at different time points after T. gondii infection using RNA sequencing (RNA-seq). We identified a total of 7, 7, 27, 45, 70, 148, 203, and 217 DEmiRNAs and 276, 355, 782, 1863, 1738, 6336, 1229, and 1680 DEcircRNAs at 1.5, 3, 6, 9, 12, 24, 36, and 48 h post infection (hpi), respectively. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses revealed that the DE transcripts were enriched in immune response, apoptosis, signal transduction, and metabolism-related pathways. These findings provide new insight into the involvement of miRNAs and circRNAs in the host response to T. gondii infection.