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.

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
|
|
|---|
|
|
|---|
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.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. List of the Phytophthora spp. and Pythium spp. used in this study
|
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 a new window] |
TABLE 2. Characteristics of all the different cDNA and genomic sequences obtained and assignation to a specific 3'UTR group
|
|
View this table: [in a new window] |
TABLE 3. List and sequences of the 3'UTR-specific reverse primers designed in this study
|
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.
|
|
|---|
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 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 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).
|
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 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).
|
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 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.
|
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.
|
|
|---|
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.
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.
Published ahead of print on 29 June 2007. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»