Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Acinetobacter baumannii Virulence Is Mediated by the Concerted Action of Three Phospholipases D

  • Julia Stahl,

    Affiliation Department of Molecular Microbiology and Bioenergetics, Institute of Molecular Biosciences, Goethe University Frankfurt, Frankfurt am Main, Germany

    ⨯
  • Holger Bergmann,

    Affiliation Department of Applied Bioinformatics, Institute for Cell Biology and Neuroscience, Goethe University Frankfurt, Frankfurt am Main, Germany

    ⨯
  • Stephan Göttig,

    Affiliation Institute of Medical Microbiology and Infection Control, University Hospital Frankfurt, Goethe University, Frankfurt am Main, Germany

    ⨯
  • Ingo Ebersberger,

    Affiliation Department of Applied Bioinformatics, Institute for Cell Biology and Neuroscience, Goethe University Frankfurt, Frankfurt am Main, Germany

    ⨯
  • Beate Averhoff

    averhoff@bio.uni-frankfurt.de

    Affiliation Department of Molecular Microbiology and Bioenergetics, Institute of Molecular Biosciences, Goethe University Frankfurt, Frankfurt am Main, Germany

    ⨯

Abstract

Acinetobacter baumannii causes a broad range of opportunistic infections in humans. Its success as an emerging pathogen is due to a combination of increasing antibiotic resistance, environmental persistence and adaptation to the human host. To date very little is known about the molecular basis of the latter. Here we demonstrate that A. baumannii can use phosphatidylcholine, an integral part of human cell membranes, as sole carbon and energy source. We report on the identification of three phospholipases belonging to the PLD superfamily. PLD1 and PLD2 appear restricted to the bacteria and display the general features of bacterial phospholipases D. They possess two PLDc_2 PFAM domains each encompassing the HxKx4Dx6GS/GGxN (HKD) motif necessary for forming the catalytic core. The third candidate, PLD3, is found in bacteria as well as in eukaryotes and harbours only one PLDc_2 PFAM domain and one conserved HKD motif, which however do not overlap. Employing a markerless mutagenesis system for A. baumannii ATCC 19606T, we generated a full set of PLD knock-out mutants. Galleria mellonella infection studies as well as invasion experiments using A549 human lung epithelial cells revealed that the three PLDs act in a concerted manner as virulence factors and are playing an important role in host cell invasion.

Citation: Stahl J, Bergmann H, Göttig S, Ebersberger I, Averhoff B (2015) Acinetobacter baumannii Virulence Is Mediated by the Concerted Action of Three Phospholipases D. PLoS ONE 10(9): e0138360. https://doi.org/10.1371/journal.pone.0138360

Editor: Eric Cascales, Centre National de la Recherche Scientifique, Aix-Marseille Université, FRANCE

Received: June 16, 2015; Accepted: August 28, 2015; Published: September 17, 2015

Copyright: © 2015 Stahl et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Data Availability: All relevant data are within the paper and its Supporting Information files.

Funding: This study was supported by a grant by the German Research Council (DFG) - FOR2251 (www.acinetobacter.de).

Competing interests: The authors have declared that no competing interests exist.

Introduction

Acinetobacter baumannii is one of the most successful and most threatening nosocomial (hospital-acquired) Gram-negative bacteria implicated in septicaemia, meningitis, wound infections, urinary tract infections and pneumonia [1,2]. In recent years, the National Nosocomial Infection Surveillance System (NNIS) revealed a substantial increase in A. baumannii-associated nosocomial pneumonia cases. Overall, they make up 5–10% of the intensive care unit (ICU)-acquired pneumonia cases in the United States [3–6]. Moreover, A. baumannii ranks 10th among the organisms causing monomicrobial bloodstream infections [1,4]. Besides its steady increase in resistance against a broad spectrum of antibiotics the pronounced resistance to desiccation of A. baumannii is also highly worrying [1,7,8]. The persistence and virulence of A. baumannii can be attributed to a broad spectrum of factors. These include biofilm formation, quorum sensing, resistance to desiccation, possession of powerful iron acquisition systems and the production of virulence factors such as phospholipases C and D [8–14]. In addition, A. baumannii displays a number of traits conveying adaptation to the human as a host. For example, it can adhere to, colonize and invade human epithelial and endothelial cells and displays human serum resistance [8–10,14–17]. However, despite the increasing threat of A. baumannii surprisingly little is known about the molecular mechanisms of host adaptation.

The utilization of carbon sources available in the host is one of the prerequisites for the persistence of pathogens in the human host. The human body is a rich source of nutrients and it has been shown that different pathogens use different carbon sources during host colonization, such as lactic acid, mucus sugars, amino acids and cholesterol, just to mention of few [18]. Phospholipids are also very abundant in the human host since they are the major building blocks of biological membranes making them also good candidates as carbon and energy source. Phosphatidylcholine (PC) is particularly abundant in eukaryotic membranes accounting for 50% of all phospholipids [19]. In the lung the prevalence even increases up to 80% [20] and PC is also the major phospholipid in tracheobronchial secretions where it has significant impact in wettability of the mucus [21]. There is strong evidence that phosphatidylcholine serves as nutrient source during lung-infections by pathogens like Pseudomonas aeruginosa [22]. Taking into account that ventilator-associated pneumonia is the most prevalent manifestation of A. baumannii infection [23], high abundance of PC in the mucus raised the question whether PC serves as nutrient source for A. baumannii in the human host.

Phospholipases (PL) are key enzymes, essential for the metabolism of PC and have been studied in a diverse spectrum of bacteria [24–26]. Three classes of phospholipases (PLA, PLC, PLD) have been defined by the site of cleavage. PLA hydrolyzes fatty acids from the glycerol backbone. Upon PLC cleavage, the phosphorylated head group is released from the phospholipid and PLD cleaves off only the head group. Both, release of the phosphorylated headgroup and release of the polar head group can affect the stability of host cell membranes. In addition, phospholipases can also interfere with cellular signalling by generating second messengers like phosphatidic acid, which are capable of modulating host immune responses [24–26]. Information about the distribution of phospholipases in different A. baumannii strains and their role in virulence is scarce. So far, only two phospholipases, one PLC and one PLD have been identified as virulence factors in A. baumannii ATCC 17978 and in the clinical isolate 98-37-09, respectively [11,12]. Mutation studies revealed an involvement in cytotoxicity and epithelial cell invasion, respectively, and suggested a potential role of the two enzymes in virulence.

In this study, we identified three pld genes in A. baumannii ATCC 19606T and performed a thorough bioinformatic and experimental characterization of the corresponding proteins. We show that pld1 and pld2 evolved through a recent gene duplication, most probably in the last common ancestor of Acinetobacter. In contrast, pld3 appears substantially older and we provide initial evidence that it may share evolutionary descent with the cardiolipin hydrolase ZUC in Drosophila. Furthermore, PLD3 shows an exchange of a highly conserved aspartate by a tyrosine in the HKD active site motif of the PLDc_2 domain, otherwise a hallmark of the PLD superfamily. For functional characterization of the PLDs by mutation studies, we developed a markerless mutagenesis system to generate a whole set of pld knock-out mutants. Bacterial competition studies, infection of Galleria mellonella and invasion of A549 human lung epithelial cells revealed that the three PLDs do not play a role in bacterial competition but act in a concerted manner in infection and invasion of eukaryotic host cells.

Materials and Methods

Organisms and Cultivation

A. baumannii ATCC 19606T was grown at 37°C in Luria Bertani medium (LB) or mineral medium (MM) with 20 mM acetate as carbon source as described earlier [27]. Kanamycin was added when appropriate to a final concentration of 50 μg/ml. E. coli DH5α was grown at 37°C in LB medium. When appropriate, kanamycin or ampicillin were added to a final concentration of 20 μg/ml and 100 μg/ml, respectively. Pseudomonas putida 548.C8 [28] was grown in LB medium with 20 μg/ml kanamycin at 30°C.

Growth Experiments

Growth experiments with A. baumannii ATCC 19606T wildtype and mutants were performed using 500 ml Erlenmeyer flasks filled with 100 ml medium (LB or MM + 20 mM acetate) and incubated shaking (130 rpm) at 37°C. Growth experiments were started by inoculation of a stationary phase culture into fresh medium to OD600 of 0.15. Growth was followed at 37°C by analyzing the optical density at 600 nm.

For growth on different carbon sources, A. baumannii was inoculated in a test tube containing 5 ml MM with either 20 mM acetate, glycerol or choline or MM with either 0.5% PC or 0.5% palmitate as sole carbon source. Growth experiments were started by 1:10 dilution of an appropriate overnight pre-culture, grown on the corresponding carbon source. Growth on acetate, glycerol and choline was monitored by increase of the OD600 over time. Growth on PC or palmitate was determined by plating serial dilutions on LB agar and CFU counting at different time points. In each case, MM without any carbon source was used as negative control, whereas growth on acetate was used as positive control.

Data Collections and Data Sources

To determine the phyletic profile of PLDs across the entire tree of life we first compiled the TOL set comprising protein collections from 1281 bacteria, 123 archaea and 204 eukaryotes that are available from the Orthology Matrix project web sites (OMA, Release 16, Mach 2014, http://cbrg-oma-test.ethz.ch/oma/archives/). These data were complemented by the set of proteins encoded in the genome of the A. baumannii type strain ATCC 19606T as provided by NCBI GenBank (http://www.ncbi.nlm.nih.gov/genome/browse/representative/). Our second data set (ACI set) was designed to provide a high-resolution presence-absence pattern of PLDs in the genus Acinetobacter. The ACI set comprises the proteins from the 14 Acinetobacter species represented in the TOL set together with the proteins from 54 additional Acinetobacter strains provided by NCBI GenBank (http://www.ncbi.nlm.nih.gov/genome/browse/representative/). The full taxon lists of both the TOL (S1 Table) set and the ACI (S2 Table) set are available as supporting information. Please note that putative products of annotated pseudogenes are not explicitly considered in the analyses.

Annotation of Functional Domain Content

We annotated Pfam-A domains [29] in the PLD proteins employing hmmscan from the HMMER3 package [30] (http://hmmer.janelia.org). Further functional features and signatures in the protein sequence were annotated with interproscan (http://www.ebi.ac.uk/Tools/pfa/iprscan5/) using default settings. Domain architectures were visualized and arranged with Domosaics v. 0.95 [31]. The presence of signal sequences was analyzed using SinglaP 4.1 [32].

Ortholog Search

For the detection of orthologs in both the TOL and the ACI set we used the targeted ortholog search tool HaMStR v. 13.2 [33] (http://sourceforge.net/projects/hamstr). We seeded the ortholog searches for the three proteins, PLD1, PLD2, and PLD3, respectively, with the corresponding protein sequences from A. baumannii ATCC 19606T. We then used an initial HaMStR search to first compile a manually curated core set of orthologs for each of the three proteins. In the curation step we verified that the domain architecture of the identified orthologs resembles that of the respective seed proteins. The final core sets comprised proteins from the following bacterial species: Paenibacillus terrae strain HPL-003, Yersinia pestis strain CO92, Chlamydia trachomatis strain B, Acinetobacter calcoaceticus strain PHEA-2, Helicobacter pylori strain PeCan4 and A. baumannii strain ATCC 19606T. Subsequently, we aligned the sequences in the three core sets with MAFFT-LINSI [34] and used the alignments to train a profile hidden Markov model (pHMM) with hmmbuild from the HMMER3 package [30]. The resulting refined pHMM served then as input for the second ortholog search in the full compilation of bacterial, archaeal and eukaryotic gene sets.

Phylogenetic Analysis

We aligned the amino acid sequences in the individual orthologous groups using MAFFT-LINSI [34]. For the resulting alignments, ProtTest 3.2 [35] identified the VT model [36] modeling rate across sites with a gamma distribution (+G) as the best-fitting substitution model. Maximum likelihood tree reconstruction was performed with RAxML 8.1.9 [37] using the PROTGAMMAVT model for amino acid sequence evolution. For statistical support, 100 nonparametric bootstrap replicates were computed. A combined tree for all three ortholog groups was reconstructed using FastTree 2.1.7 [38], and branch support was assessed with a local Shimodaira-Hasegawa test [39] as implemented in FastTree.

Markerless Mutagenesis

A sacB/kanR-cassette was amplified from the genome of A. baylyi ADP1 hptA::sacB/kanR [40] using the primer combination sacB/kanR_for and sacB/kanR_rev (S3 Table). A 1825 bp fragment spanning the pld2-gene (1632 bp), as well as 45 bp of the up- and 145 bp of the downstream region, respectively, was amplified from the genome of A. baumannii ATCC 19606T using pld2_up_for and pld2_do_rev primers. The fragment was cloned into KpnI and PstI sites of pBIISK (Invitrogen). A central 316 bp fragment of pld2 was deleted by restriction with EcoRI and HincII followed by Klenow fill-in and subsequent ligation, resulting in plasmid pBIISK_pld2_updown. The sacB/kanR cassette was inserted into the SmaI-site of this plasmid by blunt end ligation, resulting in pBIISK_sacB/kanR_pld2_updown. To generate the constructs for pld1 and pld3-mutation the sacB/kanR cassette was inserted into the EcoRV site of pBIISK (Invitrogen) by blunt end ligation leading to pBIISK_sacB/kanR. For each construct ~1500 bp of the up- and downstream region of pld1 and pld3 were amplified from the genome of A. baumannii ATCC 19606T using the following primer combinations: pld1_up_for and pld1_up_rev (pld1 upstream-region), pld1_do_for and pld1_do_rev (pld1 downstream region), pld3_up_for and pld3_up_rev (pld3 upstream region), pld3_do_for and pld3_do_rev (pld3 downstream region). Fragments containing up- and downstream regions of pld1 and pld3 were cloned into pBIISK_sacB/kanR using PstI, BamHI and NotI resulting in the vectors pBIISK_sacB/kanR_pld1_updown and pBIISK_sacB/kanR_pld3_updown.

To generate electro-competent A. baumannii cells, an overnight culture of A. baumannii ATCC 19606T was inoculated into 400 ml fresh LB medium and grown until OD600 of 0.45–0.5. The cultures were cooled down on ice for 10 minutes and washed 4 times with sterile, ice-cold mili-Q H2O, followed by one washing step in sterile, ice-cold 10% glycerol. Cell pellets were resuspended in 500 μl 10% glycerol. 50 μl of the A. baumannii cell suspensions were used as recipient for electroporation of plasmid DNA (2–3 μg DNA). Electroporation was performed at 2.5 kV, 200 Ώ and 25 μF in 2 mm pre-cooled electroporation cuvettes. After electroporation, 1 ml pre-warmed SOC medium (2% tryptone; 0.5% yeast extract; 10 mM NaCl; 2.5 mM KCl; 10 mM MgCl2; 10 mM MgSO4; 20 mM Glucose) was added. Cells were incubated for 1 hour at 37°C, and subsequent transformants were selected on LB-Agar + 50 μg/ml kanamycin. Integration of the plasmid was monitored by colony PCR using the primer combinations ctr_for + do_rev and up_for + ctr_rev (S3 Table).

Correct integrants were subjected to counter-selection overnight at 37°C in 5 ml LB medium with 10% sucrose. Appropriate dilutions of the cultures were plated onto LB agar medium with 10% sucrose and incubated at 37°C for 12 h. Single colonies were tested for kanamycin sensitivity by replica plating onto LB agar ± 50 μg/ml kanamycin. Kanamycin sensitive colonies were subjected to PCR analyzes using the primers ctr_for + ctr_rev (S3 Table). Double and triple mutants were obtained by insertion and segregation of the plasmids comprising the truncated phospholipase genes of interest into single or double pld mutants.

Interbacterial Competition Assay

The interbacterial competition between wild-type (WT) or pld-mutants and E. coli DH5α or P. putida, respectively, was analyzed as described earlier for competition between A. nosocomialis and E. coli [41,42]. A. baumannii and E. coli DH5α pET28a were grown at 37°C in 20 ml LB medium to an OD600 of 1.5–2.0. P. putida 548.C8 [28] was grown at 30°C in 20 ml LB medium to an OD600 of 1.5–2.0. Subsequently, OD600 was adjusted to 0.4 (~ 2 * 107 CFU/ml) and A. baumannii and E. coli DH5α cells or A. baumannii and P. putida cells were mixed in a 10:1 ratio. Twenty μl of this mixture were spotted onto LB agar and incubated at 37°C for 4 hours. Spots were excised from the agar plate and resuspended in 500 μl 0.9% NaCl. Serial dilutions were plated onto LB agar + 50 μg/ml kanamycin to select for surviving E. coli or P. putida cells.

Galleria mellonella Infection Assay

Galleria mellonella caterpillars (larvae of the greater wax moth) were obtained by a local provider. Caterpillars were preselected by melanization, size and movement in response to touch. A. baumannii wild-type and pld mutant cells were grown on LB medium until late exponential growth phase, washed with 0.9% NaCl and adjusted to an OD600 of 2. In each set of infection, 16 caterpillars were used to test each bacterial strain in one test. Caterpillars were whipped in 70% EtOH for external disinfection and 10 μl (~ 1 * 106 CFU/ml) of the bacterial suspension were injected into the last left proleg. As control served a set of untreated caterpillars as well as a set of caterpillars in which 10 μl of 0.9% NaCl were injected. Caterpillars were incubated at 37°C in the dark for 4 days. Caterpillars were considered as dead if they did not respond to gentle probing. All experiments were repeated at least 4 times and experiments in which 2 or more caterpillars in one of the control groups died were not considered.

Adhesion and Invasion Assays

Human lung epithelial A549 (ATCC CCL-185) cells were grown to confluence in six well plates in RPMI 1640 medium (Merck Millipore, Darmstadt, Germany) supplemented with 10% fetal calf serum (FCS) at 37°C and 5% CO2. A. baumannii was grown in LB medium overnight at 37°C. Stationary cultures were washed once in RPMI medium prior to addition of 5x107 bacteria to each well, corresponding to a multiplicity of infection (MOI) of approximately 200 bacteria per A549 cell. The plates were incubated for 5 hours at 37°C in 5% CO2 to allow adhesion and invasion of bacteria, as described earlier [43]. For analysis of invasion, A549 cells were washed three times with Dulbecco’s phosphate-buffered saline (PBS). Afterwards, 2 ml RPMI containing 500 μg/ml gentamicin were added to each well for 2 hours to eliminate all extracellular bacteria. After three washing steps with PBS, A549 cells were lyzed by addition of 800 μl distilled H2O for 15 min, to release invaded bacteria from the disrupted A549 cells. Then, 200 μl 5x PBS was added and serial dilutions were plated on LB-Agar for determination of CFU/ml. For analysis of adhesion, A549 cells were infected as described above, washed three times with PBS and lyzed with H2O to release all internal and adherent bacteria. All invasion and adhesion assays were performed in duplicates and were repeated at least three times.

Statistical Analysis

Paired t-tests were performed with GraphPad Prism 6 (GraphPad Software Inc., San Diego, CA, USA). Values are presented as mean ± SEM. Growth rate was determined with GraphPad Prism 6 using exponential growth equation.

Results

Use of Phosphatidylcholine as sole Carbon and Energy Source

Growth experiments on mineral medium (MM) containing phosphatidylcholine revealed that PC serves as carbon and energy source for A. baumannii ATCC 19606T with a growth rate of 0.378 ± 0.028 h-1 and a final cell density of ~ 1.4∙108 CFU/ml (Fig 1). To elucidate which of the PC cleavage products were used as carbon source by A. baumannii we tested growth in MM with glycerol, fatty acids and choline as sole carbon and energy sources. Neither glycerol nor choline alone could serve as sole carbon and energy source since no increase in optical density was observed after incubation for 2 days. In contrast, palmitate, which is one of the most common fatty acids in PC, was used as carbon and energy source by A. baumannii with a growth rate of 0.545 ± 0.057 h-1 and a final cell density of ~ 5∙108 CFU/ml (Fig 1). As positive control for all growth experiments served MM with 20 mM acetate, resulting in a growth rate of 0.586 ± 0.001 h-1 and a final cell density of ~ 1.1∙109 CFU/ml. As negative control served MM without carbon source. No increase in optical density as well as CFU was observed after the first transfer in MM without carbon source.

thumbnail
Download:
Fig 1. Growth of A. baumannii ATCC 19606T on PC, palmitate and acetate.

A. baumannii ATCC 19606T was grown in 5 ml MM with either 0.5% PC (triangles), 0.5% palmitate (squares) or 20 mM acetate (circles) as sole carbon source. A culture without carbon source (diamonds) was used as negative control. Growth was followed by plating appropriate dilutions on LB-agar plates at the indicated time points. The amount of colony forming units (CFU) was determined after incubation at 37°C overnight.

https://doi.org/10.1371/journal.pone.0138360.g001

Screen for Phospholipases D in the Genome of A. baumannii

In previous studies, we could show that Acinetobacter baylyi ADP1 oxidizes choline after growth on MM containing choline and identified a potential choline oxidation pathway in the genome of A. baylyi ADP1 [44]. A homologous pathway, including a choline dehydrogenase (BetA), a betaine-aldehyde dehydrogenase (BetB), a potential regulator (BetI) as well as two choline transporters (BetT1 and BetT2) was also detected in the genome of A. baumannii ATCC 19606T. Therefore, we concluded, that A. baumannii can also oxidize choline [44], raising the question whether A. baumannii exhibits PLDs, which catalyze the liberation of choline from PC. We screened the genome of A. baumannii ATCC 19606T for the presence of potential pld genes and identified three candidates: The open reading frame HMPREF0010_00607 (pld1) spans 1527 bp and encodes a hypothetical PLD1 protein with a deduced molecular mass of 58.8 kDa. HMPREF0010_03706 (pld2) spans 1632 bp and encodes a hypothetical phospholipase PLD (PLD2) with a molecular mass of 62.4 kDa. PLD2 is identical to the recently published PLD of A. baumannii 98-37-09 [12]. The third candidate, HMPREF0010_02731 (pld3), has a predicted size of 1464 bp and encodes a hypothetical 55 kDa protein.

Subsequently, we characterized the putative PLDs in silico. A screen for conserved domains and signatures in the three proteins revealed that PLD1 and PLD2 each harbor two phospholipase D-like domains PLDc_2 (PF13091). Embedded in both domains we detected the 28 amino acid long phospholipase D phosphodiesterase active site motif (PF00614; Fig 2). This motif spans the HxKx4Dx6GS/GGxN pattern that has been described as the catalytic center of the phospholipases [45]. Interestingly, the third candidate, PLD3, deviates from this domain architecture. It features only one phospholipase D-like domain and one active site motif, which however do not overlap. While the active site is located in the central part, the PLDc_2 domain is located C-terminal of PLD3. Notably, a subsequent screen of the PLD3 sequence with reduced stringency revealed a second HKD-like motif overlapping with the C-terminal PLDc_2 domain. However the typically highly conserved aspartic acid (D) is replaced in this motif by a tyrosine (Y) (Fig 2). Our further annotation revealed that none of the pld genes encodes a signal sequence indicative for protein export. Furthermore PLD1 lacks a transmembrane region, whereas PLD2 and PLD3 both possess a hydrophobic domain which is predicted to represent a membrane anchor.

thumbnail
Download:
Fig 2. Domain architecture of the three A. baumannii ATCC 19606T phospholipase D candidates.

PLD1 and PLD2 each possess PLDc_2 Pfam domains (PF13091). Each domain harbors the conserved motif characteristic for the active center of type D phospholipases. PLD3 has only a single PLD domain encompassing a HKD-like motif (HxKx4Dx6GS/GGxN) where the functionally relevant aspartate is replaced by a tyrosine. A second canonical HKD motif is present in the N-terminal half of the protein and outside the context of a PLDc_2 domain. While both, PLD2 and PLD3, appear membrane bound proteins, PLD1 lacks a transmembrane domain (TM).

https://doi.org/10.1371/journal.pone.0138360.g002

Phyletic Distribution and Evolutionary History of A. baumannii Phospholipases D

First, we analyzed the distribution of PLD candidates in the genus Acinetobacter by a targeted ortholog search in the 68 bacterial strains of the ACI set. PLD1 orthologs are present in 60 taxa, while PLD2 and PLD3 orthologs are present in all of the 68 analyzed taxa (S2 Table). This observation indicates that PLD1 can be absent in individual strains of the genus Acinetobacter due to gene loss, a loss of function, or pseudogenization—the latter can be observed for A. baumannii strain AYE and its PLD1 ortholog. Please note, that the two other PLD1 pseudogenes described by Vallenet et al. 2008 [46] are not represented in the AYE protein sets downloaded from the OMA web pages. The resulting presence/absences matrix is given in S2 Table.

Second, we characterized the phylogenetic profile and the evolutionary relationships of the three PLD candidates. To this end, we first performed a targeted ortholog search in the TOL set comprising proteins from 1,281 bacteria, 124 archaea, and 204 eukaryotes. The results are summarized in S4 Table and the detected orthologs are provided in the S1 File. Both PLD1 and PLD2 are widely distributed across bacteria with orthologs present in 498 and 436 taxa, respectively. In archaea and eukaryotes orthologs could be identified in only 10 (PLD1) and 3 (PLD2) taxa respectively. Given the sparse and scattered phyletic distribution of these candidates the orthology assignments are most likely spurious. This suggests that both proteins are confined to bacteria. The phylogenetic profile of the third protein, PLD3, differs substantially from that of the other two PLDs. Though orthologs were also found throughout the analyzed bacterial clades, PLD3 is with only 187 detected orthologs less common in bacteria than both, PLD1 and PLD2. However, and again in contrast to the other two proteins, its prevalence in the eukaryotes, and here particularly in the animals is marked (S4 Table). Analyses of the phylogenetic tree revealed that orthologs assigned to the A. baumannii PLD1 and PLD2 genes, respectively, are well mixed in the tree (Fig 3). Such an intermingled topology indicates a rather recent evolutionary separation of the corresponding A. baumannii genes. A subsequent maximum likelihood tree reconstruction with a taxon sampling tailored to shed further light on this issue confirmed this hypothesis and narrowed the gene duplication event to have happened within the Moraxellaceae lineage (Fig 4 and S5 Table). In contrast to the entwined relationships of PLD1 and PLD2 orthologs, our PLD tree also indicates that, with very few exceptions, all PLD3 orthologs from the three domains of life reside in a single clade to the exclusion of the other PLDs (Fig 3). Determination of the Pfam domain architecture revealed that the vast majority (178) of PLD3 orthologs share the presence of only a single PLDc_2 domain as a common feature (S1 Fig and S3 File). Moreover, archaean sequences are distributed across the reconstructed phylogeny suggesting several independent evolutionary origins of these sequences. In contrast, the eukaryotic sequences form, with two exceptions, a monophyletic clade and the corresponding subtree resembles by and large the eukaryotic species phylogeny. This is consistent with a single evolutionary origin of the eukaryotic PLD3 sequences. As one possibility they may share a common descent with their bacterial orthologs from an evolutionarily ancient gene in the last universal common ancestor. Alternatively, they might have been horizontally introduced into the eukaryotic lineage. In this context, it is interesting to note that the best annotated animal ortholog, the Drosophila protein ZUC (Uniprot Id Q9VKD7), is associated with the mitochondria making an introduction of this gene into the eukaryotes via the endosymbiont conceivable.

thumbnail
Download:
Fig 3. Phylogenetic tree combining the orthologous groups of PLD1, -2, and -3.

The tree demonstrates the intertwined evolutionary relationships of orthologs assigned to A. baumannii PLD1 and PLD2, respectively. In contrast, PLD3 sequences form a well separated clade suggesting an evolutionary origin distinct from that of PLD1 and PLD2. The * denotes a local branch support of 0.83.

https://doi.org/10.1371/journal.pone.0138360.g003

thumbnail
Download:
Fig 4. Evolutionary relationships of Acinetobacter baumannii PLD1 and PLD2.

The tree reveals that PLD1 and PLD2 emerged through a recent gene duplication that occurred after the split of P. aeruginosa and possibly in the last common ancestor of the genus Acinetobacter. Orthologs whose split pre-date the diversification of PLD1 and 2 were renamed to PLD1/2. Branch labels represent percent bootstrap support.

https://doi.org/10.1371/journal.pone.0138360.g004

Generation and Fitness Analyses of pld Mutants

In order to facilitate the functional characterization of the various PLDs we established a two-step procedure for markerless mutagenesis of A. baumannii. This procedure is favorable since it avoids polar effects. In the first step we performed insertion duplication mutagenesis of the target pld genes using a recombinant pBIISK plasmid, which carries a kanamycin resistance gene (kanR) as positive selection marker, a levansucrase gene (sacB) for negative selection and 1500 bp of the upstream and downstream region flanking the pld gene of interest. With this technique a complete set of phospholipase D single (Δpld1, Δpld2, Δpld3), double (Δpld1+2, Δpld1+3, Δpld2+3) and triple (Δpld1-3) mutants was generated. All mutants were verified by sequencing using the primer combination ctr_for and ctr_rev (S3 Table).

We determined whether the knockout of the phospholipases had any general impact on the growth of the mutants by comparing the growth on LB and mineral medium (MM) of all pld mutants to that of the wild-type. Notably, the WT and the mutants reached a comparable final OD600 of 3.3 to 3.5 in LB-medium with a growth rate of 0.9 to 1.0 h-1 and a maximal OD600 of 1.2 to 1.3 in MM + 20 mM acetate, with a growth rate of 0.75 to 0.80 h-1, indicating that knock out of any of the plds had no effect on general fitness. Next we raised the question whether PLDs are required for growth on PC. Therefore we tested growth of the Δpld1-3 triple mutant on PC as sole carbon source. No difference in growth rate and final cell density was observed compared to the wildtype, leading to the conclusion, that these PLDs are not essential for growth on PC.

Phospholipases D are not essential for Interbacterial Competition

Recently, it was shown that PLDs of Pseudomonas aeruginosa are implicated in bacterial competition [47]. Therefore we addressed the question whether A. baumannii ATCC 19606T outcompetes other bacteria. We used E. coli as model strain in a bacterial competition assay as described earlier for competition between A. nosocomialis and E. coli [41,42]. A. baumannii ATCC 19606T wild-type cells were mixed with kanamycin resistant E. coli in a 10:1 ratio. A mixture of LB-medium and E. coli in the same ratio was used as negative control. After incubation of the mixture for 4 hours, the colony forming units of E. coli cells were determined by plating serial dilutions onto kanamycin containing agar plates. Incubation with A. baumannii ATCC 19606T wild-type revealed 1.9∙105 ± 0.8∙105 colony forming E. coli cells (Fig 5A). Incubation of E. coli cells with LB medium as control led to the detection of 4.9∙107 ± 0.9∙107 CFU. This provides clear evidence that A. baumannii ATCC 19606T outcompetes E. coli. To verify our results we repeated the competition analyses using a kanamycin resistant Pseudomonas putida strain 548.C8 [48]. 1.4x105 ± 0.5∙105 P. putida CFU were determined when incubated with A. baumannii wild-type cells, which clearly shows that A. baumannii also outcompetes P. putida (Fig 5B).

thumbnail
Download:
Fig 5. Bacterial competition between A. baumannii ATCC 19606T and E. coli or P. putida.

A. baumannii ATCC 19606T wild-type or pld mutant cells and E. coli DH5α pET28a cells (A) or P. putida 548.C8 cells (B) were mixed at a ratio of 10:1, spotted onto an LB-agar plates and incubated for 4 hours at 37°C. To quantify the number of surviving E. coli or P. putida cells serial dilutions were plated onto kanamycin containing LB-agar and colony forming units were counted after incubation overnight at 37°C or 30°C, respectively. The standard deviation was calculated from three independent experiments. LB medium was used instead of A. baumannii as negative control.

https://doi.org/10.1371/journal.pone.0138360.g005

Apparently, the pld single mutants were not affected in outcompeting E. coli (Fig 5A). The same was observed for the triple mutant which was also not affected in outcompeting E. coli or P. putida. 6∙105 ± 2.9∙105 E. coli and 1.4∙105 ± 0.7∙105 P. putida cells were counted after incubation with the pld triple mutant (Fig 5B). As control living cell counts of A. baumannii WT and Δpld1-3 mutant were analyzed by plating dilutions on LB agar + ampicillin for selection of A. baumannii. This control supported our finding since no difference in CFU/ml was observed between wildtype and Δpld1-3 mutant after competition with either E. coli or P. putida. ~ 2.6∙108 CFU/ml of WT or Δpld1-3 mutant survived when incubated with E. coli and ~ 2.1∙108 CFU/ml survived when incubated with P. putida. Our findings are consistent with the notion, that A. baumannii PLDs are not essential for antibacterial activity against E. coli or P. putida under the conditions tested.

Identification of Phospholipases D as Virulence Factors in vivo

The larvae of the greater wax moth (Galleria mellonella) have previously been described as a suitable model system to study Acinetobacter pathogenesis in vivo [49]. In particular, strain ATCC 19606T was shown to kill G. mellonella in a time and dose dependent manner [13]. To elucidate the role of A. baumannii PLDs in the infection process, G. mellonella caterpillars were infected with 106 cells of A. baumannii ATCC 19606T wild-type or the Δpld1-3 mutant by injecting the bacterial cell suspensions (10 μl, OD600 = 2) into the hemolymph via the last left proleg. Infected caterpillars were incubated at 37°C in the dark. The A. baumannii wild-type killed the caterpillars in a time dependent manner (Fig 6). 52.8 ± 7.2%, 75.9 ± 4.8%, 84.4 ± 3.1% and 89.8 ± 2.6% of the caterpillars were dead 1, 2, 3 and 4 days after infection. In contrast infection with the pld triple mutant led to a significant reduction in the mortality of the caterpillars compared to the wild-type (paired t-test; p<0.05 for day 1; p<0.01 for days 2 to 4). Only 31.9 ± 4.9%, 35.7 ± 5.6%, 51.56 ± 6.4% and 74.2 ± 3.6% of larvae were dead after 1, 2, 3 and 4 days of infection, respectively (Fig 6). Subsequently, we addressed the question whether the observed virulence deficiency can be linked to the inactivation of individual phospholipases D. We repeated the G. mellonella infection assay with all single and double pld mutants. Interestingly, caterpillars infected with either pld single or double mutants died with a similar rate as caterpillars infected with the wild-type (Fig 7). The finding that only the pld triple mutant exhibited reduced virulence in the Galleria infection implies that the concerted action of all three phospholipases are necessary for full virulence of A. baumannii ATCC 19606T in a Galleria model.

thumbnail
Download:
Fig 6. A. baumannii Δpld1-3 mutant is less virulent in G. mellonella infection.

For infection assays 16 caterpillars were injected with approximately 106 bacteria and incubated at 37°C in the dark. Mean difference between the wild-type and the Δpld1-3 mutant in four independent experiments was significant after 1, 2, 3 and 4 days (paired t-test; p<0.05 for day 1, p<0.01 for days 2, 3 and 4). As control served an untreated set of caterpillars as well as a set of caterpillars which was treated with 10 μl of sterile 0.9% NaCl. Tests in which more than 2 caterpillars in one of the control groups died after 4 days were not considered.

https://doi.org/10.1371/journal.pone.0138360.g006

thumbnail
Download:
Fig 7. Single or double pld mutants display an unaffected virulence towards G. mellonella.

For infection assays 16 caterpillars were injected with approximately 106 bacteria and incubated at 37°C in the dark. Mean difference between wild-type and single (A) or double mutants (B) was not significant. As control served an untreated set of caterpillars as well as a set of caterpillars which was treated with 10 μl of sterile 0.9% NaCl. Tests in which more than 2 caterpillars in one of the control groups died after 4 days were not considered.

https://doi.org/10.1371/journal.pone.0138360.g007

Full Invasion Efficiency requires the concerted Action of all Three PLDs

To analyze the role of the PLDs in human host cell invasion, we tested the pld mutants for invasion efficiencies in a gentamycin protection assay, as described earlier [43]. Single pld1 and pld3 mutants had no effect on invasion efficiency, whereas a single pld2 mutant reduced the invasion rate by 32.8 ± 5.7%. A reduction in invasion by a pld2 mutant of A. baumannii had actually been described before [12]. A double knock-out of pld1+2 or pld2+3 reduced the rate even more and the most reduced rate was observed in the Δpld1-3 triple mutant (Fig 8). This finding suggests, that the concerted action of all three PLDs contributes to host cell invasion.

thumbnail
Download:
Fig 8. PLDs affect invasion efficiency of A. baumannii.

A549 cells were infected with A. baumannii wild-type or mutant cells for 5 hours. Then cells were treated with gentamycin to eliminate all external bacteria, prior to release of internal bacteria by lysis with H2O. Internal bacteria were calculated and the cell count of the wild-type was set to 100%.100% corresponds to 112.5 ± 2.5 CFU/ml Asterisks indicate statistical significance by paired t-test; *p<0.05; **p<0.01; ***p<0.001; ****p<0.0001.

https://doi.org/10.1371/journal.pone.0138360.g008

Note that the observed reduction of internal bacteria of the Δpld1-3 mutant could also be caused by a deficiency of this mutant to adhere to the A549 cells. To test this possibility we additionally performed adhesion assays. This assay revealed that the Δpld1-3 mutant is not impaired in cell association to A549 cells compared to the wild-type (105.5 ± 14.57% compared to the wildtype).

Discussion

Here, we report on the identification and characterization of three phospholipases D (PLDs) in A. baumannii ATCC 19606T: PLD1, PLD2 and PLD3. None of the three PLDs was required for growth on PC. However, our finding, that palmitate, a PC cleavage product, liberated from PC via phospholipase A (PLA) activities suggests, that a PLA might be implicated in growth on PC. This is supported by the detection of a highly conserved phospholipase A (HMPREF0010_00411) in the genome of A. baumannii.

The fact that multiple phospholipases D are present throughout the genus Acinetobacter (S2 Table) together with our finding that all PLDs play an important role in virulence and host cell invasion suggests, that PLDs are major virulence factors of A. baumannii. Two of these genes, pld1 and pld2 share a conserved domain architecture (two PLDc_2 domains) and are considerably close paralogs that emerged within contemporary Moraxellaceae and presumably in the last common ancestor of Acinetobacter. The same domain architecture is seen across a broad spectrum of bacterial orthologs suggesting that it represents the archetype of bacterial PLDs [50]. The third PLD, PLD3 differs from this scheme. Although similar in size, the encoded PLD3 protein harbors only a single PLDc_2 domain located in the C-terminal half of the protein with an active site having the functionally relevant aspartic acid residue (D) replaced by a tyrosine (Y). The functional consequence of this D to Y replacement, which is specific for the PLD3s of the genus Acinetobacter (S4 File), remains unclear and is subject to future studies. Notably, a second active site motif is located in the N-terminal half of the protein, although the corresponding PLDc_2 domain appears absent. Our finding is of particular interest since a general feature of all known members of the PLD superfamily is, that two HKD motifs are required for the catalytic activity and that all conserved residues of both HKD motifs are essential for PLD activity [24,45,51]. Even a single mutation of one aspartate residue leads to an almost complete loss of catalytic activity in mammalian PLD1 [52] or even to an insoluble protein such as in the case of Yersinia pestis murine toxin, which represents another member of the PLD superfamily [53]. It will now be interesting to see whether a non-conserved lysine residue detected next to the tyrosine residue, an aspartic acid downstream of the HKD-motif, or the tyrosine itself can take over the function of the otherwise highly conserved aspartic acid residue in the active site. In this instance, it can be speculated that the modified motif might affect functionality thereby modifying the activity of this phospholipase. Alternatively, the hydrolysis could be accomplished by the second active site motif displaying the canonical D. As, however, two fully functional copies of the active site motif are required for formation of the catalytic core, this would require (homo-) dimerization of PLD3. It was shown on the example of the PLD3 ortholog ZUC from Drosophila melanogaster, that this protein indeed is capable of exerting a double functionality. ZUC acts both as a phospholipase hydrolyzing cardiolipin and also as a nuclease, which is involved in the formation of small functional RNAs implicated in the PIWI pathway (piRNAs) [54,55]. From this aspect, it wouldn’t be too surprising if also PLD3 displays an activity beyond the cleavage of phosphatidylcholine. It will be interesting to see whether this is indeed the case, and to what extent PLD3 might even resemble the functional role as a nuclease described for its metazoan orthologs ZUC.

Each of the three PLDs was found to contribute to virulence in the eukaryotic G. mellonella infection model and invasion of human lung epithelial cells. These findings imply that the PLDs of A. baumannii are important effectors mediating virulence and invasion in eukaryotic host cells. pld2 corresponds to the recently identified pld gene in the clinical isolate 98-37-09. Disruption of this pld gene in this isolate via transposon mutagenesis resulted in reduced growth of the isolate in human serum and a decrease in invasion of BEAS-2B bronchial epithelial cells and HeLa cervical cancer cells [12]. Furthermore a reduction in bacterial burden of blood, heart and liver was found. In agreement with the findings by Jacobs et al. [12] our studies revealed that a single pld2 mutation led to a statistically significant reduction in invasion efficiency using A549 lung epithelial cells. However, only upon inactivation of all three pld genes invasion efficiency reached a minimum, indicating that the concerted action of all three PLDs was necessary for full invasion and virulence. From an evolutionary perspective, one observation in this context is particularly interesting. While the sole knockout of pld2 markedly decreases the capability for cell invasion, no effect is seen when knocking out its close paralog pld1. Given the considerably recent evolutionary separation of pld1 and pld2, most likely in the common ancestor of Acinetobacter, this indicates a neo-functionalization of either protein. One indication for such a functional change comes from the domain architecture of both proteins (Fig 2). PLD1 lacks a membrane anchor that is present in PLD2. As a consequence, both proteins are likely to differ in their cellular localization.

Future studies will show, whether the differences observed between the different pld mutants are caused by functional changes such as different substrate specificities and/or different modes of regulation. This would be in line with previous observations, which showed that PLDs not only mediate cleavage of phosphodiester bonds in phospholipids but some also hydrolyze polynucleotide backbones or exhibit phosphatase activity [24,55]. The latter might also contribute to phosphate scavenging of pathogenic bacteria. In this context, a PLD mediated release of phosphate from PC could be beneficial for invasion of the human host via wound infections. It is very interesting to note, that strain ATCC 19606T has been reported to have a reduced ability to metabolize phosphate compounds [56] and it is tempting to speculate, that PLDs might also play a role in phosphate acquisition. Whether this is indeed the case and whether there is even a link between PLDs and the phosphate (pho) regulon found in A. baumannii will be addressed in future studies. Furthermore it will be interesting to see whether the three PLDs possess an active role in the internalization process. The latter has been shown for PLDs of Neisseria gonorrhoeae and P. aeruginosa, which directly interact with host kinases recruited by the pathogens to enter epithelial cells [47,57]

Moreover, PLD activities generate bioactive molecules, such as phosphatidic acid from the hydrolyzation of phospholipids, potentially affecting eukaryotic signal cascades. The latter has already been shown in gonococci, P. aeruginosa and Streptomyces chromofuscus [47,58,59] where the PLD mediated cleavage of lysophosphatidylcholine (LPC) triggers signaling cascades via G-protein-coupled receptors in host cells. The contribution of these distinct features to the three PLDs of A. baumannii will be deciphered in future studies.

Supporting Information

S1 Fig. Pfam domain architecture of the PLD3 protein family.

The tree topology represents the one shown at S3 File. The overview reveals that most of the bacterial PLD3 members and all of the eukaryotic PLD3 orthologs share the presence of only a single PLDc_2 domain.

https://doi.org/10.1371/journal.pone.0138360.s001

(PDF)

S1 Table. TOL taxa set. Comprising 1606 taxa.

https://doi.org/10.1371/journal.pone.0138360.s002

(XLS)

S2 Table. ACI taxa set.

The set of 68 analyzed Acinetobacter taxa and the respective presence/absence of PLD orthologs. Percent identity values are given in parentheses. * denotes that the protein was not present in the protein collection of the data source.

https://doi.org/10.1371/journal.pone.0138360.s003

(XLS)

S3 Table. Primer for markerless mutagenesis.

https://doi.org/10.1371/journal.pone.0138360.s004

(DOCX)

S4 Table. Phyletic pattern of PLD orthologs across the set of 1,608 analyzed taxa.

https://doi.org/10.1371/journal.pone.0138360.s005

(DOCX)

S5 Table. Taxa set for maximum likelihood tree reconstruction of PLD1 and PLD2 orthologs.

https://doi.org/10.1371/journal.pone.0138360.s006

(DOCX)

S1 File. Orthologous and co-orthologous sequences with given protein ID, taxon name, assignment and amino acid sequence.

https://doi.org/10.1371/journal.pone.0138360.s007

(XLSX)

S2 File. Tree topology for 1180 taxa representing the tree for all three PLD orthologs.

https://doi.org/10.1371/journal.pone.0138360.s008

(NEXUS)

S3 File. Tree topology for 246 taxa representing the tree of PLD3 orthologs.

https://doi.org/10.1371/journal.pone.0138360.s009

(NEXUS)

S4 File. Alignment of HKD motif.

Alignment of HKD motif for a selection of 238 bacterial and archaeaen taxa including 68 Acinetobacter strains with the unusual HKD-like motif with given taxon name and amino acid sequence.

https://doi.org/10.1371/journal.pone.0138360.s010

(XLSX)

Acknowledgments

We thank Valerie De-Crecy Legard for providing A. baylyi ADP1 hptA::sacB/kanR for amplification of the sacB/kanR-cassette and Sara Christ for support during the work with the cell cultures. We thank Volker Müller for critical reading of the manuscript.

Author Contributions

Conceived and designed the experiments: BA JS. Performed the experiments: JS HB. Analyzed the data: BA JS HB IE SG. Contributed reagents/materials/analysis tools: BA IE SG. Wrote the paper: BA IE JS.

References

  1. 1. Dijkshoorn L, Nemec A, Seifert H. An increasing threat in hospitals: multidrug-resistant Acinetobacter baumannii. Nat Rev Microbiol. 2007;5: 939–951. pmid:18007677
  2. 2. Antunes LC, Visca P, Towner KJ. Acinetobacter baumannii: evolution of a global pathogen. Pathog Dis. 2014;71: 292–301. pmid:24376225
  3. 3. Kempf M, Rolain JM. Emergence of resistance to carbapenems in Acinetobacter baumannii in Europe: clinical impact and therapeutic options. Int J Antimicrob Agents. 2012;39: 105–114. pmid:22113193
  4. 4. Peleg AY, Seifert H, Paterson DL. Acinetobacter baumannii: emergence of a successful pathogen. Clin Microbiol Rev. 2008;21: 538–582. pmid:18625687
  5. 5. Gaynes R, Edwards JR, System NNIS. Overview of nosocomial infections caused by gram-negative bacilli. Clin Infect Dis. 2005;41: 848–854. pmid:16107985
  6. 6. Falagas ME, Bliziotis IA, Siempos II. Attributable mortality of Acinetobacter baumannii infections in critically ill patients: a systematic review of matched cohort and case-control studies. Crit Care. 2006;10: R48. pmid:16563184
  7. 7. Göttig S, Gruber TM, Higgins PG, Wachsmuth M, Seifert H, Kempf VAJ. Detection of pan drug-resistant Acinetobacter baumannii in Germany. J Antimicrob Chemother. 2014;69: 2578–2579. pmid:24833751
  8. 8. Jawad A, Seifert H, Snelling AM, Heritage J, Hawkey PM. Survival of Acinetobacter baumannii on dry surfaces: comparison of outbreak and sporadic isolates. J Clin Microbiol. 1998;36: 1938–1941. pmid:9650940
  9. 9. Clemmer KM, Bonomo RA, Rather PN. Genetic analysis of surface motility in Acinetobacter baumannii. Microbiology. 2011;157: 2534–2544. pmid:21700662
  10. 10. Dorsey CW, Beglin MS, Actis LA. Detection and analysis of iron uptake components expressed by Acinetobacter baumannii clinical isolates. J Clin Microbiol. 2003;41: 4188–4193. pmid:12958246
  11. 11. Camarena L, Bruno V, Euskirchen G, Poggio S, Snyder M. Molecular mechanisms of ethanol-induced pathogenesis revealed by RNA-sequencing. PLoS Pathog. 2010;6: e1000834. pmid:20368969
  12. 12. Jacobs AC, Hood I, Boyd KL, Olson PD, Morrison JM, Carson J, et al. Inactivation of phospholipase D diminishes Acinetobacter baumannii pathogenesis. Infect Immun. 2010;78: 1952–1962. pmid:20194595
  13. 13. Gaddy JA, Arivett BA, McConnell MJ, Lopez-Rojas R, Pachon J, Actis LA. Role of acinetobactin-mediated iron acquisition functions in the interaction of Acinetobacter baumannii strain ATCC 19606T with human lung epithelial cells, Galleria mellonella caterpillars, and mice. Infect Immun. 2012;80: 1015–1024. pmid:22232188
  14. 14. Tomaras AP, Dorsey CW, Edelmann RE, Actis LA. Attachment to and biofilm formation on abiotic surfaces by Acinetobacter baumannii: involvement of a novel chaperone-usher pili assembly system. Microbiology. 2003;149: 3473–3484. pmid:14663080
  15. 15. Lee JC, Koerten H, van den Broek P, Beekhuizen H, Wolterbeek R, van den Barselaar M, et al. Adherence of Acinetobacter baumannii strains to human bronchial epithelial cells. Res Microbiol. 2006;157: 360–366. pmid:16326077
  16. 16. Choi CH, Lee EY, Lee YC, Park TI, Kim HJ, Hyun SH, et al. Outer membrane protein 38 of Acinetobacter baumannii localizes to the mitochondria and induces apoptosis of epithelial cells. Cell Microbiol. 2005;7: 1127–1138. pmid:16008580
  17. 17. Antunes LC, Imperi F, Carattoli A, Visca P. Deciphering the multifactorial nature of Acinetobacter baumannii pathogenicity. PLoS One. 2011;6: e22674. pmid:21829642
  18. 18. Brown SA, Palmer KL, Whiteley M. Revisiting the host as a growth medium. Nat Rev Microbiol. 2008;6: 657–666. pmid:18679171
  19. 19. Vance DE. Role of phosphatidylcholine biosynthesis in the regulation of lipoprotein homeostasis. Curr Opin Lipidol. 2008;19: 229–234. pmid:18460912
  20. 20. Bernhard W, Hoffmann S, Dombrowsky H, Rau GA, Kamlage A, Kappler M, et al. Phosphatidylcholine molecular species in lung surfactant: composition in relation to respiratory rate and lung development. Am J Respir Cell Mol Biol. 2001;25: 725–731. pmid:11726398
  21. 21. Girod S, Galabert C, Lecuire A, Zahm JM, Puchelle E. Phospholipid composition and surface-active properties of tracheobronchial secretions from patients with cystic fibrosis and chronic obstructive pulmonary diseases. Pediatr Pulmonol. 1992;13: 22–27. pmid:1589308
  22. 22. Sun Z, Kang Y, Norris MH, Troyer RM, Son MS, Schweizer HP, et al. Blocking phosphatidylcholine utilization in Pseudomonas aeruginosa, via mutagenesis of fatty acid, glycerol and choline degradation pathways, confirms the importance of this nutrient source in vivo. PLoS One. 2014;9: e103778. pmid:25068317
  23. 23. McConnell MJ, Actis L, Pachon J. Acinetobacter baumannii: human infections, factors contributing to pathogenesis and animal models. FEMS Microbiol Rev. 2013;37: 130–155. pmid:22568581
  24. 24. Selvy PE, Lavieri RR, Lindsley CW, Brown HA. Phospholipase D: enzymology, functionality, and chemical modulation. Chem Rev. 2011;111: 6064–6119. pmid:21936578
  25. 25. Songer JG. Bacterial phospholipases and their role in virulence. Trends Microbiol. 1997;5: 156–161. pmid:9141190
  26. 26. Titball RW. Bacterial phospholipases. Symp Ser Soc Appl Microbiol. 1998;27: 127S–137S. pmid:9750370
  27. 27. Sand M, de Berardinis V, Mingote A, Santos H, Göttig S, Müller V, et al. Salt adaptation in Acinetobacter baylyi: identification and characterization of a secondary glycine betaine transporter. Arch Microbiol. 2011;193: 723–730. pmid:21567174
  28. 28. Berendes F, Sabarth N, Averhoff B, Gottschalk G. Construction and use of an ipb DNA module to generate Pseudomonas strains with constitutive trichloroethene and isopropylbenzene oxidation activity. Appl Environ Microbiol. 1998;64: 2454–2462. pmid:9647815
  29. 29. Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, et al. Pfam: the protein families database. Nucleic Acids Res. 2014;42: D222–D230. pmid:24288371
  30. 30. Eddy SR. Accelerated Profile HMM Searches. PLoS Comput Biol. 2011;7: e1002195. pmid:22039361
  31. 31. Moore AD, Held A, Terrapon N, Weiner J 3rd, Bornberg-Bauer E. DoMosaics: software for domain arrangement visualization and domain-centric analysis of proteins. Bioinformatics. 2014;30: 282–283. pmid:24222210
  32. 32. Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8: 785–786. pmid:21959131
  33. 33. Ebersberger I, Strauss S, von Haeseler A. HaMStR: profile hidden markov model based search for orthologs in ESTs. BMC Evol Biol. 2009;9: 157. pmid:19586527
  34. 34. Katoh K, Toh H. Recent developments in the MAFFT multiple sequence alignment program. Brief Bioinform. 2008;9: 286–298. pmid:18372315
  35. 35. Darriba D, Taboada GL, Doallo R, Posada D. ProtTest 3: fast selection of best-fit models of protein evolution. Bioinformatics. 2011;27: 1164–1165. pmid:21335321
  36. 36. Müller T, Vingron M. Modeling amino acid replacement. J Comput Biol. 2000;7: 761–776. pmid:11382360
  37. 37. Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30: 1312–1313. pmid:24451623
  38. 38. Price MN, Dehal PS, Arkin AP. FastTree 2—approximately maximum-likelihood trees for large alignments. PLoS One. 2010;5: e9490. pmid:20224823
  39. 39. Shimodaira H, Hasegawa M. Multiple comparisons of log-likelihoods with applications to phylogenetic inference. Molecular Biology and Evolution. 1999;16: 1114–1116.
  40. 40. Metzgar D, Bacher JM, Pezo V, Reader J, Doring V, Schimmel P, et al. Acinetobacter sp. ADP1: an ideal model organism for genetic analysis and genome engineering. Nucleic Acids Res. 2004;32: 5780–5790. pmid:15514111
  41. 41. Carruthers MD, Nicholson PA, Tracy EN, Munson RS Jr. Acinetobacter baumannii utilizes a type VI secretion system for bacterial competition. PLoS One. 2013;8: e59388. pmid:23527179
  42. 42. Carruthers MD, Harding CM, Baker BD, Bonomo RA, Hujer KM, Rather PN, et al. Draft Genome Sequence of the Clinical Isolate Acinetobacter nosocomialis Strain M2. Genome Announc. 2013;1.
  43. 43. Choi CH, Lee JS, Lee YC, Park TI, Lee JC. Acinetobacter baumannii invades epithelial cells and outer membrane protein A mediates interactions with epithelial cells. BMC Microbiol. 2008;8: 216. pmid:19068136
  44. 44. Sand M, Stahl J, Waclawska I, Ziegler C, Averhoff B. Identification of an osmo-dependent and an osmo-independent choline transporter in Acinetobacter baylyi: implications in osmostress protection and metabolic adaptation. Environ Microbiol. 2014;16: 1490–1502. pmid:23889709
  45. 45. Ponting CP, Kerr ID. A novel family of phospholipase D homologues that includes phospholipid synthases and putative endonucleases: identification of duplicated repeats and potential active site residues. Protein Sci. 1996;5: 914–922. pmid:8732763
  46. 46. Vallenet D, Nordmann P, Barbe V, Poirel L, Mangenot S, Bataille E, et al. Comparative analysis of Acinetobacters: three genomes for three lifestyles. PLoS One. 2008;3: e1805. pmid:18350144
  47. 47. Jiang F, Waterfield NR, Yang J, Yang G, Jin Q. A Pseudomonas aeruginosa type VI secretion phospholipase D effector targets both prokaryotic and eukaryotic cells. Cell Host Microbe. 2014;15: 600–610. pmid:24832454
  48. 48. Boeris PS, Lucchesi GI. The phosphatidylcholine synthase of Pseudomonas putida ATCC 12633 is responsible for the synthesis of phosphatidylcholine, which acts as a temporary reservoir for Al3+. Microbiology. 2012;158: 1249–1257. pmid:22343357
  49. 49. Peleg AY, Jara S, Monga D, Eliopoulos GM, Moellering RC Jr, Mylonakis E. Galleria mellonella as a model system to study Acinetobacter baumannii pathogenesis and therapeutics. Antimicrob Agents Chemother. 2009;53: 2605–2609. pmid:19332683
  50. 50. Stuckey JA, Dixon JE. Crystal structure of a phospholipase D family member. Nat Struct Biol. 1999;6: 278–284. pmid:10074947
  51. 51. Mahankali M, Alter G, Gomez-Cambronero J. Mechanism of enzymatic reaction and protein-protein interactions of PLD from a 3D structural model. Cell Signal. 2014;27: 69–81. pmid:25308783
  52. 52. Sung TC, Roper RL, Zhang Y, Rudge SA, Temel R, Hammond SM, et al. Mutagenesis of phospholipase D defines a superfamily including a trans-Golgi viral protein required for poxvirus pathogenicity. EMBO J. 1997;16: 4519–4530. pmid:9303296
  53. 53. Rudolph AE, Stuckey JA, Zhao Y, Matthews HR, Patton WA, Moss J, et al. Expression, characterization, and mutagenesis of the Yersinia pestis murine toxin, a phospholipase D superfamily member. J Biol Chem. 1999;274: 11824–11831. pmid:10207000
  54. 54. Ross RJ, Weiner MM, Lin H. PIWI proteins and PIWI-interacting RNAs in the soma. Nature. 2014;505: 353–359. pmid:24429634
  55. 55. Nishimasu H, Ishizu H, Saito K, Fukuhara S, Kamatani MK, Bonnefond L, et al. Structure and function of Zucchini endoribonuclease in piRNA biogenesis. Nature. 2012;491: 284–287. pmid:23064230
  56. 56. Peleg AY, de Breij A, Adams MD, Cerqueira GM, Mocali S, Galardini M, et al. The success of Acinetobacter species; genetic, metabolic and virulence attributes. PLoS One. 2012;7: e46984. pmid:23144699
  57. 57. Edwards JL, Apicella MA. Neisseria gonorrhoeae PLD directly interacts with Akt kinase upon infection of primary, human, cervical epithelial cells. Cell Microbiol. 2006;8: 1253–1271. pmid:16882030
  58. 58. Edwards JL, Entz DD, Apicella MA. Gonococcal phospholipase D modulates the expression and function of complement receptor 3 in primary cervical epithelial cells. Infect Immun. 2003;71: 6381–6391. pmid:14573659
  59. 59. van Dijk MC, Postma F, Hilkmann H, Jalink K, van Blitterswijk WJ, Moolenaar WH, et al. Exogenous phospholipase D generates lysophosphatidic acid and activates Ras, Rho and Ca2+ signaling pathways. Curr Biol. 1998;8: 386–392. pmid:9545198