Systems level analysis of two-component signal transduction systems in Erwinia amylovora: Role in virulence, regulation of amylovoran biosynthesis and swarming motility

  • Youfu Zhao1Email author,

    Affiliated with

    • Dongping Wang1,

      Affiliated with

      • Sridevi Nakka1,

        Affiliated with

        • George W Sundin2 and

          Affiliated with

          • Schuyler S Korban3

            Affiliated with

            BMC Genomics200910:245

            DOI: 10.1186/1471-2164-10-245

            Received: 30 December 2008

            Accepted: 26 May 2009

            Published: 26 May 2009

            Abstract

            Background

            Two-component signal transduction systems (TCSTs), consisting of a histidine kinase (HK) and a response regulator (RR), represent a major paradigm for signal transduction in prokaryotes. TCSTs play critical roles in sensing and responding to environmental conditions, and in bacterial pathogenesis. Most TCSTs in Erwinia amylovora have either not been identified or have not yet been studied.

            Results

            We used a systems approach to identify TCST and related signal transduction genes in the genome of E. amylovora. Comparative genomic analysis of TCSTs indicated that E. amylovora TCSTs were closely related to those of Erwinia tasmaniensis, a saprophytic enterobacterium isolated from apple flowers, and to other enterobacteria. Forty-six TCST genes in E. amylovora including 17 sensor kinases, three hybrid kinases, 20 DNA- or ligand-binding RRs, four RRs with enzymatic output domain (EAL-GGDEF proteins), and two kinases were characterized in this study. A systematic TCST gene-knockout experiment was conducted, generating a total of 59 single-, double-, and triple-mutants. Virulence assays revealed that five of these mutants were non-pathogenic on immature pear fruits. Results from phenotypic characterization and gene expression experiments indicated that several groups of TCST systems in E. amylovora control amylovoran biosynthesis, one of two major virulence factors in E. amylovora. Both negative and positive regulators of amylovoran biosynthesis were identified, indicating a complex network may control this important feature of pathogenesis. Positive (non-motile, EnvZ/OmpR), negative (hypermotile, GrrS/GrrA), and intermediate regulators for swarming motility in E. amylovora were also identified.

            Conclusion

            Our results demonstrated that TCSTs in E. amylovora played major roles in virulence on immature pear fruit and in regulating amylovoran biosynthesis and swarming motility. This suggested presence of regulatory networks governing expression of critical virulence genes in E. amylovora.

            Background

            Prokaryotes use their small size and metabolic diversity to dominate every conceivable niche on earth. A large part of this success comes from the evolution of elaborate sensory systems to monitor and respond to dramatic fluctuations in their environment. Primary means of signal transduction in bacteria involve two-component signal transduction systems (TCSTs) [1, 2]. The term "two-component" was coined in 1986 to describe a new class of regulatory systems found in bacteria; today, two-component systems represent major paradigms for signal transduction in prokaryotes and in lower eukaryotes [3, 4]. The prototypical two-component system consists of a histidine kinase protein (HK) containing a conserved kinase core and a response regulator protein (RR) containing a conserved regulatory domain [4, 5]. Extracellular stimuli are sensed by, and serve to modulate HK activities. The HK is autophosphorylated at a histidine residue, creating a high-energy phosphoryl group that is subsequently transferred to an aspartate residue in the RR via a reaction catalyzed by the RR itself. Phosphotransfer to the RR leads to activation of a downstream effector domain that elicits a specific response [1, 4].

            Two-component systems are distributed at varying frequencies among organisms of all domains, including Eubacteria, Archaea, and Eukarya [1, 3]. The availability of complete genome sequences has allowed for a definitive assessment of the prevalence of two-component proteins. There are 30 HKs (five of which are hybrid kinases) and 32 RRs in Escherichia coli [6]. However, the number of two-component proteins differs greatly among bacteria. Often, parasitic bacteria encode fewer signaling proteins than relatively free-living bacteria [7]. Among sequenced plant pathogenic bacteria, the number of TCST genes is quite different. For example, the xylem-limited Xylella fastidiosa has the fewest TCST genes; whereas, the relatively adaptable Pseudomonas syringae and Xanthomonas spp. have the largest number of TCST genes [8, 9]. Perhaps, the most attractive reason for studying two-component systems in bacteria is that TCST systems are used by pathogenic bacteria to control expression of virulence factors required for infection. Many such TCST systems are identified in both human and animal pathogens and in plant pathogens [7, 10]. The Agrobacterium tumefaciens VirA/VirG system, the GacA/GacS of both Pseudomonas sp. and Pectobacterium carotovora, and the RpfCG of Xanthomonas spp. and X. fastidiosa are probably the most well-known and studied TCST systems involved in virulence gene expression in plant pathogens [1012]. However, most of these studies focused on one or a few TCSTs. Only recently has genome-wide analysis of TCSTs in plant pathogenic bacteria become possible due to increased availability of whole genome sequences [8, 9, 13].

            E. amylovora is the causal agent of fire blight, a devastating necrotic disease affecting apple, pear, and other rosaceous plants. Fire blight is one of the most important bacterial plant diseases worldwide that has a significant economic impact, resulting in crop losses of millions of dollars per year [14]. As a member of the Enterobacteriaceae, E. amylovora is related to many important human and animal pathogens such as Es. coli, Salmonella enterica, Shigella flexneri, Yersinia enterocolitica, and Y. pestis. Like many other Gram-negative plant pathogenic bacteria, E. amylovora utilizes both a type III secretion (T3SS) apparatus to deliver effector proteins into host plant cells, and the extracellular polysaccharide (EPS) amylovoran to cause disease [1517]. In E. amylovora, structural components of T3SS encoded by the Hrp regulon are regulated by the two-component system HrpX and HrpY, which direct the expression of the NtrC family σ54-dependent, enhancer-binding protein HrpS [18]. Both HrpY and HrpS function in the activation of expression of the alternate sigma factor HrpL, thereby regulating various genes and operons of the Hrp regulon. Expression of hrpX and hrpS is regulated by low pH, low nutrients, and low temperature conditions, mimicking the plant apoplast [18]. The biosynthesis of amylovoran is regulated by another TCST system, the RcsCDB phosphorelay system [17]. Recently, we identified several TCST genes including hrpX, grrS, and envZ that are induced during infection of host tissue in E. amylovora, indicating that TCSTs are key players in controlling the expression of virulence factors required for infection, and thus we hypothesize that there are likely networks controlling virulence gene expression [19]. However, other TCSTs in E. amylovora have either not been identified or have not yet been studied. This prompted us to use a systems approach to explore the function of TCSTs in E. amylovora. Understanding the genetics and molecular mechanisms of E. amylovora signaling will greatly enhance the likelihood of developing novel methods of controlling the disease.

            The EnvZ-OmpR and GrrS-GrrA (also called GacSA, BarA-UvrY) are two widely-distributed and well-studied TCSTs in γ-proteobacteria. They represent those paradigms of signal transduction systems having pleiotropic effects, thus suggesting both systems are global regulators [2, 20]. EnvZ is a transmembrane sensor that predominantly responds to acidic pH conditions and changes in osmolarity, and subsequently phosphorylates OmpR. The EnvZ-OmpR system has been originally reported to govern the expression of ompC and ompF genes, encoding two major outer-membrane porins [21]. In addition to its role in porin osmoregulation, OmpR is involved in regulating various cellular components including flagellar gene expression, fatty acid transport, curli fibre formation, and cell division as a dual regulator; i.e., negative or positive [22]. In Salmonella spp., OmpR-EnvZ regulates another TCST system SsrA-SsrB, that in turn regulates the T3SS produced by Salmonella pathogenicity island 2 (SPI-2) [2326]. OmpR also negatively regulates expression of invasin in Yersinia enterocolitica and the T3SS in P. syringae [27, 28].

            GacS and GacA homologs have been identified in many Gram-negative bacterial genera, including Azotobacter, Erwinia, Escherichia, Legionella, Pectobacterium, Pseudomonas, Salmonella, Serratia, and Vibrio [12, 20, 29]. The sensor kinase GacS was initially discovered as a key regulator of virulence in the plant pathogen P. syringae [see review [12]]. The GacSA system has since been reported to regulate an array of phenotypes, including biofilm formation, alginate biosynthesis, production of toxins and extracellular enzymes [30, 31], proteases, siderophores, swarming motility [32], and type III secretion [10, 29, 3335]. Similar to EnvZ-OmpR, GacSA also has a dual regulatory function; i.e., as a positive or negative regulator. The GacS/GacA (BarA/UvrY) system positively controls the expression of one to five genes specifying small RNAs (sRNAs), thus upregulating the productions of proteins that are otherwise repressed by small RNAs (RsmA/CsrA) [32, 36]. Most gacS/gacA mutants demonstrate reduced production of virulence factors and reduced virulence in a variety of host-pathogen systems [12, 20]. An alternative scenario is that the GacSA cascade downregulates the expression of flagellar genes in P. fluorescens or Es. coli [20].

            In this study, our genome-wide analysis, using the recently closed genome sequence of E. amylovora Ea273 http://​www.​sanger.​ac.​uk/​Projects/​E_​amylovora/​ identified a total of 46 TCST genes in E. amylovora including 17 sensor kinases, 20 response regulators, three hybrid HKs, four RRs with enzymatic output domain (EAL-GGDEF proteins), and two kinases. A comparative genomic analysis of TCSTs was then conducted in related enterobacteria enabling classification of the TCSTs in E. amylovora. A total of 59 deletion mutants were generated, and their contribution to virulence, amylovoran biosynthesis, and swarming motility was characterized. Our findings suggested that TCSTs in E. amylovora played a major role in virulence and in the regulation of amylovoran biosynthesis and swarming motility. This indicated that networks of gene regulation existed in E. amylovora that respond to different environmental and host signals.

            Results

            Identification and comparison of TCSTs in E. amylovora

            We utilized two approaches to identify HKs and RRs in the genome of E. amylovora. First, using a candidate gene approach, known HK and RR sequences from Es. coli K12 and P. carotovora pv. atroseptica SCRI1043 were used to search the closed genome sequence of E. amylovora http://​www.​sanger.​ac.​uk/​Projects/​E_​amylovora/​. Second, based on conserved domains of known HKs and RRs, other putative HKs and RRs in E. amylovora were identified and confirmed by searching the complete genome sequences for proteins containing HK and RR domains using Pfam Hidden Markov Model (HMM) profiles and BLASTP. As a result, we identified 46 TCST and other related signal transduction genes in E. amylovora (Table 1 and Table S1, [see Additional file 1]), excluding two sets of chemotaxis genes (cheABRWYZ). These putative regulatory genes represented 1.36% of the genome (3367 genes) [37]. Among them, 17 are sensor kinases, three hybrid HKs, 20 DNA-or ligand-binding RRs, four RRs with enzymatic output domain (EAL-GGDEF proteins), and two putative kinases (Table 1).
            Table 1

            TCSTs and other signal transduction genes in Erwinia amylovora and mutant construction

            Gene/Operona

            HK/RR/Hybrid/kinaseb

            Mutants constructedc

            arcB, arcA

            HK, RR

            ΔarcB, ΔarcA

            baeSR

            HK, RR

            ΔbaeS, ΔbaeR, ΔbaeSR

            cpxA1R1

            HK, RR

            ΔcpxA1, ΔcpxR1, ΔcpxAR1

            cpxA2R2

            HK, RR

            ΔcpxA2, ΔcpxR2, ΔcpxAR2

            dcuSR

            HK, RR

            ΔdcuS, ΔdcuR, ΔdcuSR

            envZ/ompR

            HK, RR

            ΔenvZ, ΔompR, ΔenvZ/ompR

            grrS, grrA

            HK hybrid, RR

            ΔgrrS, ΔgrrA

            hrpXY, hrpS

            HK, RR, EBP

            ΔhrpX, ΔhrpY, ΔhrpXY, ΔhrpXYS

            kdpD, kdpE

            HK, RR

            ΔkdpD, ΔkdpE

            narQP

            HK, RR

            ΔnarQ, ΔnarP, ΔnarQP

            phoQP

            HK, RR

            ΔphoQ, ΔphoP, ΔphoQP

            phoRB

            HK, RR

            ΔphoR, ΔphoB, ΔphoRB

            pmrBA

            HK hybrid, RR

            ΔpmrB, ΔpmrA, ΔpmrBA

            rcsCDB

            HK hybrid, HPT, RR

            ΔrcsD, ΔrcsB, ΔrcsC, ΔrcsBD

            rstB, rstA

            HK, RR

            ΔrstB, ΔrstA

            yehUT

            HK, RR

            ΔyehU, ΔyehT, ΔyehUT

            yfhKA

            HK, RR

            ΔyfhK, ΔyfhA, ΔyfhKA

            Other signal transduction genes used in this study b

            eamIR

            AHL, RR

            ΔeamI, ΔeamR

            luxQ, luxP

            AI2 sensor kinase, RR

            ΔluxQ, ΔluxP

            Spk1

            Serine kinase

            Δspk1

            ybjN

            Kinase

            ΔybjN

            yciR

            GGDEF-EAL

            ΔyciR

            yddV

            GGDEF

            ΔyddV

            yegE

            GGDEF

            ΔyegE

            yoaD

            EAL

            ΔyoaD

            a Genes for one TCST system are listed separately if they are not encoded within an operon in the genome. Otherwise, TCSTs are listed as an operon if genes are present together. Mutants for ntrBC were not constructed.

            b All TCSTs are sensor kinases and their corresponding DNA-binding response regulators. We included two sets of quorum sensing genes and four RRs with enzymatic output domain (EAL-GGDEF proteins) and two kinases. EBP: NtrC family of enhancing-binding protein. HPT: Histidine phosphotransfer domain; AHL: Acyl-homoserine lactone; GGDEF: Diguanylate cyclase; EAL: Type I c-di-GMP phosphodiesterase.

            c Mutant designation see Table 3.

            When compared to TCSTs in other enterobacteria, the majority of TCSTs in E. amylovora have counterparts in related enterobacteria, including Dickeya dadantii (Erwinia chrysanthemi) 3937, P. carotovora subsp. atroseptica SCRI1043, Erwinia tasmaniensis Et1/99, and Es. coli K12 (Table 2 and Table S1, [see Additional file 1]). Only one HK and RR pair (ypdAB) in E. tasmaniensis is not present in E. amylovora Ea273 (Table 2). The number of TCSTs in E. amylovora is the lowest among other sequenced related plant pathogenic enterobacteria and Es. coli (Table S1, [see Additional file 1]). BLAST searches revealed that the HKs and RRs in E. amylovora shared the highest amino acid (aa) identity/similarity to those of E. tasmaniensis Et1/99, an epiphytic bacterium isolated from apple flowers in Australia [38], except for HrpX/Y which shared the highest aa identity/similarity to those in Erwinia pyrifoliae, a related blight pathogen of Asian pear (Table 2). It is interesting to note that HrpXY, which regulates T3SS gene expression in plant pathogenic enterobacteria, is the only TCST that is not present in other mammalian enterobacterial pathogens [39].
            Table 2

            Comparison of E. amylovora and E. tasmaniensis TCST and other signal transduction genesa

            Geneb

            Length of deduced Amino acidc

            aa Identity/Similarity (%)d

            Geneb

            Length of deduced Amino acidc

            aa Identity/Similarity (%)d

            arcB

            779

            93/95

            arcA

            238

            96/99

            baeS

            461

            89/93

            baeR

            235

            88/94

            grrS (barA)

            909

            91/94

            grrA (Eta_20780)

            219

            94/98

            cpxA1

            459

            96/98

            cpxR1

            233

            95/97

            cpxA2 (Eta_10960)

            450

            82/88

            cpxR2 (Eta_10950)

            226

            88/93

            dcuS

            546

            79/86

            dcuR

            239

            86/93

            envZ

            449

            95/97

            ompR

            239

            100/100

            hrpX

            494

            89/93

            hrpY

            213

            94/96

            kdpD

            890

            95/97

            kdpE (Eta_23200)

            226

            89/92

            luxQ (Eta_11320)

            419 (420)

            89/93

            luxP (Eta_11350)

            258

            74/89

            narQ

            330

            N/A

            narP

            209 (210)

            71/83

            ntrB

            349

            96/98

            ntrC

            469

            98/98

            phoQ

            481

            92/96

            phoP

            222

            95/98

            phoR

            437

            88/92

            phoB

            229

            96/98

            pmrB

            350

            81/89

            pmrA

            219

            94/98

            rcsC

            885 (890)

            88/93

            rcsB

            215

            99/100

            rcsD (yojN)

            987 (951)

            84/91

            eamR (sdiA)

            240

            84/92

            rstB

            426

            92/96

            rstA

            242

            93/97

            ybjN (Eta_21780)

            159

            91/96

            Spk1 (Eta_06360)

            468

            91/92

            yddV (hmsT)

            365

            81/88

            yegE

            848 (877)

            78/88

            yehU

            564 (562)

            84/90

            yehT

            239

            80/89

            yfhK

            463 (475)

            87/91

            yfhA (Eta_10100)

            444

            93/96

            yoaD (rtn)

            516 (523)

            81/89

            yciR (Eta_30620)

            818 (820)

            82/91

            a Sequences for E. amylovora TCST genes are from the whole genomic sequence of E. amylovorahttp://​www.​sanger.​ac.​uk/​Projects/​E_​amylovora/​. One TCTS system (ypdAB) in E. tasmaniensis is not present in E. amylovora Ea273. Two sets of CheABRWYZ exist in E. amylovora, but only one set in E. tasmaniensis.

            b In parentheses is gene or gene number for the corresponding gene in E. tasmaniensis.

            c The number in parentheses is the length of deduced amino acid for the corresponding protein in E. tasmaniensis which differs from those in E. amylovora.

            d All data are the highest aa similarity and identify found in NCBI database as of today except for HrpXY which share the highest aa identity/similarity to those in E. pyrifoliae (98/99%; 99/99%, respectively). N/A: not available.

            Classification of HKs and RRs in E. amylovora

            In general, HKs are highly variable in amino acid sequence length; whereas, RRs are highly conserved in their receiver domain (REC). Classification of HKs is usually based on alignment of residues surrounding the H-box that contains the conserved His residue within the histidine kinase domain (HisKA) [9, 40]. Multiple sequence alignment using Clustal X has classified the HKs of E. amylovora into four groups [40]. Among them, 15 belong to groups IA (11), IB (3), and IC (1); and one (DcuS), two (NarQ and HrpX), and one (YehU) belong to groups II, III, and IV HKs, respectively (Table S1, [see Additional file 1] and Figure S1, [see Additional file 2]). This is the first time that HrpX has been classified as a type III HK. However, RcsD (also known as YojN), containing a histidine phosphotransfer (HPt) domain, but lacking a conserved histidine residue in the HisKA domain, could not be classified.

            Analysis of the domain architecture of E. amylovora HKs, by either searching the Conserved Domain Databases (CDD) [41] or using the SMART program [42]http://​smart.​embl-heidelberg.​de/​, has revealed that all HKs in E. amylovora, except for HrpX and NtrB, contain one to three hydrophobic transmembrane (TM) helices within the N-terminal sensor region. This indicates that these HKs are periplasmic sensors. HrpX and NtrB contain two and one PAS domain (initially found in PER, ARNT, and SIM proteins), respectively, within the N-terminal sensor region, indicating that they are soluble and cytoplasmic proteins that sense intracellular signals [18]. The YehU contains a 5TMR-LYT (5 TM receptor domain, LytS-like) sensory input domain before the two TM helices, thus indicating that the YehU is an intramembrane sensor.

            Often, RRs are classified into families based on their output domains and domain combinations [43]. By searching the SMART program and conducting a structural similarity of C-terminal output domains [4, 43], it is revealed that 20 RRs contain either DNA- or ligand-binding output domains in E. amylovora (Table S1, [see Additional file 1]). Among them, eleven, four, two, and one are further classified into OmpR-, NarL-, NtrC-, and LytR-like proteins, respectively (Table S1, [see Additional file 1]). These RRs belong to four major families of RRs found in prokaryotes and account for about 64% of all RRs identified thus far [43]. In addition, four RRs (YciR, YddV, YegE, and YoaD) contain enzymatic output domains including GGDEF and EAL domains (Table 1). The GGDEF and EAL domains have diguanylate cyclase and Type I c-di-GMP phosphodiesterase activity, respectively, and these proteins are involved in regulating virulence functions in some plant pathogens [44].

            Systematic deletion of TCSTs and related signaling genes, and their roles in virulence on immature pear fruit

            To pursue comprehensive identification of two-component signaling pathways required for virulence, swarming motility, and amylovoran biosynthesis in E. amylovora, deletion mutants were generated for those genes/operons identified above, and also included hrpS and eamI. Deletions were made using the Lambda-Red recombinase technique, as previously described for Es. coli and E. amylovora [4547]. A total of 59 single-, double- and triple- stable deletion mutants were generated (Table 1 and Table 3). For these mutants, the majority of the coding region of each gene or operon was deleted and replaced by the CmR or KmR marker gene, except for the first and last 50 nt of the gene or operon, thus minimizing any polar effects of the mutation.
            Table 3

            E. amylovora TCTS mutants and their phenotypes

            Gene/Operona

            Mutant designation

            Virulence assayb

            Amylovoran production (48 h) (OD600)c

            Significance levele

            Distance of swarming motility (48 h, cm)d

            Significance levele

            WT

             

            +++

            0.08 ± 0.002

            nop

            C

            2.7 ± 0.15

            cde

            ams

            Z0118Δams

            ---

            0.002 ± 0.0005

            r

            ND

            ND

            ND

            flhD

            Z2946ΔflhD

            +++

            ND

            ND

            C

            0.7 ± 0

            w

            fliA

            Z2986ΔfliA

            +++

            ND

            ND

            C

            0.7 ± 0

            w

            arcA

            Z1306ΔarcA

            +++

            0.053 ± 0.0001

            qr

            C

            1.23 ± 0.09

            uv

            arcB

            Z0117ΔarcB

            +++

            0.067 ± 0.0002

            pq

            C

            1.83 ± 0.28

            qr

            baeR

            Z3120ΔbaeR

            +++

            0.14 ± 0.045

            klmn

            I

            1.87 ± 0.25

            qr

            baeS

            Z3119ΔbaeS

            +++

            0.13 ± 0.03

            klmn

            I

            2.07 ± 0.06

            nopq

            baeSR

            Z3219–20ΔbaeSR

            +++

            0.1 ± 0.01

            lmn

            I

            2.07 ± 0.06

            nopq

            cpxA1

            Z0686ΔcpxA1

            +++

            0.09 ± 0.014

            nop

            I

            1.40 ± 0.22

            tu

            cpxA1R1

            Z0686–87ΔcpxAR1

            +++

            0.055 ± 0.003

            pqr

            I

            1.33 ± 0.05

            uv

            cpxR1

            Z0687ΔcpxR1

            +++

            0.08 ± 0.008

            nop

            I

            1.16 ± 0.13

            v

            cpxA2

            Z3367ΔcpxA2

            +++

            0.35 ± 0.023

            hgi

            I

            2.5 ± 0.4

            efghi

            cpxA2R2

            Z3367–68ΔcpxAR2

            +++

            0.29 ± 0.07

            i

            I

            2.07 ± 0.34

            opqr

            cpxR2

            Z3368ΔcpxR2

            +++

            0.35 ± 0.058

            hgi

            I

            2.5 ± 0.14

            efgh

            dcuR

            Z2361ΔdcuR

            +++

            0.198 ± 0.01

            k

            I

            2.2 ± 0.17

            ijklm

            dcuS

            Z2362ΔdcuS

            +++

            0.171 ± 0.02

            kl

            I

            2.07 ± 0.32

            ijklm

            dcuSR

            Z2361–62ΔdcuSR

            +++

            0.256 ± 0.02

            j

            I

            2.23 ± 0.38

            mnop

            eamI

            Z2195ΔeamI

            +++

            0.256 ± 0.046

            j

            I

            2.23 ± 0.12

            ijklm

            eamR

            Z2194ΔeamR

            +++

            0.544 ± 0.1

            de

            I

            2.03 ± 0.31

            pqr

            envZ

            Z0270ΔenvZ

            +++

            1.55 ± 0.06

            c

            C

            1.03 ± 0.11

            v

            envZ/ompR

            Z0270–71 ΔenvZ/ompR

            +++

            1.7 ± 0.05

            a

            C

            1.03 ± 0.11

            v

            ompR

            Z0271ΔompR

            +++

            1.63 ± 0.06

            b

            C

            1.03 ± 0.11

            v

            grrA

            Z2198ΔgrrA

            +++

            1.58 ± 0.15

            bc

            C

            3.15 ± 0.07

            a

            grrS

            Z3742ΔgrrS

            +++

            1.56 ± 0.09

            c

            C

            3.2 ± 0.05

            a

            hrpX

            Z3962ΔhrpX

            +++

            0.69 ± 0.11

            d

            I

            1.73 ± 0.15

            rs

            hrpY

            Z3963ΔhrpY

            +++

            0.471 ± 0.05

            e

            I

            2.0 ± 0.0

            opqr

            hrpXY

            Z3962–63 ΔhrpXY

            +++

            0.45 ± 0.09

            ef

            I

            2.33 ± 0.15

            hijkl

            hrpXYS

            Z3962–64 ΔhrpXYS

            ---

            0.4 ± 0.09

            gh

            I

            2.27 ± 0.15

            hijkl

            kdpD

            Z0414ΔkdpD

            +++

            0.372 ± 0.02

            hgi

            I

            2.13 ± 0.12

            nopq

            kdpE

            Z1912ΔkdpE

            +++

            0.09 ± 0.035

            no

            I

            1.43 ± 0.19

            tu

            luxP

            Z3324ΔluxP

            +++

            0.419 ± 0.12

            g

            I

            1.87 ± 0.23

            rs

            luxQ

            Z3327ΔluxQ

            +++

            0.625 ± 0.11

            de

            I

            1.9 ± 0.17

            qr

            narP

            Z3390ΔnarP

            +++

            0.1 ± 0.001

            lmn

            I

            2.2 ± 0.0

            klmno

            narQ

            Z3389ΔnarQ

            +++

            0.08 ± 0.003

            lmn

            I

            2.4 ± 0.1

            fghij

            narQP

            Z3389–90 ΔnarQP

            +++

            0.08 ± 0.002

            nop

            I

            2.27 ± 0.25

            ghijkl

            phoB

            Z1675ΔphoB

            +++

            0.145 ± 0.003

            klmn

            I

            2.6 ± 0.2

            defg

            phoR

            Z1676ΔphoR

            +++

            0.366 ± 0.01

            hgi

            I

            2.75 ± 0.15

            cd

            phoRB

            Z1675–76 ΔphoRB

            +++

            0.312 ± 0.01

            hi

            I

            2.7 ± 0.1

            cde

            phoP

            Z2324ΔphoP

            +++

            0.15 ± 0.02

            klm

            I

            2.45 ± 0.17

            cdef

            phoQ

            Z2323ΔphoQ

            +++

            0.106 ± 0.01

            lmno

            I

            2.25 ± 0.23

            efghi

            phoQP

            Z2323–24 ΔphoQP

            +++

            0.18 ± 0.002

            jk

            I

            2.25 ± 0.15

            ijklm

            pmrA

            Z0089ΔpmrA

            +++

            0.05 ± 0.02

            pqr

            I

            2.3 ± 0.19

            hijkl

            pmrB

            Z0090ΔpmrB

            +++

            0.10 ± 0.01

            lmno

            I

            2.40 ± 0.17

            defgh

            pmrBA

            Z0089–90 ΔpmrBA

            +++

            0.07 ± 0.01

            opq

            I

            2.3 ± 0.20

            fghij

            rcsB

            Z3206ΔrcsB

            ---

            0.002 ± 0.0005

            r

            I

            2.51 ± 0.17

            efghi

            rcsC

            Z3207ΔrcsC

            ---

            0.4 ± 0.04

            fg

            I

            2.45 ± 0.2

            fghij

            rcsD (yojN)

            Z3205ΔrcsD

            ---

            0.002 ± 0.0002

            r

            I

            2.48 ± 0.28

            efghi

            rcsDB

            Z3205–06 ΔrcsBD

            ---

            0.003 ± 0.0001

            r

            I

            2.43 ± 0.11

            fghijk

            rstA

            Z2651ΔrstA

            +++

            0.359 ± 0.045

            hgi

            I

            1.6 ± 0

            st

            rstB

            Z2662ΔrstB

            +++

            0.1 ± 0.032

            lmn

            I

            2.03 ± 0.06

            opqr

            spk1

            Z1250Δspk1

            +++

            0.17 ± 0.04

            kl

            I

            2.58 ± 0.14

            cdef

            ybjN

            Z2074ΔybjN

            +++

            1.65 ± 0.12

            a

            I

            2.05 ± 0.26

            lmno

            yciR

            Z0955ΔyciR

            +++

            0.111 ± 0.0012

            lmno

            I

            1.9 ± 0.08

            pqr

            yddV

            Z1215ΔyddV

            +++

            0.043 ± 0.001

            qr

            C

            2.9 ± 0.14

            bc

            yegE

            Z3113ΔyegE

            +++

            0.22 ± 0.007

            j

            I

            2.9 ± 0.14

            bc

            yehT

            Z3996ΔyehU

            +++

            0.186 ± 0.026

            k

            I

            2.5 ± 0.2

            efghi

            yehU

            Z3995ΔyehT

            +++

            0.15 ± 0.04

            klm

            I

            2.33 ± 0.21

            fghij

            yehUT

            Z3995–96 ΔyehUT

            +++

            0.15 ± 0.004

            klm

            I

            2.4 ± 0.17

            fghij

            yfhA

            Z3498ΔyfhA

            +++

            0.07 ± 0.012

            opq

            I

            2.3 ± 0.1

            hijkl

            yfhK

            Z3495ΔyfhK

            +++

            0.06 ± 0.002

            pq

            I

            2.33 ± 0.25

            fghij

            yfhKA

            Z3495–98 ΔyfhKA

            +++

            0.062 ± 0.0005

            pq

            I

            2.17 ± 0.12

            lmno

            yoaD

            Z0355ΔyoaD

            +++

            0.22 ± 0.008

            j

            I

            2.4 ± 0.37

            ghijkl

            a ams operon encodes amylovoran biosynthesis. The entire operon (amsA-L; 15.8 kb) was deleted in the ams mutant [47]. flhD and fliA encode two regulators of flagella biosynthesis. These three mutants were used as negative controls for amylovoran biosynthesis and swarming motility, respectively.

            b Virulence assay was performed using immature pear fruits as described. +++: pathogenic as WT; ---: non-pathogenic.

            c Amylovoran production was measured 48 hr post inoculation. ND: not determined.

            d Distance of swarming motility was measured 48 hr post inoculation. Original circle was about 0.6 cm in diameter. I: irregular movement; C: circular. ND: not determined.

            e Values within the amylovoran production and distance of swarming motility columns followed by the same letter are not significantly different according to one way ANOVA and student-Newmans-Kleus test (P = 0.05).

            For prokaryotes, genes for cognate pairs of HKs and RRs are typically encoded together in a single operon [1]. Because of this unique feature for HKs and RRs, and the potential cross-talk between different HKs and RRs, these single deletion mutants generated in this study have been designed to have little or no effect on expression/function of the cognate gene in the operon. Furthermore, operon deletion mutants have been generated so that phenotypes of single gene deletion mutants could be compared with those of operon deletion mutants to determine whether or not there is cross-talk between HKs and RRs. All mutants and their designations are listed in Table 3.

            To determine the role of deletion mutants in E. amylovora virulence, an assay, routinely used to evaluate virulence of E. amylovora WT strains, of all 59 TCST mutants on immature pear fruits was conducted [19, 46, 48]. Mutants and WT strains were inoculated on immature pears as described, and disease development was assessed for up to eight days. Results showed that five mutants, including three single (rcsB, rcsC, rcsD mutant), one double (rcsBD), and one triple mutant (hrpXYS), were non-pathogenic on immature pears (Table 3). However, hrpX, hrpY, and hrpXY mutants remained to exhibit full-virulence, as that of WT, on immature pear fruit (Table 3). Previous genetic studies have shown that a mutation in rcsB renders E. amylovora nonpathogenic and abolishes amylovoran production [49]. Our results further demonstrated that the RcsCDB phosphorelay system was essential for virulence in E. amylovora [17]. In addition, all other mutants induced similar symptoms to the WT strain on immature pear fruit (Table 3).

            Phenotypic analysis of deletion mutants

            All generated mutants did not exhibit any distinctive phenotypes on rich LB medium, except for a single mutant which was mucoid (ybjN mutant). Therefore, phenotypes of the 59 deletion mutant strains were evaluated using a swarming plate assay [50, 51]. WT cells can swarm via the combined effects of flagellar motility, chemotaxis, and growth, thus creating a circular colony. Defects in cell motility, chemotaxis, or growth can produce alterations in swarming size or density [51]. For E. amylovora WT and mutant strains, bacterial suspensions were plated onto swarming agar plates containing 0.3% agar, as previously described [50, 51]. Swarming diameter and density were determined following incubation at 28°C for up to 72 hr, and three different swarming phenotypes were identified. Two mutants, grrA and grrS, exhibited substantially larger and lower density swarms than the WT (Figure 1A; Table 3, significance level a). In addition, three mutants, envZ, ompR, and envZ/ompR double-mutant, showed dramatic reduction in swarm size, but with increased density (Figure 1B; Table 3, significance level v). All these five mutants exhibited circular swarming (Figure 1A and 1B). A third group of mutants showed smaller swarms, colonies exhibited irregular circular patterns, and these included arcAB, baeRS, cpxA1R1, kdpE, luxPQ, rstAB, and yciR mutants (Figure 1C; Table 3, significance levels o to u). Although the distance of the cpxA1R1 mutants had similar significant level as that of envZ/ompR mutants at 48 hr, these mutants were included in the third group due to its irregular movements and with increased distance at 72 hr in contrast to envZ/ompR mutants which remained the same at 72 hr (Figure 2B). The rest of mutants showed even less change in swarm distance. As negative controls, all flagella-deficient mutants, flhDC and fliA, were non-motile on swarming plates (Figure 1D).
            http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-10-245/MediaObjects/12864_2008_Article_2129_Fig1_HTML.jpg
            Figure 1

            Comparison of motility on swarming plates for WT and TCST mutants. Bacterial strains were spotted on the swarming plate (0.3% agar) and incubated at 28°C. Photos were taken at one or two days post inoculation. DPI: days post inoculation. A1 to D1: WT strain; A2: grrS mutant; A3: grrA mutant; B2: envZ mutant; B3: ompR mutant; B4: envZ/ompR double mutant; C2 to C4: representative mutants with irregular movements at two days post inoculation; D2 to D4: flhD, flhC and fliA mutant, respectively. Flagella mutants were used as negative controls.

            http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-10-245/MediaObjects/12864_2008_Article_2129_Fig2_HTML.jpg
            Figure 2

            Comparison of the swarming distance of WT and TCST mutants. The diameters of the swarming circle were measured 24, 48 and 72 hrs after incubation. The experiments were repeated at least three times.

            As shown in Figure 1A, grrS/grrA mutant strains showed precocious swarming and moved faster than the WT strain, with a diameter of 1.7 cm within 24 hr compared to 1 cm for the WT strain, from an original spot of 0.6 cm in diameter (Figure 2A). The swarming distance was 3.2 and 4.8 cm for grrS/grrA mutants and 2.7 and 3.8 cm for WT strain, respectively, at 48 and 72 hr (Figure 2A). When grrS and grrA mutants were complemented with their own gene cloned from the WT strain, the hyper-motile phenotype was restored in these mutants with a swarming distance, ~2.8 cm after 48 hr of incubation, similar to that of the WT strain (data not shown). In contrast, envZ/ompR mutants showed reduced swarming even after 72 hr (Figure 2B), and colonies exhibited more dense and fuzzy appearance. The diameter for envZ/ompR mutants was about 0.8, 1.0, and 1.0 cm for 24 to72 hr, from an original spot of 0.6 cm in diameter, respectively (Figure 2B). When envZ/ompR mutants were complemented by their own gene cloned from the WT strain, the 'low-motile' phenotype was partially restored. The diameter (about 1.5 cm after 48 hr of incubation) for the complemented mutant strains was comparable to that of the WT strain (Data not shown). However, the swarming phenotype remained different between envZ/ompR mutants and those of flhDC and fliA mutants, which exhibited no movement and maintained the original size of 0.6 to 0.7 cm in diameter (Figure 1BD). In summary, these results demonstrated that EnvZ/OmpR acts as a positive regulator of swarming motility, while GrrS/GrrA is a negative regulator of swarming motility in E. amylovora. These findings demonstrated that additional TCST genes might also control swarming motility in E. amylovora.

            Regulation of amylovoran biosynthesis

            Based on the virulence assay, the RcsCDB system that regulates amylovoran biosynthesis is essential for virulence in E. amylovora. Phenotypes of the 59 deletion mutant strains were further evaluated by measuring amylovoran biosynthesis in vitro. The amylovoran concentration in the supernatant of bacterial cultures was quantitatively determined using a turbidity assay with cetylpyrimidinium chloride (CPC), as previously described [52]. In this screening, four groups of mutants exhibiting varying levels of amylovoran production in vitro were identified (Table 3). For Group I, three mutants (rcsB, rcsD, and rcsBD) exhibited lack of amylovoran production (Table 3, significance level r); For Group II, six mutants (grrA, grrS, ybjN, envZ, ompR, and envZ/ompR double-mutant) showed a dramatic increase in amylovoran production, about 20-fold, compared to that of the WT (Figure 3A; Table 3, significance levels a to c); Group III mutants demonstrated increased production of amylovoran, from 2- to 8-fold, compared to that of the WT, and these included hrpX, hrpY, hrpXY, hrpXYS, rcsC, luxPQ, and eamIR (Table 3, significance levels d to j). Group IV mutants did not show much difference in amylovoran production when compared to that of WT, and these included pmrAB, narPQ, and yfhAK (Table 3, significance levels k to p). As a negative control, an amylovoran operon deletion mutant, ams, lacked amylovoran production (Figure 3A; Table 3, significance level r) [17, 47].
            http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-10-245/MediaObjects/12864_2008_Article_2129_Fig3_HTML.jpg
            Figure 3

            TCSTs regulate amylovoran biosynthesis and gene expression. (A) Amylovoran production of E. amylovora WT and TCST mutants in vitro. Bacterial strains were grown in MBMA media with 1% sorbitol for 48 hrs at 28°C with shaking. The amount of amylovoran was measured with the CPC assay and normalized to a cell density of 1. Amylovoran operon (amsA-L) deletion mutant (Δams) was used as a negative control [47]. (B) Gene expression of the amsG gene in WT and TCST mutants in vitro. GFP intensity in WT and TCST mutants containing amsG promoter-gfp fusion plasmid was measured by flow cytometry. GFP-A: Green fluorescence protein absorbance; Count: Number of cells.

            For group I mutants, we have recently described how the RcsCDB phosphorelay system regulates amylovoran biosynthesis [17]. This study focused on group II mutants, except for the ybjN mutant, which will be reported in the future. Amylovoran production of five mutants, including grrA, grrS, envZ, ompR, and the envZ/ompR double-mutant, was partially complemented by their own gene/operon cloned from the WT strain. As shown in Figure 3A, complemented strains produced slightly less amylovoran than that of the mutant strains. To correlate amylovoran production with amylovoran biosynthesis gene expression [17], the level of expression of the amsG gene, the first gene in the ams (amylovoran) operon, was measured using GFP as a reporter in two RR mutants (grrA and ompR). GFP intensity was measured in the WT and in mutants containing the amsG promoter-GFP fusion using flow cytometry [17]. The amsG gene was expressed at a basal level in the WT strain, with a GFP intensity value of 1.7 (geometric mean), compared to a geometric mean value of 1.5 for the control vector. The geometric mean value of the GFP intensity of the amsG promoter was 26.2 and 60.5 for ompR and grrA mutants, respectively (Figure 3B). These results indicated that amylovoran production was negatively regulated by both EnvZ/OmpR and GrrA/S systems, and a regulatory network might be regulating amylovoran production in E. amylovora.

            Discussion

            In this study, 46 TCSTs have been identified in E. amylovora, a pathogen of rosaceous plants which mainly resides in the plant xylem, but can also grow epiphytically on stigmas of flowers. Compared to other plant pathogenic enterobacteria such as D. dadantii and P. carotovora subsp. atroseptica, both of which can survive not only in plants, but also in soil, or the animal counterpart Es. coli, the number of TCSTs present in E. amylovora is relatively small, thus may reflecting the particular host niche for this pathogen. Furthermore, the genome size (3.9 Mbp) of E. amylovora is also the smallest among sequenced enterobacterial plant and mammalian pathogens [53]. However, E. amylovora still maintains both enterobacterial-specific TCSTs (such as the Rcs system) [54] and plant enterobacterial-specific ones (such as HrpXY) [39].

            On the other hand, E. amylovora has almost the same number of TCSTs as that found in E. tasmaniensis, a saprophytic bacterium isolated from apple flowers in Tasmania, Australia [38]. Since both bacteria occupy the same ecological niche during colonization and can grow epiphytically on flowers, these bacteria might have developed and maintained similar TCSTs. This suggests that TCSTs may be evolutionarily maintained to cope with similar environmental and plant host signals in closely related bacteria. It is interesting to note that fire blight is endemic to North America, and it has subsequently spread to Europe and New Zealand in the 1950s and in 1917, respectively, but it has not yet been reported in Australia [14]. It is possible that these two bacteria may have not yet encountered each other, as the distribution of E. tasmaniensis outside of Australia remains unknown. Surprisingly, we found that the majority of TCSTs in these two bacteria share a high level of aa similarity and identity. These results suggest that similar TCSTs might play important roles for the survival and proliferation among closely-related bacteria in similar plant niches.

            Currently, several thousands of TCSTs have been identified in sequenced bacterial genomes [5557]. Although the basic biochemistry of TCSTs is well understood, some structural insights in phosphorylation-dependent changes of TCST domains are variable. Domain architecture has proven particularly informative for analyzing multidomain proteins involved in signal transduction and in predicting the functions of these signal transduction proteins [43, 58]. In RRs, structural characterization of DNA-binding domains has revealed several variations on the common helix-turn-helix (HTH) theme, such as NarL- and OmpR-types. Recently, many novel conserved domains have been described such as PAS, GAF, GGDEF, EAL, and HD-GYP, thereby affirming the complexity of bacterial signaling systems [43, 58]. In E. amylovora, most HKs belong to four common HK groups, and are periplasmic sensors. Moreover, most RRs in this bacterium are OmpR-, NarL-, NtrC-, and LytR-like proteins, the four most common families of DNA-binding RRs found in prokaryotes [43]. These results indicate that E. amylovora has maintained some basic signal transduction systems for the bacterium to survive. Phylogenetic and genomic analyses have revealed co-evolutionary relationships between cognate HKs and RRs. This seems obvious in E. amylovora as orphan HKs and RRs are rare.

            It is well understood that TCSTs are involved in regulating virulence gene expression in plant bacterial pathogens [10]. Previous reports on TCSTs have demonstrated the importance of TCSTs in the virulence of bacterial plant pathogens. However, data on complete or global virulence regulation networks are lacking. A recent study by Qian et al. [13] provides a useful beginning towards a better understanding of the regulatory networks involved. A genome-wide mutagenesis of all 54 RRs in X. campestris pv. campestris has revealed that two novel RRs are involved in virulence, thus facilitating future studies on signaling networks in this bacterium [13]. In this study, we have utilized a reverse genetic approach and constructed 59 HK and RR mutants in E. amylovora, which will also provide valuable tools for future global gene expression assays using microarrays to deduce signaling networks in this bacterium.

            Early studies have revealed that in E. amylovora, the Hrp T3SS, which delivers effector proteins into host plants, and the EPS amylovoran are two major virulence factors [52, 59, 60]. Previous reports have also indicated that the RcsCDB phosphorelay system regulates amylovoran biosynthesis, while the two-component system HrpXY regulates hrp-T3SS gene expression. Recently, we have further demonstrated that the Rcs system is essential for virulence in E. amylovora and may play a role in the survival of the pathogen [17]. Mutations in the Rcs system have rendered the organism non-pathogenic [17]. In this study, we have found that hrpX, hrpY and hrpXY mutants remain virulent, and could induce a spotty weak hypersensitive response (HR) on tobacco (Figure S2A, [see Additional file 2]); while, a hrpXYS triple mutant has a normal Hrp- phenotype. It is interesting to note that, in a previous report, Tn5-insertional mutants of hrpY have been reported to be non-pathogenic, and could not induce an HR on tobacco [18]. Two classes of hrpX insertional mutants have been identified, one similar to the hrpY mutant and the other that continue to cause disease and induce a spotty HR on tobacco [18]. Similar observations have been reported in Pantoea stewartii subsp. stewartii, causal agent of Stewart's wilt of corn [61]; wherein, Tn5-insertional hrpX and hrpY mutants exhibit a Hrp- phenotype and in-frame deletion hrpX mutants show reduced virulence. Since the hrpXY is transcribed as an operon, it is possible that Tn5 insertion could cause polar effects on the downstream genes such as hrpS (Figure S2B, [see Additional file 2]). Indeed, a recent study indicates that the hrp regulatory genes in P. stewartii subsp. stewartii participate in a novel regulatory loop that upregulates itself by readthrough transcription of hrpL into hrpXYS [62].

            It has been proposed that in E. amylovora, both HrpY and HrpS regulate hrpL, encoding the master regulator of T3SS, and that the effects of HrpY and HrpS are additive [18]. Subsequent studies of P. stewartii subsp. stewartii, D. dadantii, and Pantoea herbicola pv. gyposophilae have demonstrated that HrpY initially activates hrpS by binding to its promoter, and then the HrpL is activated by HrpS as well as by other regulatory genes [39, 63]. Furthermore, it has been reported that HrpY in D. dadantii acts as both a positive and negative regulator [39]. In P. stewartii subsp. stewartii, besides the HrpY binding site, additional sequences further upstream of the hrpS promoter are also required for hrpS expression, suggesting that unknown regulatory proteins may act cooperatively with HrpY [63]. Microarray studies suggest that hrpL represents only one branch of the regulatory pathways downstream of hrpRS, and a large number of genes regulated by HrpRS are hrpL-independent in P. syringae [64]. In our study, HrpX, HrpY, and HrpS also act as negative regulators, as hrpX, hrpY, hrpXY, and hrpXYS mutants produce more amylovoran. Our virulence tests suggest that the T3SS remains functional in hrpX, hrpY, and hrpXY deletion mutants, but not in the hrpXYS mutant, whereby HrpS and HrpL may be activated by other unknown regulators except that production and/or translocation of HrpN (Harpin) in tobacco is severely attenuated as showed in an HR assay (Figure S2A, [see Additional file 2]). It is possible that host signals affecting gene expression may also be different in tobacco.

            Various models have proposed that HrpX senses environmental signals in the plant apoplast or the Hrp-inducing medium to phosphorylate HrpY [18, 34, 35]. However, domain structure analysis has indicated that HrpX is a soluble cytoplasmic protein, and may sense intracellular signals. This suggests that other signaling pathways may also be involved in activating hrpXY, hrpS, or hrpL by sensing outside signals to regulate T3SS. Since the expression of HrpX, HrpS, and HrpL is regulated by low pH, also corresponding to conditions under which OmpR-EnvZ and GrrS-GrrA are activated, our results further suggest that both OmpR-EnvZ and GrrS-GrrA can regulate the hrpXY operon or hrpS either directly or indirectly as reported in other plant pathogenic bacteria [10, 28, 33, 34]. Further studies are needed to dissect the roles of OmpR-EnvZ and GrrS-GrrA in regulating T3SS.

            Several regulatory genes have been previously reported to control amylovoran biosynthesis in E. amylovora, including the Rcs system, RcsA, Lon protease, and H-NS protein [17, 50, 65]. Here, we have further identified several groups of regulators, including both negative and positive regulators. These regulators may form a network that governs the production of amylovoran under different conditions to benefit pathogen survival or pathogenesis. Regulatory cascades are also likely to occur as global regulators such as OmpR-EnvZ and GrrS-GrrA may control expression of other regulatory genes or proteins, such as hrpXY and quorum sensing systems, as reflected in the amount of amylovoran produced in these mutants. However, we cannot rule out that cross-talk between different TCSTs may further complicate this scenario. Our study indicates that regulation of amylovoran biosynthesis is highly complex and further suggests that the pathogen has developed a system to control this major virulence factor.

            Swarming is a flagella-driven form of motility for movement across solid surfaces as a group [6670]. Swarmer cells are normally hyperflagellated and require extracellular components such as EPS and surfactants that enable mass migration [71, 72]. Previous studies have identified several global regulators in Es. coli and in other bacteria, including the Rcs system, OmpR and GrrSA, known to influence flagella biosynthesis, especially the master regulator flhDC [71, 73]. In this study, we have identified both negative (GrrSA) and positive (EnvZ/OmpR) regulators of swarming motility. It is easy to accept that the GacSA system may negatively regulate flagella biosynthesis, thus rendering grrSA mutants hypermotile. It is also obvious that flagella biosynthesis is not impaired in the envZ-ompR mutants as the swarming phenotype is different between the envZ-ompR mutants and the flhDC-fliA mutants. The obvious question that arises as to why envZ-ompR mutants are non-motile although they produce prolonged flagella. A recent study in Salmonella typhimurium has reported that mutations in chemotaxis pathways are impaired for swarming motility and further revealed a role of flagellum in sensing external wetness [72]. It has been proposed that swarming requires a fluid environment generated as bacteria extract water from the underlying agar gel [74]. Flagella are designed to work in this aqueous environment; that is swarming cells move in a thin layer of fluid over the surface of the agar. Further studies have revealed that the wetting agent that draws water out of the underlying agar is an osmotic agent [74]. It is likely that in our envZ-ompR mutants, this unknown osmotic agent provides a signal for EnvZ-OmpR system that may regulate chemotaxis response or creates high osmolarity so that water can be removed from the agar, thus affecting swarming motility. Further studies are needed to clarify this hypothesis.

            Conclusion

            In summary, we have identified and classified TCSTs, and have systematically generated TCST deletion mutants in E. amylovora. The mutants generated will serve as resources for future exploration of TCSTs in E. amylovora. In-depth characterization of deletion mutants and global gene expression will be our next goal. Our current data provide experimental evidence that TCSTs, especially those global regulators such as the Rcs system, EnvZ-OmpR and GrrA-GrrS, play important roles in virulence and in regulating virulence factors such as amylovoran production and swarming motility. These findings also suggest that multiple TCSTs may form complex and highly connected circuits and signaling networks in this important pathogen.

            Methods

            Bacterial stains and culture media

            Bacterial strains and plasmids used in this study are listed in Table 4. The LB medium was used routinely for culturing E. amylovora. When necessary, the following antibiotics were added to the medium: 20 μg ml-1kanamycin and chloramphenical, and 100 μg ml-1 ampicillin. Amylovoran production was determined by growing bacteria in MBMA medium (3 g KH2PO4, 7 g K2HPO4, 1 g [NH4]2SO4, 2 ml glycerol, 0.5 g citric acid, 0.03 g MgSO4) amended with 1% sorbitol [52, 75].
            Table 4

            Bacterial strains and plasmids used in this study

            Strains or plasmids

            Relevant charactersa

            Reference or source

            E. amylovora strains

              

            Ea1189

            Wild type, isolated from apple

            81

            Z2946ΔflhD

            flhD::Km; Kmr-insertional mutant of flhD of Ea1189, Kmr

            This study

            Z2945ΔflhC

            flhC::Km; Kmr-insertional mutant of flhC of Ea1189, Kmr

            This study

            Z2986ΔfliA

            fliA::Km; Kmr-insertional mutant of fliA of Ea1189, Kmr

            This study

            Z0118Δams

            KmR-deletion mutant of ams operon (15.8 kb) of Ea1189, KmR

            47

            E. coli

              

            DH10B

            F- mcrA Δ(mrr-hsdRMS-mcrBC) Φ80lacZΔM15 ΔlacX74 recA1 endA1 araΔ139 Δ(ara, leu)7697 galU galK λ - rpsL (StrR) nupG

            Invitrogen, Carlsbad, CA, USA

            Plasmids

              

            pKD46

            Apr, PBAD gam bet exo pSC101 oriTS

            45

            pKD13

            Kmr, FRT cat FRT PS1 PS2 oriR6K rgbN

            45

            pKD3

            CmR, FRT cat FRT PS1 PS2 oriR6K rgbN

            45

            pGEM® T-easy

            Apr, PCR cloning vector

            Promega, Madison, WI, USA

            pWSK29

            Apr, cloning vector, low copy number

            82

            pFPV25

            Apr, GFP based promoter trap vector with a promoterless gfpmut3a gene

            83

            pWDP4

            721 bp KpnI-XbaI DNA fragment containing promoter sequence of amsG gene in pFPV25

            17

            pWDP5

            A 3.057 kb PCR fragment containing grrS gene in pGEM T-easy vector

            This study

            pSN4

            A 2.727 kb PCR fragment containing envZ/ompR operon in pGEM T-easy vector

            This study

            pSN5

            A 2.727 kb PCR fragment containing envZ/ompR operon in pWSK29 vector

            This study

            pSN6

            A 1.6 kb PCR fragment containing grrA gene in pGEM T-easy vector

            This study

            pSN7

            A 1.463 kb PCR fragment containing flhDC operon in pWSK29 vector

            This study

            pSN8

            A 940 bp PCR fragment containing fliA gene in pWSK29 vector

            This study

            a KmR, ApR, CmR = kanamycin, ampicillin and chloramphenical resistance, respectively.

            DNA manipulation and bioinformatic analysis

            Plasmid DNA purification, PCR amplification of genes, isolation of fragments from agarose gels, cloning, restriction enzyme digestion and T4 DNA ligation were performed using standard molecular procedures [76]. DNA sequencing was performed at the Keck Center for Functional and Comparative Genomics at the University of Illinois at Urbana-Champaign. Sequence management and contig assembly were conducted using Sequencher 4.9 software. Similarity searches of nucleic acid sequences and proteins were conducted using the BLASTN and BLASTP programs at NCBI http://​www.​ncbi.​nlm.​nih.​gov/​BLAST[77]. Protein domain organizations of the HKs and RRs were identified by searching the CDD with Reverse Specific Position BLAST [41] or the SMART program http://​smart.​embl-heidelberg.​de/​[42]. Domain limits for proteins were also derived from the graphical output of the SMART web interface. Transmembrane segments of HKs were predicted by the TMHMM2 program http://​www.​cbs.​dtu.​dk/​services/​TMHMM/​[78]. Multiple sequence alignments of proteins were carried out using Clustal ×1.80 [79]. Sequence logos were derived from alignment of proteins belonging to corresponding HK group and created by Weblogo http://​weblogo.​berkeley.​edu. Proteins without known domain were characterized by using PSI-BLAST of NCBI and methods described previously [80].

            Construction of TCST mutants in E. amylovora by Lambda-Red recombinase

            E. amylovora stable mutants were generated by using the λ phage recombinases as previously described [4547]. Briefly, E. amylovora Ea1189 was transformed with plasmid pKD46 expressing recombinases redα, β, and γ. The transformant Ea1189 (pKD46) was grown overnight at 28°C, reinoculated in LB broth containing 0.1% arabinose, and grown to exponential phase OD600= 0.8. Cells were collected, made electrocomponent, and stored at -80°C. Recombination fragments consisting of a kanamycin (kan) or chloramphenical (cat) gene with its own promoter, flanked by a 50-nucleotide (nt) homology arm, were generated by PCR using plasmid pKD13 or pKD3 as a template. The primers used for generating mutants are listed in Table S2 [see Additional file 1] (primers start with B). To confirm mutants by PCR, primers (Table S2, [see Additional file 1], primers start with Z), and internal primer pairs km1 and km2 of the kan gene, cm1 and cm2 of the cat gene were used. For the resulting mutants, the majority of the coding region of each gene was replaced by the kan or cat gene, except for the first and last 50 nt. The resulting mutants were designated and listed in Table 1 and Table 3.

            Cloning genes/operon for complementation of TCST mutants

            For complementation of selected mutants, flanking sequences of genes or operon were used to design primers to amplify fragments of gene/or operon and their promoter sequences (Table S2, [see Additional file 1], primers with gene name). Primer pairs, with or without restriction sites were used to amplify corresponding DNA fragments from E. amylovora WT strain which contains upstream and downstream sequences of the gene or operon, and PCR fragments were cloned into either pGEM T-easy vector through A-T ligation or into a low-copy number vector pWSK29. The final plasmids were listed in Table 4. For primer pair with restriction site, DNA fragments and the vector were both digested with corresponding enzymes following amplification, and ligated together. All plasmids were introduced into E. amylovora strain by electroporation. Transformants were selected on LB plates supplemented with Ap and Km or Cm. Their genotypes were confirmed by both enzymatic digestion and sequencing.

            Construction of promoter-GFP fusions for gene expression assays

            For gene expression assays, flanking sequences of the amsG ORF were used to design primers to amplify DNA fragments. Primer pairs amsG1-amsG2 with restriction sites were used to amplify 721 bp DNA fragments from E. amylovora WT strain, containing promoter sequences of amsG gene. PCR products and the promoter trapping vector pFPV25 were both digested with BamH1 and EcoRI for amsG gene. The resulting fragments were gel-purified, ligated together, and cloned to the upstream of promoter-less gfp gene. The final plasmids were designated as pWDP4 for amsG, and were confirmed by restriction enzyme digestion and sequencing.

            Immature pear virulence assays

            For E. amylovora WT and mutant strains, bacterial suspensions were grown overnight in LB broth, harvested by centrifugation, and resuspended in 0.5 × sterile phosphate buffered-saline (PBS) with bacterial cells adjusted to concentrations of ~1 × 103 to 1 × 104 colony-forming units (CFU/μl) (OD600 = 0.1 and then diluted 100 times) in PBS. Immature fruits of pear (Pyrus communis L. cv. 'Bartlett') were surface-sterilized, and pricked with a sterile needle as described previously [19, 46]. Wounded fruits were inoculated with 2 μl of cell suspensions, and incubated in a humidified chamber at 26°C. Symptoms were recorded at 2, 4, 6, and 8 days post-inoculation.

            CPC assay for determining amylovoran concentration

            The amylovoran concentration in supernatants of bacterial cultures was quantitatively determined by a turbidity assay with cetylpyrimidinium chloride (CPC), as previously described [50, 52]. Briefly, for E. amylovora WT, mutants, and complementation strains, bacterial suspensions were grown overnight in LB broth w/o appropriate antibiotics, harvested by centrifugation, and washed with PBS three times. After the final wash, the bacterial pellet was resuspended in 200 μl PBS. A total of 100 μl of bacterial suspension was inoculated into 10 ml MBMA medium with 1% sorbitol. One ml of bacterial cells was pelleted two to three days after inoculation at 28°C with shaking. Following centrifugation, 50 μl CPC at 50 mg ml-1 was added to one ml supernatant. After 10 min of incubation at room temperature, the amylovoran concentration was determined by measuring OD600 turbidity. The final concentration of amylovoran production was normalized for a cell density of 1.0. For each strain tested, the experiment was repeated at least three times.

            Swarming motility assay

            For E. amylovora WT and mutant strains, bacterial suspensions were grown overnight in LB broth w/o appropriate antibiotics, harvested by centrifugation, washed with PBS once, and resuspended in 200 μl PBS. Then, bacterial suspensions were diluted 10 × in water, and 5 μl of the diluted bacterial suspension was plated onto the center of swarming agar plates (10 g tryptone, 5 g NaCl, 3 g agar per l Liter) as previously described [50, 51]. Swarming diameters were determined following incubation at 28°C for up to three days. The experiments were repeated at least three times.

            Statistical analysis

            Amylovoran production and swarming motility data were undertaken a one-way ANOVA and Student-Newman-Keuls test to determine differences in means (P = 0.05) using SAS 9.1 program.

            GFP reporter gene assay by flow cytometry

            The BD FACSCanto flow cytometer was used to monitor the GFP intensity of WT and mutant strains containing the corresponding promoter-gfp construct [17]. For in vitro amsG gene expression, WT and mutants containing the amsG promoter-GFP fusion plasmid were grown in LB overnight, harvested, and resuspended in water. Bacterial suspensions were re-inoculated in MBMA broth with 1% sorbitol and grown at 28°C with shaking for up to three days. Bacterial cultures were then harvested by centrifugation, washed once with PBS, and then resuspended in PBS for flow cytometry assay. Flow cytometry was performed on a BD LSRII 10 parameter multilaser analyzers (BD Bioscience, San Jose, CA). For both cases, data were collected for a total of 100,000 events and statistically analyzed by gating using flow cytometry software FCS Express V3 (De Novo Software, LA, CA). A geometric mean was calculated for each sample. Each treatment was performed in triplicate and each experiment was repeated three times.

            Abbreviations

            CDD: 

            Conserved Domain Databases

            CPC: 

            cetylpyrimidinium chloride

            EAL: 

            Type I c-di-GMP phosphodiesterase

            EPS: 

            extracellular polysaccharide

            HisKA: 

            histidine kinase domain

            HK: 

            histidine kinase

            HMM: 

            Hidden Markov Model

            Hpt: 

            histidine phosphotransfer

            HTH: 

            helix-turn-helix

            GGDEF: 

            diguanylate cyclase

            PAS: 

            initially found in PER, ARNT, and SIM proteins

            REC: 

            receiver domain

            RR: 

            response regulator

            TCST: 

            two-component signal transduction system

            T3SS: 

            type III secretion system

            Declarations

            Acknowledgements

            We are grateful to Dr. Barry Wanner at Purdue University, Dr. S. Kushner at University of Georgia, and Dr. R. Valdivia at Duke University for providing the pKD46, pKD13 and pKD3, pWSK29 and pFPV25 plasmids, respectively. We would like to thank Dr. Virginia Stockwell at Oregon State University for providing part of the immature pear fruits used in this experiment. We are also grateful to Dr. Steve Beer at the Cornell University and the Sanger Institute for the genome sequence of E. amylovora strain Ea273. This work was supported by the USDA Cooperative State Research, Education and Extension Service (Hatch project # Illu-802-317) and the Agricultural Experiment Stations of Illinois and Michigan.

            Authors’ Affiliations

            (1)
            Department of Crop Sciences, University of Illinois at Urbana-Champaign
            (2)
            Department of Plant Pathology, Michigan State University
            (3)
            Department of Natural Resources and Environmental Sciences, University of Illinois at Urbana-Champaign

            References

            1. Hoch JA, Silhavy TJ, eds: Two-component signal transduction Washington DC: ASM Press 1995.
            2. Inouye M, Dutta R, eds: Histidine Kinases in Signal Transduction New York: Academic Press 2003.
            3. Alex LA, Simon MI: Protein histidine kinases and signal transduction in prokaryotes and eukaryotes. Trends Genet 1994, 10:133–138.View ArticlePubMed
            4. Stock AM, Robinson VL, Goudreau PN: Two-component signal transduction. Annu Rev Biochem 2000, 69:183–215.View ArticlePubMed
            5. Wolanin P, Thomason PA, Stock JB: Histidine protein kinases: key signal transducers outside the animal kingdom. Genome Biol 2002, 3:3013.View Article
            6. Mizuno T: Compilation of all genes encoding two-component phosphotransfer signal transducers in the genome of Escherichia coli. DNA Res 1997, 4:161–168.View ArticlePubMed
            7. Beier D, Gross R: Regulation of bacterial virulence by two-component systems. Curr Opin Microbiol 2006, 9:1–10.View Article
            8. Lavin JL, Kiil K, Resano O, Ussery DW, Oguiza JA: Comparative genomic analysis of two-component regulatory proteins in Pseudomonas syringae. BMC Genomics 2007, 8:397.View ArticlePubMed
            9. Qian W, Han Z, He C: Two-component signal transduction systems of Xanthomonas spp.: A lesson from genomics. Mol Plant-Microbe Interact 2008, 21:151–161.View ArticlePubMed
            10. Mole BM, Baltrus DA, Dangl JL, Grant SR: Global virulence regulation networks in phytopathogenic bacteria. Trends Microbiol 2007, 15:363–371.View ArticlePubMed
            11. Brencic A, Winans SC: Detection and response to signals involved in host-microbe interactions by plant associated bacteria. Microbiol Mol Biol Rev 2005, 69:155–194.View ArticlePubMed
            12. Heeb S, Haas D: Regulatory roles of the GacS/GacA two-component system in plant-associated and other gram-negative bacteria. Mol Plant-Microbe Interact 2001, 14:1351–1363.View ArticlePubMed
            13. Qian W, Han Z, Tao J, He C: Genome-scale mutagenesis and phenotypic characterization of two-component signal transductions of Xanthomonas campestris pv. campestris ATCC33913. Mol Plant-Microbe Interact 2008, 21:1128–1138.View ArticlePubMed
            14. Vanneste JL, ed: Fire blight, the disease and its causative agent, Erwinia amylovora Wallingford UK: CABI Publ 2000.View Article
            15. He SY, Nomura K, Whittam TS: Type III protein secretion mechanism in mammalian and plant pathogens. Biochem Biophys Acta 2004, 1694:181–206.View ArticlePubMed
            16. Oh CS, Beer SV: Molecular genetics of Erwinia amylovora involved in the development of fire blight. FEMS Microbiol Let 2005, 253:185–192.View Article
            17. Wang DP, Korban SS, Zhao YF: The Rcs phosphorelay system is essential for pathogenicity in Erwinia amylovora. Mol Plant Pathol 2009, 10:277–290.View ArticlePubMed
            18. Wei ZM, Kim JF, Beer SV: Regulation of hrp genes and type III protein secretion in Erwinia amylovora by HrpX/HrpY, a novel two component system, and HrpS. Mol Plant-Microbe Interact 2000, 13:1251–1262.View ArticlePubMed
            19. Zhao YF, Blumer SE, Sundin GW: Identification of Erwinia amylovora genes induced during infection of immature pear tissue. J Bacteriol 2005, 187:8088–8103.View ArticlePubMed
            20. Lapouge K, Schubert M, Allain FHT, Haas D: Gac/Rsm signal transduction pathways of γ-proteobacteria: from RNA recognition to regulation of social behaviour. Mol Microbiol 2008, 67:241–253.View ArticlePubMed
            21. Forst S, Delgado J, Inouye M: Phosphorylation of OmpR by the osmosensor EnvZ modulates expression of the ompF and ompC genes in Escherichia coli. Proc Natl Acad Sci USA 1989, 86:6052–6056.View ArticlePubMed
            22. Park D, Forst S: Co-regulation of motility, exoenzyme and antibiotic production by the EnvZ-OmpR-FlhDC-FliA pathway in Xenorhabus nematophila. Mol Microbiol 2006, 61:1397–1412.View ArticlePubMed
            23. Feng X, Oropeza R, Kenney LJ: Dual regulation by phospho-OmpR of ssrA/B gene expression in Salmonella pathogenicity island 2. Mol Microbiol 2003, 48:1131–1143.View ArticlePubMed
            24. Garmendia J, Beuzon CR, Ruiz-Albert J, Holden DW: The roles of SsrA-SsrB and OmpR-EnvZ in the regulation of genes encoding the Salmonella typhimurium SPI–2 type III secretion system. Microbiology 2003, 149:2385–2396.View ArticlePubMed
            25. Kim CC, Falkow S: Delineation of upstream signaling events in the Salmonella pathogenicity island 2 transcriptional activation pathway. J Bacteriol 2004, 186:4694–4704.View ArticlePubMed
            26. Lee AK, Detweiler CS, Falkow S: OmpR regulates the two-component system SsrA-SsrB in Salmonella pathogenicity island 2. J Bacteriol 2000, 182:771–781.View ArticlePubMed
            27. Brzostek K, Brzostkowska M, Bukowska I, Karwicka E, Raczkowska A: OmpR negatively regulates expression of invasion in Yersinia enterocolitica. Microbiology 2007, 153:2416–2425.View ArticlePubMed
            28. Xiao Y, Lan L, Yin C, Deng X, Baker D, Zhou J, Tang X: Two-component sensor RhpS promoters induction of Pseudomonas syringae type III secretion system by repressing negative regulator RhpR. Mol Plant-Microbe Interact 2007, 20:223–234.View ArticlePubMed
            29. Yang S, Peng Q, Zhang Q, Yi X, Choi CJ, Reedy RM, Charkowski AO, Yang C: Dynamic regulation of GacA in type III secretion, pectinase gene expression, pellicle formation, and pathogenicity of Dickeya dadantii (Erwinia chrysanthemi 3937). Mol Plant-Microbe Interact 2008, 21:133–142.View ArticlePubMed
            30. Lebeau A, Reverchon S, Gaubert S, Kraepiel Y, Simond-Cote E, Nasser W, Gijsegem F: The GacA global regulator is required for the appropriate expression of Erwinia chrysanthemi 3937 pathogenicity genes during plant infection. Environ Microbiol 2008, 10:545–559.View ArticlePubMed
            31. Marutani M, Taguchi F, Ogawa Y, Hossain Md M, Inagaski Y, Toyoda K, Shiraishi T, Ichinose Y: Gac two-component system in Pseudomonas syringae pv. tabaci is required for virulence but not for hypersensitive reaction. Mol Genet Genomics 2008, 279:313–322.View ArticlePubMed
            32. Cui Y, Chatterjee A, Yang H, Chatterjee AK: Regulatory network controlling extracellular proteins in Erwinia carotovora subsp. carotovora: FlhDC, the master regulator of flagellar genes, activates rsmB regulatory RNA production by affecting gacA and hexA (lrhA) expression. J Bacteriol 2008, 190:4610–4623.View ArticlePubMed
            33. Chatterjee A, Cui Y, Yang H, Collmer A, Alfano JR, Chatterjee AK: GacA, the response regulator of a two-component system, acts as a master regulator in Pseudomonas syringae pv. tomato DC3000 by controlling regulatory RNA, transcriptional activators, and alternative sigma factors. Mol Plant-Microbe Interact 2003, 16:1106–1117.View ArticlePubMed
            34. Tang X, Xiao Y, Zhou J: Regulation of the Type III secretion system in phytopathogenic bacteria. Mol Plant-Microbe Interact 2006, 19:1159–1166.View ArticlePubMed
            35. Yang S, Peng Q, Francisco MS, Wang Y, Zeng Q, Yang C: Type III secretion system genes of Dickeya dadantii 3937 are induced by plant phenolic acids. PLoS One 2008, 3:e2973.View ArticlePubMed
            36. Gonzalez N, Heeb S, Valverde C, Reimmann C, Junier T, Haas D: Genome-wide search reveals a novel GacA-regulated small RNA in Pseudomonas species. BMC Genomics 2008, 9:167.View ArticlePubMed
            37. Bocsanczy AM, Perna NT, Biehl B, Glasner JD, Cartinhour SW, Schneider DJ, DeClerck GA, Sebaihia M, Parkhill J, Bentley S, Beer SV: Contributions of the genome sequence of Erwinia amylovora to the fire blight community. Acta Hortic 2008, 793:163–170.
            38. Kube M, Migdoll AM, Muller I, Kuhl H, Beck A, Reinhardt R, Geider K: The genome of Erwinia tasmaniensis strain Et1/99, a non-pathogenic bacterium in the genus Erwinia. Environ Microbiol 2008, 10:2211–2222.View ArticlePubMed
            39. Yap MN, Yang C, Charkowski AO: The response regulator HrpY of Dickeya dadantii 3937 regulates virulence genes not linked to the hrp cluster. Mol Plant-Microbe Interact 2008, 21:304–314.View ArticlePubMed
            40. Kim D, Forst S: Genomic analysis of the histidine kinase family in bacteria and archaea. Microbiology 2001, 147:1197–1212.PubMed
            41. Marchler-Bauer A, Anderson JB, Cherukuri PF, DeWeese-Scott C, Geer LY, Gwadz M, He S, Hurwitz DI, Jackson JD, Ke Z, Lanczycki CJ, Liebert CA, et al.: CDD: A Conserved Domain Database for protein classification. Nucleic Acids Res 2005, 33:D192-D196.View ArticlePubMed
            42. Letunic I, Copley RR, Pils B, Pinkert S, Schultz J, Bork P: SMART 5: Domains in the context of genomes and networks. Nucleic Acids Res 2006, 34:D257-D260.View ArticlePubMed
            43. Galperin MY: Structural classification of bacterial response regulators: Diversity of output domains and domain combinations. J Bacteriol 2006, 188:4169–4182.View ArticlePubMed
            44. Fouhy Y, Lucey JF, Ryan RP, Dow JM: Cell-cell signaling, cyclic di-GMP turnover and regulation of virulence in Xanthomonas campestris. Res Microbiol 2006, 157:899–904.View ArticlePubMed
            45. Datsenko KA, Wanner BL: One-step inactivation of chromosomal genes in Escherichia coli K–12 using PCR products. Proc Natl Acad Sci USA 2000, 97:6640–6645.View ArticlePubMed
            46. Zhao YF, He SY, Sundin GW: The Erwinia amylovora avrRpt2EA gene contributes to virulence on pear and AvrRpt2EA is recognized by Arabidopsis RPS2 when expressed in Pseudomonas syringae. Mol Plant-Microbe Interact 2006, 19:644–654.View ArticlePubMed
            47. Zhao YF, Sundin GW, Wang DP: Construction and analysis of pathogenicity island deletion mutants of Erwinia amylovora. Can J Microbiol 55:457–464.
            48. Bogdanove AJ, Bauer DW, Beer SV: Erwinia amylovora secretes DspE, a pathogenicity factor and functional AvrE homolog, through the Hrp (type III secretion) pathway. J Bacteriol 1998, 180:2244–2247.PubMed
            49. Bereswill S, Geider K: Characterization of the rcsB gene from Erwinia amylovora and its influence on exoploysaccharide synthesis and virulence of the fire blight pathogen. J Bacteriol 1997, 179:1354–61.PubMed
            50. Hildebrand M, Aldridge P, Geider K: Characterization of hns genes from Erwinia amylovora. Mol Genet Genomics 2006, 275:310–319.View ArticlePubMed
            51. Skerker JM, Prasol MS, Perchuk BS, Biondi EG, Laub MT: Two-component signal transduction pathways regulating growth and cell cycle progression in a bacterium: A system-level analysis. PLoS Biol 2005, 3:e334.View ArticlePubMed
            52. Bellemann P, Bereswill S, Berger S, Geider K: Visualization of capsule formation by Erwinia amylovora and assays to determine amylovoran synthesis. Int J Biol Macromol 1994, 16:290–296.View ArticlePubMed
            53. Toth IK, Pritchard L, Birch PRJ: Comparative genomics reveals what makes an enterobacterial plant pathogen. Annu Rev Phytopathol 2005, 44:305–306.View Article
            54. Majdalani N, Gottesman S: The Rcs phosphorelay: a complex signal transduction system. Annu Rev Microbiol 2006, 59:379–405.View Article
            55. D'souza M, Glass EM, Syed M, Zhang Y, Rodriguez A, Maltsev N, Galperin MY: Sentra: a database of signal transduction proteins for comparative genome analysis. Nucleic Acids Res 2007, 35:D271–273.View ArticlePubMed
            56. Galperin MY: Bacterial signal transduction network in a genomic perspective. Environ Microbiol 2004, 6:552–567.View ArticlePubMed
            57. Ulrich LE, Zhulin IB: MiST: a microbial signal transduction database. Nucleic Acids Res 2007, 35:D386–390.View ArticlePubMed
            58. Galperin MY, Nikolskaya AN, Koonin EV: Novel domains of the prokaryotic two-component signal transduction systems. FEMS Microbiol Lett 2001, 203:11–21.View ArticlePubMed
            59. Nissien RM, Ytterberg AJ, Bogdanove AJ, van Wijk KJ, Beer SV: Analyses of the secretomes of Erwinia amylovora and selected hrp mutants reveal novel type III secreted proteins and an effect of HrpJ on extracellular harpin levels. Mol Plant Pathol 2007, 8:55–67.View Article
            60. Oh CS, Kin JY, Beer SV: The Hrp pathogenicity island of Erwinia amylovora and identification of three novel genes required for systemic infection. Mol Plant Pathol 2005, 6:125–138.View ArticlePubMed
            61. Merighi M, Majerczak DR, Stover EH, Coplin DL: The HrpX/Y two-component system activates hrpS expression, the first step in the regulatory cascade controlling the Hrp regulon in Pantoea stewartii subsp. stewartii. Mol Plant-Microbe Interact 2003, 16:238–248.View ArticlePubMed
            62. Merighi M, Majerczak DR, Coplin DL: A novel transcriptional autoregulatory loop enhances expression of the Pantoea stewartii subsp. stewartii Hrp Type III secretion system. FEMS Microbiol Lett 2005, 243:479–487.View ArticlePubMed
            63. Merighi M, Majerczak DR, Zianni M, Tessanne K, Coplin DL: Molecular characterization of Pantoea stewartii subsp. stewartii HrpY, a conserved response regulator of the hrp type III secretion system, and its interaction with the hrpS promoter. J Bacteriol 2006, 188:5089–5100.View ArticlePubMed
            64. Lan L, Deng X, Zhou J, Tang X: Genome-wide gene expression analysis of Pseudomonas syringae pv. tomato DC3000 reveals overlapping and distinct pathways regulated by hrpL and hrpRS. Mol Plant-Microbe Interact 2006, 19:976–987.View ArticlePubMed
            65. Kelm O, Kiecker C, Geider K, Bernhard F: Interaction of the regulator proteins RcsA and RcsB with the promoter of the operon for amylovoran biosynthesis in Erwinia amylovora. Mol Gen Genet 1997, 256:72–83.View ArticlePubMed
            66. Fraser GM, Hughes C: Swarming motility. Curr Opin Microbiol 1999, 2:630–635.View ArticlePubMed
            67. Harshey RM: Bacterial motility on a surface: many ways to a common goal. Annu Rev Microbiol 2003, 57:249–273.View ArticlePubMed
            68. Jarrell KF, McBride MJ: The surprisingly diverse ways that prokaryotes move. Nat Rev Microbiol 2008, 6:466–476.View ArticlePubMed
            69. Kaiser D: Bacterial swarming: a re-examination of cell movement patterns. Curr Biol 2007, 17:R561-R570.View ArticlePubMed
            70. Overhage J, Bains M, Brazas MD, Hancock RE: Swarming of Pseudomonas aeroginosa is a complex adaptation leading to increased production of virulence factors and antibiotic resistance. J Bacteriol 2008, 190:2671–2679.View ArticlePubMed
            71. McCarter LL: Regulation of flagella. Curr Opin Microbiol 2006, 9:180–186.View ArticlePubMed
            72. Wang Q, Suzuki A, Mariconda S, Porwollik S, Harshy RM: Sensing wetness: a new role for the bacterial flagellum. EMBO J 2005, 24:2034–2042.View ArticlePubMed
            73. Lai S, Tremblay J, Deziel E: Swarming motility: a multicellular behaviour conferring antimicrobial resistance. Environ Microbiol 2009,11(1):126–36.View ArticlePubMed
            74. Chen BG, Turner L, Berg HC: The wetting agent required for swarming in Salmonella enterica serovar Typhimurium is not a surfactant. J Bacteriol 2007, 189:8750–8753.View ArticlePubMed
            75. Torriani A: Influence of inorganic phosphate in the formation of phosphatases by Escherichia coli. Biochim Biophys Acta 1960, 38:460–479.View ArticlePubMed
            76. Sambrook J, Fritsch EF, Maniatis T: Molecular cloning: a laboratory manual Cold Spring Harbor, N.Y.: Cold Spring Harbor Laboratory Press 1989.
            77. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 1997, 25:3389–3402.View ArticlePubMed
            78. Krogh A, Larsson B, von Heijne G, Sonnhammer EL: Predicting transmembrane protein topology with a hidden Markov model: Application to complete genomes. J Mol Biol 2001, 305:567–580.View ArticlePubMed
            79. Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acid Res 1994, 22:4673–4680.View ArticlePubMed
            80. Galperin MY, Nikolskaya AN: Identification of sensory and signal-transducing domains in two-component signaling systems. Methods Enzymol 2007, 422:47–74.View ArticlePubMed
            81. Burse A, Weingart H, Ullrich MS: NorM, an Erwinia amylovora multidrug efflux pump involved in in vitro competition with other epiphytic bacteria. Appl Environ Microbiol 2004, 70:693–703.View ArticlePubMed
            82. Wang RF, Kushner SR: Construction of versatile low-copy-number vectors for cloning, sequencing and gene expression in Escherichia coli. Gene 1991, 100:195–199.View ArticlePubMed
            83. Valdivia RH, Falkow S: Fluorescence-based isolation of bacterial genes expressed within host cells. Science 1997, 277:2007–2011.View ArticlePubMed

            Copyright

            © Zhao et al. 2009

            This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://​creativecommons.​org/​licenses/​by/​2.​0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

            Advertisement