Several bacterial species from the Burkholderia cepacia complex (Bcc) are feared opportunistic pathogens that lead to debilitating lung infections with a high risk of developing fatal septicemia in cystic fibrosis (CF) patients. However, the pathogenic potential of other Bcc species is yet unknown. To elucidate clinical relevance of Burkholderia contaminans, a species frequently isolated from CF respiratory samples in Ibero-American countries, we aimed to identify its key virulence factors possibly linked with an unfavorable clinical outcome. We performed a genome-wide comparative analysis of two isolates of B. contaminans ST872 from sputum and blood culture of a female CF patient in Argentina. RNA-seq data showed significant changes in expression for quorum sensing-regulated virulence factors and motility and chemotaxis. Furthermore, we detected expression changes in a recently described low-oxygen-activated (lxa) locus which encodes stress-related proteins, and for two clusters responsible for the biosynthesis of antifungal and hemolytic compounds pyrrolnitrin and occidiofungin. Based on phenotypic assays that confirmed changes in motility and in proteolytic, hemolytic and antifungal activities, we were able to distinguish two phenotypes of B. contaminans that coexisted in the host and entered her bloodstream. Whole genome sequencing revealed that the sputum and bloodstream isolates (each representing a distinct phenotype) differed by over 1,400 mutations as a result of a mismatch repair-deficient hypermutable state of the sputum isolate. The inferred lack of purifying selection against nonsynonymous mutations and the high rate of pseudogenization in the derived isolate indicated limited evolutionary pressure during evolution in the nutrient-rich, stable CF sputum environment. The present study is the first to examine the genomic and transcriptomic differences between longitudinal isolates of B. contaminans. Detected activity of a number of putative virulence factors implies a genuine pathogenic nature of this novel Bcc species.

Centro Respiratorio Dr A Alvarez Hospital de Niños Ricardo Gutiérrez Buenos Aires Argentina

Department of Medical Microbiology 2nd Faculty of Medicine Charles University and Motol University Hospital Prague Czech Republic

Department of Medical Microbiology and Infectious Disease University of Manitoba Winnipeg Canada

Laboratory of Genomics and Bioinformatics Institute of Molecular Genetics Academy of Sciences of the Czech Republic Prague Czech Republic

School of Pharmacy and Biochemistry University of Buenos Aires Buenos Aires Argentina

Zobrazit více v PubMed

Mahenthiralingam E, Urban TA, Goldberg JB. The multifarious, multireplicon Burkholderia cepacia complex. Nat Rev Microbiol. 2005;3(2):144–56. 10.1038/nrmicro1085 . PubMed DOI

Drevinek P, Mahenthiralingam E. Burkholderia cenocepacia in cystic fibrosis: epidemiology and molecular mechanisms of virulence. Clin Microbiol Infect. 2010;16(7):821–30. 10.1111/j.1469-0691.2010.03237.x . PubMed DOI

Lipuma JJ. The changing microbial epidemiology in cystic fibrosis. Clin Microbiol Rev. 2010;23(2):299–323. 10.1128/CMR.00068-09 PubMed DOI PMC

Kalish LA, Waltz DA, Dovey M, Potter-Bynoe G, McAdam AJ, Lipuma JJ, et al. Impact of Burkholderia dolosa on lung function and survival in cystic fibrosis. Am J Respir Crit Care Med. 2006;173(4):421–5. 10.1164/rccm.200503-344OC . PubMed DOI

Blackburn L, Brownlee K, Conway S, Denton M. 'Cepacia syndrome' with Burkholderia multivorans, 9 years after initial colonization. J Cyst Fibros. 2004;3(2):133–4. 10.1016/j.jcf.2004.03.007 . PubMed DOI

Jones AM, Dodd ME, Govan JR, Barcus V, Doherty CJ, Morris J, et al. Burkholderia cenocepacia and Burkholderia multivorans: influence on survival in cystic fibrosis. Thorax. 2004;59(11):948–51. 10.1136/thx.2003.017210 PubMed DOI PMC

Vanlaere E, Baldwin A, Gevers D, Henry D, De Brandt E, LiPuma JJ, et al. Taxon K, a complex within the Burkholderia cepacia complex, comprises at least two novel species, Burkholderia contaminans sp. nov. and Burkholderia lata sp. nov. Int J Syst Evol Microbiol. 2009;59(Pt 1):102–11. 10.1099/ijs.0.001123-0 . PubMed DOI

Cunha MV, Leitão JH, Mahenthiralingam E, Vandamme P, Lito L, Barreto C, et al. Molecular analysis of Burkholderia cepacia complex isolates from a Portuguese cystic fibrosis center: a 7-year study. J Clin Microbiol. 2003;41(9):4113–20. 10.1128/JCM.41.9.4113-4120.2003 PubMed DOI PMC

Jordá-Vargas L, Degrossi J, Castañeda NC, D'Aquino M, Valvano MA, Procopio A, et al. Prevalence of indeterminate genetic species of Burkholderia cepacia complex in a cystic fibrosis center in Argentina. J Clin Microbiol. 2008;46(3):1151–2. 10.1128/JCM.01595-07 PubMed DOI PMC

Mahenthiralingam E, Baldwin A, Drevinek P, Vanlaere E, Vandamme P, LiPuma JJ, et al. Multilocus sequence typing breathes life into a microbial metagenome. PLoS One. 2006;1:e17 10.1371/journal.pone.0000017 PubMed DOI PMC

Medina-Pascual MJ, Valdezate S, Carrasco G, Villalón P, Garrido N, Saéz-Nieto JA. Increase in isolation of Burkholderia contaminans from Spanish patients with cystic fibrosis. Clin Microbiol Infect. 2015;21(2):150–6. 10.1016/j.cmi.2014.07.014 . PubMed DOI

Martina P, Bettiol M, Vescina C, Montanaro P, Mannino MC, Prieto CI, et al. Genetic diversity of Burkholderia contaminans isolates from cystic fibrosis patients in Argentina. J Clin Microbiol. 2013;51(1):339–44. 10.1128/JCM.02500-12 PubMed DOI PMC

Coutinho CP, Barreto C, Pereira L, Lito L, Melo Cristino J, Sá-Correia I. Incidence of Burkholderia contaminans at a cystic fibrosis center with an unusually high representation of Burkholderia cepacia during 15 years of epidemiological surveillance. J Med Microbiol. 2015;64:927–35. 10.1099/jmm.0.000094 . PubMed DOI

Kalferstova L, Kolar M, Fila L, Vavrova J, Drevinek P. Gene expression profiling of Burkholderia cenocepacia at the time of cepacia syndrome: loss of motility as a marker of poor prognosis? J Clin Microbiol. 2015;53(5):1515–22. 10.1128/JCM.03605-14 PubMed DOI PMC

Bloodworth RA, Selin C, López De Volder MA, Drevinek P, Galanternik L, Degrossi J, et al. Draft genome sequences of Burkholderia contaminans, a Burkholderia cepacia complex species that is increasingly recovered from cystic fibrosis patients. Genome Announc. 2015;3(4). 10.1128/genomeA.00766-15 . PubMed DOI PMC

Subramoni S, Sokol PA. Quorum sensing systems influence Burkholderia cenocepacia virulence. Future Microbiol. 2012;7(12):1373–87. 10.2217/fmb.12.118 . PubMed DOI

Gingues S, Kooi C, Visser MB, Subsin B, Sokol PA. Distribution and expression of the ZmpA metalloprotease in the Burkholderia cepacia complex. J Bacteriol. 2005;187(24):8247–55. 10.1128/JB.187.24.8247-8255.2005. PubMed DOI PMC

Kooi C, Subsin B, Chen R, Pohorelic B, Sokol PA. Burkholderia cenocepacia ZmpB is a broad-specificity zinc metalloprotease involved in virulence. Infect Immun. 2006;74(7):4083–93. 10.1128/IAI.00297-06 PubMed DOI PMC

Sokol PA, Darling P, Woods DE, Mahenthiralingam E, Kooi C. Role of ornibactin biosynthesis in the virulence of Burkholderia cepacia: characterization of pvdA, the gene encoding L-ornithine N(5)-oxygenase. Infect Immun. 1999;67(9):4443–55. PubMed PMC

Lameignere E, Malinovská L, Sláviková M, Duchaud E, Mitchell EP, Varrot A, et al. Structural basis for mannose recognition by a lectin from opportunistic bacteria Burkholderia cenocepacia. Biochem J. 2008;411(2):307–18. 10.1042/bj20071276 . PubMed DOI

Lameignere E, Shiao TC, Roy R, Wimmerova M, Dubreuil F, Varrot A, et al. Structural basis of the affinity for oligomannosides and analogs displayed by BC2L-A, a Burkholderia cenocepacia soluble lectin. Glycobiology. 2010;20(1):87–98. 10.1093/glycob/cwp151 . PubMed DOI

Sulák O, Cioci G, Lameignère E, Balloy V, Round A, Gutsche I, et al. Burkholderia cenocepacia BC2L-C is a super lectin with dual specificity and proinflammatory activity. PLoS Pathog. 2011;7(9):e1002238 10.1371/journal.ppat.1002238 PubMed DOI PMC

Huber B, Feldmann F, Köthe M, Vandamme P, Wopperer J, Riedel K, et al. Identification of a novel virulence factor in Burkholderia cenocepacia H111 required for efficient slow killing of Caenorhabditis elegans. Infect Immun. 2004;72(12):7220–30. 10.1128/IAI.72.12.7220-7230.2004 PubMed DOI PMC

Holden MT, Seth-Smith HM, Crossman LC, Sebaihia M, Bentley SD, Cerdeño-Tárraga AM, et al. The genome of Burkholderia cenocepacia J2315, an epidemic pathogen of cystic fibrosis patients. J Bacteriol. 2009;191(1):261–77. 10.1128/JB.01230-08 PubMed DOI PMC

O'Grady EP, Viteri DF, Malott RJ, Sokol PA. Reciprocal regulation by the CepIR and CciIR quorum sensing systems in Burkholderia cenocepacia. BMC Genomics. 2009;10:441 10.1186/1471-2164-10-441 PubMed DOI PMC

McKeon SA, Nguyen DT, Viteri DF, Zlosnik JE, Sokol PA. Functional quorum sensing systems are maintained during chronic Burkholderia cepacia complex infections in patients with cystic fibrosis. J Infect Dis. 2011;203(3):383–92. 10.1093/infdis/jiq054 PubMed DOI PMC

Mira NP, Madeira A, Moreira AS, Coutinho CP, Sá-Correia I. Genomic expression analysis reveals strategies of Burkholderia cenocepacia to adapt to cystic fibrosis patients' airways and antimicrobial therapy. PLoS One. 2011;6(12):e28831 10.1371/journal.pone.0028831 PubMed DOI PMC

Zlosnik JE, Mori PY, To D, Leung J, Hird TJ, Speert DP. Swimming motility in a longitudinal collection of clinical isolates of Burkholderia cepacia complex bacteria from people with cystic fibrosis. PLoS One. 2014;9(9):e106428 10.1371/journal.pone.0106428 PubMed DOI PMC

el-Banna N, Winkelmann G. Pyrrolnitrin from Burkholderia cepacia: antibiotic activity against fungi and novel activities against streptomycetes. J Appl Microbiol. 1998;85(1):69–78. 10.1046/j.1365-2672.1998.00473.x . PubMed DOI

Schmidt S, Blom JF, Pernthaler J, Berg G, Baldwin A, Mahenthiralingam E, et al. Production of the antifungal compound pyrrolnitrin is quorum sensing-regulated in members of the Burkholderia cepacia complex. Environ Microbiol. 2009;11(6):1422–37. 10.1111/j.1462-2920.2009.01870.x . PubMed DOI

Gu G, Smith L, Wang N, Wang H, Lu SE. Biosynthesis of an antifungal oligopeptide in Burkholderia contaminans strain MS14. Biochem Biophys Res Commun. 2009;380(2):328–32. 10.1016/j.bbrc.2009.01.073 . PubMed DOI

Lu SE, Novak J, Austin FW, Gu G, Ellis D, Kirk M, et al. Occidiofungin, a unique antifungal glycopeptide produced by a strain of Burkholderia contaminans. Biochemistry. 2009;48(35):8312–21. 10.1021/bi900814c PubMed DOI PMC

Gu G, Smith L, Liu A, Lu SE. Genetic and biochemical map for the biosynthesis of occidiofungin, an antifungal produced by Burkholderia contaminans strain MS14. Appl Environ Microbiol. 2011;77(17):6189–98. 10.1128/AEM.00377-11 PubMed DOI PMC

Gu G, Wang N, Chaney N, Smith L, Lu SE. AmbR1 is a key transcriptional regulator for production of antifungal activity of Burkholderia contaminans strain MS14. FEMS Microbiol Lett. 2009;297(1):54–60. 10.1111/j.1574-6968.2009.01653.x . PubMed DOI

Thomson EL, Dennis JJ. A Burkholderia cepacia complex non-ribosomal peptide-synthesized toxin is hemolytic and required for full virulence. Virulence. 2012;3(3):286–98. 10.4161/viru.19355 PubMed DOI PMC

Chapalain A, Vial L, Laprade N, Dekimpe V, Perreault J, Déziel E. Identification of quorum sensing-controlled genes in Burkholderia ambifaria. Microbiologyopen. 2013;2(2):226–42. 10.1002/mbo3.67 PubMed DOI PMC

Sass AM, Schmerk C, Agnoli K, Norville PJ, Eberl L, Valvano MA, et al. The unexpected discovery of a novel low-oxygen-activated locus for the anoxic persistence of Burkholderia cenocepacia. ISME J. 2013;7(8):1568–81. 10.1038/ismej.2013.36 PubMed DOI PMC

Lin Z, Falkinham JO, Tawfik KA, Jeffs P, Bray B, Dubay G, et al. Burkholdines from Burkholderia ambifaria: antifungal agents and possible virulence factors. J Nat Prod. 2012;75(9):1518–23. 10.1021/np300108u . PubMed DOI

Gomez-Escribano JP, Castro JF, Razmilic V, Chandra G, Andrews B, Asenjo JA, et al. The Streptomyces leeuwenhoekii genome: de novo sequencing and assembly in single contigs of the chromosome, circular plasmid pSLE1 and linear plasmid pSLE2. BMC Genomics. 2015;16:485 10.1186/s12864-015-1652-8 PubMed DOI PMC

Lieberman TD, Michel JB, Aingaran M, Potter-Bynoe G, Roux D, Davis MR, et al. Parallel bacterial evolution within multiple patients identifies candidate pathogenicity genes. Nat Genet. 2011;43(12):1275–80. 10.1038/ng.997 PubMed DOI PMC

Lieberman TD, Flett KB, Yelin I, Martin TR, McAdam AJ, Priebe GP, et al. Genetic variation of a bacterial pathogen within individuals with cystic fibrosis provides a record of selective pressures. Nat Genet. 2014;46(1):82–7. 10.1038/ng.2848 PubMed DOI PMC

Dillon MM, Sung W, Lynch M, Cooper VS. The rate and molecular spectrum of spontaneous mutations in the GC-rich multichromosome genome of Burkholderia cenocepacia. Genetics. 2015;200(3):935–46. 10.1534/genetics.115.176834 . PubMed DOI PMC

Turner KH, Wessel AK, Palmer GC, Murray JL, Whiteley M. Essential genome of Pseudomonas aeruginosa in cystic fibrosis sputum. Proc Natl Acad Sci U S A. 2015;112(13):4110–5. 10.1073/pnas.1419677112 PubMed DOI PMC

Barth AL, Pitt TL. The high amino-acid content of sputum from cystic fibrosis patients promotes growth of auxotrophic Pseudomonas aeruginosa. J Med Microbiol. 1996;45(2):110–9. . PubMed

Barth AL, Pitt TL. Auxotrophy of Burkholderia (Pseudomonas) cepacia from cystic fibrosis patients. J Clin Microbiol. 1995;33(8):2192–4. PubMed PMC

Lee H, Popodi E, Tang H, Foster PL. Rate and molecular spectrum of spontaneous mutations in the bacterium Escherichia coli as determined by whole-genome sequencing. Proc Natl Acad Sci U S A. 2012;109(41):E2774–83. 10.1073/pnas.1210309109 PubMed DOI PMC

Williams AB. Spontaneous mutation rates come into focus in Escherichia coli. DNA Repair (Amst). 2014;24:73–9. 10.1016/j.dnarep.2014.09.009 . PubMed DOI

Martina P, Feliziani S, Juan C, Bettiol M, Gatti B, Yantorno O, et al. Hypermutation in Burkholderia cepacia complex is mediated by DNA mismatch repair inactivation and is highly prevalent in cystic fibrosis chronic respiratory infection. Int J Med Microbiol. 2014;304(8):1182–91. 10.1016/j.ijmm.2014.08.011 . PubMed DOI

Sousa AM, Pereira MO. Pseudomonas aeruginosa diversification during infection development in cystic fibrosis lungs-a review. Pathogens. 2014;3(3):680–703. 10.3390/pathogens3030680 PubMed DOI PMC

Oliver A, Mena A. Bacterial hypermutation in cystic fibrosis, not only for antibiotic resistance. Clin Microbiol Infect. 2010;16(7):798–808. 10.1111/j.1469-0691.2010.03250.x . PubMed DOI

Oliver A. Mutators in cystic fibrosis chronic lung infection: Prevalence, mechanisms, and consequences for antimicrobial therapy. Int J Med Microbiol. 2010;300(8):563–72. 10.1016/j.ijmm.2010.08.009 . PubMed DOI

O'Sullivan LA, Weightman AJ, Jones TH, Marchbank AM, Tiedje JM, Mahenthiralingam E. Identifying the genetic basis of ecologically and biotechnologically useful functions of the bacterium Burkholderia vietnamiensis. Environ Microbiol. 2007;9(4):1017–34. 10.1111/j.1462-2920.2006.01228.x . PubMed DOI

Zerbino DR, Birney E. Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008;18(5):821–9. 10.1101/gr.074492.107 PubMed DOI PMC

Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, et al. Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28(12):1647–9. 10.1093/bioinformatics/bts199 PubMed DOI PMC

Benson DA, Karsch-Mizrachi I, Lipman DJ, Ostell J, Sayers EW. GenBank. Nucleic Acids Res. 2009;37(Database issue):D26–31. 10.1093/nar/gkn723 PubMed DOI PMC

Kaas RS, Leekitcharoenphon P, Aarestrup FM, Lund O. Solving the problem of comparing whole bacterial genomes across different sequencing platforms. PLoS One. 2014;9(8):e104984 10.1371/journal.pone.0104984 PubMed DOI PMC

Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9. 10.1093/molbev/mst197 PubMed DOI PMC

Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20. 10.1093/bioinformatics/btu170 PubMed DOI PMC

Kopylova E, Noé L, Touzet H. SortMeRNA: fast and accurate filtering of ribosomal RNAs in metatranscriptomic data. Bioinformatics. 2012;28(24):3211–7. 10.1093/bioinformatics/bts611 . PubMed DOI

Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41(Database issue):D590–6. 10.1093/nar/gks1219 PubMed DOI PMC

Burge SW, Daub J, Eberhardt R, Tate J, Barquist L, Nawrocki EP, et al. Rfam 11.0: 10 years of RNA families. Nucleic Acids Res. 2013;41(Database issue):D226–32. 10.1093/nar/gks1005 PubMed DOI PMC

Wu TD, Nacu S. Fast and SNP-tolerant detection of complex variants and splicing in short reads. Bioinformatics. 2010;26(7):873–81. 10.1093/bioinformatics/btq057 PubMed DOI PMC

Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5(10):R80 10.1186/gb-2004-5-10-r80 PubMed DOI PMC

Lawrence M, Huber W, Pagès H, Aboyoun P, Carlson M, Gentleman R, et al. Software for computing and annotating genomic ranges. PLoS Comput Biol. 2013;9(8):e1003118 10.1371/journal.pcbi.1003118 PubMed DOI PMC

Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550 10.1186/PREACCEPT-8897612761307401 . PubMed DOI PMC

Boratyn GM, Camacho C, Cooper PS, Coulouris G, Fong A, Ma N, et al. BLAST: a more efficient report with usability improvements. Nucleic Acids Res. 2013;41(Web Server issue):W29–33. 10.1093/nar/gkt282 PubMed DOI PMC

Winsor GL, Khaira B, Van Rossum T, Lo R, Whiteside MD, Brinkman FS. The Burkholderia Genome Database: facilitating flexible queries and comparative analyses. Bioinformatics. 2008;24(23):2803–4. 10.1093/bioinformatics/btn524 PubMed DOI PMC

Sokol PA, Ohman DE, Iglewski BH. A more sensitive plate assay for detection of protease production by Pseudomonas aeruginosa. J Clin Microbiol. 1979;9(4):538–40. PubMed PMC

Chen KC, Ravichandran A, Guerrero A, Deng P, Baird SM, Smith L, et al. The Burkholderia contaminans MS14 ocfC gene encodes a xylosyltransferase for production of the antifungal occidiofungin. Appl Environ Microbiol. 2013;79(9):2899–905. 10.1128/AEM.00263-13 PubMed DOI PMC