Leilane Oliveira Gonçalves1,2
Andrés F. Vallejo Pulido3
Fernando Augusto Siqueira Mathias2Alexandre Estevão Silvério Enes1,2Maria Gabriela Reis Carvalho2
Daniela de Melo Resende4Marta E. Polak3*†
Jeronimo C. Ruiz2*†- 1Programa de Pós-graduação em Biologia Computacional e Sistemas, Instituto Oswaldo Cruz, Fiocruz, Rio de Janeiro, Brazil
- 2Grupo Informática de Biossistemas, Instituto René Rachou, Fiocruz Minas, Belo Horizonte, Brazil
- 3Systems Immunology Group, Clinical and Experimental Sciences, Faculty of Medicine, University of Southampton, Southampton, United Kingdom
- 4Grupo Genômica Funcional de Parasitos, Instituto René Rachou, Fiocruz Minas, Belo Horizonte, Brazil
Leishmania amazonensis and Leishmania major are the causative agents of cutaneous and mucocutaneous diseases. The infections‘ outcome depends on host–parasite interactions and Th1/Th2 response, and in cutaneous form, regulation of Th17 cytokines has been reported to maintain inflammation in lesions. Despite that, the Th17 regulatory scenario remains unclear. With the aim to gain a better understanding of the transcription factors (TFs) and genes involved in Th17 induction, in this study, the role of inducing factors of the Th17 pathway in Leishmania–macrophage infection was addressed through computational modeling of gene regulatory networks (GRNs). The Th17 GRN modeling integrated experimentally validated data available in the literature and gene expression data from a time-series RNA-seq experiment (4, 24, 48, and 72 h post-infection). The generated model comprises a total of 10 TFs, 22 coding genes, and 16 cytokines related to the Th17 immune modulation. Addressing the Th17 induction in infected and uninfected macrophages, an increase of 2- to 3-fold in 4–24 h was observed in the former. However, there was a decrease in basal levels at 48–72 h for both groups. In order to evaluate the possible outcomes triggered by GRN component modulation in the Th17 pathway. The generated GRN models promoted an integrative and dynamic view of Leishmania–macrophage interaction over time that extends beyond the analysis of single-gene expression.
Introduction
Leishmaniasis is a complex of infectious diseases caused by different species of the protozoan parasite Leishmania (Kinetoplastida and Trypanosomatidae) that represent one of the major public health problems in developing countries, with 12 million cases and an incidence of 0.7–1.0 million new cases annually, according to the WHO (WHO, 2020). Depending on the Leishmania species involved, the disease can comprise two main clinical forms: visceral leishmaniasis or cutaneous leishmaniasis (Burza et al., 2018). Between both clinical forms, the former has the highest prevalence and incidence and is caused by different Leishmania species, including, but not restricted to, Leishmania major, Leishmania amazonensis, Leishmania braziliensis, Leishmania mexicana, and Leishmania peruviana (Martinez and Petersen, 2014). In humans, macrophages are the final host cells for Leishmania parasites, playing a dual role in the progression of the infection, responsible for control or exacerbation of the immune response that culminates in different clinical manifestations (Tomiotto-Pellissier et al., 2018). Infections by Leishmania sp. promote different types of immune responses, creating a microenvironment of cytokines and chemokines, which is a determinant of the establishment or not of the disease. In this context, T helper cell responses have an important role, through different T cell subsets such as Th1, Th2, Th17, and Tregs (Gonçalves-de-Albuquerque et al., 2017). The Th1 response drives an increase of CD4+ T cells with the production of interferon-γ, interleukin (IL) 2 and IL-12 that stimulate macrophages to fight against infected cells (Isnard et al., 2012; Mutiso et al., 2013). However, if the immune response is Th2 type, IL-4 and IL-10 are produced, which act in inhibiting the activation of macrophages, consequently promoting the aggravation of the disease (Isnard et al., 2012; Liu and Uzonna, 2012; Mutiso et al., 2013; de Paiva et al., 2015). Lastly, Th17 cells may play an essential role in protecting against certain extracellular pathogens and have been demonstrated to influence the balance between inflammatory and anti-inflammatory cytokines, which must be orchestrated in the course of infection to guide a successful effector response to the intracellular Leishmania protozoa (Gonçalves-de-Albuquerque et al., 2017). Particular interest in the Th17 type remains focused on its main cytokine, IL-17. Although it is also produced by other cells including CD8+ T cells and neutrophils, IL-17 is mostly produced by Th17 cells. Thus, the effects of Th17/IL-17 are still unclear, as disease-promoting and protective responses have both been attributed to their influence. In this scenario, the Th17 cell activities in the context of the different clinical forms of leishmaniasis, interaction with regulatory cytokines, and host particularities remain to be explored (Anderson et al., 2009, 17; Liu and Uzonna, 2012, 17; Gaffen et al., 2014; Bartlett and Million, 2015; Gonçalves-de-Albuquerque et al., 2017, 17; Tomiotto-Pellissier et al., 2018, 17).
The response of infection, in macrophages, occurs due to the activation of different transcription factors (TFs), such as STAT1, IRF, NF-κB, and STAT6 (Bichiou et al., 2021; Murray, 2017).
Although many reviews have discussed the cellular immune responses against the different forms of leishmaniasis (Pitta et al., 2009; Hu et al., 2010; Gonzalez-Lombana et al., 2013; Gonçalves-de-Albuquerque et al., 2017; Banerjee et al., 2018) and the roles of Th17/IL-17 in infectious and non-infectious diseases, their role in the immunopathogenesis of leishmaniasis remains unclear (Gonçalves-de-Albuquerque et al., 2017).
Some studies (Vargas-Inchaustegui et al., 2008; Boaventura et al., 2010, Castilho et al., 2010; Katara et al., 2012; Gonzalez-Lombana et al., 2013; Katara et al., 2013; Bartlett and Million, 2015; Banerjee et al., 2018) have suggested that high levels of IL-17 may play an important pro-inflammatory role in leishmaniasis development and as such should be considered an important target for the disease immunotherapy.
An alternative to better understanding the role of these elements of the immune response in the immunopathogenesis of leishmaniasis would be through the use of gene regulatory networks (GRNs). GRNs can be defined as mathematical and computational models capable of describing the logic underlying the regulatory occurrences among interacting genes while a specific cell program is operating. The fundamental challenge that GRNs address is the translation between the diverse information patterns presented by different types of biological macromolecules into biologically significant knowledge. This knowledge is directly associated with the time-dynamic processes necessarily underlying any kind of observable biological behavior (Bansal et al., 2006; Anitha et al., 2016; de Luis Balaguer et al., 2017; Banf and Rhee, 2017; de Luis Balaguer et al., 2017; Peter and Davidson, 2017; Delgado and Gómez-Vela, 2019). A key element in the modeling process of GRNs is the identification of molecular regulators associated with the activation and control of gene expression, to which TFs are considered essential. They modulate the rate of expression of genes by inhibiting or enhancing the gene’s transcription rates and, in addition, interact with regulatory elements to generate specific gene expression patterns. Additionally, TFs are themselves the products of gene expression, serving as a source of feedback for the regulation of genetic networks (Steele et al., 2009; Chandrasekaran et al., 2011; Ghedira et al., 2011; Sinha et al., 2020).
Different approaches used for modeling GRNs are described in the literature (Miyagi et al., 2002; Hardy and Robillard, 2004; Anitha et al., 2016; Livigni et al., 2016; Pennisi et al., 2016; Singh et al., 2020). These approaches include differential equations, stochastic simulations, flux balance analysis, graphical models, and Petri nets (PNs). PNs are an in silico mathematical approach that provides a visual-aided network-oriented modeling process based on nodes and edges to represent places and transitions, with visual feedback affording interpretation by a broad audience and can be applied to describing distributed systems, including biological ones (Hardy and Robillard, 2004; Livigni et al., 2016). In the context of PNs, places represent molecules (proteins, mRNA, and complexes), transitions represent possible processes (binding, phosphorylation, and inhibition, including reversibility, among others), and tokens represent concentrations.
On this basis, we decided to use GRNs to evaluate the dynamic of the infection profile correlated with the modulation of Th17 immune response in cutaneous leishmaniasis through the utilization of data from transcriptome sequencing of macrophages infected with L. major and L. amazonensis associated with data experimentally validated in the literature.
Materials and Methods
Dataset Selection
A total of 66 samples from the study SRP062278 were performed by Fernandes et al., 2016 (Fernandes et al., 2016). The data were downloaded from the Sequence Read Archive (SRA) database using SRA prefetch and converted to fastq files using fastq-dump, available in SRA toolkit version 2.10.2. All mRNA paired-end libraries with 100 bp were sequenced using Illumina Hiseq 1500 platform. The study SRP062278 comprises a time-series experiment with four different time points (4, 24, 48, and 72 h post-infection (hpi) of human macrophages infected with L. amazonensis (IFLA/BR/67/PH8), human macrophages infected with L. major (clone VI, MHOM/IL/80/Friedlin), uninfected human macrophages, and macrophages containing latex beads.
Data Preprocessing
The sequence quality of the downloaded data was assessed using the FastQC program (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/) and visualized using MultiQC reporting tool (Ewels et al., 2016).
As references, the human transcriptome data in FASTA and GTF format was obtained from the Ensembl database (ftp://ftp.ensembl.org/pub/release-99/), release 99. FASTA and GFF files with transcript sequences and annotations from L. amazonensis MHOM/BR/71973/M2269 and L. major Friedlin strain were downloaded from TriTrypDB (https://tritrypdb.org/tritrypdb/), release 46 (de Azevedo et al., 2010). The transcripts sequences of L. amazonensis MHOM/BR/71973/M2269 were used as a reference due to the absence of sequencing information of L. amazonensis IFLA/BR/67/PH8 in public databases.
Transcriptome Quantification and Differential Gene Expression Analyses
Kallisto (version 0.46.1) (Bray et al., 2016) was used to generate the transcriptome indexes for human and Leishmania sp. data. All samples were pseudoaligned against the references and quantified using Kallisto quant, using bias correction parameter and 100 as the bootstrap number. The differential gene expression analysis was performed in R using Sleuth (version 0.30.0) (Pimentel et al., 2017), Limma (version 3.42.1) (Law et al., 2018), and edgeR (version 3.28.0) (Robinson et al., 2010) packages. In order to identify differentially expressed (DE) genes, an adjusted p-value of less than 0.05 was set as a threshold to define the significance.
Network Analysis
The network analysis was performed using Graphia Pro, available in BioLayout Express3D (Theocharidis et al., 2009), where a gene correlation matrix was calculated for all DE human genes fulfilling the criteria mentioned above and considering that each column of data represents one sample and each row of data represents one gene.
A network graph of the data was generated with a Pearson’s correlation coefficient of r ≥ 0.85 using TMM normalization. The network graph was then clustered into groups of genes sharing similar expression profiles using the MCL (Li, 2003) algorithm with an inflation value of 1.7. All the clusters were analyzed considering the differences in the expression profiles between the conditions and functionally annotated using ToppGene tool (Chen et al., 2009) and David (Dennis et al., 2003) software.
With the use of the information related to the clusters functional annotation, the most relevant pathways in Leishmania infection were evaluated, and the Th17 pathway was selected (Supplementary Table 1).
Model Assembly
In order to identify components for the Th17 GRN, a systematic search in PubMed was performed. Of the returned papers, 327 were identified describing the IL-17 pathway and its TFs. The experimental findings within each listed reference paper were analyzed and categorized as the input information for the network into a) partner 1 (a gene that acts by somehow regulating a second gene); b) partner 2 (a gene that undergoes some type of induction/inhibition by partner 1); c) TFs that interact with partners 1 and 2; d) cell type in which this interaction was observed; e) type of interaction; and f) induced biological process and outcome of the interaction between partners 1 and 2 (activation of other genes, induction of cytokines or chemokines). To be included in the network diagram, all the identified interactions had to be confirmed by at least two different publications. This information was categorized as the network components: input node, transmission node, an output node, and mode of interaction (Supplementary Table 2).
In the model, edges, connecting nodes, and transitions determine the direction of information flow, representing the progress of the biological process over time (time-blocks). Each infection time in the experimental design represents a time-block that was used for the simulation. Thus, the time-blocks were defined as a) 4 hpi (time-block 0–25), b) 24 hpi (time-block 26–50), c) 48 hpi (time-block 51–75), and d) 72 hpi (time-block 76–100).
IL-17 Gene Regulatory Network Model Parametrization and In Silico Simulations
The network diagram was constructed using yEd Graph Editor following the mEPN notation, as described by Livigni et al. (2018), and saved as GraphML files (Livigni et al., 2016; O’Hara et al., 2016).
Prior to performing a simulation, a model was parameterized through the places representing the TFs at the beginning of the network. For this parametrization, the normalized read count values obtained from the RNA-seq data were used to parametrize the input tokens that correspond to TFs identified in the literature review that interact with one or more genes. This resulted in the modeling of three networks: the first considering the uninfected macrophages, the second considering macrophages infected with L. amazonensis, and the third as macrophages infected with L. major.
GraphML files were loaded into the BioLayout Express 3D, and a parser translated the diagrams drawn using the mEPN scheme. Once the files were imported, the conditions for stochastic PN (SPN) simulation were set (i.e., mode of stochasticity, number of runs, and defined time-blocks).
Each time point in the experimental design represents an interval of the time-blocks used for the simulation. All the measured read count averages for the TFs identified, at 4 hpi (0–25 time-block), 24 hpi (26–50 time-block), 48 hpi (51–75 time-block), and 72 hpi (76–100 time-block), were converted into parameterization values for each entry node in the GRN. Simulations were conducted in Graphia Pro using 100 time-blocks, 500 runs, uniform distribution in SPNs, and distribution and consumptive transition in SPN transition type.
During the simulation, some inhibitions and blocking edges were also made in the constructed models in order to observe the impact of this alteration on the outcome (Th17 response).
Statistical Analyses
In order to compare the TFs in each time point for each type of infected macrophage, the statistical analysis was performed using mixed ANOVA, with a general linear model (α = 0.05). Sphericity was evaluated using Mauchly’s test (α = 0.05) and/or cases where it was not verified; the Greenhouse–Geisser correction was performed to verify the effect of the independent variables and their interactions on the expression value of each of the factors individually. To locate the observed differences, a post-hoc analysis with Bonferroni correction was conducted, also with a significance level of 0.05.
Results
Transcriptome Quantification and Differential Gene Expression Analyses
Approximately 90% of reads were mapped against the reference transcriptomes. The samples SRR2163291 (macrophage containing beads at 48 hpi) and SRR2155160 (macrophage infected with L. major at 24 hpi) were identified as outliers and removed from the analyses due to the observation of distinct profiles in the principal component analysis.
Considering the applied cutoff to define significance, for the macrophage infected with L. major dataset when compared to uninfected macrophages, for each time-point (24, 224, 48, and 72 hpi), a total of 3,082, 1,389, 1,193, and 642 DE genes were identified, respectively.
For the dataset comprising the macrophages infected with L. amazonensis, when compared to uninfected macrophages, a total of 2,452, 1,332, 301, and 327 DE genes were found for each time-point (24, 224, 48, and 72 hpi, respectively). The results are depicted in Table 1.
Table 1 Differentially expressed genes.
The most expressed genes in infected macrophages were functionally annotated. For the biological process category, the main ontologies identified in both datasets (macrophages infected with L. major and macrophages infected with L. amazonensis) were a cellular response to the stimulus caused by DNA damage (GO:0006974), DNA repair (GO:0006281), DNA metabolic process (GO:0006259), repair (GO:0006260), and modifications in lysines (GO:0018205). Considering the category of molecular function, the main ontologies identified were related to insertion or deletion of DNA bases (GO:0032135), DNA binding (GO:0003684), and DNA-dependent ATPase activity (GO:0008094). For the cellular component category, ontologies related to chromatin (GO:0000785) and centrosome (GO:0005813) were identified.
The most overrepresented in the dataset was grouped based on semantics, and it was observed that most of the functions found were related to changes in the cell cycle, stress response, DNA repair mechanisms, protein modification processes, and diversification of immunological receptors.
Differentially Expressed Gene Clustering and Annotation
A network graph of the data was generated with a Pearson’s correlation coefficient of r ≥ 0.85 using TMM normalization in the read count table and clusterized using the MCL algorithm.
As result, a total of 57 clusters were obtained, but in only 30 clusters was it possible to identify significant alterations in the expression of a set of genes in macrophages infected with L. major and L. amazonensis when compared to uninfected macrophages (Supplementary Table 1).
The genes comprising the selected clusters were functionally annotated, which allowed the identification of many immunological pathways such as tumor necrosis factor (TNF) signaling pathway, NF-kappa B signaling pathway, cytokine–cytokine receptor interaction, Chagas disease (American trypanosomiasis), MAPK signaling pathway, toll-like receptor signaling pathway, IL-10 signaling pathway, and TH17/IL-17 signaling pathway.
Considering the pathways found in the annotation, a literature review was conducted in order to identify the most relevant pathways for Leishmania infection. In the literature review, using the search terms depicted in Supplementary Table 2, a total of 22 interactions were found considering the TFs, and a total of 73 interactions considering the genes in the IL-17 pathway were identified in at least 2 papers to be considered for the network model. All this information was categorized as network components, receiving the nomenclature of a) input node (for TFs) (n = 12); b) transmission node: genes that are activated during the course of infection (n = 22); c) output node: genes that are activated for a given immune response (n = 16); and d) mode of interaction: how these genes interact, whether through activation, binding, inhibition of phosphorylation.
Model Assembly and Parametrization
All the information categorized in the literature review was used to model the IL-17 network diagram. The components (genes) in the IL-17 GRN are represented by rectangles connected by arrows, which indicate the molecular interactions found during the literature review step. System start nodes (TFs) are represented by colored rectangles in black. The GRN output signal (Th17 response) is represented in the rectangle on the right side of the diagram (Figure 1). We observed that the flow of information starts from the activation of TFs and follows all the models of the genes and the interactions between them. All these interactions converge for the activation of cytokines, chemokines, and other ILs, such as IL-6 and IL-12.
Figure 1 Diagram of the constructed network. The network diagram was constructed using the yEd Graph Editor according to the mEPN notation, described by Livigni et al. (2018), using the information obtained from the search carried out in the literature. The components (genes) in the IL-17 gene regulatory network (GRN) are represented by rectangles connected by arrows that indicate the molecular interactions found during the literature review step. System start nodes (transcription factors) are represented by colored rectangles in black. The GRN output signal is represented in the rectangle on the right side of the diagram.
As a result of model parameterization, we identified the TFs involved in the IL-17 pathway and calculated the average of expression counts for each time-block evaluated (Table 2).
Table 2 Read count average for each transcription factor in IL-17 pathway for each time-block defined interval.
After the establishment of parameterization, we performed the network simulations. Under normal conditions, as can be seen in Figure 2, the signal intensity is about 2–3 times lower than when compared to samples of macrophages infected by L. major and L. amazonensis.
Figure 2 Accumulation token result of the initial network simulation. Under normal conditions, as can be seen in the yellow lines (uninfected macrophages), the final response signal intensity is about 2–3 times lower than when compared to samples of Leishmania major- and Leishmania amazonensis-infected macrophages (represented by the red and blue lines, respectively). To be generated, the initial information flow starts with the activation of transcription factors (input nodes). These factors interact with a set of genes (transmission node) through processes (mode of interaction). These processes are responsible for intensifying or reducing the intensity of the signal generated in response.
After the network parameterization and simulation, inhibition points on the edges between certain interactions were inserted in order to observe how the signal intensity changed based on the modifications in the model (Figure 3).
Figure 3 Examples of changes inserted into the network. Examples of representation of inhibitory and/or blocking changes that were inserted into the modeled network.
For the transcription factor STAT1 (signal transduction and activator of transcription 1), an edge change was inserted, which reduces the intensity in which phosphorylation is performed on the target gene for its activation, while in the STAT3 factor, a blocking edge was inserted, which completely inhibits the signal generated by this TF; i.e., it does not phosphorylate the target gene.
In addition to these two modifications in the TFs STAT1 and STAT3, partial inhibitory modifications were also inserted in CEBPβ (CEBPB-CCAAT/enhancer-binding protein beta) and NF-kβ (nuclear factor kappa-light-chain-enhancer of activated B cells). These modifications caused differences in the intensity of the signal generated (Figure 4), although they were not sufficient for a complete blockade of the Th17 response induction in infected macrophages when compared to uninfected macrophages.
Figure 4 Result of the network simulation with changes in transcription factors. Under normal conditions (uninfected macrophages), as can be seen in yellow lines, the intensity of the final response signal is lower than when compared to infected macrophages with Leishmania major and Leishmania amazonensis represented in red and blue lines, respectively.
After the inhibition of selected TFs, we noticed a decrease in the expression of CxCL2 and CCL7 in macrophages infected with L. major. For L. amazonensis, a decrease in CxCL8 and CXCL9 expression levels was also observed. The addition of inhibitory points in NF-kβ also causes a decrease in the potential to induce IL-12 and IL-23 levels. However, even with the inhibition of TFs, the gene expression related to the IL-17 pathway on infected macrophages did not equalize with non-infected macrophages.
For the analysis considering the TF comparisons between the time points (4, 24, 48, and 72 hpi) and across the macrophage’s type (uninfected, infected with L. major, and infected with L. amazonensis), Figure 5 shows the descriptive analysis of the dynamics of expression of the factors in the four time points.
Figure 5 Descriptive analysis of expression dynamics for the transcription factors in the four time points post-infection. The time points 1, 2, 3, and 4 represent each hour post-infection in the kinetics, respectively, 4 h post-infection, 24 h post-infection, 48 h post-infection, and 72 h post-infection. The blue lines represent uninfected macrophages, the red lines represent macrophages infected with Leishmania major, and the green lines represent macrophages infected with Leishmania amazonensis.
Significant differences (p < 0.05) in the expression of only STAT1, STAT4, TRAF6, and NFKB1 factors were observed for the different types of macrophages, at different times of infection. The expression of STAT 1 and STAT 4 was more associated in the early points of the infection (4–24 hpi) in macrophages infected with L. amazonensis and macrophages infected with L. major, respectively. An increase of expression for NF-κB was also observed in the time of 4 hpi for all infected macrophages when compared to uninfected macrophages.
Discussion
Leishmania parasites are capable of modulating the immunological response and many signaling pathways on host cells in order to promote survival and infection (Mandlik et al., 2016). The key points for a parasite infection are related to the ability to escape from the macrophage leishmanicidal mechanisms and modulate the immune system in order to survive (Ovalle-Bracho et al., 2015; Akhoundi et al., 2017). It is known that all these cellular mechanisms comprise a set of signaling events that result in the modulation of TFs, which lead to different gene expression patterns (Contreras et al., 2010; Grünebast and Clos, 2020; Lecoeur et al., 2020; Bichiou et al., 2021). Advances in transcriptome analysis techniques have enabled the identification and characterization of pathways involved in the infection process in both parasites and hosts. These approaches have been lauded due to their advanced effectiveness (Wang et al., 2009; Veras and Bezerra de Menezes, 2016).
In this study, the use of RNA-seq data from infected macrophages with L. major and L amazonensis allowed the identification of functions related to changes in cell cycle, stress response, DNA repair mechanisms, protein modification processes, and diversification of immunological receptors. These findings could be associated with alterations to macrophage suppression, favoring the parasite’s survival due to the capacity to alter signaling pathways of the macrophage in order to establish the infection (Ghedira et al., 2011; Liu and Uzonna, 2012; Ovalle-Bracho et al., 2015; Fernandes et al., 2016). Furthermore, these alterations in macrophage biological processes could differ according to the infective species (Liu and Uzonna, 2012; Ovalle-Bracho et al., 2015). Similar to Zhang and coworkers, the majority of genes found as DE was identified at initial infection time points (4 and 24 hpi) (Zhang et al., 2010).
Considering the immune response in cutaneous leishmaniasis, the presence of Th17 cells is related to the maintenance of the in situ inflammatory process due to the Th1/Th2 responses in skin lesion ulceration (Gonzalez et al., 2020). Our results showed that genes involved in the Th17 pathway are overexpressed on macrophages infected with both Leishmania species. Some studies have shown that one of the functions of the Th17 response is to kill pathogens via IL-17 production. It has also been reported that this pro-inflammatory cytokine could induce the expression of IL-6, IL-8, GM-CSF, and G-CSF (Katara et al., 2012; Katara et al., 2013; Gonçalves-de-Albuquerque et al., 2017). However, the role of Th17 in leishmaniasis is still unclear because, in CL, it is related to exacerbation of the lesions and disease pathology (Banerjee et al., 2016; Gonzalez et al., 2020).
A strategy to help clarify these questions could be the evaluation of the transcriptional profiles involved in this response. The modulation of the transcriptional profile in macrophages after the invasion by a pathogen is dynamic and changes along the course of the infection. At the initial stages, we observed the expression of genes related to immunological regulation; among these, genes that encode transcription factor STAT2, IL-2, and CXCL12 were identified as overexpressed until 24 hpi. These molecules are important for effective inflammatory response during infection, being responsible for inducing T-cell recruitment, proliferation, regulation, and production of cytokines such as type I interferons (IFN-α and IFN-β) (Steen and Gamero, 2012; Ross and Cantrell, 2018; Pol et al., 2019). Similar to other studies, an increase in the expression of CXCL10 and CCL2 was also observed (Rabhi et al., 2012).
It is known that IFN-α and IFN-β are related to the expression of nitric oxide synthase type 2 (NOS2) and production of NO by macrophages, as other studies have shown the importance of this pathway against intracellular infections (Diefenbach et al., 1998; Xin et al., 2008; Cummings et al., 2010). In this context, we hypothesize that the parasite starts the modulation of these genes until the second time-point (24 hpi), resulting in the reduction at third and fourth time-points in order to create a favorable microenvironment to establish the infection, neutralizing part of the macrophage defense machinery. In murine models, an infection caused by L. amazonensis showed that amastigotes induce a reduction of STAT-2 phosphorylation and an increase of degradation through parasite proteases (Xin et al., 2008).
Advances in transcriptome approaches allowed the identification of pathways related to the infective mechanisms in parasites and hosts. On the other hand, many of these approaches focus mainly on single genes, not in a global view of the pathogenesis process (Menezes et al., 2013; Castiglione et al., 2016; Lee et al., 2018; Veras et al., 2018). Different analytical tools are available for studying and analyzing next-generation sequencing (NGS) data; one of them was the GRN approach. The construction of these networks offers a measurement of relationships due to the reconstruction of the pathway (Huang et al., 2009; Kuleshov et al., 2016; Han et al., 2017; Jackson et al., 2020). In the network model, the addition of modifications in the TFs STAT1 and STAT3, CEBPβ (CEBPB-CCAAT/enhancer-binding protein beta), and NF-kB (nuclear factor kappa-light-chain-enhancer of activated B cells) caused differences in the intensity of the signal generated. However, they were not sufficient for a complete blockade of the Th17 response induction in infected macrophages when compared to uninfected macrophages. Approximately 30% of inhibition of the full pathway was achieved when compared to uninfected macrophages. Nonetheless, after the inhibition of selected TFs, we observed that they induced a decrease in the expression of CxCL2 and CCL7 in macrophages infected with L. major. For L. amazonensis, a decrease in CxCL8 and CXCL9 expression levels was also observed. The addition of inhibitory points in NF-kB also causes a decrease in the potential to induce IL-12 and IL-23 levels. Despite this, it is well known that these ILs were associated with a protective role against leishmaniasis (Dayakar et al., 2019). In the predictive model, this decrease could be related to missing genes in the literature review for the model construction. However, even with the inhibition of TFs, the gene expression related to the IL-17 pathway on infected macrophages did not equalize with that on non-infected macrophages.
Conclusion
The host–parasite interaction of leishmaniasis has enormous relevance in the disease outcome with many factors involved, as gene expression triggered by this interaction until cell recruitment and activation was due to cytokine and chemokine production and signalization. In this context, GRN is a suitable tool for managing this ocean of factors, specific points for evaluation, and attribute values in a foreseen scenario of dynamic interactions. The GRN computational modeling built in this study was able to predict how the signal (modulation of expression promoted by the alterations) flows through the network. Furthermore, our study was able to follow across determined kinetics the gene expression profile of macrophages infected with two different Leishmania spp. Specifically, we were able to predict the following: a) a set of TFs that can cause changes in the gene expression profile; b) a set of genes ranked as crucial for the immune modulation of the pathway (increase or decrease of Th17 final response); and c) the IL-17 pathway model for uninfected and infected macrophages, considering the time course of the infection. These promising results encourage us to keep looking for key factors present in other crucial pathways that we believe are involved in the resistance and susceptibility of macrophage infection by leishmaniasis.
Data Availability Statement
Publicly available datasets were analyzed in this study. These data can be found here: SRA accession number: SRP062278.
Author Contributions
LG acquired, analyzed, and interpreted the data and wrote the paper. AP acquired, analyzed, and interpreted the data. FM interpreted the data and wrote the paper. AE interpreted the data. MC interpreted the data, wrote the paper, and revised it critically for intellectual content. DR analyzed and interpreted the data. MP interpreted the data, wrote the paper, and revised critically for intellectual content. JR interpreted the data, wrote the paper, and revised it critically for intellectual content. All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
Funding
JR was funded by CNPq project number 310104/2018- and, INOVA-Fiocruz, Grants VPPIS-001-FIO-18-12 and VPPCB-005-FIO-20-2-42. MC was funded by INOVA-Fiocruz, Grant VPPIS-001-FIO-18-8.
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/fcimb.2022.826523/full#supplementary-material
References
Akhoundi, M., Downing, T., Votýpka, J., Kuhls, K., Lukeš, J., Cannet, A., et al. (2017). Leishmania Infections: Molecular Targets and Diagnosis. Mol. Aspects Med. 57, 1–29. doi: 10.1016/j.mam.2016.11.012
Anderson, C. F., Stumhofer, J. S., Hunter, C. A., Sacks, D. (2009). IL-27 Regulates IL-10 and IL-17 From CD4+ Cells in Nonhealing Leishmania Major Infection. J. Immunol. 183, 4619–4627. doi: 10.4049/jimmunol.0804024
Anitha, P., Anbarasu, A., Ramaiah, S. (2016). Gene Network Analysis Reveals the Association of Important Functional Partners Involved in Antibiotic Resistance: A Report on an Important Pathogenic Bacterium Staphylococcus Aureus. Gene 575, 253–263. doi: 10.1016/j.gene.2015.08.068
Banerjee, A., Bhattacharya, P., Dagur, P. K., Karmakar, S., Ismail, N., Joshi, A. B., et al. (2018). Live Attenuated Leishmania Donovani Centrin Gene–Deleted Parasites Induce IL-23–Dependent IL-17–Protective Immune Response Against Visceral Leishmaniasis in a Murine Model. J. Immunol. 200, 163–176. doi: 10.4049/jimmunol.1700674
Banerjee, A., Bhattacharya, P., Joshi, A. B., Ismail, N., Dey, R. (2016). And Nakhasi, HRole of Pro-Inflammatory Cytokine IL-17 in Leishmania Pathogenesis and in Protective Immunity by Leishmania Vaccines. L. Cell. Immunol. 309, 37–41. doi: 10.1016/j.cellimm.2016.07.004
Banf, M., Rhee, S. Y. (2017). Computational Inference of Gene Regulatory Networks: Approaches, Limitations and Opportunities. Biochim. Biophys. Acta BBA Gene Regul. Mech. 1860, 41–52. doi: 10.1016/j.bbagrm.2016.09.003
Bansal, M., Gatta, G. D., di Bernardo, D. (2006). Inference of Gene Regulatory Networks and Compound Mode of Action From Time Course Gene Expression Profiles. Bioinformatics 22, 815–822. doi: 10.1093/bioinformatics/btl003
Bartlett, H. S., Million, R. P. (2015). Targeting the IL-17–T H 17 Pathway. Nat. Rev. Drug Discovery 14, 11–12. doi: 10.1038/nrd4518
Bichiou, H., Bouabid, C., Rabhi, I., Guizani-Tabbane, L. (2021). Transcription Factors Interplay Orchestrates the Immune-Metabolic Response of Leishmania Infected Macrophages. Front. Cell. Infect. Microbiol. 11. doi: 10.3389/fcimb.2021.660415
Boaventura, V. S., Santos, C. S., Cardoso, C. R., de Andrade, J., Dos Santos, W. L. C., Clarêncio, J., et al. (2010). Human Mucosal Leishmaniasis: Neutrophils Infiltrate Areas of Tissue Damage That Express High Levels of Th17-Related Cytokines. Eur. J. Immunol. 40, 2830–2836. doi: 10.1002/eji.200940115
Bray, N. L., Pimentel, H., Melsted, P., Pachter, L. (2016). Near-Optimal Probabilistic RNA-Seq Quantification. Nat. Biotechnol. 34, 525–527. doi: 10.1038/nbt.3519
Burza, S., Croft, S. L., Boelaert, M. (2018). Leishmaniasis. Lancet Lond. Engl. 392, 951–970. doi: 10.1016/S0140-6736(18)31204-2
Castiglione, F., Tieri, P., Palma, A., Jarrah, A. S. (2016). Statistical Ensemble of Gene Regulatory Networks of Macrophage Differentiation. BMC Bioinf. 17, 506. doi: 10.1186/s12859-016-1363-4
Castilho, T. M., Goldsmith-Pestana, K., Lozano, C., Valderrama, L., Saravia, N. G., McMahon-Pratt, D. (2010). Murine Model of Chronic L. (Viannia) Panamensis Infection: Role of IL-13 in Disease. Eur. J. Immunol. 40, 2816–2829. doi: 10.1002/eji.201040384
Chandrasekaran, S., Ament, S. A., Eddy, J. A., Rodriguez-Zas, S. L., Schatz, B. R., Price, N. D., et al. (2011). Behavior-Specific Changes in Transcriptional Modules Lead to Distinct and Predictable Neurogenomic States. Proc. Natl. Acad. Sci. U. S. A. 108, 18020–18025. doi: 10.1073/pnas.1114093108
Chen, J., Bardes, E. E., Aronow, B. J., Jegga, A. G. (2009). ToppGene Suite for Gene List Enrichment Analysis and Candidate Gene Prioritization. Nucleic Acids Res. 37, W305–W311. doi: 10.1093/nar/gkp427
Contreras, I., Gómez, M. A., Nguyen, O., Shio, M. T., McMaster, R. W., Olivier, M. (2010). Leishmania-Induced Inactivation of the Macrophage Transcription Factor AP-1 Is Mediated by the Parasite Metalloprotease Gp63. PloS Pathog. 6, e1001148. doi: 10.1371/journal.ppat.1001148
Cummings, H. E., Tuladhar, R., Satoskar, A. R. (2010). Cytokines and Their STATs in Cutaneous and Visceral Leishmaniasis. J. Biomed. Biotechnol. 2010, 294389. doi: 10.1155/2010/294389
Dayakar, A., Chandrasekaran, S., Kuchipudi, S. V., Kalangi, S. K. (2019). Cytokines: Key Determinants of Resistance or Disease Progression in Visceral Leishmaniasis: Opportunities for Novel Diagnostics and Immunotherapy. Front. Immunol. 10. doi: 10.3389/fimmu.2019.00670
de Azevedo, A. F., Aurrecoechea, C., Berriman, M., Brestelli, J., Brunk, B. P., Carrington, M., et al. (2010). TriTrypDB: A Functional Genomic Resource for the Trypanosomatidae. Nucleic Acids Res. 38, D457–D462. doi: 10.1093/nar/gkp851
Delgado, F. M., Gómez-Vela, F. (2019). Computational Methods for Gene Regulatory Networks Reconstruction and Analysis: A Review. Artif. Intell. Med. 95, 133–145. doi: 10.1016/j.artmed.2018.10.006
de Luis Balaguer, M. A., Fisher, A. P., Clark, N. M., Fernandez-Espinosa, M. G., Möller, B. K., Weijers, D., et al. (2017). Predicting Gene Regulatory Networks by Combining Spatial and Temporal Gene Expression Data in Arabidopsis Root Stem Cells. Proc. Natl. Acad. Sci. U. S. A. 114, E7632–E7640. doi: 10.1073/pnas.1707566114
Dennis, G., Sherman, B. T., Hosack, D. A., Yang, J., Gao, W., Lane, H. C., et al. (2003). DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 4, R60. doi: 10.1186/gb-2003-4-9-r60
de Paiva, R. M. C., Grazielle-Silva, V., Cardoso, M. S., Nakagaki, B. N., Mendonça-Neto, R. P., Canavaci, A. M. C., et al. (2015). Amastin Knockdown in Leishmania Braziliensis Affects Parasite-Macrophage Interaction and Results in Impaired Viability of Intracellular Amastigotes. PloS Pathog. 11, e1005296. doi: 10.1371/journal.ppat.1005296
Diefenbach, A., Schindler, H., Donhauser, N., Lorenz, E., Laskay, T., MacMicking, J., et al. (1998). Type 1 Interferon (IFNalpha/beta) and Type 2 Nitric Oxide Synthase Regulate the Innate Immune Response to a Protozoan Parasite. Immunity 8, 77–87. doi: 10.1016/s1074-7613(00)80460-4
Ewels, P., Magnusson, M., Lundin, S., Käller, M. (2016). MultiQC: Summarize Analysis Results for Multiple Tools and Samples in a Single Report. Bioinformatics 32, 3047–3048. doi: 10.1093/bioinformatics/btw354
Fernandes, M. C., Dillon, L. A. L., Belew, A. T., Bravo, H. C., Mosser, D. M., El-Sayed, N. M. (2016). Dual Transcriptome Profiling of Leishmania-Infected Human Macrophages Reveals Distinct Reprogramming Signatures. mBio 7, 1–16. doi: 10.1128/mBio.00027-16
Gaffen, S. L., Jain, R., Garg, A. V., Cua, D. J. (2014). IL-23-IL-17 Immune Axis: Discovery, Mechanistic Understanding, and Clinical Testing. Nat. Rev. Immunol. 14, 585–600. doi: 10.1038/nri3707
Ghedira, K., Hornischer, K., Konovalova, T., Jenhani, A.-Z., Benkahla, A., Kel, A. (2011). Identification of Key Mechanisms Controlling Gene Expression in Leishmania Infected Macrophages Using Genome-Wide Promoter Analysis. Infect. Genet. Evol. 11, 769–777. doi: 10.1016/j.meegid.2010.10.015
Gonçalves-de-Albuquerque, S., da, C., Pessoa-e-Silva, R., Trajano-Silva, L. A. M., de Goes, T. C., de Morais, R. C. S., et al. (2017). The Equivocal Role of Th17 Cells and Neutrophils on Immunopathogenesis of Leishmaniasis. Front. Immunol. 8. doi: 10.3389/fimmu.2017.01437
Gonzalez, K., Calzada, J. E., Corbett, C. E. P., Saldaña, A., Laurenti, M. D. (2020). Involvement of the Inflammasome and Th17 Cells in Skin Lesions of Human Cutaneous Leishmaniasis Caused by Leishmania (Viannia) Panamensis. Mediators Inflamm. 2020, 9278931. doi: 10.1155/2020/9278931
Gonzalez-Lombana, C., Gimblet, C., Bacellar, O., Oliveira, W. W., Passos, S., Carvalho, L. P., et al. (2013). IL-17 Mediates Immunopathology in the Absence of IL-10 Following Leishmania Major Infection. PloS Pathog. 9, e1003243. doi: 10.1371/journal.ppat.1003243
Grünebast, J., Clos, J. (2020). Leishmania: Responding to Environmental Signals and Challenges Without Regulated Transcription. Comput. Struct. Biotechnol. J. 18, 4016–4023. doi: 10.1016/j.csbj.2020.11.058
Han, P., Gopalakrishnan, C., Yu, H., Wang, E. (2017). Gene Regulatory Network Rewiring in the Immune Cells Associated With Cancer. Genes 8, 308. doi: 10.3390/genes8110308
Hardy, S., Robillard, P. N. (2004). Modeling and Simulation of Molecular Biology Systems Using Petri Nets: Modeling Goals of Various Approaches. J. Bioinform. Comput. Biol. 02, 619–637. doi: 10.1142/S0219720004000764
Huang, D. W., Sherman, B. T., Lempicki, R. A. (2009). Bioinformatics Enrichment Tools: Paths Toward the Comprehensive Functional Analysis of Large Gene Lists. Nucleic Acids Res. 37, 1–13. doi: 10.1093/nar/gkn923
Hu, Y., Ota, N., Peng, I., Refino, C. J., Danilenko, D. M., Caplazi, P., et al. (2010). IL-17RC is Required for IL-17A- and IL-17F-Dependent Signaling and the Pathogenesis of Experimental Autoimmune Encephalomyelitis. J. Immunol. Baltim. Md. 1950 184, 4307–4316. doi: 10.4049/jimmunol.0903614
Isnard, A., Shio, M. T., Olivier, M. (2012). Impact of Leishmania Metalloprotease GP63 on Macrophage Signaling. Front. Cell. Infect. Microbiol. 2. doi: 10.3389/fcimb.2012.00072
Jackson, C. A., Castro, D. M., Saldi, G.-A., Bonneau, R., Gresham, D. (2020). Gene Regulatory Network Reconstruction Using Single-Cell RNA Sequencing of Barcoded Genotypes in Diverse Environments. eLife 9, e51254. doi: 10.7554/eLife.51254
Katara, G. K., Ansari, N. A., Singh, A., Ramesh, V., Salotra, P. (2012). Evidence for Involvement of Th17 Type Responses in Post Kala Azar Dermal Leishmaniasis (PKDL). PloS Negl. Trop. Dis. 6, e1703. doi: 10.1371/journal.pntd.0001703
Katara, G. K., Raj, A., Kumar, R., Avishek, K., Kaushal, H., Ansari, N. A., et al. (2013). Analysis of Localized Immune Responses Reveals Presence of Th17 and Treg Cells in Cutaneous Leishmaniasis Due to Leishmania Tropica. BMC Immunol. 14, 52. doi: 10.1186/1471-2172-14-52
Kuleshov, M. V., Jones, M. R., Rouillard, A. D., Fernandez, N. F., Duan, Q., Wang, Z., et al. (2016). Enrichr: A Comprehensive Gene Set Enrichment Analysis Web Server 2016 Update. Nucleic Acids Res. 44, W90–W97. doi: 10.1093/nar/gkw377
Law, C. W., Alhamdoosh, M., Su, S., Dong, X., Tian, L., Smyth, G. K., et al. (2018). RNA-Seq Analysis is Easy as 1-2-3 With Limma, Glimma and edgeR. F1000Research 5, 1–29. doi: 10.12688/f1000research.9005.3. ISCB Comm J-1408.
Lecoeur, H., Rosazza, T., Kokou, K., Varet, H., Coppée, J.-Y., Lari, A., et al. (2020). Leishmania Amazonensis Subverts the Transcription Factor Landscape in Dendritic Cells to Avoid Inflammasome Activation and Stall Maturation. Front. Immunol. 11. doi: 10.3389/fimmu.2020.01098
Lee, H. J., Georgiadou, A., Otto, T. D., Levin, M., Coin, L. J., Conway, D. J., et al. (2018). Transcriptomic Studies of Malaria: A Paradigm for Investigation of Systemic Host-Pathogen Interactions. Microbiol. Mol. Biol. Rev. MMBR 82, e00071–e00017. doi: 10.1128/MMBR.00071-17
Li, L. (2003). OrthoMCL: Identification of Ortholog Groups for Eukaryotic Genomes. Genome Res. 13, 2178–2189. doi: 10.1101/gr.1224503
Liu, D., Uzonna, J. E. (2012). The Early Interaction of Leishmania With Macrophages and Dendritic Cells and its Influence on the Host Immune Response. Front. Cell. Infect. Microbiol. 2. doi: 10.3389/fcimb.2012.00083
Livigni, A., O’Hara, L., Polak, M. E., Angus, T., Smith, L. B., Freeman, T. C. (2016). Petri Net-Based Graphical and Computational Modelling of Biological Systems. bioRxiv 13, 047043. doi: 10.1101/047043
Livigni, A., O’Hara, L., Polak, M. E., Angus, T., Wright, D. W., Smith, L. B., et al. (2018). A Graphical and Computational Modeling Platform for Biological Pathways. Nat. Protoc. 13 (4), 705–722.
Mandlik, V., Patil, S., Bopanna, R., Basu, S., Singh, S. (2016). Biological Activity of Coumarin Derivatives as Anti-Leishmanial Agents. PloS One 11, e0164585. doi: 10.1371/journal.pone.0164585
Martinez, P. A., Petersen, C. A. (2014). Chronic Infection by Leishmania Amazonensis Mediated Through MAPK ERK Mechanisms. Immunol. Res. 59, 153–165. doi: 10.1007/s12026-014-8535-y
Menezes, J. P. B., Almeida, T. F., Petersen, A. L. O. A., Guedes, C. E. S., Mota, M. S. V., Lima, J. G. B., et al. (2013). Proteomic Analysis Reveals Differentially Expressed Proteins in Macrophages Infected With Leishmania Amazonensis or Leishmania Major. Microbes Infect. 15, 579–591. doi: 10.1016/j.micinf.2013.04.005
Miyagi, P. E., Villani, E., Gustin, G. D. B., Maruyama, N., Santos Filho, D. J. (2002). Petri Net Approach for Modelling System Integration in Intelligent Buildings. J. Braz. Soc Mech. Sci. 24, 341–350. doi: 10.1590/S0100-73862002000400015
Murray, P. J. (2017). Macrophage Polarization. Annu. Rev. Physiol. 79, 541–566. doi: 10.1146/annurev-physiol-022516-034339
Mutiso, J. M., Macharia, J. C., Kiio, M. N., Ichagichu, J. M., Rikoi, H., Gicheru, M. M. (2013). Development of Leishmania Vaccines: Predicting the Future From Past and Present Experience. J. Biomed. Res. 27, 85–102. doi: 10.7555/JBR.27.20120064
O’Hara, L., Livigni, A., Theo, T., Boyer, B., Angus, T., Wright, D., et al. (2016). Modelling the Structure and Dynamics of Biological Pathways. PloS Biol. 14, e1002530. doi: 10.1371/journal.pbio.1002530
Ovalle-Bracho, C., Franco-Muñoz, C., Londoño-Barbosa, D., Restrepo-Montoya, D., Clavijo-Ramírez, C. (2015). Changes in Macrophage Gene Expression Associated With Leishmania (Viannia) Braziliensis Infection. PloS One 10, e0128934. doi: 10.1371/journal.pone.0128934
Pennisi, M., Cavalieri, S., Motta, S., Pappalardo, F. (2016). A Methodological Approach for Using High-Level Petri Nets to Model the Immune System Response. BMC Bioinf. 17, 498. doi: 10.1186/s12859-016-1361-6
Peter, I. S., Davidson, E. H. (2017). Assessing Regulatory Information in Developmental Gene Regulatory Networks. Proc. Natl. Acad. Sci. 114, 5862–5869. doi: 10.1073/pnas.1610616114
Pimentel, H., Bray, N. L., Puente, S., Melsted, P., Pachter, L. (2017). Differential Analysis of RNA-Seq Incorporating Quantification Uncertainty. Nat. Methods 14, 687–690. doi: 10.1038/nmeth.4324
Pitta, M. G. R., Romano, A., Cabantous, S., Henri, S., Hammad, A., Kouriba, B., et al. (2009). IL-17 and IL-22 are Associated With Protection Against Human Kala Azar Caused by Leishmania Donovani. J. Clin. Invest. 119, 2379–2387. doi: 10.1172/JCI38813
Pol, J. G., Caudana, P., Paillet, J., Piaggio, E., Kroemer, G. (2019). Effects of Interleukin-2 in Immunostimulation and Immunosuppression. J. Exp. Med. 217, e20191247. doi: 10.1084/jem.20191247
Robinson, M. D., McCarthy, D. J., Smyth, G. K. (2010). Edger: A Bioconductor Package for Differential Expression Analysis of Digital Gene Expression Data. Bioinformatics 26, 139–140. doi: 10.1093/bioinformatics/btp616
Ross, S. H., Cantrell, D. A. (2018). Signaling and Function of Interleukin-2 in T Lymphocytes. Annu. Rev. Immunol. 36, 411–433. doi: 10.1146/annurev-immunol-042617-053352
Singh, G. P., Jha, M., Singh, M., Naina (2020). Modeling the Mechanism Pathways of First Line Drug in Tuberculosis Using Petri Nets. Int. J. Syst. Assur. Eng. Manage. 11, 313–324. doi: 10.1007/s13198-019-00940-4
Sinha, S., Jones, B. M., Traniello, I. M., Bukhari, S. A., Halfon, M. S., Hofmann, H. A., et al. (2020). Behavior-Related Gene Regulatory Networks: A New Level of Organization in the Brain. Proc. Natl. Acad. Sci. 117, 23270–23279. doi: 10.1073/pnas.1921625117
Steele, E., Tucker, A., ‘t Hoen, P. A. C., Schuemie, M. J. (2009). Literature-Based Priors for Gene Regulatory Networks. Bioinform. Oxf. Engl. 25, 1768–1774. doi: 10.1093/bioinformatics/btp277
Steen, H. C., Gamero, A. M. (2012). The Role of Signal Transducer and Activator of Transcription-2 in the Interferon Response. J. Interferon Cytokine Res. 32, 103–110. doi: 10.1089/jir.2011.0099
Theocharidis, A., van Dongen, S., Enright, A. J., Freeman, T. C. (2009). Network Visualization and Analysis of Gene Expression Data Using BioLayout Express3D. Nat. Protoc. 4, 1535–1550. doi: 10.1038/nprot.2009.177
Tomiotto-Pellissier, F., Bortoleti, B. T., da, S., Assolini, J. P., Gonçalves, M. D., Carloto, A. C. M., et al. (2018). Macrophage Polarization in Leishmaniasis: Broadening Horizons. Front. Immunol. 9. doi: 10.3389/fimmu.2018.02529
Vargas-Inchaustegui, D. A., Xin, L., Soong, L. (2008). Leishmania Braziliensisinfection Induces Dendritic Cell Activation, ISG15 Transcription, and the Generation of Protective Immune Responses. J. Immunol. Baltim. Md. 1950 180, 7537–7545. doi: 10.4049/jimmunol.180.11.7537
Veras, P. S. T., Bezerra de Menezes, J. P. (2016). Using Proteomics to Understand How Leishmania Parasites Survive Inside the Host and Establish Infection. Int. J. Mol. Sci. 17, 1270. doi: 10.3390/ijms17081270
Veras, P. S. T., Ramos, P. I. P., de Menezes, J. P. B. (2018). In Search of Biomarkers for Pathogenesis and Control of Leishmaniasis by Global Analyses of Leishmania-Infected Macrophages. Front. Cell. Infect. Microbiol. 8. doi: 10.3389/fcimb.2018.00326
Wang, Z., Gerstein, M., Snyder, M. (2009). RNA-Seq: A Revolutionary Tool for Transcriptomics. Nat. Rev. Genet. 10, 57–63. doi: 10.1038/nrg2484
WHO, W H O (2020) Leishmaniasis (WHO). Available at: http://www.who.int/leishmaniasis/en/ (Accessed January 28, 2016).
Xin, L., Li, K., Soong, L. (2008). Down-Regulation of Dendritic Cell Signaling Pathways by Leishmania Amazonensis Amastigotes. Mol. Immunol. 45, 3371–3382. doi: 10.1016/j.molimm.2008.04.018
Keywords: Leishmania–macrophage interaction, gene regulatory network (GRN), expression profiling, Th17, immune modulation
Citation: Gonçalves LO, Pulido AFV, Mathias FAS, Enes AES, Carvalho MGR, de Melo Resende D, Polak ME and Ruiz JC (2022) Expression Profile of Genes Related to the Th17 Pathway in Macrophages Infected by Leishmania major and Leishmania amazonensis: The Use of Gene Regulatory Networks in Modeling This Pathway. Front. Cell. Infect. Microbiol. 12:826523. doi: 10.3389/fcimb.2022.826523
Received: 30 November 2021; Accepted: 09 March 2022;
Published: 14 June 2022.
Edited by:
Claudia Ida Brodskyn, Gonçalo Moniz Institute (IGM), BrazilReviewed by:
Saikat Majumder, University of Pittsburgh Institution, United StatesTaís Fontoura de Almeida, Federal University of Rio de Janeiro, Brazil
Copyright © 2022 Gonçalves, Pulido, Mathias, Enes, Carvalho, de Melo Resende, Polak and Ruiz. 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: Jeronimo C. Ruiz, jeronimo.ruiz@fiocruz.br; Marta E. Polak, m.ewa.polak@soton.ac.uk
†These authors have contributed equally to this work and share senior authorship