Genetic deletion of Sphk2 confers protection against Pseudomonas aeruginosa mediated differential expression of genes related to virulent infection and inflammation in mouse lung

doi:10.21203/rs.2.10335/v1

Download PDF

Research article

Genetic deletion of Sphk2 confers protection against Pseudomonas aeruginosa mediated differential expression of genes related to virulent infection and inflammation in mouse lung

https://doi.org/10.21203/rs.2.10335/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 16 Dec, 2019

Read the published version in BMC Genomics →

You are reading this older preprint version

Read the latest preprint version →

BACKGROUND Pseudomonas aeruginosa (PA) is an opportunistic Gram-negative bacteria that causes serious life threatening and nosocomial infections including pneumonia. PA has the ability to alter host genome to facilitate its invasion thus increasing the virulence of the organism. Sphingosine-1 phosphate (S1P), a bioactive lipid is known to play a key role in facilitating infection. Sphingosine kinases (SPHK)1&2 phosphorylate sphingosine to generate S1P in mammalian cells. We reported earlier that Sphk2-/- mice offered significant protection against lung inflammation, compared to wild type (WT) animals. Therefore, we profiled the differential expression of genes between the protected group of Sphk2-/- and the wild type controls to better understand the underlying protective mechanisms Sphk2 deletion in lung inflammatory injury. Whole transcriptome shotgun sequencing (RNA seq) was performed on mouse lung tissue using NextSeq 500 sequencing system. RESULTS: Two-way ANOVA analysis was performed and differentially expressed genes following PA infection were identified using whole transcriptome of Sphk2-/- mice and their wild WT counterparts. Pathway (PW) enrichment analyses of the RNA seq data identified several signaling pathways that are likely to play a crucial role in pneumonia caused by PA such as those involved in: 1. Immune response to PA infection and NF-κB signal transduction; 2. PKC signal transduction; 3. Impact on epigenetic regulation; 4. Epithelial sodium channel pathway; 5. Mucin expression; and 6. Bacterial infection related pathways. Our data clearly suggests a role for SPHK2 in PA-induced pneumonia through elevated expression of inflammatory genes in lung tissue. Further, validation by RT-PCR on 10 differentially expressed genes showed 100% concordance in terms of vectoral changes, and 70% of genes showed concordance for significant fold change. CONCLUSION: Using Sphk2-/- mice and differential gene expression analysis, we have shown here that S1P/SPHK2 signaling plays a key role in promoting PA pneumonia. The identified genes promote inflammation, and suppress others that naturally inhibit inflammation and host defense. We therefore propose a key role for SPHK2/S1P signaling in PA-induced lung inflammation, and could serve as a potential therapeutic target to combat PA-induced pneumonia.

Epigenetics & Genomics

Pseudomonas aeruginosa

Lung infection

Sphingosine kinase 2

Sphingolipids

Gene Profiling

Resistance to Infection

Pseudomonas aeruginosa (PA) is an aggressive Gram-negative bacillus that causes serious opportunistic infections such as pneumonia in humans, leading to significant morbidity and mortality [1-3]. However, it is interesting to note that PA is also capable of causing serious infections in plants and insects with significant correlation to virulence across the species [4, 5]. Among patients, those with cystic fibrosis (CF), chronic obstructive pulmonary disease (COPD), and on mechanical ventilation are particularly prone to develop pneumonia caused by PA infection [6]. In fact, PA plays a major role in deterioration of lung function in CF patients. An organism that can even grow in water, PA has of late been recognized to be capable of altering the host genome it infects in order to facilitate its own virulence [7-10]. It is known that PA mediated pneumonia leads to a cascade of responses in the host, starting with innate immune response followed by increased reactive oxygen species (ROS) generation and differential regulation of sphingolipid metabolic pathways [11-13]. In the sphingolipid pathway, it has been noted that sphingosine, which is normally present in respiratory tract of healthier patients, is almost absent in CF patients [14]. On the contrary, ceramides generated by acid sphingomyelinase are known to accumulate in the airway epithelium of CF patients with pneumonia [13, 15]. Among sphingolipids, sphingosine-1-phosphate (S1P), synthesized from sphingosine, by sphingosine kinases (SPHK)1 and 2, is an intercellular and intracellular bioactive lipid mediator that regulates pleotropic cellular functions under normal and pathophysiological conditions. Genetic deletion of Sphk1, but not Sphk2, in mouse caused aggravation of LPS-induced lung injury, suggesting a protective role for SPHK1/S1P signaling against endotoxemia [16]. In contrast, adenoviral overexpression of SPHK2 in wild type (WT) mouse augmented LPS-induced lung injury [16, 17]. Also, our recent study demonstrated that deletion of Sphk2, but not Sphk1, ameliorated PA-induced lung inflammation and injury in mice [11]. Using Sphk1 and Sphk2 knockout mice, we decided to unravel the key pathways selectively associated with SPHK2 signaling that play a role in PA-induced pathogenesis using differential gene expression analysis.

The infection of a host by a pathogenic microorganism initiates complex cascades of events that influence both immediate and long-term outcomes. In this study, we identified a set of PA responsive genes activated in the WT mice in comparison with Sphk2-/-. Our results show that SPHK2/S1P signaling cascade mediating PA-induced pneumonia modulates signaling events related to extracellular matrix remodeling, interleukin (IL) signaling, and complement cascade in the host lung. In addition, we also noted that genetic deletion of Sphk2-/- resisted alteration of host pulmonary genome by PA infection promoting its own virulence. The objective of this study is to identify novel pathways contributing to the pathology as well as protection of PA-induced pneumonia.

Mouse experiments and Animal Care

All experiments using animals were approved by the Institutional Animal Care and Use Committee at the University of Illinois at Chicago (protocol # 15-240). Sphk2 knockout mice were originally provided by Dr. Richard Proia (National Institutes of Health, Bethesda, MD). The knockout mice were backcrossed onto the C57BL/6 background for 8 generations. The resultant mixed background of C57BL/6 strain and the original background (F8 hybrid) was used as controls and is referred to hereafter as Wild Type (WT). All in vivo experiments were carried out with age-matched (6-8 weeks) female mice. The mice were housed in the University of Illinois Animal Care Facility.

Anesthesia and euthanasia: The mice were anesthetized using Ketamine (100mg/kg) and Xylazine (5mg/kg). The animals were sacrificed and the lung tissues collected, homogenized and whole cell lysates prepared for further analysis, RNA isolation (superior lobe of right lung), and RNA-seq studies.

Preparation of Pseudomonas aeruginosa culture

The parent strain P. aeruginosa (PA 103) used for all experiments was provided by Dr. Ruxana Sadikot (Emory University, Atlanta, GA). Preparation of the cultures and determination of colony-forming units (CFU) were carried out as described previously [11]. The bacterial concentration of PA was confirmed by plating out the diluted samples on sheep blood agar plates. For in vitro experiments, bacterial preparations heat-inactivated at 60oC for 20 min were used.

Pseudomonas aeruginosa infection of mouse lung

Age and weight-matched female WT and Sphk2-/- mice were anesthetized with ketamine as per approved protocol and were administrated a single intratracheal infusion of sterile PBS or PA103 in PBS (1x106 CFU/mouse). Three mice were used for each group. After 24 h of treatment, animals were sacrificed, whole lung tissue was collected, and processed as described below.

Sample processing and RNA-seq based gene expression analyses

Lungs were perfused with phosphate buffered saline prior to harvesting from the mice, and processed immediately. Whole lung tissues were used to isolate total RNA using microRNeasy® kit (Qiagen, Maryland, Cat no. 74004). RNA samples isolated from individual animals were separately labeled, hybridized, washed/stained and scanned according to the standard WT PLUS labeling protocol recommended by the manufacturer (Thermo Fisher Scientific, Waltham, MA).

RNA Quality control

RNA concentrations and purity were determined on a NanoDrop 1000 (Invitrogen), and RNA integrity was determined on the 2200 TapeStation system using RNA ScreenTape (Agilent, Cat. No. 5067-5576). RNA integrity number (RIN) values ranged from 7.0 to 8.4.

RNA-seq Library Preparation

Libraries were prepared with the 3’ QuantSeq mRNA-Seq Library Prep Kit REV for Illumina (Lexogen), according to manufacturer's instructions. In brief, 10-500 ug of total RNA was used to make each library. Library generation was initiated by oligo (dT) priming followed by first strand cDNA synthesis, removal of RNA and second strand cDNA synthesis using random priming and DNA polymerase. During these steps Illumina linker sequences and external barcodes were incorporated. Next the libraries were subject to final 20 cycles of PCR amplification.

RNA-seq Library Validation and Quantification

Quality of the libraries was checked on the 2200 Tape Station system using D1000 Screen Tape (Agilent, Cat. No. 5067-5582), and as expected peaks ranged from 264 to 294 bp. Libraries were quantified on the Qubit 2.0 Fluorometer with the Qubit dsDNA HS Assay Kit (Life Technologies, Cat. No. Q32854). Individual libraries were pooled in equimolar amounts and concentration of the final pool was determined by PCR quantification method using KAPA Library Quantification Kit (KAPA Biosystems). Sequencing was carried out on NextSeq 500 (Illumina), 1x75 nt reads, high output, to achieve approximately 20x106 clusters per sample.

Collected hybridization signals were processed using Genomics Suite 6.6 statistical package (Partek, Inc., Saint Louis, MO). The parameters applied for hybridization signal processing were as follows: RMA algorithm-based background correction, quantile normalization procedure, and probe set summarization [18, 19].

All processed array files were inspected for the quality metrics such as average signal present, signal intensity of species-specific housekeeping genes, relative signal intensities of labeling controls, absolute signal intensities of hybridization controls, and across-array signal distribution plots [20]. All hybridizations passed quality control according to indicated labeling and hybridization controls.

Identification of differentially expressed transcripts

In order to identify the subset of genes modulated specifically to the infection of WT and Sphk2-/- mice, we performed a two-way ANOVA using the status of PA infection and Sphk2 expression as comparison factors. We compared the following groups: Sphk2-/- PA infected (Sphk2-/- PA), Sphk2-/- control (Sphk2-/- CTRL), Wild Type PA infected (WT PA) and Wild Type control (WT CTRL). ANOVA model was based on Method of Moments [21] in combination with Fisher's Least Significant Difference (LSD) contrast (Tamhane and Dunlop, 2000). The Fisher’s contrast allowed calculation of direction and magnitude of change for all pair-wise comparisons between the treatment groups and was later validated by RT-PCR. Raw reads were aligned to reference genome using Burroughs-Wheeler Aligner Maximal Exact Matches (BWA-MEM) [22]. Gene expression was quantified using FeatureCounts [23]. Differential expression statistics (fold-change and p-value) were computed using edgeR [24, 25], generalized linear models to model the effect of genotype, infection, and their interaction. We used Globus Genomics [26] for these analyses. Calculated raw p-values were adjusted for False Discovery Rate (FDR) according to Benjamini-Hochberg (BH) correction procedure [27, 28] . Significant genes were determined based on an FDR threshold of 5% (0.05), and plotted in a heatmap. Pathway enrichment analysis on differentially expressed genes was performed using the Pathway Maps database in MetaCore. The top 35 genes, based on the interaction term FDR, were plotted in a heatmap. Additionally, we compared the significantly differentially expressed (FDR < 0.05) genes based on genotype, infection, or their interaction in a Venn diagram.

Availability of data and materials

The RNAseq datasets supporting the conclusions of this article are available in the National Center for Biotechnology Information Gene Expression Omnibus repository, with unique persistent identifier of NCBI tracking system accession number. The hyperlink to the datasets is given below.

https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE121359

Pathway enrichment analyses and data visualization

We performed pathway enrichment analyses (EA) in order to identify the biological factors driving the protective effect observed in Sphk2-/- mice with PA pneumonia. Transcripts identified as differentially expressed in KO animals in response to PA infection in two-way ANOVA test (FDR cut off of 0.05) were imported into MetaCore Genomic Analyses Tool Release 6.22 (Thomson Reuters) for analyses.

Differentially expressed genes were analyzed using the “Pathway Maps” ontology and the top 50 most enriched pathways (PW) were identified. The output of analyses using the tool contained a substantial number of individual PWs that overlap by genes, representing sub-segment of the same PWs and creating redundancy. In order to reduce duplication, we clustered nodal PWs based on their gene content in order to reduce duplication. Complete linkage hierarchical clustering on the Jaccard distance between the complete set of genes in each PW was used to identify closely related individual entities. A measure of the dissimilarity between two PWs (based on their gene sets) with scales from 0 to 1 was used; ‘0’ if the sets are exactly the same, and ‘1’ if they are completely different and have no genes in common. For the purpose of biological interpretations, we considered each cluster of closely related PWs as one unit or mega pathway (dissimilarity cut off of 0.6). We combined all associated differential genes for analyzing gene interactions and creating heatmaps as shown in the Venn diagram (Figure 1) and a dendrogram (Figure 2). The heatmaps for selected mega pathways were created by plotting z-scored normalized expression levels of differentially expressed genes (FDR< 0.05) across all experimental groups (Figures 3, 4, 5, 6, 7,8, 9&10). The z-scored normalized expression level using the color key ranging from dark blue to dark red.

Realtime RT- PCR validation of microarray results

Total RNA was isolated from mouse lung homogenate using TRIzol® reagent according to the manufacturer's instructions and purified using the RNeasy® Mini Kit according to the manufacturer's protocol (Qiagen, MD, USA). Quantitative Realtime RT-PCR was done using iQ SYBR Green Supermix using iCycler by Bio-Rad, USA.18S rRNA (sense, 5'-GTAACCCGTTGAACCCCATT-3', and antisense, 5'- CCATCCAATCGGTAGTAGCG-3') was used as external control to normalize expression [29]. All primers were designed by inspection of the genes of interest using data from PrimerBank database (Harvard University). The sequence description of mouse primers used are given in Supplementary Table 1. Negative controls, consisting of reaction mixtures containing all components but the target RNA, were included with each of the RT-PCR runs. The representative PCR mixtures for each gene were run in the absence of the RT enzyme after first being cycled to 95°C for 15 min in order to ensure that amplified products did not represent genomic DNA contamination. No PCR products were observed in the absence of reverse transcription. Direct comparison of four groups such as WT control, WT PA, Sphk2-/- control and Sphk2-/- PA was done using ANOVA test as described earlier. The level of statistical significance was set at p<0.05.

All the validation studies were done in the exact same cohort used in RNA-seq studies.

Shared and differentially expressed genes in wild type and Sphk2-/- mouse lungs with or without PA exposure

Analysis of gene expression showed that 375 genes were differentially regulated by PA infection of mouse lungs compared to the corresponding uninfected control mice. Venn diagram showing the number of genes differentially regulated in the WT and Sphk2-/- mice exposed to PA based on two-way ANOVA analysis is shown in Fig. 1. under three different categories viz: 1. Sphk2 gene knockout, 2. Exposure of the mouse to PA, and 3. Interaction of Sphk2 gene knock out and PA. The intersecting areas represent the number of genes affected by the corresponding condition. The advantage with two-way ANOVA is that the third variable of interaction between the two factors is purely dependent on interaction, thus independent of the direct effect of the other two variables. Only 2 of the 375 differentially expressed genes (DEGs) could be strictly characterized as those affected solely due to the impact of genetic deletion of Sphk2 (Fig 1). It is also interesting to note that 1496 genes were affected by the interaction between deletion of Sphk2 and PA infection of the mouse lung.

Pathway enrichment (PW) analyses and elucidation of underlying biological currents

The top 50 differentially regulated pathways were identified based on the gene expression profiles and relationship between the PWs is demonstrated in Fig 2. A list of the identified fifty PWs in the order of their position on hierarchical clustering graph and details of the clustering of similar pathways are provided in Supplementary Table 2. A dissimilarity score of 0.6 was used as cut off for the selection of clustering pathways and reducing the redundancy as described in the Methods. The data were condensed to eight cluster PWs that are differentially regulated among the four groups (WT CTRL, WT PA, Sphk2-/- CTRL and Sphk2-/- PA) studied, and heat maps show differential expression of selected genes (Figs. 3-10). The first heat map provides a depiction of the top 35 differentially expressed genes between the groups (Fig.3). The following are the most significant PW clusters: 1. Immune response to PA infection (cluster 1) and NF-κB signal transduction (cluster 3), 2. PKC signal transduction (cluster 8), 3. Impact on epigenetic regulation (cluster 10), 4. Epithelial sodium channel pathway (cluster 11a), 5. Mucin expression (cluster 11b), and 6. Bacterial infection related pathway (cluster 12) (Fig. 4-9).

Differentially regulated genes in the PW clusters

Wild type (WT) mice challenged with PA showed significant upregulation of the genes related to NF-κB pathway in contrast to similarly treated Sphk2-/- mice (Fig. 4). Of the five genes encoding for the Nfκb family, there was significant increase in expression of NF-κB1, NF-κB2, Rel A and Rel B. Genes encoding for NF-κB inhibitors alpha and beta were also significantly elevated in the WT PA group compared to the rest of the groups. Interestingly the Dtl gene, which encodes for TNF superfamily member 13b (TNFRSF13B), was down-regulated in the WT PA compared to the rest of the groups. This cytokine expressed in B cell lineage acts as a potent B cell activator stimulating their proliferation and differentiation [30].

Recently, a novel role of Protein Kinase C-δ (PKC-δ) in PA-induced phosphorylation of SPHK2 and histone acetylation was demonstrated in lung epithelial cells [11]. In this context, we noted that genes in the PKC pathway have undergone significant differential regulation following infection with PA of both WT and Sphk2-/-groups. PKC family of proteins activate target proteins by promoting phosphorylation at serine and threonine amino acid residues [31]. In the WT PA group upregulation of virulence promoting members of the PKC family such as Sele (stimulating adhesion of leucocytes) [32], and Lcp2 and Nr4a1 that promote apoptosis [33] was observed. Further, a significant down-regulation of genes essential for maintenance of normal metabolic state in the WT PA group belonging to the PKC family was observed. The down-regulated genes included Prkcz, Prkaca, Prkg2, Ppp1r14a and Plcb3. These genes were maintained in an up-regulated state in the control as well as the Sphk2-/- group challenged with PA. Prkcz is a member of the PKC family, which, unlike the classical PKC isoenzymes, shows calcium and diacylglycerol (DAG) independent activation. Prkaca, protein kinase A catalytic subunit (PKA Cα) is a member of the AGC kinase family, contributing to the control of cellular processes such as glucose metabolism and cell division [34]. Prkg2 encodes Protein Kinase, CGMP-Dependent, Type II, is a crucial regulator of intestinal secretion and bone growth [35, 36]. This protein also phosphorylates and activates CFTR on the plasma membrane. Ppp1r14a, protein phosphatase 1 regulatory subunit 14A, is an inhibitor of smooth muscle myosin phosphatase. Plcb3 encodes phospholipase C beta 3 that catalyzes the production of the second messengers such as DAG and inositol 1,4,5-triphosphate from phosphatidylinositol via G-protein-linked receptor-mediated signal transduction.

Genes affecting epigenetic pathways were also noted to be differentially regulated. PA infection causes downregulation of Hist1h2bc and Runx1t1 genes. While the former gene represents the histone cluster 1, H2bc, the latter represents Runx1 translocation partner 1. This gene encodes a member of the myeloid translocation gene family that interacts with DNA-bound transcription factors leading to recruitment of a range of corepressors thus causing transcriptional repression [37].

Further, we noted differential regulation of genes coding epithelial sodium channels (ENaC) in vertebrates. Scnn1g gene coding for the γ subunit of the ENaC was significantly down-regulated in WT PA but remained up-regulated in the rest of groups. WT PA was associated with up-regulation of Furin gene that was not clearly observed in Sphk2-/- mouse lungs. Furin is a host cell factor that enhances virulence of viral infection in cultured cells [38]. In the mucin expression pathway, Erb-B2 Receptor Tyrosine Kinase 2 gene encodes a member of the epidermal growth factor (EGF) receptor family of receptor tyrosine kinases, which was down- regulated in the WT PA group. Following identification of the PWs, we validated differential expression of specific genes by their biological impact, and further validated using real time RT-PCR.

Realtime RT- PCR validation of microarray results

Table. 1 and Fig. 11 A & B show real-time RT-PCR performed on select 10 genes based on the observations made from the RNA-seq gene expression analysis. RT-PCR confirmed all the 10 genes suggesting 100% validation in terms of vectoral changes and 70% for both vectoral changes and significance in ANOVA test.

The genes identified as maximally differentially regulated are represented in the heatmap shown in Fig. 3. The genes for RT-PCR validation were selected from those maximally differentially regulated genes and identified as significant from the various differentially regulated nodal pathways as they were related to cell migration, inflammation, epigenetic regulation of genome, and IL-6 signaling. Spata13, Hdc, Kdm6b, Gos2, Fosl2, Dgat1, Odc1, Smtnl2, Clec1a and Sox17 were therefore chosen for RT-PCR validation. In addition, these genes were also related to biologically relevant findings in our animal model (Table.1).

Western blot validation of microarray results

Western blot analyses were performed to validate protein expression of genes that were modulated in for Real-time RT-PCR. Immunoblots showed increased expression of KDM6B, SPATA13 and ODC1 in the PA exposed lungs of the WT mice compared to Sphk2-/- mouse lungs. CLEC1A expression was significantly decreased in Sphk2-/- mice exposed to PA, compared to WT PA. Western blots probed with corresponding antibodies were quantified by densitometry and normalized to the corresponding total protein. This finding correlated well with the results of the microarray and real-time RT-PCR analysis.

Biological impact of Sphk2 deletion in the animal model correlates with the differential expression of genes in mouse lungs

Our recently published data from the animal model showed that deletion of Sphk2, but not Sphk1, protected mice against PA-mediated inflammatory lung injury [11] . Following intratracheal instillation of PA (1x106 CFU/animal), infiltration of PMNs in the lungs and increased protein levels in BALF were noted in WT mice, whereas these responses were significantly lower in Sphk2-/- mice. Following PA infection, BALF concentrations of the pro-inflammatory mediators IL-6 and TNF-α were significantly elevated in WT compared to Sphk2-/- mice. Further, infection of mouse lungs with PA revealed enhanced phosphorylation of PKC δ and phospho-SPHK2 immunostaining predominantly in the nucleus of lung epithelial cells, and inhibiting PKC δ or SPHK2 activity with small molecule inhibitor attenuated PA-induced H3 and H4 histone acetylation and IL-6 secretion in lung epithelial cells. These in vivo and in vitro results confirm the PW cluster analysis.

Activation of innate host defense mechanism represents one of the initial responses by the host organism to an invading pathogen. This is essential for clearance of the pathogen and limit further damage. PA is responsible for a variety of infections in the nosocomial setting. It also causes serious pulmonary infections of CF patients, and has been reported to alter genome expression in the host cell beyond inciting a strong innate immune response [39, 40]. Here, we report the genomic host response in WT mice infected with PA and the essential role of SPHK2 in the ensuing pathogenesis that is relevant to pneumonia and sepsis. Our gene expression analysis clearly revealed a link between SPHK2 and several key downstream pathways that play a role, either directly or indirectly in PA pathogenesis.

One of the earliest attempts to study transcriptional changes following PA infection was carried out in A549 cells [39], wherein interferon regulatory factor 1 (IRF-1) was noted to be activated. In the same study, genes encoding transcription factors such as IkB-𝜶 (NF-kB Inhibitor Alpha-NF-kBIA), and epithelial-specific transcription factor (ESE-1/ELF3) were up-regulated. In addition, several inflammatory response genes, including monocyte chemotactic protein MCP-1 (CCL2) and TNF-α-induced protein A20 (TNFAIP3) were up-regulated [39].The DPH1 gene, whose product catalyzes conversion of histidine residues to diphthamide in elongation factor 2 was also up-regulated. Further, genes involved in a variety of signal transduction pathways that displayed altered expression, such as RhoB GTPase (RHOB), JAK-1 (JAK1), and c-Jun (JUN) were elevated. Our in vivo results are in accordance with the above finding. For instance, expression of Nfkbia, Nfkb1, Nfkb2, Elf3, Ccl2, Tnfaip3 genes were up-regulated in WT mice treated with PA compared to other groups. The Sphk2-/- PA group showed significant degree of protection accompanied by lack of activation of the aforementioned genes belonging to NF-kB activation pathway (Fig. 4). However, no significant changes were noted with respect to the expression of Dph-1, Rhob, Jak1 or c-Jun genes in WT PA group.

PKC family of enzymes are involved in a variety of cellular processes including proliferation, differentiation and secretion. Ichikawa et al also reported PKC family genes to be up-regulated by PA in epithelial cells (39). In another study, it has been observed that PA activates PKC α to invade middle ear epithelial cells [41]. Activation of PKC α by its phosphorylation in human middle ear epithelial cells (HMEEC) was associated with actin condensation, and blocking the PKC pathway attenuated the ability of bacteria to invade HMEECs and subsequent actin condensation. In our study we noted significant differential regulation of genes in the PKC pathway. A similarity between the two studies was noted in that infection by PA in WT mouse caused upregulation of genes promoting inflammation and apoptosis. PKC δ (Prkcd) was significantly up-regulated in WT PA compared to the rest of the group. Also, upregulation of other members of the PKC family such as Sele, which promotes adhesion of leucocytes [32], Lcp2 the product of which mediates T-cell receptor mediated signal transduction [42], Nr4a1 the product of which upon translocation from nucleus to mitochondria causes apoptosis [33] was observed. However, PA infection leads to significant downregulation of genes in the PKC pathway including Gna11 encoding the protein belonging to the family of guanine nucleotide-binding proteins (G proteins), which function as modulators or transducers in various transmembrane signaling systems and Prkcz, encoding protein kinase c , a member of the PKC family of serine/threonine kinases. In the present study, we also observed downregulation of RAS Guanyl Releasing Protein 2 gene expression and the corresponding protein activates small GTPases, including RAS and RAP1/RAS3 [43, 44].

CF patients are highly susceptible to PA infection including the chronic colonization and biofilm formation, culminating in pneumonia [45]. In order to understand the pathogenesis of infection and colonization, it is essential to understand the underlying molecular mechanisms, especially at the genomic level. Increased transcription of the mucin genes (Muc2, Muc5a) in the respiratory tract in response to the presence of bacteria followed by accumulation of viscous mucus in the airways has been described [46]. However, we did not note any differential regulation in the expression of Muc2, and Muc5a genes in our acute murine model of PA infection.

Recent studies have shed light on the role played by bacterial pathogens in reprograming host genome by influencing epigenetic factors [47]. Gene modification at epigenetic level holds the potential to alter the response of the host to future infections. PA has been shown to induce early T3SS-dependent dephosphorylation and deacetylation of histone H3 in eukaryotic cells. Epigenetic regulation of genes following PA infection is an interesting finding of ours in that it can have far reaching implications beyond immediate infection. PA infection caused downregulation of Histone1H2bc (Hist1h2bc) gene, which encodes a member of the histone 1H2B family. Histone1 family of proteins interact with linker DNA between nucleosomes aiding the compaction of chromatin into higher order structures. This protein, which has antibacterial and antifungal activity, was significantly down-regulated in the WT mice but preserved in the Sphk2-/-. This could contribute to the protection seen in the Sphk2-/-. Downregulation of Runx1T1 along with that of Hist1h2bc is important because Runx1 translocation partner 1 belongs to a family of transcriptional corepressors that interact with both transcription factors bound to promoters of target genes and with histone deacetylases (HDACs). Sphk2-/- mice infected with PA did not reveal any downregulation of the aforementioned genes in the epigenetic pathway. Interestingly, PA infection of WT mouse lung inhibited HDAC1/2 activity, and enhanced H3 and H4 histone acetylation, however, genetic deletion of Sphk2 in mice attenuated PA mediated H3 and H4 histone acetylation [11] . Further, S1P generated in the nucleus of lung epithelial cells by SPHK2 after exposure to heat-inactivated PA modulated HDAC1/2 activity that was blocked by inhibition of SPHK2 activity with a specific inhibitor, ABC 294640 [11], thus confirming the microarray and PW analysis data.

Activation of the Epithelial Sodium Channel (ENaC) by Alkaline Protease in PA treated cells has been reported [48]. Our study did not show any differential regulation of ENaC genes per se but members of the EnaC pathway genes, for example, the epithelial sodium channel Furin and endothelin1 (End 1) were up-regulated in the WT PA group. Furin plays an important role in the activation of exotoxin A, which is a major virulence factor of PA [49, 50]. Endocytosed exotoxin A is processed by Furin and is transported retrograde to the endoplasmic reticulum. The toxin activates Pseudomonas exotoxin A by specific cleavage and inhibits protein synthesis by ADP ribosylation of elongation factor 2, triggering cell death [51]. As noted earlier, WT PA group was associated with up-regulation of Furin gene unlike similarly treated Sphk2-/- mice (Fig. 7). This indicates that PA-induced activation of Furin is somehow dependent on SphK2 gene activity and could be the principal mechanism through which host combats PA infection. Further studies are in progress to validate the above hypothesis.

This study sheds light on the key role played by SPHK2 in facilitating the PA infection in the host and how the ‘hijacking’ of the host genome by the invading organism can be resisted by blocking the SPHK2/S1P pathway. PA infection caused significant upregulation of the genes related to NF-κB pathway and members of the PKC family. Also, PA infection up-regulated Furin gene, which plays a critical role in the activation of exotoxin A, a major virulence factor of PA. Translational application of the SPHK2 pathway can have far reaching implications in the therapy of PA infection in CF patients by two mechanisms such as 1) Targeting the SPHK2 pathway mediated pathology and 2) Reducing the virulence of the organism by resisting its alteration of the host genome. Further understanding of this novel mechanism can be harnessed to overcome the virulence of the organism thus opening up new avenues in the therapy of PA infections across various species.

PA: Pseudomonas aeruginosa

S1P: Sphingosine-1 phosphate

SPHK: Sphingosine kinases

WT: Wild type

ANOVA: Analysis of variance

PW: Pathway

CF: Cystic fibrosis

CFTR: Cystic fibrosis transmembrane regulator

COPD: Chronic obstructive pulmonary disease

CFU: Colony-forming units

RIN: RNA integrity number

RT-PCR: Reverse transcriptase polymerase chain reaction

DEG: Differentially expressed genes

CTRL: Control

CFU: Colony-forming units

EA: Enrichment analyses

TNFRSF13B: TNF superfamily member 13b

PKC: Protein Kinase C

ENaC: Epithelial sodium channels

EGF: Epidermal growth factor

IRF-1: Interferon regulatory factor 1

ESE-1/ELF3: Epithelial-specific transcription factor

MCP: Monocyte chemotactic protein

TNFAIP3: TNF-α-induced protein A20

RHOB: RhoB GTPase

HMEEC: Human middle ear epithelial cells

Muc: Mucin

His: Histone

HDAC: Histone deacetylases

- Ethics approval and consent to participate: All the experiments using animals were approved by the Institutional Animal Care and Use Committee at the University of Illinois at Chicago (protocol # 15-240).

- Consent for publication: Yes, the authors give consent to publish.

- Availability of data and material: The RNAseq datasets supporting the conclusions of this article are available in the National Center for Biotechnology Information Gene Expression Omnibus repository, with unique persistent identifier of NCBI tracking system accession number. The hyperlink to the datasets is given below.

https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE121359

- Competing interests: We confirm that we have read BioMed Central's guidance on competing interests and state that none of the authors have any competing interests.

- Funding: This work was supported in part by R01HD090887-01A1 from NICHD to AH and Program Project Grant P01 HL098050 from NIH/HLBI to VN.

- Authors' contributions: DLE conceptualized the study and its design, performed experiments, interpreted results and edited the manuscript. PF performed experiments, interpreted results, prepared figures and edited the manuscript. YK performed experiments, interpreted results, prepared figures and edited the manuscript. MMC analyzed data, interpreted results, prepared figures and drafted the manuscript. HH analyzed data, interpreted results, prepared figures and drafted the manuscript. SJ analyzed data, interpreted results, prepared figures and drafted the manuscript. RM analyzed data, interpreted results, prepared figures and edited the manuscript. ZA analyzed data, interpreted results, prepared figures, drafted, deposited the data in the NIH NCBI biorepository and edited the manuscript. AH conceptualized the study and its design, performed experiments, analyzed data, interpreted results, drafted and edited the manuscript. VN conceptualized the study and its design, interpreted results and edited the manuscript. All the authors read and approved the final manuscript.

- Acknowledgements: We gratefully acknowledge the help of Prof. Prasad Kanteti, PhD, in editing the manuscript and that of Ms. Uma Harijith, MA, and Ms. Alison Ha, BS, in the preparation of this manuscript.

1. Aaron SD, Vandemheen KL, Ramotar K, Giesbrecht-Lewis T, Tullis E, Freitag A, Paterson N, Jackson M, Lougheed MD, Dowson C et al: Infection with transmissible strains of Pseudomonas aeruginosa and clinical outcomes in adults with cystic fibrosis. JAMA 2010, 304(19):2145-2153.

2. Arancibia F, Bauer TT, Ewig S, Mensa J, Gonzalez J, Niederman MS, Torres A: Community-acquired pneumonia due to gram-negative bacteria and pseudomonas aeruginosa: incidence, risk, and prognosis. Arch Intern Med 2002, 162(16):1849-1858.

3. Parkins MD, Somayaji R, Waters VJ: Epidemiology, Biology, and Impact of Clonal Pseudomonas aeruginosa Infections in Cystic Fibrosis. Clin Microbiol Rev 2018, 31(4).

4. Jander G, Rahme LG, Ausubel FM: Positive correlation between virulence of Pseudomonas aeruginosa mutants in mice and insects. J Bacteriol 2000, 182(13):3843-3845.

5. Silo-Suh L, Suh SJ, Sokol PA, Ohman DE: A simple alfalfa seedling infection model for Pseudomonas aeruginosa strains associated with cystic fibrosis shows AlgT (sigma-22) and RhlR contribute to pathogenesis. Proc Natl Acad Sci U S A 2002, 99(24):15699-15704.

6. Al-Aloul M, Crawley J, Winstanley C, Hart CA, Ledson MJ, Walshaw MJ: Increased morbidity associated with chronic infection by an epidemic Pseudomonas aeruginosa strain in CF patients. Thorax 2004, 59(4):334-336.

7. Amiel E, Lovewell RR, O'Toole GA, Hogan DA, Berwin B: Pseudomonas aeruginosa evasion of phagocytosis is mediated by loss of swimming motility and is independent of flagellum expression. Infect Immun 2010, 78(7):2937-2945.

8. De Simone M, Spagnuolo L, Lore NI, Rossi G, Cigana C, De Fino I, Iraqi FA, Bragonzi A: Host genetic background influences the response to the opportunistic Pseudomonas aeruginosa infection altering cell-mediated immunity and bacterial replication. PLoS One 2014, 9(9):e106873.

9. Lagoumintzis G, Christofidou M, Dimitracopoulos G, Paliogianni F: Pseudomonas aeruginosa slime glycolipoprotein is a potent stimulant of tumor necrosis factor alpha gene expression and activation of transcription activators nuclear factor kappa B and activator protein 1 in human monocytes. Infect Immun 2003, 71(8):4614-4622.

10. Mena A, Smith EE, Burns JL, Speert DP, Moskowitz SM, Perez JL, Oliver A: Genetic adaptation of Pseudomonas aeruginosa to the airways of cystic fibrosis patients is catalyzed by hypermutation. J Bacteriol 2008, 190(24):7910-7917.

11. Ebenezer DL, Berdyshev EV, Bronova IA, Liu Y, Tiruppathi C, Komarova Y, Benevolenskaya EV, Suryadevara V, Ha AW, Harijith A et al: Pseudomonas aeruginosa stimulates nuclear sphingosine-1-phosphate generation and epigenetic regulation of lung inflammatory injury. Thorax 2019.

12. Khakimova M, Ahlgren HG, Harrison JJ, English AM, Nguyen D: The stringent response controls catalases in Pseudomonas aeruginosa and is required for hydrogen peroxide and antibiotic tolerance. J Bacteriol 2013, 195(9):2011-2020.

13. Seitz AP, Grassme H, Edwards MJ, Pewzner-Jung Y, Gulbins E: Ceramide and sphingosine in pulmonary infections. Biol Chem 2015, 396(6-7):611-620.

14. Tavakoli Tabazavareh S, Seitz A, Jernigan P, Sehl C, Keitsch S, Lang S, Kahl BC, Edwards M, Grassme H, Gulbins E et al: Lack of Sphingosine Causes Susceptibility to Pulmonary Staphylococcus Aureus Infections in Cystic Fibrosis. Cell Physiol Biochem 2016, 38(6):2094-2102.

15. Grassme H, Becker KA, Zhang Y, Gulbins E: Ceramide in bacterial infections and cystic fibrosis. Biol Chem 2008, 389(11):1371-1379.

16. Peng X, Hassoun PM, Sammani S, McVerry BJ, Burne MJ, Rabb H, Pearse D, Tuder RM, Garcia JG: Protective effects of sphingosine 1-phosphate in murine endotoxin-induced inflammatory lung injury. Am J Respir Crit Care Med 2004, 169(11):1245-1251.

17. Szczepaniak WS, Zhang Y, Hagerty S, Crow MT, Kesari P, Garcia JG, Choi AM, Simon BA, McVerry BJ: Sphingosine 1-phosphate rescues canine LPS-induced acute lung injury and alters systemic inflammatory cytokine production in vivo. Transl Res 2008, 152(5):213-224.

18. Irizarry RA, Bolstad BM, Collin F, Cope LM, Hobbs B, Speed TP: Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res 2003, 31(4):e15.

19. Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics 2003, 19(2):185-193.

20. Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics 2003, 4(2):249-264.

21. Eisenhart C: The assumptions underlying the analysis of variance. Biometrics 1947, 3(1):1-21.

22. Li H: Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv:13033997 [q-bioGN] 2013.

23. Liao Y, Smyth GK, Shi W: featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 2014, 30(7):923-930.

24. McCarthy DJ, Chen Y, Smyth GK: Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res 2012, 40(10):4288-4297.

25. Robinson MD, McCarthy DJ, Smyth GK: edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26(1):139-140.

26. Madduri RK, Sulakhe D, Lacinski L, Liu B, Rodriguez A, Chard K, Dave UJ, Foster IT: Experiences Building Globus Genomics: A Next-Generation Sequencing Analysis Service using Galaxy, Globus, and Amazon Web Services. Concurr Comput 2014, 26(13):2266-2279.

27. Benjamini YH, Y. : Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society Series B (Methodological) 1995, 57(1):289-300.

28. Benjamini YY, D.: The Control of the False Discovey Rate in Multiple Testing under Dependency. The Annals of Statistics 2001, 29(4):1165-1188.

29. Natarajan V, Ha AW, Dong Y, Reddy NM, Ebenezer DL, Kanteti P, Reddy SP, Usha Raj J, Lei Z, Maienschein-Cline M et al: Expression profiling of genes regulated by sphingosine kinase1 signaling in a murine model of hyperoxia induced neonatal bronchopulmonary dysplasia. BMC Genomics 2017, 18(1):664.

30. Romberg N, Virdee M, Chamberlain N, Oe T, Schickel JN, Perkins T, Cantaert T, Rachid R, Rosengren S, Palazzo R et al: TNF receptor superfamily member 13b (TNFRSF13B) hemizygosity reveals transmembrane activator and CAML interactor haploinsufficiency at later stages of B-cell development. J Allergy Clin Immunol 2015, 136(5):1315-1325.

31. Das J, Ramani R, Suraju MO: Polyphenol compounds and PKC signaling. Biochim Biophys Acta 2016, 1860(10):2107-2121.

32. Chase SD, Magnani JL, Simon SI: E-selectin ligands as mechanosensitive receptors on neutrophils in health and disease. Ann Biomed Eng 2012, 40(4):849-859.

33. Li QX, Ke N, Sundaram R, Wong-Staal F: NR4A1, 2, 3--an orphan nuclear hormone receptor family involved in cell apoptosis and carcinogenesis. Histol Histopathol 2006, 21(5):533-540.

34. Turnham RE, Scott JD: Protein kinase A catalytic subunit isoform PRKACA; History, function and physiology. Gene 2016, 577(2):101-108.

35. Wang R, Kwon IK, Thangaraju M, Singh N, Liu K, Jay P, Hofmann F, Ganapathy V, Browning DD: Type 2 cGMP-dependent protein kinase regulates proliferation and differentiation in the colonic mucosa. Am J Physiol Gastrointest Liver Physiol 2012, 303(2):G209-219.

36. Kawasaki Y, Kugimiya F, Chikuda H, Kamekura S, Ikeda T, Kawamura N, Saito T, Shinoda Y, Higashikawa A, Yano F et al: Phosphorylation of GSK-3beta by cGMP-dependent protein kinase II promotes hypertrophic differentiation of murine chondrocytes. J Clin Invest 2008, 118(7):2506-2515.

37. Hyde RK, Liu PP: RUNX1 repression-independent mechanisms of leukemogenesis by fusion genes CBFB-MYH11 and AML1-ETO (RUNX1-RUNX1T1). J Cell Biochem 2010, 110(5):1039-1045.

38. Tay FP, Huang M, Wang L, Yamada Y, Liu DX: Characterization of cellular furin content as a potential factor determining the susceptibility of cultured human and animal cells to coronavirus infectious bronchitis virus infection. Virology 2012, 433(2):421-430.

39. Ichikawa JK, Norris A, Bangera MG, Geiss GK, van 't Wout AB, Bumgarner RE, Lory S: Interaction of pseudomonas aeruginosa with epithelial cells: identification of differentially regulated genes by expression microarray analysis of human cDNAs. Proc Natl Acad Sci U S A 2000, 97(17):9659-9664.

40. Lory S, Ichikawa JK: Pseudomonas-epithelial cell interactions dissected with DNA microarrays. Chest 2002, 121(3 Suppl):36S-39S.

41. Mittal R, Grati M, Yan D, Liu XZ: Pseudomonas aeruginosa Activates PKC-Alpha to Invade Middle Ear Epithelial Cells. Front Microbiol 2016, 7:255.

42. Siggs OM, Miosge LA, Daley SR, Asquith K, Foster PS, Liston A, Goodnow CC: Quantitative reduction of the TCR adapter protein SLP-76 unbalances immunity and immune regulation. J Immunol 2015, 194(6):2587-2595.

43. Li L, Fan D, Wang C, Wang JY, Cui XB, Wu D, Zhou Y, Wu LL: Angiotensin II increases periostin expression via Ras/p38 MAPK/CREB and ERK1/2/TGF-beta1 pathways in cardiac fibroblasts. Cardiovasc Res 2011, 91(1):80-89.

44. Sato T, Takino JI, Nagamine K, Nishio K, Hori T: RASGRP2 Suppresses Apoptosis via Inhibition of ROS Production in Vascular Endothelial Cells. ScientificWorldJournal 2019, 2019:4639165.

45. Mathee K, Ciofu O, Sternberg C, Lindum PW, Campbell JI, Jensen P, Johnsen AH, Givskov M, Ohman DE, Molin S et al: Mucoid conversion of Pseudomonas aeruginosa by hydrogen peroxide: a mechanism for virulence activation in the cystic fibrosis lung. Microbiology 1999, 145 ( Pt 6):1349-1357.

46. Li JD, Dohrman AF, Gallup M, Miyata S, Gum JR, Kim YS, Nadel JA, Prince A, Basbaum CB: Transcriptional activation of mucin by Pseudomonas aeruginosa lipopolysaccharide in the pathogenesis of cystic fibrosis lung disease. Proc Natl Acad Sci U S A 1997, 94(3):967-972.

47. Dortet L, Lombardi C, Cretin F, Dessen A, Filloux A: Pore-forming activity of the Pseudomonas aeruginosa type III secretion system translocon alters the host epigenome. Nat Microbiol 2018, 3(3):378-386.

48. Butterworth MB, Zhang L, Heidrich EM, Myerburg MM, Thibodeau PH: Activation of the epithelial sodium channel (ENaC) by the alkaline protease from Pseudomonas aeruginosa. J Biol Chem 2012, 287(39):32556-32565.

49. Pasetto M, Antignani A, Ormanoglu P, Buehler E, Guha R, Pastan I, Martin SE, FitzGerald DJ: Whole-genome RNAi screen highlights components of the endoplasmic reticulum/Golgi as a source of resistance to immunotoxin-mediated cytotoxicity. Proc Natl Acad Sci U S A 2015, 112(10):E1135-1142.

50. Michalska M, Wolf P: Pseudomonas Exotoxin A: optimized by evolution for effective killing. Front Microbiol 2015, 6:963.

51. Inocencio NM, Moehring JM, Moehring TJ: Furin activates Pseudomonas exotoxin A by specific cleavage in vivo and in vitro. J Biol Chem 1994, 269(50):31831-31835.

Table 1.

Name of gene	Function of the corresponding protein
Spata13 Spermatogenesis Associated 13	Serves as guanine nucleotide exchange factor (GEF) for RHOA, RAC1 and CDC42 GTPases. Involved in regulation of cell migration and adhesion assembly and disassembly.
Hdc Histidine Decarboxylase	Converts L-histidine to histamine
Kdm6b Lysine Demethylase 6B	Lysine-specific demethylase which specifically demethylates di- or tri-methylated lysine 27 of histone H3
G0s2 G0/G1 Switch 2	Promotes apoptosis by binding to BCL2, thus preventing the formation of anti-apoptotic BCL2-BAX heterodimers
Fosl2 FOS Like 2, AP-1 Transcription Factor Subunit	The gene encodes leucine zipper proteins which can dimerize with proteins of the JUN family. resulting in formation of the transcription factor complex AP-1. FOS proteins act as regulators of cell proliferation, differentiation, and transformation.
Dgat1 Diacylglycerol O-Acyltransferase 1	Catalyzes conversion of diacylglycerol and fatty acyl CoA to triacylglycerol
Odc1 Ornithine Decarboxylase 1	Catalyzes the first and rate-limiting step of polyamine biosynthesis converting ornithine into putrescine. Putrescine is the precursor for the polyamines which are essential for cellular processes, ranging from DNA replication to apoptosis.
Smtnl2 Smoothelin Like 2	Smoothelin is a structural protein found exclusively in contractile smooth muscle cells
Clec1a C-Type Lectin Domain Family 1 Member A	This protein serves diverse functions, such as cell adhesion, cell-cell signaling, inflammation and immune response
Sox17 Sex Determining Region Y Box 17	Is a transcription regulator binding to target promoter DNA and bends the DNA.

Download PDF

Journal Publication

published 16 Dec, 2019

Read the published version in BMC Genomics →

Editorial decision: Major revision
02 Jul, 2019
Review #1 received at journal
01 Jul, 2019
Review #2 received at journal
01 Jul, 2019
Submission checks completed at journal
13 Jun, 2019
Reviewers invited by journal
13 Jun, 2019
Reviewer #1 agreed at journal
13 Jun, 2019
Reviewer #2 agreed at journal
13 Jun, 2019
Editor invited by journal
28 May, 2019
Editor assigned by journal
28 May, 2019
First submitted to journal
22 May, 2019

You are reading this older preprint version

Read the latest preprint version →

Genetic deletion of Sphk2 confers protection against Pseudomonas aeruginosa mediated differential expression of genes related to virulent infection and inflammation in mouse lung

Status:

Journal Publication

Version 1

Abstract

Figures

BACKGROUND

METHODS

RESULTS

DISCUSSION

CONCLUSION

ABBREVIATIONS

DECLARATIONS

REFERENCES

Table 1

Supplementary Files

Status:

Journal Publication

Version 1