Trypanosoma cruzi and Toxoplasma gondii Induce a Differential MicroRNA Profile in Human Placental Explants

Abstract

Trypanosoma cruzi and Toxoplasma gondii are two parasites than can be transmitted from mother to child through the placenta. However, congenital transmission rates are low for T. cruzi and high for T. gondii. Infection success or failure depends on complex parasite-host interactions in which parasites can alter host gene expression by modulating non-coding RNAs such as miRNAs. As of yet, there are no reports on altered miRNA expression in placental tissue in response to either parasite. Therefore, we infected human placental explants ex vivo by cultivation with either T. cruzi or T. gondii for 2 h. We then analyzed the miRNA expression profiles of both types of infected tissue by miRNA sequencing and quantitative PCR, sequence-based miRNA target prediction, pathway functional enrichment, and upstream regulator analysis of differentially expressed genes targeted by differentially expressed miRNAs. Both parasites induced specific miRNA profiles. GO analysis revealed that the in silico predicted targets of the differentially expressed miRNAs regulated different cellular processes involved in development and immunity, and most of the identified KEGG pathways were related to chronic diseases and infection. Considering that the differentially expressed miRNAs identified here modulated crucial host cellular targets that participate in determining the success of infection, these miRNAs might explain the differing congenital transmission rates between the two parasites. Molecules of the different pathways that are regulated by miRNAs and modulated during infection, as well as the miRNAs themselves, may be potential targets for the therapeutic control of either congenital Chagas disease or toxoplasmosis.

Introduction

More than one billion people worldwide are burdened by parasitic diseases (1). Of these, Chagas disease (American trypanosomiasis) and toxoplasmosis are caused by Trypanosoma cruzi (T. cruzi) and Toxoplasma gondii (T. gondii), respectively (2–4). Chagas disease is a devastating but neglected health problem in Latin America. Due to the extensive global migration of asymptomatic individuals, this infection has become an emerging disease in non-endemic countries. Congenital transmission is partially responsible for the progressive globalization of Chagas disease (5, 6). T. gondii is one of the most successful parasites on earth and is estimated to infect over one billion people worldwide (7). Importantly, both parasites can be congenitally transmitted and cause perinatal morbidity and mortality (2–4) but present different transmission rates. T. cruzi has a low transmission rate (1–12%) (6, 8) while T. gondii has a high transmission rate (22–72%) (3). Moreover, both parasites elicit a different local placental immune response that might be related to infection susceptibility (9, 10). Thus, T. cruzi and T. gondi infection is related to the expression and activation of different Toll-like receptors, which in turn mediate the secretion of different cytokines and chemokines in defense against both parasites in the placenta (9, 11).

However, the probability of congenital transmission depends on a variety of complex interactions between the pathogen and the host (4, 12). In particular, parasite factors, placental factors, and maternal and developing fetal immune systems determine infection occurrence (4, 13). In this context, both parasites display sophisticated strategies to avoid host defenses and virulence factors that increase the chance of establishing infection and long-term persistence. One of these strategies is the ability to modulate host cell gene expression (14–16) through small non-coding RNAs such as microRNAs (miRNAs) that repress mRNAs in a sequence-specific manner by either an mRNA degradation process or through mRNA translation inhibition (17–19). MiRNAs play a key role in fine tuning gene expression in multiple physiological and pathological conditions including T. cruzi (20) or T. gondii (21) infection. Interestingly, the largest miRNA cluster in humans is encoded in chromosome 19 (C19MC; 19q13.41) and is almost exclusively expressed in the placenta (22). Both C19MC-derived and non-C19MC-derived miRNAs have been associated with placental development pathologies such as pre-eclampsia and infection (23, 24). However, there is no report in the literature regarding altered miRNA expression in placental tissue in response to either parasite.

Here, we infected human placental explants (HPE) ex vivo by 2 h of incubation with either T. cruzi or T. gondii, then analyzed both miRNA expression profiles by miRNA sequencing and quantitative PCR of selected miRNAs. In addition, we used sequence-based miRNA target prediction and performed pathway functional enrichment and upstream regulator analysis of differentially expressed genes targeted by differentially expressed miRNAs (DEMs).

Materials and Methods

Parasite Culture and Harvesting

For T. cruzi, Y strain (T. cruzi II) trypomastigotes were obtained from previously infected Vero cells (ATCC^® CCL-81) grown in RPMI medium supplemented with 5% fetal bovine serum and 1% antibiotics (penicillin-streptomycin) at 37°C in a humid atmosphere with 5% CO_2. Parasites invaded the cells and replicated intracellularly as amastigotes. After 48–72 h, amastigotes transformed back to trypomastigotes and lysed the host cells. The infective trypomastigotes were separated from cellular debris by low speed centrifugation (500 × g) for 10 min. Parasites were isolated from the supernatant by centrifugation at 3500×g during 15 min, suspended in RPMI media (without fetal bovine serum, 1% (penicillin-streptomycin) (RPMI 1640^®, Biological Industries Ltd.), and quantified in a Neubauer chamber (9).

For T. gondii, semi-confluent HFF cells were infected with RH tachyzoites at a multiplicity of infection of 3 to 5 parasites per cell. After 40 h, the infected cells were washed, then monolayers were scraped from the flasks and passed through 20-, 23-, and 25-gauge needles. Tachyzoites were purified from host cell debris with a 3.0 μm Isopore filter (Merck Millipore^®) (25).

The laboratory has been certificated as a Biosafety level 2 laboratory by the Biosafety Committee (“Unidad de Prevención de Riesgo”) of the “Facultad de Medicina, Universidad de Chile” (approval # 0403/2019).

HPE Infection

Human term placentas were obtained from 3 women with uncomplicated pregnancies with vaginal or caesarean delivery. Informed consent for experimental use of the placenta was given by each patient as stipulated by the Code of Ethics of the “Servicio de Salud Metropolitana Norte” (approval number 0010/2019). Exclusion criteria consisted of the following: major fetal abnormalities, placental tumor, intrauterine infection, obstetric pathology, positive serology for Chagas disease, and any other maternal disease. Donor patients were negative for anti-T. gondii IgG/IgM antibodies. The organs were collected in cold, sterile, saline-buffered solution (PBS) and processed no more than 30 min after delivery. The dissected explants (approximately 50 mg of tissue) were washed with sterile PBS to remove the blood and co-cultivated with T. cruzi trypomastigotes or T. gondii tachyzoites (10⁵ parasites/ml) in serum free RPMI media. After 2 h of co-cultivation, explants were collected in RNAlater™ solution (ThermoFisher Scientific^®), then stored at 4°C for 24 h and at -80°C for posterior RNA isolation (9). Three independent experiments were carried out in triplicates; HPEs from each placenta were infected with either T. cruzi or T. gondii parasites. The parasite load in the HPEs was confirmed by real-time PCR as described previously by us (9, 10, 26).

RNA Extraction

Total RNA was extracted from HPE by mechanical disruption in 1.3 ml of RNA-solv^® reagent (Omega Bio-tek) and isolated using an E.Z.N.A^® total RNA kit I (Omega Bio-tek) according to manufacturer instructions. RNA was stored at -80°C until analysis. The concentration and quality of RNA was determined with a Qubit^® RNA HS Assay kit and an IQ Assay kit (Invitrogen), respectively. Only RNA samples with an IQ ≥8 were further analyzed for quality with an Agilent 2100 Bioanalyzer System (Agilent Technologies, USA) using an RNA Nano 6000 Assay Kit. RNA samples with RNA integrity numbers >5.0 were used for miRNA profiling analysis (9).

Library Construction and Sequencing

Small RNA-Seq libraries were constructed with an Illumina TruSeq Small RNA library preparation kit according to manufacturer protocols. To assess the quality of the libraries, a DNA High Sensitivity Chip was used in an Agilent 2100 Bioanalyzer (Agilent Technologies, USA). The libraries were sequenced on an Illumina NextSeq 500 platform. For each condition, three independent biological replicates were sequenced and paired-end reads were generated.

Data Analysis

Raw read quality was evaluated using the FastQC tool (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/). Raw reads were analyzed with Trim Galore Cutadapt software (27) and low-quality reads were removed (phred value < 30) in order to obtain clean reads. Clean reads with a length range of 18–35 nucleotides were chosen to perform all subsequent analyses. The sofware STAR (28) was used to align all reads to the reference human genome sequence (Hg38). The read counts per coding sequence were determined using HTSeq-count (29). To evaluate replicates, we used Principal Component Analysis, Pearson correlation and standardized median correlation analyses and box plots. The program EdgeR was used for differential expression analysis (30). Differentially expressed genes were defined as genes with p-value <0.05. Target gene prediction performed by using miRDB, psRNA target, and TargetScan sofwares.

Enrichment Analyses

miRNA set enrichment analysis was performed using the TAM 2.0 tool (http://www.lirmed.com/tam2/). KEGG pathways and functional annotation of the predicted target genes (https://www.genome.jp/kegg/kegg2.html) were also analyzed to determine the biological processes, molecular functions, cellular components, and associations with disease.

RT-qPCR

RNA enriched in small RNAs was extracted from HPEs (approximately 50 mg of tissue) by mechanical disruption in 1 ml RNAzol^® RT (Sigma-Aldrich) according to manufacturer instructions and stored at -80°C until analysis. The concentration of the isolated miRNAs was determined using a Qubit^® Quant-iT™ microRNA Assay Kit (Molecular Probes). cDNA of mature miRNAs was synthesized with a MystiCq™ microRNA cDNA Synthesis Mix Kit (Sigma-Aldrich Merck) per manufacturer guidelines. The 25 µl RT-qPCR reaction contained 12.5 µl 2× MystiCq microRNA SYBR Green qPCR Ready Mix, 0.5 µl of 10 µM MystiCq Universal PCR Primer, 0.5 µl of 10 µM of each specific MystiCq microRNA qPCR Assay Primer (Supplementary Table 1), 10.5 µl nuclease-free water, and 1 µl cDNA. All RT-qPCR reactions were performed in triplicates. RT-qPCR was performed under the following cycling conditions: initial denaturation at 95°C for 2 min, followed by 40 cycles of 95°C for 5 s and 60°C for 30 s. Gene expressions were calculated using the 2^−ΔΔCT relative expression method and normalized to snRNA U6 (RNU6-1) expression levels (31).

Results

T. cruzi and T. gondii Change the miRNA Expression Profile in HPE

The effects of T. cruzi and T. gondii on placental tissue were assayed in HPE after a 2 h challenge with 10⁵ parasites/ml. Total miRNA extracted from infected and non-infected control HPE was analyzed by miRNA Seq. Key characteristics of the obtained sequencing data are summarized in Table 1. A total of 680 and 686 DEMs were identified in T. cruzi and T. gondii infected HPE, respectively. Only 14 DEMs with a minimum 1.5-fold change in expression and a 95% probability of being differentially expressed (p ≤ 0.05) were identified in T. cruzi challenged samples (Figure 1A). In T. gondii challenged samples, the number of DEMs increased to 42 (Figure 1B). Comparison of T. cruzi infected HPE with non-infected control samples showed that five miRNAs were downregulated and nine were upregulated. In T. gondii infected tissues, 13 miRNAs were downregulated and 29 were upregulated. The Venn diagram in Figure 1C shows the miRNAs that were differentially expressed in the presence of both parasites compared to non-infected control samples and in HPE infected with either parasite. The complete list of DEMs in response to ex vivoT. cruzi and T. gondii infection is shown in Table 2.

Figure 1

Functional Annotation and KEGG Pathway Enrichment Analysis of miRNA Target Genes

To better understand the roles of the miRNAs identified in HPE in response to ex vivo infection with both parasites, the target genes of the miRNAs were identified using miRDB, psRNA target, and TargetScan. GO and KEGG enrichment analyses used to identify the biological functions of the DEMs (p < 0.05) during ex vivoT. cruzi and T. gondii infection revealed 679 best scored target genes of the 14 miRNAs from T. cruzi vs. control, 1970 best scored target genes of the 42 miRNAs from T. gondii vs. control, and 2011 best scored target genes of the 42 miRNAs from T. cruzi vs. T. gondi. The target genes of the DEMs are shown in Supplementary Table 2. Among the significantly enriched GO terms in T. cruzi vs. control samples, DEMs were significantly enriched in regulation of NFκB pathways, chondrocyte development, cell death including apoptosis, peritoneal cavity homeostasis, angiogenesis, cell cycle, megakaryocyte differentiation, Toll-like receptor signaling pathway, and immune response including innate immunity (Figure 2A). Among the significantly enriched GO terms in T. gondii vs. control samples, DEMs were significantly enriched in cell proliferation, cell migration, osteoblast differentiation, oxidative stress, lipid metabolism, regulation of stem cells including embryonic stem cells, hepatotoxicity, DNA damage response, regulation of NFκB pathways, smooth muscle proliferation, T-helper 17 cell differentiation, T-cell activation, and response to estrogen (Figure 2B).

Figure 2

In addition, we performed GO and KEGG analyses to identify different pathologies in which the T. cruzi- and T. gondii-induced DEMs were related. In T. cruzi vs. control samples, significantly enriched DEMs were related to metabolic syndrome, IgA-nephropathy, acute childhood lymphoblastic leukemia, atherosclerosis, oral lichen planus, human papilloma virus infection, psoriasis, neuropsychiatric disorders, heart diseases, pancreatic carcinoma, Löfgren’s syndrome, Mycobacterium tuberculosis infection, male infertility, and gastric carcinoma (Figure 3A). In T. gondii vs. control samples, the significantly enriched DEMs were related to ankylosing spondylitis, type 2 diabetes mellitus, hypertrophic cardiomyopathy, congenital heart disease, fetal alcohol syndrome, pulmonary hypertension, ulcerative colitis, cystic fibrosis, vascular diseases, human cytomegalovirus infection, muscular dystrophy, liver diseases, coxsackievirus infection, and diabetic retinopathy (Figure 3B).

Figure 3

Specific GO and KEGG enrichment analyses focused on the functions of DEMs in response to both parasites that were related to parasitic diseases and/or placenta pathology. These DEMs are listed in Table 3. In this context, we found that the significantly enriched DEMs were related to the regulation of apoptosis, wound healing, cardiomyocyte apoptosis, heart development, skeletal muscle cell differentiation, Toll-like receptor signaling pathway, innate immunity, epithelial to mesenchymal transition, chromatin remodeling, and nephrotoxicity (Figure 4). Moreover, we analyzed the significantly enriched DEMs related to transcription activation and found that the following transcription factors or its activators were modulated by them: MYOG, calcineurin, AP-1, TNFSF12, NFκB1, myogenin, MYOD, MYF5, MRF4, and TP53 (Figure 5).

Figure 4

Figure 5

Validation of miRNA Expression by RT-qPCR

Validation of miRNAs associated with parasite infection or pregnancy related pathologies was performed by selecting six miRNAs [miR-3074 (26), miR-127-3p (27, 28), miR-30e-3p (29), miR-512-3p (30), miR-515-5p (31), and miR-190b (32)] for confirmation by real time PCR to verify the DEM expression levels. Expression of miR-3074, miR-127-3p, and miR-30e-3p (Figures 6A–C) was analyzed in HPE in response to both parasites. miR-512-3p and miR-515-5p (Figures 6D–E) expression was determined in response to T. cruzi infection and miR-190b expression in response to T. gondii infection (Figure 6F). All selected miRNAs except for miR-30e-3p (Figure 6C) were differentially expressed. Thus, miR-3074 expression (Figure 6A) was significantly decreased (T. cruzi: 57.03 ± 19.99%, p ≤ 0.01; T. gondii: 69.84 ± 24.67%, p ≤ 0.01) with respect to the control but not the infected samples. Decreased miR-3074 expression was expected in the T. gondii infected samples but not in the T. cruzi infected samples. According to the miRNA Seq data, miR-3074 was upregulated in T. cruzi challenged samples. Similar results were observed for miR-127-3p (Figigure 6B). Expression of miR-127-3p was significantly decreased in HPE infected with either parasite (T. cruzi: 68.218 ± 16.41%, p ≤ 0.01; T. gondii: 73.13 ± 22.45%, p ≤ 0.01) compared to the control but not to infected samples; we expected an increase in miR-127-3p expression in the presence of T. gondii since in the miRNA Seq data this particular miRNA was increased (Table 3). RT-qPCR validation results for miR-512-3p, miR-515-5p, and miR-190b confirmed the miRNA Seq data. Thus, miR-512-3p expression increased (40.83 ± 22.53%, p ≤ 0.01) (Figure 6D) and miR-515-5p expression decreased (21.44 ± 8.60%, p ≤ 0.01) (Figure 6E) significantly in T. cruzi infected samples. In T. gondii infected HPE, miR-190b expression was significantly increased (59.02 ± 37.73%, p ≤ 0.01) (Figure 6F).

Figure 6

Discussion

Pathogens have evolved strategies to exploit resources from their hosts to maximize their own survival, replication, and dissemination. Thus, different kinds of pathogens (including parasites) have developed sophisticated mechanisms that include hijacking host cellular machinery to modulate host gene expression to inhibit defense responses (16, 32, 33). Both of the parasites studied here are able to manipulate host gene expression. For instance, during T. cruzi cell and tissue invasion, cell reprogramming affects cellular stress responses, host metabolism, and a significant number of transcription factors (16, 34). T. gondii also reprograms host cells, primarily targeting cell-specific transcription factors that regulate host defenses (i.e., NF-κB, interferon regulatory factor, and JAK/STAT) by regulating their intrinsic activities and expression levels (35). In addition, certain parasites including T. cruzi and T. gondii can alter host miRNA expression to favor both parasite clearance and infection (18, 19). Moreover, different strains of T. gondii can induce specific miRNAs in mice that have been proposed as biomarkers for early infection (19, 36).

Mature miRNAs regulate the expression of over 30% of fundamental genes; these are involved in key biological processes including development, cellular proliferation and differentiation, apoptosis, metabolism, and immune response (18, 19, 37); all of these determine infection success or failure. Moreover, epigenetic and genetic defects in miRNAs and their processing machinery are a common hallmark of infection and diseases that include pregnancy-specific pathologies such as preeclampsia (17, 38).

Most of the transcriptomic studies as well as those analyzing miRNA profiles have focused on a single type of cell response (34) or on tissues or organs in animal models (39, 40); no studies have focused on human tissues. The present study is the first report on the miRNA profile of the human placenta in response to T. cruzi or T. gondii infection. We identified 680 and 686 DEMs, respectively, in T. cruzi and T. gondii infected samples. T. cruzi modulated only 14 DEMs with a minimum of a 1.5-fold change in expression and a 95% probability of being differentially expressed (Figure 1A). In contrast, the number of DEMs increased to 42 in T. gondii challenged HPE (Figure 2B). Our results showed that the DEMs identified here are related to the regulation of different fundamental cellular processes (Figure 2) as well as different pathologies (Figure 3). Moreover, fundamental cellular processes related to placenta pathologies and embryonic development are affected by the presence of both parasites (Figures 4–5). It is important to point out, that tissue response to infection is relevant during disease progression. The presence of the parasites leads to tissue damage as well as immune and regulatory/repair responses, which can lead to fibrosis and tissue dysfunction as observed in chagasic cardiomyopathy (41) or Toxoplasma induced encephalitis in immune-compromised individuals (42).

Three miRNAs, miR-21, miR-146a/b, and miR-210, were overrepresented in most of the ontology terms (Table 2, Supplementary Tables 3–4). Previous studies have implicated these miRNAs in immune and inflammatory response regulation via macrophage polarization controlled through transcription factor regulation in response to signals from the microenvironment (43, 44). Concordantly, in T. cruzi-infected mice, increased miR-21 expression in the heart has been correlated with a parasitemia peak at 30 days post-infection (39). In placenta, miR-21 has been associated with trophoblast differentiation and invasion and miR-21 dysregulation leads to placental pathology (45). MiR-146a is a negative feedback regulator in TLR-4 signaling that acts by repressing TRAF6 to inhibit NFκB transcription factor activation (46, 47). In macrophages, TRAF6 mediates the induction of the pro-inflammatory cytokine IL-12, which is essential to control T. gondii infection (48). TRAF6 activation is also required for vacuole-lysosome fusion, a fundamental step during T. gondii infection (49). Our results showed that in HPEs, T. gondii and T. cruzi infection increased miR-146a expression. Our previous studies showed that both parasites modulate placental immune response differentially through TLRs and NFκB pathways in HPEs (9, 10) Interestingly, the inhibition of these pathways increased the DNA loads of both parasites in HPEs (10). Increased T. gondii infection in placental tissue is also induced by TLR-4 inhibition (9). In addition, increased levels of miR-146a have been reported in the brains of mice infected with T. gondii, moreover, miR-146a ablation affects early parasite burden and improves survival (50). It was previously reported that miR-210 is induced by damage associated molecular patterns (51). In preeclamptic placentas, miR-210 is increased (52); in the present study, miR-210 was increased in HPE infected with T. cruzi but not with T. gondii. Expression of miR-210 can be directly regulated by the specific binding of NF-κB p50 to its putative promoter (53). In this context, it is important to mention that T. cruzi, but not T. gondii, infection of HPE activates both NF-κB signaling pathways (10). Therefore, the increased level of miR-210 might be a placental response to signal transduction pathway activation.

In addition, several identified pathways, important, e.g. for chondrocyte development, megakaryocyte smooth and muscle cell differentiation, hepatotoxicity, and DNA damage response, are neither related to infection or with placental tissues (Figure 3). This can be explained be the fact, that miRNAs target multiple genes, while individual genes are targeted by multiple miRNAs. Moreover, the same miRNA regulates different genes in different tissues and organs (54, 55). Here, we chose to validate six miRNAs that were associated specifically with parasite infection and/or pregnancy related pathologies (Figure 6). Deregulation of miR-30e-3p has been reported in mice that were experimentally infected with T. gondii (40). This miRNA is also related to Chagas cardiomyopathy (39) and is upregulated in placentas with intrauterine growth restriction (56). Nonetheless, miR-30e-3p expression was unaffected by T. cruzi or T. gondii infection in HPE (Figure 6C). Increased miR-3074-5p expression has been described in placentas from recurrent miscarriages (57) and in livers from T. gondii-infected cats (42). However, miR-3074-5p expression was diminished in HPE infected either with T. cruzi or T. gondii (Figure 6A). The differences between our results and the reported data might be explained by differences in the studied organs (heart and liver versus placenta) and the complexity of the above mentioned placental pathologies. MiR-127 is a placenta-specific miRNA codified in the C14MC cluster (58) and its levels are decreased in placenta-related pathologies such as recurrent miscarriage and small-for-gestational age (59, 60); the downregulation of MiR-127 was also detected in babies infected congenitally with either parasite (3, 4). Concordantly, our results showed that HPE infection with either T. cruzi or T. gondii led to the decrease of this miRNA (Figure 6B). Moreover, a decreased expression of miR-127-3p in non-placental tissues has been reported during T. gondii infection in mice and cats (40, 42, 61), but there is no report regarding miR-127-3p expression in response to T. cruzi infection. Both miR-515-5p and miR-512-3p are placenta-specific miRNAs that are codified in the C19MC cluster (62). Decreased miR-515-5p expression is related to fetal growth restriction (63) and preeclampsia (64). Importantly, this miRNA inhibits human trophoblast differentiation by directly repressing the aromatase P450 (CYP19A1), frizzled 5 (FZD5), and glial cells missing 1 transcription factor (GCM1) genes (65). Trophoblast differentiation is part of the trophoblast epithelial turnover and it has been proposed that this mechanism is part of an antiparasitic placental response against T. cruzi infection (11, 13, 66, 67). Therefore, our reported decrease of miR-515-5p expression during ex vivoT. cruzi infection of HPE (Figure 6D) might be at least partially responsible for the parasite-induced trophoblast differentiation. In contrast, miR-512-3p was upregulated in HPE in response to T. cruzi infection (Figure 6E). Interestingly, miR-512-3p confers resistance to vesicular stomatitis virus in non-placental recipient cells (68) and represses the caspase 8 inhibitor c-FLIP (cellular FLICE-like inhibitory protein); it consequently increases caspase 8 activity (69). Caspase 8 regulates trophoblast differentiation and apoptotic cell death and is activated by T. cruzi (66). Therefore, miR-512-3p upregulation might also be a protective placental response to T. cruzi infection, as it is to viral infection. The upregulation of miR-190b in HPE during T. gondii infection was observed in the RNAseq analysis (Figure 1B, Table 3), then validated by qPCR (Figure 6F). Upregulation of miR-190b promotes cell proliferation and migration and reduces cell apoptosis in different types of cancer (70, 71). Parasites modulate apoptotic responses in infected cells to avoid rapid clearance; T. gondii is particularly capable of blocking apoptosis by different mechanisms (33). In neurons, increased miR-190b expression also increases cell viability, suppresses autophagy, and significantly decreases the levels of pro-inflammatory TNF-α, IL-6, and IL-1β cytokines (72). In this context, we have shown that T. gondii, in contrast to T. cruzi, does not induce pro-inflammatory cytokines in HPE (9). Therefore, it is postulated that the lack of pro-inflammatory cytokine secretion in response to T. gondii in HPE might be related to an increase in miR-190b expression and that, together with the modulation of the apoptotic pathway, it could allow parasite persistence and infection in the placenta.

In conclusion, the present study provides a comparative analysis of RNA sequencing-based miRNA profiles in HPE in response to ex vivoT. cruzi and T. gondii infection (Figure 7). Our findings provide new insight into the capacity of both parasites to modulate host gene expression. GO analysis revealed that the predicted targets of the DEMs were different cellular processes involved in development and immunity, and most of the identified KEGG pathways were related to chronic diseases and infection. Considering that the DEMs identified herein modulate crucial host cellular targets that participate in determining the success of infection, these miRNAs might explain the differences in congenital transmission rates. Molecules of the different pathways that are regulated by miRNAs and modulated during infection, as well as the miRNAs themselves, may be potential targets for the therapeutic control of either congenital Chagas disease or toxoplasmosis.

Figure 7

Funding

UK received a grant from the Network of the European Union, Latin America and the Caribbean Countries on Joint Innovation and Research Activities (ERANet-LAC; grant number ERANet17/HLH-0142). UK, JM, and CC received grants from the National Fund for Scientific and Technological Development (FONDECYT; grant numbers 1190341, 1170126, and 3180452, respectively).

References

• 1

  Dotters-KatzSKullerJHeineRP. Parasitic infections in pregnancy. Obstet Gynecol Surv (2011) 66:515–25. doi: 10.1097/OGX.0b013e3182385fde

  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  BłaszkowskaJGóralskaK. Parasites and fungi as a threat for prenatal and postnatal human development. Ann Parasitol (2014) 60:225–34.

  • Google Scholar
  • View reference in article

• 3

  CarlierYTruyensCDeloronPPeyronF. Congenital parasitic infections: a review. Acta Trop (2012) 121:55–70. doi: 10.1016/j.actatropica.2011.10.018

  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  KemmerlingUOsunaASchijmanAGTruyensC. Congenital Transmission of Trypanosoma cruzi: A Review About the Interactions Between the Parasite, the Placenta, the Maternal and the Fetal/Neonatal Immune Responses. Front Microbiol (2019) 10:1854. doi: 10.3389/fmicb.2019.01854

  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  SchmuñisG. Status of and cost of Chagas disease worldwide. Lancet Infect Dis (2013) 13:283–4. doi: 10.1016/S1473-3099(13)70032-X

  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  Pérez-MolinaJAMolinaI. Chagas disease. Lancet (London England) (2018) 391:82–94. doi: 10.1016/S0140-6736(17)31612-4

  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  MeloMBJensenKDCSaeijJPJ. Toxoplasma gondii effectors are master regulators of the inflammatory response. Trends Parasitol (2011) 27:487–95. doi: 10.1016/j.pt.2011.08.001

  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  RendellVRGilmanRHValenciaEGaldos-CardenasGVerasteguiMSanchezLet al. Trypanosoma cruzi-infected pregnant women without vector exposure have higher parasitemia levels: implications for congenital transmission risk. PloS One (2015) 10:e0119527. doi: 10.1371/journal.pone.0119527

  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  CastilloCMuñozLCarrilloILiempiAGallardoCGalantiNet al. Ex vivo infection of human placental chorionic villi explants with Trypanosoma cruzi and Toxoplasma gondii induces different Toll-like receptor expression and cytokine/chemokine profiles. Am J Reprod Immunol (2017) 78:1–8. doi: 10.1111/aji.12660

  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  LiempiACastilloCMedinaLRojasMMayaJDParraguezVHet al. Ex vivo infection of human placental explants with Trypanosoma cruzi and Toxoplasma gondii: Differential activation of NF kappa B signaling pathways. Acta Trop (2019) 199:105153. doi: 10.1016/j.actatropica.2019.105153

  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  CastilloCMuñozLCarrilloILiempiAMedinaLGalantiNet al. Toll-like receptor-2 mediates local innate immune response against Trypanosoma cruzi in ex vivo infected human placental chorionic villi explants. Placenta (2017) 60:40–6. doi: 10.1016/j.placenta.2017.10.005

  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  SenRNayakLDeRK. A review on host–pathogen interactions: classification and prediction. Eur J Clin Microbiol Infect Dis (2016) 35:1581–99. doi: 10.1007/s10096-016-2716-7

  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  LiempiACastilloCCarrilloIMuñozLDroguettDGalantiNet al. A local innate immune response against Trypanosoma cruzi in the human placenta: The epithelial turnover of the trophoblast. Microb Pathog (2016) 99:123–9. doi: 10.1016/j.micpath.2016.08.022

  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  De SouzaWUlisses de CarvalhoTM. Active penetration of Trypanosoma cruzi into host cells: Historical considerations and current concepts. Front Immunol (2013) 4:2. doi: 10.3389/fimmu.2013.00002

  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  HakimiMaCannellaD. Apicomplexan parasites and subversion of the host cell microRNA pathway. Trends Parasitol (2011) 27:481–6. doi: 10.1016/j.pt.2011.07.001

  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  CastilloCCarrilloILibischGJuizNSchijmanAGRobelloCet al. Host-parasite interaction: Changes in human placental gene expression induced by Trypanosoma cruzi. Parasit Vectors (2018) 11:1–13. doi: 10.1186/s13071-018-2988-0

  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  HayderHO’BrienJNadeemUPengC. MicroRNAs: crucial regulators of placental development. Reproduction (2018) 155:R259–71. doi: 10.1530/REP-17-0603

  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  ZhengYCaiXBradleyJE. MicroRNAs in parasites and parasite infection. RNA Biol (2013) 10:377–85. doi: 10.4161/rna.23716

  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  AcuñaSMFloeter-WinterLMMuxelSM. MicroRNAs: Biological Regulators in Pathogen–Host Interactions. Cells (2020) 9:113. doi: 10.3390/cells9010113

  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  FerreiraLRPFradeAFSantosRHBTeixeiraPCBaronMANavarroICet al. MicroRNAs miR-1, miR-133a, miR-133b, miR-208a and miR-208b are dysregulated in Chronic Chagas disease Cardiomyopathy. Int J Cardiol (2014) 175:409–17. doi: 10.1016/j.ijcard.2014.05.019

  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  CaiYShenJ. Modulation of host immune responses to Toxoplasma gondii by microRNAs. Parasite Immunol (2017) 39:1–6. doi: 10.1111/pim.12417

  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  SadovskyYMouilletJFOuyangYBayerACoyneCB. The function of trophomirs and other micrornas in the human placenta. Cold Spring Harb Perspect Med (2015) 5:1–16. doi: 10.1101/cshperspect.a023036

  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  DumontTMFMouilletJFBayerAGardnerCLKlimstraWBWolfDGet al. The expression level of C19MC miRNAs in early pregnancy and in response to viral infection. Placenta (2017) 53:23–9. doi: 10.1016/j.placenta.2017.03.011

  • CrossRef
  • Google Scholar
  • View reference in article

• 24

  Delorme-AxfordESadovskyYCoyneCB. The Placenta as a Barrier to Viral Infections. Annu Rev Virol (2014) 1:133–46. doi: 10.1146/annurev-virology-031413-085524

  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  NardelliSCCheF-YSilmon de MonerriNCXiaoHNievesEMadrid-AlisteCet al. The histone code of Toxoplasma gondii comprises conserved and unique posttranslational modifications. MBio (2013) 4:e00922–13. doi: 10.1128/mBio.00922-13

  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  LiempiACastilloCMedinaLGalantiNMayaJDParraguezVHet al. Comparative ex vivo infection with Trypanosoma cruzi and Toxoplasma gondii of human, canine and ovine placenta: Analysis of tissue damage and infection efficiency. Parasitol Int (2020) 76:102065. doi: 10.1016/j.parint.2020.102065

  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  LindgreenS. AdapterRemoval: Easy cleaning of next-generation sequencing reads. BMC Res Notes (2012) 5:337. doi: 10.1186/1756-0500-5-337

  • CrossRef
  • Google Scholar
  • View reference in article

• 28

  DobinADavisCASchlesingerFDrenkowJZaleskiCJhaSet al. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics (2013) 29:15–21. doi: 10.1093/bioinformatics/bts635

  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  AndersSPylPTHuberW. HTSeq-A Python framework to work with high-throughput sequencing data. Bioinformatics (2015) 31:166–9. doi: 10.1093/bioinformatics/btu638

  • CrossRef
  • Google Scholar
  • View reference in article

• 30

  RobinsonMDMcCarthyDJSmythGK. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics (2009) 26:139–40. doi: 10.1093/bioinformatics/btp616

  • CrossRef
  • Google Scholar
  • View reference in article

• 31

  AndrovicPValihrachLEllingJSjobackRKubistaM. Two-tailed RT-qPCR: a novel method for highly accurate miRNA quantification. Nucleic Acids Res (2017) 45:e144. doi: 10.1093/nar/gkx588

  • CrossRef
  • Google Scholar
  • View reference in article

• 32

  VillaresMBertheletJWeitzmanJB. The clever strategies used by intracellular parasites to hijack host gene expression. Semin Immunopathol (2020) 42:215–26. doi: 10.1007/s00281-020-00779-z

  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  MammariNHalabiMAYaacoubSChlalaHDardéM-LCourtiouxB. Toxoplasma gondii Modulates the Host Cell Responses: An Overview of Apoptosis Pathways. BioMed Res Int (2019), 2019:6152489. doi: 10.1155/2019/6152489

  • CrossRef
  • Google Scholar
  • View reference in article

• 34

  ChiribaoMLLibischGParodi-TaliceARobelloC. Early Trypanosoma cruzi infection reprograms human epithelial cells. BioMed Res Int (2014) 2014:439501. doi: 10.1155/2014/439501

  • CrossRef
  • Google Scholar
  • View reference in article

• 35

  HakimiMAOliasPSibleyLD. Toxoplasma effectors targeting host signaling and transcription. Clin Microbiol Rev (2017) 30:615–44. doi: 10.1128/CMR.00005-17

  • CrossRef
  • Google Scholar
  • View reference in article

• 36

  JiaBChangZWeiXLuHYinJJiangNet al. Plasma microRNAs are promising novel biomarkers for the early detection of Toxoplasma gondii infection. Parasit Vectors (2014) 7:433. doi: 10.1186/1756-3305-7-433

  • CrossRef
  • Google Scholar
  • View reference in article

• 37

  SalehiMSharifiM. Exosomal miRNAs as novel cancer biomarkers: Challenges and opportunities. J Cell Physiol (2018) 233:6370–80. doi: 10.1002/jcp.26481

  • CrossRef
  • Google Scholar
  • View reference in article

• 38

  MouilletJFOuyangYCoyneCBSadovskyY. MicroRNAs in placental health and disease. Am J Obstet Gynecol (2015) 213:S163–72. doi: 10.1016/j.ajog.2015.05.057

  • CrossRef
  • Google Scholar
  • View reference in article

• 39

  NavarroICFerreiraFMNakayaHIBaronMAVilar-PereiraGPereiraIRet al. MicroRNA transcriptome profiling in heart of trypanosoma cruzi-infected mice: Parasitological and cardiological Outcomes. PloS Negl Trop Dis (2015) 9:1–17. doi: 10.1371/journal.pntd.0003828

  • CrossRef
  • Google Scholar
  • View reference in article

• 40

  HuRSHeJJElsheikhaHMZhangFKZouYZhaoGHet al. Differential brain microRNA expression profiles after acute and chronic infection of mice with toxoplasma gondii oocysts. Front Microbiol (2018) 9:2316. doi: 10.3389/fmicb.2018.02316

  • CrossRef
  • Google Scholar
  • View reference in article

• 41

  MachadoFSDutraWOEsperLGollobKJTeixeiraMMFactorSMet al. Current understanding of immunity to Trypanosoma cruzi infection and pathogenesis of Chagas disease. Semin Immunopathol (2012) 34:753–70. doi: 10.1007/s00281-012-0351-7

  • CrossRef
  • Google Scholar
  • View reference in article

• 42

  CongWZhangXXHeJJLiFCElsheikhaHMZhuXQ. Global miRNA expression profiling of domestic cat livers following acute Toxoplasma gondii infection. Oncotarget (2017) 8:25599–611. doi: 10.18632/oncotarget.16108

  • CrossRef
  • Google Scholar
  • View reference in article

• 43

  CătanăCSCalinGABerindan-NeagoeI. Inflamma-miRs in aging and breast cancer: Are they reliable players? Front Med (Lausanne) (2015) 2:85. doi: 10.3389/fmed.2015.00085

  • CrossRef
  • Google Scholar
  • View reference in article

• 44

  LiHJiangTLiMQZhengXLZhaoGJ. Transcriptional regulation of macrophages polarization by microRNAs. Front Immunol (2018) 9:1175. doi: 10.3389/fimmu.2018.01175

  • CrossRef
  • Google Scholar
  • View reference in article

• 45

  ChaiwangyenWOspina-PrietoSPhotiniSMSchleussnerEMarkertURMorales-PrietoDM. Dissimilar microRNA-21 functions and targets in trophoblastic cell lines of different origin. Int J Biochem Cell Biol (2015) 68:187–96. doi: 10.1016/j.biocel.2015.08.018

  • CrossRef
  • Google Scholar
  • View reference in article

• 46

  HouJWangPLinLLiuXMaFAnHet al. MicroRNA-146a Feedback Inhibits RIG-I-Dependent Type I IFN Production in Macrophages by Targeting TRAF6, IRAK1, and IRAK2. J Immunol (2009) 183:2150–8. doi: 10.4049/jimmunol.0900707

  • CrossRef
  • Google Scholar
  • View reference in article

• 47

  TaganovKDBoldinMPChangKJBaltimoreD. NF-κB-dependent induction of microRNA miR-146, an inhibitor targeted to signaling proteins of innate immune responses. Proc Natl Acad Sci U S A (2006) 103:12481–6. doi: 10.1073/pnas.0605298103

  • CrossRef
  • Google Scholar
  • View reference in article

• 48

  MasonNJFioreJKobayashiTMasekKSChoiYHunterCA. TRAF6-dependent mitogen-activated protein kinase activation differentially regulates the production of interleukin-12 by macrophages in response to Toxoplasma gondii. Infect Immun (2004) 72:5662–7. doi: 10.1128/IAI.72.10.5662-5667.2004

  • CrossRef
  • Google Scholar
  • View reference in article

• 49

  SubausteCSAndradeRMWessendarpM. CD40-TRAF6 and autophagy-dependent anti-microbial activity in macrophages. Autophagy (2007) 3:245–8. doi: 10.4161/auto.3717

  • CrossRef
  • Google Scholar
  • View reference in article

• 50

  CannellaDBrenier-PinchartMPBraunLvanRooyenJMBougdourABastienOet al. MiR-146a and miR-155 delineate a microRNA fingerprint associated with toxoplasma persistence in the host brain. Cell Rep (2014) 6:928–37. doi: 10.1016/j.celrep.2014.02.002

  • CrossRef
  • Google Scholar
  • View reference in article

• 51

  UnluSTangSWangEMartinezITangDBianchiMEet al. Damage associated molecular pattern molecule-induced micrornas (DAMPmiRs) in human peripheral blood mononuclear cells. PloS One (2012) 7(6):e38899. doi: 10.1371/journal.pone.0038899

  • CrossRef
  • Google Scholar
  • View reference in article

• 52

  ZhuXmHanTSargentILYinGYaoYq. Differential expression profile of microRNAs in human placentas from preeclamptic pregnancies vs normal pregnancies. Am J Obstet Gynecol (2009) 7(6):e38899. doi: 10.1016/j.ajog.2008.12.045

  • CrossRef
  • Google Scholar
  • View reference in article

• 53

  ZhangYFeiMXueGZhouQJiaYLiLet al. Elevated levels of hypoxia-inducible microRNA-210 in pre-eclampsia: New insights into molecular mechanisms for the disease. J Cell Mol Med (2012) 16:249–59. doi: 10.1111/j.1582-4934.2011.01291.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 54

  PeterME. Targeting of mRNAs by multiple miRNAs: The next step. Oncogene (2010) 29:2161–4. doi: 10.1038/onc.2010.59

  • CrossRef
  • Google Scholar
  • View reference in article

• 55

  CaiYYuXHuSYuJ. A Brief Review on the Mechanisms of miRNA Regulation. Genomics Proteomics Bioinformatics (2009) 7:147–54. doi: 10.1016/S1672-0229(08)60044-3

  • CrossRef
  • Google Scholar
  • View reference in article

• 56

  WenHChenLHeJLinJ. MicroRNA expression profiles and networks in placentas complicated with selective intrauterine growth restriction. Mol Med Rep (2017) 16:6650–73. doi: 10.3892/mmr.2017.7462

  • CrossRef
  • Google Scholar
  • View reference in article

• 57

  GuYZhangXYangQWangJHeYSunZet al. Aberrant Placental Villus Expression of miR-486-3p and miR-3074-5p in Recurrent Miscarriage Patients and Uterine Expression of These MicroRNAs during Early Pregnancy in Mice. Gynecol Obstet Invest (2016) 81:112–7. doi: 10.1159/000435879

  • CrossRef
  • Google Scholar
  • View reference in article

• 58

  WilliamsZBen-DovIZEliasRMihailovicABrownMRosenwaksZet al. Comprehensive profiling of circulating microRNA via small RNA sequencing of cDNA libraries reveals biomarker potential and limitations. Proc Natl Acad Sci U S A (2013) 110:4255–60. doi: 10.1073/pnas.1214046110

  • CrossRef
  • Google Scholar

• 59

  YangQGuWWGuYYanNNMaoYYZhenXXet al. Association of the peripheral blood levels of circulating microRNAs with both recurrent miscarriage and the outcomes of embryo transfer in an in vitro fertilization process. J Transl Med (2018) 16(1):186. doi: 10.1186/s12967-018-1556-x

  • CrossRef
  • Google Scholar

• 60

  RodosthenousRSBurrisHHSandersAPJustACDereixAESvenssonKet al. Second trimester extracellular microRNAs in maternal blood and fetal growth: An exploratory study. Epigenetics (2017) 12:804–10. doi: 10.1080/15592294.2017.1358345

  • CrossRef
  • Google Scholar

• 61

  HeJJMaJWangJLXuMJZhuXQ. Analysis of miRNA expression profiling in mouse spleen affected by acute Toxoplasma gondii infection. Infect Genet Evol (2016) 37:137–42. doi: 10.1016/j.meegid.2015.11.005

  • CrossRef
  • Google Scholar

• 62

  DonkerRBMouilletJFChuTHubelCAStolzDBMorelliAEet al. The expression profile of C19MC microRNAs in primary human trophoblast cells and exosomes. Mol Hum Reprod (2012) 18:417–24. doi: 10.1093/molehr/gas013

  • CrossRef
  • Google Scholar

• 63

  HigashijimaAMiuraKMishimaHKinoshitaAJoOAbeSet al. Characterization of placenta-specific microRNAs in fetal growth restriction pregnancy. Prenat Diagn (2013) 33:214–22. doi: 10.1002/pd.4045

  • CrossRef
  • Google Scholar

• 64

  HromadnikovaIKotlabovaKOndrackovaMPirkovaPKestlerovaANovotnaVet al. Expression profile of C19MC microRNAs in placental tissue in pregnancy-related complications. DNA Cell Biol (2015) 34:437–57. doi: 10.1089/dna.2014.2687

  • CrossRef
  • Google Scholar

• 65

  ZhangMMuralimanoharanSWortmanACMendelsonCR. Primate-specific miR-515 family members inhibit key genes in human trophoblast differentiation and are upregulated in preeclampsia. Proc Natl Acad Sci (2016) 113:E7069–76. doi: 10.1073/pnas.1607849113

  • CrossRef
  • Google Scholar

• 66

  CarrilloIDroguettDCastilloCLiempiAMuñozLMayaJDJDJDet al. Caspase-8 activity is part of the BeWo trophoblast cell defense mechanisms against Trypanosoma cruzi infection. Exp Parasitol (2016) 168:9–15. doi: 10.1016/j.exppara.2016.06.008

  • CrossRef
  • Google Scholar

• 67

  LiempiACastilloCDuasoJDroguettDSandovalABarahonaKet al. Trypanosoma cruzi induces trophoblast differentiation: A potential local antiparasitic mechanism of the human placenta? Placenta (2014) 35:1035–42. doi: 10.1016/j.placenta.2014.09.017

  • CrossRef
  • Google Scholar

• 68

  Delorme-AxfordEDonkerRBMouilletJ-FChuTBayerAOuyangYet al. Human placental trophoblasts confer viral resistance to recipient cells. Proc Natl Acad Sci (2013) 110:12048–53. doi: 10.1073/pnas.1304718110

  • CrossRef
  • Google Scholar

• 69

  ChenFZhuH-HZhouL-FWuS-SWangJChenZ. Inhibition of c-FLIP expression by miR-512-3p contributes to taxol-induced apoptosis in hepatocellular carcinoma cells. Oncol Rep (2010) 23:1457–62. doi: 10.3892/or_00000784

  • CrossRef
  • Google Scholar

• 70

  DaiWHeJZhengLBiMHuFChenMet al. miR-148b-3p, miR-190b, and miR-429 regulate cell progression and act as potential biomarkers for breast cancer. J Breast Cancer (2019) 22:219–36. doi: 10.4048/jbc.2019.22.e19

  • CrossRef
  • Google Scholar

• 71

  ZhaoQLiuCCuiQLuanXWangQZhouC. miR−190b promotes colorectal cancer progression through targeting forkhead box protein P2. Exp Ther Med (2019) 19:79. doi: 10.3892/etm.2019.8175

  • CrossRef
  • Google Scholar

• 72

  PeiXLiYZhuLZhouZ. Astrocyte-derived exosomes transfer miR-190b to inhibit oxygen and glucose deprivation-induced autophagy and neuronal apoptosis. Cell Cycle (2020) 19:906–17. doi: 10.1080/15384101.2020.1731649

  • CrossRef
  • Google Scholar