Previous Article | Next Article 
Applied and Environmental Microbiology, September 2007, p. 5587-5597, Vol. 73, No. 17
0099-2240/07/$08.00+0 doi:10.1128/AEM.00721-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
Distribution and Expression of Elicitin Genes in the Interspecific Hybrid Oomycete Phytophthora alni
Renaud Ioos,1,2
Franck Panabières,3
Benoît Industri,3
Axelle Andrieux,1 and
Pascal Frey1*
INRA, Nancy-Université, UMR1136, Interactions Arbres-Microorganismes, IFR 110, F-54280 Champenoux,1
Laboratoire National de la Protection des Végétaux, Unité de Mycologie Agricole et Forestière, IFR 110, Domaine de Pixérécourt, F-54220 Malzéville,2
INRA, UMR1064, Interactions Plantes-Microorganismes et Santé Végétale, F-06930 Sophia Antipolis Cedex, France3
Received 30 March 2007/
Accepted 20 June 2007
 |
ABSTRACT
|
|---|
Phytophthora alni subsp. alni, P. alni subsp. multiformis, and P. alni subsp. uniformis are responsible for alder disease in Europe. Class I and II elicitin gene patterns of P. alni subsp. alni, P. alni subsp. multiformis, P. alni subsp. uniformis, and the phylogenetically close species P. cambivora and P. fragariae were studied through mRNA sequencing and 3' untranslated region (3'UTR)-specific PCRs and sequencing. The occurrence of multiple 3'UTR sequences in association with identical elicitin-encoding sequences in P. alni subsp. alni indicated duplication/recombination events. The mRNA pattern displayed by P. alni subsp. alni demonstrated that elicitin genes from all the parental genomes are actually expressed in this allopolyploid taxon. The complementary elicitin patterns resolved confirmed the possible involvement of P. alni subsp. multiformis and P. alni subsp. uniformis in the genesis of the hybrid species P. alni subsp. alni. The occurrence of multiple and common elicitin gene sequences throughout P. cambivora, P. fragariae, and P. alni sensu lato, not observed in other Phytophthora species, suggests that duplication of these genes occurred before the radiation of these species.
 |
INTRODUCTION
|
|---|
The oomycete Phytophthora alni (stramenopile lineage) is a recently described highly aggressive pathogen specific to alder trees (Alnus spp.) that is spreading all over Europe, especially along rivers (13). P. alni sensu lato comprises three related taxa: P. alni subsp. alni, P. alni subsp. uniformis, and P. alni subsp. multiformis (5). P. alni subsp. multiformis and P. alni subsp. uniformis are scarce in comparison to P. alni subsp. alni (2, 15). Up to now, all three sibling taxa have been exclusively isolated from alder trees, but P. alni subsp. multiformis and P. alni subsp. uniformis are reported to be significantly less aggressive to Alnus than P. alni subsp. alni (4, 30). The three taxa are phylogenetically close to P. cambivora and P. fragariae, two species that, it was previously suggested, may be P. alni subsp. alni's progenitors (3). However, it was recently inferred from nuclear and mitochondrial gene genealogies (15) and microsatellite patterns (16) that only P. alni subsp. alni is a genuine hybrid taxon, originating from hybridization between P. alni subsp. uniformis and P. alni subsp. multiformis. The status of P. alni subsp. multiformis is still questionable, while P. alni subsp. uniformis's genetic features do not fit with a hybrid origin (15, 16). Up to now, reports of other natural hybrids within the genus Phytophthora are scarce and remain confined to limited geographical areas (19, 20), whereas the allopolyploid taxon P. alni subsp. alni is currently thriving throughout Europe. Since this hybrid taxon is significantly more aggressive than its putative parents, P. alni subsp. multiformis and P. alni subsp. uniformis (4), heterosis or genetic rearrangement may be put forward as explanations for this kind of ecological advantage. Except for the observation of polymorphism in the internal transcribed spacer sequence (3), the extent of genetic rearrangement in the allopolyploid taxon P. alni subsp. alni remains unknown, and the expression of distinct genomes within this hybrid taxon has not been investigated yet.
In order to address these issues, multigenic families, such as elicitin genes, should be of great interest since they are particularly prone to duplication and recombination. Elicitin proteins are restricted to the oomycete genus Phytophthora and a few Pythium species (18, 25). They comprise a large family of polypeptides whose intrinsic function remains largely unknown and that can be divided into at least eight classes (18, 23). Elicitin genes are highly transcribed during vegetative growth, as deduced from their representation in a collection of expressed sequence tags from the broad-host-range pathogen Phytophthora parasitica (23). In particular, group I elicitins (ELI-1), characterized by the typical 98-amino-acid elicitin domain, were reported to be the most abundant secreted proteins in culture filtrates (18). Although elicitin genes are not appropriate for phylogenetic studies because of their multigenic family feature, they can be used as a tool for identification purposes since the amino acid sequence of a given elicitin may provide a signature at the species level (25). In addition, the 3' untranslated regions (3'UTRs) of class I elicitin genes are strictly conserved within an individual species but diverge between species to such an extent that sequence alignment is almost impossible (1, 6, 11; Panabières et al., unpublished results).
The aims of the present work were, by taking advantage of the high expression of elicitin genes in the Phytophthora genome, (i) to test if the multiple genomes present in the allopolyploid hybrid P. alni subsp. alni are coexpressed or not and (ii) to assess the occurrence of genetic rearrangements. We studied the occurrence and the distribution of members of class I and II elicitin genes and their expression among the different P. alni taxa. This study included the phylogenetically close species P. cambivora, P. fragariae var. fragariae, and P. fragariae var. rubi, previously suggested as P. alni subsp. alni's progenitors (3), for comparison of their elicitin gene patterns.
 |
MATERIALS AND METHODS
|
|---|
Phytophthora isolates and culture.
French isolates of Phytophthora alni sensu lato and other Phytophthora spp. were collected on naturally infected tissues and isolated on PARBHY medium (28). Foreign isolates of P. alni sensu lato and Phytophthora spp. were obtained from the Centraalbureau voor Schimmelcultures (Utrecht, The Netherlands) or from collaborative researchers (Table 1). All the cultures were kept at 10°C in the dark on V8-agar slants (21) and as small V8-agar blocks flooded with sterile distilled water. Five isolates of P. alni subsp. alni (PAA129, PAA130, PAA143, PAA151, and PAA162), three isolates of P. alni subsp. uniformis (PAU60, PAU84, and PAU89), three isolates of P. alni subsp. multiformis (PAM54, PAM71, and PAM73), two isolates of P. cambivora (PC643 and PCjc17), one isolate of P. fragariae var. fragariae (PFF309), and one isolate of P. fragariae var. rubi (PFR109), selected from different geographical locations, were used for in vitro elicitin production and mRNA studies (Table 1).
Nucleic acid manipulation.
Genomic DNA was extracted from 5-day-old cultures grown in shake culture in liquid V8 juice medium (21) at 20°C using a plant DNA extraction kit (DNeasy plant minikit; QIAGEN, Courtaboeuf, France) by following the manufacturer's instructions with slight modifications. Briefly, ca. 200 mg of fresh mycelium was harvested and mixed in a 2-ml tube with 400 µl of lysis buffer and 4 µl of the RNase A provided with the kit. The mixture was ground for 2 min with two 3-mm tungsten carbide beads at a frequency of 30 Hz, using a mixer mill grinder (Tissuelyser; QIAGEN). The ground solution was subsequently centrifuged for 5 min at 15,000 x g to compact the debris, and the resulting supernatant was treated in accordance with the manufacturer's instructions. Genomic DNA was stored at –20°C until used for PCR tests.
To enhance elicitin mRNA synthesis, oomycete cultures were grown for 3 days in shake culture in liquid elicitin secretion medium (ESM; M. Horta [Algarve University, Portugal], personal communication) at 20°C. The composition of the ESM was 0.05% (wt/vol) KH2PO4, 0.025% (wt/vol) MgSO4·7H2O, 0.1% (wt/vol) asparagine, 1 mg/liter thiamine, 0.05% (wt/vol) yeast extract, and 2% (wt/vol) glucose. The medium was sterilized by filtration through a 0.2-µm membrane. mRNAs were extracted using a QuickPrep micro-mRNA purification kit with oligo(dT) cellulose (Amersham Biosciences, Orsay, France) in accordance with the manufacturer's instructions, resuspended in 40 µl of diethyl pyrocarbonate-treated molecular biology grade water, and stored at –80°C.
Reverse transcription-PCR, cloning, and sequencing of elicitin mRNA.
Polyadenylated RNAs were reverse-transcribed using a SuperScript first-strand synthesis system (Invitrogen, Cergy Pontoise, France) with NotI-oligo(dT) [5'-ATTCGCGGCCGCAGGA(T)16-3'] (25). A PCR was performed on the cDNA template using a combination of the NotI-oligo(dT) primer and degenerate primer 1 (5'-ATGAACTTCCGCGCTCTSYTYGC-3'), initially designed from conserved sequences of class I elicitins located in the peptide signal region (25). This primer was assumed to efficiently anneal to the peptide signal region of every elicitin class I gene unraveled until now. PCRs were carried out in a 20-µl PCR mixture containing 1x polymerase buffer (Sigma-Aldrich, L'Isle d'Abeau, France), 0.9 mM MgCl2, 0.3 µM of each primer, 180 µM deoxynucleoside triphosphates, 0.6 unit of Taq DNA polymerase (Sigma-Aldrich), and 2 µl of template cDNA. Molecular biology grade water was added to 20 µl. For each isolate, PCRs were also performed using genomic DNA as a negative control to test potential amplification of genomic DNA in cDNA amplification. The PCRs were carried out with a GeneAmp 9700 thermal cycler (Applied Biosystems, Foster City, CA). The cycling profile included an initial denaturation step at 95°C for 2 min, followed by 35 cycles of denaturation, annealing, and elongation for, respectively, 20 s at 94°C, 30 s at 60°C, and 1 min at 72°C and a final extension step at 72°C for 7 min. PCR fragments were separated by a 1-h electrophoresis on a 1% agarose gel at 4 V·cm–1. Gels were stained with ethidium bromide, and images were recorded with a charge-coupled device camera and a GELDOC 2000 system (Bio-Rad, Marne-La-Coquette, France).
The PCR products generated with cDNAs were cloned for each of the 15 isolates tested using the pCR4-TOPO TA cloning kit (Invitrogen). Five microliters of the PCR product was transferred into a sterile 1.5-ml microcentrifuge tube, and the amplicons were ligated to a TOPO vector (Invitrogen) and used to transform TOP 10-competent cells (Invitrogen) according to the manufacturer's instructions. Positive clones were selected by PCR amplification of inserts with M13 sequencing primers. Positive clones were selected according to their expected PCR product sizes, corresponding to class I elicitin transcripts. PCR products were purified by centrifugation using a polyethylene glycol 8000 solution as described by Rosenthal et al. (29). Double-stranded DNA sequencing was performed by the dideoxy chain termination method using a T3-T7 sequencing kit on a CEQ 2000 XL DNA sequencer (Beckman, Fullerton, CA). Sequences were edited with Sequencher software (Gene Codes, Ann Arbor, MI) and aligned using ClustalW (33) (Table 2). The cDNA sequences were translated using Fast PCR software, version 3.6.62 (R. Kalendar, FastPCR, PCR primer design, DNA and protein tools, repeats and own database search program [www.biocenter.helsinki.fi/bi/programs/fastpcr.htm]). The isoelectric points (pI) of the deduced proteins were calculated using IEP online software (http://bioweb.pasteur.fr/seqanal/interfaces/iep.html). Multiple amino acid sequence alignments with hierarchical clustering were performed using MultAlin program, version 5.3.3 (8), with Blossum 62 as symbol comparison table. An unrooted phylogram was built using a parsimony analysis and the neighbor-joining method in PAUP* 4.0b10 (32). Bootstrap values were computed on 10,000 replicates.
View this table:
[in this window]
[in a new window]
|
TABLE 2. Characteristics of all the different cDNA and genomic sequences obtained and assignation to a specific 3'UTR group
|
|
3'UTR-specific PCR detection and sequencing.
Based on cDNA sequence alignment, a set of 11 reverse primers was designed to target 3'UTR-specific regions among elicitin-encoding sequences (Table 3). Primers were synthesized by Invitrogen and tested on the 15-isolate panel, then on the entire Phytophthora and Pythium collection listed in Table 1.
The 3'UTR PCR tests were carried out in a 20-µl mixture containing 1x polymerase buffer (Sigma-Aldrich), 1.8 mM MgCl2, 0.45 µM of degenerate primer 1, 0.45 µM of 3'UTR-specific primer, 180 µM deoxynucleoside triphosphates, 0.7 µg·µl–1 bovine serum albumin (Sigma-Aldrich), 0.6 unit of Taq DNA polymerase (Sigma-Aldrich), and 2 µl of template genomic DNA or cDNA. Molecular biology grade water was added to 20 µl. PCR parameters were as indicated above except that the annealing temperature was lowered to 58°C. PCR products were resolved by agarose gel electrophoresis as described above.
In order to verify the specificity of the 3'UTR-specific PCR assays, some of the PCR products were subsequently sequenced using degenerate primer 1 as the sequencing primer.
Nucleotide sequence accession numbers.
mRNA sequences generated in this study were deposited in GenBank under accession numbers DQ012508 to DQ012535. PCR products sequenced with degenerate primer 1 as the sequencing primer were deposited in GenBank under accession numbers EF158401 to EF158426.
 |
RESULTS
|
|---|
Cloning, sequencing, and classification of elicitin-encoding sequences.
A subset of 15 Phytophthora isolates (Table 1) was selected to study the elicitin expression pattern through mRNA sequencing. Using a combination of oligo(dT) primer and a degenerate oligonucleotide designed ahead of the 5' end of the elicitin coding sequence, reverse transcription-PCR of the mRNA extracted for each of these strains showed a strong smeared signal corresponding to a 450- to 550-bp product. For each isolate, the entire amplicon was cloned and 32 individual clones were randomly selected and checked for insert size. Inserts ranged from ca. 480 to 570 bp. Therefore, for each isolate, one to three clones of different sizes were selected for sequencing. A total of 28 different sequences were obtained with clones generated from the 15-isolate panel that were translated along with two unpublished class I elicitin sequences from P. cambivora isolate 143 (PC_cam1) and P. fragariae var. rubi isolate 486 (PFR_fra1) (F. Panabières, unpublished data).
Translation of the 28 cDNA sequences resulted in five distinct amino acid sequences (Fig. 1), which belonged to the class I (98 amino acids [aa]) and class II (99 aa) elicitins, according to the classifications of Ponchet et al. (26) and Qutob et al. (27). The unrooted phylogram based on parsimony analysis separated the sequences into three classes, corresponding to acidic, basic, and highly acidic elicitins (Fig. 2).

View larger version (37K):
[in this window]
[in a new window]
|
FIG. 1. Multiple-amino-acid sequence alignment of the elicitins characterized in this study for P. alni, P. cambivora, and P. fragariae (AE1.1, AE1.2, AE2, BE1, BE2, and HAE1) and well documented acidic ( ), basic (ß), and highly acidic (HAE) elicitins retrieved from the GenBank database: P. cryptogea Z34462, Z34459, Z34460, Z34461 (24); P. infestans AY830090 (17); P. cinnamomi AJ000071 (11); P. sojae AJ007859 (1); and P. parasitica S67432 (6). The sequence corresponding to the signal peptide is underlined. pI are indicated in parentheses, following sequence references.
|
|

View larger version (19K):
[in this window]
[in a new window]
|
FIG. 2. Phylogenetic relationships between elicitins, inferred from the sequence alignment of the mature polypeptides. The unrooted phylogram was constructed using a parsimony analysis and the neighbor-joining method, based on the multiple alignment of elicitin sequences listed in Fig. 1 and using the P. infestans inf6 sequence as an outgroup. Bootstrap values (>50%) from 10,000 replicates are indicated. AE1.1/0.2 and AE2 fall into the class I acidic elicitin group (solid line), whereas BE1 and BE2 are class I basic elicitins (dotted line). HAE1 is closely related to class II highly acidic elicitins (dashed line).
|
|
First, two amino acid sequences deduced from five cDNA clones, including PC_cam1, were identified as basic elicitins. They displayed a predicted pI of 8.22 and several key signatures such as the K13 frequently observed in basic elicitins (26). The two proteins differed only by a single A84V substitution and were therefore designated BE1 and BE2 (basic elicitin).
An alignment of the 98-aa core region corresponding to the mature protein (excluding the 20-aa peptide signal sequence) allowed the identification of two other proteins, differing by a single S40A substitution, among a set of 23 sequences. From alignment with published sequences and from their predicted pI of 4.99, they were clearly identified as acidic elicitins and named AE1 and AE2 (acidic elicitin). Extending the comparison to the entire amino acid sequence further split AE1 into two proteins, which differed by a single S17F mutation in the peptide signal sequence, and these were designated AE1.1 and AE1.2.
Last, a third type of elicitin protein was identified. The deduced peptide comprised 119 aa, including the peptide signal sequence, and displayed a predicted pI of 3.95. From alignment with known elicitin sequences, it was considered to belong to class II of highly acidic elicitins and was accordingly called HAE1 (11, 26).
Elicitin-encoding mRNAs display important 3'UTR variability.
The diversity of elicitin-encoding sequences was also investigated at the nucleotide level on the 3'UTR (Fig. 3).

View larger version (35K):
[in this window]
[in a new window]
|
FIG. 3. Sequence alignment of the regions corresponding to the 3'UTRs of the mRNAs, deduced from cDNA sequencing. 3'UTR-specific groups are indicated on the left and correspond to the clustering of identical or nearly identical sequences. Fourteen different 3'UTR groups were defined among the 28 sequences obtained from our 15-isolate panel. The different 3'UTR groups were designated according to the elicitin class, i.e., a for acidic, ha for highly acidic, and b for basic elicitins. For each 3'UTR, the boxed sequence represents the polymorphic region from which a 3'UTR-specific primer could be designed. The PFF309b 3'UTR sequence was so divergent that it could not be properly aligned. PFF309b was therefore not represented, but a specific reverse primer could be designed (a11-R).
|
|
First, the two cDNA sequences encoding the putative highly acidic elicitin, including the 3'UTR, were identical. This 3'UTR region was consequently designated ha1.
Second, sequences encoding basic elicitins were highly similar, diverging by one to four synonymous substitutions and by a C/T transversion leading to the A84V substitution in the coding region. In addition, the 3'UTRs of the different clones derived from P. alni sensu lato were strictly identical over the entire 168-bp region, with the exception of a single G/C transversion in the transcript obtained from isolate PAA162 (data not shown). This 3'UTR sequence was designated b1. The PC_cam1 sequence was more divergent, displaying 8 synonymous substitutions in the coding region and 14 substitutions, as well as a 5-bp deletion, in the 3'UTR. This 3'UTR sequence was called b2.
Last, for acidic elicitins, 11 clones encoding AE1.1 and 7 cDNA clones encoding AE1.2 were examined. They displayed extensive conservation in the coding region, as only three synonymous mutations further differentiated the two groups. The polymorphism was much higher in the 3'UTRs, which could be classified into 11 3'UTR groups, called a1 to a11. The diversity was the outcome of frequent insertions and deletions as well as single-nucleotide polymorphisms.
Overall, there was no obvious correlation between the classification of 3'UTRs and the protein sequence deduced from the coding region (Table 2). AE2 was associated with a3, a4, a9, and a11 3'UTRs. Similarly, AE1.1 was associated with a1, a2, a6, a7, a8, a9, and a10 3'UTRs, while AE1.2 was associated with a1, a2, a3', and a5 3'UTRs. As a consequence, a given 3'UTR could be associated with two different protein sequences.
Expression and distribution of the different elicitin genes.
The occurrence of the whole set of genes revealed by sequencing of cDNAs was further investigated through 3'UTR-specific PCR tests. Specific primers could be designed for 11 out of the 14 3'UTR groups (Table 3). PCR assays were carried out first with the cDNAs from the 15 isolates of the panel, then with genomic DNA extracted from a large set of 101 P. alni subsp. alni, P. alni subsp. multiformis, P. alni subsp. uniformis, P. cambivora, P. fragariae var. fragariae, and P. fragariae var. rubi isolates and other Phytophthora and Pythium species (Table 1).
Results of the 3'UTR-specific PCR tests with cDNAs for P. alni subsp. alni, P. alni subsp. multiformis, and P. alni subsp. uniformis were identical for all isolates of a given taxon, suggesting that the diversity observed in the elicitin-encoding sequences does not correspond to individual variation, but rather to a high level of elicitin complexity among the three taxa (Fig. 4).

View larger version (26K):
[in this window]
[in a new window]
|
FIG. 4. Distribution of the acidic, basic, and highly acidic elicitin genes resolved in this study for P. alni subsp. alni (Paa), P. alni subsp. multiformis (Pam), P. alni subsp. uniformis (Pau), and the phylogenetically close species P. cambivora and P. fragariae sensu lato, as inferred from sequencing of cDNA and 3'UTR-specific PCR tests conducted on genomic DNA and cDNA libraries. Coding sequences are boxed and were found to be associated with different 3'UTRs. The occurrence of a1, a6, and a9 3'UTRs is deduced only from sequencing, as no 3'UTR-specific primer could be designed from these sequences. With the exception of the HAE1 associated with the ha1 3'UTR in P. cambivora, all the elicitin genes we resolved were shown to be expressed in our 15-isolate panel by 3'UTR-specific PCR with cDNA libraries. a, only for P. cambivora isolate PCjc17; Pff, obtained for P. fragariae var. fragariae; Pfr, obtained for P. fragariae var. rubi.
|
|
In addition, the results of 3'UTR-specific PCR tests conducted on genomic DNA of the 101 Phytophthora and Pythium isolates (Table 1) were in complete agreement with those obtained on the cDNAs from the 15-isolate panel and showed that, except for P. cambivora, the elicitin gene patterns were conserved among the different isolates of each taxon (Fig. 4). Sequencing of a subset of 3'UTR-specific PCR products generated with genomic DNA extracts confirmed their assignation to specific 3'UTRs and resolved three additional 3'UTR sequences, called a2', a3', and a10' and a10'', slightly different from the original cDNA sequences a2, a3, and a10, respectively (Table 2).
The distributions of the various acidic elicitin-related 3'UTRs with cDNA and with genomic DNA were identical, demonstrating that the whole elicitin gene content was actually expressed in the current culture conditions. As a single exception, ha1 was amplified from the genomic DNA of P. cambivora isolates, whereas it was not detected in the cDNA, indicating that, for this species, this class of genes was not expressed during our vegetative-growth conditions.
Overall, P. alni subsp. alni displayed an elicitin gene pattern that combined those observed for P. alni subsp. multiformis and for P. alni subsp. uniformis, with the exception of the a8 3'UTR, which could not be observed in P. alni subsp. alni, due to unexpected cross-annealing of the a8-R PCR primer with the a1 3'UTR sequence, and of the a5 3'UTR, which could not be observed in P. alni subsp. multiformis, due to unexpected cross-annealing with the a3' 3'UTR sequence. These cross-annealing events were suggested by in silico annealing tests and further unambiguously demonstrated by sequencing the PCR products (data not shown). Additionally, the occurrence of the BE2 elicitin gene in P. alni subsp. multiformis or in P. alni subsp. uniformis could not be verified since all the sequenced products generated by b1-specific PCR corresponded to a gene encoding BE1 (Table 2). The BE2-encoding sequence may be either (i) present in the hybrid while not present in the progenitors, thus representing an autapomorphic feature (a derived characteristic unique to a given taxon or monophyletic group) that would have been generated during or after the hybridization event, or (ii) not detectable by sequencing PCR products because of underrepresentation in comparison with BE1-encoding sequences.
In contrast with P. alni sensu lato, analysis of P. cambivora and P. fragariae cDNAs revealed a patchy distribution of the various acidic elicitin-related 3'UTRs. The a8 3'UTR was detected in the cDNAs from the two strains of P. cambivora (the sequence of a6 previously obtained from these strains did not allow the design of specific primers), but only isolate PCjc17 appeared to express additional 3'UTRs, namely, a2' and a7 3'UTRs. Similarly P. fragariae var. fragariae and P. fragariae var. rubi isolates shared a2 and a9 3'UTRs (whose sequences did not allow the design of specific primers), but only the P. fragariae var. rubi isolate expressed an additional gene containing a10, while a11 appeared to be specific to the P. fragariae var. fragariae isolate.
Finally, 3'UTR-specific PCRs did not yield any positive signal with other Phytophthora or Pythium species (Table 1), confirming that the primers designed in this study are specific to P. alni subsp. alni, P. alni subsp. multiformis, P. alni subsp. uniformis, P. cambivora, and P. fragariae and that these five taxa share particular evolutionary relationships.
 |
DISCUSSION
|
|---|
Taken as a whole and inferred both from sequence data and from 3'UTR-specific PCR amplifications, the elicitin gene family is more diverse in P. alni subsp. alni and P. alni subsp. multiformis than in the diploid or nearly diploid taxa P. alni subsp. uniformis, P. cambivora, P. fragariae var. fragariae, and P. fragariae var. rubi. This complexity is consistent with the hybrid status for P. alni subsp. alni but also supports the hypothesis that P. alni subsp. multiformis is probably also an ancient allopolyploid taxon (15). In addition, most of the mRNA patterns observed in the allopolyploid hybrid P. alni subsp. alni are a composite of the mRNA patterns from its putative progenitors, P. alni subsp. multiformis and P. alni subsp. uniformis. In this respect, we showed that, at least for the elicitin genes, the different genomes are currently expressed in the hybrid P. alni subsp. alni. However, additional data would be needed to demonstrate that the entire distinct genome sets are still coexpressed in P. alni subsp. alni. Garcia-Olmedo et al. (12) and Volkov et al. (34) showed that, in the course of evolution of protein-coding genes in allopolyploid plants, the alleles of one parental species are mainly transcriptionally active, whereas the alleles from the other parent are gradually transformed into pseudogenes. By contrast, the results presented here suggest that gene silencing is not yet observed for this highly expressed family of elicitin genes in the allopolyploid hybrid P. alni subsp. alni. These results strengthen the hypothesis that this taxon is of recent origin and still evolving (3), which is in good accordance with its recent emergence as an aggressive alder pathogen (14).
More unexpected, from an evolutionary point of view, is the association of given 3'UTRs, namely, those of a1 and a2, with two distinct coding regions in the hybrid P. alni subsp. alni and in P. alni subsp. multiformis, since a higher selection pressure is assumed to have been exerted on the coding regions. In this respect, it may be hypothesized that recent (P. alni subsp. alni) or likely more ancient (P. alni subsp. multiformis) hybridization events may have been followed by recombination with coding sequences from paralogs present in the putative parental species. In polyploid plants and animals, chromosomal reorganization and gene silencing generally occur rapidly and may be so extensive that the genome is no longer structured as an allopolyploid (31). P. alni subsp. alni still appears structured as an allopolyploid taxon since previous studies demonstrated that P. alni subsp. alni combined the alleles of its progenitors, P. alni subsp. uniformis and P. alni subsp. multiformis, for a series of single-copy genes (15) and for microsatellite loci (16). However, this study demonstrates that, similar to the additivity observed in the ribosomal DNA internal transcribed spacer (3), at least for the allopolyploid taxon P. alni subsp. alni, genetic recombination also occurred between paralogs of the elicitin genes.
Except for the observed new combinations between coding sequences and 3'UTRs, the genomic elicitin pattern of P. alni subsp. alni combined as expected those of the putative parental taxa, P. alni subsp. multiformis and P. alni subsp. uniformis. In addition, while P. alni subsp. uniformis shared several 3'UTRs with P. cambivora, P. alni subsp. alni and P. alni subsp. multiformis possessed private sequences, e.g., a3, a4, and b1. These sequences were found neither in P. cambivora nor in P. fragariae. These findings confirm the close relationship between P. cambivora and P. alni subsp. uniformis (3, 15). Conversely, both P. cambivora and P. fragariae displayed specific elicitin sequences not found in P. alni subsp. alni, in agreement with previous results rejecting the hypothesis that these species may have been the putative parents of the hybrid P. alni subsp. alni (15).
Furthermore, P. alni subsp. alni, P. alni subsp. multiformis, P. alni subsp. uniformis, P. cambivora, and P. fragariae seem to have retained identical genes independently of the speciation. This phenomenon has been partially shown for P. cactorum and P. pseudotsugae (26), but cactorein and pseudotsugaein genes possess distinct 3'UTR sequences (F. Panabières, unpublished data). Jiang et al. (18) demonstrated that the main diversification events of elicitin genes occurred before Phytophthora radiation and that elicitin genes of a given clade, such as the one analyzed in the present work, are under purifying selection. In this respect, the duplication of elicitin genes, creating paralogs, could also explain the multiplicity of divergent 3'UTR sequences for a given elicitin gene, with the assumption of a lower selection pressure on this part of the gene. It is likely that the diversity of elicitin genes in P. cambivora, P. fragariae, and P. alni reflects duplication events prior to the radiation of these species from their common ancestor. It would also indicate that the radiation of these species is of particularly recent origin (7), which seems to be supported by the cross-amplification of microsatellite markers in all three species (16). However, the unexpected conservation of several 3'UTR regions for these different species (e.g., ha1, b2, a8) could also be the outcome of reticulation or introgression events. Gene transfer after speciation should not be ruled out as a possibility to explain the elicitin gene patterns observed within this clade.
 |
ACKNOWLEDGMENTS
|
|---|
This work was supported by grants from INRA, the Agence de l'Eau Rhin-Meuse, and the Institut Fédérateur de Recherches 110.
We thank the colleagues listed in Table 1 for sharing Phytophthora sp. isolates. We are grateful to Everett Hansen for his critical reading of an early version of the manuscript and to Christine Delaruelle for her help in sequencing.
 |
FOOTNOTES
|
|---|
* Corresponding author. Mailing address: INRA, Nancy-Université, UMR1136, Interactions Arbres-Microorganismes, F-54280 Champenoux, France. Phone: 33 383 394056. Fax: 33 383 394069. E-mail: frey{at}nancy.inra.fr 
Published ahead of print on 29 June 2007. 
 |
REFERENCES
|
|---|
- Becker, J., S. Nagel, and R. Tenhaken. 2000. Cloning, expression and characterization of protein elicitors from the soybean pathogenic fungus Phytophthora sojae. J. Phytopathol. 148:161-167.[CrossRef]
- Brasier, C. M. 2003. The hybrid alder Phytophthora: genetic status, pathogenicity, distribution, and competitive survival, p. 39-54. In J. Gibbs, C. Van Dijk, and J. Webber (ed.), Phytophthora disease of alder in Europe. Bulletin 126. Forestry Commission, Edinburgh, United Kingdom.
- Brasier, C. M., D. E. L. Cooke, and J. M. Duncan. 1999. Origin of a new Phytophthora pathogen through interspecific hybridization. Proc. Natl. Acad. Sci. USA 96:5878-5883.[Abstract/Free Full Text]
- Brasier, C. M., and S. Kirk. 2001. Comparative aggressiveness of standard and variant hybrid alder Phytophthora, Phytophthora cambivora and other Phytophthora species on bark of Alnus, Quercus and other woody hosts. Plant Pathol. 50:218-229.[CrossRef]
- Brasier, C. M., S. A. Kirk, J. Delcan, D. E. L. Cooke, T. Jung, and W. A. Man in't Veld. 2004. Phytophthora alni sp. nov. and its variants: designation of emerging heteroploid hybrid pathogens spreading on Alnus trees. Mycol. Res. 108:1172-1184.[CrossRef][Medline]
- Colas, V., S. Conrod, P. Venard, H. Keller, P. Ricci, and F. Panabières. 2001. Elicitin genes expressed in vitro by certain tobacco isolates of Phytophthora parasitica are down regulated during compatible interactions. Mol. Plant-Microbe Interact. 14:326-335.[Medline]
- Cooke, D. E. L., A. Drenth, J. M. Duncan, G. Wagels, and C. M. Brasier. 2000. A molecular phylogeny of Phytophthora and related oomycetes. Fungal Genet. Biol. 30:17-32.[CrossRef][Medline]
- Corpet, F. 1988. Multiple sequence alignment with hierarchical clustering. Nucleic Acids Res. 16:10881-10890.[Abstract/Free Full Text]
- Delcán, J., and C. M. Brasier. 2001. Oospore viability and variation in zoospore and hyphal tip derivatives of the hybrid alder phytophthoras. For. Pathol. 31:65-83.
- De Merlier, D., A. Chandelier, N. DeBruxelles, L. Noldus, F. Laurent, E. Dufays, H. Classens, and M. Cavelier. 2005. Characterization of alder Phytophthora isolates from Wallonia and development of SCAR primers for their specific detection. J. Phytopathol. 153:99-107.[CrossRef]
- Duclos, J., A. Fauconnier, A. C. Coelho, A. Bollen, A. Cravador, and E. Godfroit. 1998. Identification of a gene cluster in Phytophthora cinnamomi. DNA Seq. 9:231-237.[Medline]
- Garcia-Olmedo, F., P. Carbonero, C. Aragancillo, and G. Salcedo. 1978. Loss of redundant gene expression after polyploidization in plants. Experentia 34:332-333.[CrossRef]
- Gibbs, J., C. Van Dijk, and J. Webber. 2003. Phytophthora disease of alder in Europe. Bulletin 126. Forestry Commission, Edinburgh, United Kingdom.
- Gibbs, J. N. 1995. Phytophthora root disease of alder in Britain. EPPO Bull. 25:661-664.
- Ioos, R., A. Andrieux, B. Marçais, and P. Frey. 2006. Genetic characterization of the natural hybrid species Phytophthora alni as inferred from nuclear and mitochondrial DNA analyses. Fungal Genet. Biol. 43:511-529.[CrossRef][Medline]
- Ioos, R., B. Barrès, A. Andrieux, and P. Frey. 2007. Characterization of microsatellite markers in the hybrid Phytophthora alni subsp. alni and cross-amplification with related taxa. Mol. Ecol. Notes 7:133-137.[CrossRef]
- Jiang, R. H., A. L. Dawe, R. Weide, M. van Staveren, S. Peters, D. L. Nuss, and F. Govers. 2005. Elicitin genes in Phytophthora infestans are clustered and interspersed with various transposon-like elements. Mol. Genet. Genomics 273:20-32.[CrossRef][Medline]
- Jiang, R. H. Y., B. M. Tyler, S. C. Whisson, A. R. Hardham, and F. Govers. 2006. Ancient origin of elicitin gene clusters in Phytophthora genomes. Mol. Biol. Evol. 23:338-351.[Abstract/Free Full Text]
- Man in't Veld, W. A., A. W. A. M. de Cock, and R. C. Summerbell. 2007. Natural hybrids of resident and introduced Phytophthora species proliferating on new hosts. Eur. J. Plant Pathol. 117:25-33.[CrossRef]
- Man in't Veld, W. A., W. J. Veenbaas-Rijk, E. Ilieva, A. W. A. M. De Cock, P. J. M. Bonants, and R. Pieters. 1998. Natural hybrids of Phytophthora nicotianae and P. cactorum demonstrated by isozyme analysis and random amplified polymorphic DNA. Phytopathology 88:922-929.[CrossRef]
- Miller, P. M. 1955. V-8 juice agar as a general purpose medium for fungi and bacteria. Phytopathology 45:461-462.
- Nagy, Z. A., J. Bakonyi, and T. Ersek. 2003. Standard and Swedish variant types of the hybrid alder Phytophthora attacking alder in Hungary. Pest Manag. Sci. 59:484-492.[CrossRef][Medline]
- Panabières, F., J. Amselem, E. Galiana, and J. Y. Le Berre. 2005. Gene identification in the oomycete pathogen Phytophthora parasitica during in vitro vegetative growth through expressed sequence tags. Fungal Genet. Biol. 42:611-623.[CrossRef][Medline]
- Panabières, F., A. Marais, J. Y. Le Berre, I. Penot, D. Fournier, and P. Ricci. 1995. Characterization of a gene cluster of Phytophthora cryptogea which codes for elicitins, proteins inducing a hypersensitive-like response in tobacco. Mol. Plant-Microbe Interact. 8:996-1003.[Medline]
- Panabières, F., M. Ponchet, V. Allasia, L. Cardin, and P. Ricci. 1997. Characterization of border species among Pythiaceae: several Pythium isolates produce elicitins, typical proteins from Phytophthora spp. Mycol. Res. 101:1459-1468.[CrossRef]
- Ponchet, M., F. Panabières, M. L. Milat, J. L. Montillet, L. Suty, C. Triantaphylides, Y. Tirilly, and J. P. Blein. 1999. Are elicitins cryptogams in plant-oomycete communications? Cell. Mol. Life Sci. 56:1020-1047.[CrossRef][Medline]
- Qutob, D., E. Huitema, M. Gijzen, and S. Kamoun. 2003. Variation in structure and activity among elicitins from Phytophthora sojae. Mol. Plant Pathol. 4:119-124.[CrossRef]
- Robin, C., M. L. Desprez-Loustau, G. Capron, and C. Delatour. 1998. First record of Phytophthora cinnamomi on cork and holm oak in France and evidence of pathogenicity. Ann. Sci. For. 55:869-883.[CrossRef]
- Rosenthal, A., O. Coutelle, and M. Craxton. 1993. Large-scale production of DNA sequencing templates by microtitre format PCR. Nucleic Acids Res. 21:173-174.[Free Full Text]
- Santini, A., G. P. Barzanti, and P. Capretti. 2003. Susceptibility of some mesophilic hardwoods to alder Phytophthora. J. Phytopathol. 151:406-410.[CrossRef]
- Soltis, D. E., and P. S. Soltis. 1999. Polyploidy: recurrent formation and genome evolution. Trends Ecol. Evol. 9:348-352.
- Swofford, D. L. 2002. PAUP*. Phylogenetic analysis using parsimony (*and other methods), version 4.0b10. Sinauer Associates, Sunderland, MA.
- Thompson, J. D., D. G. Higgins, and T. J. Gibson. 1994. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 22:4673-4680.[Abstract/Free Full Text]
- Volkov, R. A., N. V. Borisjuk, I. I. Panchuk, D. Schweizer, and V. Hemleben. 1999. Elimination and rearrangement of parental rDNA in the allotetraploid Nicotiana tabacum. Mol. Biol. Evol. 16:311-320.[Abstract]
Applied and Environmental Microbiology, September 2007, p. 5587-5597, Vol. 73, No. 17
0099-2240/07/$08.00+0 doi:10.1128/AEM.00721-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.