- Open Access
Host-parasite interaction: changes in human placental gene expression induced by Trypanosoma cruzi
Parasites & Vectorsvolume 11, Article number: 479 (2018)
Chagas disease is caused by Trypanosoma cruzi, a parasite endemic to Latin America. Most infections occur in children by vector or congenital transmission. Trypanosoma cruzi establishes a complexity of specific molecular parasite-host cell interactions to invade the host. However, most studies have been mainly focused on the interaction between the parasite and different cell types, but not on the infection and invasion on a tissue level. During congenital transmission, T. cruzi must cross the placental barrier, composed of epithelial and connective tissues, in order to infect the developing fetus. Here we aimed to study the global changes of transcriptome in the placental tissue after a T. cruzi challenge.
Strong changes in gene expression profiling were found in the different experimental conditions, involving the reprogramming of gene expression in genes involved in the innate immune response.
Trypanosoma cruzi induces strong changes in genes involved in a wide range of pathways, especially those involved in immune response against infections.
Chagas disease is a zoonotic disease caused by Trypanosoma cruzi, a parasite endemic to Latin America. Most infections occur in children by vector or congenital transmission. The prevalence of Chagas disease in pregnant women in Latin America ranges between 5–40% depending on the geographical area and the rate of congenital transmission is estimated to be 1–12% . In addition, due to population mobility, Chagas disease has been increasingly detected in other non-endemic countries and continents (where the vector does not exist) such as the USA, Canada, Australia, Europe and Asia [2, 3]. Congenital transmission, in spite of its low transmission rates, is partially responsible for the progressive globalization of the disease [2, 4, 5]. Importantly, congenital infection is responsible for an estimated 22% of new infections in 2010, making this form of transmission epidemiologically relevant .
The parasite presents a complex life-cycle that occurs in both vertebrate and invertebrate hosts, where three major developmental stages are observed: epimastigotes, trypomastigotes and amastigotes. Trypomastigotes constitute the extracellular infective form in mammals where they are able to infect a wide range of nucleated mammalian cells . Interestingly, T. cruzi has co-evolved with mammals to establish a complexity of specific molecular parasite-host cell interactions to invade host cells and tissues, to evade the host immune system and to undergo intracellular replication . Key steps in parasite infection include its host cell penetration and replication of the protozoa in the cytoplasm of infected cells. The application of oligonucleotide and cDNA microarray technologies in the study of host-parasite interactions have permitted rapid and unbiased examination of changes in expression of a large number of genes at the level of transcription [9, 10]. However, these studies have been mainly focused on the interaction between the parasite and different cell types, but not on the infection and invasion on a tissue level.
During congenital transmission, T. cruzi must cross the placental barrier in order to infect the developing fetus [3, 11]. This anatomical barrier is formed by the trophoblast, a two-layer epithelium which is in direct contact with maternal blood, the fetal connective tissue (villous stroma), the endothelium of fetal vessels and the basal laminae that support the epithelia [3, 12]. Interestingly, the congenital transmission rate for T. cruzi is low [4, 13] and it has been proposed that the placenta might play an important role avoiding parasite infection .
The study of gene expression profiles during infection constitutes a very powerful tool to analyze global responses of several kinds of cells and tissues, allowing the identification of new genes and/or pathways implicated in the establishment of the infection and pathogenesis as well as possible local tissue responses [3, 10]. Therefore, here we aimed to study the global changes of transcriptome in the placental tissue after T. cruzi challenge.
Trypomastigotes from T. cruzi (Y Strain, T. cruzi II) were obtained from previously infected Vero cells (ATCC® CCL-81) grown in RPMI medium supplemented with 5% fetal bovine serum (FBS) and antibiotics (penicillin-streptomycin) at 37 °C in a humid atmosphere at 5% CO2. Parasites invaded the cells and replicated intracellularly as amastigotes, after 48–72 h; amastigotes transformed back to trypomastigotes and lysed host cells. The infective trypomastigotes were separated from cellular debris by low speed centrifugation (500× g). From the supernatant, the parasites were isolated by centrifugation at 3500× g, suspended in RPMI media (without FBS, 1% antibiotics; RPMI 1640, Biological Industries Ltd., Kibbutz Beit Haemek, Israel) and quantified in a Neubauer Chamber [14, 15].
HPE were obtained from healthy mothers with uncomplicated pregnancies by cesarean delivery. Placentas were processed in a class II laminar flow hood immediately after delivery. The maternal and fetal surfaces were discarded and villous tissue was obtained from the central part of the cotyledons. The dissected explants were washed with sterile PBS in order to get rid of the blood and co-cultivated with T. cruzi trypomastigotes in serum free RPMI media. HPE were challenged with 105 or 106 parasites/ml, since these concentrations have been proposed to correlate with low or high parasitaemia, respectively . For validation experiments, LPS (10 ng/ml) was used as positive control. After 2 or 24 h of infection (in order to study early and late placental responses [16,17,18], explants were collected in RNA later solution (Thermo Fisher Scientific, Waltham, Massachusetts, USA), stored at 4 °C for 24 h and at -80 °C for posterior RNA isolation .
RNA purification and microarray experiment
Total RNA was isolated with a Purelink RNA isolation kit (Thermo Fisher Scientific) according to the manufacturer’s instructions. RNA integrity was analyzed with a Bioanalyzer 2100 (Agilent Technologies, Santa Clara, California, USA) obtaining RNA integrity numbers (RIN) above 8 for all samples (on a scale based on an rRNA 28S/18S ratio where a RIN of 1 corresponds to a totally degraded RNA and 10 to a totally non-degraded RNA). RNA concentration was quantified by spectrophotometry (Nanodrop, Thermo Fisher Scientific). One hundred nanograms of total RNA was reverse-transcribed into cDNA, then transcribed to cRNA and Cy3-labeled with a Low Input Quick Amp-One Color Labeling Kit (Agilent Technologies). The labeled cRNA was purified with an illustra RNAspin Mini Isolation Kit (GE Healthcare, Little Chalfont, UK) and the total yield was measured with a Qubit RNA HS Kit (Thermo Fisher Scientific). Hybridization, washing, assembling of the chips, and scanning were performed according to the manufacturers’ instructions. Briefly, labeled samples were hybridized with SurePrint G3 Human GE 8x60K chips for 17 h at 60 °C in an Agilent hybridization oven at 10× rpm. Posterior washing, stabilization and drying procedures were performed according to Agilent’s Low Input Quick Amp Labeling Kit instructions .
Chips were scanned with an Agilent microarray scanner G2565BA; the software Agilent Feature Extraction (version 9.5.1), was used for quality control, data filtering and data normalization. Extracted data from the SurePrintG3 8x60K chips were analyzed using GeneSpring GX 13.0 software. Genes showing a 2-fold change in their expression (or more) with P ≤ 0.05 were considered differentially expressed using ANOVA and Benjamini-Hochberg false discovery rate correction for multiple testing. Analysis of interaction networks between upregulated genes in each experimental group was performed with Cytoscape network visualization and integration software and GeneMania open-source gene function prediction service plug-in (http://www.genemania.org/) and visualized by the corresponding Cytoscape software version 3.0.2 plug-ins  The weighting of the network attributes was set to Gene Ontology (GO)-based weighting for biological processes. Gene set enrichment analysis (GSEA) was performed with GSEA 3.0 software (Broad Institute, Cambridge, Massachussets, USA) . Each gene set permutation was performed 1000 times, analyses were based on the GO pathways database (http://geneontology.org/)  and a normalized enrichment score (NES) was obtained for each gene set. An enrichment map from data obtained with GSEA was generated with the Enrichment Map plug-in for Cytoscape 3.0 . NES and false discovery rate (FDR) q-value were considered as parameters for the analysis. The NES value allows the comparison of analysis results across gene sets because it considers differences in gene sets sizes corrected by the size of the expression dataset; the FDR represents the estimated probability that a gene set with a given NES represents a false positive finding .
One hundred picograms of total RNA was retro-transcribed to cDNA using an M-MLV Reverse Transcriptase system with Oligo(dT) primers (Thermo Fisher Scientific). For real-time reactions, 10 μl of Sensifast qPCR Master Mix (Bioline, London, UK) was mixed with 100 mM forward and reverse primers and 2 μl of cDNA. Samples were analyzed in an ABI 7300 real-time PCR system (Applied Biosystems) with an initial step of 3 minutes at 95°C for polymerase activation, followed by 40 cycles at 95°C for 5 seconds for denaturation, and 60°C for 15 seconds for annealing/extension.333. Results were analyzed against human GAPDH as a housekeeping gene and expressed using the ΔΔCT method (Pfaffl, 2001). The sequences of primers can be found in Table 1. The results were expressed as the mean ± SD. The significance of differences was evaluated using ANOVA followed by Dunnett’s post-hoc test as indicated.
Trypanosoma cruzi changes the gene expression profile in HPE
The effect of the parasite on placental tissue was assayed in HPE after challenges with a low (105 parasites/ml) or a high (106 parasites/ml) concentration of trypomastigotes for 2 or 24 h.
Total RNA extracted from infected and non-infected control HPE, was labeled and hybridized to a Human GE 60K Microarray, which allows the evaluation of the gene expression profile of 26,083 different human genes. Genes showing at least a 2-fold change in their expression and a 95% probability of being differentially expressed (P ≤ 0.05) were significantly regulated during parasite challenge. Figure 1 shows the total number of significant differentially expressed genes between infected and non-infected control HPE. A low parasite concentration induces the downregulation of 431 and 1474 genes as well as the upregulation of 210 and 469 genes after 2 and 24 h of parasite challenge, respectively (Fig. 1a). After a high parasite concentration challenge, 157 and 722 genes were downregulated, and 342 and 454 were upregulated after 2 and 24 h, respectively (Fig. 1b). Major changes occurred after 24 h of parasite challenge with the lowest parasite concentration. A selection of the most upregulated and downregulated genes (fold change range between 34.70 and 71.43) is shown in Table 2. Among the most upregulated genes are those involved in immune response such as CXCL9, TLR-7, TLR-8, CD46, C1qTNF3, HLA-DQB1 and CCL20; genes involved in extracellular matrix (ECM) remodeling (ADAM12, ADAMTSL3, MMP10) and related to pregnancy. The list of most downregulated genes also includes genes related to immunity such as LBP, CD14, DCD and IL-6.
The Venn diagrams in Fig. 2a show that 19 genes are upregulated in the four different experimental conditions compared to control non-infected samples, which are also shown in the corresponding heatmap (Fig. 2b) listed in Table 3. Contrarily, only 5 genes are downregulated in the same conditions (Fig. 3a, b), which are listed in Table 3. Most of the upregulated genes are related to pregnancy processes.
Trypanosoma cruzi alters a wide range of biological processes in HPE
GO and pathway analysis were performed using GeneSpringGX 13.0 software (Agilent Technologies), comparing the different experimental conditions described above. Our results indicate that a wide range of biological processes are altered at the different conditions in presence of the parasite (Table 4). The different biological processes detected include immune response, pregnancy related processes and signaling. In order to understand the relationships between those biological processes, we performed a GSEA analysis of gene sets at different times and parasite load challenges. We analyzed biological processes pathways based on gene ontology results (Fig. 4). The biggest cluster is composed of pathways related with immune response, followed by development morphogenesis cluster, regulation of metabolic processes and signal transduction genes, metabolic processes, homeostasis, response to stimulus, cell death and endocytosis (Fig. 4).
Several gene sets grouped within the main cluster (immune response), are enriched at 2 or 24 h post-infection after parasite challenges of 105 trypomastigotes/ml but not of 106 parasites. Thus, immune system process pathways are positively regulated against 105 parasites at 2 h (NES: 1.57; FDR q-value = 0.0013) and 24 h (NES: 2.75; FDR q-value = 0.0018), downregulated with 106 trypomastigotes/ml at 2 h (NES: -2.1; FDR q-value = 0.0087) but not at 24 h (NES: 1.44; FDR q-value = 0.18). Regulation of immune response is positively regulated with 105 parasites at 2 h (NES: 1,50; FDR q-value = 0.022) and 24 h (NES: 1.98; FDR q-value = 0.01) but negatively with 106 parasites at 2 h (NES: -1,77; FDR q-value = 0.067) and without significant changes at 24 h (NES: -0,83; FDR q-value not significant).
Gene Ontology process regulation of cytokine production is upregulated 2 h post-infection (NES: 1.36; FDR q-value = 0.023) but not in the other experimental groups. The inflammatory response pathway is upregulated 24 h post-infection against 105 trypomastigotes (NES: 1.65; FDR q-value = 0.049) but downregulated with 106 parasites 2 h post-infection (NES: 1.99; FDR q-value = 0.023). Amongst the changes of pathways related to other processes, it can be highlighted the upregulation of cell proliferation pathway 24 h post-infection (NES: 2.18; FDR q-value = 0.003) and cellular response to hormone stimulus (2 h 105 Trypos NES: 2.53; FDR q-value = 0.005).
To understand the nature of the interaction between upregulated or downregulated genes, we performed a gene interaction analysis using the GeneMANIA plug-in for Cytoscape 3.0 software [24, 25]. For each experimental condition, we analyzed co-expression, co-localization, physical interactions, genetic interactions shared protein domains and pathways amongst all the differentially expressed genes (Fc ≥ 2) in both upregulated or downregulated gene lists. The relative weight of each process between all interacting genes is depicted in Table 5. In all experimental conditions, co-expression (a category where two genes have similar expression levels) is the predominating interaction, co-localization (genes expressed in the same tissue) the second in the upregulated conditions; however, in downregulated groups physical interactions (when two gene products are found to interact in protein-protein interaction studies) is the second most common interaction. Shared protein domains, genetic and pathways interaction represent marginal interactions in all groups. A circular layout of the interaction network for each is shown in Additional file 1: Figure S1.
Trypanosoma cruzi induces differential expression of pathogen pattern recognition receptors in HPE
Based on the results of the microarray analysis showing an activation of local immune processes, we decided to validate by RT-qPCR PRR activation and cytokine production. Explants were co-incubated with 105 T. cruzi trypomastigotes for 2 h and the expression of Toll-like receptors (TLRs) and NOD-like receptors (NLRs) was assayed. Trypanosoma cruzi trypomastigotes induce statistically significant increases of TLR-2 (96.13 ± 61.6%; F(2, 9) = 5.409; P ≤ 0.01, Fig. 5a), TLR-4 (47.56 ± 27.99%; F(2, 45) = 2.173; P ≤ 0.0001, Fig. 5b), TLR-7 (57.29 ± 24.59%; F(2, 2) = 6.048; P ≤ 0.05, Fig. 5c) and TLR-9 (61.56 ± 5.11%; F(2, 6) = 1.630; P ≤ 0.01, Fig. 5d) expression compared to non-infected HPE. However, T. cruzi does not increase significantly NOD-1 and NOD-2 receptors (Fig. 5d, e).
Trypanosoma cruzi increases pro-inflammatory and immune-modulating cytokines in HPE
HPE were co-incubated with 105 T. cruzi trypomastigotes for 2 h as well as in the presence and absence of LPS as positive controls. T. cruzi induces significant increases of the pro-inflammatory cytokines IL-1β (5542.55 ± 1090.11%; F(2, 19) = 64.91; P ≤ 0.0001, Fig. 6a), IL-6 (94.70 ± 40.38%; F(2, 15) = 8.849; P ≤ 0.01, Fig. 6b), IL-12α (77.08 ± 33.01%; F(2, 19) = 64.91; F(2, 15) = 16.13; P ≤ 0.01, Fig. 6c), IFNγ (329.29 ± 162.22%; F(2, 17) = 7.729; P ≤ 0.05, Fig. 6d) and of the immune-modulating cytokines IL-10 (303.34 ± 104.28%; F(2, 20) = 14.87; P ≤ 0.001, Fig. 6e) and TGF β (329.29 ± 162.22%; F(2, 11) = 1.684; P ≤ 0.05, Fig. 6f) but not of IL-17 (Fig. 6g).
The interaction between the host and pathogens, including T. cruzi, is the most important factor in determining whether an infection is successful. Host-parasite interaction includes invasion of the host through primary barriers (such as the placental barrier), evasion of host defenses, pathogen replication in the host, and immunological capacity of the host to control or eliminate the pathogen . Importantly, infected organisms are capable of sensing the intrusion by pathogens and react by triggering host defenses [17, 26]. On the other hand, the parasite is equipped with multiple tools to establish a long-term relationship with the infected host. Tissue infection in particular is relevant during disease progression. The presence of the parasite provokes tissue damage as well as immune and reparatory responses, which can lead to fibrosis and tissue dysfunction as observed in chagasic cardiomyopathy . Considering the temporary existence of the placenta, the effect on parasite infection on this particular tissue is relevant to understand the physiopathology of congenital transmission in order to obtain tools for diagnosis, prognosis and treatment of the disease.
Previous studies about transcriptomics related to T. cruzi and Chagas disease have been focused on a single type cell response [10, 28] or on tissues or organs in animal models [24, 29] but not on human tissues. Here, we describe for the first time, the transcriptomics of an ex vivo human placental tissue model in response to challenges with the parasite.
As expected, T. cruzi modifies an ample range of biological processes during tissue invasion and infection. As described before, the parasite dramatically changes the gene expression in single cells . However, in tissue and organ samples a more complex change in gene expression can be expected since they are composed of different cell types or tissues. For instance, in HPE epithelial cells derived from the trophoblast and fetal capillaries as well as fibroblasts and macrophages in the fetal connective tissue can be found, between others . In addition, ECM components, that are synthetized by the resident cells are also present in tissue and organ samples [18, 19]. Similar results have been obtained in animal models, where important changes in murine myocardium metabolic pathways  and in murine placental response  are described. The increase of gene expression of proteases involved in ECM-remodeling (Table 2) is in concordance with our previous results showing that the parasite increased expression and activity of matrix metalloproteases (MMP-2 and MMP-9) in human placenta . The profound changes in genes involved in signaling agrees with numerous previous studies that showed that T. cruzi activates or inhibits several signal transduction pathways [14, 28, 30].
The effect of T. cruzi on the immune system responses are particularly relevant since they are our main defense against the pathogen. The parasite dramatically changes host genes involved in the immune response. The expression of an important number of genes of innate immunity is increased. Thus, genes related to complement regulation and function such as CD46 and C1q are upregulated. We have previously shown, that during ex vivo infection of HPE, T. cruzi calreticulin (TcCRT) acts as a virulence factor since it binds maternal classical complement component C1q and increases parasite infectivity . On the other hand, TLRs are also increased, particularly TLR-7 and TLR-8 which are increased over 20-fold. The validation experiments show that both mentioned TLRs as well as TLR-2, TLR-4 and TLR-9 are significantly increased (Fig. 5). However, TLR-2, the TLR whose inhibition increases parasite infection in HPE as well as parasite-induced tissue damage  showed no significant increase in the microarray analysis (Additional file 1: Figure S1).
A similar contradictory result was obtained with IL-6; a high parasite concentration decreases the expression of this cytokine more than 60-fold (Table 2). However, a low parasite concentration does not change IL-6 expression (Additional file 1: Figure S1). However, our RT-qPCR data show a significant increase of IL-6 (Fig. 5) that is in concordance with the increase of IL-6 protein in the culture media of HPE after parasite challenge in the same condition , suggesting regulation at post-transcriptional levels.
Another important group of genes that change their expression are those related to pregnancy. Most of the 19 upregulated genes in the four different experimental conditions compared to control samples are related to pregnancy processes (Table 3) and to the maintenance and development of the fetus such as pregnancy specific beta-1-glycoproteins, GH2 (growth hormone 2), CSH1 and CSH2 (chorionic somatomammotropin hormone 1 and 2). Given that the placenta is the sole interface between mother and fetus and that this organ not only protects the fetus from infection but also regulates important metabolic and other physiological processes , it appears easily explainable that different pregnancy related processes are affected.
Trypanosoma cruzi induces strong changes in genes involved in a wide range of pathways, especially those involved in immune response against infections.
Moscatelli G, Moroni S, Garcia-Bournissen F, Ballering G, Bisio M, Freilij H, et al. Prevention of congenital Chagas through treatment of girls and women of childbearing age. Mem Inst Oswaldo Cruz. 2015;110:507–9.
Schmunis G. Epidemiology of Chagas disease in non-endemic countries: the role of international migration. Mem Inst Oswaldo Cruz. 2007;102:75–85.
Liempi A, Castillo C, Carrillo I, Munoz L, Droguett D, Galanti N, et 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.
Perez-Molina JA, Perez AM, Norman FF, Monge-Maillo B, Lopez-Velez R. Old and new challenges in Chagas disease. Lancet Infect Dis. 2015;15:1347–56.
Castillo C, Munoz L, Carrillo I, Liempi A, Medina L, Galanti N, et 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.
Kaplinski M, Jois M, Galdos-Cardenas G, Rendell VR, Shah V, Do RQ, et al. Sustained domestic vector exposure is associated with increased Chagas cardiomyopathy risk but decreased parasitemia and congenital transmission risk among young women in Bolivia. Clin Infect Dis. 2015;61:918–26.
Barrias ES, de Carvalho TM, De Souza W. Trypanosoma cruzi: entry into mammalian host cells and parasitophorous vacuole formation. Front Immunol. 2013;4:186.
Walker DM, Oghumu S, Gupta G, McGwire BS, Drew ME, Satoskar AR. Mechanisms of cellular invasion by intracellular parasites. Cell Mol Life Sci. 2014;71:1245–63.
Shigihara T, Hashimoto M, Shindo N, Aoki T. Transcriptome profile of Trypanosoma cruzi-infected cells: simultaneous up- and down-regulation of proliferation inhibitors and promoters. Parasitol Res. 2008;102:715–22.
Chiribao ML, Libisch G, Parodi-Talice A, Robello C. Early Trypanosoma cruzi infection reprograms human epithelial cells. Biomed Res Int. 2014;2014:439501.
Carlier Y, Truyens C, Deloron P, Peyron F. Congenital parasitic infections: a review. Acta Trop. 2012;121:55–70.
Arora N, Sadovsky Y, Dermody TS, Coyne CB. Microbial vertical transmission during human pregnancy. Cell Host Microbe. 2017;21:561–7.
Rendell VR, Gilman RH, Valencia E, Galdos-Cardenas G, Verastegui M, Sanchez L, et al. Trypanosoma cruzi-infected pregnant women without vector exposure have higher parasitemia levels: implications for congenital transmission risk. PLoS One. 2015;10:e0119527.
Castillo C, Villarroel A, Duaso J, Galanti N, Cabrera G, Maya JD, et al. Phospholipase C gamma and ERK1/2 mitogen activated kinase pathways are differentially modulated by Trypanosoma cruzi during tissue invasion in human placenta. Exp Parasitol. 2013;133:12–7.
Liempi A, Castillo C, Duaso J, Droguett D, Sandoval A, Barahona K, et al. Trypanosoma cruzi induces trophoblast differentiation: a potential local antiparasitic mechanism of the human placenta? Placenta. 2014;35:1035–42.
Fretes RE, Kemmerling U. Mechanism of Trypanosoma cruzi placenta invasion and infection: the use of human chorionic villi explants. J Trop Med. 2012;2012:614820.
Castillo C, Munoz L, Carrillo I, Liempi A, Gallardo C, Galanti N, et 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. https://doi.org/10.1111/aji.12660.
Duaso J, Rojo G, Cabrera G, Galanti N, Bosco C, Maya JD, et al. Trypanosoma cruzi induces tissue disorganization and destruction of chorionic villi in an ex vivo infection model of human placenta. Placenta. 2010;31:705–11.
Castillo C, Lopez-Munoz R, Duaso J, Galanti N, Jaña F, Ferreira J, et al. Role of matrix metalloproteinases 2 and 9 in ex vivo Trypanosoma cruzi infection of human placental chorionic villi. Placenta. 2012;33:991–7.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005;102:15545–50.
Juiz NA, Torrejon I, Burgos M, Fernanda Torres AM, Duffy T, Cayo NM, et al. Alterations in placental gene expression of pregnant women with chronic Chagas disease. Am J Pathol. 2018;188:1345–53.
Merico D, Isserlin R, Stueker O, Emili A, Bader GD. Enrichment map: a network-based method for gene-set enrichment visualization and interpretation. PLoS One. 2010;5:e13984.
Juiz NA, Solana ME, Acevedo GR, Benatar AF, Ramirez JC, da Costa PA, et al. Different genotypes of Trypanosoma cruzi produce distinctive placental environment genetic response in chronic experimental infection. PLoS Negl Trop Dis. 2017;11:e0005436.
Montojo J, Zuberi K, Rodriguez H, Bader GD, Morris Q. GeneMANIA: Fast gene network construction and function prediction for Cytoscape. F1000Res. 2014;3:153.
Sen R, Nayak L, De RK. A review on host-pathogen interactions: classification and prediction. Eur J Clin Microbiol Infect Dis. 2016;35:1581–99.
Nagajyothi F, Machado FS, Burleigh BA, Jelicks LA, Scherer PE, Mukherjee S, et al. Mechanisms of Trypanosoma cruzi persistence in Chagas disease. Cell Microbiol. 2012;14:634–43.
Li Y, Shah-Simpson S, Okrah K, Belew AT, Choi J, Caradonna KL, et al. Transcriptome remodeling in Trypanosoma cruzi and human cells during intracellular infection. PLoS Pathog. 2016;12:e1005511.
Girones N, Carbajosa S, Guerrero NA, Poveda C, Chillon-Marinas C, Fresno M. Global metabolomic profiling of acute myocarditis caused by Trypanosoma cruzi infection. PLoS Negl Trop Dis. 2014;8:e3337.
Maeda FY, Cortez C, Yoshida N. Cell signaling during Trypanosoma cruzi invasion. Front Immunol. 2012;3:361.
Castillo C, Ramirez G, Valck C, Aguilar L, Maldonado I, Rosas C, et al. The interaction of classical complement component C1 with parasite and host calreticulin mediates Trypanosoma cruzi infection of human placenta. PLoS Negl Trop Dis. 2013;7:e2376.
Mayhew TM. Turnover of human villous trophoblast in normal pregnancy: what do we know and what do we need to know? Placenta. 2014;35:229–40.
This work was supported by ERANET-LAC grant ELAC2014/HID-0328 (to UK and AGS), REDES130118 (to UK and CR), UREDES URC-024/16 (to UK) and FONDECYT 3180452 (to ChC).
Availability of data and materials
The datasets generated during and analyzed during the current study are available in the GEO repository under the accession number GSE113155.
Ethics approval and consent to participate
This study was conducted in accordance with Ethic Committee of the Faculty of Medicine, University of Chile (No. 041-2011). Each patient signed an informed consent for experimental use.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Interaction networks in differentially expressed genes from each experimental group. HPCVE were incubated during 2 and 24 h with 105 or 106 T. cruzi trypomastigotes. Interaction networks from differentially expressed genes (FC ≥ 2) compared with uninfected control with GeneMANIA function prediction service plug-in in Cytoscape software. Co-expression, co-localization, physical interactions, genetic interactions shared protein domains and pathways are shown and each color represents specific interactions according to legend. (PDF 707 kb)