- Open Access
Identification of differentially expressed genes and signaling pathways with Candida infection by bioinformatics analysis
European Journal of Medical Research volume 27, Article number: 43 (2022)
Opportunistic Candida species causes severe infections when the human immune system is weakened, leading to high mortality.
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.
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.
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.
Candida species colonizes human mucosal surfaces as commensals, which can turn to pathogenic behavior under certain conditions, especially in immunocompromised patients . 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 . 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 . The major virulence factors of these pathogens were the Candida peptide and the extracellular aspartic proteases of the Candida peptide family . 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 . 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 . 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 .
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.
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.
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 . 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 . 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 . 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  by Network Analyst (https://www.networkanalyst.ca/NetworkAnalyst/home.xhtml).
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.
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).
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).
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).
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).
Candida species is one of the most common pathogens of invasive fungal infections among hospitalized patients . 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 . 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 . 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 . TH1 cell responses and IFNγ productions were quite important for the fungicidal activities of both neutrophils and macrophages .
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 . 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 . 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 . Primarily produced by activated macrophages, TNF signals are transmitted through two different cell surface receptors, TNF-R1 and TNF-R2 . 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 .
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 . 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 . 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 . 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/.
Gene expression omnibus
Differentially expressed genes
Kyoto Encyclopedia of Genes and Genomes
- C. albicans :
- C. glabrata :
- C. tropicalis :
- C. parapsilosis :
- C. krusei :
National Center for Biotechnology Information
Maximal Clique Centrality
Pattern recognition receptors
Pathogen-associated molecular patterns
C-type lectin receptors
Reactive oxygen species
Tumor necrosis factor
NF-kappa B essential modulator
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.
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.
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.
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.
Pfaller MA, Castanheira M. Nosocomial candidiasis: antifungal stewardship and the importance of rapid diagnosis. Med Mycol. 2016;54(1):1–22.
Enoch DA, Yang H, Aliyu SH, Micallef C. The changing epidemiology of invasive fungal infections. Methods Mol Biol. 2017;1508:17–65.
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.
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.
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.
Pristov KE, Ghannoum MA. Resistance of Candida to azoles and echinocandins worldwide. Clin Microbiol Infect. 2019;25(7):792–8.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Ji F, Sadreyev RI. RNA-seq: basic bioinformatics analysis. Curr Protocols Mol Biol. 2018;124(1):e68.
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.
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.
Pathakumari B, Liang G, Liu W. Immune defence to invasive fungal infections: a comprehensive review. Biomed Pharmacother. 2020;130:110550.
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.
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.
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.
Burgener SS, Schroder K. Neutrophil extracellular traps in host defense. Cold Spring Harbor Perspect Biol. 2020;12(7):a037028.
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.
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.
Kawasaki T, Kawai T. Toll-like receptor signaling pathways. Front Immunol. 2014;5:461.
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.
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.
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.
Chen G, Goeddel DV. TNF-R1 signaling: a beautiful pathway. Science (New York, NY). 2002;296(5573):1634–5.
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.
Mitchell S, Vargas J, Hoffmann A. Signaling via the NFκB system. Wiley Interdiscipl Rev Syst Biol Med. 2016;8(3):227–41.
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.
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.
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.
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.
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.
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.
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.
Mauxion F, Faux C, Séraphin B. The BTG2 protein is a general activator of mRNA deadenylation. EMBO J. 2008;27(7):1039–48.
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.
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.
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.
We acknowledge GEO database for providing their platforms and contributors for uploading their meaningful datasets.
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).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Significant enrichment of GO terms for Candida albicans (top 5 according to P value).
Significant enrichment of GO terms for Candida glabrata (top 5 according to P value).
Significant enrichment of GO terms for Candida parapsilosis (top 5 according to P value).
Significant enrichment of GO terms for Candida tropicalis (top 5 according to P value).
Top 10 significantly enriched KEGG pathways of Candida albicans (according to P value).
Top 10 significantly enriched KEGG pathways of Candida glabrata (according to P value).
Top 10 significantly enriched KEGG pathways of Candida parapsilosis (according to P value).
Top 10 significantly enriched KEGG pathways of Candida tropicalis (according to P value).
Construction of the PPI network from four groups of Candida. Candida albicans (A), Candida glabrata (B), Candida parapsilosis (C) and Candida tropicalis (D).
Top 10 hub genes of four groups of Candida. Candida albicans (A), Candida glabrata (B), Candida parapsilosis (C) and Candida tropicalis (D).
About this article
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
- High-throughput sequencing
- Differentially expressed genes
- Signaling pathways
- Bioinformatics analysis