Erica Silberstein
Charles C. Chung
Alain Debrabant- 1Laboratory of Emerging Pathogens, Office of Blood Research and Review, Center for Biologics Evaluation and Research, Food and Drug Administration, Silver Spring, MD, United States
- 2High-performance Integrated Virtual Environment Team, Office of Biostatistics and Epidemiology, Center for Biologics Evaluation and Research, Food and Drug Administration, Silver Spring, MD, United States
Vertical transmission of Trypanosoma cruzi (T. cruzi) become a globalized health problem accounting for 22% of new cases of Chagas disease (CD). Congenital infection is now considered the main route of CD spread in non-endemic countries where no routine disease testing of pregnant women is implemented. The main mechanisms that lead to fetal infection by T. cruzi remain poorly understood. Mother-to-child transmission may occur when bloodstream trypomastigotes interact with the syncytiotrophoblasts (SYNs) that cover the placenta chorionic villi. These highly specialized cells function as a physical barrier and modulate immune responses against pathogen infections. To model the human placenta environment, we have previously used a three-dimensional (3D) cell culture system of SYNs that exhibits differentiation characteristics comparable to placental trophoblasts. Further, we have shown that 3D-grown SYNs are highly resistant to T. cruzi infection. In this work, we used RNA sequencing and whole transcriptome analysis to explore the immunological signatures that drive SYNs’ infection control. We found that the largest category of differentially expressed genes (DEGs) are associated with inflammation and innate immunity functions. Quantitative RT-PCR evaluation of selected DEGs, together with detection of cytokines and chemokines in SYNs culture supernatants, confirmed the transcriptome data. Several genes implicated in the Toll-like receptors signaling pathways were upregulated in 3D-grown SYNs. In fact, TLR2 blockade and TLR3/7 knockdown stimulated T. cruzi growth, suggesting that these molecules play a significant role in the host cell response to infection. Ingenuity Pathway Analysis of DEGs predicted the activation of canonical pathways such as S100 protein family, pathogen induced cytokine storm, wound healing, HIF1α signaling and phagosome formation after T. cruzi exposure. Our findings indicate that SYNs resist infection by eliciting a constitutive pro-inflammatory response and modulating multiple defense mechanisms that interfere with the parasite’s intracellular life cycle, contributing to parasite killing and infection control.
1. Introduction
Chagas disease (CD) is a lifelong zoonotic disease caused by the protozoan parasite Trypanosoma cruzi (T. cruzi) that affects six to seven million people worldwide (Perez-Molina and Molina, 2018). CD is endemic in Latin America where the parasite is mainly transmitted to humans and domestic animals by infected triatomine bugs. T. cruzi infection can also be transmitted through blood transfusion, organ transplantation, and from mother-to-baby. In non-endemic areas, T. cruzi congenital infection has become the main active route of CD spread (Carlier et al., 2019; Rios et al., 2020) with cases reported in the US, Canada, Europe and Japan (Imai et al., 2019; Irish et al., 2022; Navarro et al., 2022).
An estimated 43,000 T. cruzi-infected women of reproductive age live the United States where the vertical transmission rate is approximately 1–5% (Edwards et al., 2019; Irish et al., 2022). Congenital infections can occur both during the acute and chronic phases of maternal infection, but transmission rates increase in pregnant women with high parasitemia (Carlier et al., 2020; Klein et al., 2021). Most infected newborns do not present clinical signs at birth although in exceptional cases, they may experience severe symptoms that can lead to death (Matthews et al., 2022). Early diagnosis and therapeutic interventions are crucial to prevent cardiac or gastrointestinal complications later in life (Carlier et al., 2019; Perez-Zetune et al., 2020).
During pregnancy, the placenta forms a physical barrier between the mother and developing fetus. A single layer of multinucleated cells, the syncytiotrophoblasts (SYNs), covers the placenta free-floating chorionic villous trees that are in direct contact with maternal blood. The SYNs not only allow the exchange of nutrients and waste but also produce pregnancy hormones and immunoregulatory molecules that prevent pathogen fetal infections (Arora et al., 2017; Ander et al., 2019; Kemmerling et al., 2019).
Previous studies, based on histological analysis of infected human placentas and placental explants, indicate that T. cruzi induces destruction and detachment of SYNs (Duaso et al., 2010; Keogh et al., 2021; Mezzano et al., 2022). Other investigators have found parasitic lesions in placental tissues at the marginal zone (Fernandez-Aguilar et al., 2005). These observations suggest that T. cruzi may reach umbilical and fetal capillaries either by infecting SYNs or through invasion of the marginal sinus epithelial cells (Carlier et al., 2020; Avalos-Borges et al., 2022).
The main mechanisms and immunoregulatory factors that modulate T. cruzi fetal infection have yet to be determined. Several animal and ex vivo models have been used to study CD congenital transmission (Avalos-Borges et al., 2022). For instance, modulation of the NF-κB and Toll-like receptors (TLRs) pathways, as well as the role of the placental barrier integrity in parasite infection have been evaluated in chorionic villi explants (Diaz-Lujan et al., 2016; Liempi et al., 2016; Castillo et al., 2017a,b) and cord blood cells (Ait Djebbara et al., 2023). Yet, the ability to reproduce parasite infection in an environment that resembles the complex architecture of the human placenta continues to be a challenge for the field (Barrila et al., 2018; Torres-Vargas et al., 2018; Ander et al., 2019).
To simulate the human placenta environment, we and others have used the rotating wall vessel bioreactor (RWV) platform to culture SYNs (McConkey et al., 2016; Corry et al., 2017; Silberstein et al., 2021). When grown under the dynamic culture conditions provided by this system, SYNs differentiate into 3D tissue-like spheroids that form syncytia, produce placental-specific hormones and show a transcriptome profile strikingly similar to primary human SYNs (McConkey et al., 2016). Furthermore, and consistent with the intrinsic protective function of the placenta (Hoo et al., 2020), we have shown that 3D-grown SYNs are highly resistant to T. cruzi infection (Silberstein et al., 2021).
In this study, we sought to explore T. cruzi-host cell interactions and determine the immunological hallmark of 3D-grown SYNs in response to parasite infection using RNA sequencing followed by whole transcriptome analysis and functional studies. Our data suggest that SYNs’ exposure to T. cruzi results in the activation of multiple defense mechanisms related to parasite sensing, internalization and phagocytosis, which can ultimately culminate with killing of the parasite by the host cell.
2. Materials and methods
2.1. Cells
2D-cultured JEG-3 cells (2D SYNs; ATCC® HTB-36™, American Type Culture Collection, VA), human brain microvascular endothelial cells [HBMECs; (Coyne et al., 2007)], LLC-MK2 cells (ATCC® CCL-7™; American Type Culture Collection, VA), and three-dimensional (3D) cultures of JEG-3 cells (3D SYNs) were grown as described previously (Silberstein et al., 2021). Experiments with 3D SYNs and 3D HBMECs were carried out between day 20–22 and day 4–6 after culture initiation, respectively. The number of viable cells was determined using the CellTiter-Fluor™ (for 2D SYNs) or the CellTiter-Glo™ 3D (for 3D SYNs) cell viability assays following manufacturer’s instructions (Promega, WI).
2.2. Trypanosoma cruzi propagation and infection
Trypomastigotes of the T. cruzi strain Colombiana expressing nanoluciferase [TcCOL-NLuc (Silberstein et al., 2018)] were harvested from culture supernatants of infected LLC-MK2 cells and the number of parasites was determined using a Cellometer K2 Fluorescent Viability Cell Counter, following manufacturer’s instructions (Nexcelom Bioscience, MA). Infection of 2D SYNs was performed in 96-well plates (seeded at 1×104 cells /well) or T75 culture flasks using a ratio of 5 parasites/cell (MOI = 5). For cultures grown in 3D, spheroids were plated in 96-well plates (100 μL/well) or 6 well plates (2 mL/well) and cells were infected at MOI = 5. T. cruzi growth was monitored in 96-well plates using the Nano-Glo® Luciferase Assay following manufacturer’s instructions (Promega, WI) as previously described (Silberstein et al., 2021).
2.3. RNA sequencing and data processing
Total RNA was extracted from harvested 3D and 2D-cultured SYNs, and 3D-cultured HBMECs 18 h post-infection, using the PureLink™ RNA mini kit (Invitrogen, CA). RNA integrity number (RIN) and concentration for each sample was determined with the Agilent 2,100 Bioanalyzer (Agilent Technologies, CA). All RIN values were within the acceptable range (RIN ≥ 9). Two to four replicates from cells exposed and not exposed to T. cruzi were used to perform transcriptome sequencing at Novogene1 according to their protocols. Briefly, oligo d(T) beads were employed to perform mRNA enrichment, which were then fragmented randomly and converted into cDNA. Poly(A) tails were added and enriched by PCR to generate a cDNA library using NEBNext® Ultra™ RNA Library Prep Kit for Illumina® (NEB, MA). Illumina HiSeq 2,500 was used to perform paired-end sequencing of 150 bp reads. The read quality was analyzed with fastp (Chen et al., 2018) trimming polyG tail, polyX in 3′ ends and filtered reads shorter than length 20. Quality passed reads were aligned to the human reference genome (GRCh38.p13) using HISAT2 v2.1.0 (Kim et al., 2019). The mapped reads were counted with featureCounts (Liao et al., 2014) for genes annotated in the GCF_000001405.39 annotation file from NCBI.
2.4. Differential gene expression and pathway analysis
The differential gene expression (DGE) analysis was conducted for the 3D SYNs (n = 4) and 2D SYNs (n = 2), and for the 3D HBMECs (n = 3) and 2D SYNs (n = 2), unexposed or exposed to the parasite for 18 h using the DESeq2 R package (Love et al., 2014). The alpha and the fold change thresholds were set to a Benjamin–Hochberg-adjusted p-value <0.05 and 1.5, respectively for significance in the DGE test. Read counts were normalized using the median of ratios method. DEG cluster analysis was performed using the pheatmap R package.
The log2 (fold-change) expression values of DEGs and p-adjusted values were uploaded into the Ingenuity Pathway Analysis software [IPA; QIAGEN (Kramer et al., 2014)] to determine activated or inhibited canonical pathways. Z-scores were used to predict pathway activation (z-score ≥ 2) or inhibition (z-score ≤ 2). Statistical significance was calculated using a right-tailed Fisher Exact test and p--values were adjusted for multiple hypothesis testing using the Benjamin–Hochberg method. An absolute B-H p-value <0.05 was considered significant.
2.5. Quantitative RT-PCR
Following purification, total RNA was reverse- transcribed into cDNA using the iScript gDNA clear cDNA synthesis kit (Bio-Rad Laboratories, CA) and 0.5 to 1 μg of RNA template. Quantitative PCR was carried out using the SsoAdvanced™-Universal SYBR® Green Supermix (Bio-Rad Laboratories, CA) in a CFX96 touch real-time PCR detection system (Bio-Rad Laboratories, CA). Eighty-four genes associated with host immune defenses were profiled using the Innate and Adaptative Immune Responses PrimePCR Array (Bio-Rad Laboratories, CA). The cycling protocol included an initial activation step at 95°C for 2 min, followed by 40 cycles of 5 s at 95°C and 30 s at 60°C. Finally, a melt curve analysis was run for 5 s/step between 65 and 95°C with 0.5°C increments. Data were analyzed with CFX Manager Software. Gene expression was determined by the ∆Ct method with samples normalized to human GAPDH (Livak and Schmittgen, 2001).
2.6. Bio-Plex assay for cytokines profiling
Cytokine levels in culture supernatants of 3D SYNs (unexposed or exposed to T. cruzi) were measured using a Bio-Plex Pro Human Cytokine 17-plex Assay (Bio-Rad Laboratories, CA), according to the manufacturer’s instructions. Briefly, 50 μL of coupled magnetic beads mixture was added to each well. After three washes, 50 μL of undiluted culture supernatants and serial dilutions of standards were added to duplicate wells. The plate was incubated on a microplate orbital shaker at 850 rpm for 2 h and next washed three times before adding 50 μL of detection antibodies conjugated with biotin. Following a 30 min incubation, wells were washed again three times, and 50 μL streptavidin-phycoerythin solution was added. After washing three times, beads were resuspended in assay buffer and incubated for 10 min. All incubations were carried out at room temperature. Data were acquired using a Luminex 200 instrument (Luminex Corporation, TX) and analyzed with the BioPlex Pro software (Bio-Rad Laboratories, CA). Concentration of each analyte in pg/ml was calculated from corresponding standard curves.
2.7. TLR2 blockade
2D SYNs (7.5×103 cells) were seeded into each well of 96 well plates. After 1 h incubation with 5 mg/mL anti-hTLR2 (InvivoGen, CA) or RecombiMAb human IgG1 (LALA-PG) isotype control (Bio X Cell, NH), cells were infected with TcCOL-NLuc trypomastigotes using a ratio of 5 parasites/cell. After an overnight incubation, free parasites were removed, and cells were washed three times with culture media containing 2% fetal bovine serum (FBS). Next, 100 μL of fresh culture media containing 5 mg/mL of antibodies were added per well and parasite growth was determined at 48- and 72 h post-infection by measuring nanoluciferase (NLuc) activity in triplicate wells as described (Silberstein et al., 2021). Cells from four replicate wells were pooled and total RNA was extracted at identical time points. Quantitative RT-PCR to measure cytokines mRNA was carried out using the SsoAdvanced™-Universal SYBR® Green Supermix (Bio-Rad Laboratories, CA) and primer sets specific to human IL6 and TNF-α (Bio-Rad Laboratories, CA). Gene expression was determined by the ∆Ct method with samples normalized to human GAPDH (Livak and Schmittgen, 2001).
2.8. siRNA mediated knockdown of TLR3 and TLR7
On-target plus non-targeting siRNA control pool and on-target plus human TLR3 and TLR7 siRNA smartpools were purchased from Horizon (Lafayette, CO). Transfections were performed according to the manufacturer’s instructions. In brief, 2D SYNS (1×104) were seeded into each well of 96-well plate using antibiotic-free complete medium. Each transfection was carried out in three replicate wells using a 100 nM final concentration of siRNA (50 nM of TLR3 siRNA and 50 nM of TLR7 siRNA) mixed with 0.3 μL/well of DharmaFECT 3 transfection reagent (Horizon, CO). After 24 h, cells were infected with TcCOL-NLuc trypomastigotes using a ratio of 5 parasites/cell. Free parasites were removed after 24 h post-infection and cells were washed three times with culture media. Parasite growth was determined at 48- and 72 h post-infection by measuring nanoluciferase (NLuc) activity as described (Silberstein et al., 2021). Cells from triplicate wells were pooled and total RNA was extracted at identical time points. Quantitative RT-PCR to verify TLRs silencing and to measure cytokines mRNA, was carried out using the SsoAdvanced™-Universal SYBR® Green Supermix (Bio-Rad Laboratories, CA) and primer sets specific to human TLR3, TLR7, IRF3, IRF7, IFN-α1 and IFN-β1 (Bio-Rad Laboratories, CA). Gene expression was determined by the ∆Ct method with samples normalized to human GAPDH (Livak and Schmittgen, 2001).
2.9. Statistical analysis
All bar graphs are displayed as means ± SD. Statistical analysis of differences between mean values of groups was determined by the unpaired two-tailed Student’s t-test, using GraphPad Prism software version 9.0.0. A p value <0.05 was considered significant.
3. Results
3.1. Transcriptomic profiling reveals distinct gene expression changes in 3D-grown SYNs compared to 2D-grown
To identify the possible mechanisms responsible for the resistance of 3D-grown SYNs to T. cruzi infection, we performed high throughput RNA sequencing (RNA-Seq) on four independent SYNs cultures grown in 3D and two independent 2D SYNs cultures, unexposed or exposed to the parasite for 18 h (Figure 1A). We used a ratio of 5 parasites/cell, conditions representative of a chronic infection characterized by moderate to low parasitemia (Carlier et al., 2020).
Figure 1. Alterations in global gene expression profiles of 3D and 2D-grown SYNs. 2D or 3D SYNs were infected with T. cruzi at MOI = 5 and harvested at 18 hpi for transcriptomic analysis by RNA-Seq. Differentially expressed genes were identified using the DEseq2 software. The top signaling pathways activated in 3D SYNs were predicted through Ingenuity Pathway Analysis. Functional studies were carried out to verify the impact of selected genes in the regulation of the identified biological pathways. (A) Experimental workflow for the whole transcriptome analysis of 2D and 3D-grown SYNs cultures unexposed or exposed to T. cruzi. (B) Principal component analysis (PCA) of RNA-Seq data of 2D and 3D-cultured SYNs unexposed or exposed to T. cruzi. (C) Hierarchical cluster analysis of normalized read counts of differentially expressed genes between 3D and 2D SYNs unexposed or exposed to T. cruzi. (D) Number of significantly up- and down-regulated genes in 3D SYNs with respect to 2D SYNs (fold-change ≥1.5 and FDR < 0.05).
When we compared the transcriptome profile of T. cruzi -exposed 2D SYNs with unexposed 2D cultures, and T. cruzi -exposed 3D SYNs with unexposed 3D cultures, we observed few changes in the abundance of transcripts. These results correlate with previous findings that showed a weak modulation of host cell transcription in response to trypomastigotes invasion (Vaena de Avalos et al., 2002; Perumal et al., 2023).
With the goal of finding the main factors that drive 3D SYNs resistance to the parasite, we next studied the transcriptome profile of 3D SYNs compared to 2D SYNs, in both unexposed and T. cruzi exposed cultures. Interestingly, we found many transcriptional changes in 3D SYNs compared to 2D SYNs. As illustrated by the principal component analysis (PCA), two groups could be clearly discriminated with samples closely clustering together in each group (Figure 1B, 3D vs. 2D, unexposed or T. cruzi-exposed).
We further examined gene expression data using hierarchical clustering analysis of the DEGs. The clustering patterns correlate with the PCA observations, showing the genes with similar patterns of expression, and a strong contrast between 3D and 2D-cultured SYNs (Figure 1C, 3D vs. 2D, unexposed or T. cruzi-exposed). When we compared 3D vs. 2D SYNs cultures unexposed to T. cruzi (Figure 1D), we found 7,323 DEGs (4,334 upregulated genes and 2,989 downregulated genes). After T. cruzi exposure the number of DEGs increased to 8,318 genes, with 4,817 and 3,501 genes being up and downregulated, respectively.
To confirm the RNA-Seq data observations, we profiled the expression of selected genes by quantitative RT-PCR using arrays targeting molecules associated with host immune defenses (Figure 2). We observed upregulation of CASP1, IL6, CXCL8, IRF7, NOS2, TLR2, TLR3 and TLR7, which validates our whole transcriptome analysis findings. However, transcriptional levels of these genes remained nearly unchanged after exposure to T. cruzi. These data suggest that 3D SYNs constitutively express high levels of molecules involved in the activation of innate immune responses and correlate with the capacity of placental trophoblasts to control microbial infections (Ander et al., 2019; Megli et al., 2021; Semmes and Coyne, 2022).
Figure 2. Validation of RNA-Seq data by qRT-PCR. mRNA quantification of selected genes carried out by qRT-PCR in unexposed or T. cruzi- exposed 3D SYNs. Violin-plot distributions showing log2 fold changes in gene expression in 3D SYNs compared to 2D SYNs as determined by the ∆Ct method with samples normalized to GAPDH. Dotted lines represent first and third quartiles, and solid lines indicate the median (values shown over plots). The dashed blue line indicates the assay cut off value (0.58). The qRT-PCR experiment was repeated two times with similar results. The table shows the log2 fold change in expression of 3D SYNs vs. 2D SYNs obtained through RNA sequencing analysis. FC: Fold change expression; U: unexposed; TC: exposed to T. cruzi.
3.2. Innate immune response signaling pathways are constitutively activated in 3D-cultured SYNs
To characterize 3D SYNs responses following T. cruzi exposure, we analyzed the significant changes in gene expression using Ingenuity Pathway Analysis (IPA).
We first evaluated the canonical pathways activated in cultures that were not exposed to T. cruzi, comparing 3D SYNs to 2D SYNs (Figure 3A). The top 15 more significantly activated pathways included pathogen induced cytokine storm, wound healing, role of JAK family kinases in IL-6 type cytokine signaling, regulation of the epithelial mesenchymal transition by growth factors, S100 protein family, calcium signaling and role of MAPK signaling in inhibiting the pathogenesis of Influenza. Consistent with previous reports (Ander et al., 2019; Semmes and Coyne, 2022), our findings indicate that 3D SYNs not only constitutively express molecules associated with innate immunity, but they also exhibit wound repair and calcium-dependent signal transduction activities.
Figure 3. Pathway analysis reveals activation of processes associated with cellular immune defenses, pathogen-influenced signaling and cellular stress/injury in 3D-cultured SYNs. (A) Top 15 significant canonical pathways identified by IPA in 3D SYNs compared to 2D SYNs in cultures not exposed to T. cruzi. (B) Top significant canonical pathways identified by IPA in 3D SYNs compared to 2D SYNs in T. cruzi-exposed cultures. The red-colored bars represent the -log (B-H p-value) calculated for each pathway. (C) Heat-map showing the activation z-scores of the top canonical pathways in unexposed and T. cruzi-exposed cultures. Asterisks indicate the significance of each pathway comparing 3D and 2D SYNs (* B-H p-value <0.05; ** B-H p-value <0.01; *** B-H p-value <0.001). Only pathways with z-score ≥ 2 and B-H p-value ≥0.05 were considered significant. U: unexposed; TC: T. cruzi-exposed. ns: no significant B-H p-values. #: Pathways enriched in parasite-exposed 3D SYNs compared to 2D SYNs.
When we analyzed the top canonical pathways enriched in cultures after T. cruzi exposure (comparing 3D SYNs to 2D SYNs, Figure 3B), we observed activation of additional pathways such as G-protein coupled receptor signaling, HIF1α signaling and phagosome formation. Of note, S100 protein family signaling was the most significantly activated pathway followed by pathogen induced cytokine storm, calcium signaling and wound healing.
We next conducted a comparison analysis between unexposed and T. cruzi-exposed cultures to determine if there were differences in the z-scores and significance values (B-H p-values) between the top enriched pathways (Figure 3C; Supplementary Table 1). Based on the statistical analysis, pathways such as S100 protein family, pathogen induced cytokine storm, calcium signaling, and wound healing showed similar activation when cultures were unexposed or exposed to T. cruzi. However, G-protein coupled receptor signaling, HIF1α signaling and phagosome formation pathways were highly enriched in parasite-exposed 3D SYNs compared to 2D SYNs.
To initiate the culture of human SYNs in 3D, trophoblasts are co-cultured with human brain microvascular endothelial cells (HBMECs) attached to microcarrier beads. It was previously demonstrated that 3D-grown SYNs induce the dissociation of HBMECs from beads (McConkey et al., 2016) and form a layer that coats the surface of the spheroids (McConkey et al., 2016; Silberstein et al., 2021). Since some HBMECs may perhaps remain attached to beads after 20 days in culture, we also carried out RNA Seq analysis of HBMEC cultures grown in 3D. Next, we run an IPA core analysis to determine the canonical pathways activated in 2D SYNs compared to 3D HBMECs, unexposed or T. cruzi-exposed. Importantly, a comparison analysis of z-scores across all groups (3D vs. 2D SYNs and 3D HBMECs vs. 2D SYNs, Supplementary Figure 1) showed that the 3D SYNs top enriched canonical pathways where downregulated or not significant in 3D-grown HBMECs. Therefore, it’s unlikely that the presence of some residual HBMECs RNA in the 3D SYNs RNA preparation may possibly impact our transcriptome analysis results.
3.3. Trypanosoma cruzi modulates the gene expression landscape of 3D-cultured SYNs
To better understand the different processes triggered in 3D SYNs after exposure to T. cruzi, we generated heatmaps of the most significantly regulated DEGs (log2 fold change ≥ ± 1.5; p-value <0.05) in the top enriched canonical pathways (Figure 4; Supplementary Figure 2; Supplementary Tables 2–6).
Figure 4. Toll-like receptors and cytokines are common genes in the top activated canonical pathways. Heat-maps of top differentially expressed genes for (A) S100 family, (B) pathogen induced cytokine storm, (C) wound healing, (D) HIF1α signaling and (E) phagosome formation pathways. Only genes with log2 fold change ≥ ± 1.5 and p-value ≥0.05 were considered significant. ns: no significant log2 fold change or B-H p-values. U: unexposed; TC: T. cruzi-exposed. (F) Levels of IL6, CCL2, TNF-α, GMCSF and IFN-γ determined using a Bio-Plex Pro Human Cytokine Screening Panel in culture supernatants of unexposed (blue) and T. cruzi-exposed (red) 3D-grown SYNs. Boxes represent first and third quartiles with a line at the median. Whiskers indicate maximum and minimum data values. CA: cytosolic adapter protein; C: cytokines and cytokine receptors; E: enzymes; GF: growth factors; TR: transcription regulators.
We first examined the gene expression patterns of molecules belonging to the S100 family signaling pathway, a family of calcium-binding proteins that are activators of immune functions including inflammation and also play an important role in calcium homeostasis (Xia et al., 2017; Singh and Ali, 2022). Of note, T. cruzi internalization is reported to be a Ca2+ dependent process (Fernandes and Andrews, 2012; Ferri and Edreira, 2021). We found that the most up-regulated molecules were CXCR6 and KISS1R (G-protein coupled receptors), IL10RA (cytokine receptor) and STAT3 (transcription regulator) while genes involved in cytoskeleton assembly (S100A6 and S100A9), calcium transport (CACNG7) and CXCL8 (cytokine) were down-regulated (Figure 4A; Supplementary Figure 2A; Supplementary Table 2).
Next, we studied the genes that were present in our data sets which are associated with pathogen induced cytokine storm signaling, a process characterized by uncontrolled inflammation caused by cytokine excess (Karki and Kanneganti, 2021). We observed up-regulation of CXCL2 (chemokine), and NGFR (TNF receptor). Expression of IL17C (cytokine) was down-regulated whereas transcription levels of other innate immunity-related molecules and factors (TLR3, TLR7, IRF7, NOS2, CASP1) only showed modest changes (Figure 4B; Supplementary Figure 2B; Supplementary Table 3).
We also found changes in the expression of several genes associated with the wound healing pathway, an important process required for T. cruzi internalization into host cells (Ferri and Edreira, 2021). Specifically, KRT16 and COL8A2 (structural proteins), EGF (growth factor), and IL11 (cytokine) were upregulated. In addition, a marked down-regulation of structural constituents VIM and LAMA3 were observed following T. cruzi exposure (Figure 4C; Supplementary Figure 2C; Supplementary Table 4).
Consistent with previous findings in human cardiomyocytes (Venturini et al., 2023), HIF1α signaling was also activated in 3D-grown SYNs exposed to the parasite. Transcripts encoding FLT1 (growth factor receptor), VEGFA (growth factor), HSPA6 (heat shock protein), HK2 (hexokinase), MMP17 (peptidase) and PRKCG (kinase) were the most enriched. PRKD1 (kinase) and RAC2 (GTPase) were among the most down-regulated genes (Figure 4D; Supplementary Figure 2D; Supplementary Table 5).
Phagocytosis is a fundamental process in immunity against pathogens (Rosales and Uribe-Querol, 2017) and T. cruzi exploits this mechanism to facilitate host cell invasion (Fernandes and Andrews, 2012; Ferri and Edreira, 2021). In the phagosome formation pathway, we identified genes such as GPRC5A (G-protein coupled receptor), FCGR2B (Fc receptor), PLA2G4D and PLD6 (phospholipases), LBP (lipopolysaccharide-binding protein) and CARD9 (caspase adaptor protein) whose transcript levels were greater in T. cruzi-exposed 3D SYNs than in unexposed cultures. We also found down-regulation of BDKRB2(G-protein coupled receptor), ITGA3, MRC1 and PLA2R1 (transmembrane receptors), and MYL7(calcium-binding protein) (Figure 4E; Supplementary Figure 2E; Supplementary Table 6).
To determine the commonly expressed genes across the top predicted canonical pathways we compared the lists generated by the IPA core analysis using a multiple list comparator tool2 (Supplementary Tables 6, 7). Among the genes shared by two or more pathways, we found molecules that are known to be involved in the innate immunity response to T. cruzi (Macaluso et al., 2023) such as TLR3 and TLR7 (shared by the pathogen induced cytokine storm and phagosome formation signaling pathways) and NOS2, STA3 and IL6 (shared by the S100 family, pathogen induced cytokine storm and HIF1α signaling pathways).
To further study the immunological profile of 3D SYNs, we measured cytokine and chemokine levels in culture supernatants using a multianalyte screening panel. In agreement with previous reports (Ander et al., 2018; Megli et al., 2021), we found that high levels (>1,000 pg/mL) of IL6, CCL2 and TNF-α were released to the culture media (Figure 4F, left panel). GM-CSF, IL5 and IFN-γ were also secreted but at lower levels (< 300 pg/mL; Figure 4F, right panel). Notably, constitutive expression of these molecules was not significantly altered following exposure to T. cruzi.
Together, the findings described above (Figures 3 4) strongly suggest that exposure of 3D SYNs to T. cruzi leads to activation of processes related to parasite sensing, internalization and phagocytosis (Fernandes and Andrews, 2012; Ferri and Edreira, 2021).
3.4. Blockade of TLR2 boosts T. cruzi growth
Toll-like receptors play an important role in the elimination of microbes through recognition of pathogen-associated patterns and induction of inflammatory cytokines and type I interferons (Fitzgerald and Kagan, 2020; Motomura et al., 2023). T. cruzi molecules such as surface glycoinositolphospholipids (GPIs) and parasite DNA/RNA sequences can be sensed by TLRs (Macaluso et al., 2023). In fact, several investigators have shown the significance and involvement of TLRs in the control of T. cruzi infection (Caetano et al., 2011; Castillo et al., 2017b; Queiroga et al., 2021; Ait Djebbara et al., 2023).
Because we found that molecules associated with the TLR2 pathway (IL6, CXCL8 and CXCL2), and TLR3 and TLR7 were differentially regulated in 3D SYNs compared to 2D SYNs (Figure 4; Supplementary Figure 2), we explored their function as mediators of trophoblasts’ immunity to T. cruzi. For these experiments, we could not employ 3D SYNs because they resist parasite infection (Silberstein et al., 2021). Therefore, we used 2D SYNs as working model since they not only express TLRs (although to lower levels compared to 3D SYNs) but also support T. cruzi growth.
To study the effect of TLR2 blockade on parasite infection, 2D SYNs were pre-treated with anti-hTLR2 antibodies or isotype control antibody (5 μg/mL each) and next infected with TcCOL-NLuc trypomastigotes. Parasite growth was monitored at 48 and 72 h post infection (hpi) (Figure 5A) by measuring nanoluciferase activity as described previously (Silberstein et al., 2021). We found that TLR2 inhibition significantly promoted parasite growth as indicated by the increase in nanoluciferase activity in cells treated with anti-TLR2 antibody at both 48 and 72 hpi. Since TLR2 activation triggers expression of IL6 and TNF-α, we explored the downstream effects of TLR2 inhibition by measuring their mRNA levels by qRT-PCR (Figure 5B). The most reduced cytokine in anti-TLR2 antibody treated 2D SYNs was IL6 (fourfold down-regulation), at 48 hpi while a marked decrease in TNF-α was found at 72 hpi.
Figure 5. Inhibition of TLRs significantly decreases expression of cytokines and interferons promoting T. cruzi growth. 2D SYNs were treated with anti-hTLR2 (TLR2 blockade) or transfected with TLR3/TLR7 siRNA mix, followed by T. cruzi infection at ratio of 5 parasites/cell. (A) T. cruzi growth over time after TLR2 blockade or (C) TLR3/TLR7 siRNA transfection. Bars show the nanoluciferase activity expressed in relative luminescence units (RLU) measured at 48 and, 72 hpi and represent the mean ± SD of three technical replicates. The experiment was repeated two times with similar results. Asterisks indicate statistically significant differences calculated using an unpaired t-test. *: p ≤ 0.05; **: p ≤ 0.01. (B) qRT-PCR mRNA quantification after TLR2 blockade or (D) TLR3/TLR7 siRNA transfection. Fold changes in gene expression compared to cells treated with isotype control antibodies (B) or to cells transfected with non-targeting siRNA control (D), were determined by the ∆Ct method with samples normalized to GAPDH. Violin-plots show the log2 fold change expression vs. control (mean ± SD) of two technical replicates. Dotted lines represent first and third quartiles, and solid lines indicate the median (values shown below or above the plots). The dashed blue line shows the assay cut off value (0.58). The experiment was repeated two times with similar results.
3.5. TLR3 and TLR7 siRNA silencing stimulates T. cruzi expansion
We also investigated the outcome of TL3 and TLR7 double knockdown on T. cruzi infection. To test this, 2D SYNs were transfected with TLR3/TLR7 siRNAs or non-targeting siRNA (control). Gene silencing was verified by measuring TLR3 and TLR7 mRNA quantities using qRT-PCR (Supplementary Figure 3). Twenty-four hours post-transfection, cells were infected with TcCOL-NLuc trypomastigotes and parasite multiplication was monitored as explained above. We found that knockdown of TLR3 and TLR7 also promoted parasite growth as evidenced by the significant increase in nanoluciferase activity in double knock-downed cells (Figure 5C). Given that TLR3 and TLR7 downstream signaling involves activation of interferon regulatory factors (IRFs), which in turn trigger synthesis of type I interferon genes, we measured IRF3, IRF7, IFN-α and IFN-β mRNA levels by qRT-PCR (Figure 5D). The knockdown effect on cytokine mRNAs expression was observed for IRF3, IFN-α and IFN-β at 72 h post-transfection, with a remarkable 2.5-fold reduction in IFN-α.
These data (Figure 5) suggest that inhibition of TLR2, TLR3 and TLR7 causes a reduction in the expression of immunoregulatory molecules, which results in increased parasite growth. Furthermore, our findings highlight the significant role played by TLRs in SYNs’ innate response and resistance to T. cruzi infection.
4. Discussion
Congenital T. cruzi infection is now the main route of CD transmission in non-endemic countries where pregnant women are not routinely tested for the disease (Antinori et al., 2017). A complex parasite–host cell interplay occurs at the maternal-fetal interface when bloodstream trypomastigotes interact with the SYNs (Blaszkowska and Goralska, 2014). Understanding the molecular basis of T. cruzi invasion and replication in SYNs could provide insights into the specific molecules and host signaling pathways that modulate parasite transmission from mother-to-fetus.
Although composed by different trophoblasts populations (SYNs, cytotrophoblast, and extravillous trophoblasts), maternal macrophages and Hofbauer cells (Semmes and Coyne, 2022), human chorionic villi explants have been previously used to examine the transcriptome changes induced after T. cruzi challenge (Castillo et al., 2018). However, the signaling pathways activated in monotypic cultures of SYNs in the context of parasite infection, have not been investigated. Using 3D-grown cultures, and through RNA sequencing and whole transcriptome analysis, we characterized the SYN’s response to T. cruzi and confirmed their key role as immune sensors of parasite infection.
Consistent with previous work, our data revealed significant activation of signaling pathways associated with the innate immune response (Castillo et al., 2018; Barbosa et al., 2022). Notably, we observed that pathways such as S100 protein family, pathogen induced cytokine storm, calcium signaling, and wound healing were similarly enriched in unexposed or T. cruzi exposed 3D SYNs compared to 2D SYNs (Figure 3C). Indeed, induction of these canonical pathways was verified through the detection of cytokines and chemokines (Figure 4F) in the culture media, as reported by others (Vekemans et al., 2000; Arora et al., 2017; Corry et al., 2017; Ander et al., 2019). We also found overexpression of caspase 1 (CASP1, Figures 2, 4B), a key driver of placental antimicrobial defenses through the inflammasome signaling pathway (Megli et al., 2021). These findings suggest that 3D SYNs exhibit a pro-inflammatory state that shapes their resistance to T. cruzi and other microbial infections.
During its life cycle inside the host-cell, T. cruzi hijacks different cell functions to establish a productive infection. The main steps of parasite invasion include host cell recognition and adhesion, internalization, formation of and escape of parasites from the phagolysosome, multiplication of amastigotes in the cytosol where they transform into trypomastigotes that can infect new cells (Ferri and Edreira, 2021). Based on our data, it is likely that 3D SYNs resist T. cruzi infection by modulating the expression of host molecules linked to the top canonical pathways identified in this work. These proteins may play a role in inhibiting the parasite’s intracellular cycle at different stages.
For instance, T. cruzi can be recognized by cell surface TLRs and some G protein coupled receptor such as bradykinin B2 (BDKRB2), which in turn trigger innate defense mechanisms in the host cell (Schmitz et al., 2014; Castillo et al., 2017b; Macaluso et al., 2023). Notably, we observed up-regulation of TLR2 (Figure 2), and BDKRB2 (Figures 4A,E; associated with the S100 protein family pathway) in 3D SYNs. Blockade of TLR2 resulted in more T. cruzi replication (Figure 5A) likely due to a reduction in TLR2-mediated signaling and cytokines production (Figure 6). Interestingly, we have not detected significant overexpression of TLR4, reported to induce IFN-γ and TNF-α production in human cord blood cells after stimulation a recombinant T. cruzi macrophage infectivity potentiator (TcMIP) (Ait Djebbara et al., 2023). Likewise, we have not observed an increase in the level of nuclear factor kappa B (NF-κB) RNA transcripts, a molecule shown to be activated in placental explants via TLRs signaling pathways (Liempi et al., 2019). However, we detected significant up-regulation of the AP-1 transcription subunits FOSL1 and JUN in unexposed (log2 FC 3.5 and 4.2, respectively) and T. cruzi-exposed (log2 FC 3.6 and 4.2, respectively) SYNs (Supplementary Table 9). Activation of AP-1 leads to production of proinflammatory cytokines such us IL6 and TNF-α (Duan et al., 2022) and could explain the detection of these cytokines in supernatants of T. cruzi exposed SYNs (Figure 4F).
Figure 6. A proposed model of immune response mechanisms against T. cruzi mediated by Toll-like receptors in 3D-grown SYNs. T. cruzi invasion begins with binding to the SYNs extracellular matrix followed by TLR2 recognition and adhesion to the cell (1). This interaction promotes parasite internalization (2) into the parasitophorous vacuole and subsequent formation of the phagolysosome, where DNA and RNA digestion products are recognized by TLR3 and 7 (3). Parasite sensing by TLRs leads to synthesis of cytokines, chemokines and INFα/β (4), which in turn induce expression of NOS2 (5). Oxidation of arginine by NOS2 (6) generates nitric oxide which kills the parasite inside the phagolysosome or in the cytosol (7).
Another crucial event required for T. cruzi cell invasion is the mobilization of calcium (Ca2+) for lysosome recruitment. When trypomastigotes wound host cells to access the cytoplasm, they trigger Ca2+ influx from intracellular stores and the extracellular environment (Fernandes and Andrews, 2012; Ferri and Edreira, 2021). In agreement with former studies that evaluated the alterations of gene expression in placentas of pregnant women with chronic CD (Juiz et al., 2018), we found enrichment of the calcium signaling pathway (Figure 3). In addition, and as mentioned above, we also observed activation of the S100 protein family pathway, which includes Ca2+ binding and Ca2+ voltage-gated channel molecules. S100 proteins regulate diverse cellular processes including intracellular Ca2+ homeostasis (Sreejit et al., 2020) and their expression may be associated with HBV intrauterine transmission (Zhao et al., 2020). It is possible that SYNs alter their Ca2+ intracellular storage and mobilization limiting parasite invasion. However, the specific Ca2+ binding and transport proteins along with their role in the context of placenta infection by T. cruzi, remain to be investigated.
In line with the capacity of SYNs to experience a continuous and dynamic repair process (Castillo et al., 2018; Hoo et al., 2020), we found that the wound healing signaling pathway was significantly enriched in 3D SYNs. Because damage of the host cell plasma membrane is required for parasite entry (Brigido et al., 2017), our results suggest that SYNs’ constitutive healing mechanism can be actively inhibiting this important process. Additionally, we detected down-regulation of the expression of laminin-subunit alpha 3 (LAMA3) and vimentin (VIM) along with a marked up-regulation of cytokeratin 16 (KRT16) (Figure 4C; Supplementary Figure 2C). It’s been shown that Tc-85 proteins (members of the gp85/trans-sialidase superfamily) are involved in cell adhesion and parasite cell entry by interaction with molecules such as laminin (extracellular matrix constituent) and the intermediate filaments proteins cytokeratin and vimentin (Teixeira et al., 2015). Thus, downregulation of LAMA3 and VIM may possibly constitute another mechanism by which SYNs restrict parasite colonization while KRT16 over expression could be needed to maintain cell shape and integrity.
T. cruzi exploits several internalization pathways (including phagocytosis) for cell invasion (Ferri and Edreira, 2021). As reviewed by Rosales et al. (Rosales and Uribe-Querol, 2017), phagocytosis is a fundamental process in immunity which concludes with the formation of phagolysosomes where parasites are eventually digested and killed (Cardoso et al., 2015). In fact, we found significant activation of the phagosome formation pathway following exposure to the parasite (Figures 3B,C). We also observed upregulation of phospholipases PLD6 and PLA2G4D, enzymes known to be expressed in phagocytes (Griffiths and Mayorga, 2007; Dabral and van den Bogaart, 2021) and TLR3 and 7, pattern recognition receptors that detect nucleic acids from lysed parasites in the phagolysosome (Cardoso et al., 2015). Furthermore, we show that silencing of TLR3 and 7 increased T. cruzi infection (Figure 5C). Surprisingly, we have not found significant overexpression of TLR9 in 3D SYNs compared to 2D SYNs although it was reported that T. cruzi can also be sensed by intracellular TLR9 (Bafica et al., 2006; Gravina et al., 2013). Together, these findings indicate that altering phagocytosis-related pathways, including increased expression of TLR3 and 7, could be one more strategy to control T. cruzi colonization of 3D SYNs.
Previous reviews have explored the crosstalk between hypoxia-inducible factor (HIF-1α) and nitric oxide (NO) signaling pathways (Olson and van der Vliet, 2011; Pappas et al., 2023). Besides inducing the expression of proteins implicated in glucose metabolism such as HK2, HIF-1α can activate NOS2 and trigger NO production. Consistent with recent reports, our transcriptome analysis studies unveiled activation of the HIF-1α pathway (Figures 3B,C) along with overexpression of NOS2 (Figures 2 4D; Supplementary Figure 2D) and HK2 (Figure 4D; Supplementary Figure 2D) in parasite-exposed 3D SYNs (Triquell et al., 2018; Silberstein et al., 2021; Venturini et al., 2023). Upregulation of NOS2 may increase arginine oxidation which in turn could generate high NO levels to kill T. cruzi (Cardoso et al., 2015). Our results suggest that activation of the HIF-1α pathway by SYNs could also contribute to their resistance to T. cruzi.
It is important to note that even though 3D SYNs have characteristics similar to human syncytiotrophoblasts, including the ability to resist microbial infections (McConkey et al., 2016; Corry et al., 2017; Silberstein et al., 2021), this cell-based system may not completely reproduce in vivo responses to T. cruzi infection. In fact, the maternal-fetal interface is a complex environment where the chorionic villi are associated with the decidua containing different immune cell populations which can also impact immunity signaling pathways in SYNs (Semmes and Coyne, 2022). While SYNs have been reported to be a route of parasite transmission to fetal tissues, it has also been shown that T. cruzi invasion can occur through infection of non-trophoblastic epithelial cells located at the placenta marginal sinus (Fernandez-Aguilar et al., 2005; Carlier et al., 2020). Our 3D cell culture system addresses only the role of SYNs in the control of T. cruzi vertical transmission, and therefore our observations may specifically apply to these trophoblastic cell type. Another limitation of this work is that the functional analysis experiments only evaluated the role of TLRs in controlling parasite growth in the context of a low parasite exposure representing the scenario of a chronic infection in pregnant women. Further studies using higher parasite loads and focused on the functionality of the activated signaling pathways and molecules predicted by IPA, may provide additional insights into the multiple resistance processes that limit T. cruzi infection of SYNs.
Our data support a model by which 3D SYNs may resist infection by modulating TLR signaling pathways (Figure 6). T. cruzi invasion begins with parasite binding to laminin (LAMA3) within the extracellular matrix. Next, parasites are recognized by TLR2 and adhere to the cell. This interaction promotes internalization into the parasitophorous vacuole and subsequent formation of the phagolysosome where some parasites are digested by enzymes and others escape into the cytoplasm. The DNA and RNA digestion products are recognized by TLR3 and 7 inside the phagolysosome. Parasite sensing by TLRs leads to synthesis of cytokines, chemokines and INFα/β, which in turn induce expression of NOS2. Oxidation of arginine by NOS2 generates nitric oxide which could kill the parasite either inside the phagolysosome or in the cytosol. Thus, downregulation of laminin (LAMA3) and overexpression of TLRs and NOS2, could potentially impact parasite growth in SYNs at different stages of the parasite intracellular life cycle.
Collectively, our studies show that increased TLRs signaling and resulting inflammatory responses may together play an important role in SYNs’ resistance to T. cruzi. Understanding the immunological landscape of the placenta, including the contribution of SYNs and how innate immunity is regulated, may contribute to the development of new therapeutics to reduce the risk of congenital CD.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: NCBI – PRJNA994533.
Author contributions
ES: Conceptualization, Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing. CC: Data curation, Formal analysis, Software, Writing – review & editing. AD: Conceptualization, Formal analysis, Supervision, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was supported by the intramural grants from the U.S. Food and Drug Administration to AD. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. The findings and conclusions in this article have not been formally disseminated by the Food and Drug Administration and should not be construed to represent any Agency determination and policy.
Acknowledgments
We wish to thank Luis Santana-Quintero (FDA High-performance Integrated Virtual Environment) for the transfer and analysis of the original RNA sequencing data, Sreenivas Gannavaram for valuable advice and suggestions, David Acosta for technical support, and FDA-CBER Facility for Biotechnology Resources staff for assisting with the RNA sequencing.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2023.1256385/full#supplementary-material
Footnotes
References
Ait Djebbara, S., McHeik, S., Percier, P., Segueni, N., Poncelet, A., and Truyens, C. (2023). The macrophage infectivity potentiator of Trypanosoma cruzi induces innate IFN-gamma and TNF-alpha production by human neonatal and adult blood cells through TLR2/1 and TLR4. Front. Immunol. 14:1180900. doi: 10.3389/fimmu.2023.1180900
Ander, S. E., Diamond, M. S., and Coyne, C. B. (2019). Immune responses at the maternal-fetal interface. Sci. Immunol. 4:eaat6114. doi: 10.1126/sciimmunol.aat6114
Ander, S. E., Rudzki, E. N., Arora, N., Sadovsky, Y., Coyne, C. B., and Boyle, J. P. (2018). Human placental Syncytiotrophoblasts restrict Toxoplasma gondii attachment and replication and respond to infection by producing immunomodulatory chemokines. MBio 9:e01678-17. doi: 10.1128/mBio.01678-17
Antinori, S., Galimberti, L., Bianco, R., Grande, R., Galli, M., and Corbellino, M. (2017). Chagas disease in Europe: a review for the internist in the globalized world. Eur. J. Intern. Med. 43, 6–15. doi: 10.1016/j.ejim.2017.05.001
Arora, N., Sadovsky, Y., Dermody, T. S., and Coyne, C. B. (2017). Microbial vertical transmission during human pregnancy. Cell Host Microbe 21, 561–567. doi: 10.1016/j.chom.2017.04.007
Avalos-Borges, E. E., Rios, L. E., Jimenez-Coello, M., Ortega-Pacheco, A., and Garg, N. J. (2022). Animal models of Trypanosoma cruzi congenital transmission. Pathogens 11:1172. doi: 10.3390/pathogens11101172
Bafica, A., Santiago, H. C., Goldszmid, R., Ropert, C., Gazzinelli, R. T., and Sher, A. (2006). Cutting edge: TLR9 and TLR2 signaling together account for MyD88-dependent control of parasitemia in Trypanosoma cruzi infection. J. Immunol. 177, 3515–3519. doi: 10.4049/jimmunol.177.6.3515
Barbosa, C. G., Gomez-Hernandez, C., da Silva, M. V., Rezende-Oliveira, K., Ferreira, P. T. M., de Oliveira, A. C. M., et al. (2022). Congenital transmission of Mexican strains of Trypanosoma cruzi TcIa: interaction between parasite and human placental explants. Parasitology 149, 418–426. doi: 10.1017/S0031182021002018
Barrila, J., Crabbe, A., Yang, J., Franco, K., Nydam, S. D., Forsyth, R. J., et al. (2018). Modeling host-pathogen interactions in the context of the microenvironment: three-dimensional cell culture comes of age. Infect. Immun. 86:e00282-18. doi: 10.1128/IAI.00282-18
Blaszkowska, J., and Goralska, K. (2014). Parasites and fungi as a threat for prenatal and postnatal human development. Ann. Parasitol. 60, 225–234.
Brigido, R., Tavares, P. C. B., Santos, M. A. D., Santos, J. G., Souza, M. A., Goulart, I. M. B., et al. (2017). Trypanosoma cruzi modulates gene expression of plasma membrane repair-related proteins. Acta Trop. 174, 153–157. doi: 10.1016/j.actatropica.2016.06.008
Caetano, B. C., Carmo, B. B., Melo, M. B., Cerny, A., dos Santos, S. L., Bartholomeu, D. C., et al. (2011). Requirement of UNC93B1 reveals a critical role for TLR7 in host resistance to primary infection with Trypanosoma cruzi. J. Immunol. 187, 1903–1911. doi: 10.4049/jimmunol.1003911
Cardoso, M. S., Reis-Cunha, J. L., and Bartholomeu, D. C. (2015). Evasion of the immune response by Trypanosoma cruzi during acute infection. Front. Immunol. 6:659. doi: 10.3389/fimmu.2015.00659
Carlier, Y., Altcheh, J., Angheben, A., Freilij, H., Luquetti, A. O., Schijman, A. G., et al. (2019). Congenital Chagas disease: updated recommendations for prevention, diagnosis, treatment, and follow-up of newborns and siblings, girls, women of childbearing age, and pregnant women. PLoS Negl. Trop. Dis. 13:e0007694. doi: 10.1371/journal.pntd.0007694
Carlier, Y., Schijman, A., and Kemmerling, U. (2020). Placenta, Trypanosoma cruzi, and congenital Chagas disease. Curr. Trop. Med. Rep. 7, 172–182. doi: 10.1007/s40475-020-00210-5
Castillo, C., Carrillo, I., Libisch, G., Juiz, N., Schijman, A., Robello, C., et al. (2018). Host-parasite interaction: changes in human placental gene expression induced by Trypanosoma cruzi. Parasit. Vectors 11:479. doi: 10.1186/s13071-018-2988-0
Castillo, C., Munoz, L., Carrillo, I., Liempi, A., Gallardo, C., Galanti, N., et al. (2017a). 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. 78:e12660. doi: 10.1111/aji.12660
Castillo, C., Munoz, L., Carrillo, I., Liempi, A., Medina, L., Galanti, N., et al. (2017b). Toll-like receptor-2 mediates local innate immune response against Trypanosoma cruzi in ex vivo infected human placental chorionic villi explants. Placenta 60, 40–46. doi: 10.1016/j.placenta.2017.10.005
Chen, S., Zhou, Y., Chen, Y., and Gu, J. (2018). Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34, i884–i890. doi: 10.1093/bioinformatics/bty560
Corry, J., Arora, N., Good, C. A., Sadovsky, Y., and Coyne, C. B. (2017). Organotypic models of type III interferon-mediated protection from Zika virus infections at the maternal-fetal interface. Proc. Natl. Acad. Sci. U. S. A. 114, 9433–9438. doi: 10.1073/pnas.1707513114
Coyne, C. B., Kim, K. S., and Bergelson, J. M. (2007). Poliovirus entry into human brain microvascular cells requires receptor-induced activation of SHP-2. EMBO J. 26, 4016–4028. doi: 10.1038/sj.emboj.7601831
Dabral, D., and van den Bogaart, G. (2021). The roles of phospholipase a(2) in phagocytes. Front. Cell Dev. Biol. 9:673502. doi: 10.3389/fcell.2021.673502
Diaz-Lujan, C., Triquell, M. F., Castillo, C., Hardisson, D., Kemmerling, U., and Fretes, R. E. (2016). Role of placental barrier integrity in infection by Trypanosoma cruzi. Acta Trop. 164, 360–368. doi: 10.1016/j.actatropica.2016.09.021
Duan, T., Du, Y., Xing, C., Wang, H. Y., and Wang, R. F. (2022). Toll-like receptor signaling and its role in cell-mediated immunity. Front. Immunol. 13:812774. doi: 10.3389/fimmu.2022.812774
Duaso, J., Rojo, G., Cabrera, G., Galanti, N., Bosco, C., Maya, J. D., et al. (2010). Trypanosoma cruzi induces tissue disorganization and destruction of chorionic villi in an ex vivo infection model of human placenta. Placenta 31, 705–711. doi: 10.1016/j.placenta.2010.05.007
Edwards, M. S., Stimpert, K. K., Bialek, S. R., and Montgomery, S. P. (2019). Evaluation and Management of Congenital Chagas Disease in the United States. J. Pediatric Infect. Dis. Soc. 8, 461–469. doi: 10.1093/jpids/piz018
Fernandes, M. C., and Andrews, N. W. (2012). Host cell invasion by Trypanosoma cruzi: a unique strategy that promotes persistence. FEMS Microbiol. Rev. 36, 734–747. doi: 10.1111/j.1574-6976.2012.00333.x
Fernandez-Aguilar, S., Lambot, M. A., Torrico, F., Alonso-Vega, C., Córdoba, M., Suarez, E., et al. (2005). Placental lesions in human Trypanosoma cruzi infection. Rev. Soc. Bras. Med. Trop. 38, 84–86.
Ferri, G., and Edreira, M. M. (2021). All roads Lead to cytosol: Trypanosoma cruzi multi-strategic approach to invasion. Front. Cell. Infect. Microbiol. 11:634793. doi: 10.3389/fcimb.2021.634793
Fitzgerald, K. A., and Kagan, J. C. (2020). Toll-like receptors and the control of immunity. Cells 180, 1044–1066. doi: 10.1016/j.cell.2020.02.041
Gravina, H. D., Antonelli, L., Gazzinelli, R. T., and Ropert, C. (2013). Differential use of TLR2 and TLR9 in the regulation of immune responses during the infection with Trypanosoma cruzi. PLoS One 8:e63100. doi: 10.1371/journal.pone.0063100
Griffiths, G., and Mayorga, L. (2007). Phagosome proteomes open the way to a better understanding of phagosome function. Genome Biol. 8:207. doi: 10.1186/gb-2007-8-3-207
Hoo, R., Nakimuli, A., and Vento-Tormo, R. (2020). Innate immune mechanisms to protect against infection at the human decidual-placental Interface. Front. Immunol. 11:2070. doi: 10.3389/fimmu.2020.02070
Imai, K., Misawa, K., Osa, M., Tarumoto, N., Sakai, J., Mikita, K., et al. (2019). Chagas disease: a report of 17 suspected cases in Japan, 2012-2017. Trop. Med. Health 47:38. doi: 10.1186/s41182-019-0168-3
Irish, A., Whitman, J. D., Clark, E. H., Marcus, R., and Bern, C. (2022). Updated estimates and mapping for prevalence of Chagas disease among adults, United States. Emerg. Infect. Dis. 28, 1313–1320. doi: 10.3201/eid2807.212221
Juiz, N. A., Torrejon, I., Burgos, M., Torres, A. M. F., Duffy, T., Cayo, N. M., et al. (2018). Alterations in placental gene expression of pregnant women with chronic Chagas disease. Am. J. Pathol. 188, 1345–1353. doi: 10.1016/j.ajpath.2018.02.011
Karki, R., and Kanneganti, T. D. (2021). The 'cytokine storm': molecular mechanisms and therapeutic prospects. Trends Immunol. 42, 681–705. doi: 10.1016/j.it.2021.06.001
Kemmerling, U., Osuna, A., Schijman, A. G., and Truyens, C. (2019). 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. 10:1854. doi: 10.3389/fmicb.2019.01854
Keogh, A., Ferguson, W., Nolder, D., Doyle, E., and McEntagart, N. (2021). Vertical transmission of Trypanosoma cruzi in a non-endemic country: histology of the infected placenta. Pediatr. Dev. Pathol. 24, 559–563. doi: 10.1177/10935266211012192
Kim, D., Paggi, J. M., Park, C., Bennett, C., and Salzberg, S. L. (2019). Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat. Biotechnol. 37, 907–915. doi: 10.1038/s41587-019-0201-4
Klein, M. D., Proano, A., Noazin, S., Sciaudone, M., Gilman, R. H., and Bowman, N. M. (2021). Risk factors for vertical transmission of Chagas disease: a systematic review and meta-analysis. Int. J. Infect. Dis. 105, 357–373. doi: 10.1016/j.ijid.2021.02.074
Kramer, A., Green, J., Pollard, J. Jr., and Tugendreich, S. (2014). Causal analysis approaches in ingenuity pathway analysis. Bioinformatics 30, 523–530. doi: 10.1093/bioinformatics/btt703
Liao, Y., Smyth, G. K., and Shi, W. (2014). featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30, 923–930. doi: 10.1093/bioinformatics/btt656
Liempi, A., Castillo, C., Carrillo, I., Munoz, L., Droguett, D., Galanti, N., et al. (2016). A local innate immune response against Trypanosoma cruzi in the human placenta: the epithelial turnover of the trophoblast. Microb. Pathog. 99, 123–129. doi: 10.1016/j.micpath.2016.08.022
Liempi, A., Castillo, C., Medina, L., Rojas, M., Maya, J. D., Parraguez, V. H., et al. (2019). Ex vivo infection of human placental explants with Trypanosoma cruzi and Toxoplasma gondii: differential activation of NF kappa B signaling pathways. Acta Trop. 199:105153. doi: 10.1016/j.actatropica.2019.105153
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262
Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15:550. doi: 10.1186/s13059-014-0550-8
Macaluso, G., Grippi, F., Di Bella, S., Blanda, V., Gucciardi, F., Torina, A., et al. (2023). A review on the immunological response against Trypanosoma cruzi. Pathogens 12. doi: 10.3390/pathogens12020282
Matthews, S., Tannis, A., Puchner, K. P., Bottazzi, M. E., Cafferata, M. L., Comande, D., et al. (2022). Estimation of the morbidity and mortality of congenital Chagas disease: a systematic review and meta-analysis. PLoS Negl. Trop. Dis. 16:e0010376. doi: 10.1371/journal.pntd.0010376
McConkey, C. A., Delorme-Axford, E., Nickerson, C. A., Kim, K. S., Sadovsky, Y., Boyle, J. P., et al. (2016). A three-dimensional culture system recapitulates placental syncytiotrophoblast development and microbial resistance. Sci. Adv. 2:e1501462. doi: 10.1126/sciadv.1501462
Megli, C., Morosky, S., Rajasundaram, D., and Coyne, C. B. (2021). Inflammasome signaling in human placental trophoblasts regulates immune defense against Listeria monocytogenes infection. J. Exp. Med. 218. doi: 10.1084/jem.20200649
Mezzano, L., Moran, J. P., Moreira-Espinoza, M. J., Triquell, M. F., Mezzano, J., Diaz-Lujan, C. M., et al. (2022). Chagas disease affects the human placental barrier's turnover dynamics during pregnancy. Mem. Inst. Oswaldo Cruz 117:e210304. doi: 10.1590/0074-02760210304
Motomura, K., Hara, M., Ito, I., Morita, H., and Matsumoto, K. (2023). Roles of human trophoblasts' pattern recognition receptors in host defense and pregnancy complications. J. Reprod. Immunol. 156:103811. doi: 10.1016/j.jri.2023.103811
Navarro, M., Reguero, L., Subira, C., Blazquez-Perez, A., and Requena-Mendez, A. (2022). Estimating Chagas disease prevalence and number of underdiagnosed, and undertreated individuals in Spain. Travel Med. Infect. Dis. 47:102284. doi: 10.1016/j.tmaid.2022.102284
Olson, N., and van der Vliet, A. (2011). Interactions between nitric oxide and hypoxia-inducible factor signaling pathways in inflammatory disease. Nitric Oxide 25, 125–137. doi: 10.1016/j.niox.2010.12.010
Pappas, G., Wilkinson, M. L., and Gow, A. J. (2023). Nitric oxide regulation of cellular metabolism: adaptive tuning of cellular energy. Nitric Oxide 131, 8–17. doi: 10.1016/j.niox.2022.11.006
Perez-Molina, J. A., and Molina, I. (2018). Chagas disease. Lancet 391, 82–94. doi: 10.1016/S0140-6736(17)31612-4
Perez-Zetune, V., Bialek, S. R., Montgomery, S. P., and Stillwaggon, E. (2020). Congenital Chagas disease in the United States: the effect of commercially priced Benznidazole on costs and benefits of maternal screening. Am. J. Trop. Med. Hyg. 102, 1086–1089. doi: 10.4269/ajtmh.20-0005
Perumal, N., White, B., Sanchez-Valdez, F., and Tarleton, R. L. (2023). cGAS-STING pathway activation during Trypanosoma cruzi infection leads to tissue-dependent parasite control. J. Immunol. ji2300373. doi: 10.4049/jimmunol.2300373 [Epub ahead of print].
Queiroga, T. B. D., Pereira, N. S., da Silva, D. D., Andrade, C. M., de Araujo Junior, R. F., Brito, C., et al. (2021). Virulence of Trypanosoma cruzi strains is related to the differential expression of innate immune receptors in the heart. Front. Cell. Infect. Microbiol. 11:696719. doi: 10.3389/fcimb.2021.696719
Rios, L., Campos, E. E., Menon, R., Zago, M. P., and Garg, N. J. (2020). Epidemiology and pathogenesis of maternal-fetal transmission of Trypanosoma cruzi and a case for vaccine development against congenital Chagas disease. Biochim. Biophys. Acta Mol. basis Dis. 1866:165591. doi: 10.1016/j.bbadis.2019.165591
Rosales, C., and Uribe-Querol, E. (2017). Phagocytosis: a fundamental process in immunity. Biomed. Res. Int. 2017, 9042851–9042818. doi: 10.1155/2017/9042851
Schmitz, V., Almeida, L. N., Svensjo, E., Monteiro, A. C., Kohl, J., and Scharfstein, J. (2014). C5a and bradykinin receptor cross-talk regulates innate and adaptive immunity in Trypanosoma cruzi infection. J. Immunol. 193, 3613–3623. doi: 10.4049/jimmunol.1302417
Semmes, E. C., and Coyne, C. B. (2022). Innate immune defenses at the maternal-fetal interface. Curr. Opin. Immunol. 74, 60–67. doi: 10.1016/j.coi.2021.10.007
Silberstein, E., Kim, K. S., Acosta, D., and Debrabant, A. (2021). Human placental trophoblasts are resistant to Trypanosoma cruzi infection in a 3D-culture model of the maternal-fetal Interface. Front. Microbiol. 12:626370. doi: 10.3389/fmicb.2021.626370
Silberstein, E., Serna, C., Fragoso, S. P., Nagarkatti, R., and Debrabant, A. (2018). A novel nanoluciferase-based system to monitor Trypanosoma cruzi infection in mice by bioluminescence imaging. PLoS One 13:e0195879. doi: 10.1371/journal.pone.0195879
Singh, P., and Ali, S. A. (2022). Multifunctional role of S100 protein family in the immune system: an update. Cells 11:2274. doi: 10.3390/cells11152274
Sreejit, G., Flynn, M. C., Patil, M., Krishnamurthy, P., Murphy, A. J., and Nagareddy, P. R. (2020). S100 family proteins in inflammation and beyond. Adv. Clin. Chem. 98, 173–231. doi: 10.1016/bs.acc.2020.02.006
Teixeira, A. A., de Vasconcelos Vde, C., Colli, W., Alves, M. J., and Giordano, R. J. (2015). Trypanosoma cruzi binds to cytokeratin through conserved peptide motifs found in the laminin-G-like domain of the gp85/trans-sialidase proteins. PLoS Negl. Trop. Dis. 9:e0004099. doi: 10.1371/journal.pntd.0004099
Torres-Vargas, J., Jimenez-Coello, M., Guzman-Marin, E., Acosta-Viana, K. Y., Yadon, Z. E., Gutierrez-Blanco, E., et al. (2018). Quantitative and histological assessment of maternal-fetal transmission of Trypanosoma cruzi in guinea pigs: an experimental model of congenital Chagas disease. PLoS Negl. Trop. Dis. 12:e0006222. doi: 10.1371/journal.pntd.0006222
Triquell, M. F., Diaz-Lujan, C., Romanini, M. C., Ramirez, J. C., Paglini-Oliva, P., Schijman, A. G., et al. (2018). Nitric oxide synthase and oxidative-nitrosative stress play a key role in placental infection by Trypanosoma cruzi. Am. J. Reprod. Immunol. 80:e12852. doi: 10.1111/aji.12852
Vaena de Avalos, S., Blader, I. J., Fisher, M., Boothroyd, J. C., and Burleigh, B. A. (2002). Immediate/early response to Trypanosoma cruzi infection involves minimal modulation of host cell transcription. J. Biol. Chem. 277, 639–644. doi: 10.1074/jbc.M109037200
Vekemans, J., Truyens, C., Torrico, F., Solano, M., Torrico, M. C., Rodriguez, P., et al. (2000). Maternal Trypanosoma cruzi infection upregulates capacity of uninfected neonate cells to produce pro- and anti-inflammatory cytokines. Infect. Immun. 68, 5430–5434. doi: 10.1128/IAI.68.9.5430-5434.2000
Venturini, G., Alvim, J. M., Padilha, K., Toepfer, C. N., Gorham, J. M., Wasson, L. K., et al. (2023). Cardiomyocyte infection by Trypanosoma cruzi promotes innate immune response and glycolysis activation. Front. Cell. Infect. Microbiol. 13:1098457. doi: 10.3389/fcimb.2023.1098457
Xia, C., Braunstein, Z., Toomey, A. C., Zhong, J., and Rao, X. (2017). S100 proteins as an important regulator of macrophage inflammation. Front. Immunol. 8:1908. doi: 10.3389/fimmu.2017.01908
Keywords: congenital Chagas disease, Trypanosoma cruzi, 3D-trophoblasts, resistance to infection, transcriptome analysis, Toll-like receptors
Citation: Silberstein E, Chung CC and Debrabant A (2023) The transcriptome landscape of 3D-cultured placental trophoblasts reveals activation of TLR2 and TLR3/7 in response to low Trypanosoma cruzi parasite exposure. Front. Microbiol. 14:1256385. doi: 10.3389/fmicb.2023.1256385
Edited by:
Carlos Robello, Universidad de la República, UruguayReviewed by:
Yves Carlier, Université libre de Bruxelles, BelgiumUlrike Kemmerling, University of Chile, Chile
Copyright © 2023 Silberstein, Chung and Debrabant. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Alain Debrabant, YWxhaW4uZGVicmFiYW50QGZkYS5oaHMuZ292