Microbiology
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Microbiology 155 (2009), 174-185; DOI  10.1099/mic.0.022103-0
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary Table
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Xu, J.
Right arrow Articles by Hegedus, D. D.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Xu, J.
Right arrow Articles by Hegedus, D. D.
Agricola
Right arrow Articles by Xu, J.
Right arrow Articles by Hegedus, D. D.
Microbiology 155 (2009), 174-185; DOI  10.1099/mic.0.022103-0
© 2009 Society for General Microbiology

Comparative EST analysis of a Zoophthora radicans isolate derived from Pieris brassicae and an isogenic strain adapted to Plutella xylostella

J. Xu1,2, D. Baldwin1, C. Kindrachuk1 and D. D. Hegedus1

1 Agriculture and Agri-Food Canada, Saskatoon, SK S7N OX2, Canada
2 Institute of Microbiology, College of Life Science, Zhejiang University, Hangzhou 310029, PR China

Correspondence
Dwayne D. Hegedus
Hegedusd{at}agr.gc.ca


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Zoophthora radicans is an entomopathogenic fungus with the potential to be used as an insect biological control agent. To better understand the mechanisms used by Z. radicans to infect different hosts, we generated expressed sequence tag (EST) datasets from a Z. radicans strain originally isolated from Pieris brassicae, and an isogenic strain passaged through Plutella xylostella. In total, 1839 ESTs were generated which clustered into 466 contigs and 433 singletons to provide a set of 899 unique sequences. Approximately 85 % of the ESTs were significantly similar (E≤e–03) to other fungal genes, of which 69.6 % encoded proteins with a reported function. Proteins involved in protein synthesis and metabolism were encoded by 38.3 % of the ESTs, while 26.3 % encoded proteins involved in cell-cycle regulation, DNA synthesis, protein fate, transport, cell defence, transcription and RNA synthesis, and 4.9 % encoded proteins associated with cellular transport, signal transduction, control of cellular organization and cell-wall degradation. Several proteinases, including aspartic proteinases, trypsins, trypsin-like serine proteases and metalloproteases, with the potential to degrade insect cuticle were expressed by the two isolates.


Abbreviations: EST, expressed sequence tag; Hsp70, 70 kDa heat-shock protein

The GenBank/EMBL/DDBJ accession numbers for the sequences reported in this paper are available in Tables 2–4 and supplementary Table S1.

A table of annotations and functional categories of Zoophthora radicans ESTs from both original and new host-passage isolates is available with the online version of this paper.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Zoophthora radicans (Brefeld) Batko (Zygomycetes: Entomophthorales) is an entomopathogenic fungus capable of infecting and killing a variety of insects. Unlike other entomophthoralean fungi, like Pandora delphacis and Pandora neoaphidis (Xu & Feng, 2000Down) with narrow host ranges, Z. radicans is best known for its broad geographical and host range. While the fungus is pathogenic on species from more than 10 insect orders, including Coleoptera, Lepidoptera, Diptera, Homoptera, Hymenoptera, Heteroptera, Orthoptera, Plecopetera, Psocoptera and Thysanoptera (Humber, 1992Down; Papierok et al., 1984Down; Furlong & Pell, 2001Down; Wraight et al., 2003Down), individual isolates exhibit differing degrees of host specificity. Generally, isolates exhibit greater virulence against the species from which they were originally isolated than against species from other insect orders (Glare, 1988Down). Repeated passage through an insect host is known to enhance the virulence of entomopathogenic fungi toward a specific host, while continual propagation on artificial medium leads to reduced virulence (Brownbridge et al., 2001Down; Quesada-Moraga & Vey, 2003Down); however, the mechanisms underlying these phenomena remain unclear. To identify factors involved in host specificity, we studied a strain of Z. radicans from an isolate originally associated with the cabbage butterfly, Pieris brassicae, and an isogenic strain that was adapted to be highly virulent toward the diamond back moth, Plutella xylostella, through serial infection and reisolation (Li et al., 2004aDown). Previous investigations found that the most obvious difference between the Pieris brassicae isolate and the isolates adapted to Plutella xylostella was the profile of their respective extracellular proteases (Xu et al., 2006Down).

Penetration, dissemination and subsequent death of the insect host involve several biochemical and physical activities, and therefore many genes are associated with fungal pathogenesis. The infection by Z. radicans starts when the primary or capillary conidia that are forcibly discharged from conidiophores land on the cuticle of a suitable host. These germinate and differentiate to form appressoria by extension and swelling of the germ tube tip. Penetration is brought about by a combination of hydrolytic enzymes and mechanical forces exerted by the protruding penetration peg beneath the appressorium. Once the cuticle has been breached, host defence responses are induced or enabled. Some entomophthoralean species initially produce protoplasts lacking the cell-wall components recognized by insect haemocytes, and proliferate with impunity (Glare & Milner, 1991Down; Pell et al., 2001Down; Shah & Pell, 2003Down).

The limited information on genes responsible for pathogenicity and virulence restricts our ability to make targeted improvements to these fungi, and thus a better understanding of the pathogenicity mechanisms and of host adaptation is required. Expressed sequence tag (EST) studies with Metarhizium anisopliae, Conidiobolus coronatus and Beauveria bassiana generated around 1700, 2000 and 10 000 ESTs, respectively, and suggested they share many of the same pathogenicity genes (Freimoser et al., 2003aDown, bDown; Cho et al., 2006aDown, bDown). Genes encoding enzymes with the potential to degrade host cuticle or to disable antimicrobial peptides, including multiple chitinases, subtilisins, trypsins, metalloproteases, aspartyl proteases, phospholipases, lipases and esterases, were identified. These EST data, combined with those from phytopathogenic fungi such as Blumeria graminis (Thomas et al., 2001Down), Gibberella zeae (Trail et al., 2003Down), Sclerotinia sclerotiorum (Li et al., 2004bDown), Magnaporthe grisea (Ebbole et al., 2004Down), Mycosphaerella graminicola (Keon et al., 2005Down), Fusarium verticillioides (Brown et al., 2005Down), Phakopsora pachyrhizi (Posada-Buitrago & Frederick, 2005Down) and a saprophyte, Aspergillus nidulans (Jeong et al., 2000Down; Sims et al., 2004Down), enrich our understanding of pathogenic mechanisms and will help to delineate key factors for fungal adaptation to new insect hosts. In this study, we generated ESTs from a Z. radicans strain originally isolated from Pieris brassicae, and from an isogenic strain adapted to Plutella xylostella to identify candidate pathogenicity genes and to compare gene expression patterns in different insect hosts.


    METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Fungal isolates.
An isolate of Zoophthora radicans, ARSEF 1342 (Pieris brassicae, Pieridae: Lepidoptera) was obtained from the USDA-ARS collection of entomopathogenic fungal cultures (USDA, Ithaca, NY, USA). ARSEF 1342 was originally isolated from Pieris brassicae in Poland in 1983 and was maintained at very low temperature and rarely subcultured. We reactivated the strain through Pieris brassicae, and conidia were used to infect Plutella xylostella larvae and subsequently reisolated from diseased cadavers. This reisolated strain was designated R1 and the serial infection and reisolation process was repeated another four times, at which time the strain was highly adapted to this host (Li et al., 2004aDown). All isolates were routinely grown on SEMA (80 % Sabouraud dextrose agar, 11.5 % fresh milk, 8.5 % egg yolk in 1 l distilled water) at 10 °C in a photoperiod of 12 h light : 12 h dark (L12 : D12) and subcultured every 3 months.

Bioassay and data analysis.
For bioassays, all strains were grown on SEMA at 15 °C (L12 : D12) for 7 days. The resultant fungal colony was macerated with a needle and inoculated into a 100 ml flask containing 30 ml Sabouraud dextrose broth (SDB; 4 % glucose, 1 % peptone, 1 % yeast extract and 12.5 µg tetracycline ml–1). The liquid culture was incubated at 28 °C for 24 h with agitation (150 r.p.m.) and then poured into a 200 ml flask containing 80 ml SDB and incubated for an additional 48 h under the same regimen. Subsequently, the liquid culture was distributed in 10 ml aliquots onto 90 mm Petri dishes containing only 2 % agar. Excess liquid was removed using sterile filter paper and mycelial mats were used for inoculation when the fungus had sporulated uniformly.

Plutella xylostella were maintained on Brassica napus plants in a constant environment room under standard temperature and photoperiodic conditions (23~25 °C; L12 : D12). Pupae were collected and stored at 4 °C until sufficient quantities were obtained and then transferred into an oviposition cage for eclosion. Adults were allowed to lay ova on a perforated plastic bag containing a fresh cabbage leaf, and the bag was subsequently used to encase a plant during incubation. Second to third instar Plutella xylostella larvae were inoculated using the conidial shower method (Xu & Feng, 2000Down). Petri plates containing larvae on a cabbage leaf were placed beneath an inverted sporulating culture to allow the conidia to discharge directly onto the larval epidermis. Leaves and their resident larvae were transferred to Petri dishes covered with a porous membrane and were incubated at 20 °C at high relative humidity (90 %). Live and dead larvae were recorded daily. Dead larvae were removed and examined individually for verification of infection after overnight incubation.

All data were corrected for control mortality and analysed by using time–concentration–mortality modelling (Nowierski et al., 1996Down; Xu & Feng, 2000Down). The model includes I concentrations and J times of observation. The cumulative mortality probability (pij) related to the concentration Ci (i=1, 2, ..., I) at time tj (j=1, 2, ..., J) can be expressed as:

Down

Formula 001
where β is the slope to describe the concentration effect, and {tau}j is the parameter(s) for the time effect of Ci during the period from start to the jth observation, (t1, t2, ..., tj-1, tj, tj+1,..., tJ). To guarantee that the observed mortality probability is independent of time, the true mortality induced with Ci at the interval [tj–1, tj] must be considered. This true mortality, qij, is called the conditional mortality probability (Robertson & Preisler, 1992Down) and is given as follows:Down

Formula 002
where β is equal to that in equation 1 and {gamma}j describes the conditional time effect of Ci at the interval [tj–1, tj]. The independence between time intervals of observation allowed the fitting of data to equation 2 by approaching the binomial response variable to the maximum-likelihood equation (Nowierski et al., 1996Down; Xu & Feng, 2000Down), yielding {gamma}j and β. The values of {tau}j were then calculated using the formula (Robertson & Preisler, 1992Down) below:Down

Formula 003
The procedures, including modelling, estimation of time- and concentration-effect parameters for both conditional and cumulative models, test for goodness of fit, and estimation of virulence indices (LC20 and LT20) using the parameters, were conducted using the DPS data processing system software (Tang & Feng, 2002Down).

cDNA library construction.
Fungal mycelia were collected after 3 days of growth on SEMA at room temperature and ground in liquid nitrogen. Total RNA was isolated by dispensing 200 mg (wet wt) ground mycelia into a 1.5 ml micro-centrifuge tube containing 600 µl extraction buffer [0.1 M NaCl, 2 % SDS, 50 mM Tris/HCl (pH 9.0), 10 mM EDTA (pH 8.0)] and 600 µl phenol/chloroform/isoamyl alcohol (25 : 24 : 1, by vol.). The solution was mixed using a vortex for 30 s, transferred to a new tube and extracted once with 250 µl chloroform/isoamyl alcohol, and the RNA was precipitated with 0.1 vols 3 M sodium acetate (pH 5.2) and 2.5 vols 95 % ethanol. The RNA pellet was washed with 70 % ethanol, dried for 5 min and resuspended in 50 µl RNase-free ddH2O. Poly(A)+ mRNA was purified on an oligo(dT)-cellulose column using an mRNA purification kit (Amersham Pharmacia Biotech). The cDNA library was constructed using the ZAP cDNA synthesis kit (Stratagene) and amplified in Escherichia coli XL1-Blue. In vivo mass excision of pBluescript clones was performed according to the manufacturer's instructions (Stratagene). Colonies were randomly picked and put into wells of a 384-well microtitre plate with 80 µl freezing broth [1 % tryptone, 0.5 % yeast extract, 1 % NaCl, 4.5 % glycerol, 36 mM K2HPO4 . 3H2O, 13.2 mM KH2PO4, 1.7 mM sodium citrate, 0.41 mM MgSO4 . 7H2O, 6.8 mM (NH4)2SO4] with 100 µg carbenicillin ml–1. The plate was incubated for 16 h at 37 °C and stored at –80 °C.

DNA sequencing and EST analysis.
Bacteria in the 384-well plate were transferred using a pin tool to a 4x96-well flat bottom block with each well containing 1.25 ml LB medium plus 100 µg carbenicillin ml–1. Plasmid DNA was isolated using the DirectPrep 96 Miniprep kit (Qiagen). cDNA inserts were sequenced from both the 5' and 3' termini using T3 (forward) and T7 (reverse) primers, respectively, on an ABI Prism 3700 DNA sequencer using Big Dye Terminator chemistry (Applied Biosystems). The resultant sequences were processed for sequence quality by selecting those having a Phred score (Ewing et al., 1998Down) of more than 20 and the vector sequences were removed. Overlapping forward and reverse sequences were aligned using Sequencher 4.1 (GeneCodes) and combined to form a single EST contig. The sequences were translated and compared to international databases (non-redundant GenBank CDS translations, PDB, Swiss-Prot, PIR, PRF) using the gapped BLASTX algorithm (Altschul et al., 1997Down).


    RESULTS AND DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Virulence variation of Z. radicans original strain adapted to Plutella xylostella
All strains were pathogenic toward Plutella xylostella larvae when applied using the spore shower method; however, the new strains adapted to Plutella xylostella were more virulent than the original strain derived from Pieris brassicae. The number of larvae inoculated and the cumulative mortality after exposure to varying spore doses are provided in Table 1Down. The deaths attributed to infection by the original and adapted strains R1, R3 and R5 occurred mostly on days 3–5, days 2–4, days 2–4 and days 3–5, respectively. All cadavers displayed typical symptoms for Z. radicans infection upon microscopic examination after overnight incubation in a humidified chamber. The cumulative mortality at the end of the observation period (6 days) ranged from 6.3 to 15.2 % for the original strain, 7.5–33.3 % for R1, 5.1–28.9 % for R3 and 5.5–17.6 % for R5 as the spore concentration increased. The cumulative mortality of the control larvae was 2.0–3.9 % on day 6 (Table 1Down).


View this table:
[in this window]
[in a new window]

 
Table 1. Virulence of a Z. radicans isolate derived from Pieris brassicae and isolates successively passaged through Plutella xylostella against Plutella xylostella larvae

 
The time–concentration–mortality data were adjusted using the background mortality records and then fitted to the conditional mortality probability model (equation 2), yielding parameters for conditional concentration effects (β) of the four strains assayed (Table 1Up). The t statistics for all parameters estimated were significant (P<0.01). The Hosmer–Lemeshow test for the heterogeneity of the goodness of the fit (Nowierski et al., 1996Down) was insignificant for all strains (ARSEF 1342, C=4.45, P=0.81; R1-1342, C=10.11, P=0.26; R3-1342, C=2.93, P=0.94; R5-1342, C=7.02, P=0.53), indicating that the data for the Z. radicans isolates against Plutella xylostella fit the model well. The time–concentration–mortality relationships of four strains against Plutella xylostella based on the fitting of the cumulative and conditional mortality probabilities from equations 1 and 2 are shown in Fig. 1Down.


Figure 1
View larger version (82K):
[in this window]
[in a new window]

 
Fig. 1. Cumulative mortality probabilities based on the time–concentration–mortality of a Z. radicans isolate (ARSEF 1342) from Pieris brassicae (R0) or isolates successively passaged through Plutella xylostella (R1, R3 and R5) against Plutella xylostella larvae.

 
The slope, β, for concentration effect of the original strain (0.22±0.07) was not significantly different from that for R1 (0.37±0.08; t=1.89, P>0.01) whereas that of R3 and R5 (1.21±.0.06 and 1.17±0.10, respectively) was significantly different from that for original strain (t=13.17 and 11.19, respectively, P<0.01), indicating that virulence increases with more passages through Plutella xylostella. In addition, according to the value of parameter {gamma}j estimated by modelling, the mortality peak is the third day for the original strain and the second day for all strains adapted to Plutella xylostella, indicating that the infection time was reduced for strains adapted to Plutella xylostella.

Based on time–concentration–mortality modelling, the LC20 (the lethal concentration required to cause 20 % mortality) on days 1–6 after exposure was estimated as 14.7, 14.5, 9.0, 7.1, 6.0 and 5.5 spores mm–2, respectively, for the original strain; 9.6, 5.0, 4.2, 3.6, 3.1 and 2.9 spores mm–2 for R1; 4.6, 2.9, 2.8, 2.5, 2.4 and 2.2 spores mm–2 for R3; and 5.2, 3.7, 3.2, 2.8, 2.6 and 2.6 spores mm–2 for R5. Thus the adapted strains tended to exhibit higher infectivity to the new host than the original isolate as the number of host passages increased.

cDNA library analysis and characterization of the EST data
Two cDNA libraries were constructed for Z. radicans EST analysis. One library (ZRA) was constructed using mycelia from a strain isolated from the original host (Pieris brassicae) and a second library (ZRB) was constructed using mycelia of an isogenic strain that had been passaged and reisolated from an alternate host (Plutella xylostella). The libraries from the strains isolated from Pieris brassicae and Plutella xylostella had primary titres of 5.18x105 and 7.39x105 plaques, respectively, with mean insert sizes of 1.35 kbp and ranging from 0.65 to 3.5 kbp. A combined total of 1536 randomly selected cDNA clones from both Z. radicans cDNA libraries were sequenced, yielding 1839 forward and reverse ESTs of satisfactory quality. For the Pieris brassicae isolate, a total of 844 ESTs of satisfactory quality were generated, representing 455 unique sequences consisting of 220 singletons and 235 contigs with an abundance of 2–40 times in the dataset. For the strain isolated from Plutella xylostella, a total of 995 ESTs with satisfactory quality were generated, representing 533 unique sequences consisting of 249 singletons and 284 contigs with an abundance of 2–29 times. In total, 899 unique Z. radicans sequences consisting of 433 singletons and 466 contigs were obtained of which 837 (284 contigs and 553 singletons) were similar to other sequences (E≤e–03 and 95.2 % with E≤e–05) in the NCBI ‘nr’ and Saccharomyces Genome databases. For this group, 69.65 % of the sequences encoded proteins with a known function with the most similar orthologues being from fungi (54.6 %) and animal (31.9 %) and to a lesser extent from plants (8.1 %) and bacteria (5.4 %).

Functional classification of expressed Z. radicans genes
We compared our EST data to several databases using BLASTX with a threshold E-value of ≤10–3 and grouped the annotated genes into 14 categories based on the Saccharomyces cerevisiae Functional Catalogue developed by the Munich Information Center for Protein Sequences (MIPS) (http://mips.gsf.de/projects/funcat) and the Gene Ontology consortium (http://www.geneontology.org/GO.doc.shtml). These results and corresponding GenBank accession numbers for each unique sequence are available in supplementary Table S1 (available with the online version of this paper). As shown in Fig. 2Down, most of the proteins were involved in protein synthesis (22.44 %) and metabolism (15.88 %), with the others involved in cell cycle and DNA synthesis (5.27 %), protein fate (5.53 %), transport facilitation (4.89 %), cell defence (4.18 %), transcription and RNA synthesis (3.43 %), energy (2.96 %), cellular transport (2.64 %), signal transduction (1.41 %), control of cellular organization (0.71 %) and cell-wall degradation (0.13 %). Approximately 30.35 % of transcripts were similar to hypothetical proteins or proteins for which a function has yet to be determined, indicating that a great deal of functional analysis of genes from entomopathogenic fungi remains to be done.


Figure 2
View larger version (21K):
[in this window]
[in a new window]

 
Fig. 2. Functional classification of ESTs from a Z. radicans isolate from Pieris brassicae (ZRA; black bars) and an isogenic isolate adapted to Plutella xylostella (ZRB; grey bars).

 
Highly expressed Z. radicans genes
Since both cDNA libraries were constructed using mycelia grown under the same conditions, the relative expression level of ESTs could be used to identify genes important for Z. radicans growth, and possibly those that are differentially expressed in strains isolated from different hosts (Table 2Down). Genes encoding ribosomal proteins and translational machinery, such as elongation factors, were abundantly expressed. The longer form of elongation factor 1{alpha} (EX148848) was encoded by 69 cDNAs and represented by 58 % of the ESTs in the original Pieris brassicae isolate and 42 % in the strain passaged through Plutella xylostella. Genes encoding other ribosomal proteins, including ribosomal protein L6 (EX14295), L7A (EX149294), L8 (EX148837), L9 (EX148813), L11 (EX149252), L13 (EX149271), L15 (EX149263), L19 (EX149259), L21 (EX149285), S4 (EX148831), S6 (EX148816) and S-phase-specific ribosomal protein (EX149289), were also abundant in both libraries. High levels of expression of genes associated with protein synthesis indicated that the Z. radicans isolates require new protein synthesis for rapid growth, which is similar to other entomophthoralean fungi. Conidiobolus coronatus has been reported to express a much higher proportion of genes encoding ribosomal proteins and enzymes involved in intermediate metabolism that facilitate rapid growth compared to Metarhizium anisopliae (Freimoser et al., 2003aDown).


View this table:
[in this window]
[in a new window]

 
Table 2. Most prevalent Z. radicans EST clusters

 
Genes involved in energy metabolism, amino acid metabolism, carbohydrate metabolism, lipid/fatty acid and isoprenoid metabolism were also highly expressed. cDNAs encoding proteins associated with energy metabolism, including ATP synthase β chain (EX148879), ATP synthase H+ transporter (EX149352) and ADP/ATP translocase (EX148881), were found 13–16 times in the libraries. cDNAs encoding enzymes associated with amino acid metabolism, such as adenosylhomocysteinase (EX148709) and aconitate hydratase (EX148707), carbohydrate metabolism, such as Gapd-prov protein (EX148721), and lipid/fatty acid and isoprenoid metabolism, such as phosphoglycerate kinase (EX148742), were more abundant in the library from the original Pieris brassicae isolate than in the strain passaged through Plutella xylostella. Although the EST dataset is small, this observation may indicate that the original isolate may derive nutrients from carbohydrate, amino acid and lipid/fatty acid sources, and may correspond to our observation that its growth is more rapid on SEMA medium than the isolate from Plutella xylostella.

cDNAs encoding proteins involved in the stress response, such as those of the 70 kDa heat-shock protein (Hsp70) family, including the PrBiP precursor (EX148912), Hsp70 protein 2 (EX149384) and ribosome-associated Hsp70-like protein 3 (EX149538), were highly represented in both libraries. cDNA encoding Hsp70-like protein 3 was found in the library from the Plutella xylostella isolate. These proteins are presumably involved in protecting the fungus from various environmental stresses (Da Silva et al., 1999Down).

Other transcripts, such as those encoding histone H4 (EX148794), cyclophilin A (EX149321 and EX149323), polyubiquitin (EX148866) and trypsin-like serine protease (EX149338), were also found to be highly represented in both libraries. cDNAs encoding histone H4, which is associated with DNA synthesis and processing, were more prevalent in the library from the Pieris brassicae isolate than from the strain passaged through Plutella xylostella.

Proteinase genes expressed in Z. radicans isolates
Several aspects of the fungus–insect interaction, including spore attachment and germination on the cuticle, secretion of cuticle-degrading enzymes and penetration of the integument, recognition by the insect defence system and enablement of such, come to bear on the success or failure of the pathogen to infect the host. Many virulence-associated genes have been identified from human fungal pathogens, such as Candida albicans (Goldman et al., 2003Down; Rodier et al., 1999Down), and plant-pathogenic fungi, such as Magnaporthe grisea (Viaud et al., 2002Down) and Ustilago maydis (Lee & Kronstad, 2002Down). Only a few virulence genes have been identified from entomopathogenic fungi, such as Metarhizium anisopliae (Freimoser et al., 2003aDown), Conidiobolus coronatus (Freimoser et al., 2003bDown) and Beauveria bassiana (Cho et al., 2006aDown, bDown).

Among the most widely studied pathogenicity and virulence factors for entomopathogenic fungi are proteinases as they are involved in degradation of insect cuticle and impairment of host defence responses (Joshi & St Leger, 1999Down; Screen & St Leger, 2000Down; Freimoser et al., 2003aDown, bDown; Bagga et al., 2004Down; Cho et al., 2006aDown, bDown; Xu et al., 2006Down). Both isolates expressed genes encoding secreted proteinases (containing a signal peptide), including several aspartic, serine and metalloproteases (Table 3Down); however, neither strain expressed genes encoding cysteine proteinases. Aspartic proteinases function in acidic environments and have been found in plant-pathogenic fungi, such as Botrytis cinerea, which release oxalic acid (ten Have et al., 2004Down). We identified two types of aspartic proteinases (EX149341 and EX149324), both of which were found in the library from the isolate from Plutella xylostella. Serine and metalloproteinases are common among the entomopathogenic fungi and are considered to be important for cuticle penetration (Xu et al., 2006Down; St Leger et al., 1994Down; Freimoser et al., 2003aDown, bDown). The library from the Plutella xylostella isolate contained cDNAs encoding a trypsin (EX149336 and EX149337), a trypsin-like serine protease (EX149338), an extracellular trypsin protease (EX149340) and a serine-type protease (EX149333), whereas the library from the Pieris brassicae isolate contained cDNAs encoding a cuticle-degrading proteinase similar to CDEP-1 (EX148865), an Astryp1-like protease (EX148859) and a trypsin-like serine protease (EX149338). Astryp1 is an acid proteinase and may function in the same environment as aspartic proteinases. Both isolates expressed genes encoding trypsins, but not subtilisins and chymotrypsins, which are common in other entomopathogenic fungi such as Metarhizium anisopliae (Freimoser et al., 2003aDown; Bagga et al., 2004Down) and Conidiobolus coronatus (Freimoser et al., 2003bDown). In this study, we found that the Plutella xylostella isolate also expressed genes encoding metalloproteinases (EX149342, EX149343, EX149344, EX149346 and EX149322). Previously, the original Z. radicans isolate from Pieris brassicae was shown to produce a metalloproteinase when grown in submerged medium for 30 days (Xu et al., 2006Down). Many researchers have reported that proteinases, including trypsin-like serine, metallo- and aspartic proteinases, are involved in virulence (Xu et al., 1998Down). St Leger et al. (1994)Down found that the subtilisin-like proteinase (Pr1), the trypsin-like serine (Pr2) and a thermolysin-like metalloproteinase in Metarhizium anisopliae are vital for digestion of cockroach cuticle. In addition, aspartic and metalloproteinases seem to play an important role in invasion by Candida albicans and could be considered as potential pathogenicity factors (Rodier et al., 1999Down).


View this table:
[in this window]
[in a new window]

 
Table 3. cDNAs encoding proteinases identified in the original and newly adapted isolates of Z. radicans

 
Other candidate genes associated with fungal virulence and pathogenesis
Attachment of the spore and germination on the cuticle are the first steps for infection by Z. radicans. These processes require nutrient acquisition for germ-tube extension and secretion of substances for anchoring the appressorium to the insect cuticle. Both isolates expressed isocitrate lyase (EX149145), suggesting the glyoxylate cycle is operating in both fungi (Table 4Down). Isocitrate lyase has been associated with pathogenicity and virulence in several pathogenic fungi as it is a key enzyme in the glyoxylate bypass pathway, which circumvents two decarboxylation steps of the citric acid cycle, allowing for more timely assimilation of carbon substrates. The mammalian pathogen Candida albicans also requires the glyoxylate cycle for assimilation of lipids (Lorenz & Fink, 2001Down). When phagocytosed by macrophages and neutrophils in mice, Candida albicans upregulates principal enzymes of the glyoxylate cycle, isocitrate lyase and malate synthase, and mutants lacking isocitrate lyase are markedly less virulent than the wild-type. The glyoxylate cycle is also required for virulence of plant-pathogenic fungi such as Magnaporthe grisea (Wang et al., 2003Down) and Leptosphaeria maculans (Idnurm & Howlett, 2002Down). Wang et al. (2003)Down found that Magnaporthe grisea employs isocitrate lyase to catabolize plant lipids for turgor generation and appressorium formation prior to infection. Although L. maculans enters the host through stomata without the aid of an appressorium, isocitrate lyase is also essential for successful colonization on Brassica napus (Idnurm & Howlett, 2002Down). However, isocitrate lyase is not essential for disease development in an isolate of Saccharomyces cerevisiae (Goldstein & McCusker, 2001Down) or in the opportunistic animal pathogen Cryptococcus neoformans (Rude et al., 2002Down). Whether isocitrate lyase is involved in disease development in entomopathogenic fungi remains to be resolved.


View this table:
[in this window]
[in a new window]

 
Table 4. Z. radicans cDNAs representing proteins associated with pathogenicity in original and newly adapted isolates

 
Glucosidase has also been implicated as having a role in the infection process and we found one EST (EX149127) in the database from the Plutella xylostella isolate of Z. radicans. Glucosidases are involved in the secretion of adhesive mucilage at the end of the germ tube; this mucilage is crucial for pathogenicity in Metarhizium anisopliae (Freimoser et al., 2003bDown).

After fungi breach the physical barrier provided by the cuticle, they encounter a battery of insect cellular and humoral defences. The current view is that components on the surface of the pathogen, such as 1,3-β-glucan from fungi or peptidoglycan from bacteria, are recognized by specific receptors on haemocytes which in turn trigger a variety of defence responses, including prophenoloxidase activation, phagocytosis, nodule formation and encapsulation. Thus subtle changes in the cell wall could modify adhesion and recognition by the defence system (Navarro-Garcia et al., 2001Down; Goldman et al., 2003Down). Both isolates expressed genes encoding proteins involved in the synthesis of 1,3-β-glucan (EX148729, EX149121, EX149132 and EX149372). The Plutella xylostella isolate expressed a gene encoding thiol peroxidase (EX149380), which is considered to be an antioxidant protein involved in fungal pathogenicity (Wu et al., 1997Down) and may protect the pathogen from host oxidative defences (Iwanaga & Kawabata, 1998Down). In addition, both strains expressed a gene encoding superoxide dismutase (SOD) (EX149194 and EX148763), an antioxidant enzyme which catalyses the dismutation of superoxide radical anions to di-oxygen and hydrogen peroxide. Cu/Zn SOD has been shown to play a role in protecting cells against oxygen toxicity by decomposing superoxide radical anions, and has been proposed to be a virulence determinant in some pathogenic organisms. For example, Candida albicans strains lacking Cu/Zn SOD were more susceptible to macrophage attack and showed attenuated virulence in mice (Hwang et al., 2002Down).

Both Z. radicans strains expressed genes encoding cyclophilin A (EX149321 and EX149323), a member of the cyclophilin family of enzymes that have cistrans isomerase activity. They have been implicated in cellular processes, including stress response, signalling, cell-cycle control, RNA splicing and facilitating viral replication (Howard et al., 2003Down). Cyclophilin A is present in bacteria, fungi, plants and animals, and is required for virulence of both plant-pathogenic fungi, such as Magnaporthe grisea (Viaud et al., 2002Down), and human-pathogenic fungi, such as Cryptococcus neoformans (Wang et al., 2001Down).

The database from the Plutella xylostella isolate also contained ESTs representing GTP-binding proteins (EX149327, EX148886 and EX149365), which are required for pathogenicity, budding growth and mating in Ustilago maydis (Lee & Kronstad, 2002Down). The Pieris brassicae isolate expressed a chitin synthase (EX148735) that was similar to that associated with virulence of Wangiella dermatitidis against mice (Wang et al., 2001Down). While entomopathogenic fungi, such as Metarhizium anisopliae and Conidiobolus coronatus, secrete toxins (Freimoser et al., 2003aDown, bDown), we did not identify genes encoding enzymes involved in their synthesis within the limited Z. radicans EST dataset. Cuticle degradation and evasion of insect defences appear to be important pathogenicity and virulence factors for Z. radicans.

Conclusions
Although the EST datasets presented here represent a small portion of the expressed Z. radicans genes, the analysis has provided insight into the mechanisms required for host infection. Around 80 % of the unique ESTs were found only in libraries from one isolate and encoded proteins involved in all functional categories, indicating that much more information can still be gathered from this approach. The study found that the highly expressed genes encoded proteins involved in basic metabolism, such as ATP-generating enzymes, histones and some ribosomal proteins. These have also been found to be highly expressed in Beauveria bassiana (Cho et al., 2006aDown, bDown). We found that genes encoding potential pathogenicity and virulence factors, such as a trypsin-like serine protease and superoxide dismutase, were also highly expressed. In addition, genes encoding several ribosomal proteins such as L11, L13, L15, L19 and L21 and S-phase-specific ribosomal protein, genes involved in the stress response (ribosome-associated Hsp70-like protein) and metabolism (phosphoglycerate kinase and aconitate hydratase) were expressed in the two isolates and may be required for adaptation to different environmental stress, nutrient availability or the defence response associated with the different hosts. The knowledge gained through this study has contributed to the understanding of entomophthoralean fungal pathogenesis and will assist in the development of highly virulent and host-specific isolates.


    ACKNOWLEDGEMENTS
 
This work was supported with funding from the Chinese Scholarship Council Foundation and the National Natural Science Foundation of China (grant No. 30470063).

Edited by: S. D. Harris


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Altschul, S. F., Madden, T. L., Schäffer, A. A., Zhang, J., Zhang, Z., Miller, W. & Lipman, D. J. (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25, 3389–3402.[Abstract/Free Full Text]

Bagga, S., Hu, G., Screen, S. E. & St Leger, R. J. (2004). Reconstructing the diversification of subtilisins in the pathogenic fungus Metarhizium anisopliae. Gene 324, 159–169.[CrossRef][Medline]

Brown, D. W., Cheung, F., Proctor, R. H., Butchko, R. A., Zheng, L., Lee, Y., Utterback, T., Smith, S., Feldblyum, T. & other authors (2005). Comparative analysis of 87,000 expressed sequence tags from the fumonisin-producing fungus Fusarium verticillioides. Fungal Genet Biol 42, 848–861.[CrossRef][Medline]

Brownbridge, M., Costa, S. & Jaronski, S. T. (2001). Effects of in vitro passage of Beauveria bassiana on virulence to Bemisia argentifolii. J Invertebr Pathol 77, 280–283.[CrossRef][Medline]

Cho, E. M., Liu, L., Farmerie, W. & Keyhani, N. O. (2006a). EST analysis of cDNA libraries from the entomopathogenic fungus Beauveria (Cordyceps) bassiana. I. Evidence for stage-specific gene expression in aerial conidia, in vitro blastospores and submerged conidia. Microbiology 152, 2843–2854.[Abstract/Free Full Text]

Cho, E. M., Boucias, D. & Keyhani, N. O. (2006b). EST analysis of cDNA libraries from the entomopathogenic fungus Beauveria (Cordyceps) bassiana. II. Fungal cells sporulating on chitin and producing oosporein. Microbiology 152, 2855–2864.[Abstract/Free Full Text]

Da Silva, S. P., Borges-Walmsley, M. I., Pereira, I. S., Soares, C. M. A., Walmsley, A. R. & Felipe, M. S. S. (1999). Differential expression of an hsp70 gene during transition from the mycelial to the infective yeast form of the human pathogenic fungus Paracoccidioides brasiliensis. Mol Microbiol 31, 1039–1050.[CrossRef][Medline]

Ebbole, D. J., Jin, Y., Thon, M., Pan, H., Bhattarai, E., Thomas, T. & Dean, R. (2004). Gene discovery and gene expression in the rice BLAST fungus, Magnaporthe grisea: analysis of expressed sequence tags. Mol Plant Microbe Interact 17, 1337–1347.[Medline]

Ewing, B., Hillier, L., Wendl, M. C. & Green, P. (1998). Base-calling of automated sequencer traces using Phred. I. Accuracy assessment. Genome Res 8, 175–185.[Abstract/Free Full Text]

Freimoser, F. M., Screen, S., Bagga, S., Hu, G. & St Leger, R. J. (2003a). Expressed sequence tag (EST) analysis of two subspecies of Metarhizium anisopliae reveals a plethora of secreted proteins with potential activity in insect hosts. Microbiology 149, 239–247.[Abstract/Free Full Text]

Freimoser, F. M., Screen, S., Hu, G. & St Leger, R. J. (2003b). EST analysis of genes expressed by the zygomycete pathogen Conidiobolus coronatus during growth on insect cuticle. Microbiology 149, 1893–1900.[Abstract/Free Full Text]

Furlong, M. J. & Pell, J. K. (2001). Horizontal transmission of entomopathogenic fungi by the diamondback moth. Biol Control 22, 288–299.[CrossRef]

Glare, T. R. (1988). Variation in the entomopathogenic fungal genus Zoophthora, with special reference to Z. radicans. PhD dissertation, Australian National University.

Glare, T. R. & Milner, R. J. (1991). Ecology of entomopathogenic fungi. In Handbook of Applied Mycology, Vol. 2: Humans, Animals, and Insects, pp. 547–612. Edited by D. K. Arora, L. Ajello & K. G. Mukerji. New York: Marcel Dekker.

Goldman, G. H., dos Reis Marques, E., Duarte Ribeiro, D. C., de Souza Bernardes, L. A., Quiapin, A. C., Vitorelli, P. M., Savoldi, M., Semighini, C. P., de Oliveira, R. C. & other authors (2003). Expressed sequence tag analysis of the human pathogen Paracoccidioides brasiliensis yeast phase: identification of putative homologues of Candida albicans virulence and pathogenicity genes. Eukaryot Cell 2, 34–48.[Abstract/Free Full Text]

Goldstein, A. L. & McCusker, J. H. (2001). Development of Saccharomyces cerevisiae as a model pathogen: a system for the genetic identification of gene products required for survival in the mammalian host environment. Genetics 159, 499–513.[Abstract/Free Full Text]

Howard, B. R., Vajdos, F. F., Li, S., Sundquist, W. I. & Hill, C. P. (2003). Structural insights into the catalytic mechanism of cyclophilin A. Nat Struct Biol 10, 475–481.[CrossRef][Medline]

Humber, R. A. (1992). Collection of Entomopathogenic Fungal Cultures: Catalog of Strains. United States Department of Agriculture Publication ARS-100. Ithaca, NY: Agricultural Research Service.

Hwang, C. S., Rhie, G., Oh, J. H., Huh, W. K., Yim, H. S. & Kang, S. O. (2002). Copper- and zinc-containing superoxide dismutase (Cu/ZnSOD) is required for the protection of Candida albicans against oxidative stresses and the expression of its full virulence. Microbiology 148, 3705–3713.[Abstract/Free Full Text]

Idnurm, A. & Howlett, B. J. (2002). Isocitrate lyase is essential for pathogenicity of the fungus Leptosphaeria maculans to canola (Brassica napus). Eukaryot Cell 1, 719–724.[Abstract/Free Full Text]

Iwanaga, S. & Kawabata, S. (1998). Evolution and phylogeny of defense molecules associated with innate immunity in horseshoe crab. Front Biosci 3, D973–D984.[Medline]

Jeong, H. Y., Han, D. M., Jahng, K. Y. & Chae, K. S. (2000). The rpl16a gene for ribosomal protein L16A identified from expressed sequence tags is differentially expressed during sexual development of Aspergillus nidulans. Fungal Genet Biol 31, 69–78.[CrossRef][Medline]

Joshi, L. & St Leger, R. J. (1999). Cloning, expression, and substrate specificity of MeCPA, a zinc carboxypeptidase that is secreted into infected tissues by the fungal entomopathogen Metarhizium anisopliae. J Biol Chem 274, 9803–9811.[Abstract/Free Full Text]

Keon, J., Antoniw, J. R., Skinner, J. H. & Hammond-Kosack, K. (2005). Analysis of expressed sequence tags from the wheat leaf blotch pathogen Mycosphaerella graminicola (anamorph Septoria tritici). Fungal Genet Biol 42, 376–389.[CrossRef][Medline]

Lee, N. & Kronstad, J. W. (2002). Ras2 controls morphogenesis, pheromone response, and pathogenicity in the fungal pathogen Ustilago maydis. Eukaryot Cell 1, 954–966.[Abstract/Free Full Text]

Li, J., Xu, J. H. & Feng, M. G. (2004a). Infectivity of a Pieris brassicae-derived Zoophthora radicans isolate and its new host-passage isolates against Plutella xylostella in relation to the phenoloxidase activity in the new host hemolymph after infection. Acta Entomol Sin 47, 567–572.

Li, R., Rimmer, R., Buchwaldt, L., Sharpe, A. G. & Swartz, G. S. (2004b). Interaction of Sclerotinia sclerotiorum with a resistant Brassica napus cultivar: expressed sequence tag analysis identifies genes associated with fungal pathogenesis. Fungal Genet Biol 41, 735–753.[CrossRef][Medline]

Lorenz, M. C. & Fink, G. R. (2001). The glyoxylate cycle is required for fungal virulence. Nature 412, 83–86.[CrossRef][Medline]

Navarro-Garcia, F., Sánchez, M., Nombela, C. & Pla, J. (2001). Virulence genes in the pathogenic yeast Candida albicans. FEMS Microbiol Rev 25, 245–268.[Medline]

Nowierski, R. M., Zeng, Z., Jaronski, S., Delgado, F. & Swearingen, W. (1996). Analysis and modeling of time–dose–mortality of Melanoplus sanguinipes, Locusta migratoria migratorioides, and Schistocerca gregaria (Orthoptera: Acrididae) from Beauveria, Metarhizium, and Paecilomyces isolates from Madagascar. J Invertebr Pathol 67, 236–252.[CrossRef][Medline]

Papierok, B., Valadao, L., Torres, B. & Arnault, M. (1984). Contribution à l'étude de la spécificité parasitaire du champignon entomopathogéne Zoophthora radicans (Zygomycetes, Entomophthorales). Entomophaga 29, 109–119.[CrossRef]

Pell, J. K., Eilenberg, J., Hajek, A. E. & Steinkraus, D. C. (2001). Biology, ecology and pest management potential of Entomophthorales. In Fungi as Biocontrol Agents: Progress, Problems and Potential, pp. 71–153. Edited by T. M. Butt, C. Jackson & N. Magan. Wallingford: CAB International.

Posada-Buitrago, M. L. & Frederick, R. D. (2005). Expressed sequence tag analysis of the soybean rust pathogen Phakopsora pachyrhizi. Fungal Genet Biol 42, 949–962.[Medline]

Quesada-Moraga, E. & Vey, A. (2003). Intra-specific variation in virulence and in vitro production of macromolecular toxins active against locust among Beauveria bassiana strains and effects of in vivo and in vitro passage on these factors. Biocontrol Sci Technol 13, 323–340.[CrossRef]

Robertson, J. L. & Preisler, H. K. (1992). Pesticide Bioassays with Arthropods. Boca Raton, FL: CRC Press.

Rodier, M. H., Moudni, B. E., Kauffmann-Lacroix, C., Daniault, G. & Jacquemin, J. L. (1999). A Candida albicans metallopeptidase degrades constitutive proteins of extracellular matrix. FEMS Microbiol Lett 177, 205–210.[CrossRef][Medline]

Rude, T. H., Toffaletti, D. L., Cox, G. M. & Perfect, J. R. (2002). Relationship of the glyoxylate pathway to the pathogenesis of Cryptococcus neoformans. Infect Immun 70, 5684–5694.[Abstract/Free Full Text]

Screen, S. E. & St Leger, R. J. (2000). Cloning, expression, and substrate specificity of a fungal chymotrypsin. J Biol Chem 275, 6689–6694.[Abstract/Free Full Text]

Shah, P. A. & Pell, J. K. (2003). Entomopathogenic fungi as biological control agents. Appl Microbiol Biotechnol 61, 413–423.[Medline]

Sims, A. H., Robson, G. D., Hoyle, D. C., Oliver, S. G., Turner, G., Prade, R. A., Russell, H. H., Dunn-Coleman, N. S. & Gent, M. E. (2004). Use of expressed sequence tag analysis and cDNA microarrays of the filamentous fungus Aspergillus nidulans. Fungal Genet Biol 41, 199–212.[CrossRef][Medline]

St Leger, R. J., Bidochka, M. J. & Roberts, D. W. (1994). Isoforms of the cuticle-degrading Pr1 proteinase and production of a metalloproteinase by Metarhizium anisopliae. Arch Biochem Biophys 313, 1–7.[CrossRef][Medline]

Tang, Q. Y. & Feng, M. G. (2002). Practical Statistics and DPS Data Processing System. Beijing, China: Science Press.

ten Have, A., Dekkers, E., Kay, J., Phylip, L. H. & Kan, J. A. L. (2004). An aspartic proteinase gene family in the filamentous fungus Botrytis cinerea contains members with novel features. Microbiology 150, 2475–2489.[Abstract/Free Full Text]

Thomas, S. W., Rasmussen, S. W., Glaring, M. A., Rouster, J. A., Christiansen, S. K. & Oliver, R. P. (2001). Gene identification in the obligate fungal pathogen Blumeria graminis by expressed sequence tag analysis. Fungal Genet Biol 33, 195–211.[CrossRef][Medline]

Trail, F., Xu, J. R., Miguel, P. S., Halgren, R. G. & Kistlere, H. C. (2003). Analysis of expressed sequence tags from Gibberella zeae (anamorph Fusarium graminearum). Fungal Genet Biol 38, 187–197.[CrossRef][Medline]

Viaud, M. C., Balhadere, P. V. & Talbot, N. J. (2002). A Magnaporthe grisea cyclophilin acts as a virulence determinant during plant infection. Plant Cell 14, 917–930.[Abstract/Free Full Text]

Wall, D. M., Duffy, P. S., DuPont, C., Prescott, J. F. & Meijer, W. G. (2005). Isocitrate lyase activity is required for virulence of the intracellular pathogen Rhodococcus equi. Infect Immun 73, 6736–6741.[Abstract/Free Full Text]

Wang, Z., Zheng, L., Liu, H. B., Wang, Q. F., Hauser, M., Kauffman, S., Becker, J. M. & Szaniszlo, P. J. (2001). WdChs2p, a class I chitin synthase, together with WdChs3p (Class III) contributes to virulence in Wangiella (Exophiala) dermatitidis. Infect Immun 69, 7517–7526.[Abstract/Free Full Text]

Wang, Z. Y., Thornton, C. R., Kershaw, M. J., Li, D. B. & Talbot, N. J. (2003). The glyoxylate cycle is required for temporal regulation of virulence by the plant pathogenic fungus Magnaporthe grisea. Mol Microbiol 47, 1601–1612.[CrossRef][Medline]

Wraight, S. P., Galaini-Wraight, S., Carruthers, R. I. & Robertsa, D. W. (2003). Zoophthora radicans (Zygomycetes: Entomophthorales) conidia production from naturally infected Empoasca kraemeri and dry-formulated mycelium under laboratory and field conditions. Biol Control 28, 60–77.[CrossRef]

Wu, G. S., Shortt, B. J., Lawrence, E. B., Leon, J., Fitzsimmons, K. C., Levine, E. B., Raskin, I. & Shah, D. W. (1997). Activation of host defense mechanisms by elevated production of H2O2 in transgenic plants. Plant Physiol 115, 427–435.[Abstract]

Xu, J. H. & Feng, M. G. (2000). The time-dose-mortality modeling and virulence indices for the two entomophthoralean species, Pandora delphacis and P. neoaphidis, against the green peach aphid, Myzus persicae. Biol Control 17, 29–34.[CrossRef]

Xu, J. H., Feng, M. G. & Xu, Q. (1998). The role of extracellular proteases of entomopathogenic fungi in host-pathogenic process and their assay techniques. Chinese J Biol Control 14, 123–130.

Xu, J. H., Baldwin, D., Kindrachuk, C. & Hegedus, D. D. (2006). Serine proteases and metalloproteases associated with pathogenesis but not host specificity in the Entomophthoralean fungus, Zoophthora radicans. Can J Microbiol 52, 550–559.[CrossRef][Medline]

Received 9 July 2008; revised 19 September 2008; accepted 22 September 2008.



This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary Table
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Xu, J.
Right arrow Articles by Hegedus, D. D.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Xu, J.
Right arrow Articles by Hegedus, D. D.
Agricola
Right arrow Articles by Xu, J.
Right arrow Articles by Hegedus, D. D.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
INT J SYST EVOL MICROBIOL MICROBIOLOGY J GEN VIROL
J MED MICROBIOL ALL SGM JOURNALS
Copyright © 2009 Society for General Microbiology.