Comparative analysis of secreted protein evolution using expressed sequence tags from four poplar leaf rusts (Melampsora spp.)

The Fungal Effector Mlp37347 Alters Plasmodesmata Fluxes and Enhances Susceptibility to Pathogen

Melampsora larici-populina (Mlp) is a devastating pathogen of poplar trees, causing the defoliating poplar leaf rust disease. Genomic studies have revealed that Mlp possesses a repertoire of 1184 small secreted proteins (SSPs), some of them being characterized as candidate effectors. However, how they promote virulence is still unclear. This study investigates the candidate effector Mlp37347’s role during infection. We developed a stable Arabidopsis transgenic line expressing Mlp37347 tagged with the green fluorescent protein (GFP). We found that the effector accumulated exclusively at plasmodesmata (PD). Moreover, the presence of the effector at plasmodesmata favors enhanced plasmodesmatal flux and reduced callose deposition. Transcriptome profiling and a gene ontology (GO) analysis of transgenic Arabidopsis plants expressing the effector revealed that the genes involved in glucan catabolic processes are up-regulated. This effector has previously been shown to interact with glutamate decarboxylase 1 (GAD1), and in silico docking analysis supported the strong binding between Mlp37347 and GAD1 in this study. In infection assays, the effector promoted Hyalonoperospora arabidopsidis growth but not bacterial growth. Our investigation suggests that the effector Mlp37347 targets PD in host cells and promotes parasitic growth.

Antivirulence and avirulence genes in human pathogenic fungi

Opportunistic commensal and environmental fungi can cause superficial to systemic diseases in humans. But how did these pathogens adapt to infect us and how does host-pathogen co-evolution shape their virulence potential? During evolution toward pathogenicity, not only do microorganisms gain virulence genes, but they also tend to lose non-adaptive genes in the host niche. Additionally, virulence factors can become detrimental during infection when they trigger host recognition. The loss of non-adaptive genes as well as the loss of the virulence potential of genes by adaptations to the host has been investigated in pathogenic bacteria and phytopathogenic fungi, where they are known as antivirulence and avirulence genes, respectively. However, these concepts are nearly unknown in the field of pathogenic fungi of humans. We think that this unnecessarily limits our view of human-fungal interplay, and that much could be learned if we applied a similar framework to aspects of these interactions. In this review, we, therefore, define and adapt the concepts of antivirulence and avirulence genes for human pathogenic fungi. We provide examples for analogies to antivirulence genes of bacterial pathogens and to avirulence genes of phytopathogenic fungi. Introducing these terms to the field of pathogenic fungi of humans can help to better comprehend the emergence and evolution of fungal virulence and disease.

Identification of Genes Underlying the Resistance to Melampsora larici-populina in an R Gene Supercluster of the Populus deltoides Genome

Identification of the particular genes in an R genes supercluster underlying resistance to the rust fungus Melampsora larici-populina in poplar genome remains challenging. Based on the de novo assembly of the Populus deltoides genome, all of the detected major genetic loci conferring resistance to M. larici-populina were confined to a 3.5-Mb region on chromosome 19. The transcriptomes of the resistant and susceptible genotypes were sequenced for a timespan from 0 to 168 hours postinoculation. By mapping the differentially expressed genes to the target genomic region, we identified two constitutive expression R genes and one inducible expression R gene that might confer resistance to M. larici-populina. Nucleotide variations were predicted based on the reconstructed haplotypes for each allele of the candidate genes. We also confirmed that salicylic acid was the phytohormone mediating signal transduction pathways, and PR-1 was identified as a key gene inhibiting rust reproduction. Finally, quantitative reverse transcription PCR assay revealed consistent expressions with the RNA-sequencing data for the detected key genes. This study presents an efficient approach for the identification of particular genes underlying phenotype of interest by the combination of genetic mapping, transcriptome profiling, and candidate gene sequences dissection. The identified key genes would be useful for host resistance diagnosis and for molecular breeding of elite poplar cultivars exhibiting resistance to M. larici-populina infection. The detected R genes are also valuable for testing whether the combination of individual R genes can induce durable quantitative resistance.

Comparative genomic and transcriptomic analyses reveal different pathogenicity-related genes among three eucalyptus fungal pathogens

Ceratocystis fimbriata is an important plant pathogen known to cause Ceratocystis Wilt (CW), a prevalent fungal disease known to affect Eucalyptus spp. plantations in Brazil. To better understand the molecular mechanisms related to pathogenicity in eucalyptus, we generated a high-quality assembly and annotation of the Ce. fimbriata LPF1912 isolate (LPF1912) genome, as well as the first transcriptome of LPF1912 from 16 eucalyptus clones at three infection incubation periods (12, 18, and 24 h). The LPF1912 genome assembly contains 805 scaffolds, totaling 31.8 Mb, with 43% of the genome estimated to be coding sequence comprised of 7,390 protein-coding genes of which 626 (8.5%) were classified as secreted proteins, 120 ribosomal RNAs, and 532 transfer RNAs. Comparative genomic analysis among three eucalyptus fungal pathogens (Ce. fimbriata, Ce. eucalypticola, and Calonectria pseudoreteaudii), showed high similarity in the proteome (21.81%) and secretome (52.01%) of LPF1912 and Ce. eucalypticola. GO annotation of pathogenicity-related genes of LPF1912 and Ce. eucalypticola, revealed enrichment in cell wall degrading enzymes (CWDEs), and lipid/cutin metabolism for Ca. pseudoreteaudii. Additionally, a transcriptome analysis between resistant and susceptible eucalyptus clones to CW infection indicated that a majority (11) of LPF1912 differentially expressed genes had GO terms associated with enzymatic functions, such as the polygalacturonase gene family, confirming the crucial role of CWDEs for Ce. fimbriata pathogenicity. Finally, our genomic and transcriptomic analysis approach provides a better understanding of the mechanisms involved in Ce. fimbriata pathogenesis, as well as a framework for further studies.

Mitotic Recombination and Rapid Genome Evolution in the Invasive Forest Pathogen Phytophthora ramorum

Alien species are often successful invaders in new environments, despite the introduction of a few isolates with a reduced genetic pool. This is called the genetic paradox of invasion. We found two mechanisms by which the invasive forest pathogen causing sudden oak and sudden larch death can evolve. Extensive mitotic recombination producing runs of homozygosity generates genotypic diversity even in the absence of sexual reproduction, and rapid turnover of genes in the non-core, or nonessential portion of genome not shared by all isolates, allows pathogenicity genes to evolve rapidly or be eliminated while retaining essential genes. Mitotic recombination events occur in genomic hot spots, resulting in similar ROH patterns in different isolates or groups; one ROH, independently generated in two different groups, was enriched in pathogenicity genes and may be a target for selection. This provides important insights into the evolution of invasive alien pathogens and their potential for adaptation and future persistence.

Invasive alien species often have reduced genetic diversity and must adapt to new environments. Given the success of many invasions, this is sometimes called the genetic paradox of invasion. Phytophthora ramorum is invasive, limited to asexual reproduction within four lineages, and presumed clonal. It is responsible for sudden oak death in the United States, sudden larch death in Europe, and ramorum blight in North America and Europe. We sequenced the genomes of 107 isolates to determine how this pathogen can overcome the invasion paradox. Mitotic recombination (MR) associated with transposons and low gene density has generated runs of homozygosity (ROH) affecting 2,698 genes, resulting in novel genotypic diversity within the lineages. One ROH enriched in effectors was fixed in the NA1 lineage. An independent ROH affected the same scaffold in the EU1 lineage, suggesting an MR hot spot and a selection target. Differences in host infection between EU1 isolates with and without the ROH suggest that they may differ in aggressiveness. Non-core regions (not shared by all lineages) had signatures of accelerated evolution and were enriched in putative pathogenicity genes and transposons. There was a striking pattern of gene loss, including all effectors, in the non-core EU2 genome. Positive selection was observed in 8.0% of RxLR and 18.8% of Crinkler effector genes compared with 0.9% of the core eukaryotic gene set. We conclude that the P. ramorum lineages are diverging via a rapidly evolving non-core genome and that the invasive asexual lineages are not clonal, but display genotypic diversity caused by MR.

The Rust Fungus Melampsora larici-populina Expresses a Conserved Genetic Program and Distinct Sets of Secreted Protein Genes During Infection of Its Two Host Plants, Larch and Poplar

Mechanisms required for broad-spectrum or specific host colonization of plant parasites are poorly understood. As a perfect illustration, heteroecious rust fungi require two alternate host plants to complete their life cycles. Melampsora larici-populina infects two taxonomically unrelated plants, larch, on which sexual reproduction is achieved, and poplar, on which clonal multiplication occurs, leading to severe epidemics in plantations. We applied deep RNA sequencing to three key developmental stages of M. larici-populina infection on larch: basidia, pycnia, and aecia, and we performed comparative transcriptomics of infection on poplar and larch hosts, using available expression data. Secreted protein was the only significantly overrepresented category among differentially expressed M. larici-populina genes between the basidial, the pycnial, and the aecial stages, highlighting their probable involvement in the infection process. Comparison of fungal transcriptomes in larch and poplar revealed a majority of rust genes were commonly expressed on the two hosts and a fraction exhibited host-specific expression. More particularly, gene families encoding small secreted proteins presented striking expression profiles that highlight probable candidate effectors specialized on each host. Our results bring valuable new information about the biological cycle of rust fungi and identify genes that may contribute to host specificity.

Positively Selected Effector Genes and Their Contribution to Virulence in the Smut Fungus Sporisorium reilianum

Plants and fungi display a broad range of interactions in natural and agricultural ecosystems ranging from symbiosis to parasitism. These ecological interactions result in coevolution between genes belonging to different partners. A well-understood example is secreted fungal effector proteins and their host targets, which play an important role in pathogenic interactions. Biotrophic smut fungi (Basidiomycota) are well-suited to investigate the evolution of plant pathogens, because several reference genomes and genetic tools are available for these species. Here, we used the genomes of Sporisorium reilianum f. sp. zeae and S. reilianum f. sp. reilianum, two closely related formae speciales infecting maize and sorghum, respectively, together with the genomes of Ustilago hordei, Ustilago maydis, and Sporisorium scitamineum to identify and characterize genes displaying signatures of positive selection. We identified 154 gene families having undergone positive selection during species divergence in at least one lineage, among which 77% were identified in the two investigated formae speciales of S. reilianum. Remarkably, only 29% of positively selected genes encode predicted secreted proteins. We assessed the contribution to virulence of nine of these candidate effector genes in S. reilianum f. sp. zeae by deleting individual genes, including a homologue of the effector gene pit2 previously characterized in U. maydis. Only the pit2 deletion mutant was found to be strongly reduced in virulence. Additional experiments are required to understand the molecular mechanisms underlying the selection forces acting on the other candidate effector genes, as well as the large fraction of positively selected genes encoding predicted cytoplasmic proteins.

Infection assays in Arabidopsis reveal candidate effectors from the poplar rust fungus that promote susceptibility to bacteria and oomycete pathogens

Fungi of the Pucciniales order cause rust diseases which, altogether, affect thousands of plant species worldwide and pose a major threat to several crops. How rust effectors-virulence proteins delivered into infected tissues to modulate host functions-contribute to pathogen virulence remains poorly understood. Melampsora larici-populina is a devastating and widespread rust pathogen of poplar, and its genome encodes 1184 identified small secreted proteins that could potentially act as effectors. Here, following specific criteria, we selected 16 candidate effector proteins and characterized their virulence activities and subcellular localizations in the leaf cells of Arabidopsis thaliana. Infection assays using bacterial (Pseudomonas syringae) and oomycete (Hyaloperonospora arabidopsidis) pathogens revealed subsets of candidate effectors that enhanced or decreased pathogen leaf colonization. Confocal imaging of green fluorescent protein-tagged candidate effectors constitutively expressed in stable transgenic plants revealed that some protein fusions specifically accumulate in nuclei, chloroplasts, plasmodesmata and punctate cytosolic structures. Altogether, our analysis suggests that rust fungal candidate effectors target distinct cellular components in host cells to promote parasitic growth.

Evolutionary dynamics of interactions between plants and their enemies: comparison of herbivorous insects and pathogens

Plants colonized land over 400 million years ago. Shortly thereafter, organisms began to consume terrestrial plant tissue as a nutritional resource. Most plant enemies are plant pathogens or herbivores, and they impose natural selection for plants to evolve defenses. These traits generate selection pressures on enemies. Coevolution between terrestrial plants and their enemies is an important element of the evolutionary history of both groups. However, coevolutionary studies of plant-pathogen interactions have tended to focus on different research topics than plant-herbivore interactions. Specifically, studies of plant-pathogen interactions often adopt a "gene-for-gene" conceptual framework. In contrast, studies of plants and herbivores often investigate escalation or elaboration of plant defense and herbivore adaptations to overcome it. The main exceptions to the general pattern are studies that focus on small, sessile herbivores that share many features with plant pathogens, studies that incorporate both herbivores and pathogens into a single investigation, and studies that test aspects of Thompson's geographic mosaic theory for coevolution. We discuss the implications of these findings for future research.

The prediction of a pathogenesis-related secretome of Puccinia helianthi through high-throughput transcriptome analysis

Background

Many plant pathogen secretory proteins are known to be elicitors or pathogenic factors,which play an important role in the host-pathogen interaction process. Bioinformatics approaches make possible the large scale prediction and analysis of secretory proteins from the Puccinia helianthi transcriptome. The internet-based software SignalP v4.1, TargetP v1.01, Big-PI predictor, TMHMM v2.0 and ProtComp v9.0 were utilized to predict the signal peptides and the signal peptide-dependent secreted proteins among the 35,286 ORFs of the P. helianthi transcriptome.

Results

908 ORFs (accounting for 2.6% of the total proteins) were identified as putative secretory proteins containing signal peptides. The length of the majority of proteins ranged from 51 to 300 amino acids (aa), while the signal peptides were from 18 to 20 aa long. Signal peptidase I (SpI) cleavage sites were found in 463 of these putative secretory signal peptides. 55 proteins contained the lipoprotein signal peptide recognition site of signal peptidase II (SpII). Out of 908 secretory proteins, 581 (63.8%) have functions related to signal recognition and transduction, metabolism, transport and catabolism. Additionally, 143 putative secretory proteins were categorized into 27 functional groups based on Gene Ontology terms, including 14 groups in biological process, seven in cellular component, and six in molecular function. Gene ontology analysis of the secretory proteins revealed an enrichment of hydrolase activity. Pathway associations were established for 82 (9.0%) secretory proteins. A number of cell wall degrading enzymes and three homologous proteins specific to Phytophthora sojae effectors were also identified, which may be involved in the pathogenicity of the sunflower rust pathogen.

Conclusions

This investigation proposes a new approach for identifying elicitors and pathogenic factors. The eventual identification and characterization of 908 extracellularly secreted proteins will advance our understanding of the molecular mechanisms of interactions between sunflower and rust pathogen and will enhance our ability to intervene in disease states.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1577-0) contains supplementary material, which is available to authorized users.

The Hemileia vastatrix effector HvEC-016 suppresses bacterial blight symptoms in coffee genotypes with the SH 1 rust resistance gene

A number of genes that confer resistance to coffee leaf rust (S_H 1-S_H 9) have been identified within the genus Coffea, but despite many years of research on this pathosystem, the complementary avirulence genes of Hemileia vastatrix have not been reported. After identification of H. vastatrix effector candidate genes (HvECs) expressed at different stages of its lifecycle, we established an assay to characterize HvEC proteins by delivering them into coffee cells via the type-three secretion system (T3SS) of Pseudomonas syringae pv. garcae (Psgc). Employing a calmodulin-dependent adenylate cyclase assay, we demonstrate that Psgc recognizes a heterologous P. syringae T3SS secretion signal which enables us to translocate HvECs into the cytoplasm of coffee cells. Using this Psgc-adapted effector detector vector (EDV) system, we found that HvEC-016 suppresses the growth of Psgc on coffee genotypes with the S_H 1 resistance gene. Suppression of bacterial blight symptoms in S_H 1 plants was associated with reduced bacterial multiplication. By contrast, HvEC-016 enhanced bacterial multiplication in S_H 1-lacking plants. Our findings suggest that HvEC-016 may be recognized by the plant immune system in a S_H 1-dependent manner. Thus, our experimental approach is an effective tool for the characterization of effector/avirulence proteins of this important pathogen.

Effector-Mining in the Poplar Rust Fungus Melampsora larici-populina Secretome

The poplar leaf rust fungus, Melampsora larici-populina has been established as a tree-microbe interaction model. Understanding the molecular mechanisms controlling infection by pathogens appears essential for durable management of tree plantations. In biotrophic plant-parasites, effectors are known to condition host cell colonization. Thus, investigation of candidate secreted effector proteins (CSEPs) is a major goal in the poplar-poplar rust interaction. Unlike oomycetes, fungal effectors do not share conserved motifs and candidate prediction relies on a set of a priori criteria established from reported bona fide effectors. Secretome prediction, genome-wide analysis of gene families and transcriptomics of M. larici-populina have led to catalogs of more than a thousand secreted proteins. Automatized effector-mining pipelines hold great promise for rapid and systematic identification and prioritization of CSEPs for functional characterization. In this review, we report on and discuss the current status of the poplar rust fungus secretome and prediction of candidate effectors from this species.

Candidate Effector Proteins of the Rust Pathogen Melampsora larici-populina Target Diverse Plant Cell Compartments

Rust fungi are devastating crop pathogens that deliver effector proteins into infected tissues to modulate plant functions and promote parasitic growth. The genome of the poplar leaf rust fungus Melampsora larici-populina revealed a large catalog of secreted proteins, some of which have been considered candidate effectors. Unraveling how these proteins function in host cells is a key to understanding pathogenicity mechanisms and developing resistant plants. In this study, we used an effectoromics pipeline to select, clone, and express 20 candidate effectors in Nicotiana benthamiana leaf cells to determine their subcellular localization and identify the plant proteins they interact with. Confocal microscopy revealed that six candidate effectors target the nucleus, nucleoli, chloroplasts, mitochondria, and discrete cellular bodies. We also used coimmunoprecipitation (coIP) and mass spectrometry to identify 606 N. benthamiana proteins that associate with the candidate effectors. Five candidate effectors specifically associated with a small set of plant proteins that may represent biologically relevant interactors. We confirmed the interaction between the candidate effector MLP124017 and TOPLESS-related protein 4 from poplar by in planta coIP. Altogether, our data enable us to validate effector proteins from M. larici-populina and reveal that these proteins may target multiple compartments and processes in plant cells. It also shows that N. benthamiana can be a powerful heterologous system to study effectors of obligate biotrophic pathogens.

The effect on the transcriptome of Anemone coronaria following infection with rust (Tranzschelia discolor)

In order to understand plant/pathogen interaction, the transcriptome of uninfected (1S) and infected (2I) plant was sequenced at 3'end by the GS FLX 454 platform. De novo assembly of high-quality reads generated 27,231 contigs leaving 37,191 singletons in the 1S and 38,393 in the 2I libraries. ESTcalc tool suggested that 71% of the transcriptome had been captured, with 99% of the genes present being represented by at least one read. Unigene annotation showed that 50.5% of the predicted translation products shared significant homology with protein sequences in GenBank. In all 253 differential transcript abundance (DTAs) were in higher abundance and 52 in lower abundance in the 2I library. 128 higher abundance DTA genes were of fungal origin and 49 were clearly plant sequences. A tBLASTn-based search of the sequences using as query the full length predicted polypeptide product of 50 R genes identified 16 R gene products. Only one R gene (PGIP) was up-regulated. The response of the plant to fungal invasion included the up-regulation of several pathogenesis related protein (PR) genes involved in JA signaling and other genes associated with defense response and down regulation of cell wall associated genes, non-race-specific disease resistance1 (NDR1) and other genes like myb, presqualene diphosphate phosphatase (PSDPase), a UDP-glycosyltransferase 74E2-like (UGT). The DTA genes identified here should provide a basis for understanding the A. coronaria/T. discolor interaction and leads for biotechnology-based disease resistance breeding.

Fungal effectors and plant susceptibility

Plants can be colonized by fungi that have adopted highly diverse lifestyles, ranging from symbiotic to necrotrophic. Colonization is governed in all systems by hundreds of secreted fungal effector molecules. These effectors suppress plant defense responses and modulate plant physiology to accommodate fungal invaders and provide them with nutrients. Fungal effectors either function in the interaction zone between the fungal hyphae and host or are transferred to plant cells. This review describes the effector repertoires of 84 plant-colonizing fungi. We focus on the mechanisms that allow these fungal effectors to promote virulence or compatibility, discuss common plant nodes that are targeted by effectors, and provide recent insights into effector evolution. In addition, we address the issue of effector uptake in plant cells and highlight open questions and future challenges.

High-resolution transcript profiling of the atypical biotrophic interaction between Theobroma cacao and the fungal pathogen Moniliophthora perniciosa

Witches' broom disease (WBD), caused by the hemibiotrophic fungus Moniliophthora perniciosa, is one of the most devastating diseases of Theobroma cacao, the chocolate tree. In contrast to other hemibiotrophic interactions, the WBD biotrophic stage lasts for months and is responsible for the most distinctive symptoms of the disease, which comprise drastic morphological changes in the infected shoots. Here, we used the dual RNA-seq approach to simultaneously assess the transcriptomes of cacao and M. perniciosa during their peculiar biotrophic interaction. Infection with M. perniciosa triggers massive metabolic reprogramming in the diseased tissues. Although apparently vigorous, the infected shoots are energetically expensive structures characterized by the induction of ineffective defense responses and by a clear carbon deprivation signature. Remarkably, the infection culminates in the establishment of a senescence process in the host, which signals the end of the WBD biotrophic stage. We analyzed the pathogen's transcriptome in unprecedented detail and thereby characterized the fungal nutritional and infection strategies during WBD and identified putative virulence effectors. Interestingly, M. perniciosa biotrophic mycelia develop as long-term parasites that orchestrate changes in plant metabolism to increase the availability of soluble nutrients before plant death. Collectively, our results provide unique insight into an intriguing tropical disease and advance our understanding of the development of (hemi)biotrophic plant-pathogen interactions.

Compositions of fungal secretomes indicate a greater impact of phylogenetic history than lifestyle adaptation

Background

Since the first fungal genome sequences became available, investigators have been employing comparative genomics to understand how fungi have evolved to occupy diverse ecological niches. The secretome, i.e. the entirety of all proteins secreted by an organism, is of particular importance, as by these proteins fungi acquire nutrients and communicate with their surroundings.

Results

It is generally assumed that fungi with similar nutritional lifestyles have similar secretome compositions. In this study, we test this hypothesis by annotating and comparing the soluble secretomes, defined as the sets of proteins containing classical signal peptides but lacking transmembrane domains of fungi representing a broad diversity of nutritional lifestyles. Secretome size correlates with phylogeny and to a lesser extent with lifestyle. Plant pathogens and saprophytes have larger secretomes than animal pathogens. Small secreted cysteine-rich proteins (SSCPs), which may comprise many effectors important for the interaction of plant pathogens with their hosts, are defined here to have a mature length of ≤ 300 aa residues, at least four cysteines, and a total cysteine content of ≥5%. SSCPs are found enriched in the secretomes of the Pezizomycotina and Basidiomycota in comparison to Saccharomycotina. Relative SSCP content is noticeably higher in plant pathogens than in animal pathogens, while saprophytes were in between and closer to plant pathogens. Expansions and contractions of gene families and in the number of occurrences of functional domains are largely lineage specific, e.g. contraction of glycoside hydrolases in Saccharomycotina, and are only weakly correlated with lifestyle. However, within a given lifestyle a few general trends exist, such as the expansion of secreted family M14 metallopeptidases and chitin-binding proteins in plant pathogenic Pezizomycotina.

Conclusions

While the secretomes of fungi with similar lifestyles share certain characteristics, the expansion and contraction of gene families is largely lineage specific, and not shared among all fungi of a given lifestyle.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-722) contains supplementary material, which is available to authorized users.

The haustorial transcriptomes of Uromyces appendiculatus and Phakopsora pachyrhizi and their candidate effector families

Haustoria of biotrophic rust fungi are responsible for the uptake of nutrients from their hosts and for the production of secreted proteins, known as effectors, which modulate the host immune system. The identification of the transcriptome of haustoria and an understanding of the functions of expressed genes therefore hold essential keys for the elucidation of fungus-plant interactions and the development of novel fungal control strategies. Here, we purified haustoria from infected leaves and used 454 sequencing to examine the haustorial transcriptomes of Phakopsora pachyrhizi and Uromyces appendiculatus, the causal agents of soybean rust and common bean rust, respectively. These pathogens cause extensive yield losses in their respective legume crop hosts. A series of analyses were used to annotate expressed sequences, including transposable elements and viruses, to predict secreted proteins from the assembled sequences and to identify families of candidate effectors. This work provides a foundation for the comparative analysis of haustorial gene expression with further insights into physiology and effector evolution.

The genome sequence and effector complement of the flax rust pathogen Melampsora lini

Rust fungi cause serious yield reductions on crops, including wheat, barley, soybean, coffee, and represent real threats to global food security. Of these fungi, the flax rust pathogen Melampsora lini has been developed most extensively over the past 80 years as a model to understand the molecular mechanisms that underpin pathogenesis. During infection, M. lini secretes virulence effectors to promote disease. The number of these effectors, their function and their degree of conservation across rust fungal species is unknown. To assess this, we sequenced and assembled de novo the genome of M. lini isolate CH5 into 21,130 scaffolds spanning 189 Mbp (scaffold N50 of 31 kbp). Global analysis of the DNA sequence revealed that repetitive elements, primarily retrotransposons, make up at least 45% of the genome. Using ab initio predictions, transcriptome data and homology searches, we identified 16,271 putative protein-coding genes. An analysis pipeline was then implemented to predict the effector complement of M. lini and compare it to that of the poplar rust, wheat stem rust and wheat stripe rust pathogens to identify conserved and species-specific effector candidates. Previous knowledge of four cloned M. lini avirulence effector proteins and two basidiomycete effectors was used to optimize parameters of the effector prediction pipeline. Markov clustering based on sequence similarity was performed to group effector candidates from all four rust pathogens. Clusters containing at least one member from M. lini were further analyzed and prioritized based on features including expression in isolated haustoria and infected leaf tissue and conservation across rust species. Herein, we describe 200 of 940 clusters that ranked highest on our priority list, representing 725 flax rust candidate effectors. Our findings on this important model rust species provide insight into how effectors of rust fungi are conserved across species and how they may act to promote infection on their hosts.

Using transcription of six Puccinia triticina races to identify the effective secretome during infection of wheat

Wheat leaf rust, caused by the basidiomycete Puccinia triticina, can cause yield losses of up to 20% in wheat producing regions. During infection, the fungus forms haustoria that secrete proteins into the plant cell and effect changes in plant transcription, metabolism, and defense. It is hypothesized that new races emerge as a result of overcoming plant resistance via changes in the secreted effector proteins. To understand gene expression during infection and find genetic differences associated with races, RNA from wheat leaves infected with six different rust races, at 6 days post inoculation, was sequenced using Illumina. As P. triticina is an obligate biotroph, RNA from both the host and fungi were present and separated by alignment to the P. triticina genome and a wheat EST reference. A total of 222,571 rust contigs were assembled from 165 million reads. An examination of the resulting contigs revealed 532 predicted secreted proteins among the transcripts. Of these, 456 were found in all races. Fifteen genes were found with amino acid changes, corresponding to putative avirulence effectors potentially recognized by 11 different leaf rust resistance (Lr) genes. Twelve of the potential avirulence effectors have no homology to known genes. One gene had significant similarity to cerato-platanin, a known fungal elicitor, and another showed similarity to fungal tyrosinase, an enzyme involved in melanin synthesis. Temporal expression profiles were developed for these genes by qRT-PCR and show that the genes expression patterns were consistent between races from infection initiation to just prior to spore eruption.

Annotation of a hybrid partial genome of the coffee rust (Hemileia vastatrix) contributes to the gene repertoire catalog of the Pucciniales

Coffee leaf rust caused by the fungus Hemileia vastatrix is the most damaging disease to coffee worldwide. The pathogen has recently appeared in multiple outbreaks in coffee producing countries resulting in significant yield losses and increases in costs related to its control. New races/isolates are constantly emerging as evidenced by the presence of the fungus in plants that were previously resistant. Genomic studies are opening new avenues for the study of the evolution of pathogens, the detailed description of plant-pathogen interactions and the development of molecular techniques for the identification of individual isolates. For this purpose we sequenced 8 different H. vastatrix isolates using NGS technologies and gathered partial genome assemblies due to the large repetitive content in the coffee rust hybrid genome; 74.4% of the assembled contigs harbor repetitive sequences. A hybrid assembly of 333 Mb was built based on the 8 isolates; this assembly was used for subsequent analyses. Analysis of the conserved gene space showed that the hybrid H. vastatrix genome, though highly fragmented, had a satisfactory level of completion with 91.94% of core protein-coding orthologous genes present. RNA-Seq from urediniospores was used to guide the de novo annotation of the H. vastatrix gene complement. In total, 14,445 genes organized in 3921 families were uncovered; a considerable proportion of the predicted proteins (73.8%) were homologous to other Pucciniales species genomes. Several gene families related to the fungal lifestyle were identified, particularly 483 predicted secreted proteins that represent candidate effector genes and will provide interesting hints to decipher virulence in the coffee rust fungus. The genome sequence of Hva will serve as a template to understand the molecular mechanisms used by this fungus to attack the coffee plant, to study the diversity of this species and for the development of molecular markers to distinguish races/isolates.

Patterns of genomic variation in the poplar rust fungus Melampsora larici-populina identify pathogenesis-related factors

Melampsora larici-populina is a fungal pathogen responsible for foliar rust disease on poplar trees, which causes damage to forest plantations worldwide, particularly in Northern Europe. The reference genome of the isolate 98AG31 was previously sequenced using a whole genome shotgun strategy, revealing a large genome of 101 megabases containing 16,399 predicted genes, which included secreted protein genes representing poplar rust candidate effectors. In the present study, the genomes of 15 isolates collected over the past 20 years throughout the French territory, representing distinct virulence profiles, were characterized by massively parallel sequencing to assess genetic variation in the poplar rust fungus. Comparison to the reference genome revealed striking structural variations. Analysis of coverage and sequencing depth identified large missing regions between isolates related to the mating type loci. More than 611,824 single-nucleotide polymorphism (SNP) positions were uncovered overall, indicating a remarkable level of polymorphism. Based on the accumulation of non-synonymous substitutions in coding sequences and the relative frequencies of synonymous and non-synonymous polymorphisms (i.e., PN/PS ), we identify candidate genes that may be involved in fungal pathogenesis. Correlation between non-synonymous SNPs in genes encoding secreted proteins (SPs) and pathotypes of the studied isolates revealed candidate genes potentially related to virulences 1, 6, and 8 of the poplar rust fungus.

Duplications and losses in gene families of rust pathogens highlight putative effectors

Rust fungi are a group of fungal pathogens that cause some of the world's most destructive diseases of trees and crops. A shared characteristic among rust fungi is obligate biotrophy, the inability to complete a lifecycle without a host. This dependence on a host species likely affects patterns of gene expansion, contraction, and innovation within rust pathogen genomes. The establishment of disease by biotrophic pathogens is reliant upon effector proteins that are encoded in the fungal genome and secreted from the pathogen into the host's cell apoplast or within the cells. This study uses a comparative genomic approach to elucidate putative effectors and determine their evolutionary histories. We used OrthoMCL to identify nearly 20,000 gene families in proteomes of 16 diverse fungal species, which include 15 basidiomycetes and one ascomycete. We inferred patterns of duplication and loss for each gene family and identified families with distinctive patterns of expansion/contraction associated with the evolution of rust fungal genomes. To recognize potential contributors for the unique features of rust pathogens, we identified families harboring secreted proteins that: (i) arose or expanded in rust pathogens relative to other fungi, or (ii) contracted or were lost in rust fungal genomes. While the origin of rust fungi appears to be associated with considerable gene loss, there are many gene duplications associated with each sampled rust fungal genome. We also highlight two putative effector gene families that have expanded in Cqf that we hypothesize have roles in pathogenicity.

Mosaic genome structure of the barley powdery mildew pathogen and conservation of transcriptional programs in divergent hosts

Barley powdery mildew, Blumeria graminis f. sp. hordei (Bgh), is an obligate biotrophic ascomycete fungal pathogen that can grow and reproduce only on living cells of wild or domesticated barley (Hordeum sp.). Domestication and deployment of resistant barley cultivars by humans selected for amplification of Bgh isolates with different virulence combinations. We sequenced the genomes of two European Bgh isolates, A6 and K1, for comparative analysis with the reference genome of isolate DH14. This revealed a mosaic genome structure consisting of large isolate-specific DNA blocks with either high or low SNP densities. Some of the highly polymorphic blocks likely accumulated SNPs for over 10,000 years, well before the domestication of barley. These isolate-specific blocks of alternating monomorphic and polymorphic regions imply an exceptionally large standing genetic variation in the Bgh population and might be generated and maintained by rare outbreeding and frequent clonal reproduction. RNA-sequencing experiments with isolates A6 and K1 during four early stages of compatible and incompatible interactions on leaves of partially immunocompromised Arabidopsis mutants revealed a conserved Bgh transcriptional program during pathogenesis compared with the natural host barley despite ~200 million years of reproductive isolation of these hosts. Transcripts encoding candidate-secreted effector proteins are massively induced in successive waves. A specific decrease in candidate-secreted effector protein transcript abundance in the incompatible interaction follows extensive transcriptional reprogramming of the host transcriptome and coincides with the onset of localized host cell death, suggesting a host-inducible defense mechanism that targets fungal effector secretion or production.

Unequal recombination and evolution of the mating-type (MAT) loci in the pathogenic fungus Grosmannia clavigera and relatives

Sexual reproduction in fungi is regulated by the mating-type (MAT) locus where recombination is suppressed. We investigated the evolution of MAT loci in eight fungal species belonging to Grosmannia and Ophiostoma (Sordariomycetes, Ascomycota) that include conifer pathogens and beetle symbionts. The MAT1-2 idiomorph/allele was identified from the assembled and annotated Grosmannia clavigera genome, and the MAT locus is flanked by genes coding for cytoskeleton protein (SLA) and DNA lyase. The synteny of these genes is conserved and consistent with other members in Ascomycota. Using sequences from SLA and flanking regions, we characterized the MAT1-1 idiomorph from other isolates of G. clavigera and performed dotplot analysis between the two idiomorphs. Unexpectedly, the MAT1-2 idiomorph contains a truncated MAT1-1-1 gene upstream of the MAT1-2-1 gene that bears the high-mobility-group domain. The nucleotide and amino acid sequence of the truncated MAT1-1-1 gene is similar to its homologous copy in the MAT1-1 idiomorph in the opposite mating-type isolate, except that positive selection is acting on the truncated gene and the alpha(α)-box that encodes the transcription factor has been deleted. The MAT idiomorphs sharing identical gene organization were present in seven additional species in the Ophiostomatales, suggesting that the presence of truncated MAT1-1-1 gene is a general pattern in this order. We propose that an ancient unequal recombination event resulted in the ancestral MAT1-1-1 gene integrated into the MAT1-2 idiomorph and surviving as the truncated MAT1-1-1 genes. The α-box domain of MAT1-1-1 gene, located at the same MAT locus adjacent to the MAT1-2-1 gene, could have been removed by deletion after recombination due to mating signal interference. Our data confirmed a 1:1 MAT/sex ratio in two pathogen populations, and showed that all members of the Ophiostomatales studied here including those that were previously deemed asexual have the potential to reproduce sexually. This ability can potentially increase genetic variability and can enhance fitness in new, ecological niches.

Methodologies to decipher the cell secretome

The cell secretome is a collection of proteins consisting of transmembrane proteins (TM) and proteins secreted by cells into the extracellular space. A significant portion (~13-20%) of the human proteome consists of secretory proteins. The secretory proteins play important roles in cell migration, cell signaling and communication. There is a plethora of methodologies available like Serial Analysis of Gene Expression (SAGE), DNA microarrays, antibody arrays and bead-based arrays, mass spectrometry, RNA sequencing and yeast, bacterial and mammalian secretion traps to identify the cell secretomes. There are many advantages and disadvantages in using any of the above methods. This review aims to discuss the methodologies available along with their potential advantages and disadvantages to identify secretory proteins. This review is a part of a Special issue on The Secretome. This article is part of a Special Issue entitled: An Updated Secretome.

The rust transferred proteins-a new family of effector proteins exhibiting protease inhibitor function

Only few fungal effectors have been described to be delivered into the host cell during obligate biotrophic interactions. RTP1p, from the rust fungi Uromyces fabae and U. striatus, was the first fungal protein for which localization within the host cytoplasm could be demonstrated directly. We investigated the occurrence of RTP1 homologues in rust fungi and examined the structural and biochemical characteristics of the corresponding gene products. The analysis of 28 homologues showed that members of the RTP family are most likely to occur ubiquitously in rust fungi and to be specific to the order Pucciniales. Sequence analyses indicated that the structure of the RTPp effectors is bipartite, consisting of a variable N-terminus and a conserved and structured C-terminus. The characterization of Uf-RTP1p mutants showed that four conserved cysteine residues sustain structural stability. Furthermore, the C-terminal domain exhibits similarities to that of cysteine protease inhibitors, and it was shown that Uf-RTP1p and Us-RTP1p are able to inhibit proteolytic activity in Pichia pastoris culture supernatants. We conclude that the RTP1p homologues constitute a rust fungi-specific family of modular effector proteins comprising an unstructured N-terminal domain and a structured C-terminal domain, which exhibit protease inhibitory activity possibly associated with effector function during biotrophic interactions.

Toward understanding the role of aphid effectors in plant infestation

In recent years, immense progress has been made toward understanding the functions of effectors from a range of plant pathogens, such as oomycetes, fungi, bacteria, and nematodes. Like plant pathogens, aphids form close associations with host plants, featuring signal exchange between the two organisms. While feeding and probing, aphids deliver effector proteins mixed with saliva directly into the host-stylet interface. With the increasing availability of aphid genome and transcriptome sequence data, aphid effector biology is emerging as a new and exciting area of research. In this review, we provide an overview of recent advances in the aphid effector biology field and highlight some of the current questions.

Defining the predicted protein secretome of the fungal wheat leaf pathogen Mycosphaerella graminicola

The Dothideomycete fungus Mycosphaerella graminicola is the causal agent of Septoria tritici blotch, a devastating disease of wheat leaves that causes dramatic decreases in yield. Infection involves an initial extended period of symptomless intercellular colonisation prior to the development of visible necrotic disease lesions. Previous functional genomics and gene expression profiling studies have implicated the production of secreted virulence effector proteins as key facilitators of the initial symptomless growth phase. In order to identify additional candidate virulence effectors, we re-analysed and catalogued the predicted protein secretome of M. graminicola isolate IPO323, which is currently regarded as the reference strain for this species. We combined several bioinformatic approaches in order to increase the probability of identifying truly secreted proteins with either a predicted enzymatic function or an as yet unknown function. An initial secretome of 970 proteins was predicted, whilst further stringent selection criteria predicted 492 proteins. Of these, 321 possess some functional annotation, the composition of which may reflect the strictly intercellular growth habit of this pathogen, leaving 171 with no functional annotation. This analysis identified a protein family encoding secreted peroxidases/chloroperoxidases (PF01328) which is expanded within all members of the family Mycosphaerellaceae. Further analyses were done on the non-annotated proteins for size and cysteine content (effector protein hallmarks), and then by studying the distribution of homologues in 17 other sequenced Dothideomycete fungi within an overall total of 91 predicted proteomes from fungal, oomycete and nematode species. This detailed M. graminicola secretome analysis provides the basis for further functional and comparative genomics studies.

mRNA-Seq analysis of the Pseudoperonospora cubensis transcriptome during cucumber (Cucumis sativus L.) infection

Pseudoperonospora cubensis, an oomycete, is the causal agent of cucurbit downy mildew, and is responsible for significant losses on cucurbit crops worldwide. While other oomycete plant pathogens have been extensively studied at the molecular level, Ps. cubensis and the molecular basis of its interaction with cucurbit hosts has not been well examined. Here, we present the first large-scale global gene expression analysis of Ps. cubensis infection of a susceptible Cucumis sativus cultivar, ‘Vlaspik’, and identification of genes with putative roles in infection, growth, and pathogenicity. Using high throughput whole transcriptome sequencing, we captured differential expression of 2383 Ps. cubensis genes in sporangia and at 1, 2, 3, 4, 6, and 8 days post-inoculation (dpi). Additionally, comparison of Ps. cubensis expression profiles with expression profiles from an infection time course of the oomycete pathogen Phytophthora infestans on Solanum tuberosum revealed similarities in expression patterns of 1,576–6,806 orthologous genes suggesting a substantial degree of overlap in molecular events in virulence between the biotrophic Ps. cubensis and the hemi-biotrophic P. infestans. Co-expression analyses identified distinct modules of Ps. cubensis genes that were representative of early, intermediate, and late infection stages. Collectively, these expression data have advanced our understanding of key molecular and genetic events in the virulence of Ps. cubensis and thus, provides a foundation for identifying mechanism(s) by which to engineer or effect resistance in the host.

A comprehensive analysis of genes encoding small secreted proteins identifies candidate effectors in Melampsora larici-populina (poplar leaf rust)

The obligate biotrophic rust fungus Melampsora larici-populina is the most devastating and widespread pathogen of poplars. Studies over recent years have identified various small secreted proteins (SSP) from plant biotrophic filamentous pathogens and have highlighted their role as effectors in host-pathogen interactions. The recent analysis of the M. larici-populina genome sequence has revealed the presence of 1,184 SSP-encoding genes in this rust fungus. In the present study, the expression and evolutionary dynamics of these SSP were investigated to pinpoint the arsenal of putative effectors that could be involved in the interaction between the rust fungus and poplar. Similarity with effectors previously described in Melampsora spp., richness in cysteines, and organization in large families were extensively detailed and discussed. Positive selection analyses conducted over clusters of paralogous genes revealed fast-evolving candidate effectors. Transcript profiling of selected M. laricipopulina SSP showed a timely coordinated expression during leaf infection, and the accumulation of four candidate effectors in distinct rust infection structures was demonstrated by immunolocalization. This integrated and multifaceted approach helps to prioritize candidate effector genes for functional studies.

Using hierarchical clustering of secreted protein families to classify and rank candidate effectors of rust fungi

Rust fungi are obligate biotrophic pathogens that cause considerable damage on crop plants. Puccinia graminis f. sp. tritici, the causal agent of wheat stem rust, and Melampsora larici-populina, the poplar leaf rust pathogen, have strong deleterious impacts on wheat and poplar wood production, respectively. Filamentous pathogens such as rust fungi secrete molecules called disease effectors that act as modulators of host cell physiology and can suppress or trigger host immunity. Current knowledge on effectors from other filamentous plant pathogens can be exploited for the characterisation of effectors in the genome of recently sequenced rust fungi. We designed a comprehensive in silico analysis pipeline to identify the putative effector repertoire from the genome of two plant pathogenic rust fungi. The pipeline is based on the observation that known effector proteins from filamentous pathogens have at least one of the following properties: (i) contain a secretion signal, (ii) are encoded by in planta induced genes, (iii) have similarity to haustorial proteins, (iv) are small and cysteine rich, (v) contain a known effector motif or a nuclear localization signal, (vi) are encoded by genes with long intergenic regions, (vii) contain internal repeats, and (viii) do not contain PFAM domains, except those associated with pathogenicity. We used Markov clustering and hierarchical clustering to classify protein families of rust pathogens and rank them according to their likelihood of being effectors. Using this approach, we identified eight families of candidate effectors that we consider of high value for functional characterization. This study revealed a diverse set of candidate effectors, including families of haustorial expressed secreted proteins and small cysteine-rich proteins. This comprehensive classification of candidate effectors from these devastating rust pathogens is an initial step towards probing plant germplasm for novel resistance components.

454-pyrosequencing of Coffea arabica leaves infected by the rust fungus Hemileia vastatrix reveals in planta-expressed pathogen-secreted proteins and plant functions in a late compatible plant-rust interaction

Coffee (Coffea arabica L.), one of the key export and cash crops in tropical and subtropical countries, suffers severe losses from the rust fungus Hemileia vastatrix. The transcriptome of H. vastatrix was analysed during a compatible interaction with coffee to obtain an exhaustive repertoire of the genes expressed during infection and to identify potential effector genes. Large-scale sequencing (454-GS-FLEX Titanium) of mixed coffee and rust cDNAs obtained from 21-day rust-infected leaves generated 352 146 sequences which assembled into 22 774 contigs. In the absence of any reference genomic sequences for Coffea or Hemileia, specific trinucleotide frequencies within expressed sequence tags (ESTs) and blast homology against a set of dicots and basidiomycete genomes were used to distinguish pathogen from plant sequences. About 30% (6763) of the contigs were assigned to H. vastatrix and 61% (13 951) to C. arabica. The majority (60%) of the rust sequences did not show homology to any genomic database, indicating that they were potential novel fungal genes. In silico analyses of the 6763 H. vastatrix contigs predicted 382 secreted proteins and identified homologues of the flax rust haustorially expressed secreted proteins (HESPs) and bean rust transferred protein 1 (RTP1). These rust candidate effectors showed conserved amino-acid domains and conserved patterns of cysteine positions suggestive of conserved functions during infection of host plants. Quantitative reverse transcription-polymerase chain reaction profiling of selected rust genes revealed dynamic expression patterns during the time course of infection of coffee leaves. This study provides the first valuable genomic resource for the agriculturally important plant pathogen H. vastatrix and the first comprehensive C. arabica EST dataset.

Meta-analytic approach to the accurate prediction of secreted virulence effectors in gram-negative bacteria

Background

Many pathogens use a type III secretion system to translocate virulence proteins (called effectors) in order to adapt to the host environment. To date, many prediction tools for effector identification have been developed. However, these tools are insufficiently accurate for producing a list of putative effectors that can be applied directly for labor-intensive experimental verification. This also suggests that important features of effectors have yet to be fully characterized.

Results

In this study, we have constructed an accurate approach to predicting secreted virulence effectors from Gram-negative bacteria. This consists of a support vector machine-based discriminant analysis followed by a simple criteria-based filtering. The accuracy was assessed by estimating the average number of true positives in the top-20 ranking in the genome-wide screening. In the validation, 10 sets of 20 training and 20 testing examples were randomly selected from 40 known effectors of Salmonella enterica serovar Typhimurium LT2. On average, the SVM portion of our system predicted 9.7 true positives from 20 testing examples in the top-20 of the prediction. Removal of the N-terminal instability, codon adaptation index and ProtParam indices decreased the score to 7.6, 8.9 and 7.9, respectively. These discrimination features suggested that the following characteristics of effectors had been uncovered: unstable N-terminus, non-optimal codon usage, hydrophilic, and less aliphathic. The secondary filtering process represented by coexpression analysis and domain distribution analysis further refined the average true positive counts to 12.3. We further confirmed that our system can correctly predict known effectors of P. syringae DC3000, strongly indicating its feasibility.

Conclusions

We have successfully developed an accurate prediction system for screening effectors on a genome-wide scale. We confirmed the accuracy of our system by external validation using known effectors of Salmonella and obtained the accurate list of putative effectors of the organism. The level of accuracy was sufficient to yield candidates for gene-directed experimental verification. Furthermore, new features of effectors were revealed: non-optimal codon usage and instability of the N-terminal region. From these findings, a new working hypothesis is proposed regarding mechanisms controlling the translocation of virulence effectors and determining the substrate specificity encoded in the secretion system.

The poplar-poplar rust interaction: insights from genomics and transcriptomics

Poplars are extensively cultivated worldwide, and their susceptibility to the leaf rust fungus Melampsora larici-populina leads to considerable damages in plantations. Despite a good knowledge of the poplar rust life cycle, and particularly the epidemics on poplar, the perennial status of the plant host and the obligate biotrophic lifestyle of the rust fungus are bottlenecks for molecular investigations. Following the completion of both M. larici-populina and Populus trichocarpa genome sequences, gene families involved in poplar resistance or in rust fungus virulence were investigated, allowing the identification of key genetic determinants likely controlling the outcome of the interaction. Specific expansions of resistance and defense-related genes in poplar indicate probable innovations in perennial species in relation with host-pathogen interactions. The genome of M. Larici-populina contains a strikingly high number of genes encoding small secreted proteins (SSPs) representing hundreds of candidate effectors. Transcriptome analyses of interacting partners in compatible and incompatible interactions revealed conserved set of genes involved in poplar defense reactions as well as timely regulated expression of SSP transcripts during host tissues colonisation. Ongoing functional studies of selected candidate effectors will be achieved mainly on the basis of recombinant protein purification and subsequent characterisation.

Melampsora larici-populina transcript profiling during germination and timecourse infection of poplar leaves reveals dynamic expression patterns associated with virulence and biotrophy

Melampsora larici-populina is responsible for poplar leaf rust disease and causes severe epidemics in poplar plantations in Europe. The poplar rust genome has been recently sequenced and, in order to find the genetic determinants associated with its biotrophic lifestyle, we generated a whole-genome custom oligoarray and analyzed transcript profiles of M. larici-populina during the infection timecourse in poplar leaves. Different stages were investigated during the asexual development of the rust fungus, including resting and germinating urediniospores and seven in planta stages in the telial host. In total, 76% of the transcripts were detected during leaf infection as well as in urediniospores, whereas 20% were only detected in planta, including several transporters and many small secreted proteins (SSP). We focused our analysis on gene categories known to be related to plant colonization and biotrophic growth in rust pathogens, such as SSP, carbohydrate active enzymes (CAZymes), transporters, lipases, and proteases. Distinct sets of SSP transcripts were expressed all along the infection process, suggesting highly dynamic expression of candidate rust effectors. In contrast, transcripts encoding transporters and proteases were mostly expressed after 48 h postinoculation, when numerous haustoria are already formed in the leaf mesophyll until uredinia formation, supporting their role in nutrient acquisition during biotrophic growth. Finally, CAZymes and lipase transcripts were predominantly expressed at late stages of infection, highlighting their importance during sporulation.

Gene discovery in EST sequences from the wheat leaf rust fungus Puccinia triticina sexual spores, asexual spores and haustoria, compared to other rust and corn smut fungi

Background

Rust fungi are biotrophic basidiomycete plant pathogens that cause major diseases on plants and trees world-wide, affecting agriculture and forestry. Their biotrophic nature precludes many established molecular genetic manipulations and lines of research. The generation of genomic resources for these microbes is leading to novel insights into biology such as interactions with the hosts and guiding directions for breakthrough research in plant pathology.

Results

To support gene discovery and gene model verification in the genome of the wheat leaf rust fungus, Puccinia triticina (Pt), we have generated Expressed Sequence Tags (ESTs) by sampling several life cycle stages. We focused on several spore stages and isolated haustorial structures from infected wheat, generating 17,684 ESTs. We produced sequences from both the sexual (pycniospores, aeciospores and teliospores) and asexual (germinated urediniospores) stages of the life cycle. From pycniospores and aeciospores, produced by infecting the alternate host, meadow rue (Thalictrum speciosissimum), 4,869 and 1,292 reads were generated, respectively. We generated 3,703 ESTs from teliospores produced on the senescent primary wheat host. Finally, we generated 6,817 reads from haustoria isolated from infected wheat as well as 1,003 sequences from germinated urediniospores. Along with 25,558 previously generated ESTs, we compiled a database of 13,328 non-redundant sequences (4,506 singlets and 8,822 contigs). Fungal genes were predicted using the EST version of the self-training GeneMarkS algorithm. To refine the EST database, we compared EST sequences by BLASTN to a set of 454 pyrosequencing-generated contigs and Sanger BAC-end sequences derived both from the Pt genome, and to ESTs and genome reads from wheat. A collection of 6,308 fungal genes was identified and compared to sequences of the cereal rusts, Puccinia graminis f. sp. tritici (Pgt) and stripe rust, P. striiformis f. sp. tritici (Pst), and poplar leaf rust Melampsora species, and the corn smut fungus, Ustilago maydis (Um). While extensive homologies were found, many genes appeared novel and species-specific; over 40% of genes did not match any known sequence in existing databases. Focusing on spore stages, direct comparison to Um identified potential functional homologs, possibly allowing heterologous functional analysis in that model fungus. Many potentially secreted protein genes were identified by similarity searches against genes and proteins of Pgt and Melampsora spp., revealing apparent orthologs.

Conclusions

The current set of Pt unigenes contributes to gene discovery in this major cereal pathogen and will be invaluable for gene model verification in the genome sequence.

Proteome analysis of wheat leaf rust fungus, Puccinia triticina, infection structures enriched for haustoria

Puccinia triticina (Pt) is a representative of several cereal-infecting rust fungal pathogens of major economic importance world wide. Upon entry through leaf stomata, these fungi establish intracellular haustoria, crucial feeding structures. We report the first proteome of infection structures from parasitized wheat leaves, enriched for haustoria through filtration and sucrose density centrifugation. 2-D PAGE MS/MS and gel-based LC-MS (GeLC-MS) were used to separate proteins. Generated spectra were compared with a partial proteome predicted from a preliminary Pt genome and generated ESTs, to a comprehensive genome-predicted protein complement from the related wheat stem rust fungus, Puccinia graminis f. sp. tritici (Pgt) and to various plant resources. We identified over 260 fungal proteins, 16 of which matched peptides from Pgt. Based on bioinformatic analyses and/or the presence of a signal peptide, at least 50 proteins were predicted to be secreted. Among those, six have effector protein signatures, some are related and the respective genes of several seem to belong to clusters. Many ribosomal structural proteins, proteins involved in energy, general metabolism and transport were detected. Measuring gene expression over several life cycle stages of ten representative candidates using quantitative RT-PCR, all were shown to be strongly upregulated and four expressed solely upon infection.