. $41.... r.‘ {47.3%}. :J. G. . h . .. ,iwmhhu Le... x I: i‘zu‘... . acts-3:3” .3 Fr... (3...! .1? Al. 553...}...3 1.... 1:11.", ‘ VF , 1.: I. . . 1" i. 3... 5.1.54. 2.» i“ & 3%». .3... a} in.) .3 S.‘......4.¢.. .. int: .. .{ .! Wag}. . s 3.37.. ,...n. t .h. a ‘ .. MI ’32 . 1'. )n .HHAHII... $.41“ .5‘ 1.9.. .LJ...‘ .75.: R. “(0:15! ‘ .Iihfl x b 9...: 3.7. 3, Rh. 7 x...u....,.a 1...: niggzo. 41.. S a 2.?! :3 «$5,... .. :1...» fix, iii]. .51) vii-1‘! 3!: 3231‘ gawaiywu. . _ 1 9. ... . Jill-(«41" s #3.?» a??? 3%. 3| . This is to certify that the dissertation entitled IDENTIFICATION AND ANALYSIS OF INDUCED GENES FROM ERWINIA AMYLOVORA AND MALUS X DOMESTICA DURING FIRE BLIGHT INFECTION presented by Sara E. Blumer—Schuette has been accepted towards fulfillment of the requirements for the Ph.D. degree in Plant Pathology W’— Major Professor’s Signature Dan. [5' 200C Date MSU is an Affirmative Action/Equal Opportunity Institution unq-—.—.-.-.—.—.-.—-.--.-.—-.-.-.-.-.-.-.-.--.-.-.----.-.-.-.g.-.---.--.-.-.-.--.-.-.-.-.-.-----.--.-v--.- LIBRARY Michigan State University PLACE IN RETURN BOX to remove this checkout from your record. To AVOID FINES return on or before date due. MAY BE RECALLED with earlier due date if requested. DATE DUE DATE DUE DATE DUE “@103 9099 2/05 pleIRC/DateDuejnddpJ IDENTIFICATION AND ANALYSIS OF INDUCED GENES FROM ER WINIA AMYLOVORA AND MAL US X DOMESTICA DURING FIRE BLIGHT INFECTION By Sara E. Blumer—Schuette A DISSERTATION Submitted to Michigan State University in partial fitlfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY Department of Plant Pathology 2006 ABSTRACT IDENTIFICATION AND ANALYSIS OF INDUCED GENES FROM ER WINIA AMYLOVORA AND MALUS X DOMESTICA DURING FIRE BLIGHT INFECTION By Sara E. Blumer-Schuette Fire blight, caused by the bacterial pathogen Erwim‘a amylovora, is a difficult disease to manage due to the virulent nature of the pathogen, disease susceptibility of most popular apple varieties, and the lack of known resistance genes in Malus spp. The goals of my research are to identify and characterize genes fi‘om E. amylovora that are induced during infection and to identify genes in Malus spp. that are associated with disease resistance. I utilized an in viva expression technology screen to identify in planta upregulated E. amylovora genes, two of which had homology to pseudopilins from a type II secretion operon and an endo-polygalacturonase (peh). Deletion mutants of both the peh gene and outhA were constructed to flirther analyze the contribution of type II secretion to the E. amylovora-Malus spp. interaction. This work resulted in the first report of a polygalacturonase enzyme influencing virulence of E. amylovora. Some apple cultivars, such as Red Delicious, exhibit tolerance to fire blight. Because resistance to fire blight appears to be a quantitative trait in apple, I chose to identify genes associated with resistance using suppressive subtractive hybridization. I generated 183 unique expressed sequence tags (ESTs) fiom cultivar Red Delicious during infection that are absent or repressed in the fire blight sensitive cultivar Gala. Temporal expression analysis identified 21 of the ESTs which are induced in Red Delicious of which eight show differential expression in other apple cultivars of varying resistance to fire blight. Further functional characterization of these genes will shed new light on signaling cascades that lead to a successful resistance response; in addition these genes can be converted into molecular markers for mapping onto linkage maps of apple for use in marker assisted selection of E. amylovora resistance in apple breeding projects. Copyright by SARA E. BLUMER-SCHUETTE 2006 To my husband Eric, who never stopped smiling. ACKNOWLEDGMENTS I first need to acknowledge my major advisor Dr. George Sundin for his effort, patience and support during my graduate studies. I appreciate and admire the continual optimism and adventurous attitude that George brings to the Sundin lab. I thank my committee members, Dr. Dennis Fulbright, Dr. Sheng Yang He and Dr. Steve vanNocker for their suggestions and input on my research and also for their comments and reviews for my dissertation. I would like to thank Dr. Youfu Zhao for collaborating with me on the first chapter of my dissertation, and for the discussions and help during my studies as well. I thank Gayle McGhee for kindly donated the gus-reporter strain, pGCMO used in chapters one and two and is always willing to talk about life, MUG and pears. I would like to thank Gail Ehret for always managing to find the apple trees we needed and for taking me out into the field. Even though the trips were as much about the Blue Bird as they were about work, I learned a great deal fiom Gail about plant pathology while out in the field. I also want to thank the other students from the Sundin lab past and present: Lisa Renick, Andrea Cogal, Lindsay Tripplett, Jessica Koczan, Mike Weigand, Erin Lizotte and Matt Berry for their camaraderie, research discussions and general mayhem when appropriate. None of this would have been possible if it wasn’t for my parents Ronald and Susan Blumer who have supported and encouraged me throughout my studies and have patiently waited ten years to see this day. I also am lucky to have supportive in-laws, Ernest and Sandra Schuette who opened their home to me. vi TABLE OF CONTENTS LIST OF FIGURES .................................................................................. ix LIST OF TABLES .................................................................................... x LITERATURE REVIEW ............................................................................ 1 Epidemiology of fire blight ..................................................................... 2 Molecular mechanisms of E. amylovora pathogenicity ..................................... 4 Virulence factors of E. amylovora ................................................................................. 6 Type H secretion in phytopathogenic bacteria ................................................ 8 Exoproteins secreted by type II secretion ................................................... 11 Host response to fire blight infection ......................................................... 13 Using plant genomics to complement genetic analysis of resistance loci ............... 16 Resistance gene analogues ..................................................................... 19 Candidate gene analysis with crop-specific RGAs ......................................... 20 Comparative genomics for R gene discovery ............................................... 23 Candidate gene analysis using crop-specific EST data .................................... 24 Synthesis of genomic and genetic data: Broad spectrum disease resistance. . . . . ....25 Concluding Remarks ........................................................................... 26 CHAPTER 1 Identification of Erwinia amylovora genes induced during infection of immature pear tissue .................................................................................................. 28 Abstract .......................................................................................... 29 Introduction ...................................................................................... 30 Materials and Methods ......................................................................... 34 Results ............................................................................................ 43 Discussion ....................................................................................... 64 CHAPTER 2 Taking a page from macerogenic bacteria: Functional characterization of an Erwinia amylovora endo-polygalacturonase and its effect on virulence .............................. 72 Abstract .......................................................................................... 73 Introduction ...................................................................................... 74 Materials and Methods ......................................................................... 78 Results ............................................................................................ 88 Discussion ...................................................................................... 1 05 CHAPTER 3 Genetic response of resistance-related genes to E. amylovora infection in Malus X domestica cultivars with varying levels of resistance ......................................... 111 Abstract ......................................................................................... 1 12 Materials and Methods ....................................................................... 117 vii Results .......................................................................................... 124 Discussion ...................................................................................... 141 APPENDIX ........................................................................................ 147 REFERENCES .................................................................................... 1 5 5 viii LIST OF FIGURES Figure 1: Overview of the IV ET screen for E. amylovora genes induced during infection of immature pear discs .............................................................................. 45 Figure 2: Symptoms and bacterial growth of Erwinia amylovora WT strains and hrpA and dspE mutants in immature pear ............................................................... 48 Figure 3: Symptoms and growth of Erwinz’a amylovora WT Ea1189 and corresponding hothngA and mltEEA mutants in immature pear ............................................... 64 Figure 4: Alignment of type H secretion systems from Y. enterocolitica (ytsl) and E. amylovora (outga) ................................................................................... 89 Figure 5: Areas of necrosis on inoculated immature pears ..................................... 95 Figure 6: Ruthenium red stain for polygalacturonase activity ............................... 102 Figure 7: Western blot ofEa1189(pSEB31) subcellular fractions .......................... 102 Figure 8: Supernatant proteins separated out on SDS-polyacrylamide gel ................ 104 Figure 9: Representation of the suppression subtractive hybridization library construction ......................................................................................... 119 Figure 10: Representative radiographs of differentially screened EST clones ............ 126 Figure 11: Representative cDNA nylon macroarray hybridization ......................... 129 Figure 12: Semi-quantitative RT-PCR with cultivars Empire, Fuji and Red Delicious ................................................. i .......................................... 138 Figure 13: PCR amplifying gene fragments corresponding to the EST of the 14-3-3 and unknown gene with genomic DNA of Empire, Fuji, Gala and Red Delicious ............ 140 ix LIST OF TABLES TABLE 1. Bacterial strains, plasmids and primers used in this study ........................ 35 TABLE 2. Selected list of Erwim'a amylovora genes induced during infection of immature pear tissue ................................................................................ 52 TABLE 3. Expression of NET clones after inoculation of immature pear fi'uit ........... 63 TABLE 4. Bacterial strains, plasmids and primers ............................................. 80 TABLE 5. Protein homology of E. amylovora Out operon .................................... 91 TABLE 6. Fold induction of promoter-uidA fusion constructs inoculated in minimal medium ............................................................................................... 97 TABLE 7. Fold induction of promoters upon infection of apple cultivar ‘Gala’ ......... 100 TABLE 8. Oligonucleotide primers used in this study ....................................... 123 TABLE 9. Classification of ESTs recovered from an SSH experiment subtracting cDNA of apple cultivar Gala fi'om Red Delicious following inoculation with E. amylovora.... 125 TABLE 10. Red Delicious ESTs with two-fold or higher induction after inoculation with E. amylovora ........................................................................................ 128 TABLE 11. ESTs that are induced to higher levels in apple cultivar Red Delicious versus Empire following inoculation with E. amylovora .............................................. 132 TABLE 12. ESTs that are induced to higher levels in apple cultivar Red Delicious versus Fuji following inoculation with E. amylovora ................................................. 133 TABLE 13. Constitutively expressed ESTs of apple cultivar Red Delicious and analysis of expression of these genes in the cultivars Empire and Fuji ................................ 136 TABLE 14. Complete list of sequenced BSTs fi'om an SSH library constructed by subtracting cDNA pools from the apple cultivars Red Delicious and Gala after inoculation with E. amylovora ................................................................... 147 LITERATURE REVIEW Fire blight is a necrogenic disease of rosaceous plants caused by the gram negative bacterium, Erwinia amylovora. In the United States, fire blight can be a devastating disease on non-native hosts such as domestic apple, pear and quince due to the lack of resistance in popular cultivars and limited control options. Due to the limited amount of control options available, there is interest in the molecular mechanisms of fire blight in order to generate novel control measures. Initial genetic research of the pathogen E. amylovora has identified major pathogenicity factors such as the secreted effector protein DspE, a type III secretion system and production of the exopolysaccharide amylovoran (18, 21, 26, 82, 127). The discovery of these pathogenicity factors relied mainly on the creation and mapping of non-pathogenic mutants using phage insertional mutagenesis (239). Research investigating the host response to fire blight infection has demonstrated the induction of known pathogenesis related (PR), antioxidant and phenylpropanoid pathway genes (31, 35, 153, 248, 249). Currently, no resistance genes for fire blight have been identified; however absolute resistance does exist in Malus species related to domestic apple and certain cultivars of domestic apple show elevated levels of resistance to fire blight. Due to the increased ease and availability of genomic and genetic techniques, I chose to look at the fire blight pathosystem on a larger genomic scale in order to identify additional virulence factors from E. amylovora and resistance response genes from Malus X domestica cultivar Red Delicious. Epidemiology and control of fire blight Fire blight was first described on rosaceous fruit trees in New York around 1780, and is thought to be indigenous to North America. The disease moved westward along with the movement of European settlers bringing nursery stock from the eastern United States, until fire blight was present throughout the entire United States and part of Canada. As apple and pear nursery stock were exported from the United States, fire blight spread worldwide to Europe, the Middle East, and New Zealand (32). Fire blight is disseminated in orchards by either insects or wind and by the use of contaminated tools and epidemics can develop when both humidity and temperature are high during bloom. In late spring to early summer, over-wintering cankers will start to exude bacterial ooze produced by E. amylovora that can be spread to flowering trees by pollinating insects such as bees or by wind and or rain. Once E. amylovora is transferred to a blossom, the bacteria can survive briefly as epiphytes, multiplying on the surface of stigmas (241). Afier this brief epiphytic phase, E. amylovora enters the plant through nectarthodes located in the floral cup. The blossom typically becomes necrotic at this point, and this mode of infection is referred to as blossom blight. Bacteria that enter through the nectarthodes may also infect immature fruit that develop from an infected blossom, creating secondary sources of inoculum as bacterial ooze exudes from the infected fruit (119). Again, wind and rain can carry droplets of ooze that contain embedded E. amylovora to wounds on stems and leaves created by wind or mechanical damage. Infection through the stems and leaves results in the movement of E. amylovora to the vascular tissue, where surrounding parenchyma tissue collapses and disrupts water and solute transport. The disruption of the vascular system contributes to the formation of shoot blight which is typified by a shepard’s crook that quickly becomes necrotic and also may exude bacterial ooze (241). Movement of E. amylovora through the vascular system can continue down into the rootstock, killing the vascular system in a susceptible rootstock which will eventually kill the entire tree, this is commonly referred to as rootstock blight (168). Shoot blight will typically form cankers in woody stem tissue that are the primary source of over-wintering inoculum mentioned above. There are also cases of E. amylovora surviving in budwood that could also act as primary sources of inoculum (241). Control of fire blight is often through the use of cultural practices, copper sprays and antibiotics. Cultural practices are comprised of the selection of resistant rootstocks and scions, reduction of inoculum by pruning of diseased tissue and cankers from the orchard and sanitation of tools used for pruning. In respect to chemical control, the application of copper sprays can only be used before leaf set, due to phytotoxicity, which limits copper sprays to only controlling primary E. amylovora populations in the early summer (182). Streptomycin is the most commonly used antibiotic for fire blight control, which is used to control blossom blight and to reduce inoculum levels after violent rain storms that cause tissue damage (157). Development of streptomycin resistance has reduced the effectiveness of streptomycin applications for controlling fire blight. Streptomycin resistance in E. amylovora populations has developed through two independent means. The first report of streptomycin resistant E. amylovora in California was due to a single point mutation in the streptomycin binding site of ribosomal protein 812 (48). In addition to the isolation of ribosomal point mutations, some streptomycin resistant E. amylovora populations in Califomia carry plasmid pEa8.7 that contains the streptomycin resistance genes strA-strB. Streptomycin resistance in Michigan however has arisen due to the horizontal transfer of plasmid pEA34 which harbors the transposon Tn5393 that expresses the streptomycin resistance genes strA-strB under the control of a promoter located in the insertion element 181133 (49). The transposon Tn5393 has also jumped to the ubiquitous plasmid pEa29 and to the chromosome of some field isolates and as a result, fewer streptomycin resistant E. amylovora are found that harbor the plasmid pEa34 (155, 158). Molecular mechanisms of E. amylovora pathogenicity As mentioned before, three major pathogenicity factors have been identified for E. amylovora. First is the capability to secrete effector proteins using the type III secretion system. Type III secretion utilizes a needle-like structure to translocate proteins across the bacterial inner and outer membranes into the host cell. This secretion system is highly conserved among both animal and plant pathogenic bacteria and is independent of Sec-mediated protein secretion (39). In E. amylovora the type III secretion system is organized into hip (hypersensitive response and pathogencity) operons (27, 127, 253). Regulation of the hrp operons ensures that the type III secretion machinery is only expressed during conditions that mimic the apoplast of the host such as low nutrient content and low pH. The two-component sensor hrpXY locus is responsible for the first known signaling step for hip expression. The sensor protein, HrpX is theorized to be anchored to the inner membrane and responds to favorable hrp-inducing stimuli by phosphorylating Her. Activated response regulator protein Her contains a helix-tum- helix motif that is then able to bind upstream of the hrpL gene to induce transcription along with sigma factor 54 and another enhancer protein HrpS (252). HrpL is the main activator of transcription of the structural hrp operons and effector protein genes that contain hrp boxes upstream of the gene (125). Interestingly enough, recent evidence suggests that HrpL downregulates flagellar biosynthesis due to the observed increased motility and flagella formation in hrpL mutants (44) All bacteria that posses functional type III secretion systems share a core set of nine conserved hrp genes that were renamed as hrc (hypersensitive response and gonserved) genes. This includes hrcC which encodes the outer membrane pore forming protein, hch which encodes a lipoprotein, five hrc genes encoding inner membrane proteins, an ATPase homologue hrcN and the gene hch that encodes a protein that may be secreted (125). Another hm gene that are known to form the type III secretion apparatus in E. amylovora is hipA which was demonstrated to encode the pilin that forms the type 1H secretion pilus in E. amylovora, mutants in hrpA were non-pathogenic and known type III secreted proteins remained in the cytoplasm (118). Besides the type III secretion system, a secreted effector protein, DspE is the second pathogenicity factor of E. amylovora. Originally, dspE mutants were observed as disease specific mutations, with no effect on the HR when the mutants were inoculated on tobacco (193). It was later demonstrated that dspE mutants could elicit an HR on tobacco especially in highly virulent strains of E. amylovora and also on soybean leaves. In addition, the transfer of the dspEF locus to Pseudomonas syringae pv. glycinea rendered it avirulent on soybean (26). DspE was demonstrated to rely on its chaperone DspF for secretion and that secretion was of DspE used type III secretion (25, 82). Later studies determined that the chaperone DspF interacts with and protects DspE from degradation (83). Genetic analysis determined that the upstream sequence of dspE contains a hrp box, and expression of dspE is HrpL-dependant and coordinated with expression of the neighboring hrp operons (82). Amylovoran, the main exopolysaccharide produced by E. amylovora is the third major pathogenicity factor. Biosynthetic genes for amylovoran production are organized into the ams operon that consists of 12 genes. Regulation of the ams operon is by a two- component regulatory system RcsB/C (2, 19, 123), and the global regulator protein H-NS which inhibits the production of both amylovoran and levan (103). Amylovoran is comprised mainly of a repeating galactose backbone with galactose and glucuronic acid side chains (84). The composition of sugar monomers in amylovoran can change amongst E. amylovora isolates from differing host tissues, such as E. amylovora isolates from Rubus spp. that are unable to infect apple or pear (148). Mutants in key genes of the amylovoran biosynthetic operon were non-pathogenic and unable to move or increase cell numbers in planta (18, 28). This exopolysaccharide is thought to obstruct the host vascular tissue, have a protective role for E. amylovora from oxidative burst and desiccation and also a possible role for masking outer membrane proteins from the host (84,128,148) Virulence factors of E. amylovora Virulence factors also contribute to infection by E. amylovora including additional effector proteins HrpN and AerptZEA (253, 268), the iron-binding siderophore desferrioxamine (60), and genes required for metabolism of sorbitol (3), sucrose (29) and galactose (161). The effector protein HrpN was the first type III secreted protein from E. amylovora to be described and insertional mutants were non- pathogenic on immature pear fruit (253). Irnmunogold labelling of HrpN protein demonstrated that HrpN is secreted into the host apoplast, but not delivered into the host cells (193). Recently, HrpN3937 from Erwinia chrysanthemi was determined to act as an aggregation factor and is required for the successful formation of the pellicle (262). It remains to be determined if HrpN from E. amylovora is involved in the formation of a biofilm in planta. Recently AerptZEA, an additional effector protein contributing to virulence was described by Zhao et al. (267). Deletion mutagenesis of aerptZEA affected virulence and bacterial cell count in immature pear fruits. Additionally, the transfer of aerptZEA to Pseudomonas syringae pv. tomato DC3000 increased its ability to infect Arabidopsis rpsZ mutants, indicating that AerptZEA is recognized by RPS2 (267). The ability to scavenge for iron in nutrient-lirniting environments is also required for E. amylovora virulence under certain conditions. The siderophore desferrioxamine is secreted into the extracellular milieu to complex iron for uptake by the ferrioxamine receptor FoxR. Siderophores are also important for protection of bacterial cells against the host oxidative burst (73). Insertional mutants of desferrioxamine dfoA and ferrioxamine receptor foxR were hindered in growth under iron limiting conditions. In addition, both mutants were reduced in their ability to colonize and infect apple blossoms (60). Erwim'a amylovora also secretes other virulence factors that affect its interaction with host plants. A metalloprotease (PrtA) which is accompanied by its own type I secretion machinery was identified in the supernatant fi'action from a highly proteolytic strain of E. amylovora (E8). Although this metalloprotease did not contribute greatly to virulence of E. amylovora, the ability of prtA mutants to escape out of xylem vessels and into xylem parenchyma and movement through parenchyma was reduced (265). Another component of E. amylovora exopolysaccharide is levan, a 6-2, 6-fructan produced by the secreted enzyme levansucrase. Upon secretion to the extracellular milieu, levansucrase produces levan from extracellular sucrose and releases glucose as a by-product. Mutants that do not produce levan are not significantly altered in virulence, although one mutant exhibited reduced virulence under high sucrose conditions (20, 85). Levansucrase is not modified in any manner for its secretion and a putative porin downstream of the lsc gene is thought to have a role in its secretion (69). Type II secretion in phytopathogenic bacteria E. amylovora also possesses a full sec transport system which supports other secretion systems such as type II, type IV and type V (auto-transporter) secretion (188, 238, 266). Type II secretion (T28) is an important virulence factor in Gram-negative animal pathogens and phytopatho gens such as Pectobacterium carotovorum, E. chrysanthemi, Ralstonia solanacearum, and Xanthomonas campestris. The T28 operon was identified first from Klebsiella oxytoca by transferring a plasmid encoding the T28 operon from K. oxytoca to Escherichia coli to demonstrate that T28 proteins were required for the secretion of the enzyme pullulanase (215). Transfer of a cosmid containing the T28 from E. chrysanthemi to E. coli was also used to demonstrate the ability of the T 28 to secrete various extracellular enzymes (138). Genome sequencing has made it possible to identify the presence of T28 genes in many gram negative bacteria other than pathogenic bacteria (50). However, the comparison of known functional T28 systems highlights 12 to 15 proteins expressed from the operon are required for secretion across the outer membrane (215, 219, 220). Although the gene naming convention has differed slightly amongst species possessing a T28 system, the majority of the gene products follow the Klebsiella gene designation as A through 0 plus 8 (50, 219, 220), with proteins CDEFGHIJKLMNO acting as the core of 12 essential T28 components (50). Here the general nomenclature of the T28 operon, “general secretion pathway” (Gsp) will be followed Overall, the core T28 proteins are used to facilitate the transport of folded proteins from the periplasm across the outer membrane. The T28 protein GspD is a member of the secretin protein family, which includes proteins fiom T388 and type IV pilus generation. This group of related proteins function to form pores in the outer membrane for protein secretion (238). In some cases where a GspS homologue is present in the T28 operon, an extension of the GspD C-terminus is required for GspS interaction, stabilizing GspD multirners in the outer membrane (54, 183, 228). Electron microscopy analysis of K. oxytoca GspD (PulD) multirners with and without Gsp8 (PulS) show a stacked complex with an apparent occlusion in the center of the complex, which is hypothesized to be the N-terminal of GspD (183, 184). This evidence, combined with the observed channel activity of multimeric GspD/ GspS complexes after the application of an electrical current now indicates that GspD may act as a voltage-induced gated channel rather than a passive pore for protein transport (183). There are examples, however, of species that do not posses a GspS homologue in their T28 operon, such as Pseudomonas aeruginosa and X. campestris, where their GspD homologue is still able to localize to the outer membrane (228), however the absence of a GspS homologue does not mean that a pilot protein coordinating the insertion of GspD multimers into the outer membrane does not exist (219). A likely source for energy in the secreton is GspE, a T28 protein theorized to be localized in the cytoplasm and demonstrated to act as an ATPase in Vibrio cholerae (43). GspE has been demonstrated to interact with the inner membrane bound GspL in V. cholerae (EpsL) and E. chrysanthemi (OutLEch) (221), and that conformational changes occur with both OutEEch and OutLEc}l in E. Chrysanthemi (198). GspL has also been shown to interact with another inner membrane bound protein, GspM, through immunoprecipitation experiments and GspE-GspL-GspM tripartite complexes are thought to involved in the assembly of the entire secreton or just the transport and assembly of a T28 pilus that would bring T28 secreted proteins to the GspD/GspS complex for secretion (43, 198). Other critical components of T28 are inner membrane bound with some degree of periplasmic domains, such as the major pseudopilin, GspG and minor pseudopilins GspH-J, named for their homology to type IV pilins. Pseudopilins are processed by GspO and are capable of forming protein-protein interactions amongst each other in vitro (67, 145). GspG homologues are able to form pilus-like structures when overexpressed (70, 106, 223); this has led to speculations of pseudopilins forming a pilus that polymerizes to push secreted enzymes through the GspD pore (220). 10 Exoproteins secreted by type II secretion Specificity of exoproteins secreted through the T28 has been observed in K. oxytoca, P. carotovorum and E. chrysanthemi. In all three systems, mutations of GspC or GspD could not be complemented by homologous proteins from other systems (173, 197, 228). This led to the hypothesis that they determine the specificity of secreted proteins. The N—terminus of OutDEch from E. chrysanthemi was determined through a protein deletion series to interact with secreted endogenous proteins, and was not able to interact with T28 secreted proteins from P. carotovorum (228). In the system of K. oxytoca, however, the N-terminal of OutDEch from E. chrysanthemi fused to a truncated PulD protein was still able to secrete pullulanase, indicating that the requirements for secretion specificity may differ fiom species to species (95). Therefore, the contribution of GspC to secretion specificity must not be overlooked. GspC and its homologues are inner membrane bound proteins with large periplasm-associated domains (34). Based on amino acid sequence analysis of GspC homologues, a PDZ motif is present near the C-terminal of GspC homologues from K. oxytoca (PulC), P. carotovorum (OutCEca) and E. chrysanthemi (OutCEch). Certain mutations of the PDZ motif of PulC with the addition or deletion of amino acids closest to the C-teminus rendered the T28 secreton non-fimctional (196). When domains from OutCEca and OutCEch were swapped to form fusion proteins, the PDZ motif was shown to be involved in the secretion of species specific exoproteins in E. chrysanthemi (34). The search for a type II secretion signal in exoproteins has had some success, with regions responsible for secretion across the outer membrane determined for pullulanase 11 of K. oxytoca, PelC of E. chrysanthemi and PehA of P. carotovorum (76, 137, 190). Recent work with PehA from P. carotovorum indicates that the secretion signal needed for species dependant secretion is a three-dimensional motif (1 89). No easily identifiable T28 signal has been determined by amino acid sequence comparisons as of yet. Exoproteins secreted by phytopathogenic bacteria include proteinases, amylase, pectinases, cellulases and polygalacturonases (50). Polygalacturonases present in phytopathogenic bacteria fall into two categories, endo- and exo-polygalacturonases. The former catalyze the random hydrolysis of pectic acid forming oligogalacturonates while the latter catalyze the terminal hydrolysis of pectic acid forming galacturonic acid monomers (115). While endo-polygalacturonase activity has not been detected in E. chrysanthemi (109), endo-polygalacturonase activity is a known virulence factor in P. carotovorum (216), R. solanacearum (61, 107) and Agrobacterium vitis (210). Analysis of the secretome of E. chrysanthemi identified 14 T28 secreted proteins, including 11 pectinases, cellulase, rhamnogalacturonan lyase, a novel esterase and a novel Avr-like protein (122). Sequence analysis of Ralstonia solanacearum has identified 6 exoproteins including three polygalacturonases, a pectin methylesterase, and two cellulases (139). Previously, E. amylovora was described not to possess any cell wall degradation abilities (226); however, recent studies have described the presence of a type II secretion gene, outF, based on Southern hybridization to an outFEca homologue and weak CelA activity (207). 12 Host response to fire blight infection Interest in the genetic mechanisms of the host response to fire blight infection has increased due to the limited control options available to growers. Venisse et al. demonstrated that wild type E. amylovora and an ams mutant were capable of inducing the oxidative burst in pear and a susceptible apple cultivar’s leaves, whereas the type HI secretion mutant was unable to induce an oxidative burst (248, 249). Furthermore, in another study, Venisse et al. were able to link the secretion of two effector proteins, HrpN and DspE to the elicitation of the oxidative burst of pear (247). When the levels of expression of phenylpropanoid genes were observed in a resistant and susceptible cultivar of apple, it appeared that effector proteins were also down-regulating some of the flavanol biosynthesis genes in the susceptible cultivar (249). Based on these studies, the elicitation of an oxidative burst by the host appears to aid infection by E. amylovora, and the lack of down-regulation of flavanol biosynthesis genes appears to be related to resistance (248, 249). Other studies concentrated on pathogenesis-related (PR) genes and their expression after chemical elicitation or infection in apple. Use of a salicylic acid analogue, acibenzolar-S-methyl (ASM) was demonstrated to induce PR-I, PR-2, PR-8 in cultivar Jonathan, PR-Z and PR-IO isoforrns in cultivar Golden Delicious (35, 153, 269). When apple cultivar Gala was treated with ASM, there was no noticeable induction of PR—Ia, PR-2, PR-5 or PR-8, the authors explain that this may be because in the previous studies, apple seedlings were used compared to the use of shoots on woody tissue for the Gala experiment. The induction of PR-Z, PR-5 and PR-8 were the strongest however after inoculation with E. amylovora 48 and 96 hours post inoculation (31). 13 Plant proteins interacting with the effector protein DspE have been identified fi'om Malus X domestica using a yeast two-hybrid system. Four of the proteins, DIPM 1 to 4 are serine threonine kinases with a leucine rich repeat (LRR) motif (159). This is similar to other receptor-like kinase (RLK)-LRR resistance (R) genes such as Xa21 fi'om rice, rng from barley and FLSZ from Arabidopsis (37, 55, 134). The site of interaction between DIPMs and DspE is the intercellular RLK domain, in contrast to other RLK- LRR family members where the interacting site is the extracellular LRR domain. The RLK-LRR DIPMs were present in all hosts of fire blight tested, both resistant and susceptible, alluding that the DIPM proteins are not R genes that lead to a successful resistance response (159). Cloning and mapping of RGAs fi'om apple species is another way to identify resistance-related genes. Classic resistance genes encode proteins with a nucleotide binding site (NBS) and a LRR domain that participates in protein-protein interactions (13). These resistance genes are now hypothesized to be guard proteins that monitor changes to other proteins that are the true targets of pathogen effectors (181). Resistance gene analogues from apple in two studies were cloned from sequences amplified based on degenerate primers in the nucleotide binding site (NBS) (14, 135). In one study, after the cloning of 27 sequences from RGAs, 18 of those sequences were successfully converted into molecular markers and mapped the locations of the RGAs on linkage maps (14). The technique NBS-profiling was proven useful by identifying and mapping of 23 RGAs at the same time in F1 progeny, adding to the number of RGA markers available for mapping (42). In a recent study, expressed sequence tags (ESTs) 14 homologous to RGAs were converted into expressed sequence tagged site (E-STS) markers, facilitating the mapping of ESTs onto linkage maps (175). Groups have also sought to map quantitative trait loci (QTLs) in apple and pear for fire blight resistance, to aid in breeding efforts. The first QTLs identified for fire blight resistance were located in a resistant cultivar of European pear. Four QTLs were identified, two on linkage group two, one on linkage group four and another on linkage group nine. Two molecular makers, one co-localizing with the QTL on linkage group two-b and the other co-localizing with linkage group four are resistance gene analogue markers, indicating that nucleotide binding site (NBS) type resistance genes may be involved in fire blight resistance (65). Calenge et al. were able to identify a major QTL in apple that accounts for 34.3 to 46.6 % of the resistance phenotype. This QTL was located on linkage group seven, with no RGA markers co-localizing with the QTL (41). This indicates a high probability that the main genetic component of fire blight resistance is not R gene mediated. Khan et al. also identified a major QTL for fire blight resistance on linkage group seven that accounted for 37.5% to 38.6% of the resistance phenotype. Both QTLs identified on linkage group seven were in roughly the same area and were observed using two different strains of E. amylovora for infected, therefore this QTL is likely to be a major fire blight resistance locus (41, 124). Furthermore, areas of the pear QTLs on linkage groups two and nine were determined to be involved in digenic interactions with other linkage groups, thus in apple, the contribution of the pear QTLs are in combination with other loci (41, 65). Although a major QTL for fire blight resistance has been determined, the exact nature of that resistance is still unknown. Since no currently cloned RGAs are located in 15 that region of linkage group seven, there are no mapped candidate genes for resistance in apple. The identification of the DspE interacting apple proteins is an insight into the potential innate immunity response of Rosaceous plants, however RLK-LRR proteins are unlikely to be involved in resistance to fire blight, either. In order to identify other genes involved in the resistance response to fire blight, attention should be given to genes other than the traditional R gene types. Using plant genomics to complement genetic analysis of resistance loci Major food crops such as rice, wheat, barley, maize, potatoes and soybeans still suffer losses due to plant pathogens on a global scale, despite the best efforts of breeding programs and the availability of pesticides (108, 235). The search for durable resistance is not a novel concept in plant biology and the overall goal of disease resistance breeding programs is to identify chromosomal regions conferring this type of resistance (121, 163). Identification of useful disease resistance QTLs are further complicated by rapid evolution of some pathogens, by overcoming released resistant cultivars in a matter of years (163, 233). In the search for the ultimate resistance gene, a great wealth of genomic material pertaining to plant pathogen interactions has been generated for pathogens and plants in model and crop species. Currently, there are four complete plant genome sequences available. The first plant genome sequence to be released was the model plant Arabidopsis thaliana in 2000 (110), followed by the genome sequence of two rice types (Oryza sativa L ssp. indica and 0. sativa L ssp. japonica) in 2002 (89, 263) and most recently, the genome of the model tree Papulus trichocarpa (245). These genome sequences are invaluable for establishing l6 areas of synteny with other related plants and using the genome sequences for molecular marker development and choosing candidate genes from the collinear area. With increased ability to obtain large scale expression profiling of pathogen-induced genes, more studies are also seeking to identify the global gene expression response of crop species to various pathogens. Expressed sequence tag libraries are available for rice infected with rice blast (Magnaporthe grisea), potato leaves infected with late blight (Phytophthora infestans), wheat infected with Fusarium head blight (F usarium graminearum), apple infected with apple scab (Venturia inaequalis) and nodule tissues of the legume model species Lotus japonicus (10, 114, 178, 213). In addition, there are a great deal of subtractive libraries produced from infected tissue, aimed at identifying ESTs in infected tissue that are only involved in the resistance response in: rice (144, 258), wheat (94, 129, 146), potato (242), barley (177), tomato (87), apple (58) and coffee (74). Resistance (R) genes were first described by Flor in his seminal work that established the gene for gene theory where an incompatible response (disease resistance) occurs when an R gene “recognizes” an avirulence gene fi'om a pathogen (reviewed in 55). More recently, a “guard hypothesis” has been proposed where the R genes are monitoring the status of a particular protein, and upon a change to that protein mediated by a pathogen effector protein, the R gene initiates a signaling cascade (55). The identification of a large number of R genes involved in gene for gene interactions has allowed for the classification of types of R genes and their conserved motifs. The most common R genes were those with a nucleotide-binding site (NBS) and leucine rich repeats (LRR) that are split into two types, those with coiled coil (CC) 17 domains at their N-terminal or those with the animal-like 1011 and interleukin receptors (TIR). The conservation of motifs in the NBS-LRR regions of these R genes has allowed researchers to clone resistance gene analogues (RGAs) from numerous plant species (156). In addition to the identification of R genes, characterization of disease response genes in the model system, A. thaliana, is useful for a reference point in which to identify similar mechanisms in crop species. The majority of the disease-related (DR) genes from Arabidopsis are involved in the signaling pathways that coordinate defense and were identified by mutant analysis (97). Defense signaling pathways are classified into defense against biotrophs or necrotrophs for simplicity, although the true picture of how a plant perceives a pathogen is assuredly more complicated (88). Very briefly, signaling cascades for defense against biotrophic pathogens typically occur through CC-NBS-LRR or TIR-NBS-LRR R genes that initiate an oxidative burst, activate MAP kinases, induce salicylic acid production and the induction of PR gene expression through NPR] (97). Detection of a necrotroph will initiate jasmonic acid and ethylene responsive pathways that induce gene expression through ethylene and jasmonic acid induced transcription factor ERFl, and jasmonic acid induced transcription factors RAP2.6 and JIN 1. In addition, host detection of a necrotroph will also activate a MAP kinase pathway that leads to the downregulation of the salicylic response (88, 97). As useful as genetic and genomic data are for assessing the molecular pathways and interactions leading up to a resistance response, the question of how that information can be applied is not often answered. Bridging gaps in between basic lmowledge of plant pathogen interactions and applied knowledge of resistance phenotypes in available 18 gerrnplasm is vital to identifying genes that coordinate durable resistance. In this literature review I will examine the coordination of large scale genomics projects and resistance phenotype analysis and the promise of these methods in yielding crops with durable resistance. Resistance gene analogues Approximately 149 NBS-LRR gene homologues were identified from Arabidopsis after the genome sequence was released. A physical map of the locations of these RGAs shows areas of 43 RGA clusters on the chromosomes, the result of localized gene duplication (162). While evolutionary insights leading to the duplication and distribution of RGAs are critical for the understanding of the development and maintenance of disease resistance in plants, the function of the majority of RGAs is unknown; function has only been determined for twelve NBS-LRR genes (105). As such, once an RGA is found to contribute to disease resistance, it is then known as an R gene. The mapping of RGAs to areas of QTL localization for disease resistance could help to identify the potential function of these genes. Mapping of RGAs to QTL regions has also been conducted in major crops such as rice, soybean, potato, rapeseed, and barley (reviewed in 195). RGAs are also currently being used as candidate genes in disease resistance QTL analysis in rice (202, 255, 256), potato (108, 233), and crucifers (231,236) 19 Candidate gene analysis with crop-specific RGAs Candidate gene analysis involves the identification of previously characterized defense-related genes, such as those identified already in Arabidopsis (13, 55), or the identification of RGAs before disease resistance QTL analysis. The candidate genes are then converted into molecular markers for use in the QTL analysis and for placement onto the genetic map. Those candidate genes that co-localize with the disease resistance QTLs can be firrther analyzed without the need for map based cloning (reviewed in 195). In addition to the potential for an RGA to be identified as a ftmctional R gene, RGAs may also be used as molecular markers for other resistance loci, since RGA clusters have mapped to resistance loci as is the case with Arabidopsis (l, 234). Therefore, caution must be exercised on assumptions from co-localizing candidate genes and the function of the candidate gene and phenotype must be confirmed using other methods such as Northern analysis or the use of transgenic plants for phenotype analysis. In the case of the Rng resistance gene to stem rust in barley, RGAs mapped near Rng, but attempts to use synteny with rice failed to identify any candidate genes in that area (12) and a non-traditional R gene, Rng was later cloned using positional cloning and was confirmed as being absent in the rice genome (37). RGAs were used in the identification and cloning of a resistance gene RB from Solanum bulbocasta after other R genes that were bred into domestic potato (Solanum tuberosum) failed (233). Cloning of the RB resistance gene from S. bulbocasta used both traditional map-based cloning techniques, along with long-range PCR targeting the identified candidate gene RGAs from genomic DNA of resistant S. bulbocasta. The BAC carrying the RB allele was unstable in E. coli; therefore the complete RB gene was 20 amplified from genomic DNA and cloned into a binary vector for plant transformation and validation as the gene conferring resistance in previous QTL analysis (174). Insertion of RB as a transgene into a late blight-susceptible potato variety conferred resistance to all races of P. infestans, including one race that had overcome all previous R genes (233). In maize, a combination of RGA mapping and transposon mutational analysis was used to clone a resistance gene from the 1p] locus. RGAs had previously been mapped to the maize physical map and one RGA, PIC20 had mapped to the rpI locus which is responsible for maize resistance to common rust (Puccinia sorghi). The detection of an insertion into the R gene RpI-D was detected by the PIC20 probe, and led to the subsequent cloning and characterization of RpI-D (53). Candidate gene analysis in rice has been useful for the identification of rice RGAs and defense-related genes that co-localize with QTLs for bacterial blight resistance (Xanthomonas oryzae pv. oryzae) and rice blast (Magnaprothe grisea). One study compiled a group of 118 molecular markers based on RGAs and previously determined defense-related (DR) genes (202). In total, six of the candidate gene markers co- localized with QTLs for bacterial blight resistance: five were RGAs and one was a defense-related gene, oxalate oxidase. A higher number of candidate genes with functional characterization co-localized with previously determined blast resistance QTLs (202). Due to the low number of defense-related markers associating with the bacterial blast resistance QTLs, caution must be exercised with choosing potential candidate genes for screening, the higher number of DR markers that associated with the fungal (blast) 21 resistance QTLs suggests that the authors selected genes that were more likely to be involved in fungal, not bacterial disease resistance in rice. A second candidate gene approach to QTL analysis determined which RGAs and defense-related genes co-localized with rice blast resistance QTL in advanced backcross populations of the rice cultivar Vandana (256). Advanced backcross populations were used so that lines with promising resistance already possessed the positive agronomic traits from the susceptible parent, and could be released with minimal further selection and breeding (237). Seven markers were associated with blast resistance when a single blast isolate was used for infection; including two candidate genes and one RGA. A higher number of candidate genes were associated with resistance when multiple blast isolates were used for inoculation; there were four candidate genes and two RGAs. In addition, one of the candidate genes, a NBS-LRR gene, was also determined to have a digenic effect with two other NBS-LRR genes and the RGA homologue (256). The increased number of candidate genes identified as being significant to resistance during multiple isolate pressures is interesting, since this more closely mimics natural field conditions that new cultivars face. Particularly since there were digenic effects detected between NBS-LRR genes, supporting the hypothesis that R gene diversity is driven in part, by pathogen diversity (105). Overall, the localization of the RGAs with rice resistance QTLs is interesting; however it requires further characterization of the RGAs to determine their involvement in the resistance response. 22 Comparative genomics for R gene discovery Comparative genomics is also a potentially useful method to identify collinear resistance genes between two related species where the fimction of the R gene is known in one of the species. The identification of an R gene in tomato; the model Solanaceous plant, against Fusarium wilt was used to identify the potato homologue based on extensive synteny between the genetic maps of tomato and potato. By aligning the genetic maps of the 12 resistance locus from tomato and the R3 resistance locus from potato, the identification of the collinear R gene, R3a was determined in potato (108). This was the first known report of an R gene being cloned based on synteny. Since Arabidopsis is also a crucifer, RGAs derived from the Arabidopsis genome are of interest to groups working with agriculturally important crucifers. Because no R genes have been cloned from Brassica napus or its progenitors, Brassica rapa and Brassica oleracea, previously identified RGAs from Arabidopsis (33) were selected as R- EST clones and mapped onto a B. napus genetic linkage map in order to determine areas of collinearity, which were found to exhibit conserved gene content and order (231). Similarly, for B. rapa, markers based on known open reading frames (ORFs) from Arabidopsis were used to map areas of collinearity between clubroot (Plasmodiophora brassicae) resistance QTLs and the corresponding segment of an A. thaliana chromosome that is known to have R gene clusters (236). Overall for B. napus, the inheritance of R genes on chromosomal segments from the Brassica progenitors have contributed to R gene duplication, however these areas are still collinear to Arabidopsis. Once areas are identified in B. napus as being related to resistance and in the case of B. 23 rapa QTLs, researchers can take advantage of the genomic data from Arabidopsis to identify candidate genes (231, 236). The areas of synteny determined above for B. napus were used for comparative mapping with Arabidopsis against the Lle QTL for blackleg (Leptosphaeria maculans) (154). Lle had previously been mapped so as to narrow the QTL down to a locus, however the identity of the gene is unknown; using EST markers that co-localize with this locus has helped narrow the candidate genes down to seven possible RGAs. The EST markers that the Lle locus co-segregates with are a WD-40 repeat family protein and a nucleoporin family protein (154). Even though the WD-4O repeat family protein wasn’t considered a candidate gene, coincidentally one of the QTLs for clubroot resistance in B. rapa also encompasses a WD-40 repeat family protein sequence (236). Even though there are no reports of WD-4O proteins being involved in resistance, this would be a potential candidate gene based on position. Candidate gene analysis using crop-specific EST data Attempts at using comparative genomic analysis for the identification of resistance genes against Fusarium head blight (FHB) proved difficult for wheat, due in part to a lack of collinearity between the QTL in wheat and the corresponding area in rice. The breakdown of synteny at that region was due to rearrangements of the chromosomes in wheat and rice (140). As such, a previous attempt at using synteny between wheat and rice to identify candidate R genes or DR genes failed, in part because there were no annotated R or DR genes present in that chromosomal area of rice. A new approach to find candidate genes for the F HB resistance QTL utilized the vast amount of 24 EST sequences available for wheat that are mapped onto chromosomes. A low stringency tBLASTx search returned three ESTs with homology to R genes from other cereal species. One of the ESTs; homologous to a barley R gene firnctional against rust, was found to be polymorphic between parents used for analysis, and is a good candidate for the R gene responsible for the QTL (227). Synthesis of genomic and genetic data: broad spectrum disease resistance Overall synthesis of existing information on disease resistance QTLs, cDNA libraries and the genome sequence can also yield positional candidates for disease resistance. By combining QTL data for resistance to firngal, bacterial and viral pathogens against a genomic map, areas of multiple QTL congregation and RGA clustering are hypothesized to be involved in broad spectrum resistance (255). RGAs were identified by nucleotide homology to other known NBS-LRR sequences and amino acid similarity to a conserved domain. In addition to mapping QTLs, RGAs and known R genes, a digital Northern approach was used to analyze the occurrence of ESTs fiom inoculated leaf tissue by using full length cDNAs to perform BLASTn searches against the EST libraries. Full length cDNA with significantly different numbers of EST matches were placed on the genomic map in order to determine if they also co-localize with broad spectrum disease resistance regions. Four chromosomal areas were identified by the concurrent mapping of QTLs along with other genetic and genomic data. Positional candidate genes from each of these areas were determined based on homology to annotated Arabidopsis proteins with homology to the translated cDNA sequences from 25 that particular chromosomal segment (255). Further narrowing of the areas on the chromosome that QTLs encompass is still required for better introgression of these broad spectrum disease resistance areas into new cultivars, but this was a first important step in the melding of current breeding and genomic data for rice for disease resistance. Concluding Remarks Since the majority of the molecular mechanisms of E. amylovora virulence are still elucidated by mutational analysis, I chose to use an in vivo expression technology (IVET) approach to identify genes upregulated in E. amylovora during colonization and virulence. This work is reported on in chapter one and pertains to the diversity of upregulated genes in E. amylovora that are required for cellular processes during virulence beyond just plant-microbe interactions. In the NET screen, genes from a type II secretion operon and an endo-polygalacturonase were identified. Since type II secretion was previously though not to be important for E. amylovora virulence, I chose to characterize both the type II operon and endo-polygalacturonase and to determine their contribution to virulence in chapter two. Finally, with the increased amount of genomic data available for plant systems, scientists are able to apply knowledge about disease response in model plant systems to crop plant systems. The use of RGAs and defense related genes as molecular markers and candidate genes for QTL analysis is helping to identify new introgressible R genes for agricultural use. However, in a slow growing crop species such as apple, the ability to identify resistance loci from QTLs is not easy and is time consuming. In order to identify discrete genes that are involved in the resistance response to E. amylovora, I chose to use a suppression subtractive 26 hybridization library combined with expression analysis with cDNA macroarrays. Chapter three concerns these genes that are upregulated in the moderately resistance cultivar Red Delicious, and are part of the resistance response. 27 CHAPTER 1 Identification of Erwinia amylovora Genes Induced During Infection of Immature Pear Tissue 28 ABSTRACT I used a modified in vivo expression technology (IVET) system to identify E. amylovora genes that are activated during infection of immature pear tissue, and identified 394 unique pear fruit-induced (pfi) genes on the basis of sequence similarity to known genes. Known virulence genes, including hrp/hrc components of the type HI secretion system, the major effector gene dspE, type H secretion, levansucrase lsc, and regulators of levansucrase and amylovoran biosynthesis were upregulated during pear tissue infection. Known virulence factors previously identified in E. (Pectobacterium) carotovora and Pseudomonas syringae were identified for the first time in E. amylovora and included HecA hemagglutin family adhesion, Peh polygalacturonase, new effector HothoCEA, and membrane-bound lytic murein transglycosylase MltEEA. An insertional mutation within hothoCEA did not result in reduced virulence; however, an mltEgA knockout mutant was reduced in virulence and growth in immature pears. 29 INTRODUCTION Erwinia amylovora is the causative agent of fire blight, a devastating necrotic disease affecting apple, pear and other rosaceous plants. Entry of the bacterium into plants can occur via flower blossoms, actively growing young shoots or through wounds. Upon entry, the fire blight pathogen moves through intercellular spaces towards the xylem, and also the cortical parenchyma (247). Symptoms often appear as water-soaked tissue that rapidly wilts and becomes necrotic, leading to the characteristic “shepherd’s crook”. As a member of the Enterobacteriaceae, E. amylovora is related to many important human and animal pathogens such as Escherichia coli, Yersinia pestis, Y. enterocolitica, Salmonella enterica and Shigellaflexneri. Like many other Gram-negative plant pathogenic bacteria, E. amylovora produces a type IH Hrp secretion (TT88) apparatus that delivers effector proteins into host plants (125). The TTSS in E. amylovora controls the ability of E. amylovora to cause disease in susceptible host plants and to elicit the hypersensitive response (HR) in resistant and non- host plants. Most hrp genes have been found to encode proteins involved in gene regulation or in assembly of the TTSS apparatus (4, 99, 125). The TTSS of E. amylovora secretes several virulence proteins, including HrpA, HrpN, HrpW and disease-specific protein DspA/E (hereafter referred to as DspE) (25, 26, 82, 125, 126, 251, 252). The HrpA protein is the major structural protein of a pilus called the Hrp pilus, which is the extracellular part of the TTSS (117). DspE, HrpN, and HrpW proteins are effector proteins of the TTSS and are believed to be injected directly into host cells (25, 26). 3O Additional E. amylovora virulence factors that contribute to pathogenesis and plant colonization include the exopolysaccharides amylovoran and levan, iron- scavenging siderophore desferrioxamine, metalloprotease PrtA, multidrug efflux pump AcrAB, and carbohydrate metabolism genes specifically involved in the utilization of sorbitol, sucrose and galactose (3, 29, 38, 161, 264). Transcriptional regulators of the amylovoran and levan biosynthetic operons have also been identified (22, 51, 264) and are required for the expression of the biosynthetic machinery for the exopolysaccharides (19, 68, 123, 264). E. amylovora pathogenesis is also subject to global regulation by the small regulatory RNA rsmB which functions by titrating and countering the activity of the repressor protein RsmA; this system is reported to positively regulate exopolysaccharide production, motility and pathogenicity (147). In addition, E. amylovora strains contain a ubiquitous nonconjugative plasmid of 28-30 kb designated pEA29; laboratory-derived plasmid-cured strains exhibit a reduction in virulence (155). pEA29 encodes several potential virulence genes including a thiamine-biosynthetic operon that is proposed to influence amylovoran production (155). Genetic analysis of virulence genes in E. amylovora have been performed mostly through the production and screening of mutants. Additionally, most of the genes discovered so far have been identified from mutant screening under controlled conditions. However, it is not feasible to mimic all of the nutrient and defense conditions in vitro to characterize all the genes from E. amylovora required for infection and colonization of plants. There is a need then, for a hi gh-throughput method of screening for genes that are involved in virulence and grth in planta of E. amylovora. 31 In the last decade, many gene expression technologies including in vivo expression technology (IVET) have been developed to identify gene expression profiles of organisms during interactions with various host environments (8, 101, 150). IVET screening theoretically scans the entire genome, and through the use of appropriate environmental conditions and different strategies, can yield large numbers of potentially important genes (201). IV ET screens have identified genes upregulated upon infection with enteric human and animal pathogens such as Salmonella enterica, Shigella flexneri, and Y. enterocolitica (15, 150, 171, 261). IVET systems have also been used to identify genes expressed during plant infection by Xanthomonas campestris, E. chrysanthemi, Pseudomonas syringae and Ralstonia solanacearum (24, 36, 96, 172, 187, 254), phyllosphere colonization by P. syringae (151), and saprophytic colonization by P. fluorescens (200, 225). Like many plant pathogenic bacteria, E. amylovora can infect different host tissues at different stages of disease development. E. amylovora not only infects blossoms, leaves, and succulent shoots, but also immature finits of susceptible hosts. The bacterium also grows epiphytically on stigmas and endophytically inside plant tissue. The maintenance of large numbers of apple trees for study of E. amylovora pathogenesis is quite difficult due to the extensive greenhouse and growth chamber space required. As an alternative, many researchers have utilized immature pear fruits to study E. amylovora infection (19, 26, 85). Immature pear infection is initiated through a wound inoculation; wound colonization is a frequently-utilized mechanism of E. amylovora infection in nature (247). Immature pear assays, either using intact pear fruits or pear slices, have 32 been used successfully to analyze virulence effects of several E. amylovora genes (26, 85, 125) Although key virulence factors contributing to fire blight have been identified, little knowledge is available of the global host-regulated genes of E. amylovora during infection. To gain a better understanding of the molecular mechanism governing E. amylovora—host plant interactions, we undertook a comprehensive genome-wide examination of gene expression patterns during host infection to uncover pathogenesis strategies of the organism and to lay the groundwork for future studies examining the expression and function of critical virulence genes during infection of different host tissues and survival within the host. Several known virulence and pathogenesis factors were identified using this modified IVET screen, along with new potential virulence genes that were previously only described in other bacterial pathosystems. We also confirmed that infection of immature pear tissue by E. amylovora required the major pathogenicity factors of the bacterium. 33 MATERIALS AND METHODS Bacterial strains and plasmids. The bacterial strains and plasmids utilized in this study are listed in Table 1. Erwinia amylovora wild type (WT) and mutant strains and E. coli strains were grown in Luria Bertani (LB) medium at 28°C and 37°C, respectively. Antibiotics were added to the culture medium at the following concentrations: rifampicin, 100 ug/ml; kanamycin, 30 rig/ml; gentamicin, long/ml and arnpicillin, 100 rig/ml. Oligonucleotide primers used for polymerase chain reaction (PCR) and sequencing in this study are also listed in Table 1. DNA manipulation and sequence analysis. Plasmid DNA purification, PCR amplification of genes, isolation of fiagments from agarose gels, restriction enzyme digestion, T4 DNA ligation, and Southern hybridization were performed using standard molecular procedures (217). Chromosomal DNA was isolated using a genomic DNA purification kit (Qiagen, Valencia, CA). TAIL-PCR was performed using the degenerate primers ADl, AD2 and AD3 as described previously (142) and PCR products from secondary and tertiary nested PCR were used for sequencing. DNA sequencing was performed at the Genomic Technology Support Facility at Michigan State University. The Oligonucleotide primer Aj1585 corresponding to the 5’ end of uidA gene was used for sequencing fiagment inserts cloned into the pGCMO plasmid. Sequence management and contig assembly were conducted using DNAStar sofiware (DNAStar Inc., Madison, WI, USA). Database searches were conducted using the BLAST programs at NCBI (www.ncbi.nlm.nih.gov/BLAST). Percent similarity was also calculated using the 34 BLAST program (7). Amino acid alignments were done with ClustalW, v. 1.83 (European Bioinformatics Institute, Cambridge, UK). TABLE 1. Bacterial strains, plasmids and primers used in this study Strains, plasmids, primers Relevant characters or sequences (5’—3 ’)" Reference or source Erwinia amylovora EallO EallO' Eal 189 CFBPl43O M52 (Ea dspA) Ea l 10 hrpA ZY C 1 -3 (Ba hothoCEA) ZYE3-l l E. coli DHlOB 817-1 817-1 2. pir Plasmids pBluescript II SK(+) pGem32f+ pCAMl40 pCAM140-MC8 le9l8GT pBSL15 pGCMO Wild type, isolated from apple Eal 10, cured of pEA29 Wild type, isolated from apple Wild type, isolated fiom Crataegus CFBP1430, dspA::uidA-l801acZAM15 AlacX74 recAl endAl araAl39 A(ara, leu)7697 galU galK A. - rpsL (StrR) nqu recA pro hst RP4-2-Tc::Mu-Km::Tn7 it -pir lysogen of 817-1 ApR, cloning vector ApR, cloning vector SmR, SpR, ApR, R6k origin, an588gusA4O ApR, R6K origin, pCAM140 derivative without mini-TnS, contains the multiple cloning site of pBluescript II SK (+) xylE-GmR fusion cassettes-containing plasmid flanked by inverted repeats of the pUC19 MCS Km cassette flanked by inverted repeats of the pUC18 MCS GmR cassette with downstream transcriptional terminator and gusA with upstream translational stop codons in pGem32f+ (155) (155) (38) (82) (82) (l 17) This study This study Invitrogen (253) (253) Stratagene Promega (253) (3 8) (224) (8) This study 35 TABLE 1 (con’t) Strains, plasmids, Relevant characters or sequences (5 ’—3’)“ Reference or source primers pZYFZ 570 bp dspE promoter in opposite orientation This study relative to uidA in pGCMO pZYF 8 570 bp dspE promoter in correct orientation This study relative to uidA in pGCMO Primersb Aj 1388 CCCAAGCTTGGTGCGCCAGGAGAGTTGTTG (Hind HI) Aj 1389 AAAACTGCAGTGATTGATTGACGGACCAGTATTATTATC (PstI) Aj 1390 CCGGAATTCCGAATTGACATAAGCCTGTTCGG (EcoRI) Aj 1391 CGGGGTACCTGGACGCGGCCGATCACCTGGCCG'ITG (Kpnl) Aj 15 85 GATAATAATACTGGTCCGTCAATC Aj 1565 CGGTI'I‘ACAAGCATAAAGCTGGGCAACGGCC DspEl TCCCCCGGGCAGTGAGGGGGGGCAGACTTITI’I‘TTAACC (Smal) DspE2 TCCCCCGGGTATCTTCGCCGCTGCCACCTTTCACCATTG (SmaI) PtoCl TCCCCGCGGGCGGGCTGTTGGTC’IT GCT CT (SacII) PtoC2 TGCT CT AGACTCT GGCAAAATTCAACTGA (XbaI) PtoC3 CCGGAA'ITCCATGGCAGGGACCCGCAGTTTG (EcoRI) PtoC4 CCGCTCGAGGGCTGATGGCGGGTTAGTCTGTCG (Xhol) MltEl TCCCCGCGGTGAATAGTGCGTGGCGTGATGTGC (SacII) MltE2 TGCT CT AGATTAATCATTGCAATCGCCT CGTC (XbaI) MltE3 CCGGAATTCTI‘ACCAGCACGTGCAGACAAAACA (EcoRI) MltE4 CCGCTCGAGCCGGATGGATCTGGTGAGGGGCGC (Xhol) ADl NTCGASTWTSGWGTT AD2 NGTCGASWGANAWGAA AD3 WGTGNAGWAN CAN AGA R R R R . . . . . . . . 0 Km , Ap , Gm , SpR , Sm = kanamycrn, amprcrlhn, gentarrucrn, spectmomycrn, and streptomycin resistance, respectively. Underlined nucleotides are restriction sites added and the restriction enzymes are indicated at the end of primers. Mixed nucleotides: 8 = C+G; W = A+T; N = A+T+C+G. Immature pear infection assays. Immature pears are routinely used to examine the pathogenicity of naturally-occurring isolates or bacterial mutants of E. amylovora (26). In order to confirm that infection of immature pear required major pathogenicity factors as previously reported (26), we inoculated wounded immature pear fruits with E. amylovora M52 (CFBP143O dspE) and Eal 10 hrpA mutants and monitored for symptom development and in planta bacterial growth. Bacterial suspensions of all strains were 36 grown overnight in LB broth , harvested by centrifugation, and resuspended in 0.5X sterile phosphate buffered-saline (PBS) with the cells adjusted to approximately 1 x 104 colony-forming units (CFU)/ul (OD600 = 0.1 and then diluted 100 times) in PBS. Immature pears (Pyms communis L. cv. ‘Bartlett’) were surface sterilized with 10% bleach and pricked with a sterile needle as described previously (155). Wounded pears were inoculated with 2 ul of cell suspensions and incubated in a humidified chamber at 28°C. Symptoms were recorded at 2, 4, 6, and 8 days post inoculation. For bacterial population studies, the pear tissue surrounding the inoculation site was excised by using a #4 cork borer as described previously (26) and homogenized in 0.5 ml of 0.5 X PBS. Bacterial growth within the pear tissue was monitored by dilution-plating of the ground material on LB medium amended with the appropriate antibiotics. For each strain tested, fruits were assayed in triplicate, and each experiment was repeated. Construction of the genomic library of transcriptional fusions to uidA. We used E. amylovora Ea110' (cured of the ubiquitous plasmid pEA29) as the source of chromosomal DNA for the IV ET experiments. We excluded pEA29 genes from this study because an analysis of the expression of pEA29-encoded genes during infection will be presented in a separate report (McGhee and Sundin; unpublished). To create a library of transcriptional fusions, chromosomal DNA from E. amylovora Eal 10' was partially digested with HaeIII, and fiagments between 800 bp and 2 kb in length were separated by electrophoresis and gel-purified. The purified fragments were ligated into pGCMO prepared by SmaI digestion and transformed into WT E. amylovora Eal 10 (containing pEA29) by electroporation. The use of WT strain Eal 10 was necessary 37 because the ubiquitous pEA29 plasmid contributes to E. amylovora virulence (155). The 6.2-kb pGCMO reporter vector was constructed by cloning the aacCI gene (confening resistance to gentamicin) into the EcoRI and KpnI sites and the promoter-less uidA (B-glucuronidase) reporter gene into PstI and HindIII sites of pGem3zf through multiple cloning steps (Figure 1A). The aacCI gene was amplified from plasmid pX1918GT by PCR using the primer pair Aj1390 and Aj139l, whereas the promoter-less uidA gene was amplified from plasmid pCAMl4O using the primer pair Aj1388 and Aj1389. A transcriptional terminator sequence, also from pX1918G, was located immediately downstream of the aacCI gene, and we added translational termination codons in all three reading flames upstream of the uidA gene. The pGCMO vector was first digested with SmaI and the ends were dephosphorylated with calf intestinal alkaline phosphatase (CIAP) and checked for self-ligation before ligation with E. amylovora chromosomal DNA fragments. After ligation, DNA was introduced into Eal 10 by electroporation, and transformants growing on LB medium amended with gentamicin and arnpicillin were randomly collected and plasmids were recovered. The randomness of the inserts in the IV ET collection was confirmed by checking insert size from 30 random colonies through restriction digestion and PCR (data not shown). I As a control, we cloned a 570-bp fiagment containing the dspE promoter into pGCMO. The fragment was amplified by PCR from strain Ea110 using the primer pair DspEl and DspE2. The resulting 570-bp product was cleaved with SmaI and ligated into pGCMO in both orientations. The resulting plasmids were designated as pZYF2 (dspE'wzszCMO, dspE promoter in opposite orientation to uidA) and pZYF8 (dspEf°’::pGCMO, dspE promoter in correct orientation to uidA), respectively, and each 38 plasmid was introduced into strain Ea110 by electroporation. Screening of the E. amylovora IVET library using a GUS-based microtiter plate assay. An in vivo microtiter plate assay was developed for screening of the E. amylovora IVET library (Figure 13). Briefly, approximately 19,200 transformants in strain EallO were randomly collected and initially screened for GUS activity on LB plates containing 5-bromo-4-chloro-3-indolyl-B-D-glucuronide (Xgluc). After incubation at 25 °C for 48 h, bacteria were transferred individually using a 48-pin colony transfer apparatus and inoculated onto immature pear discs (3mm) in 96-well microtiter plates. Intact pears were surface sterilized using 10% bleach for 10 min and rinsed three times with sterile water. Discs were cut from pears using a #2 cork borer and immediately immersed into microtiter plate wells containing 25 ul 0.5X PBS buffer to avoid oxidation. The microtiter plates were then covered with AirPore tape (Qiagen, Valencia, CA) after inoculation and incubated in a humidity chamber at 25 °C for 48h. After incubation, a qualitative GUS assay was performed as described below. Transformants showing GUS activity on pear tissue but not on LB plates were selected and re-screened on LB plates containing Xgluc and re-inoculated onto pear discs in 96-well microtiter plates. Confirmed differentially-expressing transformants were again selected and stored at -70 °C in glycerol stocks for firrther analysis. Plasmids were isolated from the consistent differentially expressed transformants and were end sequenced to identify the genes or promoter regions. Transformants showing GUS activities on both LB plates and on pear tissues were assumed to contain constitutively expressed firsions and were not analyzed further in this study. 39 r] y :1 Construction of hothaCEA and mltEEA mutants. For the construction of hothoCEA and mltEEA mutants, the sequences of the putative ORFs defined by the corresponding clones were determined and used to design primers to amplify fragments of the genes and its upstream and downstream sequences. Primer pairs PtoCl and PtoC2, PtoC3 and PtoC4 were used to amplify 590 bp and 670 bp fragments from E. amylovora strain Eal 189 corresponding to the upstream and downstream of hothoCEA gene, respectively. Primer pairs M1031 and MltE2, MltE3 and MltE4 were used to amplify 700 bp and 560 bp fragments from E. amylovora strain Eal 189 corresponding to the upstream and downstream of mltEEA gene, respectively. The two fragments for each ORF were cloned into pBluescript-H SK(+ ) through multiple cloning steps with corresponding restriction enzyme digestion (SacII and XbaI; EcoRI and Xhol, respectively). The whole fragment was excised using SacII and Xhol, gel purified and cloned into the suicide vector pCAM- MCS (38) digested with the same enzymes. The resulting plasmids were digested with SmaI and ligated with a 1.2 kb fragment of the aph gene (conferring kanamycin resistance) released from plasmid pBSL15. The final plasmids were designated as pZYC8 and pZYE8, respectively and introduced into E. amylovora strain Eal 189 by electroporation. Transconjugants resistant to kanamycin (Km) were selected. To further exclude mutants resulting from single crossover events, transformants were selected on LB plates supplemented with Km and onto LB with arnpicillin (Ap). Km-resistant and Ap-sensitive colonies were selected and their genotypes were confirmed by hybridization or PCR analysis. 40 GUS Assays. The B-glucuronidase (GUS) reporter gene (uidA) on pGCMO was used to monitor promoter activity of IVET clones both in vitro and in viva. Qualitative GUS activity of IVET clones was monitored visually by the development of a blue color within 48 h of cells on LB medium containing 1 mM Xgluc. Qualitative GUS activity of IVET clones grown on pear slices in microtiter plates after 48 h at 25 °C was also monitored visually by adding 10 ul of 20mM Xgluc into the wells followed by incubation for 30 min at 37 °C. The development of a blue color indicated GUS activity. To monitor the expression of IVET clones in pear tissue, quantitative GUS activity of bacteria in either culture or in pear tissue was determined as described previously (26, 116) using 4-methylumbelliferyl-B-D-glucuronide (MUG) as a substrate and 0.2 M Na2CO3 as stop buffer. Briefly, E. amylovora strains containing the IVET clones were grown on LB medium, resuspended in 0.5 X PBS, and inoculated in immature pear fi'uits as described above. At 0, 24, and 48 h post inoculation, the pear tissue surrounding the inoculation site was excised using a #4 cork borer and homogenized in 0.5 ml 0.5 x PBS. Forty microliters of homogenate was mixed with 160 pl of GUS extraction buffer. Reactions were stopped by NazCO3 addition and fluorescence was measured using a SAFIRE fluorometer (TECAN Boston, Medford, MA). Bacterial cell numbers in the sample were estimated by dilution plating, and GUS activity (umol of 4-methylumbelliferone (MU) produced per min) was normalized per 109 CPU (26). Three replicate fi'uits for each strain were tested and the experiment was repeated. GenBank Accession numbers. Nucleotide sequence data reported for the hothngA 41 and mltEEA genes were deposited in the GenBank database under the accession no. AY887538 and AY887539. 42 RESULTS Development of an I VET system and identification of E. amylovora upregulated genes during immature pear infection Our immature pear infection results clearly demonstrated that infection of immature pear by E. amylovora required major pathogenicity factors (Figure 2). At 48 h afier inoculation, E. amylovora strains CFBP 1430 and Eal 10 produced water-soaking symptoms in pears with visible bacterial ooze (data not shown). Two different strains were used in this initial experiment because of the availability of mutants of these strains. This work also confirmed that growth and symptom development of WT strains CFBP1430 and EallO were similar during immature pear inoculation (Figure 2). Four days after inoculation, inoculated immature pears showed necrotic lesions and bacterial ooze formation (Figure 2A). Afier eight days, the entire pears showed necrosis, turning black with copious ooze production at the inoculation site (Figure 2A). In contrast, disease symptoms were not observed on immature pears inoculated with either the E. amylovora hrpA or dspE mutants (Figure 2A). Disease symptoms caused by WT strains on immature pear were correlated with high levels of bacterial growth in pear tissue during the four days post inoculation (Figure 2B); however, the CF BP1430 dspE mutant M52 grew only slightly in pears, representing an approximately 105-fold reduction relative to the WT CFBP1430 strain. Populations of the hrp/t mutant declined quickly after inoculation, indicating that the hrpA mutant was not able to survive in immature pear (Figure 2B). Verification that infection of immature pear required major pathogenicity factors 43 of E. amylovora facilitated the development of a simple, high-throughput IVET system to identify E. amylovora genes induced during colonization and infection. We used cores of immature pear tissue in a microtiter plate format, a system that was conducive to handling large numbers of samples. To develop an immature pear fruit assay, we used the B-glucuronidase gene uidA as a reporter (116) in the vector pGCMO which was constructed as described in Materials and Methods (Figure 1A). The vector was verified with a control construct containing the dspE promoter in both orientations (Table l). The dspE promoter was previously reported to be strongly induced during immature pear infection (26). GUS activity was not observed after two days gowth in LB medium for either Ea110(pZYF8) (dspE promoter in correct orientation to uidA) or Ea110(pZYF2) (dspE promoter in opposite orientation to uidA). However, GUS activity was observed for strain Eal 10(pZYF 8) in qualitative assays two days following inoculation onto immature pear discs but not following inoculation of strain Eal'10(pZYF2) (data not shown). GUS activity was not observed for the WT Ea110 strain containing the empty pGCMO vector on either LB medium or in pear discs. To identify E. amylovora genes expressed during colonization and infection of pear discs, we constructed a library of 0.8 to 2-kb fragnents of genomic DNA of Eal 10' (cured of PEA29) in pGCMO and introduced the library into WT Ea110 by electroporation. In order to screen for differentially-expressed promoter fusions, we developed an in-planta pear disc microtiter plate assay (Figure 13). Strain Eal 10' containing library clones were first gown on LB/Xgluc medium for two days, visually pGCMO I. Ligate Ea110- DNA laacC uidA I into pGCMO Sma l v 2. Transform library 0 0 directly into Ea 110 O O V 3. Screen transformants on LB media for GUS activity 4. lnoculate all clones on immature pear disks 5. After 2—day incubation, screen transformants again for GUS activity FIGURE 1. Overview of the IV ET screen for E. amylovora genes induced during infection of immature pear discs. (A) Schematic map of the IVET vector pGCMO. The 6.2-kb pGCMO vector was constructed by cloning the aacCI gene (conferring resistance to gentamicin) into the EcoRI and KpnI sites and the promoter-less uidA (B- glucuronidase) reporter gene into PstI and HindIII sites of pGem32f through multiple cloning steps. The (D symbol represents translational terminator codons in all three reading frames upstream of the uidA gene; and the symbol 0 represents a transcriptional terminator sequence immediately downstream of the aacCI gene. The Smal site was used for ligation of random chromosomal inserts. (B) A library (19,200 clones) of Smal chromosomal DNA fragments (0.8 to 2 kb) from E. amylovora was constructed in pGCMO, transformed into E. amylovora Eal 10 and screened individually for GUS activity on LB medium amended with Xgluc. A 96-well microplate containing slices of pear tissue was inoculated with EalOO containing random IV ET fusion clones and incubated for 48h at 25°C. Clones exhibiting GUS activity on pear discs but not on LB+Xgluc medium were selected and the plasmids were recovered for further analysis. 45 monitored for GUS activity, and then inoculated onto pear discs in nricrotiter plates (Figure 1B). GUS activity was qualitatively detected after two days of incubation at 25°C . Only clones that showed high GUS activity in pear discs but no GUS activity on LB plates were recognized as pear-upregulated clones. Those differentially-expressed clones were screened again on LB/Xgluc plates and pear discs to confirm the results. A total of 19,200 transcriptional fusion clones were screened on both LB/Xgluc medium and pear discs, and 498 clones (2.5%) were repeatedly found to differentially express GUS activity on pear discs in this qualitative assay. Sequence analysis of E. amylovora genes upregulated in immature pear tissue. We determined the sequence of the inserts from the 498 clones and identified the putative genes induced following BLAST searches of the non-redundant GenBank database. Of the 498 inserts sequenced, a total of 55 genes were identified two or more times and 12 clones contained either an intragenic sequence or with the putative gene present in the incorrect orientation. Although it is possible that these 12 clones may contain cryptic promoter sequences as has been shown in a previous study with P. fluorescens (66), we separated the clones from the others in the current study and did not subject them to further analysis. Thus, a total of 394 unique putative pear-inducible genes were identified, and these pear fruit-induced (pfi) genes could be divided into nine putative functional goups, including host-microbe interactions (3.8%), stress response (5.3%), regulatory (11.9%), cell surface (8.9%), transport (13.5%), mobile elements - phage (1.0%), metabolism (20.3%), nutrient acquisition and synthesis (15.5%), and unknown or hypothetical proteins (19.8%). 46 Figure 2. Symptoms and bacterial gowth of Erwinia amylovora WT strains and hrpA and dspE mutants in immature pear. (A) Symptoms caused by Erwinia amylovora Eal 10, CFBP1430 and corresponding hrpA and dspE mutants in immature pear. DPI: days post inoculation. (B) Bacterial gowth of Erwinia amylovora WT EAllO, CFBP1430 and hrpA and dspE mutants during infection of immature pears. The gowth of bacterial strains was monitored at 0, 1, 2, 3, 4 days after inoculation. Data points represent the means of three replicates i SE. Similar results were obtained in a second independent experiment. 47 Log CFU/ g tissue CFBP1430 dspE' 4 DPI 8 DPI Eal 10 hrp/1' Eal 10 hrpA' i CFBP1430 dspE 11 10‘ 9. 8‘ +Ea11_0 +hrpA 7‘ +CFBP1430 —l—- dspE 6< 5. 4| _. 3. 1. 0 1 2 3 Days post inoculation 48 The majority of the putative gene-products identified as inducible during infection of pear tissue shared high amino acid similarity with that of proteins from Yersinia spp., Salmonella spp., E. coli, Shigella spp., and Erwinia spp. (Table 2). Several known virulence factors previously reported in E. amylovora such as the TTSS genes hrpGF, hrpL, hrpX and genes encoding known or new effector proteins DspE and HothoCEA were upregulated during pear infection (Table 2). Other known E. amylovora virulence genes identified as upregulated in this study were levansucrase (lsc), regulator of levansucrase (rlsA), amylovoran regulator (rcsA), and zinc-binding metalloprotease (prtA) (Table 2). In addition, genes encoding polygalacturonase (peh), a hemaggulatinin- farnily adhesion (hecA), and membrane-bound lytic murein transglycosylase (mltE) were identified for the first time in E. amylovora Eal 10 (Table 2). Peh and HecA are important virulence factors in E. chrysanthemi (205), and MltE plays a role in the virulence of P. syringae (24). A total of 54 upregulated genes identified were homologs of genes identified in IVET studies performed with other bacterial plant or animal pathogens (Table 2). Type H secretion system (T288) genes similar to Yersinia enterocolitica ytsIIJ (112), a known virulence factor, and one of five major protein secretion systems in many pathogenic bacteria, were identified for the first time in E. amylovora. YtsllJ are known type 4 pilin-like proteins (pseudopilins). Interestingly, the type H secretion system is dependent on the general secretory pathway (GSP), i. e. the Sec secretion pathway (112). In our study, the major preprotein translocase SecA (ATPase), molecular chaperone SecB and membrane proteins SecDF of the GSP were also induced during infection of pears along with the T288 (Table 2). Furthermore, the peh gene, encoding an enzyme that is 49 known to be secreted by the T288, was also up-regulated (Table 2), indicating that a functional T288 is present in E. amylovora and to a geater extent could contribute to the virulence of this pathogen. Transport genes (pfi 16 to 51) including general, ion, sugar, amino acid, peptide, and nucleotide transport proteins were induced in pear tissues (Table 2). Some of the transporters may belong to the type I secretion system that is known to be involved in secreting toxins, proteases and lipases and are potential virulence factors in E. amylovora. Cell surface proteins including inner, periplasmic, and outer membrane proteins, lipoproteins, flagella and polysaccharide proteins were also induced during pear tissue infection (pfi 52 to 72; Table 2). These membrane proteins may be involved in protein secretion and membrane maintenance. The sensor component (envZ, pfi 94) of a two- component regulatory system and cogiate outer membrane protein genes that this system regulates ompA (pfi 65) and ompC (pfi 64) were also differentially expressed in pears compared to LB medium. Under unfavorable conditions such as nutritional stress or exposure to a host defense response, bacterial pathogens respond by over-expressing stress response genes. Several stress response genes (pfi 75 to 90) were identified in our screen (Table 2). These genes included DNA repair or protection (mutS, recA, sulA), carbon starvation, heat or phage shock and antioxidant genes (such as grpE and ahpC). These results suggest that pear tissue at least initially is not a favorable habitat for E. amylovora gowth and /or that DNA damage and the neutralization of plant-derived reactive oxygen species are involved in virulence and in planta gowth. The sensor component of a two component regulatory system, grrS (pfi 93), was 50 identified as upregulated in this study. GrrS is a homolog of GacS which, along with GacA, globally regulate a network of virulence functions in E. carotovora, including the production of quorum-sensing sigraling molecules (72). Besides amylovoran and levansucrase regulators (rcsA and rlsA), and genes encoding sensor component of a two- component regulatory system (grrS and envZ), our screen identified fliZ, a positive regulator of the flagellar biosynthetic operon in enterobacteria, as upregulated. Other regulatory genes (pfi 96 to 126) and genes involved in sequence mobilility (pfi 73 to 74) and metabolism and nutrient scavenging (pfi 127 to 157) were identified in our screen and listed in Table 2. We recovered several metabolic genes that are potential precursors for the siderophore desferrioxamine biosynthesis in this study (pfi 153 to 157). It is probable that under unfavorable conditions, the bacterium itself adjusts and overcomes nutrient and iron deficiencies. The large number of unknown or hypothetical proteins identified in this IV ET screen (78 genes, 19.8%) indicates the firture possibilities of characterizing novel virulence traits in E. amylovora and assigring functions to these proteins. A complete genome sequence of E. amylovora is expected soon. When an annotated genome sequence is released, we will make a listing available upon request of the gene numbers of the unknown or hypothetical proteins identified in this study. 51 Ga .88 8381.2 8:88.. EC. 538538 888.5% m8 8 mNmooclmZ A5085 amp—Sm :0: :3 Etta.“ urban—9030 Emma £5.98» 8.x 2 2»: 5m mug—Gee omfifiamofiS—oaoxm 8882833» 882.8% x&& w— Nvmoovlmz 888qu 889? ton—meg 839.39 93845 Emma 55me $3 5 champ—ab com :83 aom 5341.2 2.8983 cocoa magi: 8:828 cm< 8:88 8322 So 2 team—3.: «2822:. .E2. 2:5: came. .8. £2239 mnwmvmad 83852—808 wamvfintoEN 888188 SEEM .. SLR n— ?NV .3; Namovflmz emu—xmoobwmeflu 5238 our cgcnocfinEo—Z Marni 8338» “.38 3 an 5.x ”2325. <8: eseo oxzéeaawaeoéeoboeom 88888 858.8 woo: 2 Ge .oom 5292 ea_oeoa\=§.=_wm§o€oo€o oases :8 58.888 88‘ 2 3.893 8.8.63. .85: .3 38.88 2— 3 5755— 2:895 $8 .me ~mv§<Exa:u:=8> .888 838.8888m .. V38 mv on— Swizz 5 88:58:33 3288-883»; magma 8.8.3.88. One .8 E852 885388-38 .88 .8 .885 8 33 .on 588.82 5228 888388 8:888. RE 28 822 888 8:88.. 828.: an $32.18 Z 888qu 85836593 ~8888§8-8w3w $5 mm £2 :32 2 529a 8582: 83% .888: 8:88: .88 88.5888 3588 R 529?. ton—mug Um< Boa 9:88 520685538885 888.82 888.888 853985888582 3858-85 888 588: 838 on 5895 wfiufin o2 Sela—Z -omontAH egg—atom cor—88mg 0% 83.;— 88888~8M8 88888883883 mm? mm and .E Swaoslmm 38.8.88 8228.8 8.36 8:8 888:5.— 858: 888288 .8228; cm Nwmam<<< 838.53 2585 in 828888” 238% mm 3235 .888 858528.. 8878.582 952:: 8.88288 888882 8% mm .98 88%“ § 5.82 m 5228 8.582 952:2 .88 588888 8:8 R 3H— 5qu £888 ton—28b oEmSoflawAEHoE 58.285383 8288888” 8 EMS cm mvvchIm Z 83538 8.28on 89858453 Eatafiimb 828.88%. 388% am $38an 888:8 828-8 53.88-85 8.888 85.8; 9.8 mm *6.vach 828355 82:0 888E888 8.88888888l 3.8 um 82?? 9888.88,: 885988 22: 32 -1— 358in 8 889:8 85858-8? 88 28.8888 888888 8 Ga .88 828182 8883 323283.28? 83838388 8888838 :83 2 mnmnmwlmz £205 mnxm teams“: BEA—892m 88.2888 8-8888~8m M38 vm E89? 82 SoocImN tong Baconuafiéoficobsmbugfi 995-0% 888888885885 888888D 2.3886 mm 83o3<< tame“: .8: Eu 8882 285888 8880 888888888 8.838% N88 mm $8832 5888 85838-92 82.88 888.28 8.58 88 88.5888 .58 s 888.882 .0: .oo< .o: b95385 455800 Emmomlfimflm 88.8830 88.88% Km 383 N 39:. 53 2&meth 280V 8088 285808 $80 2885 8.228 33.88%. _._ V88 3 :33 69:0 53% 0 522a 8882. 530 :8 38:38 05% E mmmocVImZ 356.68 .33 8208 28.5808 385 02:85 53.235388 SBEVEBM Fix mo Ge .5 $8232 528 0833ch 3:58 53.23538 8%st .3 85m 8 38 and $83th 8088 28.5808 8:8 98.85 8.288 Sfincéum _. “NEE S gob 88 Sou. 5m? honmwofiu $80828 8088 235808 888 8088 Eek 8858398 583% 52 00 NM 333m Z 8208 883808 888 829$ .888 BE EEEEESO 58838” 38 am 3 2:32 5205 BE?» 263.: .688 3:28 3:6,: mm A $3 .fim 828...? 53.8 08538 025.5 583538 £828 E? S 208 .8m 3: Sela—Z 8808 8898o8 9885 8288 38:285. #18 em 8% .8m ommwmvflmz 8088 235808 o_£o=v=_-own88- 58,92 ms: 538: :82 8.9:”? SE :8: 52.28: as: p 8.8%?— 39.5 88232 83%: NEE 828$ :8 58:88: Q2: 2. mfimmvdz 888:: 82: 2885 8.228 $8853: .825: Q can...— - 8.3529 2392 hmi‘wNImZ E205 gown—3.5% 033:5 :8 339.2393 amok Nb Emvfidz SE 532: :25” 888 8:08.: :8 58.58:: 8% K 585-? 0:". .885 8:8 888 8:82,: t2§§8m§ Di 2 Eammlmz 2»; .888 $855 2.8“: .85»: 58:3: 2% 8 383-? <8: see: wagsouéfim :3. :8 5:38:98: 52: we mmevov1m 2 H858:— 5 833.5 3.8% 852m: .. 53 no wovgdz => 88.2: .3 an 532.. 88:88. 830 :8 £8,238: _. 5:: co :uotoao: .o: .oo< I .o: 39530.5 x5980 gummamom Ema—m =55»ch :m2mb : & 9.8“: N 5mg 55 mace“; Z zd £8qu @22st cowobmz 2:an £32m» MEM as. 8531? 58ch >83st actuate“ BEE Ea BEE—am :8 uSutmfinm .E\ wo— egmoeumz $8562 “285:5 233$ seen rename; ”can :3me .289: 8. $38.? 508a rename 2:55 nag «ES; em 8% 2: $281.2 5065 E332 383822“ bfiééfi ”.53 552.: #38”: 8. mnmmahlmz Baa—:3“ .aaotatomaab bun—um Much @328. SR 3??? Qacsohammm E nabsmmfi X: 23 .fim c8852 «m3 8.2:? 389388“ 2653 :8 sfitmfifi $33 8. 882 833? 38358.3 25-5: 3&52 £238: «35 N2 $821.2 6:5 .22st 383585 Efiéufimsm 9% z: Smammlmz 52.5 9.332 25-05. ”ESE :8 35.23% .22 2: once—51ml 83052 :0 «3:8— omfia 8m 2.28% 8385 tmauukuammim v3: 3 :3st 05:250.: 05 mo EmmoEoSEo—z 22.31% £88m >53:on :03w8 82v 83052 82 388:3»: naheufsoal ENE ma SmmelmZ Samoan.— =ohoqc omcntaxouE zoom 8.28:» 532255; «8% 3 ohmmeImZ .8329: 9:30 5.58:5: Bacon—w 822$ §m§§§w «3% ca wovmovflmz £8qu >85:on 38am «Ema 335:2? 033:5 3me 8.:me Na 3 Co; 5m 3331A? 5205 833 bra—cfimoREmBmaafl 05235 Mama £5.28. *Nhs a a3; 3<< 08:3 888 8.52:? 3:38“ Eco—w 258$ cuiuasbm 3.29:3; .. ”EM mo < 3889 w~wm¢< 553.8 $35583 8338 Eon 2:200 9833.5» 52.32% _. 5.2 Na cave—<15 onEo: mm». 38ch 33 mm— anocnmz 38:3 <5 .8208?&< £3 35$» . is: a. 30 5mm NNNNE <35 88:2. a$=§¢¢< :8 sfiEfifi 35 NN_ Navmvwlm Z £88m 8E8“: Eugene—Ymmcoh gums?» nuzosomohi Sums NS ONSNSOINN 2:5. 8223: £822 2&5 seomésksofigé @2395 _N_ mwbmoocoN—N omniv— anoa cacaofitoctom §=a~m5> .SBBSQNV waVS>V om— Scmovdz mfitomafix 053083 a 25 958.5 $2.3 SSEQ 330%» a: NNwoNodz 82: 35.5 fig 352“: . an? N: homemmlmz 2.53 833w“: emu—cameraman.“ Egan—oi ameSMEma 9.395638% R NNVm >2 96M 55 8893 ©334<< EOE—smog 28888-05 5 3:2st 09.38% 5:23:3So 332°st MQE 0: £33th commmoaxu 28m 82 .3 “Sn—swag 03325 9.28:» SENSE?” .. 335;. m: ONNVNéufi seen 55.. m5— 835 §§§§m 235m 3. N. .332 85% we 3.. saamom .55»: £35 .22: N: 383152 ESSA .8568 Ema—9:0: 258nm Magma 852m» MAE N: fiwmmVIAZ 3:268 .agmatomuab :82? 95 3.225 £3.85?” Ema :— vaowmlmz 528a boa—swag 839. 235302 :8 Sautmfinm $3.: a: $5953» .0: .oo< I .o: bmsoSoi x5980 utmumom 5.35 22.2836 amzmw NE 9.8"; N 339 57 m _ «womlmz N 8535983 2:93 :8 uSutmfiam $3.6m mm _ movmcv1mz 82.23330 Sgaq—oaooaamonm Eatagfig czmzcgum SEQ cm— w 53‘an omgowofimsow 8398:9668030 Manna Ssrmx E mm— Nmaacclmz ommaowoéxzov oafiamonqééfinoga «2me 352m» 3% v2 oSEImz omaéamofifigawosam :8 asutmfifi Ex 3 chfiumz 9F: usufia 033on fig 33%: 329: N2 mg hwowlmz ommcowofizgoc Sma 028226 3.58.5 «bank 35me 1K NNOnE G— mmmmovflwz omfiommcg SEEM 03.3on Mani £5».me *QNVNQai on _ caged—Z £88m mung—>85 aficiounbo—z “68% SEEN» Dues 9; SmNV .cm mvmanImZ U 5805 £35585 “Bowman ”Santa—o2 8.22:» uzchEm. #032: mi mmaoclmz omfiomagafiofi 89.2320 Mama 55me x23 5: an .oommwm ”851.2 8805:“ 825:6 Mean esrfi _. ER 9: 63 5mm £8932 z 88.53%; Eurasia sagas» 2&3 n: so: 50> £03.32 aeaégbog 28558< ”can 552“: Es 3: 69 6mm wemmhlmz 520 8:80 .E 322:»? <29 :8 pingam .93» m: an .8: NESME 33m: “Ease“--fia 322% EU 82%» §~§§£ 33$ N: $832 89:5 E cowoqfigeoge: §§53§ sagessm as»: :2 woficeumz seen 95. ”ER 5:23. 9% o: Avmv in 8331.? E8388 5 08536306 88,3?— uotSzm §m=c£~cm Mme: m2 Ava .2 232$ €33 gauge»; 25.5 fig :3me 95 ”2 $3 .2 38.32 aéugées 2653 “53 35$» 229: 5 $3 .2 356.32 aefiagzseoaea $5.5 Etafiib §§§§ NE 92 cé ,2 3832 8283: £825 €0-38 2653 mean 35$: 38% 9: as as. mama was 32.5 3:338 «ENE §.§§8§ .931 v2 aversion .o: .oo< I .o: $325.5 Ema—00 92:“QO 535 :5:qu0 ewzmb E 383 N 35¢ 58 “Hypothetical and unknown genes were not listed. bGene name designations were based on originally-reported gene products that shared high similarity to pfi clone sequences. Asterisks afier the gene name indicated that two or more clones were identified for the same gene in the experiment. cPredicted proteins or functions based on similar proteins identified using BlastX searches. Genes identified to be induced in other bacteria-plant or animal systems using similar IVET screen or other in vivo expression systems. Rso: Ralstonia solanacearum; Ech: Erwinia Chrysanthemi; Pst: Pseudomonas syringae pv. tomato; Pfl: Pseudomonas fluorescens; Sen: Salmonella enterica; Sty: Salmonella typhimurium; Sau: Staphylococcus aureus, Apl: Actinobacillus pleuropneumoniae, Vch: Vibrio cholerae, Eco: Escherichia coli. Number indicated the reference. 59 Quantitative expression analysis of selected pear up-regulated genes To verify that the pear-upregulated gene promoters identified using the qualitative IV ET assay are induced in pear, quantitative GUS activity for six strains containing pfi promoter constructs and for promoter constructs containing the dspE promoter (in both directions) was monitored 24 and 48 hours after infection of immature pear. The positive control dspE f“ promoter in pZYF 8 was highly induced in pear after infection at both 24 and 48 hr post inoculation (Table 3), whereas the negative control dspE" promoter in pZYF2 showed very low GUS expression (Table 3). Most of the pfi clones tested showed varying degrees of induction of promoter activity at 24 h and 48 h after infection of immature pear (Table 3). The pfi 43 clone (oppA) was found to be highly induced at both 24 and 48 h after inoculation; whereas pfi 5 (hothoCEA), pfi 9 (ytsIIJ), pfi 91 (rlsA), and pfi 93 (grrS) were only induced at 48 h post inoculation. The clone containing hypothetical protein (sav2932), on the other hand, was strongly induced at 24 h post inoculation with expression tailing off at 48 h (Table 3). 6O TABLE 3. Expression of IV ET clones afier inoculation of immature pear fi'uit pfi clone Gene 0 hpi" 24 hpi" 48 hpi" or strain homolog GUS GUS Fold GUS Fold activityb activityb inductionc activityb inductionc PZYF2 dspE r" d u.l.f u.l.f O u.l.f O pZYF8 dspE for e 1.0 i 2.1 282.0 :1: 145.5 282.0 118.3 :t 18.4 118.3 pfi5 hothoC 10.8 :t 6.7 11.8 :1: 4.4 1.1 90.1 i 43.3 8.4 pfi9 ytsIIJ 37.1 d: 2.0 74.7 :1: 34.9 2.0 128.0 d: 48.8 3.5 pfi43 oppA 1.0 d: 0.9 91.4 :t 27.4 91.4 138.9 :1: 25.2 138.9 pfi9I rlsA 1.0 a 1.3 11.1." o 19.2 :1: 12.6 19.2 pfi93 grrS 77.7 i 2.3 30.2 i 7.0 0.4 106.3 i 35.8 2.4 sav2932g 107.4 :1: 3.5 852.2 d: 380.3 7.9 86.2 i 40.4 0.8 “Hours post inoculation (hpi). bGUS activity is shown in mol of 4-methylumbelliferyl produced min1 109 CFU". Data represented the mean of three measurements i standard error. Similar results were obtained in a second independent experiment. CF old induction is shown as GUS activity at 24 or 48 hpi/GUS activity at 0 hpi. ddspE promoter in opposite direction of uidA gene. edspE promoter in correct direction of uidA gene. u.l.: Under detection limit of vector control. gsav2932 encodes a hypothetical protein and is not listed in Table 2. Construction and analysis of knockout mutants Although the IVET experiments were conducted in the E. amylovora Eal 10' background and an Ea110 hrpA mutant was available, we were unsuccessful in subsequent attempts to construct chromosomal knockout mutants in this strain. Thus we utilized the strain Ea1189 for mutant construction. Although the similarity of genetic backgrounds of these strains is currently unknown, the virulence of the two strains is similar (data not shown) and the overall genome diversity of E. amylovora is relatively low and therefore, we hypothesize that the expression of promoters identified as upregulated in Ea110 would be comparable to that in Ea1189. We chose the hothoCEA and mltEEA genes as candidates for insertional mutagenesis using allele marker exchange 61 to investigate the potential roles of those genes in virulence. The full sequences of these genes and their corresponding upstream and downstream sequences were obtained by various methods including fully sequencing available clone inserts or by recovering additional flanking DNA sequences using thermal asymmetric interlaced (TAIL-) PCR. Sequence analysis showed that the deduced amino acid sequence of the hothngA gene shared 77% similarity with that of hothonsr gene from P. syringae pv. tomato (data not shown). The deduced amino acid sequence of E. amylovora mltEgA gene showed 75% similarity with that of mltEyp gene from Yersinia pestis (data not shown). As described in Materials and Methods, we were able to generate insertional mutants for hothngA and mltEgA genes and tested the effect of these mutations on pathogenesis and bacterial growth in immature pear. The hothoCEA mutant ZYC1-3 was fully virulent like the WT Ea1189 strain upon infection of immature pear (Figure 3A). The symptoms caused by the hothngA mutant also progressed similarly to the wild type strain and caused typical symptoms, i. e. necrotic lesions and the production of bacterial ooze (Figure 3A). In addition, we were unable to detect significant difference in bacterial grth upon infection of immature pear for the hothngA mutant (Figure 3B). In contrast, the mltEEA knockout mutant (ZYE3-11) was slightly reduced in virulence and bacterial growth (Figure 3A and 33). Although the mutant caused similar symptoms as the wild type strain, symptom progression was slightly reduced (Figure 3A). There was no difference for bacterial growth of the mltEEA mutant compared to that of the wild type strain two days after inoculation, however, cell counts of the mltEgA mutant were 5 to 10 fold less than that of the WT strain at three and four days afier inoculation (Figure 3B). 62 WT ZYC1-3 fi ZYE3-ll W B 12 11‘ 10‘ i 0 9‘ a —o—E31189 .3 81 —-—zvc1-3 an . +ZYE3-11 D 7 in U 6‘ M 6 1-3 5‘ 4 3 v r x 0 1 2 3 4 Days post-inoculation Figure 3. Symptoms and growth of Erwinia amylovora WT Ea1189 and corresponding hothoCEA and mltEEA mutants in immature pear. (A) Symptoms caused by Erwinia amylovora Ea1189 and hothoCEA (ZYC1-3) and mltEEA (ZYE3-11) mutants in immature pear. WT: Wild type strain; W: water control. DPI: days post inoculation. (B) Bacterial growth of Erwinia amylovora WT Ea1189, and hothoCEA (ZYC1-3) and mltEEA (ZYE3- ll) mutants during infection of immature pears. The growth of strains was monitored at 0, 1, 2, 3, 4 days afier inoculation. Data points represent the means of three replicates 1- SE. Similar results were obtained in a second independent experiment. 63 DISCUSSION We utilized a simplified qualitative IVET approach to scan the E. amylovora genome and recovered 394 miique chromosomal genes with increased expression during infection of pear fruit tissue. As expected, this study highlighted the importance of type III secretion in E. amylovora pathogenesis with the recovery of genes encoding regulatory and structural components of the Hrp type III secretion system and effector proteins. While we did not recover all of the currently known hrp-regulated genes in E. amylovora, our results are similar to those of other IVET studies with plant pathogenic bacteria. For example, IVET studies of E. chrysanthemi and P. syringae pv. tomato identified two and eight hrp-regulated genes, respectively (24, 261). These findings validated our approach and suggested that a detailed analysis of the genes recovered in this study would further reveal additional determinants involved in the pathogenesis of the fire blight bacterium. The dspEF operon, encoding the major effector and pathogenicity factor DspE and its cognate chaperone DspF, was recovered multiple times in our analysis and shown by quantitative expression analysis to be highly expressed during pear infection (Table 3). The importance of DspE and its homologs to plant pathogenesis is well known in a number of pathosystems (26, 77, 143, 170, 243) although the elucidation of the function(s) of this large protein remains unknown. DspE was recently shown to contribute to the suppression of salicylic-acid-mediated basal immunity (57); effector suppression of the host defense response is rapidly becoming recognized as an important strategy of bacterial plant pathogenesis (5). We identified a new putative effector 64 HothoCEA in this study, an ortholog of HothoC from P. syringae pv. tomato (224). As with many effectors from P. syringae, a knockout mutant of HothoCEA in E. amylovora Ea1189 was not reduced in virulence presumably due to functional redundancy with other effectors in the E. amylovora genome. The other known E. amylovora effectors HrpN and HrpW were not identified as upregulated in this study; although the roles of hrpN and hrpW in the pathogenicity of E. amylovora were reported to differ during infection of immature pear fruit (126, 252). It is tempting to speculate that additional effector proteins may exist in the genome of E. amylovora and contribute to the virulence of the bacterium. The importance of type H secretion in E. amylovora pathogenesis was also highlighted with the identification of the upregulation of genes of the ytsIIJ operon and components of the general secretion pathway. Type H secretion is a c00perative process initially dependent upon the secretion of enzymes into the periplasm by the general secretion pathway followed by targeted secretion through the type II apparatus (9, 220). In Y. enterocolitica, the Ytsl protein secretion apparatus is unique to highly pathogenic species, is important for virulence in a mouse model, and shares homology with type II secretion clusters from E. chrysanthemi and E. carotovora (17, 112). Peh (polygalacturonase), an enzyme thought to be secreted by the TZSS, was also upregulated and recovered in our IV ET screen (122). While the importance of polygalacturonase to virulence in soft-rotting Erwinia spp. is well known (122), the role of cell wall degrading enzymes in E. amylovora pathogenesis is currently still unknown. In addition, the upregulation of MltE, a specialized cell-wall-degrading enzyme was interesting in that the function of this enzyme is to generate localized openings of the bacterial 65 peptidoglycan envelope for export of bulky materials including possibly toxins and fimbrial proteins and to allow the efficient assembly and anchoring of supramolecular transport complexes such as T2SS and TTSS in the cell envelope (63, 130). As in P. syringae (24), we found that E. amylovora MltE made a small contribution to virulence. We identified three additional upregulated enzymes in our IV ET assay that are potentially secreted from the cell including levansucrase Lsc and the adhesin-like protein HecA which belongs to a class of external virulence factors that is widely distributed among plant and animal pathogens. HecA from E. chrysanthemi contributes to attachment, aggregation, and epidermal cell killing and is thought to be involved in the earliest stages of E. chrysanthemi pathogenesis (212). Levansucrase, an enzyme that directs the synthesis of levan from sucrose, has a known effect on the virulence of E. amylovora during pear seedling infection (85). The PrtA metalloprotease contributes to E. amylovora virulence in an apple leaf infection assay and is apparently dependent upon the type I Prt machinery for secretion (122, 265). These results demonstrate the importance of TTSS, T2SS and of other external virulence factors in E. amylovora infection of fruit tissue. A total of 5.3% of the IVET genes identified were placed in the functional category of stress response including genes involved in the response to reactive oxygen species, both heat and cold shock, and carbon and sulfate starvation. E. amylovora apparently induces an initial host defense response early after infection (247, 248); the bacterium is capable of surviving this plant oxidative burst with the initial plant cell death and nutrient leakage thought to provide the impetus for fiirther spreading of the pathogen within the plant. The role of individual proteins in oxidative stress survival is currently 66 unlmown in E. amylovora, however, the alkyl hydroperoxide reductase AhpC is a known virulence factor in several plant pathogenic bacteria contributing to protection from oxidative stress from an active plant defense response (169). We recovered a multitude of transporters frmctioning in the uptake of iron, sugars, amino acids, and inorganic ions. The induction of these systems during infection indicates that E. amylovora elaborates various factors as needed to colonize host tissues. Iron availability is critical to most bacterial pathogens and the siderophore desferrioxamine is a virulence factor in E. amylovora (73). We recovered three upregulated proteins involved in iron transport or storage. It is probable that under unfavorable conditions, the bacterium itself adjusts and overcomes nutrient and iron deficiencies. Several upregulated transport proteins recovered were ABC transporters which is potentially significant because ABC transporters both directly and indirectly affect virulence of bacterial pathogens (80). While most of the transporters were involved in uptake, the multidrug resistance proteins EmrB (pfi 31) was also upregulated and presumably functions in the efflux of plant-derived toxins encountered during infection. The role of the AcrAB efflux pump in E. amylovora virulence and tolerance of phytoalexins including phloretin, naringenin, and quercetin was recently reported (3 8). Thus, it is possible that many of these ABC transporters are involved in the virulence of E. amylovora. In conjunction with the number of upregulated transporters found, a large proportion of the genes identified in this study were involved in metabolism (20.3%) and nutrient acquisition (16%). These fi'equencies may be associated with the host tissue (immature pear fruit) chosen for study; however, a number of genes we identified were 67 also identified in other IVET studies involving E. coli, P. fluorescens, P. syringae, R. solanacearum, and V. cholerae (Table 2; 24, 36, 150, 261). About 12% of the genes identified in this study were involved in regulation, which is a similar ratio to that identified in an IVET examination of E. chrysanthemi infection (261). Previously-known E. amylovora transcription factors that were upregulated included RcsA, an activator (along with RcsB) of amylovoran production (251) and RlsA, an activator of levan production (264) along with the capsular polysaccharide export protein KpsC. This further confirms that the production of both amylovoran and levan in E. amylovora is induced during infection. Another important regulator, GrrS (global response regulator sensor in a two-component regulatory system), is a homolog of GacS which, along with GacA, globally regulate a network that controls exoenzyme and secondary metabolite (toxin) production in Pseudomonas spp., virulence functions in E. carotovora, and also regulates the production of quorum-sensing signalling molecules (45, 72, 203). GacA/GacS-regulated networks also function by positively controlling the transcription of small regulatory RNAs, transcriptional activators, and alternative sigma factors such as HrpL (45, 100). In E. amylovora, the small regulatory RNA rsmB titrates the repressor RsmA in a system that affects exopolysaccharide production and therefore, pathogenicity (147). EnvZ is the sensor component of the OmpR/EnvZ two component regulatory system that is very important in regulating various cellular components such as outer memberane proteins OmpC and OmpA, which is also upregulated in this study. In Salmonella spp., OmpR-EnvZ regulates another two component system SsrA-SsrB that in turn regulates the type III secretion system produced by Salmonella pathogenicity island 68 2 (Spi-2; 133). EnvZ is a transmembrane sensor that predominantly responds to acidic pH conditions and subsequently phosphorylates OmpR which functions as a transcriptional activator in the expression of the ssrAB genes (79). SsrA is second sensor protein that is responsive to acidic pH and also detects low osmolarity conditions and the absence of Ca2+ ions, all environmental conditions within macrophages where the Spi-2 type III secretion system is exclusively expressed (79). In E. amylovora, the structural components of the TTSS encoded by the Hrp regulon are regulated by the two component system HrpX and Her, which direct the expression of the GSA-dependent, enhancer- binding protein HrpS (251). Both Her and HrpS function in activating the expression of the alternate sigma factor HrpL, thereby regulating the various genes and operons of the Hrp regulon which contains HrpL-dependent promoter sequences (251). The expression of HrpX and HrpS is regulated by low pH, low nutrients, and low temperature conditions, mimicking the plant apoplast, but also representing conditions that suggest a two component regulatory system such as OmpR-EnvZ could further regulate the hrpXY operon despite no direct evidence to support this claim. Interestingly, both hrpX and hrpL, along with EnvZ were found to be upregulated during infection of immature pear in this study (Table 2). Among the bacterial cell surface and transmembrane upregulated proteins, three flagellar proteins FliG, FliM, and FlgN were upregulated. The trait of motility is not required for E. amylovora pathogenesis, however, motility does increase blossom infectivity, particularly at lower cell concentrations (16). Furthermore, a homolog of Y. pestis FliZ, a positive regulator for the flagellar biosynthetic operon and an alternative sigma factor, was also found to be upregulated in our study. In Salmonella enterica 69 serovar Typhimurium, F liZ upregulates HilA, which in turn activates production of several invasion proteins encoded within the Salmonella pathogenicity island 1 (113). Finally, the contribution of cell shape to virulence was also highlighted by the recovery of an E. coli RodA homolog; E. amylovora mutants with TnPhoA insertions within the rodA-ppr operon were previously reported to be avirulent (165). In summary, our IVET screen successfully identified a variety of genes upregulated during fruit infection by E. amylovora. We utilized a modified IVET method in this study which is different from many other IVET studies in that we did not impose a rigorous selection step, i.e. one that necessitates rescue of an essential phenotype, in our gene identification work. An advantage of our approach is that we screened clones individually which would remove potential competition among firsion clones in vivo if inoculated as a pool. The successful identification of a large number of known virulence genes of E. amylovora in this study further validated our approach. However, because of the qualitative nature of the gene identification step, through B-glucoronidase staining and visualization of gene expression on pear slices and agar medium, it is possible that this methodology may have resulted in some artifacts. Nevertheless, the main goal of this work, as in other IV ET analyses, was to identify potentially important genes in the E. amylovora infection process that could be subjected to further detailed studies to clearly delineate the role of these genes in pathogenesis. 70 We further confirmed that the TTSS and its major effector protein DspE are essential for full virulence in E. amylovora during infection of immature pear. We also found a complete and functional T2SS and its potential secreted proteins in E. amylovora for the first time. We identified a new putative effector, external virulence factors such as HecA which were previously unknown in E. amylovora, and discovered a number of putative regulatory proteins that may influence the regulation of virulence factors on a global level and eventually contribute to the virulence of the bacterium. We can now ask questions concerning the comparative regulation of critical genes identified in this study during infection of other host tissues, particularly blossoms and shoots. It is possible that E. amylovora the pathogen may utilize differential virulence strategies depending upon the host tissue encountered. Of interest to us also is the expression profile of these same genes during infection of highly susceptible versus fire blight-tolerant apple varieties. 71 CHAPTER 2 Taking a page from macerogenic bacteria: Functional characterization of an Erwinia amylovora endo—polygalacturonase and its effect on virulence 72 ABSTRACT In a previous in vivo expression technology (IVET) screen of induced Erwinia amylovora genes during infection (Chapter 1), two genes associated with type II secretion (T2SS), the pseudopilins outIJ and an endo-polygalacturonase peh were identified. Here I describe the contribution of TZSS to virulence, using deletion mutants of outDEA and peh. In particular, E. amylovora peh was reduced in virulence during apple leaf inoculations when cells were not introduced into the main vascular tissue of leaves. The activity of both the T2SS operon promoter and peh promoters were also determined under differing conditions using flourometric methods. Additionally, the peh gene was confirmed to encode a protein with polygalacturonase activity. His-tagged polygalactuonase was present in the supernatant of E. amylovora, demonstrating secretion to the extracellular milleu. Overall, this is the first report of a polygalacturonase enzyme influencing virulence of E. amyolvora. 73 INTRODUCTION Erwinia amylovora is a Gram-negative phytopathogenic bacterium that infects members of the Rosaceae family. This necrogenic pathogen causes tissue necrosis and additional symptoms including water soaking, wilt, and cankers. Infection can occur either by entry through nectarthodes in the base of flowers, or through wounds on young, actively growing tissue. Entry into the plant via wounding is followed by multiplication in the cortical parenchyma, and the expression of disease promoting genes (246). Ultrastructural analyses have shown bacteria-dependent plasmolysis by 24 hours post infection, of xylem parenchyma. As the infection progressed, the extensive plasmolysis appeared to lower the turgor pressure around xylem vessels resulting in the bursting of cells (90). Later studies have demonstrated that infection and symptom development depends on the secretion of extracellular proteins from E. amylovora. Major pathogenicity factors of E. amylovora include a type III secretion system (T3SS), the type 1H effector DspA/E, and the exopolysaccharide amylovoran (11, 18, 26, 27, 82, 127, 239). Secretion and translocation of DspE and an additional effector HrpN through the T3SS causes tissue collapse in susceptible hosts, and a hypersensitive response in resistant hosts and nonhosts (246). Amylovoran is hypothesized to protect E. amylovora from the oxidative burst in host tissues and also to mask the presence of E. amylovora from its host (128). Virulence factors also contribute to infection by E. amylovora including the iron- binding siderophore desferrioxamine (60) and genes conferring the utilization of sorbitol, the major sugar present in apple (3), galactose (161) and sucrose (29). E. amylovora also 74 secretes other virulence factors that affect interaction with its host, including the recently- described type III effector AerptZEa (267), metalloprotease (PrtA) (265) and the levan biosynthetic enzyme levansucrase (20, 85). E. amylovora also possesses a full sec transport system which supports other secretion systems such as type II, type IV and type V (auto-transporter) secretion (188, 238, 266). Type H secretion (T28) is an important virulence factor in gram negative animal pathogens and phytopathogens such as Pectobacterium carotovorum, Erwinia chrysanthemi, Ralstonia solanacearum, and Xanthomonas campestris. The TZS operon was identified first from Klebsiella oxytoca by transferring a plasmid encoding the T28 operon from K. oxytoca to Escherichia coli to demonstrate that T2S proteins were required for the secretion of the enzyme pullulanase (215). Transfer of a cosmid containing the T2S from E. chrysanthemi to E. coli was also used to demonstrate the ability of the T2S to secrete various extracellular enzymes (138). Genome sequencing has made it possible to identify the presence of T28 genes in many gram negative bacteria other than pathogenic bacteria (50). However, the comparison of known functional T28 systems highlights 12 to 15 proteins expressed from the operon are required for secretion across the outer membrane (215, 219, 220). Although the gene naming convention has differed slightly amongst species possessing a T2S system, the majority of the gene products follow the Klebsiella gene designation as A through 0 plus S (50, 219, 220), with proteins CDEFGHUKLMNO acting as the core of 13 essential T28 components (50). Here the general nomenclature of the T28 operon, “general secretion pathway” (Gsp) will be followed. 75 Exoproteins secreted by phytopathogenic bacteria include proteinases, amylase, pectinases, cellulases and polygalacturonases (50). Polygalacturonases present in phytopathogenic bacteria fall into two categories, endo- and exo-polygalacturonases. The former catalyze the random hydrolysis of pectic acid forming oligogalacturonates while the latter catalyze the terminal hydrolysis of pectic acid forming galacturonic acid monomers (115). While endo-polygalacturonase activity has not been detected in E. chrysanthemi (109), endo-polygalacturonase activity is a known virulence factor in P. carotovorum (216), R. solanacearum (61, 107) and Agrobacerium vitis (210). Analysis of the secretome of E. chrysanthemi identified 14 T28 secreted proteins, including 11 pectinases, cellulase, rhamnogalacturonan lyase, a novel esterase and a novel Avr-like protein (122). Sequence analysis of Ralstonia solanacearum has identified 6 exoproteins including three polygalacturonases, a pectin methylesterase, and two cellulases (139). Previously, E. amylovora was described not to possess any cell wall degradation abilities (226); however, recent studies have described the presence of a type II secretion gene, outF, based on Southern hybridization to an outhca homologue and weak CelA activity (207). We have previously described the upregulation of outIJEA and peh expression upon colonization and infection of immature pear tissue (266). The pseudopilins outIJEA share the highest homology to ytsIIJ from the human pathogen Yersinia enterocolitica. In addition, the identification of peh, an endopolygalacturonase, indicates that E. amylovora may use a polysaccharide degrading enzyme during host infection. In order to determine the contribution of T2SS to virulence, deletion mutants of the porin gene outDEA and polygalacturonase peh were created. Additionally, the expression of uidA 76 fused to the T2SS operon promoter and peh promoter was determined under differing conditions. 77 METHODS AND MATERIALS Bacterial strains and plasmids. The bacterial strains and plasmids used in this study are listed in Table 4. E. amylovora and E. coli were cultured on Luria-Bertani (LB) medium at 28°C and 37°C, respectively, unless otherwise indicated. Minimal media was comprised of 1XM9 salts according to Sambrook et al. (217) with the addition of thiamine (0.02% w/v), niacin (0.02% w/v), and 2% (w/v) of one of the following: fructose, galactose, glucose, sorbitol or sucrose or 0.25% (w/v) pectin derived from apple fruit (82, 155). Liquid SOC medium used for transformed cell recovery was made according to Sambrook et al. (217). Antibiotics were used at the following concentrations in medium: ampicillin (Ap), 100 rig/m1; carbenicillin (Cb), 50 ug/ml; chlorarnphenicol (Cm), 25 [Lg/ml; gentimicin (Gm), 10 jig/ml and kanamycin (Kn) 25 jig/m1. Oligonucleotide primers used for PCR were purchased from Invitrogen (Carlsbad, CA) and are listed in Table 4. Molecular biology techniques. Genomic and plasmid DNA isolation, gel purification of DNA, and PCR fragment clean up were performed using Qiagen kits (Qiagen Inc., Valencia, CA) according to manufacture’s recommendations. PCR amplification, restriction enzyme digestion, T4 DNA ligation, and cloning were all performed according to standard molecular biology protocols (217). Sequencing of DNA was performed by the Research Technology Support Facility (RTSF) at Michigan State University. Database searches were conducted using BLAST programs (7) at the National Center for Biotechnology Information (NCBI) and the E. amylovora genome sequence at the Sanger 78 Institute (http://www.sanger.ac.uk/Projects/E_arnylovora/). Amino acid alignments were conducted using ClustalW (v1.83) from the European Bioinforrnatics Institute (http://www.ebi.ac.uk/services/index.html). Signal peptide determination was done by the SignalP v3.0 program (http://www.cbs.dtu.dk/services/SignalP/) using both neural networks and hidden Markov model predictors of signal peptides (71). Construction of the peh and outD deletion mutants. Construction of stable E. amylovora mutants was as previously described (267) using the A phage Red recombinase system (56). Eall89(pKD46) was grown up overnight and diluted 1/10 in fresh LB medium supplemented with ampicillin and 10 mM L-arabinose and grown to an ODwo of ~ 0.6. The cells were then pelleted (5,000X g) and washed three times in 10% glycerol at 4°C. After the final wash, the cells were resuspended in l/ 100 volume. The primer sets OutDmut-F/ OutDmut-R and Pehmut-F/ Pehmut-R were used to generate PCR fiagments that were transformed into 100 pl of competent Ea1189(pKD46). Cells were resuspended in 500 pl of SOC medium without antibiotics and allowed to recover overnight at 28°C. Transformants were then plated on solid LB medium supplemented with Cm and Ap. Colonies that were resistant to both Cm and Ap were selected for PCR screening using primer pairs that amplified from the cat cassette to the up and downstream sequences of outD (le/outD-F and Cm2/outD-R) or peh (le/peh-F and Cm2/peh-R). The primer pairs outD-F/outD-R or peh-F/peh-R were used to confirm the insertion of the cat cassette. 79 TABLE 4. Bacterial strains, plasmids and primers Strains, plasmids Source or and primers reference Erwinia amylovora Eal 10 Wild type, isolated from apple Ea1189 Wild type, isolated from apple (3 8) 8B] As Ea1189 but also outD::Cmr This study SB2 As Ea1189 but also peh::Cmr This study Escherichia coli DHlOB F mcrA A(mrr-hstMS-mchC) 801acZAM15 Invitrogen AlacX74 recAl endAl araA139 A(ara, leu)7697 galU galK 7t - rpsL (Str') nqu TGl Stratagene Plasrrrids pBluescriptII(+) Ap’, cloning vector Invitrogen pGEMT-easy Ap’, PCR cloning vector Promega pQE60 Ap’, Qiagen pBBRl ~MCS2 Km', cloning vector derived from pBBRlMCS (131) pKD46 Ap’, PBAD gam bet exo pSClOl oriTS pKD3 Ap’, FRT cat FRT PSI PS2 oriR6K rng pGCMO GrnR cassette with downstream transcriptional terminator (266) and gusA with upstream translational stop codons in pGem32f+ pZYF2 570bp region of dspEF promoter in SmaI site of pGCMO, (266) in opposite orientation relative to uidA pZYF 8 570bp region of dspEF promoter in SmaI site of pGCMO, (266) in correct orientation relative to uidA pSEB28 863bp region of outD“ promoter cloned upstream of uidA This study in pGCMO pSEB29 546bp region of peh promoter cloned upstream of uidA in This study pGCMO pSEB30 477bp region of IS] 133 promoter cloned upstream of uidA This study in pGCMO pSEB3l 1248bp NcoI/ BglII fragment of peh ORF in pQE60 This study pSEB36 1251bp SacII/XhoI peh ORF in pBBRl-MCS2 This study pSEB37 193 8bp SacII/Xhol outD ORF in pBBRl-MCSZ This study Primers OutDmut-F ATGAAGAAGAGATCCCCCCAATCCGCTACCAGCATCCGGCGGCT GTTACC GCGA TT GT GT A GGCT GGA GCT OutDmut-R TTACGCT’I'ITIT CCGGTAGAAATCAGCGATGTGCTT CT GTATTTC CACCA CCA T GGT CCA TA T GAA TA T CCT CC Pehmut-F ATGACTATTCTAACCAATTGTITATTAAGGTTT'ITATACCTCGGA ATATC GCGA TT GT GT A GGCT GGA GCT Pehmut-R TTAC'ITATCGAT'I‘TGAACGTTGTTAACGTTATAT'I'I‘TGTCGCAGG ATCAA CCA T GGT CCA TA T GAA TA TCCT CC outD-F CAACGGATCATCATTCGTCACC outD-R GCACGTI'ATGCTGTTGTACC peh-F GCCCCCTGTTAAACCCGTACTCAATAAC 80 TABLE 4 @on’t) Strains, plasmids Source or and primers reference peh-R CCTGGCCACGGTGATAGAAG'I‘TTTG le TTATACGCAAGGCGACAAGG Cm2 GATCTTCCGTCACAGGTAGG Orflprom-F GTCT CCGGTACCTTGCGGATAATACCTACCGCAACCT G (Kpnl) Orflprom-RZ CGCACTGCAGTTTCTTTGITCCTTATGATGTCTCCG (PstI) Pehprom-F GTCTCCGGTACCCCGCCATTGGCCGTCTTITAATCGCCG (KpnI) Pehprom-R2 CGCACT GCAGTTGAGATCCT CT GTGTAACT GG (PstI) dspEF-F TCCCCCGGGCAGTGAGGGGGGGCAGACTTTITI'ITAACC (SmaI) dspEF-R TCCCCCGGGTATCI'TCGCCGCTGCCACCT'ITCACCATTG (SmaI) ISl l33-F CCGGTACCCGTCGCGTGA'ITGGCTGG (Kpnl) Aj 1655 CGCACT GCAGGAAGCGCGGAGGTGGCT C (PstI) PehpQE-F CATGCCATGGATGACTATTCTAACCAA’ITGT'ITATTAAGG (NcoI) PehpQE-R GAAGATCTCTTATCGATI'TGAACGITGTTAACG (BglII) PehMCSZ-F ACT GATCCGCGGATGACTA'ITCT AACCAATTGTITATTAAGG (SacII) PehMCSZ-R TCTCGAGTTACI'I‘ATCGA'I'ITGAACG’ITGTTAACG (XhoI) OutDMCSZ-F ACTGATCCGCGGATGAAGAAGAGATCCCCCCAATCCGCTACC (SacII) OutDMCSZ-R TCTCG_1_A§TTACGCTI‘TI'ITCCGGTAGAAATCAGCGATGTGC (XhoI) Creation of uidA reporter strains. Fusions of the putative type H out operon promoter, peh promoter and dspEF promoter to the uidA gene were constructed using pGCMO. The PCR primer pairs Orflprom-F/Orfl prom-R, Pehprom-F/ Pehprom-R and 181133- F/ajl655 were used to amplify the upstream sequence fi'om the type II operon, peh and IS1133 respectively. Purified upstream sequences were inserted into the KpnI and PstI sites in pGCMO and ligated overnight at 14°C overnight using T4 ligase fiom Invitrogen (Carlsbad, CA). Promoter constructs were transformed into chemically competent E. coli DHlOB, for screening using PCR and restriction enzyme analysis. Confirmed orflpzszCMO, pehp::pGCMO and IS1133p::pGCMO were named pSEB28, pSEB29 and pSEB30, respectively. Positive colonies were stored in LB plus 10% glycerol at -80°C, 81 and the respective promoter fusion constructs were transformed into electrocompetent Ea1189. Fluorometric assays with promoter fusion constructs. The activity of the promoters fused in front of the uidA gene was measured using the fluorescent substrate 4- methylumbelliferyl-B-D-glucuronide (MUG) as was done for the E. amylovora IVET screen (266). The promoter fusion constructs pSEB28, 29 and 30 and pZYF2 and 8 were grown in l><106 CFU/ml) was then used to 83 dip sterile scissors into and cut either across the midrib or roughly 1 cm into interveinal tissue of the second emerged leaf. Pear infection assays used immature pear (Pyrus communis L. cv. ‘Bartlett’) according to previous studies for the analysis of E. amylovora mutants (26, 155, 267). Immature pears of similar size were sterilized in 10% bleach solution for 15 minutes and washed three times with sterile dHZO. After drying in a sterile flow hood, the pears were pricked with a sterile 30 1/2-gauge needle and 3 ul of ~5><105 CFU/ml E. amylovora was pipetted onto the wound. Infected pears were kept on sterile wet paper towels in sealed containers in a 28°C incubation chamber. For statistical analysis of aggressiveness of the E. amylovora outD mutant SBl, peh mutant SB2 and wild type Ea1189 the method of Cabrefiga and Montesinos was used (40). E. amylovora cultures used for inoculation were grown overnight in liquid LB medium supplemented with the appropriate antibiotics. Cultures were then spun down at 5000 x g and washed with 0.5x PBS twice and diluted to an 013600 of 0.8 (~5x109 CFU/ml). Cultures were then serially diluted 1/ 10 in 0.5>< PBS to obtain cultures that were logro 7.5, 6.5, 5.5, 4.5, 3.5, 2.5 and 1.5. Six immature pears were inoculated 4 times opposite of each other for each concentration by pipetting 1011.1 of the bacterial dilution on each wounding site. The presence or absence of necrosis and/or ooze was recorded for each inoculation site every day following inoculation at approximately the same time as inoculations. These experiments were repeated three times. For the calculation of the median effective dose, the hyperbolic equation as described by Cabrefiga and Montesinos (40) was applied: y=ymax(x/(X+Kx)) (1) 84 where y is the incidence of infections, ymx is the maximum of infection incidence, x is the concentration of the culture and K, is the median effective dose (EDso). For the cultures used, the EDso was solved for using data collected on day 3 of infection. Aggressiveness of the E. amylovora cultures was determined using a modified Gompertz equation (40): y = K — exp(-Bg° exr)[-rg(t - to)]) (2) where K is the curve asymptote signifying the maximum incidence of disease, 38, accounts for the origin of the curve, and to is days of delay until symptoms are observed. For the calculation of aggressiveness, data from all days post inoculation were used for the culture diluted to 5><103 CFU/ml. The statistical program IMP version 6.0 (SAS Institue, Inc., Cary, NC) was used to fit the data collected to the hyperbolic and modified Gompertz equations. Protein expression of Peh. Protein expression constructs were built to express Peh with a carboxyl-terminal 6x histidine tag. The primer set PehpQE-F/ PehpQE-R was used to PCR amplify a 1251 bp fragment of peh, which is missing the stop codon. The peh fragment was inserted into the NcoI and BglII sites of pQE6O and transformed into E. coli TGl for screening using PCR amplification and restriction enzyme digestion. Confirmed peh::pQE6O constructs were named pSEB3l and transformed into Ea1189, both TGl(pSEB31) and Eal 189(pSEB31) were stored in LB freezing medium and stored at -80°C. A fast mini-preparation of protein fractions was prepared according to (52) and used for subsequent Western blots. Eall89(pSEB31) and TG1(pSEB3l) overnight 85 cultures were diluted 1:5 into fresh LB medium supplemented with 50 jig/ml carbenicillin and allowed to recover for one hour, after which IPTG was added to a final concentration of 1mM. Cells were harvested 48 hours post-induction, and centrifuged at 13,000 rpm to separate the supernatant from cells. The pellet was used for a fast fractionation protocol (52) and separated into periplasmic and cytoplasmic/membrane fractions. Purification of E. amylovora Supernatant Protein. For the visualization of E. amylovora secreted proteins and the ruthenium red staining assay, supernatant was harvested fiom cultures grown on solid l>