Chenglong Rao1†
Chan Mao2†Yupei Xia1†Meijuan Zhang1Zhiqiang Hu1Siqi Yuan1Wenbo Yang1Jingmin Yan1Ling Deng1Xiaolian Cai3
Xuhu Mao1
Qian Li1*Yaling Liao1*- 1Department of Clinical Microbiology and Immunology, College of Pharmacy and Medical Laboratory, Army Medical University (Third Military Medical University), Chongqing, China
- 2Department of Pharmacy, Second Affiliated Hospital, Army Medical University (Third Military Medical University), Chongqing, China
- 3Department of Cardiology, First Affiliated Hospital, Army Medical University (Third Military Medical University), Chongqing, China
Burkholderia pseudomallei is a zoonotic pathogen that usually affects patients' lungs and causes serious melioidosis. The interaction of B. pseudomallei with its hosts is complex, and cellular response to B. pseudomallei infection in humans still remains to be elucidated. In this study, transcriptomic profiling of B. pseudomallei-infected human lung epithelial A549 cells was performed to characterize the cellular response dynamics during the early infection (EI) stage. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed by using the online databases DAVID 6.8 and KOBAS 3.0. Real-time quantitative PCR and western blot were used for validation experiments. Compared with the negative control group (NC), a set of 36 common genes varied over time with a cut-off level of 1.5-fold change, and a P-value < 0.05 was identified. Bioinformatics analysis indicated that the PERK-mediated unfolded protein response (UPR) was enriched as the most noteworthy biological process category, which was enriched as a branch of UPR in the signaling pathway of protein processing in the endoplasmic reticulum. Other categories, such as inflammatory responses, cell migration, and apoptosis, were also focused. The molecular chaperone Bip (GRP78), PERK, and PERK sensor-dependent phosphorylation of eIF2α (p-eIF2α) and ATF4 were verified to be increasing over time during the EI stage, suggesting that B. pseudomallei infection activated the PERK-mediated UPR in A549 cells. Collectively, these results provide important initial insights into the intimate interaction between B. pseudomallei and lung epithelial cells, which can be further explored toward the elucidation of the cellular mechanisms of B. pseudomallei infections in humans.
Introduction
Melioidosis is a potentially fatal infectious disease in tropical and subtropical countries worldwide caused by the Gram-negative bacteria Burkholderia pseudomallei (Wiersinga et al., 2012). Humans and animals usually acquire this disease through broken skin, inhalation, or ingestion of the Tier-1 (top tier) select agent B. pseudomallei (Wiersinga et al., 2006). It is estimated that there are ~165,000 human melioidosis cases and ~89,000 deaths (54%) worldwide per year (Limmathurotsakul et al., 2016). Due to the globalization of the world's tourism and trade, the epidemic areas are expanding (Wiersinga et al., 2012, 2018; Limmathurotsakul et al., 2016). In addition, there is no licensed vaccine for the prevention of disease (Titball et al., 2017). It is particularly important and prominent to make a profound study on the pathogenic mechanism of B. pseudomallei.
B. pseudomallei is a facultative intracellular bacterium, and it can adhere to and invade a number of mammalian cell types (Jones et al., 1996; Williams et al., 2014) and persist in vivo for many years (Mays and Ricketts, 1975). During the early infection (EI) stage, B. pseudomallei can escape from the endocytic vacuole into the cytosol, next modulate different cellular responses, and evade intracellular killing in infected cells (Willcocks et al., 2016). Previous studies have shown that B. pseudomallei modulated host iron homeostasis to facilitate iron availability and intracellular survival (Schmidt et al., 2018). It could also evade LC3-associated phagocytosis (Gong et al., 2011) and induce macrophages pyroptosis (Bast et al., 2014). Our previous research demonstrated that B. pseudomallei could interfere with the progress of Rab32 GTPase-mediated B. pseudomallei-containing phagosomes maturation and escape into the cytosol (Hu et al., 2019). B. pseudomallei could evade autophagy by regulating ATG10 in A549 cells (Li et al., 2015). Moreover, B. pseudomallei acquired the evolutionary ability to subvert autophagy by hijacking host lipid metabolism for intracellular survival (Tang et al., 2020). However, the interaction of B. pseudomallei with its hosts is more complex, and cellular response to B. pseudomallei infection in humans still remains incomplete.
Transcriptome profiling is an effective strategy for understanding the molecular events in host–pathogen interactions in recent years (Tuanyok et al., 2006; Chin et al., 2010). Previous transcriptome studies about B. pseudomallei have shown that it could up- or down-regulate numerous genes to adapt rapidly to the intracellular environment in human macrophage-like U937 cells during the EI stage (Chieng et al., 2012). B. pseudomallei could up-regulate two component signal transduction systems and a denitrification enzyme pathway for biofilm production and virulence (Chin et al., 2015; Wong et al., 2015). The small colony variant (SCV) could up-regulate many virulence and survival factors pre- and post-exposed to A549 cells (Al-Maleki et al., 2020). Previous studies are limited by a narrow dynamic expression range, and little is known about the dynamic changes of host response to B. pseudomallei infection.
Human lung epithelial cells are particularly susceptible following exposure by inhalation (Zueter et al., 2016). Multinucleated giant cell (MNGC) formation is an important characteristic feature of the intracellular life cycle of B. pseudomallei during the late dissemination stage (Whiteley et al., 2017). The persistent infection of B. pseudomallei promotes the fusion of an infected mononuclear cell with one or more neighboring cells (Kespichayawattana et al., 2000; McNally and Anderson, 2011). Identifying the differentially expressed genes (DEGs) of host cells and digging their inter-relationships, such as enriched biological processes and pathways, will help to provide insight into the combined influence of B. pseudomallei on body function. Therefore, we divided the EI stage according to no more than 50% MNGC formation and dynamically investigated the global transcriptional response of human lung epithelial cell A549 infected with B. pseudomallei using microarray analysis during the EI stage to get a better understanding of genes association. Analysis of the DEGs revealed that the response genes of A549 cells to B. pseudomallei infection were mainly involved in unfolded protein response (UPR). UPR affects host's proteins synthesis and maturation, and the finding may provide a new direction for future work to reveal the pathogenic mechanism of B. pseudomallei in manipulating UPR signaling.
Materials and Methods
Cell Lines and Bacterial Strains
Human lung epithelial cell line A549 (ATCC, CCL_185) was grown in Dulbecco's Modified Eagle Medium (DMEM, Gibco, 11965-092) containing 10% fetal bovine serum (Gibco, 10099-141) without the addition of antibiotics at 37°C with an atmosphere containing 5% CO2. B. pseudomallei C006 (BPC006), as a representative of clinical virulent strain, which was the first sequencing strain in China, was used for bacterial infection in this study (Fang et al., 2012). The genome sequence of BPC006 was very close to the international standard strain K96243, and the mapped reads were ~94% (Supplementary Table 1). BPC006 was cultured in Luria-Bertani (LB) broth at 37°C with shaking at 200 rpm, and the overnight culture was washed twice in phosphate-buffered saline (PBS) and adjusted to an appropriate concentration by measurement of the optical density at 600 nm for infection. Live B. pseudomallei was handled under standard laboratory conditions (biosafety containment level 3).
Immunofluorescence and MNGC Formation
A549 cells were seeded in 12-well plates with 1 × 105 cells per well and grown on glass coverslips (NEXT, China) overnight. Subsequently, the cells were infected with B. pseudomallei at a multiplicity of infection (MOI) of 10 and incubated for 6 or 12 h. The cells were then fixed with 4% paraformaldehyde (PFA, Electron Microscopy Sciences) for 10 min and permeabilized with 0.05% Triton X-100 (Sigma, T8787). B. pseudomallei-infected fixed cells were incubated with a rabbit polyclonal antibody of B. pseudomallei (1:500) from our previous research (Hu et al., 2019) and subsequently with the Alexa Fluor 568 anti-rabbit secondary antibody (Abcam, UK) (1:2,000). The cells were stained with Actin-Tracker Green (Biotechnology, Shanghai, China). Nuclei were counterstained with DAPI (Life Technologies, 300 nM). Cover glasses were mounted on glass slides using a fluorescence mounting medium (DakoCytomation). Images were acquired on a fluorescence microscope (Zeiss, Germany) using a 40× (60× for B. pseudomallei) lens objective.
MNGCs were defined as cells containing at least three nuclei. Counting statistics of MNGCs was conducted according to the method outlined in previous reports (Whiteley et al., 2017). In brief, MNGC formation efficiency (as a percentage) was determined with a 20× objective using the following formula: (N within multinucleated giant cells / total N) × 100, where N is the number of nuclei. A minimum of 2,300 nuclei were counted for each condition. Three independent replicates were performed at each time point in the test.
Intracellular Survival of Bacteria
Intracellular survival of B. pseudomallei in A549 was estimated as previously described with some modifications (Li et al., 2015). In the initial 2 h infection with B. pseudomallei, the bacteria and cells were co-incubated for 1 h, and then the cells were overlaid with DMEM containing 250 μg/ml kanamycin for 1 h to kill extracellular bacteria. Subsequently, the cells were washed with PBS and maintained with 20 μg/ml of kanamycin in DMEM for the next 1, 4, 8, 12, 16, and 20 h. After incubation, infected A549 cells were washed three times with PBS and lysed with 1 ml of 0.1% Triton X-100 (Sigma, T8787). Cell lysates were collected and serially diluted 10-fold in PBS, and aliquots of 10 μl were plated onto LB agar to assess viable bacterial counts. We performed the experiment three replicates independently.
GeneChip™ Microarray Assay
The A549 RNA extracted from 6 to 12 h post-infection (hpi) was isolated by means of the RNeasy Micro Kit (Qiagen) according to the manufacturer's protocol for transcriptome analyses using GeneChip Human Gene 2.0 ST Arrays (Affymetrix). The control wells at 0 hpi without B. pseudomallei were used as reference group for normalization. Sample preparation for microarray hybridization was carried out as described in the Affymetrix GeneChip® Whole Transcript (WT) Sense Target Labeling Assay manual (Affymetrix, Inc., Santa Clara, CA, USA). In brief, 300 ng of total RNA was used to generate double-stranded cDNA. First, cRNA was synthesized (WT cDNA Synthesis and Amplification Kit, Affymetrix), purified, and reverse transcribed into single-stranded DNA. Purified ssDNA was then fragmented and labeled with biotin (WT Terminal Labeling Kit, Affymetrix). Finally, 2.3 μg DNA was hybridized to GeneChip Human Gene 2.0 ST Arrays (Affymetrix) for 16 h at 45°C in a rotating chamber. Hybridized arrays were washed and stained in the Affymetrix Washing Station FS450 using Hyb, Wash & Stain Kit (Affymetrix), and the fluorescent signals were measured in the Affymetrix GeneChip® Scanner 3000-7G. Sample processing was performed at the Affymetrix Service Provider and Core Facility, “KFB-Center of Excellence for Fluorescent Bioanalytics” (Regensburg, Germany; http://www.kfb-regensburg.de). For microarray data analysis, the RMA algorithm in the Affymetrix GeneChip Expression Console software 4.0.1 was used to summarize probe signals. Then, they were exported to Microsoft Excel, and average signal values, comparison fold changes, and significance P-values were calculated. Probe sets with a fold change above 1.5-fold and a Student's t-test P-value lower than 0.05 were considered as significantly regulated according to the criteria in a previous report (Manalo et al., 2005).
Functional Annotation and Enrichment Analysis of DEGs
Gene Ontology (GO) is a functional annotation of the DEGs. The mRNAs exhibiting differential expression were entered into the integrated discovery online server (http://david.abcc.ncifcrf.gov) to be annotated and visualized. The analyses included classifications of cell constituents and molecular function, as well as biological processes, with a confidence level of 95%. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis using the KEGG Orthology-Based Annotation System (KOBAS; http://kobas.cbi.pku.edu.cn/home.do) was also conducted, and the signaling pathways were verified using the Fisher's test and an false discovery rate (FDR) of ≤0.05. Significant DEGs were used to populate human pathways in KOBAS, with statistical significance determined by Fisher's exact test or hypergeometric test.
Validation of Microarray Assay
To validate the data generated from GO analysis and pathway analysis, 11 DEGs were selected from different functional categories for real-time quantitative PCR (RT-qPCR) analysis. Supplementary Table 2 shows the respective primer sequences of selected mRNA transcripts. β-Actin and GAPDH were used as reference genes for normalization. The total RNA used for RT-PCR was extracted by TRIzol (Invitrogen Life Technologies) according to the manufacturer's recommendation. Next, complementary DNA (cDNA) was retro-transcribed from 1 μg of total RNA using PrimeScript™ RT reagent kit with gDNA Eraser (Takara, Dalian, China). CFX96 Touch Real-Time PCR detection system (Bio-Rad, USA) was used for PCR detection. The reaction system (total volume: 25 μl) consisted of the following: 12.5 μl SYBR Premix Ex Taq II (Takara, Dalian, China), 1 μl 10 mmol/L upstream primer and 1 μl 10 mmol/L downstream primer, 2 μl of ~50 ng/μl cDNA template, and supplemented with ddH2O to a final reaction volume of 25 μl. The PCR program was set as 95°C for 2 min and then 40 cycles of 95°C for 5 s and 60°C for 30 s. To confirm that only one PCR product was amplified, the PCR product was subjected to dissociation curve analysis at the end of each PCR reaction. Each assay was conducted in three replicates. The method of 2−ΔΔCT was used to calculate the relative expression levels of mRNAs in cells after B. pseudomallei infection, which were expressed as the relative fold change in the expression level in infected cells divided by that in control cells.
Western Blot Analysis
The cells from 6 to 12 hpi were collected and lysed in RIPA lysis buffer (50 mM Tris pH 7.4, 150 mM NaCl, 1% Triton X-100, 1% sodium deoxycholate, 0.1% SDS) containing protease and phosphatase inhibitor on ice for 30 min. Cell suspensions were centrifuged at 12,000g for 15 min at 4°C, and then cell supernatants were collected and incubated at 100°C for 10 min. Protein concentration was determined by BCA Protein Assay according to the instructions of the supplier (Beyotime, Beijing, China). Equal amounts of protein samples (40 μg) were loaded on 10% SDS-PAGE and electrophoretically transferred onto PVDF membranes (Roche, Switzerland). A commercial protein marker was used for identification of protein size. The membranes were blocked with 5% (w/v) nonfat dry milk for 2 h, followed by an overnight incubation at 4°C with BiP (CST 3177S, 1:1,000), PERK (CST 12703, 1:1,000), Phospho-eIF2α (Ser51, CST 3398, 1:1,000), eIF2α (CST 5324S, 1:1,000), and ATF4 (CST 11815, 1:1,000) antibody, respectively, in which p-eIF2α was an activated form of eIF2α in PERK-mediated UPR. Membranes were subsequently incubated with HRP-linked anti-rabbit IgG secondary antibody (CST 7074S, 1:5,000) at the room temperature of 25°C for 2 h. The bound antibodies were monitored with chemiluminescence using an electrogenerated chemiluminescence detection system (ChemiDoc XRS System, Bio-Rad, USA), and the density of each protein band was quantified using Image Lab software 6.0 (Bio-Rad Laboratories). GAPDH (CST 8884, 1:1,000) and β-actin (CST 3700, 1:5,000) were used as loading controls, and three independent replicates were performed.
Statistical Analysis
Statistical analysis of gene expression was performed by means of Student's t-test. SPSS software for Windows (IBM, V.20.0) was used. Statistical significance was stated in case of P-values being lower than 0.05.
Results
MNGC Formation
To observe the cell morphology, cytoskeleton protein actin was stained with Actin-Tracker Green, and confocal microscope was employed. A time series assessment of A549 phenotype challenged with B. pseudomallei was shown in Figure 1A. Three independent replicates were performed for each time point in the test. There were no observable microscopic reaction lesions at 0 hpi. By 6 hpi, MNGCs had appeared with a percentage of 5.67 ± 1.15, and most of the infected cells did not fuse at this period, maintaining normal cell morphology with clear cell boundary and intact outline. However, at 12 hpi, the MNGCs had developed much larger, and the percentage increased up to 53.33 ± 4.935 (Figure 1B). Actually, the infection of B. pseudomallei was a continuous process. The early and late stages of B. pseudomallei infection were not strictly defined. We divide the early stage of infection based on the formation of MNGCs no more than 50% for the convenience of research. Therefore, the results indicated that the early stage of infection could be determined before 12 hpi.
Figure 1. MNGC formation and B. pseudomallei intracellular growth in A549 cells. (A) Confocal microscope shows the formation of MNGCs induced by B. pseudomallei at different infected times. Red is the B. pseudomallei indicated by Alexa Fluor 568, blue is the nucleus indicated by DAPI, and green is the cytoskeleton protein actin indicated by Actin-Tracker Green. (B) The percentage of MNGC formation was shown at 0, 6, and 12 hpi. (C) Intracellular growth of B. pseudomallei in A549 cell lines was assayed at different time points (1–20 h) following the initial 2 h infection with a MOI of 10. Three replicates were used for each time point. Bars represent the standard deviation of three replicates.
B. pseudomallei Intracellular Growth
Taking the infected time as the abscissa axis and the number of viable bacteria as the ordinate axis, growth kinetics curve was drawn in Figure 1C. The ability of B. pseudomallei to survive and replicate intracellularly demonstrated a slow increase from 1 to 12 hpi, with a growth rate of about 9 × 104 cfu/h. However, during 12–20 hpi, the growth of intracellular B. pseudomallei was rapid, with a high rate of about 4 × 105 cfu/h. Combined with intracellular B. pseudomallei in Figure 1A, the intracellular number at 6 hpi was much lower than that at 12 hpi. It should be mentioned that the morphology of the majority of infected cells at 6 hpi still remained normal. On the basis of all the above results, it could be deduced that the time of 6 hpi was similar to the lag phase of intracellular B. pseudomallei, that the bacteria were in the internalization stage, and that the time of 12 hpi was similar to the intracellular multiplication stage. We further choose the two time points for subsequent dynamic transcriptome research to appropriately reveal transcriptomic profiles of A549 with B. pseudomallei infection during the EI stage.
Microarray Analyses of the EI Stage
To investigate the host cellular factor involved in B. pseudomallei infection, we analyzed the mRNA by means of Affymetrix Human 2.0 Gene Chips at two time points (6 and 12 hpi) during the early stage of infection. Cells without infection (0 hpi) were used as the reference group for normalization. Applying a cut-off level of 1.5-fold change and a P-value < 0.05 on the mRNA level, we found that the cellular transcriptome profiles of duplicates were highly similar but diverged between time points, and that the number of DEGs at 6 hpi was greater than that at 12 hpi (Figure 2A). Figure 2B shows that the venn diagrams depict the overlap of DEGs from the time points of 6 and 12 hpi during the early stage of infection. A set of 306 genes at 6 hpi (Supplementary Table 3) and 103 genes at 12 hpi (Supplementary Table 2) compared with the uninfected groups (0 hpi) were observed. Thirty-six genes were involved in all two categories that indicated the 36 DEGs varied over time during the EI stage. The hierarchical clustering was performed for 306 genes (93 genes up-regulated and 213 genes down-regulated) at 6 hpi and 103 genes (45 genes up-regulated and 58 genes down-regulated) at 12 hpi, using the Pearson correlation as a similarity matrix. The heat-maps were generated to visualize the quantitative differences in the expression levels of DEGs at 6 and 12 hpi compared with 0 hpi. Clustered samples in the columns and the highly up-regulated (red) and down-regulated (blue) genes in the rows were shown in Figures 2C,D, respectively.
Figure 2. Differentially expressed genes (DEGs) at 6 and 12 hpi. (A) Changes of B. pseudomallei regulated cellular genes over time compared with NC (0 hpi) during the early stage of infection. Counts of genes up-regulated or down-regulated are shown in red or light blue dots, respectively, and a cut-off level of 1.5-fold change and a P-value < 0.05 were applied. (B) Venn diagrams depict the overlap of DEGs from each pairwise comparison between the time points, 6/0 hpi and 12/0 hpi. (C,D) Hierarchical clustering of DEGs at 6 and 12 hpi. Genes with high expression level are in red; genes with low expression level are in blue.
Cellular Processes and Pathways Regulated by B. pseudomallei Infection
To explore cellular responses to B. pseudomallei infection and characterize DEG functions, we predicted enriched biological processes for up-regulated gene set and down-regulated gene set against the GO database, respectively, as a function of time following B. pseudomallei infection (Figure 3). No matter it was at 6 or 12 hpi, the most noteworthy enrichment biological process category was the PERK-mediated UPR with P-value = 2.48 × 10−5 at 6 hpi and P-value = 2.75 × 10−4 at 12 hpi in the up-regulated gene set, including three common DEGs of EIF2S1, HSPA5, also known as GRP78 or BiP, and ATF4. The expression levels of these three DEGs at later time points were comparable to or higher than those at earlier time points (Table 1), suggesting that PERK-mediated UPR was induced and became stronger over time during the EI stage.
Figure 3. Cellular responses to B. pseudomallei infection analysis. Gene Ontology enrichment analysis was performed on up-regulated gene set and down-regulated gene set at 6 and 12 hpi, respectively. The significantly enriched Gene Ontology processes were ranked by their enrichment scores, and the top 20 processes were shown for each gene set. Blue and red panels represent down-regulated and up-regulated gene sets, respectively.
Table 1. Differentially expressed genes of A549 infected with Burkholderia pseudomallei at 6 and 12 hpi in contrast to control cells (0 hpi).
Certainly, many enrichment biological processes of DEGs in both up- and down-regulated gene sets pointed to inflammatory response, which included the processes of cellular response to tumor necrosis factor (TNF), cellular response to lipopolysaccharide, positive regulation of interleukin-8 and interleukin-6 production, and cellular response to interleukin-1. The DEGs of SPINK5, MAP2K6, IL20RB, IL1A, and ALOX5AP down-regulated over time (↓) and THBS1, TNFAIP3, IL6, GKN2, F3, and NFκB2 up-regulated over time (↑) were involved in this category. Other enrichment cellular processes that varied over time were cell migration including the DEGs of ANXA9↓, CLDN7↓, ICAM1↑, UBASH3B↑, SAMD9L↓, and CTGF↑ and apoptosis including the DEGs of GADD45A↑, FOSL1↑, MT1X↑, BCL2L15↓, and IGFBP3↑. The functional category and DEGs that varied over time were listed in Table 1.
Moreover, we performed KEGG pathway classification for the DEGs in the early stage of infection, and it showed that the 10 top pathways were activated (Supplementary Table 5), including the signaling pathway of protein processing in the endoplasmic reticulum (ER) (Figure 4A). Differentially expressed proteins of BiP (also known as HSPA5/GRP78), eIF2α, ATF4, Hsp70 (including HSPA1A, HSPA1B, and HSPA2), HERPUD1, and EDEM2 were enriched in this signaling pathway, as shown in Figure 4B. What caught our attention was that the signaling pathway of protein processing in the ER was mainly induced by the branch of PERK-mediated UPR. The other top enrichment pathways were mainly involved in inflammatory response, including TNF signaling pathway and some viral-like pathogens infection signaling pathways, which were consistent with the results of GO enrichment analysis. In general, our results inferred that B. pseudomallei triggered the integration response of the host cell, but PERK-mediated UPR in the pathogenesis of B. pseudomallei should be paid more attention.
Figure 4. Pathway analysis of DEGs from GeneChip™ microarray assay. (A) The 10 most enriched KEGG pathways were displayed in scatter charts. Different colors indicate different P-value ranges. Important pathways were outlined in red bubbles. The different sizes of the dot indicate the number of target genes in the pathway. The pathway of protein processing in the endoplasmic reticulum was significantly enriched, P < 0.001. (B) Comparison of single gene expression changes in the pathway of protein processing in the endoplasmic reticulum. DEGs were colored in red (up-regulated) or green (down-regulated), and white or gray represents no significant difference in gene expression. Each colored box was divided into three parts, and the left, middle, and right parts represent the expression state of cells without infection, 6, and 12 hpi, respectively.
Validation of Microarray Data and UPR
To validate the expression results observed in microarray, we further randomly selected 11 genes for verification by means of RT-qPCR. Selected genes were, respectively, involved in PERK-mediated UPR (eIF2α and ATF4), inflammatory responses (IL1A, IL6, TNFAIP3, NFκB2, and GKN2), cell migration (ANXA9, CLDN7, and ICAM1), and apoptosis (BCL2L15) in Table 1. As summarized in Figure 5, our RT-qPCR experiments indicated that the mRNA expression level of the 11 selected genes was consistent with the expression of microarray (Figures 5A–K), with the Pearson correlation r = 0.93 and P-value < 0.0001 (Figure 5L), which indicate that the microarray data were reliable.
Figure 5. Expression level of host genes in A549 cells verified using RT-qPCR. (A–K) The relative expression level of each mRNA transcript was represented as the n-fold change relative to NC (0 hpi). EIF2α and ATF4 were verified in PERK-mediated unfolded protein response; IL1A, IL6, TNFAIP3, NFκB2, and GKN2 were verified in inflammatory responses; ANXA9, CLDN7, and ICAM1 were verified in cell migration; and BCL2L15 were verified in apoptosis. All data are presented as mean ± SD (n = 3 in each time point). (L) The correlation of expression level between GeneChip™ microarray and RT-qPCR was high and significant (Pearson's r = 0.93; P < 0.0001).
The PERK-mediated UPR is one branch of UPR in the signaling pathway of protein processing in the ER (Bettigole and Glimcher, 2015). It will be activated to maintain protein-folding homeostasis in the ER. Concomitantly, mRNA translation was transiently attenuated through the phosphorylation of eIF2α leading to ATF4 activation (Hotamisligil, 2010). To investigate B. pseudomallei-activated PERK arm, BiP, PERK, eIF2α/p-eIF2α, and ATF4 were selected for verification via western blotting analysis, and the result was shown in Figure 6. With B. pseudomallei infection prolonged in the EI stage, the expression level of the ER chaperone proteins BiP, PERK, and ATF4 was verified to be significantly increased. The relative expression of BiP, PERK, and ATF4 was 2.20, 3.82, and 1.68 at 6 hpi, respectively, and then it went up to 2.58, 7.60, and 3.38 at 12 hpi, respectively. There was no significant change about the expression level of eIF2α, but the phosphorylation of eIF2α, the key to the activation of PERK-mediated UPR (Pavio et al., 2003), was verified to be significantly increased with the relative expression from 1.40 to 1.74. Our present data suggested that B. pseudomallei infection induced PERK-mediated UPR in human lung epithelial cells.
Figure 6. Expression of BiP (GRP78), PERK, eIF2α/p-eIF2α, and ATF4 was increased in B. pseudomallei-infected A549 cells. (A) Western blot analysis of DEGs involved in UPR challenged with B. pseudomallei. GAPDH and β-actin were used as loading controls, and all experiments were performed three times. (B) The relative expression level of DEGs in every group. Cells without infection (0 hpi) was the control group. A t-test was performed to compare the means, and the asterisks represent the P-value as follows * <0.05, n.s = non-significant, and P-value > 0.05 is considered not statistically significant.
Discussion
The intracellular lifestyle of B. pseudomallei in host cells has been reported to impact the severity of melioidosis, ranging from acute fatal sepsis to chronic infection with or without clinical symptoms (Pilatz et al., 2006; Lazar Adler et al., 2009; Allwood et al., 2011). However, the pathogenic mechanism of B. pseudomallei is still not fully understood. Elucidating the cellular response to B. pseudomallei infection will aid in unraveling the pathogenic mechanisms of the pathogen. We previously reported that B. pseudomallei could modulate autophagy (Li et al., 2015) and Rab32-mediated phagocytic vesicle transport to evade intracellular killing in host cells (Hu et al., 2019). Although previous studies have been carried out to determine the host transcriptional profiling of B. pseudomallei infection (Perumal Samy et al., 2015; Aschenbroich et al., 2019; Al-Maleki et al., 2020), limited by a narrow dynamic expression range, little is known about the dynamic changes of host response to B. pseudomallei infection. To further understand the cellular response dynamics caused by B. pseudomallei, we conducted a transcriptomic profiling of B. pseudomallei-infected A549 cells at 6 and 12 hpi during the EI stage. Our results generated a number of interesting observations. Upon B. pseudomallei infection, a set of 36 common genes that varied over time from 6 to 12 hpi were found. These DEGs are consistent with previous studies on the infection of B. pseudomallei, which involved inflammatory response (Wongprompitak et al., 2009; Perumal Samy et al., 2015), cell migration (Williams et al., 2014), apoptotic process (Hseu et al., 2014), and so on. The results indicated that these continuously regulated DEGs might be a sign of B. pseudomallei infection and played a key role in pathogenicity. Interestingly, the PERK-mediated UPR was enriched as the most noteworthy biological process category.
UPR is a cytoprotective response that prevents cytotoxicity effects caused by cellular accumulation of unfolded or misfolded proteins, which also invokes innate immune signaling in response to invading microorganisms (Celli and Tsolis, 2015). It is reported that various intracellular bacteria trigger the UPR to their advantage of replication during infection, such as Mycobacterium tuberculosis (Seimon et al., 2010), Brucella melitensis (Smith et al., 2013), and Listeria monocytogenes (Pillich et al., 2012). In our study, the expression of BiP, an ER chaperone protein, was increased as B. pseudomallei infection progresses, and B. pseudomallei infection increased PERK expression (Figure 6). This might be due to the accumulation of unfolded protein in the ER, resulting in cells needing more BiP to compensate for the dissociation of BiP from the transmembrane protein located in the ER, thereby helping the unfolding and assembly of protein structures (Bertolotti et al., 2000). Activated PERK could phosphorylate its target of the translation initiation factor eIF2α (Harding et al., 1999). Consistent with previous results, we confirmed that p-eIF2α increased with the time of B. pseudomallei infection. The p-eIF2α induced the activation of ATF4, which led to the transient weakening of gene translation involved in unfolded protein reaction (Choi and Song, 2019). The results showed that B. pseudomallei infection increased ATF4 expression (Figure 6). In conclusion, our present data suggested that B. pseudomallei induced PERK-mediated UPR in human lung epithelial cells. It was reported that UPR was exploited by energy and metabolically deficient pathogens to increase ATP and metabolites required for development, certain microbial agents, including viruses and bacteria, induce UPR (Celli and Tsolis, 2015), and it may be required for Brucella replication (Smith et al., 2013). Whether B. pseudomallei could hijack PERK-mediated UPR for intracellular survival and replication remains to be delved deeper. Moreover, the expression levels of up-regulated genes at later time points were comparable to or higher than those at earlier time points (Figure 6), suggesting that PERK-mediated UPR to B. pseudomallei infection became stronger over time. Although the UPR plays a major role in certain microbial infectivity, its role in B. pseudomallei pathogenesis is unknown. Other enrichment cellular processes, such as inflammatory responses, cell migration, and apoptosis, were enriched too. In general, our results inferred that B. pseudomallei triggered an integrated host cell response and the role of the UPR should be paid more attention to in the pathogenesis of B. pseudomallei.
Inflammatory response was also enriched according to microarray assay, as seen in Figures 3, 4. Microarray technology for transcriptome yields a large amount of data, and it is important to validate differential expression by independent methods, which conformed to the results of microarray assay. RT-qPCR analysis of DEGs included in inflammatory response, such as IL1A, IL6, TNFAIP3, NFκB2, and GKN2 (Figure 5), was therefore conducted. B. pseudomallei increased the mRNA expression of TNFAIP3 and IL-6 in human lung epithelial cell, which was consistent with the current research that proinflammatory mediators of TNF-α and IL-6 increased in the liver of mice following infection with B. pseudomallei and in serum samples of melioidosis patients (Ulett et al., 2000; Koosakulnirand et al., 2018). Activated UPR transcription factors can induce inflammatory response, which play an important role in the pathogenesis of inflammation (Bettigole and Glimcher, 2015). Previous studies showed that intracellular bacteria, such as M. tuberculosis, could induce p38 MAPK activation by TLR2 and TLR4, which in turn caused the accumulation of misfolded or unfolded TNF-α in the ER, and further induces the induction of macrophage ER stress (Oh et al., 2018; Choi and Song, 2019). Similarly, the MAPK signaling pathway was found in the process of B. pseudomallei infection in our manuscript (Figure 4A), and this might also be a possible way for B. pseudomallei to activate inflammatory response.
Additionally, our data yield that B. pseudomallei exposure could active the p-eIF2α/ATF4 pathway that was verified by western blotting in Figure 6. Accumulating evidences suggested that the p-eIF2α/ATF4 pathway mediated the induction of a gene expression program that referred to integrated stress response (ISR). ISR is a common adaptive pathway to restore cellular homeostasis, which could be triggered by integrating various types of stress signals, including ER stress, amino acid deprivation, virus infection, and oxidative stress (Harding et al., 1999; Zhang et al., 2002; Yerlikaya et al., 2008), while it is not clear whether this response benefits the host or the pathogen. However, considering the steady growth of intracellular B. pseudomallei at the early stage of infection in Figure 1C, it is likely that B. pseudomallei has evolved to modulate the ISR to its advantage during infection, but it still needs sufficient experimental evidence to reveal to what extent B. pseudomallei uses to exploit this intracellular niche to promote their replication. Additionally, the current evidence that bacteria may regulate UPR comes from the research on bacterial toxins. For instance, subtilase SubAB from Shiga-toxigenic Escherichia coli (STEC) can specifically degrade the ER chaperone BiP and further activated IRE1, ATF6, and PERK, the three known UPR signaling arms (Morinaga et al., 2008). B. pseudomallei has a strong virulence system (Wiersinga et al., 2006), whether UPR can be activated by it is not clear. Future work is likely to reveal to which novel mechanisms are employed by B. pseudomallei to manipulate UPR signaling.
Conclusions
In conclusion, this study has elucidated the dynamic DEG profiles following B. pseudomallei infection of A549 cells during the EI stage. Interestingly, a set of 36 common genes varied over time were found, and the PERK-mediated UPR was enriched as the most noteworthy biological process category. We have verified the DEGs with RT-qPCR and further confirmed that B. pseudomallei induced the expression of Bip (GRP78), PERK, p-eIF2α, and ATF4, the key proteins involved in the UPR. Overall, our analysis suggests that B. pseudomallei infection induced the PERK-mediated UPR in protein processing in the ER. Therefore, the data reported here are valuable resources for practical use to further characterize the mechanisms that are potentially involved in the intracellular living of B. pseudomallei.
Data Availability Statement
The datasets generated for this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
CR, CM, and YX are fully responsible for the study, including the study design, data collection, data analysis and writing. MZ, SY, WY and ZH help with experiment implementation and data collection. JY and XC provide experimental technical support. LD provides research funding management and reimbursement support. XM provides guidance and supporting contribution. QL and YL bear the fund, design, writing and critically reviewing the intellectual content of the article. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Natural Science Foundation of China (NSFC, Grant No. 31970135 and 81971907) and Army Medical University Youth Development Foundation (2020XQN06).
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.
Acknowledgments
This manuscript has been released as a pre-print at Research Square (Rao et al., 2020). We thank Gminix (Shanghai, China) for its assistance with the GeneChip™ microarray assay. Affymetrix microarray data are available as Excel files that we upload in the additional materials (Supplementary Table 6). NCBI accession numbers are not available due to a damaged CEL file.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2020.585203/full#supplementary-material
Supplementary Table 1. Analysis of transcriptome sequencing reads mapped to the K96243 genome.
Supplementary Table 2. Primer sequences for analysis of gene expression using RT-qPCR.
Supplementary Table 3. Summary of significant differentially expressed genes at 6 hpi.
Supplementary Table 4. Summary of significant differentially expressed genes at 12 hpi.
Supplementary Table 5. The enrichment of the top ten pathways.
Supplementary Table 6. Affymetrix microarray raw data.
Abbreviations
MNGCs, multinucleated giant cells; EI, early infection stage; LI, late infection stage; RT-qPCR, real-time quantitative PCR; IL-1A, interleukin-1A; IL-6, interleukin-6; TNF-α, tumor necrosis factor-alpha; p-eIF2α, phosphorylation of the alpha subunit of eukaryotic initiation factor 2; ATF4, cyclic AMP-dependent transcription factor 4; BPC006, B. pseudomallei C006 strain; DEG, differentially expressed gene; GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; WHO, World Health Organization.
References
Allwood, E. M., Devenish, R. J., Prescott, M., Adler, B., and Boyce, J. D. (2011). Strategies for intracellular survival of Burkholderia pseudomallei. Front Microbiol. 2:170. doi: 10.3389/fmicb.2011.00170
Al-Maleki, A. R., Vellasamy, K. M., Mariappan, V., Venkatraman, G., Tay, S. T., and Vadivelu, J. (2020). Transcriptome analysis of Burkholderia pseudomallei SCV reveals an association with virulence, stress resistance and intracellular persistence. Genomics 112, 501–512. doi: 10.1016/j.ygeno.2019.04.002
Aschenbroich, S. A., Lafontaine, E. R., Lopez, M. C., Baker, H. V., and Hogan, R. J. (2019). Transcriptome analysis of human monocytic cells infected with Burkholderia species and exploration of pentraxin-3 as part of the innate immune response against the organisms. BMC Med Genomics 12:127. doi: 10.1186/s12920-019-0575-7
Bast, A., Krause, K., Schmidt, I. H., Pudla, M., Brakopp, S., Hopf, V., et al. (2014). Caspase-1-dependent and -independent cell death pathways in Burkholderia pseudomallei infection of macrophages. PLoS Pathog 10:e1003986. doi: 10.1371/journal.ppat.1003986
Bertolotti, A., Zhang, Y., Hendershot, L. M., Harding, H. P., and Ron, D. (2000). Dynamic interaction of BiP and ER stress transducers in the unfolded-protein response. Nat. Cell Biol. 2, 326–332. doi: 10.1038/35014014
Bettigole, S. E., and Glimcher, L. H. (2015). Endoplasmic reticulum stress in immunity. Ann. Rev. Immunol. 33, 107–138. doi: 10.1146/annurev-immunol-032414-112116
Celli, J., and Tsolis, R. M. (2015). Bacteria, the endoplasmic reticulum and the unfolded protein response: friends or foes? Nat. Rev. Microbiol. 13, 71–82. doi: 10.1038/nrmicro3393
Chieng, S., Carreto, L., and Nathan, S. (2012). Burkholderia pseudomallei transcriptional adaptation in macrophages. BMC Gen. 13:328. doi: 10.1186/1471-2164-13-328
Chin, C. Y., Hara, Y., Ghazali, A. K., Yap, S. J., Kong, C., Wong, Y. C., et al. (2015). Global transcriptional analysis of Burkholderia pseudomallei high and low biofilm producers reveals insights into biofilm production and virulence. BMC Gen. 16:471. doi: 10.1186/s12864-015-1692-0
Chin, C. Y., Monack, D. M., and Nathan, S. (2010). Genome wide transcriptome profiling of a murine acute melioidosis model reveals new insights into how Burkholderia pseudomallei overcomes host innate immunity. BMC Gen. 11:672. doi: 10.1186/1471-2164-11-672
Choi, J. A., and Song, C. H. (2019). Insights Into the Role of Endoplasmic Reticulum Stress in Infectious Diseases. Front. Immunol. 10:3147. doi: 10.3389/fimmu.2019.03147
Fang, Y., Huang, Y., Li, Q., Chen, H., Yao, Z., Pan, J., et al. (2012). First genome sequence of a Burkholderia pseudomallei Isolate in China, strain BPC006, obtained from a melioidosis patient in Hainan. J. Bacteriol. 194, 6604–6605. doi: 10.1128/JB.01577-12
Gong, L., Cullinane, M., Treerat, P., Ramm, G., Prescott, M., Adler, B., et al. (2011). The Burkholderia pseudomallei type III secretion system and BopA are required for evasion of LC3-associated phagocytosis. PLoS ONE 6:e17852. doi: 10.1371/journal.pone.0017852
Harding, H. P., Zhang, Y., and Ron, D. (1999). Protein translation and folding are coupled by an endoplasmic-reticulum-resident kinase. Nature 397, 271–274. doi: 10.1038/16729
Hotamisligil, G. S. (2010). Endoplasmic reticulum stress and the inflammatory basis of metabolic disease. Cell 140, 900–917. doi: 10.1016/j.cell.2010.02.034
Hseu, Y. C., Sung, J. C., Shieh, B. S., and Chen, S. C. (2014). Burkholderia pseudomallei infection induces the expression of apoptosis-related genes and proteins in mouse macrophages. J. Microbiol. Immunol. Infect. 47, 394–398. doi: 10.1016/j.jmii.2013.03.011
Hu, Z. Q., Rao, C. L., Tang, M. L., Zhang, Y., Lu, X. X., Chen, J. G., et al. (2019). Rab32 GTPase, as a direct target of miR-30b/c, controls the intracellular survival of Burkholderia pseudomallei by regulating phagosome maturation. PLoS Pathog. 15:e1007879. doi: 10.1371/journal.ppat.1007879
Jones, A. L., Beveridge, T. J., and Woods, D. E. (1996). Intracellular survival of Burkholderia pseudomallei. Infect. Immun. 64, 782–790. doi: 10.1128/IAI.64.3.782-790.1996
Kespichayawattana, W., Rattanachetkul, S., Wanun, T., Utaisincharoen, P., and Sirisinha, S. (2000). Burkholderia pseudomallei induces cell fusion and actin-associated membrane protrusion: a possible mechanism for cell-to-cell spreading. Infect. Immun. 68, 5377–5384. doi: 10.1128/IAI.68.9.5377-5384.2000
Koosakulnirand, S., Phokrai, P., Jenjaroen, K., Roberts, R. A., Utaisincharoen, P., and Dunachie, S. J. (2018). Immune response to recombinant Burkholderia pseudomallei. PLoS ONE 13:e0198906. doi: 10.1371/journal.pone.0198906
Lazar Adler, N. R., Govan, B., Cullinane, M., Harper, M., Adler, B., and Boyce, J. D. (2009). The molecular and cellular basis of pathogenesis in melioidosis: how does Burkholderia pseudomallei cause disease? FEMS Microbiol. Rev. 33, 1079–1099. doi: 10.1111/j.1574-6976.2009.00189.x
Li, Q., Fang, Y., Zhu, P., Ren, C.-Y., Chen, H., Gu, J., et al. (2015). Burkholderia pseudomallei survival in lung epithelial cells benefits from miRNA-mediated suppression of ATG10. Autophagy 11, 1293–1307. doi: 10.1080/15548627.2015.1058474
Limmathurotsakul, D., Golding, N., and Dance, D. A. (2016). Predicted global distribution of Burkholderia pseudomallei and burden of melioidosis. Nat. Microbiol. 1:15008. doi: 10.1038/nmicrobiol.2015.8
Manalo, D. J., Rowan, A., Lavoie, T., Natarajan, L., Kelly, B. D., Ye, S. Q., et al. (2005). Transcriptional regulation of vascular endothelial cell responses to hypoxia by HIF-1. Blood 105, 659–669. doi: 10.1182/blood-2004-07-2958
Mays, E. E., and Ricketts, E. A. (1975). Melioidosis: recrudescence associated with bronchogenic carcinoma twenty-six years following initial geographic exposure. Chest 68, 261–263. doi: 10.1378/chest.68.2.261
McNally, A. K., and Anderson, J. M. (2011). Macrophage fusion and multinucleated giant cells of inflammation. Adv. Exp. Med. Biol. 713, 97–111. doi: 10.1007/978-94-007-0763-4_7
Morinaga, N., Yahiro, K., Matsuura, G., Moss, J., and Noda, M. (2008). Subtilase cytotoxin, produced by Shiga-toxigenic Escherichia coli, transiently inhibits protein synthesis of Vero cells via degradation of BiP and induces cell cycle arrest at G1 by downregulation of cyclin D1. Cell. Microbiol. 10, 921–929. doi: 10.1111/j.1462-5822.2007.01094.x
Oh, S. M., Lim, Y. J., Choi, J. A., Lee, J., Cho, S. N., Go, D., et al. (2018). TNF-α-mediated ER stress causes elimination of Mycobacterium fortuitum reservoirs by macrophage apoptosis. FASEB J. 32, 3993–4003. doi: 10.1096/fj.201701407R
Pavio, N., Romano, P. R., Graczyk, T. M., Feinstone, S. M., and Taylor, D. R. (2003). Protein synthesis and endoplasmic reticulum stress can be modulated by the hepatitis C virus envelope protein E2 through the eukaryotic initiation factor 2alpha kinase PERK. J. Virol. 77, 3578–3585. doi: 10.1128/JVI.77.6.3578-3585.2003
Perumal Samy, R., Manikandan, J., Pachiappan, A., Ooi, E. E., Aw, L. T., Stiles, B. G., et al. (2015). Gene Microarray Analyses of Daboia russelli russelli Daboiatoxin Treatment of THP-1 Human Macrophages Infected with Burkholderia pseudomallei. Curr. Mol. Med. 15, 961–974. doi: 10.2174/1566524016666151123114123
Pilatz, S., Breitbach, K., Hein, N., Fehlhaber, B., Schulze, J., Brenneke, B., et al. (2006). Identification of Burkholderia pseudomallei genes required for the intracellular life cycle and in vivo virulence. Infect. Immun. 74, 3576–3586. doi: 10.1128/IAI.01262-05
Pillich, H., Loose, M., Zimmer, K. P., and Chakraborty, T. (2012). Activation of the unfolded protein response by Listeria monocytogenes. Cell. Microbiol. 14, 949–964. doi: 10.1111/j.1462-5822.2012.01769.x
Rao, C., Mao, C., Xia, Y., Zhang, M., Hu, Z., Yuan, S., et al. (2020). Transcriptome analysis reveals dynamic changes of inflammation and stress responses during different infected stages with Burkholderia pseudomallei. Res. Square [Preprint]. doi: 10.21203/rs.2.22365/v1
Schmidt, I. H. E., Gildhorn, C., Böning, M. A. L., Kulow, V. A., and Steinmetz, I. (2018). Burkholderia pseudomallei modulates host iron homeostasis to facilitate iron availability and intracellular survival. PLoS Negl. Trop. Dis. 12:e0006096. doi: 10.1371/journal.pntd.0006096
Seimon, T. A., Kim, M. J., Blumenthal, A., Koo, J., Ehrt, S., Wainwright, H., et al. (2010). Induction of ER stress in macrophages of tuberculosis granulomas. PLoS ONE 5:e12772. doi: 10.1371/journal.pone.0012772
Smith, J. A., Khan, M., Magnani, D. D., Harms, J. S., Durward, M., Radhakrishnan, G. K., et al. (2013). Brucella induces an unfolded protein response via TcpB that supports intracellular replication in macrophages. PLoS pathogens 9:e1003785. doi: 10.1371/journal.ppat.1003785
Tang, M., Hu, Z., Rao, C., Chen, J., Yuan, S., Zhang, J., et al. (2020). Burkholderia pseudomallei interferes with host lipid metabolism via NR1D2-mediated PNPLA2/ATGL suppression to block autophagy-dependent inhibition of infection. Autophagy 2020, 1–16. doi: 10.1080/15548627.2020.1801270
Titball, R. W., Burtnick, M. N., Bancroft, G. J., and Brett, P. (2017). Burkholderia pseudomallei and Burkholderia mallei vaccines: Are we close to clinical trials? Vaccine 35, 5981–5989. doi: 10.1016/j.vaccine.2017.03.022
Tuanyok, A., Tom, M., Dunbar, J., and Woods, D. E. (2006). Genome-wide expression analysis of Burkholderia pseudomallei infection in a hamster model of acute melioidosis. Infect. Immun. 74, 5465–5476. doi: 10.1128/IAI.00737-06
Ulett, G. C., Ketheesan, N., and Hirst, R. G. (2000). Cytokine gene expression in innately susceptible BALB/c mice and relatively resistant C57BL/6 mice during infection with virulent Burkholderia pseudomallei. Infect. Immun. 68, 2034–2042. doi: 10.1128/IAI.68.4.2034-2042.2000
Whiteley, L., Meffert, T., Haug, M., and Weidenmaier, C. (2017). Entry, intracellular survival, and multinucleated-giant-cell-forming activity of Burkholderia pseudomallei in human primary phagocytic and nonphagocytic cells. Infect Immun. 85:e00468-17. doi: 10.1128/IAI.00468-17
Wiersinga, W. J., Currie, B. J., and Peacock, S. J. (2012). Melioidosis. N. Engl. J. Med. 367, 1035–1044. doi: 10.1056/NEJMra1204699
Wiersinga, W. J., van der Poll, T., White, N. J., Day, N. P., and Peacock, S. J. (2006). Melioidosis: insights into the pathogenicity of Burkholderia pseudomallei. Nat. Rev. Microbiol. 4, 272–282. doi: 10.1038/nrmicro1385
Wiersinga, W. J., Virk, H. S., Torres, A. G., Currie, B. J., Peacock, S. J., Dance, D. A. B., et al. (2018). Melioidosis. Nat. Rev. Dis. Prim. 4:17107. doi: 10.1038/nrdp.2017.107
Willcocks, S. J., Denman, C. C., Atkins, H. S., and Wren, B. W. (2016). Intracellular replication of the well-armed pathogen Burkholderia pseudomallei. Curr. Opin. Microbiol. 29, 94–103. doi: 10.1016/j.mib.2015.11.007
Williams, N. L., Morris, J. L., Rush, C. M., and Ketheesan, N. (2014). Migration of dendritic cells facilitates systemic dissemination of Burkholderia pseudomallei. Infect. Immun. 82, 4233–4240. doi: 10.1128/IAI.01880-14
Wong, J., Chen, Y., and Gan, Y. H. (2015). Host cytosolic glutathione sensing by a membrane histidine kinase activates the type vi secretion system in an intracellular bacterium. Cell Host Microbe 18, 38–48. doi: 10.1016/j.chom.2015.06.002
Wongprompitak, P., Sirisinha, S., and Chaiyaroj, S. C. (2009). Differential gene expression profiles of lung epithelial cells exposed to Burkholderia pseudomallei and Burkholderia thailandensis during the initial phase of infection. Asian Pacif J. Allerg. Immunol. 27, 59–70.
Yerlikaya, A., Kimball, S. R., and Stanley, B. A. (2008). Phosphorylation of eIF2alpha in response to 26S proteasome inhibition is mediated by the haem-regulated inhibitor (HRI) kinase. Biochem J. 412, 579–588. doi: 10.1042/BJ20080324
Zhang, P., McGrath, B. C., Reinert, J., Olsen, D. S., Lei, L., Gill, S., et al. (2002). The GCN2 eIF2alpha kinase is required for adaptation to amino acid deprivation in mice. Mol. Cell. Biol. 22, 6681–6688. doi: 10.1128/MCB.22.19.6681-6688.2002
Keywords: Burkholderia pseudomallei, transcriptomics, unfolded protein response (UPR), protein processing in endoplasmic reticulum, differential expression genes (DEGs)
Citation: Rao C, Mao C, Xia Y, Zhang M, Hu Z, Yuan S, Yang W, Yan J, Deng L, Cai X, Mao X, Li Q and Liao Y (2020) Transcriptome Analysis Reveals Unfolded Protein Response Was Induced During the Early Stage of Burkholderia pseudomallei Infection in A549 Cells. Front. Genet. 11:585203. doi: 10.3389/fgene.2020.585203
Received: 27 July 2020; Accepted: 09 November 2020;
Published: 08 December 2020.
Edited by:
Juan Caballero, Universidad Autónoma de Querétaro, MexicoReviewed by:
Vanitha Mariappan, University of Malaya, MalaysiaSorujsiri Chareonsudjai, Khon Kaen University, Thailand
Copyright © 2020 Rao, Mao, Xia, Zhang, Hu, Yuan, Yang, Yan, Deng, Cai, Mao, Li and Liao. 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: Qian Li, bGlxaWFuamFuZSYjeDAwMDQwOzE2My5jb20=; Yaling Liao, NDI1MzQ0NTgyJiN4MDAwNDA7cXEuY29t
†These authors have contributed equally to this work