Multiple horizontal gene transfer events and domain fusions have created novel regulatory and metabolic networks in the oomycete genome

Genome-Wide Identification and Evolutionary Analyses of SrfA Operon Genes in Bacillus

A variety of secondary metabolites contributing to plant growth are synthesized by bacterial nonribosomal peptide synthases (NRPSs). Among them, the NRPS biosynthesis of surfactin is regulated by the SrfA operon. To explore the molecular mechanism for the diversity of surfactins produced by bacteria within the genus Bacillus, we performed a genome-wide identification study focused on three critical genes of the SrfA operon-SrfAA, SrfAB and SrfAC-from 999 Bacillus genomes (belonging to 47 species). Gene family clustering indicated the three genes can be divided into 66 orthologous groups (gene families), of which a majority comprised members of multiple genes (e.g., OG0000009 had members of all three SrfAA, SrfAB and SrfAC genes), indicating high sequence similarity among the three genes. Phylogenetic analyses also found that none of the three genes formed monophyletic groups, but were usually arranged in a mixed manner, suggesting the close evolutionary relationship among the three genes. Considering the module structure of the three genes, we propose that self-duplication, especially tandem duplications, might have contributed to the initial establishment of the entire SrfA operon, and further gene fusion and recombination as well as accumulated mutations might have continuously shaped the different functional roles of SrfAA, SrfAB and SrfAC. Overall, this study provides novel insight into metabolic gene clusters and operon evolution in bacteria.

Draft genome assemblies for tree pathogens Phytophthora pseudosyringae and Phytophthora boehmeriae

Species of Phytophthora, plant pathogenic eukaryotic microbes, can cause disease on many tree species. Genome sequencing of species from this genus has helped to determine components of their pathogenicity arsenal. Here, we sequenced genomes for two widely distributed species, Phytophthora pseudosyringae and Phytophthora boehmeriae, yielding genome assemblies of 49 and 40 Mb, respectively. We identified more than 270 candidate disease promoting RXLR effector coding genes for each species, and hundreds of genes encoding candidate plant cell wall degrading carbohydrate active enzymes (CAZymes). These data boost genome sequence representation across the Phytophthora genus, and form resources for further study of Phytophthora pathogenesis.

Aphanomyces astaci mtDNA: insights into the pathogen's differentiation and its genetic diversity from other closely related oomycetes

The causative agent of crayfish plague, Aphanomyces astaci (Saprolegniales, Oomycota), is one of the 100 world's worst invasive alien species and represents a major threat to freshwater crayfish species worldwide. A better understanding of the biology and epidemiology of A. astaci relies on the application of efficient tools to detect the pathogen and assess its genetic diversity. In this study, we validated the specificity of two recently developed PCR-based approaches used to detect A. astaci groups. The first relies on the analysis of mitochondrial ribosomal rnnS (small) and rnnL (large) subunit sequences and the second, of sequences obtained by using genotype-specific primers designed from A. astaci whole genome sequencing. For this purpose, we tested the specificity against 76 selected isolates, including other oomycete species and the recently described species Aphanomyces fennicus, which, when used in nrITS-based specific tests for A. astaci, is known to result in a false positive. Under both approaches, we were able to efficiently and accurately identify A. astaci and its genetic groups in both pure cultures and clinical samples. We report that sequence analysis of the rnnS region alone is sufficient for the identification of A. astaci and a partial characterization of haplogroups. In contrast, the rnnL region alone is not sufficiently informative for A. astaci identification as other oomycete species present sequences identical to those of A. astaci.

Structural and catalytic analysis of two diverse uridine phosphorylases in Phytophthora capsici

Uridine phosphorylase (UP) is a key enzyme of pyrimidine salvage pathways that enables the recycling of endogenous or exogenous-supplied pyrimidines and plays an important intracellular metabolic role. Here, we biochemically and structurally characterized two evolutionarily divergent uridine phosphorylases, PcUP1 and PcUP2 from the oomycete pathogen Phytophthora capsici. Our analysis of other oomycete genomes revealed that both uridine phosphorylases are present in Phytophthora and Pythium genomes, but only UP2 is seen in Saprolegnia spp. which are basal members of the oomycetes. Moreover, uridine phosphorylases are not found in obligate oomycete pathogens such as Hyaloperonospora arabidopsidis and Albugo spp. PcUP1 and PcUP2 are upregulated 300 and 500 fold respectively, within 90 min after infection of pepper leaves. The crystal structures of PcUP1 in ligand-free and in complex with uracil/ribose-1-phosphate, 2'-deoxyuridine/phosphate and thymidine/phosphate were analyzed. Crystal structure of this uridine phosphorylase showed strict conservation of key residues in the binding pocket. Structure analysis of PcUP1 with bound ligands, and site-directed mutagenesis of key residues provide additional support for the "push-pull" model of catalysis. Our study highlights the importance of pyrimidine salvage during the earliest stages of infection.

Activity and Phylogenetics of the Broadly Occurring Family of Microbial Nep1-Like Proteins

Necrosis- and ethylene-inducing peptide 1 (Nep1)-like proteins (NLP) have an extremely broad taxonomic distribution; they occur in bacteria, fungi, and oomycetes. NLPs come in two forms, those that are cytotoxic to eudicot plants and those that are noncytotoxic. Cytotoxic NLPs bind to glycosyl inositol phosphoryl ceramide (GIPC) sphingolipids that are abundant in the outer leaflet of plant plasma membranes. Binding allows the NLP to become cytolytic in eudicots but not monocots. The function of noncytotoxic NLPs remains enigmatic, but the expansion of NLP genes in oomycete genomes suggests they are important. Several plant species have evolved the capacity to recognize NLPs as molecular patterns and trigger plant immunity, e.g., Arabidopsis thaliana detects nlp peptides via the receptor-like protein RLP23. In this review, we provide a historical perspective from discovery to understanding of molecular mechanisms and describe the latest developments in the NLP field to shed light on these fascinating microbial proteins.

Genome-wide sequencing and metabolic annotation of Pythium irregulare CBS 494.86: understanding Eicosapentaenoic acid production

BACKGROUND

Pythium irregulare is an oleaginous Oomycete able to accumulate large amounts of lipids, including Eicosapentaenoic acid (EPA). EPA is an important and expensive dietary supplement with a promising and very competitive market, which is dependent on fish-oil extraction. This has prompted several research groups to study biotechnological routes to obtain specific fatty acids rather than a mixture of various lipids. Moreover, microorganisms can use low cost carbon sources for lipid production, thus reducing production costs. Previous studies have highlighted the production of EPA by P. irregulare, exploiting diverse low cost carbon sources that are produced in large amounts, such as vinasse, glycerol, and food wastewater. However, there is still a lack of knowledge about its biosynthetic pathways, because no functional annotation of any Pythium sp. exists yet. The goal of this work was to identify key genes and pathways related to EPA biosynthesis, in P. irregulare CBS 494.86, by sequencing and performing an unprecedented annotation of its genome, considering the possibility of using wastewater as a carbon source.

RESULTS

Genome sequencing provided 17,727 candidate genes, with 3809 of them associated with enzyme code and 945 with membrane transporter proteins. The functional annotation was compared with curated information of oleaginous organisms, understanding amino acids and fatty acids production, and consumption of carbon and nitrogen sources, present in the wastewater. The main features include the presence of genes related to the consumption of several sugars and candidate genes of unsaturated fatty acids production.

CONCLUSIONS

The whole metabolic genome presented, which is an unprecedented reconstruction of P. irregulare CBS 494.86, shows its potential to produce value-added products, in special EPA, for food and pharmaceutical industries, moreover it infers metabolic capabilities of the microorganism by incorporating information obtained from literature and genomic data, supplying information of great importance to future work.

Comparative genomic analysis of the 'pseudofungus' Hyphochytrium catenoides

Eukaryotic microbes have three primary mechanisms for obtaining nutrients and energy: phagotrophy, photosynthesis and osmotrophy. Traits associated with the latter two functions arose independently multiple times in the eukaryotes. The Fungi successfully coupled osmotrophy with filamentous growth, and similar traits are also manifested in the Pseudofungi (oomycetes and hyphochytriomycetes). Both the Fungi and the Pseudofungi encompass a diversity of plant and animal parasites. Genome-sequencing efforts have focused on host-associated microbes (mutualistic symbionts or parasites), providing limited comparisons with free-living relatives. Here we report the first draft genome sequence of a hyphochytriomycete ‘pseudofungus’; Hyphochytrium catenoides. Using phylogenomic approaches, we identify genes of recent viral ancestry, with related viral derived genes also present on the genomes of oomycetes, suggesting a complex history of viral coevolution and integration across the Pseudofungi. H. catenoides has a complex life cycle involving diverse filamentous structures and a flagellated zoospore with a single anterior tinselate flagellum. We use genome comparisons, drug sensitivity analysis and high-throughput culture arrays to investigate the ancestry of oomycete/pseudofungal characteristics, demonstrating that many of the genetic features associated with parasitic traits evolved specifically within the oomycete radiation. Comparative genomics also identified differences in the repertoire of genes associated with filamentous growth between the Fungi and the Pseudofungi, including differences in vesicle trafficking systems, cell-wall synthesis pathways and motor protein repertoire, demonstrating that unique cellular systems underpinned the convergent evolution of filamentous osmotrophic growth in these two eukaryotic groups.

Rethinking the evolution of eukaryotic metabolism: novel cellular partitioning of enzymes in stramenopiles links serine biosynthesis to glycolysis in mitochondria

Background

An important feature of eukaryotic evolution is metabolic compartmentalization, in which certain pathways are restricted to the cytosol or specific organelles. Glycolysis in eukaryotes is described as a cytosolic process. The universality of this canon has been challenged by recent genome data that suggest that some glycolytic enzymes made by stramenopiles bear mitochondrial targeting peptides.

Results

Mining of oomycete, diatom, and brown algal genomes indicates that stramenopiles encode two forms of enzymes for the second half of glycolysis, one with and the other without mitochondrial targeting peptides. The predicted mitochondrial targeting was confirmed by using fluorescent tags to localize phosphoglycerate kinase, phosphoglycerate mutase, and pyruvate kinase in Phytophthora infestans, the oomycete that causes potato blight. A genome-wide search for other enzymes with atypical mitochondrial locations identified phosphoglycerate dehydrogenase, phosphoserine aminotransferase, and phosphoserine phosphatase, which form a pathway for generating serine from the glycolytic intermediate 3-phosphoglycerate. Fluorescent tags confirmed the delivery of these serine biosynthetic enzymes to P. infestans mitochondria. A cytosolic form of this serine biosynthetic pathway, which occurs in most eukaryotes, is missing from oomycetes and most other stramenopiles. The glycolysis and serine metabolism pathways of oomycetes appear to be mosaics of enzymes with different ancestries. While some of the noncanonical oomycete mitochondrial enzymes have the closest affinity in phylogenetic analyses with proteins from other stramenopiles, others cluster with bacterial, plant, or animal proteins. The genes encoding the mitochondrial phosphoglycerate kinase and serine-forming enzymes are physically linked on oomycete chromosomes, which suggests a shared origin.

Conclusions

Stramenopile metabolism appears to have been shaped through the acquisition of genes by descent and lateral or endosymbiotic gene transfer, along with the targeting of the proteins to locations that are novel compared to other eukaryotes. Colocalization of the glycolytic and serine biosynthesis enzymes in mitochondria is apparently necessary since they share a common intermediate. The results indicate that descriptions of metabolism in textbooks do not cover the full diversity of eukaryotic biology.

Electronic supplementary material

The online version of this article (10.1186/s12862-017-1087-8) contains supplementary material, which is available to authorized users.

Metabolic Diversity and Novelties in the Oomycetes

The eukaryotic microbes called oomycetes include many important saprophytes and pathogens, with the latter exhibiting necrotrophy, biotrophy, or obligate biotrophy. Understanding oomycete metabolism is fundamental to understanding these lifestyles. Genome mining and biochemical studies have shown that oomycetes, which belong to the kingdom Stramenopila, secrete suites of carbohydrate- and protein-degrading enzymes adapted to their environmental niches and produce unusual lipids and energy storage compounds. Despite having limited secondary metabolism, many oomycetes make chemicals for communicating within their species or with their hosts. Horizontal and endosymbiotic gene transfer events have diversified oomycete metabolism, resulting in biochemical pathways that often depart from standard textbook descriptions by amalgamating enzymes from multiple sources. Gene fusions and duplications have further shaped the composition and expression of the enzymes. Current research is helping us learn how oomycetes interact with host and environment, understand eukaryotic diversity and evolution, and identify targets for drugs and crop protection chemicals.

Systematic Search for Evidence of Interdomain Horizontal Gene Transfer from Prokaryotes to Oomycete Lineages

While most commonly associated with prokaryotes, horizontal gene transfer (HGT) can also have a significant influence on the evolution of microscopic eukaryotes. Systematic analysis of HGT in the genomes of the oomycetes, filamentous eukaryotic microorganisms in the Stramenopiles-Alveolates-Rhizaria (SAR) supergroup, has to date focused mainly on intradomain transfer events between oomycetes and fungi. Using systematic whole-genome analysis followed by phylogenetic reconstruction, we have investigated the extent of interdomain HGT between bacteria and plant-pathogenic oomycetes. We report five putative instances of HGT from bacteria into the oomycetes. Two transfers were found in Phytophthora species, including one unique to the cucurbit pathogen Phytophthora capsici. Two were found in Pythium species only, and the final transfer event was present in Phytopythium and Pythium species, the first reported bacterium-inherited genes in these genera. Our putative transfers included one protein that appears to be a member of the Pythium secretome, metabolic proteins, and enzymes that could potentially break down xenobiotics within the cell. Our findings complement both previous reports of bacterial genes in oomycete and SAR genomes and the growing body of evidence suggesting that interdomain transfer from prokaryotes into eukaryotes occurs more frequently than previously thought. IMPORTANCE Horizontal gene transfer (HGT) is the nonvertical inheritance of genetic material by transfer between different species. HGT is an important evolutionary mechanism for prokaryotes and in some cases is responsible for the spread of antibiotic resistance from resistant to benign species. Genome analysis has shown that examples of HGT are not as frequent in eukaryotes, but when they do occur they may have important evolutionary consequences. For example, the acquisition of fungal genes by an ancestral Phytophthora (plant destroyer) species is responsible for the large repertoire of enzymes in the plant-degrading arsenal of modern-day Phytophthora species. In this analysis, we set out to systematically search oomycete genomes for evidence of interdomain HGT (transfer of bacterial genes into oomycete species). Our results show that interdomain HGT is rare in oomycetes but has occurred. We located five well-supported examples, including one that could potentially break down xenobiotics within the cell.

PsAAT3, an oomycete-specific aspartate aminotransferase, is required for full pathogenicity of the oomycete pathogen Phytophthora sojae

Pathogen nutrient acquisition and metabolism are critical for successful infection and colonization. However, the nutrient requirements and metabolic pathways related to pathogenesis in oomycete pathogens are unknown. In this study, we bioinformatically identified Phytophthora sojae aspartate aminotransferases (AATs), which are key enzymes that coordinate carbon and nitrogen metabolism. We demonstrated that P. sojae encodes more AATs than the analysed fungi. Some of the AATs contained additional prephenate dehydratase and/or prephenate dehydrogenase domains in their N-termini, which are unique to oomycetes. Silencing of PsAAT3, an infection-inducible expression gene, reduced P. sojae pathogenicity on soybean plants and affected the growth under N-starving condition, suggesting that PsAAT3 is involved in pathogen pathogenicity and nitrogen utilisation during infection. Our results suggest that P. sojae and other oomycete pathogens may have distinct amino acid metabolism pathways and that PsAAT3 is important for its full pathogenicity.

The secreted proteins of Achlya hypogyna and Thraustotheca clavata identify the ancestral oomycete secretome and reveal gene acquisitions by horizontal gene transfer

Saprotrophic and parasitic microorganisms secrete proteins into the environment to breakdown macromolecules and obtain nutrients. The molecules secreted are collectively termed the "secretome" and the composition and function of this set of proteins varies depending on the ecology, life cycle, and environment of an organism. Beyond the function of nutrient acquisition, parasitic lineages must also secrete molecules to manipulate their host. Here, we use a combination of de novo genome and transcriptome sequencing and bioinformatic identification of signal peptides to identify the putative secreted proteome of two oomycetes, the facultative parasite Achlya hypogyna and free-living Thraustotheca clavata. By comparing the secretomes of these saprolegnialean oomycetes with that of eight other oomycetes, we were able to characterize the evolution of this protein set across the oomycete clade. These species span the last common ancestor of the two major oomycete families allowing us to identify the ancestral secretome. This putative ancestral secretome consists of at least 84 gene families. Only 11 of these gene families are conserved across all 10 secretomes analyzed and the two major branches in the oomycete radiation. Notably, we have identified expressed elicitin-like effector genes in the saprotrophic decomposer, T. clavata. Phylogenetic analyses show six novel horizontal gene transfers to the oomycete secretome from bacterial and fungal donor lineages, four of which are specific to the Saprolegnialeans. Comparisons between free-living and pathogenic taxa highlight the functional changes of oomycete secretomes associated with shifts in lifestyle.

A multilocus timescale for oomycete evolution estimated under three distinct molecular clock models

BACKGROUND

Molecular clock methodologies allow for the estimation of divergence times across a variety of organisms; this can be particularly useful for groups lacking robust fossil histories, such as microbial eukaryotes with few distinguishing morphological traits. Here we have used a Bayesian molecular clock method under three distinct clock models to estimate divergence times within oomycetes, a group of fungal-like eukaryotes that are ubiquitous in the environment and include a number of devastating pathogenic species. The earliest fossil evidence for oomycetes comes from the Lower Devonian (~400 Ma), however the taxonomic affinities of these fossils are unclear.

RESULTS

Complete genome sequences were used to identify orthologous proteins among oomycetes, diatoms, and a brown alga, with a focus on conserved regulators of gene expression such as DNA and histone modifiers and transcription factors. Our molecular clock estimates place the origin of oomycetes by at least the mid-Paleozoic (~430-400 Ma), with the divergence between two major lineages, the peronosporaleans and saprolegnialeans, in the early Mesozoic (~225-190 Ma). Divergence times estimated under the three clock models were similar, although only the strict and random local clock models produced reliable estimates for most parameters.

CONCLUSIONS

Our molecular timescale suggests that modern pathogenic oomycetes diverged well after the origin of their respective hosts, indicating that environmental conditions or perhaps horizontal gene transfer events, rather than host availability, may have driven lineage diversification. Our findings also suggest that the last common ancestor of oomycetes possessed a full complement of eukaryotic regulatory proteins, including those involved in histone modification, RNA interference, and tRNA and rRNA methylation; interestingly no match to canonical DNA methyltransferases could be identified in the oomycete genomes studied here.

Horizontal transfer and evolution of prokaryote transposable elements in eukaryotes

Horizontal transfer (HT) of transposable elements (TEs) plays a key role in prokaryotic evolution, and mounting evidence suggests that it has also had an important impact on eukaryotic evolution. Although many prokaryote-to-prokaryote and eukaryote-to-eukaryote HTs of TEs have been characterized, only few cases have been reported between prokaryotes and eukaryotes. Here, we carried out a comprehensive search for all major groups of prokaryotic insertion sequences (ISs) in 430 eukaryote genomes. We uncovered a total of 80 sequences, all deriving from the IS607 family, integrated in the genomes of 14 eukaryote species belonging to four distinct phyla (Amoebozoa, Ascomycetes, Basidiomycetes, and Stramenopiles). Given that eukaryote IS607-like sequences are most closely related to cyanobacterial IS607 and that their phylogeny is incongruent with that of their hosts, we conclude that the presence of IS607-like sequences in eukaryotic genomes is the result of several HT events. Selection analyses further suggest that our ability to detect these prokaryote TEs today in eukaryotes is because HT of these sequences occurred recently and/or some IS607 elements were domesticated after HT, giving rise to new eukaryote genes. Supporting the recent age of some of these HTs, we uncovered intact full-length, potentially active IS607 copies in the amoeba Acanthamoeba castellani. Overall, our study shows that prokaryote-to-eukaryote HT of TEs occurred at relatively low frequency during recent eukaryote evolution and it sets IS607 as the most widespread TE (being present in prokaryotes, eukaryotes, and viruses).

Repeated loss of an anciently horizontally transferred gene cluster in Botrytis

At least five of the six genes of the bikaverin secondary metabolic gene cluster were shown to have undergone horizontal transfer (HGT) from a Fusarium donor to the Botrytis lineage. Of these five, two enzyme-encoding genes are found as pseudogenes in B. cinerea whereas two regulatory genes and the transporter remain intact. To reconstruct the evolutionary events leading to decay of this gene cluster and infer a more precise timing of its transfer, we examined the genomes of nine additional broadly sampled Botrytis species. We found evidence that a Botrytis ancestor acquired the entire gene cluster through an ancient HGT that occurred before the diversification of the genus. During the subsequent evolution and diversification of the genus, four of the 10 genomes appear to have lost the gene cluster, while in the other six the cluster is in various stages of degeneration. Across the Botrytis genomes, the modes of gene decay in the cluster differed between enzyme-encoding genes, which had higher rates of transition to or retention of pseudogenes and were universally inactivated, and regulatory genes (particularly the non-pathway-specific regulator bik4), which more frequently appeared intact. Consistent with these results, the regulatory genes bik4 and bik5 showed stronger evidence of transcriptional expression than other bikaverin genes under multiple conditions in B. cinerea. These results could be explained by pleiotropy in the bikaverin regulatory genes either through rewiring or their interaction with more central pathways or by constraints on the order of gene loss driven by the intrinsic toxicity of the pathway. Our finding that most of the bikaverin pathway genes have been lost or pseudogenized in these Botrytis genomes suggests that the incidence of HGT of gene cluster-encoded metabolic pathways might be higher than what is possible to be inferred from isolated genome analyses.

Distinctive expansion of potential virulence genes in the genome of the oomycete fish pathogen Saprolegnia parasitica

Oomycetes in the class Saprolegniomycetidae of the Eukaryotic kingdom Stramenopila have evolved as severe pathogens of amphibians, crustaceans, fish and insects, resulting in major losses in aquaculture and damage to aquatic ecosystems. We have sequenced the 63 Mb genome of the fresh water fish pathogen, Saprolegnia parasitica. Approximately 1/3 of the assembled genome exhibits loss of heterozygosity, indicating an efficient mechanism for revealing new variation. Comparison of S. parasitica with plant pathogenic oomycetes suggests that during evolution the host cellular environment has driven distinct patterns of gene expansion and loss in the genomes of plant and animal pathogens. S. parasitica possesses one of the largest repertoires of proteases (270) among eukaryotes that are deployed in waves at different points during infection as determined from RNA-Seq data. In contrast, despite being capable of living saprotrophically, parasitism has led to loss of inorganic nitrogen and sulfur assimilation pathways, strikingly similar to losses in obligate plant pathogenic oomycetes and fungi. The large gene families that are hallmarks of plant pathogenic oomycetes such as Phytophthora appear to be lacking in S. parasitica, including those encoding RXLR effectors, Crinkler's, and Necrosis Inducing-Like Proteins (NLP). S. parasitica also has a very large kinome of 543 kinases, 10% of which is induced upon infection. Moreover, S. parasitica encodes several genes typical of animals or animal-pathogens and lacking from other oomycetes, including disintegrins and galactose-binding lectins, whose expression and evolutionary origins implicate horizontal gene transfer in the evolution of animal pathogenesis in S. parasitica.

Fish are an increasingly important source of animal protein globally, with aquaculture production rising dramatically over the past decade. Saprolegnia is a fungal-like oomycete and one of the most destructive fish pathogens, causing millions of dollars in losses to the aquaculture industry annually. Saprolegnia has also been linked to a worldwide decline in wild fish and amphibian populations. Here we describe the genome sequence of the first animal pathogenic oomycete and compare the genome content with the available plant pathogenic oomycetes. We found that Saprolegnia lacks the large effector families that are hallmarks of plant pathogenic oomycetes, showing evolutionary adaptation to the host. Moreover, Saprolegnia harbors pathogenesis-related genes that were derived by lateral gene transfer from the host and other animal pathogens. The retrotransposon LINE family also appears to be acquired from animal lineages. By transcriptome analysis we show a high rate of allelic variation, which reveals rapidly evolving genes and potentially adaptive evolutionary mechanisms coupled to selective pressures exerted by the animal host. The genome and transcriptome data, as well as subsequent biochemical analyses, provided us with insight in the disease process of Saprolegnia at a molecular and cellular level, providing us with targets for sustainable control of Saprolegnia.

Decay of genes encoding the oomycete flagellar proteome in the downy mildew Hyaloperonospora arabidopsidis

Zoospores are central to the life cycles of most of the eukaryotic microbes known as oomycetes, but some genera have lost the ability to form these flagellated cells. In the plant pathogen Phytophthora infestans, genes encoding 257 proteins associated with flagella were identified by comparative genomics. These included the main structural components of the axoneme and basal body, proteins involved in intraflagellar transport, regulatory proteins, enzymes for maintaining ATP levels, and others. Transcripts for over three-quarters of the genes were up-regulated during sporulation, and persisted to varying degrees in the pre-zoospore stage (sporangia) and motile zoospores. Nearly all of these genes had orthologs in other eukaryotes that form flagella or cilia, but not species that lack the organelle. Orthologs of 211 of the genes were also absent from a sister taxon to P. infestans that lost the ability to form flagella, the downy mildew Hyaloperonospora arabidopsidis. Many of the genes retained in H. arabidopsidis were also present in other non-flagellates, suggesting that they play roles both in flagella and other cellular processes. Remnants of the missing genes were often detected in the H. arabidopsidis genome. Degradation of the genes was associated with local compaction of the chromosome and a heightened propensity towards genome rearrangements, as such regions were less likely to share synteny with P. infestans.

Dynamics and innovations within oomycete genomes: insights into biology, pathology, and evolution

The eukaryotic microbes known as oomycetes are common inhabitants of terrestrial and aquatic environments and include saprophytes and pathogens. Lifestyles of the pathogens extend from biotrophy to necrotrophy, obligate to facultative pathogenesis, and narrow to broad host ranges on plants or animals. Sequencing of several pathogens has revealed striking variation in genome size and content, a plastic set of genes related to pathogenesis, and adaptations associated with obligate biotrophy. Features of genome evolution include repeat-driven expansions, deletions, gene fusions, and horizontal gene transfer in a landscape organized into gene-dense and gene-sparse sectors and influenced by transposable elements. Gene expression profiles are also highly dynamic throughout oomycete life cycles, with transcriptional polymorphisms as well as differences in protein sequence contributing to variation. The genome projects have set the foundation for functional studies and should spur the sequencing of additional species, including more diverse pathogens and nonpathogens.

Reconstruction of oomycete genome evolution identifies differences in evolutionary trajectories leading to present-day large gene families

The taxonomic class of oomycetes contains numerous pathogens of plants and animals but is related to nonpathogenic diatoms and brown algae. Oomycetes have flexible genomes comprising large gene families that play roles in pathogenicity. The evolutionary processes that shaped the gene content have not yet been studied by applying systematic tree reconciliation of the phylome of these species. We analyzed evolutionary dynamics of ten Stramenopiles. Gene gains, duplications, and losses were inferred by tree reconciliation of 18,459 gene trees constituting the phylome with a highly supported species phylogeny. We reconstructed a strikingly large last common ancestor of the Stramenopiles that contained ∼10,000 genes. Throughout evolution, the genomes of pathogenic oomycetes have constantly gained and lost genes, though gene gains through duplications outnumber the losses. The branch leading to the plant pathogenic Phytophthora genus was identified as a major transition point characterized by increased frequency of duplication events that has likely driven the speciation within this genus. Large gene families encoding different classes of enzymes associated with pathogenicity such as glycoside hydrolases are formed by complex and distinct patterns of duplications and losses leading to their expansion in extant oomycetes. This study unveils the large-scale evolutionary dynamics that shaped the genomes of pathogenic oomycetes. By the application of phylogenetic based analyses methods, it provides additional insights that shed light on the complex history of oomycete genome evolution and the emergence of large gene families characteristic for this important class of pathogens.

Mechanisms and evolution of virulence in oomycetes

Many destructive diseases of plants and animals are caused by oomycetes, a group of eukaryotic pathogens important to agricultural, ornamental, and natural ecosystems. Understanding the mechanisms underlying oomycete virulence and the genomic processes by which those mechanisms rapidly evolve is essential to developing effective long-term control measures for oomycete diseases. Several common mechanisms underlying oomycete virulence, including protein toxins and cell-entering effectors, have emerged from comparing oomycetes with different genome characteristics, parasitic lifestyles, and host ranges. Oomycete genomes display a strongly bipartite organization in which conserved housekeeping genes are concentrated in syntenic gene-rich blocks, whereas virulence genes are dispersed into highly dynamic, repeat-rich regions. There is also evidence that key virulence genes have been acquired by horizontal transfer from other eukaryotic and prokaryotic species.

The aspartic proteinase family of three Phytophthora species

BACKGROUND

Phytophthora species are oomycete plant pathogens with such major social and economic impact that genome sequences have been determined for Phytophthora infestans, P. sojae and P. ramorum. Pepsin-like aspartic proteinases (APs) are produced in a wide variety of species (from bacteria to humans) and contain conserved motifs and landmark residues. APs fulfil critical roles in infectious organisms and their host cells. Annotation of Phytophthora APs would provide invaluable information for studies into their roles in the physiology of Phytophthora species and interactions with their hosts.

RESULTS

Genomes of Phytophthora infestans, P. sojae and P. ramorum contain 11-12 genes encoding APs. Nine of the original gene models in the P. infestans database and several in P. sojae and P. ramorum (three and four, respectively) were erroneous. Gene models were corrected on the basis of EST data, consistent positioning of introns between orthologues and conservation of hallmark motifs. Phylogenetic analysis resolved the Phytophthora APs into 5 clades. Of the 12 sub-families, several contained an unconventional architecture, as they either lacked a signal peptide or a propart region. Remarkably, almost all APs are predicted to be membrane-bound.

CONCLUSIONS

One of the twelve Phytophthora APs is an unprecedented fusion protein with a putative G-protein coupled receptor as the C-terminal partner. The others appear to be related to well-documented enzymes from other species, including a vacuolar enzyme that is encoded in every fungal genome sequenced to date. Unexpectedly, however, the oomycetes were found to have both active and probably-inactive forms of an AP similar to vertebrate BACE, the enzyme responsible for initiating the processing cascade that generates the Aβ peptide central to Alzheimer's Disease. The oomycetes also encode enzymes similar to plasmepsin V, a membrane-bound AP that cleaves effector proteins of the malaria parasite Plasmodium falciparum during their translocation into the host red blood cell. Since the translocation of Phytophthora effector proteins is currently a topic of intense research activity, the identification in Phytophthora of potential functional homologues of plasmepsin V would appear worthy of investigation. Indeed, elucidation of the physiological roles of the APs identified here offers areas for future study. The significant revision of gene models and detailed annotation presented here should significantly facilitate experimental design.

A domain-centric analysis of oomycete plant pathogen genomes reveals unique protein organization

Oomycetes comprise a diverse group of organisms that morphologically resemble fungi but belong to the stramenopile lineage within the supergroup of chromalveolates. Recent studies have shown that plant pathogenic oomycetes have expanded gene families that are possibly linked to their pathogenic lifestyle. We analyzed the protein domain organization of 67 eukaryotic species including four oomycete and five fungal plant pathogens. We detected 246 expanded domains in fungal and oomycete plant pathogens. The analysis of genes differentially expressed during infection revealed a significant enrichment of genes encoding expanded domains as well as signal peptides linking a substantial part of these genes to pathogenicity. Overrepresentation and clustering of domain abundance profiles revealed domains that might have important roles in host-pathogen interactions but, as yet, have not been linked to pathogenicity. The number of distinct domain combinations (bigrams) in oomycetes was significantly higher than in fungi. We identified 773 oomycete-specific bigrams, with the majority composed of domains common to eukaryotes. The analyses enabled us to link domain content to biological processes such as host-pathogen interaction, nutrient uptake, or suppression and elicitation of plant immune responses. Taken together, this study represents a comprehensive overview of the domain repertoire of fungal and oomycete plant pathogens and points to novel features like domain expansion and species-specific bigram types that could, at least partially, explain why oomycetes are such remarkable plant pathogens.

The kinome of Phytophthora infestans reveals oomycete-specific innovations and links to other taxonomic groups

BACKGROUND

Oomycetes are a large group of economically and ecologically important species. Its most notorious member is Phytophthora infestans, the cause of the devastating potato late blight disease. The life cycle of P. infestans involves hyphae which differentiate into spores used for dispersal and host infection. Protein phosphorylation likely plays crucial roles in these stages, and to help understand this we present here a genome-wide analysis of the protein kinases of P. infestans and several relatives. The study also provides new insight into kinase evolution since oomycetes are taxonomically distant from organisms with well-characterized kinomes.

RESULTS

Bioinformatic searches of the genomes of P. infestans, P. ramorum, and P. sojae reveal they have similar kinomes, which for P. infestans contains 354 eukaryotic protein kinases (ePKs) and 18 atypical kinases (aPKs), equaling 2% of total genes. After refining gene models, most were classifiable into families seen in other eukaryotes. Some ePK families are nevertheless unusual, especially the tyrosine kinase-like (TKL) group which includes large oomycete-specific subfamilies. Also identified were two tyrosine kinases, which are rare in non-metazoans. Several ePKs bear accessory domains not identified previously on kinases, such as cyclin-dependent kinases with integral cyclin domains. Most ePKs lack accessory domains, implying that many are regulated transcriptionally. This was confirmed by mRNA expression-profiling studies that showed that two-thirds vary significantly between hyphae, sporangia, and zoospores. Comparisons to neighboring taxa (apicomplexans, ciliates, diatoms) revealed both clade-specific and conserved features, and multiple connections to plant kinases were observed. The kinome of Hyaloperonospora arabidopsidis, an oomycete with a simpler life cycle than P. infestans, was found to be one-third smaller. Some differences may be attributable to gene clustering, which facilitates subfamily expansion (or loss) through unequal crossing-over.

CONCLUSION

The large sizes of the Phytophthora kinomes imply that phosphorylation plays major roles in their life cycles. Their kinomes also include many novel ePKs, some specific to oomycetes or shared with neighboring groups. Little experimentation to date has addressed the biological functions of oomycete kinases, but this should be stimulated by the structural, evolutionary, and expression data presented here. This may lead to targets for disease control.

Runaway repeats force expansion of the Phytophthora infestans genome

Sequencing of the potato late blight pathogen Phytophthora infestans provides insight into the structure and evolution of its genome.

Sequencing of the genome of the potato late blight pathogen Phytophthora infestans provides insight into genome structure and evolution within this genus of plant pathogenic oomycetes.