Computational and comparative analyses of 150 full-length cDNA sequences from the oomycete plant pathogen Phytophthora infestans

Phytophthora infestans is a devastating phytopathogenic oomycete that causes late blight on tomato and potato. Recent genome sequencing efforts of P. infestans and other Phytophthora species are generating vast amounts of sequence data providing opportunities to unlock the complex nature of pathogenesis. However, accurate annotation of Phytophthora genomes will be a significant challenge. Most of the information about gene structure in these species was gathered from a handful of genes resulting in significant limitations for development of ab initio gene-calling programs. In this study, we collected a total of 150 bioinformatically determined near full-length cDNA (FLcDNA) sequences of P. infestans that were predicted to contain full open reading frame sequences. We performed detailed computational analyses of these FLcDNA sequences to obtain a snapshot of P. infestans gene structure, gauge the degree of sequence conservation between P. infestans genes and those of Phytophthora sojae and Phytophthora ramorum, and identify patterns of gene conservation between P. infestans and various eukaryotes, particularly fungi, for which genome-wide translated protein sequences are available. These analyses helped us to define the structural characteristics of P. infestans genes using a validated data set. We also determined the degree of sequence conservation within the genus Phytophthora and identified a set of fast evolving genes. Finally, we identified a set of genes that are shared between Phytophthora and fungal phytopathogens but absent in animal fungal pathogens. These results confirm that plant pathogenic oomycetes and fungi share virulence components, and suggest that eukaryotic microbial pathogens that share similar lifestyles also share a similar set of genes independently of their phylogenetic relatedness.     

The repertoire of transfer RNA genes is tuned to codon usage bias in the genomes of Phytophthora sojae and Phytophthora ramorum

In all, 238 and 155 transfer (t)RNA genes were predicted from the genomes of Phytophthora sojae and P. ramorum, respectively. After omitting pseudogenes and undetermined types of tRNA genes, there remained 208 P. sojae tRNA genes and 140 P. ramorum tRNA genes. There were 45 types of tRNA genes, with distinct anticodons, in each species. Fourteen common anticodon types of tRNAs are missing altogether from the genome in the two species; however, these appear to be compensated by wobbling of other tRNA anticodons in a manner which is tied to the codon bias in Phytophthora genes. The most abundant tRNA class was arginine in both P. sojae and P. ramorum. A codon usage table was generated for these two organisms from a total of 9,803,525 codons in P. sojae and 7,496,598 codons in P. ramorum. The most abundant codon type detected from the codon usage tables was GAG (encoding glutamic acid), whereas the most numerous tRNA gene had a methionine anticodon (CAT). The correlation between the frequencies of tRNA genes and the codon frequencies in protein-coding genes was very low (0.12 in P. sojae and 0.19 in P. ramorum); however, the correlation between amino acid tRNA gene frequency and the corresponding amino acid codon frequency in P. sojae and P. ramorum was substantially higher (0.53 in P. sojae and 0.77 in P. ramorum). The codon usage frequencies of P. sojae and P ramorum were very strongly correlated (0.99), as were tRNA gene frequencies (0.77). Approximately 60% of orthologous tRNA gene pairs in P sojae and P. ramorum are located in regions that have conserved synteny in the two species.     

Gene duplication event in family 12 glycosyl hydrolase from Phytophthora spp

A total of 18 paralogs of xyloglucan-specific endoglucanases (EGLs) from the glycosyl hydrolase family 12 were identified and characterized in Phytophthora sojae and Phytophthora ramorum. These genes encode predicted extracellular enzymes, with sizes ranging from 189 to 435 amino acid residues, that would be capable of hydrolyzing the xyloglucan component of the host cell wall. In two cases, four and six functional copies of these genes were found in tight succession within a region of 5 and 18 kb, respectively. The overall gene copy number and relative organization appeared well conserved between P. sojae and P. ramorum, with apparent synteny in this region of their respective genomes. Phylogenetic analyses of Phytophthora endoglucanases of family 12 and other known members of EGL 12, revealed a close relatedness with a fairly conserved gene sub-family containing, among others, sequences from the fungi Emericella desertorum and Aspergillus aculeatus. This is the first report of family 12 EGLs present in plant pathogenic eukaryotes.     

Targeted gene mutation in Phytophthora spp

The genus Phytophthora belongs to the oomycetes and is composed of plant pathogens. Currently, there are no strategies to mutate specific genes for members of this genus. Whole genome sequences are available or being prepared for Phytophthora sojae, P. ramorum, P. infestans, and P. capsici and the development of molecular biological techniques for functional genomics is encouraged. This article describes the adaptation of the reverse-genetic strategy of targeting induced local lesions in genomes (TILLING) to isolate gene-specific mutants in Phytophthora spp. A genomic library of 2,400 ethylnitrosourea (ENU) mutants of P. sojae was created and screened for induced point mutations in the genes encoding a necrosisinducing protein (PsojNIP) and a Phytophthora-specific phospholipase D (PsPXTM-PLD). Mutations were detected in single individuals and included silent, missense, and nonsense changes. Homozygous mutant isolates carrying a potentially deleterious missense mutation in PsojNIP and a premature stop codon in PsPXTM-PLD were identified. No phenotypic effect has yet been found for the homozygous mutant of PsojNIP. For those of PsPXTM-PLD, a reduction in growth rate and an appressed mycelial growth was observed. This demonstrates the feasibility of target-selected gene disruption for Phytophthora spp. and adds an important tool for functional genomic investigation.     

Ancient origin of elicitin gene clusters in Phytophthora genomes

The genus Phytophthora belongs to the oomycetes in the eukaryotic stramenopile lineage and is comprised of over 65 species that are all destructive plant pathogens on a wide range of dicotyledons. Phytophthora produces elicitins (ELIs), a group of extracellular elicitor proteins that cause a hypersensitive response in tobacco. Database mining revealed several new classes of elicitin-like (ELL) sequences with diverse elicitin domains in Phytophthora infestans, Phytophthora sojae, Phytophthora brassicae, and Phytophthora ramorum. ELIs and ELLs were shown to be unique to Phytophthora and Pythium species. They are ubiquitous among Phytophthora species and belong to one of the most highly conserved and complex protein families in the Phytophthora genus. Phylogeny construction with elicitin domains derived from 156 ELIs and ELLs showed that most of the diversified family members existed prior to divergence of Phytophthora species from a common ancestor. Analysis to discriminate diversifying and purifying selection showed that all 17 ELI and ELL clades are under purifying selection. Within highly similar ELI groups there was no evidence for positively selected amino acids suggesting that purifying selection contributes to the continued existence of this diverse protein family. Characteristic cysteine spacing patterns were found for each phylogenetic clade. Except for the canonical clade ELI-1, ELIs and ELLs possess C-terminal domains of variable length, many of which have a high threonine, serine, or proline content suggesting an association with the cell wall. In addition, some ELIs and ELLs have a predicted glycosylphosphatidylinositol site suggesting anchoring of the C-terminal domain to the cell membrane. The eli and ell genes belonging to different clades are clustered in the genomes. Overall, eli and ell genes are expressed at different levels and in different life cycle stages but those sharing the same phylogenetic clade appear to have similar expression patterns.     

Expressed peptide tags: an additional layer of data for genome annotation

While genome sequencing is becoming ever more routine, genome annotation remains a challenging process. Identification of the coding sequences within the genomic milieu presents a tremendous challenge, especially for eukaryotes with their complex gene architectures. Here, we present a method to assist the annotation process through the use of proteomic data and bioinformatics. Mass spectra of digested protein preparations of the organism of interest were acquired and searched against a protein database created by a six-frame translation of the genome. The identified peptides were mapped back to the genome, compared to the current annotation, and then categorized as supporting or extending the current genome annotation. We named the classified peptides Expressed Peptide Tags (EPTs). The well-annotated bacterium Rhodopseudomonas palustris was used as a control for the method and showed a high degree of correlation between EPT mapping and the current annotation, with 86% of the EPTs confirming existing gene calls and less than 1% of the EPTs expanding on the current annotation. The eukaryotic plant pathogens Phytophthora ramorum and Phytophthora sojae, whose genomes have been recently sequenced and are much less well-annotated, were also subjected to this method. A series of algorithmic steps were taken to increase the confidence of EPT identification for these organisms, including generation of smaller subdatabases to be searched against, and definition of EPT criteria that accommodates the more complex eukaryotic gene architecture. As expected, the analysis of the Phytophthora species showed less correlation between EPT mapping and their current annotation. While approximately 76% of Phytophthora EPTs supported the current annotation, a portion of them (7.7% and 12.9% for P. ramorum and P. sojae, respectively) suggested modification to current gene calls or identified novel genes that were missed by the current genome annotation of these organisms.     

Ancient isolation and independent evolution of the three clonal lineages of the exotic sudden oak death pathogen Phytophthora ramorum

The genus Phytophthora includes some of the most destructive plant pathogens affecting agricultural and native ecosystems and is responsible for a number of recent emerging and re-emerging infectious diseases of plants. Sudden oak death, caused by the exotic pathogen P. ramorum , has caused extensive mortality of oaks and tanoaks in Northern California, and has brought economic losses to US and European nurseries as well due to its infection of common ornamental plants. In its known range, P. ramorum occurs as three distinct clonal lineages. We inferred the evolutionary history of P. ramorum from nuclear sequence data using coalescent-based approaches. We found that the three lineages have been diverging for at least 11% of their history, an evolutionarily significant amount of time estimated to be on the order of 165 000 to 500 000 years. There was also strong evidence for historical recombination between the lineages, indicating that the ancestors of the P. ramorum lineages were members of a sexually reproducing population. Due to this recombination, the ages of the lineages varied within and between loci, but coalescent analyses suggested that the European lineage may be older than the North American lineages. The divergence of the three clonal lineages of P. ramorum supports a scenario in which the three lineages originated from different geographic locations that were sufficiently isolated from each other to allow independent evolution prior to introduction to North America and Europe. It is thus probable that the emergence of P. ramorum in North America and Europe was the result of three independent migration events.     

Cross-species global proteomics reveals conserved and unique processes in Phytophthora sojae and Phytophthora ramorum

Phytophthora ramorum and Phytophthora sojae are destructive plant pathogens. P. sojae has a narrow host range, whereas P. ramorum has a wide host range. A global proteomics comparison of the vegetative (mycelium) and infective (germinating cyst) life stages of P. sojae and P. ramorum was conducted to identify candidate proteins involved in host range, early infection, and vegetative growth. Sixty-two candidates for early infection, 26 candidates for vegetative growth, and numerous proteins that may be involved in defining host specificity were identified. In addition, common life stage proteomic trends between the organisms were observed. In mycelia, proteins involved in transport and metabolism of amino acids, carbohydrates, and other small molecules were up-regulated. In the germinating cysts, up-regulated proteins associated with lipid transport and metabolism, cytoskeleton, and protein synthesis were observed. It appears that the germinating cyst catabolizes lipid reserves through the beta-oxidation pathway to drive the extensive protein synthesis necessary to produce the germ tube and initiate infection. Once inside the host, the pathogen switches to vegetative growth in which energy is derived from glycolysis and utilized for synthesis of amino acids and other molecules that assist survival in the plant tissue.     

一种土壤中大豆疫霉菌分离新方法*

由于腐霉菌的干扰,土壤中大豆疫霉菌的分离十分困难。利用大豆疫霉菌的致病性和大豆对病原菌的选择作用排除腐霉菌,我们建立了一种简单、有效的土壤中大豆疫霉菌的分离方法。该方法用不含抗大豆疫霉根腐病基因的大豆叶碟诱钓大豆疫霉菌的游动孢子,将诱钓叶碟直接接种不含抗大豆疫霉菌基因的大豆植株,再对病株进行选择性或非选择性分离获得大豆疫霉菌。此方法能十分有效地排除腐霉菌干扰和细菌的污染,直接获得纯化菌株。应用该方法我们在以前未报道有大豆疫霉根腐病发生的山东、河南、安徽、江苏和浙江分离到大豆疫霉菌。     

Common infection strategies of pathogenic eukaryotes

Pathogenic eukaryotes belong to several distinct phylogenetic lineages and have evolved the ability to colonize a range of hosts, including animals and plants. Pathogenic lifestyles have evolved repeatedly in eukaryotes, indicating that unique molecular processes are involved in host infection. However, evidence is now emerging that divergent eukaryotic pathogens might share common mechanisms of pathogenicity. The results from recent studies demonstrate that Plasmodium falciparum and Phytophthora infestans use equivalent host-targeting signals to deliver virulence adhesins and avirulence gene products into human and plant cells, respectively. Remodelling of host cells by different eukaryotic pathogens might therefore share some common features.     

The malarial host-targeting signal is conserved in the Irish potato famine pathogen

Animal and plant eukaryotic pathogens, such as the human malaria parasite Plasmodium falciparum and the potato late blight agent Phytophthora infestans, are widely divergent eukaryotic microbes. Yet they both produce secretory virulence and pathogenic proteins that alter host cell functions. In P. falciparum, export of parasite proteins to the host erythrocyte is mediated by leader sequences shown to contain a host-targeting (HT) motif centered on an RxLx (E, D, or Q) core: this motif appears to signify a major pathogenic export pathway with hundreds of putative effectors. Here we show that a secretory protein of P. infestans, which is perceived by plant disease resistance proteins and induces hypersensitive plant cell death, contains a leader sequence that is equivalent to the Plasmodium HT-leader in its ability to export fusion of green fluorescent protein (GFP) from the P. falciparum parasite to the host erythrocyte. This export is dependent on an RxLR sequence conserved in P. infestans leaders, as well as in leaders of all ten secretory oomycete proteins shown to function inside plant cells. The RxLR motif is also detected in hundreds of secretory proteins of P. infestans, Phytophthora sojae, and Phytophthora ramorum and has high value in predicting host-targeted leaders. A consensus motif further reveals E/D residues enriched within approximately 25 amino acids downstream of the RxLR, which are also needed for export. Together the data suggest that in these plant pathogenic oomycetes, a consensus HT motif may reside in an extended sequence of approximately 25-30 amino acids, rather than in a short linear sequence. Evidence is presented that although the consensus is much shorter in P. falciparum, information sufficient for vacuolar export is contained in a region of approximately 30 amino acids, which includes sequences flanking the HT core. Finally, positional conservation between Phytophthora RxLR and P. falciparum RxLx (E, D, Q) is consistent with the idea that the context of their presentation is constrained. These studies provide the first evidence to our knowledge that eukaryotic microbes share equivalent pathogenic HT signals and thus conserved mechanisms to access host cells across plant and animal kingdoms that may present unique targets for prophylaxis across divergent pathogens.     

Emergence of the sudden oak death pathogen Phytophthora ramorum

The recently emerged plant pathogen Phytophthora ramorum is responsible for causing the sudden oak death epidemic. This review documents the emergence of P. ramorum based on evolutionary and population genetic analyses. Currently infection by P. ramorum occurs only in Europe and North America and three clonal lineages are distinguished: EU1, NA1 and NA2. Ancient divergence of these lineages supports a scenario in which P. ramorum originated from reproductively isolated populations and underwent at least four global migration events. This recent work sheds new light on mechanisms of emergence of exotic pathogens and provides crucial insights into migration pathways.     

Genomewide analysis of phospholipid signaling genes in Phytophthora spp.: novelties and a missing link

Phospholipids are cellular membrane components in eukaryotic cells that execute many important roles in signaling. Genes encoding enzymes required for phospholipid signaling and metabolism have been characterized in several organisms, but only a few have been described for oomycetes. In this study, the genome sequences of Phytophthora sojae and P. ramorum were explored to construct a comprehensive genomewide inventory of genes involved in the most universal phospholipid signaling pathways. Several genes and gene families were annotated, including those encoding phosphatidylinositol synthase (PIS), phosphatidylinositol (phosphate) kinase (PI[P]K), diacylglycerol kinase (DAG), and phospholipase D (PLD). The most obvious missing link is a gene encoding phospholipase C (PLC). In all eukaryotic genomes sequenced to date, PLC genes are annotated based on certain conserved features; however, these genes seem to be absent in Phytophthora spp. Analysis of the structural and regulatory domains and domain organization of the predicted isoforms of PIS, PIK, PIPK, DAG, and PLD revealed many novel features compared with characterized representatives in other eukaryotes. Examples are transmembrane proteins with a C-terminal catalytic PLD domain, secreted PLD-like proteins, and PIPKs that have an N-terminal G-protein-coupled receptor-transmembrane signature. Compared with other sequenced eukaryotes, the genus Phytophthora clearly has several exceptional features in its phospholipid-modifying enzymes.     

Mosaic origin of the heme biosynthesis pathway in photosynthetic eukaryotes

Heme biosynthesis represents one of the most essential metabolic pathways in living organisms, providing the precursors for cytochrome prosthetic groups, photosynthetic pigments, and vitamin B(12). Using genomic data, we have compared the heme pathway in the diatom Thalassiosira pseudonana and the red alga Cyanidioschyzon merolae to those of green algae and higher plants, as well as to those of heterotrophic eukaryotes (fungi, apicomplexans, and animals). Phylogenetic analyses showed the mosaic character of this pathway in photosynthetic eukaryotes. Although most of the algal and plant enzymes showed the expected plastid (cyanobacterial) origin, at least one of them (porphobilinogen deaminase) appears to have a mitochondrial (alpha-proteobacterial) origin. Another enzyme, glutamyl-tRNA synthase, obviously originated in the eukaryotic nucleus. Because all the plastid-targeted sequences consistently form a well-supported cluster, this suggests that genes were either transferred from the primary endosymbiont (cyanobacteria) to the primary host nucleus shortly after the primary endosymbiotic event or replaced with genes from other sources at an equally early time, i.e., before the formation of three primary plastid lineages. The one striking exception to this pattern is ferrochelatase, the enzyme catalyzing the first committed step to heme and bilin pigments. In this case, two red algal sequences do not cluster either with the other plastid sequences or with cyanobacterial sequences and appear to have a proteobacterial origin like that of the apicomplexan parasites Plasmodium and Toxoplasma. Although the heterokonts also acquired their plastid via secondary endosymbiosis from a red alga, the diatom has a typical plastid-cyanobacterial ferrochelatase. We have not found any remnants of the plastidlike heme pathway in the nonphotosynthetic heterokonts Phytophthora ramorum and Phytophthora sojae.     

Comparative analysis of Phytophthora genes encoding secreted proteins reveals conserved synteny and lineage-specific gene duplications and deletions

Comparative analysis of two Phytophthora genomes revealed overall colinearity in four genomic regions consisting of a 1.5-Mb sequence of Phytophthora sojae and a 0.9-Mb sequence of P. ramorum. In these regions with conserved synteny, the gene order is largely similar; however, genome rearrangements also have occurred. Deletions and duplications often were found in association with genes encoding secreted proteins, including effectors that are important for interaction with host plants. Among secreted protein genes, different evolutionary patterns were found. Elicitin genes that code for a complex family of highly conserved Phytophthora-specific elicitors show conservation in gene number and order, and often are clustered. In contrast, the race-specific elicitor gene Avrlb-1 appeared to be missing from the region with conserved synteny, as were its five homologs that are scattered over the four genomic regions. Some gene families encoding secreted proteins were found to be expanded in one species compared with the other. This could be the result of either repeated gene duplications in one species or specific deletions in the other. These different evolutionary patterns may shed light on the functions of these secreted proteins in the biology and pathology of the two Phytophthora spp.     

A functional screen to characterize the secretomes of eukaryotic pathogens and their hosts in planta

Complex suites of proteins that are secreted by plants and phytopathogens into the plant apoplast play crucial roles in surveillance, assault, defense, and counter-defense. High-throughput genome-scale strategies are being developed to better understand the nature of these "secretomes" and the identity of pathogen-derived effector proteins that subvert plant defenses and promote pathogenicity. Although combined bioinformatic and experimental approaches recently have provided comprehensive coverage of secreted proteins from bacterial phytopathogens, far less is known about the secretomes and batteries of effectors of eukaryotic phytopathogens; notably fungi and oomycetes. The yeast secretion trap (YST) represents a potentially valuable technique to simultaneously target pathogen and host secretomes in infected plant material. A YST screen, using a new vector system, was applied to study the interaction between tomato (Solanum lycopersicum) and the oomycete Phytophthora infestans, revealing sets of genes encoding secreted proteins from both pathogen and host. Most of those from the oomycete had no identifiable function and were detectable in planta only during pathogenesis, underlining the value of YST as a tool to identify new candidate effectors and pathogenicity factors. In addition, the majority of the P. infestans proteins had homologs in the genomes of the related oomycetes R. sojae and P. ramorum.     

A multi-locus phylogeny for Phytophthora utilizing markers derived from complete genome sequences.

Phytophthora species are devastating plant pathogens in both agricultural and natural environments. Due to their significant economic and environmental impact, there has been increasing interest in Phytophthora genetics and genomics, culminating in the recent release of three complete genome sequences (P. ramorum, P. sojae, and P. infestans). In this study, genome and other large sequence databases were used to identify over 225 potential genetic markers for phylogenetic analyses. Here, we present a genus-wide phylogeny for 82 Phytophthora species using seven of the most informative loci (approximately 8700 nucleotide sites). Our results support the division of the genus into 10 well-supported clades. The relationships among these clades were rigorously evaluated using a number of phylogenetic methods. This is the most comprehensive study of Phytophthora relationships to date, and many newly discovered species have been included. A more resolved phylogeny of Phytophthora species will allow for better interpretations of the overall evolutionary history of the genus.     

Recognition of Phytophthora infestans Avr4 by potato R4 is triggered by C-terminal domains comprising W motifs

Oomycete RXLR-dEER effector proteins are rapidly evolving proteins with the selective pressure targeted predominantly at their C-terminal ends. The majority of RXLR-dEER proteins have recognizable motifs of 21–30 amino acids in the C-terminal domain that are named after conserved amino acid residues at fixed positions within the respective motifs. In this article, it is reported that the Phytophthora infestans RXLR-dEER protein Avr4 contains three W motifs and one Y motif in its C-terminal domain. Agroinfection assays using constructs encoding modified forms of PiAvr4 have shown that the region containing the W2 motif, in combination with either the W1 or W3 motif, triggers a necrotic response in potato plants carrying the resistance gene R4 . By mining the superfamily of avirulence homologues (Avh) deduced from three sequenced Phytophthora genomes, several Avh proteins were identified as homologues of PiAvr4: six in P. infestans , one in P. ramorum and seven in P. sojae . One very close homologue of PiAvr4 was cloned from the sibling species, P. mirabilis. This species is not pathogenic on potato but, similar to PiAvr4, PmirAvh4 triggered a necrotic response on potato clones carrying R4 , but not on clones lacking R4 . Genes encoding RXLR-dEER effectors are often located in regions showing genome rearrangements. Alignment of the genomic region harbouring PiAvr4 with syntenic regions in P. sojae and P. ramorum revealed that PiAvr4 is located on a 100-kb indel block and is surrounded by transposable elements.     

A multi-locus phylogeny for Phytophthora utilizing markers derived from complete genome sequences

Phytophthora species are devastating plant pathogens in both agricultural and natural environments. Due to their significant economic and environmental impact, there has been increasing interest in Phytophthora genetics and genomics, culminating in the recent release of three complete genome sequences (P. ramorum, P. sojae, and P. infestans). In this study, genome and other large sequence databases were used to identify over 225 potential genetic markers for phylogenetic analyses. Here, we present a genus-wide phylogeny for 82 Phytophthora species using seven of the most informative loci (approximately 8700 nucleotide sites). Our results support the division of the genus into 10 well-supported clades. The relationships among these clades were rigorously evaluated using a number of phylogenetic methods. This is the most comprehensive study of Phytophthora relationships to date, and many newly discovered species have been included. A more resolved phylogeny of Phytophthora species will allow for better interpretations of the overall evolutionary history of the genus.