Bacterial adaptation during chronic infection revealed by independent component analysis of transcriptomic data
© Yang et al; licensee BioMed Central Ltd. 2011
Received: 5 June 2011
Accepted: 18 August 2011
Published: 18 August 2011
Bacteria employ a variety of adaptation strategies during the course of chronic infections. Understanding bacterial adaptation can facilitate the identification of novel drug targets for better treatment of infectious diseases. Transcriptome profiling is a comprehensive and high-throughput approach for characterization of bacterial clinical isolates from infections. However, exploitation of the complex, noisy and high-dimensional transcriptomic dataset is difficult and often hindered by low statistical power.
In this study, we have applied two kinds of unsupervised analysis methods, principle component analysis (PCA) and independent component analysis (ICA), to extract and characterize the most informative features from transcriptomic dataset generated from cystic fibrosis (CF) Pseudomonas aeruginosa isolates. ICA was shown to be able to efficiently extract biological meaningful features from the transcriptomic dataset and improve clustering patterns of CF isolates. Decomposition of the transcriptomic dataset by ICA also facilitates gene identification and gene ontology enrichment.
Our results show that P. aeruginosa employs multiple patient-specific adaption strategies during the early stage infections while certain essential adaptations are evolved in parallel during the chronic infections.
Bacterial infections are one of the major causes of mortality among human and animals in the world . Understanding adaptation of bacterial pathogens to the dynamic and hostile environment is crucial for improvement of therapies of infectious diseases. Bacteria associated with chronic infections in patients suffering from e.g. AIDS, burn wound sepsis, diabetes and cystic fibrosis (CF) are ideal objects for studying bacterial adaptation.
In airways of CF patients, mucus forms a stationary and thickened gel adhering to the epithelial lining fluid of the airway surfaces, which affects the mucociliary escalator and results in impaired clearance of inhaled microbes . CF patients suffer from chronic and recurrent respiratory tract infections which eventually lead to lung failure followed by death. Pseudomonas aeruginosa is one of the major pathogens for CF patients and is the principal cause of mortality and morbidity in CF patients . Early P. aeruginosa infection in CF patients is characterized by a diverse of P. aeruginosa strains which have similar phenotypes as those of environmental isolates [4, 5]. In contrast, adapted dominant epidemic strains are often identified from patients chronically infected with P. aeruginosa from different CF centers [4, 6, 7]. Once it gets adapted, P. aeruginosa can persist for several decades in the respiratory tracts of CF patients, overcoming host defense mechanisms as well as intensive antibiotic therapies .
As P. aeruginosa has been sequenced, transcriptome profiling (e.g. microarray analysis and RNA-Seq) becomes a convenient approach for characterizing biological differences among different P. aeruginosa clinical isolates from CF patients. Transcriptome profiling enables researchers to measure genome-wide gene expressions in a high-throughput manner thus can provide valuable information for P. aeruginosa adaptation during infections. However, the interpretation of transcriptomic data is a great challenge for researchers due to the complexity and noise. Clinical strains isolated from different patients have adapted to distinct host environments since patients vary in their ages, infection histories and medical treatments (e.g. different kinds of antibiotics and their dosages). Therefore, researchers need to reduce dimensionality and extract the underlying features from the multi-variable transcriptomic dataset.
Principle component analysis (PCA) is a classic projection method which is widely used to accomplish the above mentioned tasks . PCA transforms a number of correlated variables into a smaller number of uncorrelated variables called principal components (PC). The first PC captures as much of the variability in the data as possible, and each succeeding PCs capture as much of the remaining variability as possible. However, the constraint of mutual orthogonality of components implied in classical PCA methods may not be appropriate for the biological systems. Recently, independent component analysis (ICA), which decomposes input data into statistically independent components, was shown to be able to classify gene expressions into biologically meaningful groups and relate them to specific biological processes . ICA has been successfully applied by different research groups to analyze transcriptomic data from yeast, cancer, Alzheimer samples and is shown to be more powerful at feature extraction than PCA and other traditional methods for microarray data analysis [11–13]. In a study by Zhang et al., ICA was used to extract specific gene expression patterns of normal and tumor tissues, which can serve as biomarkers for molecular diagnosis of human cancer type . Yet to the best of our knowledge, there have been no reports of application of ICA to the study of bacterial transcriptomic data from chronic infections.
In this study, we applied ICA to project the transcriptomic data of 26 CF P. aeruginosa isolates into independent components. P. aeruginosa genes are unsupervisedly clustered into non-mutually exclusive groups. Each retrieved independent component is considered as a putative adaptation process, which is revealed by the functional annotations of genes that give heavy loadings to the component.
ICA improved clustering patterns of P. aeruginosamicroarray data
ICA identified significant genes for adaptation of P. aeruginosato the CF airways
Latent variables related to specific adaptation
Functions of selected enriched genes by ICA
Early stage isolates from 1973
Type III secretion
Late stage isolates
Antimicrobial peptide tolerance
Late stage isolates
Potassium uptake system
Late stage isolates
ICA revealed common adaptations shared by a group of P. aeruginosa CF isolates. IC14 revealed that the early stage isolates from 1973 had higher expression level of genes involved in type III secretion and exoenzyme activities than other isolates (Figure 4 and Additional file 1, Table S1). More importantly, IC6, IC10 and IC18 revealed adaptations shared by the late stage isolates. IC6 mainly identified antimicrobial peptide resistance related arn and pmr genes (PA3552-PA3559 and PA4773-PA4782) (Figure 4 and Additional file 1, Table S1). IC10 mainly identified the alginate biosynthesis regulatory algU (PA0762), mucA (PA0763), mucB (PA0764), mucC (PA0765) and algR (PA5261) genes; the potassium uptake kdp genes (PA1632-PA1635) (Figure 4 and Additional file 1, Table S1) and the quorum sensing genes (PA1430-PA1431) (Figure 4 and Additional file 1, Table S1). IC18 mainly identified alginate biosynthesis alg genes (PA3540-PA3551) and flagellum and type IV pilus biogenesis genes (Figure 4 and Additional file 1, Table S1).
Besides common adaptations shared by a group of P. aeruginosa CF isolates, the ICA also showed that P. aeruginosa CF isolates from early infection stage employed multiple patient-specific strategies of adaptation in the CF airways. IC2 revealed that the early stage B12-4 and B12-7 isolates induced the expression of genes related to MexAB-OprM efflux system, iron uptake as well as citronellol/leucine catabolism (Figure 4 and Additional file 1, Table S1). IC4 revealed that the early stage B6-0 and B6-4 isolates up-regulated expression of LPS biosynthesis wbp genes (PA3146-PA3159) and down-regulated expression of genes involved in the flagellum biogenesis (Figure 4 and Additional file 1, Table S1). IC16 revealed that the early stage CF114-1973 isolate up-regulated the expression of genes involved in fimbrial biogenesis while down-regulated expression of the PA0632-PA0639 genes (Figure 4 and Additional file 1, Table S1). IC20 revealed that the late stage CF66-2008 isolate up-regulated the expression of the LPS biosynthesis wbp genes (PA5448-PA5454) (Figure 4 and Additional file 1, Table S1).
ICA enhanced identification of co-regulated genes for adaptation of P. aeruginosato the CF airways
Understanding the bacterial adaptation is a great challenge for scientists and medical doctors to battle infectious diseases. Bacterial cells have a high level of mutation rate and can adapt to the dynamic host environments by selecting mutants which are more fit to the condition. Thus, a systematic investigation of the whole gene expression profiles of clinical isolates would be needed for modern diagnostic and treatment of infectious diseases. Fortunately, the rapid progress of DNA sequencing projects has made genome sequences of most of the pathogenic bacteria available now. And this has brought DNA microarray technique as a conventional and high-throughput tool for researchers. However, how to properly and accurately analyze the microarray data and extract useful information is another obstacle for using DNA microarray technique.
In the study here, we have analyzed DNA microarray dataset generated from 26 P. aeruginosa strains. ICA was shown to be an efficient approach to identify patient-specific adaptations of P. aeruginosa isolates. First of all, ICA decomposes and extracts genes from the microarray dataset simultaneously. Thus, co-regulated genes are more easily identified (Figure 6). Secondly, unlike conventional clustering approaches which group genes based on their expression levels, ICA grouped genes independent of expression levels but in a more biologically meaningful manner.
ICA shows that P. aeruginosa clinical isolates employ multiple patient-specific adaption strategies during the early stage infection. Most of these early stage adaptive changes are involved in modification of cell surface molecules and appendages. IC4 reveals that B6-0 and B6-4 isolates enhanced the expression of B-band lipopolysaccharide (LPS) biosynthesis genes while reduced the expression of flagellum biogenesis genes. The B-band LPS is a well known virulence factor which confers P. aeruginosa resistance to phagocytosis and serum-mediated killing [17–20]. Loss of flagellum as well as flagellum-mediated motility is documented to render P. aeruginosa CF isolates an advantage in the context of immune evasion [21–23]. IC16 reveals that CF114-1973 isolate enhanced the expression of the cupA fimbrial gene cluster and the type IV pilus biogenesis cluster. The gene products of these two clusters are required for P. aeruginosa adherence and biofilm formation [24–28]. Interestingly, IC16 also reveals the increased expression of pprB gene in CF114-1973, which was recently reported as a new regulatory element controlling the cupE gene expression and transition between planktonic and community lifestyles in P. aeruginosa .
ICA facilitates enrichment of co-regulated genes of P. aeruginosa CF isolates. For example, IC6 groups the two antimicrobial peptide resistance related gene clusters (arn and pmr) together. IC18 groups alginate biosynthesis gene cluster PA3540-PA3551 and flagellum biogenesis gene cluster PA1077-PA1086 together. These two gene clusters are impossible to be grouped together by other approaches since they are not localized adjacently in the genome and have different expression levels (one up-regulated and one down-regulated). And this grouping is biologically meaningful since it is well known that alginate regulator inhibits flagellum synthesis gene expression [30–32]. Many genes which encode hypothetical proteins are grouped in IC6, IC10 and IC18. It will be interesting to investigate whether these genes are functionally related with the annotated genes identified in the same ICs.
Since ICA can reveal patient-specific adaptations of P. aeruginsoa isolates, it is possible to design patient-specific therapies based on these adaptations. For example, combination of iron chelators and efflux pump inhibitors might be used to inhibit the growth of B12-4 and B12-7, which have high expression levels of genes involved in efflux pump and iron uptake systems . Ligands with high affinity to pili can be used to inhibit adhesion and biofilm formation of the CF114-1973 isolate .
In conclusion, the ICA is shown to be able to extract the most essential features from the complex multiple variant microarray dataset and identify significant genes contribute to these features. Our results show that P. aeruginosa employ a diverse set of patient-specific adaption strategies during the early stage infections while certain essential evolutionary events occurred in parallel during the chronic infections in CF infections. The ICA has a great potential in studying large-scale datasets acquired from omics research from different areas.
P. aeruginosa clinical isolates
The P. aeruginosa strains were isolated from 6 CF patients with long-term chronic infection and 3 CF patients who were intermittently colonized or recently chronically infected and who were attending the Danish CF Center, Rigshospitalet, Copenhagen. P. aeruginosa PAO1  was used as a reference strain.
Transcriptomic profiles of clinical isolates were obtained using the Affymetrix P. aeruginosa gene chip (Santa Clara, CA) [5, 8]. Triplicate experiments were performed for each strain. The microarray raw datasets are accessible at NCBI's Gene Expression Omnibus (GEO) with series accession number GSE31227.
Mathematical model of gene regulation by ICA
In this ICA model, the columns of A= [a 1 , a 2 ,..., a n ] are the n × n latent vectors of the gene microarray data. Each column of Ais associated with a specific gene expression mode. Scontains the n × m gene signatures where the rows of Sare statistically independent to each other. The gene profiles in Xare considered to be a linear mixture of statistically independent components Scombined by an unknown mixing matrix A. Once latent variable matrix Ahas been obtained, the corresponding elementary modes can be identified to extract information for classification.
This work was supported by a grant from the Danish Research Council for Independent Research (09-073917) to L.Y.
- Demuth A, Aharonowitz Y, Bachmann TT, Blum-Oehler G, Buchrieser C, Covacci A, Dobrindt U, Emody L, van der Ende A, Ewbank J, et al: Pathogenomics: an updated European Research Agenda. Infect Genet Evol. 2008, 8 (3): 386-393. 10.1016/j.meegid.2008.01.005.PubMedView ArticleGoogle Scholar
- Worlitzsch D, Tarran R, Ulrich M, Schwab U, Cekici A, Meyer KC, Birrer P, Bellon G, Berger J, Weiss T, et al: Effects of reduced mucus oxygen concentration in airway Pseudomonas infections of cystic fibrosis patients. J Clin Invest. 2002, 109 (3): 317-325.PubMedPubMed CentralView ArticleGoogle Scholar
- Govan JR, Deretic V: Microbial pathogenesis in cystic fibrosis: mucoid Pseudomonas aeruginosa and Burkholderia cepacia. Microbiol Rev. 1996, 60 (3): 539-574.PubMedPubMed CentralGoogle Scholar
- Jelsbak L, Johansen HK, Frost AL, Thogersen R, Thomsen LE, Ciofu O, Yang L, Haagensen JA, Hoiby N, Molin S: Molecular epidemiology and dynamics of Pseudomonas aeruginosa populations in lungs of cystic fibrosis patients. Infect Immun. 2007, 75 (5): 2214-2224. 10.1128/IAI.01282-06.PubMedPubMed CentralView ArticleGoogle Scholar
- Rau MH, Hansen SK, Johansen HK, Thomsen LE, Workman CT, Nielsen KF, Jelsbak L, Hoiby N, Yang L, Molin S: Early adaptive developments of Pseudomonas aeruginosa after the transition from life in the environment to persistent colonization in the airways of human cystic fibrosis hosts. Environ Microbiol. 2010, 12 (6): 1643-1658.PubMedGoogle Scholar
- Romling U, Fiedler B, Bosshammer J, Grothues D, Greipel J, von der Hardt H, Tummler B: Epidemiology of chronic Pseudomonas aeruginosa infections in cystic fibrosis. J Infect Dis. 1994, 170 (6): 1616-1621. 10.1093/infdis/170.6.1616.PubMedView ArticleGoogle Scholar
- Smith EE, Buckley DG, Wu Z, Saenphimmachak C, Hoffman LR, D'Argenio DA, Miller SI, Ramsey BW, Speert DP, Moskowitz SM, et al: Genetic adaptation by Pseudomonas aeruginosa to the airways of cystic fibrosis patients. Proc Natl Acad Sci USA. 2006, 103 (22): 8487-8492. 10.1073/pnas.0602138103.PubMedPubMed CentralView ArticleGoogle Scholar
- Yang L, Jelsbak L, Marvig RL, Damkiaer S, Workman CT, Rau MH, Hansen SK, Folkesson A, Johansen HK, Ciofu O, et al: Evolutionary dynamics of bacteria in a human host environment. Proc Natl Acad Sci USA. 2011, 108 (18): 7481-7486. 10.1073/pnas.1018249108.PubMedPubMed CentralView ArticleGoogle Scholar
- Raychaudhuri S, Stuart JM, Altman RB: Principal components analysis to summarize microarray experiments: application to sporulation time series. Pac Symp Biocomput. 2000, 455-466.Google Scholar
- Kong W, Vanderburg CR, Gunshin H, Rogers JT, Huang X: A review of independent component analysis application to microarray gene expression data. Biotechniques. 2008, 45 (5): 501-520. 10.2144/000112950.PubMedPubMed CentralView ArticleGoogle Scholar
- Lee SI, Batzoglou S: Application of independent component analysis to microarrays. Genome Biol. 2003, 4 (11): R76-10.1186/gb-2003-4-11-r76.PubMedPubMed CentralView ArticleGoogle Scholar
- Saidi SA, Holland CM, Kreil DP, MacKay DJ, Charnock-Jones DS, Print CG, Smith SK: Independent component analysis of microarray data in the study of endometrial cancer. Oncogene. 2004, 23 (39): 6677-6683. 10.1038/sj.onc.1207562.PubMedView ArticleGoogle Scholar
- Kong W, Mou X, Liu Q, Chen Z, Vanderburg CR, Rogers JT, Huang X: Independent component analysis of Alzheimer's DNA microarray gene expression data. Mol Neurodegener. 2009, 4 (1): 5-10.1186/1750-1326-4-5.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang XW, Yap YL, Wei D, Chen F, Danchin A: Molecular diagnosis of human cancer type by gene expression profiles and independent component analysis. Eur J Hum Genet. 2005, 13 (12): 1303-1311. 10.1038/sj.ejhg.5201495.PubMedView ArticleGoogle Scholar
- Hyvarinen A, Oja E: Independent component analysis: algorithms and applications. Neural Netw. 2000, 13 (4-5): 411-430. 10.1016/S0893-6080(00)00026-5.PubMedView ArticleGoogle Scholar
- Smyth GK: limma: Linear Models for Microarray Data. Edited by: Gentleman R, Carey VJ, Huber W, Irizarry RA, Dudoit S. 2005, Bioinformatics and Computational Biology Solutions using R and Bioconductor NY: SpringerGoogle Scholar
- Dasgupta T, de Kievit TR, Masoud H, Altman E, Richards JC, Sadovskaya I, Speert DP, Lam JS: Characterization of lipopolysaccharide-deficient mutants of Pseudomonas aeruginosa derived from serotypes O3, O5, and O6. Infect Immun. 1994, 62 (3): 809-817.PubMedPubMed CentralGoogle Scholar
- Cryz SJ, Pitt TL, Furer E, Germanier R: Role of lipopolysaccharide in virulence of Pseudomonas aeruginosa. Infect Immun. 1984, 44 (2): 508-513.PubMedPubMed CentralGoogle Scholar
- Engels W, Endert J, Kamps MA, van Boven CP: Role of lipopolysaccharide in opsonization and phagocytosis of Pseudomonas aeruginosa. Infect Immun. 1985, 49 (1): 182-189.PubMedPubMed CentralGoogle Scholar
- Hancock RE, Mutharia LM, Chan L, Darveau RP, Speert DP, Pier GB: Pseudomonas aeruginosa isolates from patients with cystic fibrosis: a class of serum-sensitive, nontypable strains deficient in lipopolysaccharide O side chains. Infect Immun. 1983, 42 (1): 170-177.PubMedPubMed CentralGoogle Scholar
- 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. 10.1128/IAI.00144-10.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang Z, Louboutin JP, Weiner DJ, Goldberg JB, Wilson JM: Human airway epithelial cells sense Pseudomonas aeruginosa infection via recognition of flagellin by Toll-like receptor 5. Infect Immun. 2005, 73 (11): 7151-7160. 10.1128/IAI.73.11.7151-7160.2005.PubMedPubMed CentralView ArticleGoogle Scholar
- Mahenthiralingam E, Speert DP: Nonopsonic phagocytosis of Pseudomonas aeruginosa by macrophages and polymorphonuclear leukocytes requires the presence of the bacterial flagellum. Infect Immun. 1995, 63 (11): 4519-4523.PubMedPubMed CentralGoogle Scholar
- Vallet I, Olson JW, Lory S, Lazdunski A, Filloux A: The chaperone/usher pathways of Pseudomonas aeruginosa: identification of fimbrial gene clusters (cup) and their involvement in biofilm formation. Proc Natl Acad Sci USA. 2001, 98 (12): 6911-6916. 10.1073/pnas.111551898.PubMedPubMed CentralView ArticleGoogle Scholar
- O'Toole GA, Kolter R: Flagellar and twitching motility are necessary for Pseudomonas aeruginosa biofilm development. Mol Microbiol. 1998, 30 (2): 295-304. 10.1046/j.1365-2958.1998.01062.x.PubMedView ArticleGoogle Scholar
- Klausen M, Heydorn A, Ragas P, Lambertsen L, Aaes-Jorgensen A, Molin S, Tolker-Nielsen T: Biofilm formation by Pseudomonas aeruginosa wild type, flagella and type IV pili mutants. Mol Microbiol. 2003, 48 (6): 1511-1524. 10.1046/j.1365-2958.2003.03525.x.PubMedView ArticleGoogle Scholar
- Barken KB, Pamp SJ, Yang L, Gjermansen M, Bertrand JJ, Klausen M, Givskov M, Whitchurch CB, Engel JN, Tolker-Nielsen T: Roles of type IV pili, flagellum-mediated motility and extracellular DNA in the formation of mature multicellular structures in Pseudomonas aeruginosa biofilms. Environ Microbiol. 2008, 10 (9): 2331-2343. 10.1111/j.1462-2920.2008.01658.x.PubMedView ArticleGoogle Scholar
- Ruer S, Stender S, Filloux A, de Bentzmann S: Assembly of fimbrial structures in Pseudomonas aeruginosa: functionality and specificity of chaperone-usher machineries. J Bacteriol. 2007, 189 (9): 3547-3555. 10.1128/JB.00093-07.PubMedPubMed CentralView ArticleGoogle Scholar
- Giraud C, Bernard CS, Calderon V, Yang L, Filloux A, Molin S, Fichant G, Bordi C, de Bentzmann S: The PprA-PprB two-component system activates CupE, the first non-archetypal Pseudomonas aeruginosa chaperone-usher pathway system assembling fimbriae. Environ Microbiol. 2011, 13 (3): 666-683. 10.1111/j.1462-2920.2010.02372.x.PubMedView ArticleGoogle Scholar
- Garrett ES, Perlegas D, Wozniak DJ: Negative control of flagellum synthesis in Pseudomonas aeruginosa is modulated by the alternative sigma factor AlgT (AlgU). J Bacteriol. 1999, 181 (23): 7401-7404.PubMedPubMed CentralGoogle Scholar
- Tart AH, Wolfgang MC, Wozniak DJ: The alternative sigma factor AlgT represses Pseudomonas aeruginosa flagellum biosynthesis by inhibiting expression of fleQ. J Bacteriol. 2005, 187 (23): 7955-7962. 10.1128/JB.187.23.7955-7962.2005.PubMedPubMed CentralView ArticleGoogle Scholar
- Tart AH, Blanks MJ, Wozniak DJ: The AlgT-dependent transcriptional regulator AmrZ (AlgZ) inhibits flagellum biosynthesis in mucoid, nonmotile Pseudomonas aeruginosa cystic fibrosis isolates. J Bacteriol. 2006, 188 (18): 6483-6489. 10.1128/JB.00636-06.PubMedPubMed CentralView ArticleGoogle Scholar
- Liu Y, Yang L, Molin S: Synergistic activities of an efflux pump inhibitor and iron chelators against Pseudomonas aeruginosa growth and biofilm formation. Antimicrob Agents Chemother. 2010, 54 (9): 3960-3963. 10.1128/AAC.00463-10.PubMedPubMed CentralView ArticleGoogle Scholar
- Wu HY, Zhang XL, Pan Q, Wu J: Functional selection of a type IV pili-binding peptide that specifically inhibits Salmonella Typhi adhesion to/invasion of human monocytic cells. Peptides. 2005, 26 (11): 2057-2063. 10.1016/j.peptides.2005.03.035.PubMedView ArticleGoogle Scholar
- Holloway BW, Morgan AF: Genome organization in Pseudomonas. Annu Rev Microbiol. 1986, 40: 79-105. 10.1146/annurev.mi.40.100186.000455.PubMedView ArticleGoogle Scholar