Skip to main content

Identification of differentially expressed genes and signaling pathways with Candida infection by bioinformatics analysis

Abstract

Background

Opportunistic Candida species causes severe infections when the human immune system is weakened, leading to high mortality.

Methods

In our study, bioinformatics analysis was used to study the high-throughput sequencing data of samples infected with four kinds of Candida species. And the hub genes were obtained by statistical analysis.

Results

A total of 547, 422, 415 and 405 differentially expressed genes (DEGs) of Candida albicans, Candida glabrata, Candida parapsilosis and Candida tropicalis groups were obtained, respectively. A total of 216 DEGs were obtained after taking intersections of DEGs from the four groups. A protein–protein interaction (PPI) network was established using these 216 genes. The top 10 hub genes (FOSB, EGR1, JUNB, ATF3, EGR2, NR4A1, NR4A2, DUSP1, BTG2, and EGR3) were acquired through calculation by the cytoHubba plug-in in Cytoscape software. Validated by the sequencing data of peripheral blood, JUNB, ATF3 and EGR2 genes were  significant statistical significance.

Conclusions

In conclusion, our study demonstrated the potential pathogenic genes in Candida species and their underlying mechanisms by bioinformatic analysis methods. Further, after statistical validation, JUNB, ATF3 and EGR2 genes were attained, which may be used as potential biomarkers with Candida species infection.

Background

Candida species colonizes human mucosal surfaces as commensals, which can turn to pathogenic behavior under certain conditions, especially in immunocompromised patients [1]. Opportunistic Candida species includes Candida albicans (C. albicans), Candida krusei (C. krusei), Candida tropicalis (C. tropicalis), Candida glabrata (C. glabrata), and Candida parapsilosis (C. parapsilosis). Among them, Candida albicans is the major pathogen isolated [2, 3], accounting for about 50% of all candidiasis cases with a total mortality rate of 43% [4, 5]. However, in the last decade, the proportion of Candida has been changed, with a reduction of Candida albicans and an increase in Candida parapsilosis, Candida tropicalis, and Candida glabrata [2, 6]. Non-albicans Candida such as Candida glabrata, Candida parapsilosis, and Candida tropicalis have now been identified as frequent human pathogens [7]. The number of patients and the geographic location determined that the overall distribution of the species Candida albicans was more common in Australia, Japan, Korea, Hong Kong, and Malaysia, whereas Candida tropicalis is more common in Pakistan and India [2, 8].

Candidiasis is a common bloodstream infection in hospitals around the world, causing high morbidity and mortality [2, 9]. In recent years, the emergence of Candida resistant strains brought about the further risk of clinical infection [10]. The major virulence factors of these pathogens were the Candida peptide and the extracellular aspartic proteases of the Candida peptide family [11]. Despite the introduction of intensive care facilities and modern antifungal drugs, the results of progress in curing Candida infections over the past decades have been disappointing [10, 12]. At present, adjuvant immunotherapy can further reduce the morbidity and mortality caused by Candida infection [13]. Therefore, understanding how host defense pathways participate in candidiasis is crucial for determining new targets for immunotherapy.

Sequencing technology has been applied to find targets for immunotherapies. High-throughput sequencing has become an important method for studying genomics, epigenomics, and transcriptome [14]. At present, some studies applied high-throughput sequencing to the identification and characterization of clinical microbiology. And when clinical data are more complex and contains multiple species, this technology is more reliable than normal sequencing [15].

To study the differentially expressed genes (DEGs) and signaling pathways related to Candida infection, this study applied bioinformatics analysis to analyze the gene expression profiles of human whole blood infected by four common Candida species including Candida albicans, Candida parapsilosis, Candida glabrata, and Candida tropicalis.

Methods

Data sources

The Gene Expression Omnibus (GEO, http://www.ncbi.nlm.nih.gov/geo/) database, attached to the National Center for Biotechnology Information (NCBI), was used to store gene expression datasets, series, and platform records. The gene expression profiles of GSE114174, GSE114175, GSE114177, GSE114178, and GSE114179 provided by Philipp Kämmer from the GEO database were downloaded. These gene expression profiles were attained by high-throughput sequencing. The number of infection and control samples is shown in Table 1. As shown in Table 1, The number of whole blood sample infected with Candida albicans, Candida glabrata, Candida parapsilosis and Candida tropicalis are 15, 15, 15 and 15, respectively. As a control group, GSE114179 includes three whole blood sample. These whole blood samples were donated by German volunteers. Additionally, the sources of samples infected with different Candida species were from Germany, America and Netherlands. Further details of the data sources from the GEO database for this study are shown in Table 1.

Table 1 Details of the data sources from Gene Expression Omnibus (GEO) for this study

Data possessing and identification of DEGs

With R software, the data of GSE114174, GSE114175, GSE114177, GSE114178 and GSE114179 were batch-corrected and standardized using the affy software package, and the screening and identification of DEGs were carried out using the limma software package. Additionally, the ggplot2 package was applied to draw a volcano map, whereas the pheatmap package was employed to make a heatmap to visualize the DEGs. Using P < 0.05 and |log2FC|≥ 1 as critical values, the gene expression profiles in infected samples and uninfected samples were compared.

Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of DEGs

Gene Ontology was the main bioinformatics means used to annotate genes and analyze their biological processes. It includes the biological process (BP), molecular function (MF), and cellular component (CC), which reflects the conceptual category of gene product function, the biological processes of DNA metabolism and known differences between different organisms [16]. The KEGG is a database that can be used to understand biological systems and advanced functions from large-scale sequencing data generated by high-throughput sequencing technology [17]. In order to illustrate the biological function of genes and signal pathways involved in the vivo and cells, differentially expressed genes were annotated based on GO and KEGG analysis. The clusterProfler package in R software was used to performed the two analyses.

Construction of protein–protein interaction (PPI) network and identification of hub genes

The DEGs of four transcription profile data (GSE114174, GSE114175, GSE114177 and GSE114178) were overlapped by using online tools (https://bioinfogp.cnb.csic.es/tools/venny/index.html). The construction of the PPI network was carried out using the Search Tool for the Retrieval of Interacting Genes (STRING) database (http://string-db.org/), which aimed at providing important estimates and integrations of protein–protein interactions, including functional and physical association [18]. Moreover, the cytoHubba plug-in in Cytoscape 3.7.2 software was applied to obtain the top 10 hub genes which ranked by Maximal Clique Centrality (MCC).

Verification of intersection hub genes and construction of intersection gene–miRNA interaction

After identifying the intersection hub genes from four groups of data, the data of GSE42630 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE42630) obtained after the bioinformatics analysis were used for verification. GraphPad Prism 8.0 software was used for statistical analysis. The Gene–miRNA interaction was constructed using the miRTarBase v8.0 database [19] by Network Analyst (https://www.networkanalyst.ca/NetworkAnalyst/home.xhtml).

Results

Identification of DEGs associated with Candida infection

A total of 547, 422, 415 and 405 DEGs of Candida albicans, Candida glabrata, Candida parapsilosis and Candida tropicalis groups were separately obtained. Table 2 shows the detailed results of DEGs between the infection and control groups. The heat map of the DEGs proved that these DEGs could clearly distinguish between control samples and samples infected by C. albicans (Fig. 1B), Candida glabrata (Fig. 1D), Candida parapsilosis (Fig. 1F) and Candida tropicalis (Fig. 1H), respectively. Figure 1A, C, E and G, respectively, shows the volcano plot of the DEGs for C. albicans, C. glabrata, C. parapsilosis and C. tropicalis.

Table 2 Detailed results of differentially expressed genes (DEGs) between infection and control group
Fig. 1
figure 1

Heat map and volcano map of DEGs from four groups of Candida species. AB Candida albicans, CD Candida glabrata, EF Candida parapsilosis, and GH Candida tropicalis. In heat maps, gene expression data are converted into a data matrix. Each column represents the genetic data of a sample, and each row represents a gene. The color of each cell represents the expression level, and there are references to expression levels in different colors in the upper right corner of the figure. In volcano maps, red dots indicated up-regulated genes. The green or blue dots indicated down-regulated genes. Black dots indicated the rest of the genes with no significant expression change. The threshold was set as followed: P < 0.05 and |log2FC|≥ 1. FC: fold change

Go enrichment analysis of DEGs associated with Candida infection

Figures 2, 3 and 4 represent the results of biological processes, cellular components and molecular function in the GO enrichment analysis. For the biological process, both the C. tropicalis and C. parapsilosis groups were significantly enriched in T cell activation, whereas the C. albicans, C. tropicalis and C. parapsilosis groups were significantly enriched in leukocyte differentiation and regulation of leukocyte activation. With regard to the cellular component, the C. albicans, C. glabrata, C. parapsilosis and C. tropicalis groups were significantly enriched in the cytosolic ribosome and cytosolic part. Regarding the molecular function, both C. tropicalis and C. parapsilosis groups were significantly enriched in cytokine activity and DNA-binding transcription activator activity, RNA polymerase II-specific. The top 5 terms of significant enrichment of GO analysis for Candida albicans, Candida glabrata, Candida parapsilosis and Candida tropicalis in the Additional files 1, 2, 3, 4 (Tables S1–S4).

Fig. 2
figure 2

The biological processes (BP) of GO analysis from four groups of Candida species. A Candida albicans, B Candida glabrata, C Candida parapsilosis, and D Candida tropicalis

Fig. 3
figure 3

The cellular components (CC) of GO analysis from four groups of Candida species. A Candida albicans, B Candida glabrata, C Candida parapsilosis, and D Candida tropicalis

Fig. 4
figure 4

The molecular function (MF) of GO analysis from four groups of Candida species. A Candida albicans, B Candida glabrata, C Candida parapsilosis, and D Candida tropicalis

KEGG pathway enrichment analysis of DEGs associated with Candida infection

KEGG pathway analysis was conducted to identify the biological functions of DEGs. The results of the analysis are shown in Fig. 5. Five significant enrichment pathways were found simultaneously in C. albicans, C. glabrata, C. parapsilosis and C. tropicalis, which included the NF-kappa B signaling pathway, TNF signaling pathway, viral protein interaction with cytokine and cytokine receptor, salmonella infection and osteoclast differentiation. Additionally, the Toll-like receptor signaling pathway and MAPK signaling pathway were found simultaneously in C. glabrata, C. parapsilosis and C. tropicalis. The top 10 terms of significantly enriched KEGG pathways for Candida albicans, Candida glabrata, Candida parapsilosis and Candida tropicalis in the Additional files 5, 6, 7, 8 (Tables S5–S8).

Fig. 5
figure 5

The analysis of KEGG pathways from four groups of Candida species.  A Candida albicans, B Candida glabrata, C Candida parapsilosis, and D Candida tropicalis

Construction of PPI network and hub genes identification

The overlapping results of DEGs of four groups of Candida species are shown in Fig. 6A. The PPI network and the top 10 hub genes of the C. albicans, C. glabrata, C. parapsilosis, and C. tropicalis groups are, respectively, shown in Additional files 9, 10 (Figs. S1 and S2). A total of 216 DEGs were obtained after taking the intersection of the DEGs from the four groups. The 216 genes were uploaded to STRING to construct a PPI network, consisting of 51 nodes and 96 edges (Fig. 6C). The top 10 hub genes obtained included FOSB, EGR1, JUNB, ATF3, EGR2, NR4A1, NR4A2, DUSP1, BTG2, and EGR3 (Fig. 6D).

Fig. 6
figure 6

The intersection results and relevant analysis of DEGs from four groups of Candida species. A The overlapping results of DEGs of four groups of Candida species. Crossed regions indicate co-expressed DEGs. B The gene-miRNA network of overlapping DEGs based on miRTarBase v8.0 database. C The PPI network of overlapping DEGs. D The top 10 hub genes of PPI network

Construction of intersection gene and miRNA interaction and verification of intersection hub genes

The Gene-miRNA interaction of intersection genes is shown in Figure 6B. According to Degree, the top 10 genes interacted with miRNA were: BTG2, PMAIP1, GLUL, HS3ST1, DDX3X, TFRC, ORC6, SGK1, ZFP36 and RPL23A. In consideration of the rigorousness of this study, the data from the GSE42630 were used to verify the 10 intersection hub genes obtained. With P < 0.05 as the standard, the analytic results of the genes (JUNB, ATF3, EGR2) were found to be statistically significant (Fig. 7).

Fig. 7
figure 7

Verification of hub genes associated with Candida infection. P-value < 0.05 is considered to be statistically significant

Discussion

Candida species is one of the most common pathogens of invasive fungal infections among hospitalized patients [20]. Bioinformatics analysis can quantitatively analyze the gene expression of Candida species and identify the differentially expressed genes generated in cells, tissues or organisms that were infected with Candida species and uninfected. It is extremely important to understand the molecular mechanism of genetic regulation of Candida species and to better treat and diagnose diseases [21]. We downloaded the gene expression profiles relevant to the Candida infection from the GEO database. The obtained DEGs were used to perform GO enrichment analysis and KEGG pathway analysis, construct the PPI network, and gene–miRNA interaction network and identify the top 10 hub genes.

The GO enrichment analysis indicated that an increase in pattern recognition receptor activity may enhance the host’s defense against Candida species. The study has shown that when infected with Candida species, the first step to develop an immune response to Candida species is the recognition of invasive fungi [22]. At present, studies have found that PRRs of Candida albicans include TLRs, CLRs, NLRs and RIG-I-like receptors (RLRs) [22, 23]. Additionally, in our study, the leukocyte differentiation, regulation of leukocyte activation and T cells activation increased resistance to the invasion of Candida species, which was consistent with previous studies [24,25,26]. In addition to killing Candida by the production of reactive oxygen species (ROS) and phagocytosis, activated neutrophils can also release neutrophil extracellular traps (NETs), that capture Candida conidia and hyphae and contain the antimicrobial proteins to inhibit fungal growth [27, 28]. TH17 cell responses played an important role in mucosal host defenses against Candida by producing IL-17 and IL-22. These cytokines recruit and activate neutrophils, activate epithelial cells and release antifungal β-defensins cooperatively [29]. TH1 cell responses and IFNγ productions were quite important for the fungicidal activities of both neutrophils and macrophages [30].

In the current study, certain cellular signaling pathways identified using KEGG analysis were closely associated with Candida fungal infections. In addition to recognizing a wide variety of microbial products including lipoproteins, flagellin, lipopolysaccharides and bacterial DNA, signal transduction through TLRs also led to the production of inflammatory mediators [31]. Previous studies have shown that TLR4 mediated the induction of pro-inflammatory cytokines after Candida stimulation, whereas TLR2’s recognition of Candida primarily led to the release of anti-inflammatory cytokines [32]. TLRs initiate downstream signaling that culminates in the activation of nuclear factor kappa B (NF-κB), mitogen activating protein (MAP) kinases, or Interferon regulatory factors (IRFs), to regulate the expression of type I IFNs, cytokines and chemokines that ultimately should protect the host from infection by pathogens [33]. Primarily produced by activated macrophages, TNF signals are transmitted through two different cell surface receptors, TNF-R1 and TNF-R2 [34]. A number of experimental studies have revealed that the TNF-R1 activates most of the biological activity of TNF. The binding of TNF to TNF-R1 initiates downstream signaling that culminates in the activation of NF-kB and c-Jun, two major transcription factors [35, 36]. The NF-kappa B pathway is divided into two different but interacting pathways: the classical NF-kappa B essential modulator (NEMO)—a dependent pathway and the alternate NEMO—an independent pathway. While the classical NF-κB signaling pathway, induced by TNF-α, IL-1, or by-products of bacterial and viral infections, is mainly associated with inflammatory, proliferative, and survival responses, the activation of the noncanonical pathway results in a chemokine expression. Taken together, detecting these pathways may be helpful to predict the progression of Candidiasis [37].

Protein–protein interactions and acquired networks are very important in most biological functions and processes, as most proteins seem to activate their functions through interactions [38]. The hub genes screened through the PPI network are closely related to the potential molecular mechanism of Candida infection in humans. Therefore, a total of 10 hub genes were selected in this study. Keeping in mind the rigor of this study, the data of GSE42630 was used to verify the 10 hub genes. Through verification, among the 10 hub genes, the three hub genes—JUNB, ATF3, and EGR2—were significant statistical significance. Studies have shown that JUNB possesses an important effect during the growth of Treg cells, as it promotes IL-2 signal transduction [39]. Therefore, during Candida infection, the regulation of JUNB may affect Treg cells in resisting Candida infection. The research of Rynes et al. found that ATF3 can maintain the homeostasis of the metabolism and immune system [40, 41]. The loss of ATF3 can cause chronic inflammation. Through NF-B/Relish in the ATF3 mutant, the overactive pro-inflammatory and stress signals caused by Jun N-terminal kinase and FOXO can remove the regulation of the important genes in immune defense [42]. In Candida infection, the main role of ATF3 is to inhibit inflammation. Studies have found that the adaptive immune response is regulated by EGR2 and EGR3 by uncoupling and expanding the time of T cell differentiation. EGR2 binds to and controls the expression of proliferation regulating genes (Myc and Myb), differentiation inhibitors (Bcl6, Id3) and inhibits transcription factors (Zeb2, RORa, RORc, and Bhlhe40) required for effector functions. EGR2 and EGR3 are upstream regulation factors of CD4 and CD8 T cells, which are essential for optimal response under limited immunopathology [43, 44]. In our research, EGR2 was a key adjustment factor. While its impact has been confirmed, the effect of EGR3 was not reflected clearly.

BTG2 is an archetype member of the BTG/Tob antiproliferative protein family, and its expression is related to various cellular processes, for instance, the generation cycle, divergence, or apoptosis of cells. BTG2 may act as a regulatory factor of the intracellular signal transduction cascade [45, 46]. BTG2 expression is induced through a p53-dependent mechanism, and the function of BTG2 may be related to cell cycle control and DNA damage reaction [47, 48].

Overall, bioinformatics analysis can be used to study the complex underlying molecular mechanisms related to diseases. In this study, the hub genes related to Candida infection were identified. However, further experimentation is required to verify these predicted results from bioinformatics analysis. This study has some limitations. Firstly, the number of sample of each group of Candida species was 18, therefore, the sample size of the study was relatively small. Secondly, the research sample did not eliminate factors like gender, whether the gene expression profiles were infected with other diseases, or whether drugs were used. This may have affected the factors of Candida infection in the gene expression.

Availability of data and materials

The dataset(s) supporting the conclusions of this article is available in the Gene Expression Omnibus (GEO) database, http://www.ncbi.nlm.nih.gov/geo/.

Abbreviations

GEO:

Gene expression omnibus

DEGs:

Differentially expressed genes

GO:

Gene Ontology

KEGG:

Kyoto Encyclopedia of Genes and Genomes

PPI:

Protein–protein interaction

C. albicans :

Candida albicans

C. glabrata :

Candida glabrata

C. tropicalis :

Candida tropicalis

C. parapsilosis :

Candida parapsilosis

C. krusei :

Candida krusei

NCBI:

National Center for Biotechnology Information

MF:

Molecular function

BP:

Biological process

CC:

Cellular component

MCC:

Maximal Clique Centrality

PRRs:

Pattern recognition receptors

PAMPs:

Pathogen-associated molecular patterns

TLRs:

Toll-like receptors

CLRs:

C-type lectin receptors

NLRs:

NOD-like receptors

RLRs:

RIG-I-like receptors

ROS:

Reactive oxygen species

DCs:

Dendritic cells

TH:

T helper

IL-17:

Interleukin-17

IFNγ:

Interferon-γ

TNF:

Tumor necrosis factor

NEMO:

NF-kappa B essential modulator

NIK:

NF-κB-inducing kinase

References

  1. Sardi JCO, Scorzoni L, Bernardi T, Fusco-Almeida AM, Mendes Giannini MJS. Candida species: current epidemiology, pathogenicity, biofilm formation, natural antifungal products and new therapeutic options. J Med Microbiol. 2013;62(Pt 1):10–24.

    CAS  PubMed  Google Scholar 

  2. Lamoth F, Lockhart SR, Berkow EL, Calandra T. Changes in the epidemiological landscape of invasive candidiasis. J Antimicrob Chemother. 2018;73(suppl_1):i4-i13.

  3. Xiao Z, Wang Q, Zhu F, An Y. Epidemiology, species distribution, antifungal susceptibility and mortality risk factors of candidemia among critically ill patients: a retrospective study from 2011 to 2017 in a teaching hospital in China. Antimicrob Resist Infect Control. 2019;29(8):89.

    Google Scholar 

  4. Drgona L, Khachatryan A, Stephens J, Charbonneau C, Kantecki M, Haider S, et al. Clinical and economic burden of invasive fungal diseases in Europe: focus on pre-emptive and empirical treatment of Aspergillus and Candida species. Eur J Clin Microbiol Infect Dis. 2014;33(1):7–21.

    CAS  PubMed  Google Scholar 

  5. Pfaller MA, Castanheira M. Nosocomial candidiasis: antifungal stewardship and the importance of rapid diagnosis. Med Mycol. 2016;54(1):1–22.

    PubMed  Google Scholar 

  6. Enoch DA, Yang H, Aliyu SH, Micallef C. The changing epidemiology of invasive fungal infections. Methods Mol Biol. 2017;1508:17–65.

    CAS  PubMed  Google Scholar 

  7. Whaley SG, Berkow EL, Rybak JM, Nishimoto AT, Barker KS, Rogers PD. Azole antifungal resistance in Candida albicans and emerging non-albicans Candida species. Front Microbiol. 2016;7:2173.

    PubMed  Google Scholar 

  8. Wang H, Xu YC, Hsueh PR. Epidemiology of candidemia and antifungal susceptibility in invasive Candida species in the Asia-Pacific region. Future Microbiol. 2016;11:1461–77.

    CAS  PubMed  Google Scholar 

  9. Israel S, Amit S, Israel A, Livneh A, Nir-Paz R, Korem M. The epidemiology and susceptibility of Candidemia in Jerusalem, Israel. Front Cell Infect Microbiol. 2019;11(9):352.

    Google Scholar 

  10. Pristov KE, Ghannoum MA. Resistance of Candida to azoles and echinocandins worldwide. Clin Microbiol Infect. 2019;25(7):792–8.

    CAS  PubMed  Google Scholar 

  11. Rapala-Kozik M, Bochenska O, Zajac D, Karkowska-Kuleta J, Gogol M, Zawrotniak M, et al. Extracellular proteinases of Candida species pathogenic yeasts. Mol Oral Microbiol. 2018;33(2):113–24.

    CAS  PubMed  Google Scholar 

  12. Smeekens SP, Ng A, Kumar V, Johnson MD, Plantinga TS, van Diemen C, et al. Functional genomics identifies type I interferon pathway as central for host defense against Candida albicans. Nat Commun. 2013;4:1342.

    PubMed  Google Scholar 

  13. Nami S, Aghebati-Maleki A, Morovati H, Aghebati-Maleki L. Current antifungal drugs and immunotherapeutic approaches as promising strategies to treatment of fungal diseases. Biomed Pharmacother. 2019;110:857–68.

    CAS  PubMed  Google Scholar 

  14. Churko JM, Mantalas GL, Snyder MP, Wu JC. Overview of high throughput sequencing technologies to elucidate molecular pathways in cardiovascular diseases. Circ Res. 2013;112(12):1613–23.

    CAS  PubMed  Google Scholar 

  15. Deurenberg RH, Bathoorn E, Chlebowicz MA, Couto N, Ferdous M, García-Cobos S, et al. Application of next generation sequencing in clinical microbiology and infection prevention. J Biotechnol. 2017;243:16–24.

    CAS  PubMed  Google Scholar 

  16. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25(1):25–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  17. Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2012;40(Database issue):D109–14.

  18. Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, et al. STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015;43(Database issue):D447–52.

  19. Chou CH, Shrestha S, Yang CD, Chang NW, Lin YL, Liao KW, et al. miRTarBase update 2018: a resource for experimentally validated microRNA-target interactions. Nucleic Acids Res. 2018;46(D1):D296–302.

    CAS  PubMed  Google Scholar 

  20. Antinori S, Milazzo L, Sollima S, Galli M, Corbellino M. Candidemia and invasive candidiasis in adults: a narrative review. Eur J Intern Med. 2016;34:21–8.

    PubMed  Google Scholar 

  21. Ji F, Sadreyev RI. RNA-seq: basic bioinformatics analysis. Curr Protocols Mol Biol. 2018;124(1):e68.

    Google Scholar 

  22. Gow NA, van de Veerdonk FL, Brown AJ, Netea MG. Candida albicans morphogenesis and host defence: discriminating invasion from colonization. Nat Rev Microbiol. 2011;10(2):112–22.

    PubMed  PubMed Central  Google Scholar 

  23. Netea MG, Joosten LA, van der Meer JW, Kullberg BJ, van de Veerdonk FL. Immune defence against Candida fungal infections. Nat Rev Immunol. 2015;15(10):630–42.

    CAS  PubMed  Google Scholar 

  24. Pathakumari B, Liang G, Liu W. Immune defence to invasive fungal infections: a comprehensive review. Biomed Pharmacother. 2020;130:110550.

    CAS  PubMed  Google Scholar 

  25. Castellano-Gonzalez G, Clancy LE, Gottlieb D. Prospects for adoptive T-cell therapy for invasive fungal disease. Curr Opin Infect Dis. 2017;30(6):518–27.

    CAS  PubMed  Google Scholar 

  26. Bär E, Gladiator A, Bastidas S, Roschitzki B, Acha-Orbea H, Oxenius A, LeibundGut-Landmann S. A novel Th cell epitope of Candida albicans mediates protection from fungal infection. J Immunol. 2012;188(11):5636–43.

    PubMed  Google Scholar 

  27. McCormick A, Heesemann L, Wagener J, Marcos V, Hartl D, Loeffler J, et al. NETs formed by human neutrophils inhibit growth of the pathogenic mold Aspergillus fumigatus. Microbes Infect. 2010;12(12–13):928–36.

    CAS  PubMed  Google Scholar 

  28. Burgener SS, Schroder K. Neutrophil extracellular traps in host defense. Cold Spring Harbor Perspect Biol. 2020;12(7):a037028.

    CAS  Google Scholar 

  29. Eyerich S, Wagener J, Wenzel V, Scarponi C, Pennino D, Albanesi C, et al. IL-22 and TNF-α represent a key cytokine combination for epidermal integrity during infection with Candida albicans. Eur J Immunol. 2011;41(7):1894–901.

    CAS  PubMed  Google Scholar 

  30. Kandel-Kfir M, Garcia-Milan R, Gueta I, Lubitz I, Ben-Zvi I, Shaish A, et al. IFNγ potentiates TNFα/TNFR1 signaling to induce FAT10 expression in macrophages. Mol Immunol. 2020;117:101–9.

    CAS  PubMed  Google Scholar 

  31. Kawasaki T, Kawai T. Toll-like receptor signaling pathways. Front Immunol. 2014;5:461.

    PubMed  PubMed Central  Google Scholar 

  32. Gasparoto TH, Tessarolli V, Garlet TP, Torres SA, Garlet GP, da Silva JS, et al. Absence of functional TLR4 impairs response of macrophages after Candida albicans infection. Med Mycol. 2010;48(8):1009–17.

    CAS  PubMed  Google Scholar 

  33. Kawai T, Akira S. The role of pattern-recognition receptors in innate immunity: update on Toll-like receptors. Nat Immunol. 2010;11(5):373–84.

    CAS  PubMed  Google Scholar 

  34. Horiuchi T, Mitoma H, Harashima S, Tsukamoto H, Shimoda T. Transmembrane TNF-alpha: structure, function and interaction with anti-TNF agents. Rheumatology (Oxford). 2010;49(7):1215–28.

    CAS  Google Scholar 

  35. Chen G, Goeddel DV. TNF-R1 signaling: a beautiful pathway. Science (New York, NY). 2002;296(5573):1634–5.

    CAS  Google Scholar 

  36. Shah A, Plaza-Sirvent C, Weinert S, Buchbinder JH, Lavrik IN, Mertens PR, et al. YB-1 mediates TNF-induced pro-survival signaling by regulating NF-κB activation. Cancers. 2020;12(8):2188.

    CAS  PubMed Central  Google Scholar 

  37. Mitchell S, Vargas J, Hoffmann A. Signaling via the NFκB system. Wiley Interdiscipl Rev Syst Biol Med. 2016;8(3):227–41.

    CAS  Google Scholar 

  38. Athanasios A, Charalampos V, Vasileios T, Ashraf GM. Protein–protein interaction (PPI) network: recent advances in drug discovery. Curr Drug Metab. 2017;18(1):5–10.

    CAS  PubMed  Google Scholar 

  39. Katagiri T, Yamazaki S, Fukui Y, Aoki K, Yagita H, Nishina T, et al. JunB plays a crucial role in development of regulatory T cells by promoting IL-2 signaling. Mucosal Immunol. 2019;12(5):1104–17.

    CAS  PubMed  Google Scholar 

  40. Sha H, Zhang D, Zhang Y, Wen Y, Wang Y. ATF3 promotes migration and M1/M2 polarization of macrophages by activating tenascin-C via Wnt/β-catenin pathway. Mol Med Rep. 2017;16(3):3641–7.

    CAS  PubMed  Google Scholar 

  41. Bambouskova M, Gorvel L, Lampropoulou V, Sergushichev A, Loginicheva E, Johnson K, et al. Electrophilic properties of itaconate and derivatives regulate the IκBζ-ATF3 inflammatory axis. Nature. 2018;556(7702):501–4.

    CAS  PubMed  PubMed Central  Google Scholar 

  42. Rynes J, Donohoe CD, Frommolt P, Brodesser S, Jindra M, Uhlirova M. Activating transcription factor 3 regulates immune and metabolic homeostasis. Mol Cell Biol. 2012;32(19):3949–62.

    CAS  PubMed  PubMed Central  Google Scholar 

  43. Miao T, Symonds ALJ, Singh R, Symonds JD, Ogbe A, Omodho B, et al. Egr2 and 3 control adaptive immune responses by temporally uncoupling expansion from T cell differentiation. J Exp Med. 2017;214(6):1787–808.

    CAS  PubMed  PubMed Central  Google Scholar 

  44. Morita K, Okamura T, Inoue M, Komai T, Teruya S, Iwasaki Y, et al. Egr2 and Egr3 in regulatory T cells cooperatively control systemic autoimmunity through Ltbp3-mediated TGF-β3 production. Proc Natl Acad Sci U S A. 2016;113(50):E8131–40.

    CAS  PubMed  PubMed Central  Google Scholar 

  45. Mauxion F, Faux C, Séraphin B. The BTG2 protein is a general activator of mRNA deadenylation. EMBO J. 2008;27(7):1039–48.

    CAS  PubMed  PubMed Central  Google Scholar 

  46. Bai Y, Qiao L, Xie N, Shi Y, Liu N, Wang J. Expression and prognosis analyses of the Tob/BTG antiproliferative (APRO) protein family in human cancers. PLoS ONE. 2017;12(9):e0184902.

    PubMed  PubMed Central  Google Scholar 

  47. Yuniati L, Scheijen B, van der Meer LT, van Leeuwen FN. Tumor suppressors BTG1 and BTG2: beyond growth control. J Cell Physiol. 2019;234(5):5379–89.

    CAS  PubMed  Google Scholar 

  48. Hwang SS, Lim J, Yu Z, Kong P, Sefik E, Xu H, Harman CCD, Kim LK, Lee GR, Li HB, Flavell RA. mRNA destabilization by BTG1 and BTG2 maintains T cell quiescence. Science. 2020;367(6483):1255–60.

    CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We acknowledge GEO database for providing their platforms and contributors for uploading their meaningful datasets.

Funding

This work was supported by the National Key Research and Development Program of China (2021YFC2009404), Guangzhou Planned Project of Science and Technology (202102010030) and Guangzhou General Science and Technology Project of Health and Family Planning (20191A011113).

Author information

Authors and Affiliations

Authors

Contributions

XG conceived and designed the experiments. GZ and LX analyzed the data and made the figures and tables. GZ, LX, JS, XC, XY, YL, YG participated in the writing, reading, and revising of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Xu-Guang Guo.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Significant enrichment of GO terms for Candida albicans (top 5 according to P value).

Additional file 2: Table S2.

Significant enrichment of GO terms for Candida glabrata (top 5 according to P value).

Additional file 3: Table S3.

Significant enrichment of GO terms for Candida parapsilosis (top 5 according to P value).

Additional file 4: Table S4.

Significant enrichment of GO terms for Candida tropicalis (top 5 according to P value).

Additional file 5: Table S5.

Top 10 significantly enriched KEGG pathways of Candida albicans (according to P value).

Additional file 6: Table S6.

Top 10 significantly enriched KEGG pathways of Candida glabrata (according to P value).

Additional file 7: Table S7.

Top 10 significantly enriched KEGG pathways of Candida parapsilosis (according to P value).

Additional file 8: Table S8.

Top 10 significantly enriched KEGG pathways of Candida tropicalis (according to P value).

Additional file 9: Figure S1.

Construction of the PPI network from four groups of Candida. Candida albicans (A), Candida glabrata (B), Candida parapsilosis (C) and Candida tropicalis (D).

Additional file 10: Figure S2.

Top 10 hub genes of four groups of Candida. Candida albicans (A), Candida glabrata (B), Candida parapsilosis (C) and Candida tropicalis (D).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zhu, GD., Xie, LM., Su, JW. et al. Identification of differentially expressed genes and signaling pathways with Candida infection by bioinformatics analysis. Eur J Med Res 27, 43 (2022). https://doi.org/10.1186/s40001-022-00651-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40001-022-00651-w

Keywords

  • Candida
  • High-throughput sequencing
  • Differentially expressed genes
  • Signaling pathways
  • Bioinformatics analysis