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

High-throughput RNA sequencing reveals differences between the transcriptomes of the five spore forms of Puccinia striiformis f. sp. tritici, the wheat stripe rust pathogen

The devastating wheat stripe (yellow) rust pathogen, Puccinia striiformis f. sp. tritici (Pst), is a macrocyclic and heteroecious fungus. Pst produces urediniospores and teliospores on its primary host, wheat, and pycniospores and aeciospores are produced on its alternate hosts, barberry (Berberis spp.) or mahonia (Mahonia spp.). Basidiospores are developed from teliospores and infect alternate hosts. These five spore forms play distinct roles in Pst infection, disease development, and fungal survival, etc. However, the specific genes and mechanisms underlying these functional differences are largely unknown. In this study, we performed, for the first time in rust fungi, the deep RNA sequencing to examine the transcriptomic shift among all five Pst spore forms. Among a total of 29,591 identified transcripts, 951 were specifically expressed in basidiospores, whereas 920, 761, 266, and 110 were specific for teliospores, pycniospores, aeciospores, and urediniospores, respectively. Additionally, transcriptomes of sexual spores, namely pycniospores and basidiospores, showed significant differences from those of asexual spores (urediniospores, teliospores, and aeciospores), and transcriptomes of urediniospores and aeciospores were more similar to each other than to the three other spore forms. Especially, the basidiospores and pycniospores which infected the berberis shows wide differences in the cell wall degrading-enzymes and mating and pheromone response genes. Besides, we also found that there are 6234 differential expressed genes between the urediniospores and pycniospores, while only have 3 genes have alternative splicing enents, suggesting that differential genes expression may make more contribution than AS. This comprehensive transcriptome profiling can substantially improve our understanding of the developmental biology of the wheat stripe rust fungus.

Molecular Mechanisms of the Co-Evolution of Wheat and Rust Pathogens

Wheat (Triticum spp.) is a cereal crop domesticated >8000 years ago and the second-most-consumed food crop nowadays. Ever since mankind has written records, cereal rust diseases have been a painful awareness in antiquity documented in the Old Testament (about 750 B.C.). The pathogen causing the wheat stem rust disease is among the first identified plant pathogens in the 1700s, suggesting that wheat and rust pathogens have co-existed for thousands of years. With advanced molecular technologies, wheat and rust genomes have been sequenced, and interactions between the host and the rust pathogens have been extensively studied at molecular levels. In this review, we summarized the research at the molecular level and organized the findings based on the pathogenesis steps of germination, penetration, haustorial formation, and colonization of the rusts to present the molecular mechanisms of the co-evolution of wheat and rust pathogens.

Comparative Genome Analyses of Plant Rust Pathogen Genomes Reveal a Confluence of Pathogenicity Factors to Quell Host Plant Defense Responses

Switchgrass rust caused by Puccinia novopanici (P. novopanici) has the ability to significantly affect the biomass yield of switchgrass, an important biofuel crop in the United States. A comparative genome analysis of P. novopanici with rust pathogen genomes infecting monocot cereal crops wheat, barley, oats, maize and sorghum revealed the presence of larger structural variations contributing to their genome sizes. A comparative alignment of the rust pathogen genomes resulted in the identification of collinear and syntenic relationships between P. novopanici and P. sorghi; P. graminis tritici 21–0 (Pgt 21) and P. graminis tritici Ug99 (Pgt Ug99) and between Pgt 21 and P. triticina (Pt). Repeat element analysis indicated a strong presence of retro elements among different Puccinia genomes, contributing to the genome size variation between ~1 and 3%. A comparative look at the enriched protein families of Puccinia spp. revealed a predominant role of restriction of telomere capping proteins (RTC), disulfide isomerases, polysaccharide deacetylases, glycoside hydrolases, superoxide dismutases and multi-copper oxidases (MCOs). All the proteomes of Puccinia spp. share in common a repertoire of 75 secretory and 24 effector proteins, including glycoside hydrolases cellobiohydrolases, peptidyl-propyl isomerases, polysaccharide deacetylases and protein disulfide-isomerases, that remain central to their pathogenicity. Comparison of the predicted effector proteins from Puccinia spp. genomes to the validated proteins from the Pathogen–Host Interactions database (PHI-base) resulted in the identification of validated effector proteins PgtSR1 (PGTG_09586) from P. graminis and Mlp124478 from Melampsora laricis across all the rust pathogen genomes.

Physical separation of haplotypes in dikaryons allows benchmarking of phasing accuracy in Nanopore and HiFi assemblies with Hi-C data

Background

Most animals and plants have more than one set of chromosomes and package these haplotypes into a single nucleus within each cell. In contrast, many fungal species carry multiple haploid nuclei per cell. Rust fungi are such species with two nuclei (karyons) that contain a full set of haploid chromosomes each. The physical separation of haplotypes in dikaryons means that, unlike in diploids, Hi-C chromatin contacts between haplotypes are false-positive signals.

Results

We generate the first chromosome-scale, fully-phased assembly for the dikaryotic leaf rust fungus Puccinia triticina and compare Nanopore MinION and PacBio HiFi sequence-based assemblies. We show that false-positive Hi-C contacts between haplotypes are predominantly caused by phase switches rather than by collapsed regions or Hi-C read mis-mappings. We introduce a method for phasing of dikaryotic genomes into the two haplotypes using Hi-C contact graphs, including a phase switch correction step. In the HiFi assembly, relatively few phase switches occur, and these are predominantly located at haplotig boundaries and can be readily corrected. In contrast, phase switches are widespread throughout the Nanopore assembly. We show that haploid genome read coverage of 30–40 times using HiFi sequencing is required for phasing of the leaf rust genome, with 0.7% heterozygosity, and that HiFi sequencing resolves genomic regions with low heterozygosity that are otherwise collapsed in the Nanopore assembly.

Conclusions

This first Hi-C based phasing pipeline for dikaryons and comparison of long-read sequencing technologies will inform future genome assembly and haplotype phasing projects in other non-haploid organisms.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02658-2.

Host Adaptation and Virulence in Heteroecious Rust Fungi

Rust fungi (Pucciniales, Basidiomycota) are obligate biotrophic pathogens that cause rust diseases in plants, inflicting severe damage to agricultural crops. Pucciniales possess the most complex life cycles known in fungi. These include an alternation of generations, the development of up to five different sporulating stages, and, for many species, the requirement of infecting two unrelated host plants during different parts of their life cycle, termed heteroecism. These fungi have been extensively studied in the past century through microscopy and inoculation studies, providing precise descriptions of their infection processes, although the molecular mechanisms underlying their unique biology are poorly understood. In this review, we cover recent genomic and life cycle transcriptomic studies in several heteroecious rust species, which provide insights into the genetic tool kits associated with host adaptation and virulence, opening new avenues for unraveling their unique evolution.

A Chromosome-Scale Assembly of the Wheat Leaf Rust Pathogen Puccinia triticina Provides Insights Into Structural Variations and Genetic Relationships With Haplotype Resolution

Despite the global economic importance of the wheat leaf rust pathogen Puccinia triticina (Pt), genomic resources for Pt are limited and chromosome-level assemblies of Pt are lacking. Here, we present a complete haplotype-resolved genome assembly at a chromosome-scale for Pt using the Australian pathotype 64-(6),(7),(10),11 (Pt64; North American race LBBQB) built upon the newly developed technologies of PacBio and Hi-C sequencing. PacBio reads with ∼200-fold coverage (29.8 Gb data) were assembled by Falcon and Falcon-unzip and subsequently scaffolded with Hi-C data using Falcon-phase and Proximo. This approach allowed us to construct 18 chromosome pseudomolecules ranging from 3.5 to 12.3 Mb in size for each haplotype of the dikaryotic genome of Pt64. Each haplotype had a total length of ∼147 Mb, scaffold N 50 of ∼9.4 Mb, and was ∼93% complete for BUSCOs. Each haplotype had ∼29,800 predicted genes, of which ∼2,000 were predicted as secreted proteins (SPs). The investigation of structural variants (SVs) between haplotypes A and B revealed that 10% of the total genome was spanned by SVs, highlighting variations previously undetected by short-read based assemblies. For the first time, the mating type (MAT) genes on each haplotype of Pt64 were identified, which showed that MAT loci a and b are located on two chromosomes (chromosomes 7 and 14), representing a tetrapolar type. Furthermore, the Pt64 assembly enabled haplotype-based evolutionary analyses for 21 Australian Pt isolates, which highlighted the importance of a haplotype resolved reference when inferring genetic relationships using whole genome SNPs. This Pt64 assembly at chromosome-scale with full phase information provides an invaluable resource for genomic and evolutionary research, which will accelerate the understanding of molecular mechanisms underlying Pt-wheat interactions and facilitate the development of durable resistance to leaf rust in wheat and sustainable control of rust disease.

Effector Biology of Biotrophic Plant Fungal Pathogens: Current Advances and Future Prospects

The interaction of fungal pathogens with their host requires a novel invading mechanism and the presence of various virulence-associated components responsible for promoting the infection. The small secretory proteins, explicitly known as effector proteins, are one of the prime mechanisms of host manipulation utilized by the pathogen to disarm the host. Several effector proteins are known to translocate from fungus to the plant cell for host manipulation. Many fungal effectors have been identified using genomic, transcriptomic, and bioinformatics approaches. Most of the effector proteins are devoid of any conserved signatures, and their prediction based on sequence homology is very challenging, therefore by combining the sequence consensus based upon machine learning features, multiple tools have also been developed for predicting apoplastic and cytoplasmic effectors. Various post-genomics approaches like transcriptomics of virulent isolates have also been utilized for identifying active consortia of effectors. Significant progress has been made in understanding biotrophic effectors; however, most of it is underway due to their complex interaction with host and complicated recognition and signaling networks. This review discusses advances, and challenges in effector identification and highlighted various features of the potential effector proteins and approaches for understanding their genetics and strategies for regulation.

Long-Read-Based de novo Genome Assembly and Comparative Genomics of the Wheat Leaf Rust Pathogen Puccinia triticina Identifies Candidates for Three Avirulence Genes

Leaf rust, caused by Puccinia triticina (Pt), is one of the most devastating diseases of wheat, affecting production in nearly all wheat-growing regions worldwide. Despite its economic importance, genomic resources for Pt are very limited. In the present study, we have used long-read sequencing (LRS) and the pipeline of FALCON and FALCON-Unzip (v4.1.0) to carry out the first LRS-based de novo genome assembly for Pt. Using 22.4-Gb data with an average read length of 11.6 kb and average coverage of 150-fold, we generated a genome assembly for Pt104 [strain 104-2,3,(6),(7),11; isolate S423], considered to be the founding isolate of a clonal lineage of Pt in Australia. The Pt104 genome contains 162 contigs with a total length of 140.5 Mb and N50 of 2 Mb, with the associated haplotigs providing haplotype information for 91% of the genome. This represents the best quality of Pt genome assembly to date, which reduces the contig number by 91-fold and improves the N50 by 4-fold as compared to the previous Pt race1 assembly. An annotation pipeline that combined multiple lines of evidence including the transcriptome assemblies derived from RNA-Seq, previously identified expressed sequence tags and Pt race 1 protein sequences predicted 29,043 genes for Pt104 genome. Based on the presence of a signal peptide, no transmembrane segment, and no target location to mitochondria, 2,178 genes were identified as secreted proteins (SPs). Whole-genome sequencing (Illumina paired-end) was performed for Pt104 and six additional strains with differential virulence profile on the wheat leaf rust resistance genes Lr26, Lr2a, and Lr3ka. To identify candidates for the corresponding avirulence genes AvrLr26, AvrLr2a, and AvrLr3ka, genetic variation within each strain was first identified by mapping to the Pt104 genome. Variants within predicted SP genes between the strains were then correlated to the virulence profiles, identifying 38, 31, and 37 candidates for AvrLr26, AvrLr2a, and AvrLr3ka, respectively. The identification of these candidate genes lays a good foundation for future studies on isolating these avirulence genes, investigating the molecular mechanisms underlying host-pathogen interactions, and the development of new diagnostic tools for pathogen monitoring.

Breeding Wheat for Durable Leaf Rust Resistance in Southern Africa: Variability, Distribution, Current Control Strategies, Challenges and Future Prospects

Leaf or brown rust of wheat caused by Puccinia triticina (Pt) is one of the most damaging diseases globally. Considerable progress has been made to control leaf rust through crop protection chemicals and host plant resistance breeding in southern Africa. However, frequent changes in the pathogen population still present a major challenge to achieve durable resistance. Disease surveillance and monitoring of the pathogen have revealed the occurrence of similar races across the region, justifying the need for concerted efforts by countries in southern Africa to develop and deploy more efficient and sustainable strategies to manage the disease. Understanding the genetic variability and composition of Pt is a pre-requisite for cultivar release with appropriate resistance gene combinations for sustainable disease management. This review highlights the variability and distribution of the Pt population, and the current control strategies, challenges and future prospects of breeding wheat varieties with durable leaf rust resistance in southern Africa. The importance of regular, collaborative and efficient surveillance of the pathogen and germplasm development across southern Africa is discussed, coupled with the potential of using modern breeding technologies to produce wheat cultivars with durable resistance.

Progress on Molecular Genetics and Manipulation of Rust Fungi

Among the thousands of rust species described, many are known for their devastating effects on their hosts, which include major agriculture crops and trees. Hence, for over a century, these basidiomycete pathogenic fungi have been researched and experimented with. However, due to their biotrophic nature, they are challenging organisms to work with and, needing their hosts for propagation, represent pathosystems that are not easily experimentally accessible. Indeed, efforts to perform genetics have been few and far apart for the rust fungi, though one study performed in the 1940s was famously instrumental in formulating the gene-for-gene hypothesis describing pathogen-host interactions. By taking full advantage of the molecular genetic tools developed in the 1980s, research on many plant pathogenic microbes thrived, yet similar work on the rusts remained very challenging though not without some successes. However, the genomics era brought real breakthrough research for the biotrophic fungi and with innovative experimentation and the use of heterologous systems, molecular genetic analyses over the last 2 decades have significantly advanced our insight into the function of many rust fungus genes and their role in the interaction with their hosts. This has allowed optimizing efforts for resistance breeding and the design and testing of various novel strategies to reduce the devastating diseases they cause.

The progress of leaf rust research in wheat

Leaf rust (also called brown rust) in wheat, caused by fungal pathogen Puccinia triticina Erikss. (Pt) is one of the major constraints in wheat production worldwide. Pt is widespread with diverse population structure and undergoes rapid evolution to produce new virulent races against resistant cultivars that are regularly developed to provide resistance against the prevailing races of the pathogen. Occasionally, the disease may also take the shape of an epidemic in some wheat-growing areas causing major economic losses. In the recent past, substantial progress has been made in characterizing the sources of leaf rust resistance including non-host resistance (NHR). Progress has also been made in elucidating the population biology of Pt and the mechanisms of wheat-Pt interaction. So far, ∼80 leaf rust resistance genes (Lr genes) have been identified and characterized; some of them have also been used for the development of resistant wheat cultivars. It has also been shown that a gene-for-gene relationship exists between individual wheat Lr genes and the corresponding Pt Avr genes so that no Lr gene can provide resistance unless the prevailing race of the pathogen carries the corresponding Avr gene. Several Lr genes have also been cloned and their products characterized, although no Avr gene corresponding a specific Lr gene has so far been identified. However, several candidate effectors for Pt have been identified and functionally characterized using genome-wide analyses, transcriptomics, RNA sequencing, bimolecular fluorescence complementation (BiFC), virus-induced gene silencing (VIGS), transient expression and other approaches. This review summarizes available information on different aspects of the pathogen Pt, genetics/genomics of leaf rust resistance in wheat including cloning and characterization of Lr genes and epigenetic regulation of disease resistance.

Comparative transcriptomics of Gymnosporangium spp. teliospores reveals a conserved genetic program at this specific stage of the rust fungal life cycle

BACKGROUND

Gymnosporangium spp. are fungal plant pathogens causing rust disease and most of them are known to infect two different host plants (heteroecious) with four spore stages (demicyclic). In the present study, we sequenced the transcriptome of G. japonicum teliospores on its host plant Juniperus chinensis and we performed comparison to the transcriptomes of G. yamadae and G. asiaticum at the same life stage, that happens in the same host but on different organs.

RESULTS

Functional annotation for the three Gymnosporangium species showed the expression of a conserved genetic program with the top abundant cellular categories corresponding to energy, translation and signal transduction processes, indicating that this life stage is particularly active. Moreover, the survey of predicted secretomes in the three Gymnosporangium transcriptomes revealed shared and specific genes encoding carbohydrate active enzymes and secreted proteins of unknown function that could represent candidate pathogenesis effectors. A transcript encoding a hemicellulase of the glycoside hydrolase 26 family, previously identified in other rust fungi, was particularly highly expressed suggesting a general role in rust fungi. The comparison between the transcriptomes of the three Gymnosporangium spp. and selected Pucciniales species in different taxonomical families allowed to identify lineage-specific protein families that may relate to the biology of teliospores in rust fungi. Among clustered gene families, 205, 200 and 152 proteins were specifically identified in G. japonicum, G. yamadae and G. asiaticum, respectively, including candidate effectors expressed in teliospores.

CONCLUSIONS

This comprehensive comparative transcriptomics study of three Gymnosporangium spp. identified gene functions and metabolic pathways particularly expressed in teliospores, a stage of the life cycle that is mostly overlooked in rust fungi. Secreted protein encoding transcripts expressed in teliospores may reveal new candidate effectors related to pathogenesis. Although this spore stage is not involved in host plant infection but in the production of basidiospores infecting plants in the Amygdaloideae, we speculate that candidate effectors may be expressed as early as the teliospore stage for preparing further infection by basidiospores.

Advances in understanding obligate biotrophy in rust fungi

Contents Summary 1190 I. Introduction 1190 II. Rust fungi: a diverse and serious threat to agriculture 1191 III. The different facets of rust life cycles and unresolved questions about their evolution 1191 IV. The biology of rust infection 1192 V. Rusts in the genomics era: the ever-expanding list of candidate effector genes 1195 VI. Functional characterization of rust effectors 1197 VII. Putting rusts to sleep: Pucciniales research outlooks 1201 Acknowledgements 1202 References 1202 SUMMARY: Rust fungi (Pucciniales) are the largest group of plant pathogens and represent one of the most devastating threats to agricultural crops worldwide. Despite the economic importance of these highly specialized pathogens, many aspects of their biology remain obscure, largely because rust fungi are obligate biotrophs. The rise of genomics and advances in high-throughput sequencing technology have presented new options for identifying candidate effector genes involved in pathogenicity mechanisms of rust fungi. Transcriptome analysis and integrated bioinformatics tools have led to the identification of key genetic determinants of host susceptibility to infection by rusts. Thousands of genes encoding secreted proteins highly expressed during host infection have been reported for different rust species, which represents significant potential towards understanding rust effector function. Recent high-throughput in planta expression screen approaches (effectoromics) have pushed the field ahead even further towards predicting high-priority effectors and identifying avirulence genes. These new insights into rust effector biology promise to inform future research and spur the development of effective and sustainable strategies for managing rust diseases.

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.

Host-induced gene silencing of an important pathogenicity factor PsCPK1 in Puccinia striiformis f. sp. tritici enhances resistance of wheat to stripe rust

Rust fungi are devastating plant pathogens and cause a large economic impact on wheat production worldwide. To overcome this rapid loss of resistance in varieties, we generated stable transgenic wheat plants expressing short interfering RNAs (siRNAs) targeting potentially vital genes of Puccinia striiformis f. sp. tritici (Pst). Protein kinase A (PKA) has been proved to play important roles in regulating the virulence of phytopathogenic fungi. PsCPK1, a PKA catalytic subunit gene from Pst, is highly induced at the early infection stage of Pst. The instantaneous silencing of PsCPK1 by barley stripe mosaic virus (BSMV)-mediated host-induced gene silencing (HIGS) results in a significant reduction in the length of infection hyphae and disease phenotype. These results indicate that PsCPK1 is an important pathogenicity factor by regulating Pst growth and development. Two transgenic lines expressing the RNA interference (RNAi) construct in a normally susceptible wheat cultivar displayed high levels of stable and consistent resistance to Pst throughout the T3 to T4 generations. The presence of the interfering RNAs in transgenic wheat plants was confirmed by northern blotting, and these RNAs were found to efficiently down-regulate PsCPK1 expression in wheat. This study addresses important aspects for the development of fungal-derived resistance through the expression of silencing constructs in host plants as a powerful strategy to control cereal rust diseases.

Comparative transcriptome analysis and identification of candidate effectors in two related rust species (Gymnosporangium yamadae and Gymnosporangium asiaticum)

BACKGROUND

Rust fungi constitute the largest group of plant fungal pathogens. However, a paucity of data, including genomic sequences, transcriptome sequences, and associated molecular markers, hinders the development of inhibitory compounds and prevents their analysis from an evolutionary perspective. Gymnosporangium yamadae and G. asiaticum are two closely related rust fungal species, which are ecologically and economically important pathogens that cause apple rust and pear rust, respectively, proved to be devastating to orchards. In this study, we investigated the transcriptomes of these two Gymnosporangium species during the telial stage of their lifecycles. The aim of this study was to understand the evolutionary patterns of these two related fungi and to identify genes that developed by selection.

RESULTS

The transcriptomes of G. yamadae and G. asiaticum were generated from a mixture of RNA from three biological replicates of each species. We obtained 49,318 and 54,742 transcripts, with N50 values of 1957 and 1664, for G. yamadae and G. asiaticum, respectively. We also identified a repertoire of candidate effectors and other gene families associated with pathogenicity. A total of 4947 pairs of putative orthologues between the two species were identified. Estimation of the non-synonymous/synonymous substitution rate ratios for these orthologues identified 116 pairs with Ka/Ks values greater than1 that are under positive selection and 170 pairs with Ka/Ks values of 1 that are under neutral selection, whereas the remaining 4661 genes are subjected to purifying selection. We estimate that the divergence time between the two species is approximately 5.2 Mya.

CONCLUSION

This study constitutes a de novo assembly and comparative analysis between the transcriptomes of the two rust species G. yamadae and G. asiaticum. The results identified several orthologous genes, and many expressed genes were identified by annotation. Our analysis of Ka/Ks ratios identified orthologous genes subjected to positive or purifying selection. An evolutionary analysis of these two species provided a relatively precise divergence time. Overall, the information obtained in this study increases the genetic resources available for research on the genetic diversity of the Gymnosporangium genus.

Divergent and convergent modes of interaction between wheat and Puccinia graminis f. sp. tritici isolates revealed by the comparative gene co-expression network and genome analyses

BACKGROUND

Two opposing evolutionary constraints exert pressure on plant pathogens: one to diversify virulence factors in order to evade plant defenses, and the other to retain virulence factors critical for maintaining a compatible interaction with the plant host. To better understand how the diversified arsenals of fungal genes promote interaction with the same compatible wheat line, we performed a comparative genomic analysis of two North American isolates of Puccinia graminis f. sp. tritici (Pgt).

RESULTS

The patterns of inter-isolate divergence in the secreted candidate effector genes were compared with the levels of conservation and divergence of plant-pathogen gene co-expression networks (GCN) developed for each isolate. Comprative genomic analyses revealed substantial level of interisolate divergence in effector gene complement and sequence divergence. Gene Ontology (GO) analyses of the conserved and unique parts of the isolate-specific GCNs identified a number of conserved host pathways targeted by both isolates. Interestingly, the degree of inter-isolate sub-network conservation varied widely for the different host pathways and was positively associated with the proportion of conserved effector candidates associated with each sub-network. While different Pgt isolates tended to exploit similar wheat pathways for infection, the mode of plant-pathogen interaction varied for different pathways with some pathways being associated with the conserved set of effectors and others being linked with the diverged or isolate-specific effectors.

CONCLUSIONS

Our data suggest that at the intra-species level pathogen populations likely maintain divergent sets of effectors capable of targeting the same plant host pathways. This functional redundancy may play an important role in the dynamic of the "arms-race" between host and pathogen serving as the basis for diverse virulence strategies and creating conditions where mutations in certain effector groups will not have a major effect on the pathogen's ability to infect the host.

Prediction of the in planta Phakopsora pachyrhizi secretome and potential effector families

Asian soybean rust (ASR), caused by the obligate biotrophic fungus Phakopsora pachyrhizi, can cause losses greater than 80%. Despite its economic importance, there is no soybean cultivar with durable ASR resistance. In addition, the P. pachyrhizi genome is not yet available. However, the availability of other rust genomes, as well as the development of sample enrichment strategies and bioinformatics tools, has improved our knowledge of the ASR secretome and its potential effectors. In this context, we used a combination of laser capture microdissection (LCM), RNAseq and a bioinformatics pipeline to identify a total of 36 350 P. pachyrhizi contigs expressed in planta and a predicted secretome of 851 proteins. Some of the predicted secreted proteins had characteristics of candidate effectors: small size, cysteine rich, do not contain PFAM domains (except those associated with pathogenicity) and strongly expressed in planta. A comparative analysis of the predicted secreted proteins present in Pucciniales species identified new members of soybean rust and new Pucciniales- or P. pachyrhizi-specific families (tribes). Members of some families were strongly up-regulated during early infection, starting with initial infection through haustorium formation. Effector candidates selected from two of these families were able to suppress immunity in transient assays, and were localized in the plant cytoplasm and nuclei. These experiments support our bioinformatics predictions and show that these families contain members that have functions consistent with P. pachyrhizi effectors.

Comparative Analysis Highlights Variable Genome Content of Wheat Rusts and Divergence of the Mating Loci

Three members of the Puccinia genus, Pucciniatriticina (Pt), Pstriiformis f.sp. tritici (Pst), and Pgraminis f.sp. tritici (Pgt), cause the most common and often most significant foliar diseases of wheat. While similar in biology and life cycle, each species is uniquely adapted and specialized. The genomes of Pt and Pst were sequenced and compared to that of Pgt to identify common and distinguishing gene content, to determine gene variation among wheat rust pathogens, other rust fungi, and basidiomycetes, and to identify genes of significance for infection. Pt had the largest genome of the three, estimated at 135 Mb with expansion due to mobile elements and repeats encompassing 50.9% of contig bases; in comparison, repeats occupy 31.5% for Pst and 36.5% for Pgt We find all three genomes are highly heterozygous, with Pst [5.97 single nucleotide polymorphisms (SNPs)/kb] nearly twice the level detected in Pt (2.57 SNPs/kb) and that previously reported for Pgt Of 1358 predicted effectors in Pt, 784 were found expressed across diverse life cycle stages including the sexual stage. Comparison to related fungi highlighted the expansion of gene families involved in transcriptional regulation and nucleotide binding, protein modification, and carbohydrate degradation enzymes. Two allelic homeodomain pairs, HD1 and HD2, were identified in each dikaryotic Puccinia species along with three pheromone receptor (STE3) mating-type genes, two of which are likely representing allelic specificities. The HD proteins were active in a heterologous Ustilago maydis mating assay and host-induced gene silencing (HIGS) of the HD and STE3 alleles reduced wheat host infection.

Investigating Gene Function in Cereal Rust Fungi by Plant-Mediated Virus-Induced Gene Silencing

Cereal rust fungi are destructive pathogens, threatening grain production worldwide. Targeted breeding for resistance utilizing host resistance genes has been effective. However, breakdown of resistance occurs frequently and continued efforts are needed to understand how these fungi overcome resistance and to expand the range of available resistance genes. Whole genome sequencing, transcriptomic and proteomic studies followed by genome-wide computational and comparative analyses have identified large repertoire of genes in rust fungi among which are candidates predicted to code for pathogenicity and virulence factors. Some of these genes represent defence triggering avirulence effectors. However, functions of most genes still needs to be assessed to understand the biology of these obligate biotrophic pathogens. Since genetic manipulations such as gene deletion and genetic transformation are not yet feasible in rust fungi, performing functional gene studies is challenging. Recently, Host-induced gene silencing (HIGS) has emerged as a useful tool to characterize gene function in rust fungi while infecting and growing in host plants. We utilized Barley stripe mosaic virus-mediated virus induced gene silencing (BSMV-VIGS) to induce HIGS of candidate rust fungal genes in the wheat host to determine their role in plant-fungal interactions. Here, we describe the methods for using BSMV-VIGS in wheat for functional genomics study in cereal rust fungi.

Effectors from Wheat Rust Fungi Suppress Multiple Plant Defense Responses

Fungi that cause cereal rust diseases (genus Puccinia) are important pathogens of wheat globally. Upon infection, the fungus secretes a number of effector proteins. Although a large repository of putative effectors has been predicted using bioinformatic pipelines, the lack of available high-throughput effector screening systems has limited functional studies on these proteins. In this study, we mined the available transcriptomes of Puccinia graminis and P. striiformis to look for potential effectors that suppress host hypersensitive response (HR). Twenty small (<300 amino acids), secreted proteins, with no predicted functions were selected for the HR suppression assay using Nicotiana benthamiana, in which each of the proteins were transiently expressed and evaluated for their ability to suppress HR caused by four cytotoxic effector-R gene combinations (Cp/Rx, ATR13/RPP13, Rpt2/RPS-2, and GPA/RBP-1) and one mutated R gene-Pto(Y207D). Nine out of twenty proteins, designated Shr1 to Shr9 (suppressors of hypersensitive response), were found to suppress HR in N. benthamiana. These effectors varied in the effector-R gene defenses they suppressed, indicating these pathogens can interfere with a variety of host defense pathways. In addition to HR suppression, effector Shr7 also suppressed PAMP-triggered immune response triggered by flg22. Finally, delivery of Shr7 through Pseudomonas fluorescens EtHAn suppressed nonspecific HR induced by Pseudomonas syringae DC3000 in wheat, confirming its activity in a homologous system. Overall, this study provides the first evidence for the presence of effectors in Puccinia species suppressing multiple plant defense responses.

Draft Genome of the Wheat Rust Pathogen (Puccinia triticina) Unravels Genome-Wide Structural Variations during Evolution

Leaf rust is one of the most important diseases of wheat and is caused by Puccinia triticina, a highly variable rust pathogen prevalent worldwide. Decoding the genome of this pathogen will help in unraveling the molecular basis of its evolution and in the identification of genes responsible for its various biological functions. We generated high quality draft genome sequences (approximately 100- 106 Mb) of two races of P. triticina; the variable and virulent Race77 and the old, avirulent Race106. The genomes of races 77 and 106 had 33X and 27X coverage, respectively. We predicted 27678 and 26384 genes, with average lengths of 1,129 and 1,086 bases in races 77 and 106, respectively and found that the genomes consisted of 37.49% and 39.99% repetitive sequences. Genome wide comparative analysis revealed that Race77 differs substantially from Race106 with regard to segmental duplication (SD), repeat element, and SNP/InDel characteristics. Comparative analyses showed that Race 77 is a recent, highly variable and adapted Race compared with Race106. Further sequence analyses of 13 additional pathotypes of Race77 clearly differentiated the recent, active and virulent, from the older pathotypes. Average densities of 2.4 SNPs and 0.32 InDels per kb were obtained for all P. triticina pathotypes. Secretome analysis demonstrated that Race77 has more virulence factors than Race 106, which may be responsible for the greater degree of adaptation of this pathogen. We also found that genes under greater selection pressure were conserved in the genomes of both races, and may affect functions crucial for the higher levels of virulence factors in Race77. This study provides insights into the genome structure, genome organization, molecular basis of variation, and pathogenicity of P. triticina The genome sequence data generated in this study have been submitted to public domain databases and will be an important resource for comparative genomics studies of the more than 4000 existing Puccinia species.

Characterization of protein kinase PsSRPKL, a novel pathogenicity factor in the wheat stripe rust fungus

As in other eukaryotes, protein kinases (PKs) are generally evolutionarily conserved and play major regulatory roles in plant pathogenic fungi. Many PKs have been proven to be important for pathogenesis in model fungal plant pathogens, but little is currently known about their roles in the pathogenesis of cereal rust fungi, devastating pathogens in agriculture worldwide. Here, we report on an in planta highly induced PK gene PsSRPKL from the wheat stripe rust fungus Puccinia striiformis f. sp. tritici (Pst), one of the most important cereal rust fungi. PsSRPKL belongs to a group of PKs that are evolutionarily specific to cereal rust fungi. It shows a high level of intraspecies polymorphism in the kinase domains and directed green fluorescent protein chimers to plant nuclei. Overexpression of PsSRPKL in fission yeast induces aberrant cell morphology and a decreased resistance to environmental stresses. Most importantly, PsSRPKL is proven to be an important pathogenicity factor responsible for fungal growth and responses to environmental stresses, therefore contributing significantly to Pst virulence in wheat. We hypothesize that cereal rust fungi have developed specific PKs as pathogenicity factors for adaptation to their host species during evolution. Thus, our findings provide significant insights into pathogenicity and virulence evolution in cereal rust fungi.

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.

Effector proteins of rust fungi

Rust fungi include many species that are devastating crop pathogens. To develop resistant plants, a better understanding of rust virulence factors, or effector proteins, is needed. Thus far, only six rust effector proteins have been described: AvrP123, AvrP4, AvrL567, AvrM, RTP1, and PGTAUSPE-10-1. Although some are well established model proteins used to investigate mechanisms of immune receptor activation (avirulence activities) or entry into plant cells, how they work inside host tissues to promote fungal growth remains unknown. The genome sequences of four rust fungi (two Melampsoraceae and two Pucciniaceae) have been analyzed so far. Genome-wide analyses of these species, as well as transcriptomics performed on a broader range of rust fungi, revealed hundreds of small secreted proteins considered as rust candidate secreted effector proteins (CSEPs). The rust community now needs high-throughput approaches (effectoromics) to accelerate effector discovery/characterization and to better understand how they function in planta. However, this task is challenging due to the non-amenability of rust pathosystems (obligate biotrophs infecting crop plants) to traditional molecular genetic approaches mainly due to difficulties in culturing these species in vitro. The use of heterologous approaches should be promoted in the future.

Transcriptome analysis of poplar rust telia reveals overwintering adaptation and tightly coordinated karyogamy and meiosis processes

Most rust fungi have a complex life cycle involving up to five different spore-producing stages. The telial stage that produces melanized overwintering teliospores is one of these and plays a fundamental role for generating genetic diversity as karyogamy and meiosis occur at that stage. Despite the importance of telia for the rust life cycle, almost nothing is known about the fungal genetic programs that are activated in this overwintering structure. In the present study, the transcriptome of telia produced by the poplar rust fungus Melampsora larici-populina has been investigated using whole genome exon oligoarrays and RT-qPCR. Comparative expression profiling at the telial and uredinial stages identifies genes specifically expressed or up-regulated in telia including osmotins/thaumatin-like proteins (TLPs) and aquaporins that may reflect specific adaptation to overwintering as well numerous lytic enzymes acting on plant cell wall, reflecting extensive cell wall remodeling at that stage. The temporal dynamics of karyogamy was followed using combined RT-qPCR and DAPI-staining approaches. This reveals that fusion of nuclei and induction of karyogamy-related genes occur simultaneously between the 25 and 39 days post inoculation time frame. Transcript profiling of conserved meiosis genes indicates a preferential induction right after karyogamy and corroborates that meiosis begins prior to overwintering and is interrupted in Meiosis I (prophase I, diplonema stage) until teliospore germination in early spring.

Analysis of Phakopsora pachyrhizi transcript abundance in critical pathways at four time-points during infection of a susceptible soybean cultivar using deep sequencing

BACKGROUND

Phakopsora pachyrhizi, the causal agent responsible for soybean rust, is among the top hundred most virulent plant pathogens and can cause soybean yield losses of up to 80% when appropriate conditions are met. We used mRNA-Seq by Illumina to analyze pathogen transcript abundance at 15 seconds (s), 7 hours (h), 48 h, and 10 days (d) after inoculation (ai) of susceptible soybean leaves with P. pachyrhizi to gain new insights into transcript abundance in soybean and the pathogen at specific time-points during the infection including the uredinial stage.

RESULTS

Over three million five hundred thousand sequences were obtained for each time-point. Energy, nucleotide metabolism, and protein synthesis are major priorities for the fungus during infection and development as indicated by our transcript abundance studies. At all time-points, energy production is a necessity for P. pachyrhizi, as indicated by expression of many transcripts encoding enzymes involved in oxidative phosphorylation and carbohydrate metabolism (glycolysis, glyoxylate and dicarboxylate, pentose phosphate, pyruvate). However, at 15 sai, transcripts encoding enzymes involved in ATP production were highly abundant in order to provide enough energy for the spore to germinate, as observed by the expression of many transcripts encoding proteins involved in electron transport. At this early time-point, transcripts encoding proteins involved in RNA synthesis were also highly abundant, more so than transcripts encoding genes involved in DNA and protein synthesis. At 7 hai, shortly after germination during tube elongation and penetration, transcripts encoding enzymes involved in deoxyribonucleotide and DNA synthesis were highly abundant. At 48 hai, transcripts encoding enzymes involved in amino acid metabolism were highly abundant to provide for increased protein synthesis during haustoria maturation. During sporulation at 10 dai, the fungus still required carbohydrate metabolism, but there also was increased expression of transcripts encoding enzymes involved in fatty acid metabolism.

CONCLUSION

This information provides insight into molecular events and their timing throughout the life cycle of the P. pachyrhizi, and it may be useful in the development of new methods of broadening resistance of soybean to soybean rust.

Host-induced gene silencing of wheat leaf rust fungus Puccinia triticina pathogenicity genes mediated by the Barley stripe mosaic virus

Rust fungi are devastating plant pathogens and several Puccinia species have a large economic impact on wheat production worldwide. Disease protection, mostly offered by introgressed host-resistance genes, is often race-specific and rapidly overcome by newly-emerging virulent strains. Extensive new genomic resources have identified vital pathogenicity genes but their study is hampered because of the biotrophic life styles of rust fungi. In cereals, Barley stripe mosaic virus (BSMV)-induced RNAi has emerged as a useful tool to study loss-of-function phenotypes of candidate genes. Expression of pathogen-derived gene fragments in this system can be used to obtain in planta-generated silencing of corresponding genes inside biotrophic pathogens, a technique termed host-induced gene silencing (HIGS). Here we test the effectiveness of BSMV-mediated HIGS in the wheat leaf rust fungus Puccinia triticina (Pt) by targeting three predicted pathogenicity genes, a MAPK, a cyclophilin, and a calcineurin regulatory subunit. Inoculation of BSMV RNAi constructs generated fungal gene-specific siRNA molecules in systemic leaves of wheat plant. Subsequent Pt inoculation resulted in a suppressed disease phenotype and a reduction in endogenous transcript levels of the targeted fungal genes indicating translocation of siRNA molecules from host to fungal cells. Efficiency of this host-generated trans-specific RNAi was enhanced by using BSMV silencing vectors defective in coat protein coupled with introducing fungal gene sequences simultaneously in sense and antisense orientation. The disease suppression indicated the likely involvement of these fungal genes in pathogenicity. This study demonstrates that BSMV-mediated in planta-generated RNAi is an effective strategy for functional genomics in rust fungi.

Endogenous silencing of Puccinia triticina pathogenicity genes through in planta-expressed sequences leads to the suppression of rust diseases on wheat

Rust fungi are destructive plant pathogens. The draft genomes of several wheat-infecting species have been released and potential pathogenicity genes identified through comparative analyses to fungal pathogens that are amenable to genetic manipulation. Functional gene analysis tools are needed to understand the infection process of these obligate parasites and to confirm whether predicted pathogenicity genes could become targets for disease control. We have modified an Agrobacterium tumefaciens-mediated in planta-induced transient gene silencing (PITGS) assay for use in Triticum spp. (wheat), and used this assay to target predicted wheat leaf rust fungus, Puccinia triticina (Pt) pathogenicity genes, a MAP kinase (PtMAPK1), a cyclophilin (PtCYC1) and calcineurin B (PtCNB), to analyze their roles in disease. Agroinfiltration effectively delivered hairpin silencing constructs in wheat, leading to the generation of fungal gene-specific siRNA molecules in infiltrated leaves, and resulting in up to 70% reduction in transcription of the endogenous target genes in superinfected Pt. In vivo silencing caused severe disease suppression, compromising fungal growth and sporulation, as viewed by confocal microscopy and measured by reductions in fungal biomass and emergence of uredinia. Interestingly, using the same gene constructs, suppression of infection by Puccinia graminis and Puccinia striiformis was also achieved. Our results show that A. tumefaciens-mediated PITGS can be used as a reverse-genetics tool to discover gene function in rust fungi. This proof-of-concept study indicates that the targeted fungal transcripts might be important in pathogenesis, and could potentially be used as promising targets for developing RNA interference-based resistance against rust fungi.

Conserved loci of leaf and stem rust fungi of wheat share synteny interrupted by lineage-specific influx of repeat elements

Background

Wheat leaf rust (Puccinia triticina Eriks; Pt) and stem rust fungi (P. graminis f.sp. tritici; Pgt) are significant economic pathogens having similar host ranges and life cycles, but different alternate hosts. The Pt genome, currently estimated at 135 Mb, is significantly larger than Pgt, at 88 Mb, but the reason for the expansion is unknown. Three genomic loci of Pt conserved proteins were characterized to gain insight into gene content, genome complexity and expansion.

Results

A bacterial artificial chromosome (BAC) library was made from P. triticina race 1, BBBD and probed with Pt homologs of genes encoding two predicted Pgt secreted effectors and a DNA marker mapping to a region of avirulence. Three BACs, 103 Kb, 112 Kb, and 166 Kb, were sequenced, assembled, and open reading frames were identified. Orthologous genes were identified in Pgt and local conservation of gene order (microsynteny) was observed. Pairwise protein identities ranged from 26 to 99%. One Pt BAC, containing a RAD18 ortholog, shares syntenic regions with two Pgt scaffolds, which could represent both haplotypes of Pgt. Gene sequence is diverged between the species as well as within the two haplotypes. In all three BAC clones, gene order is locally conserved, however, gene shuffling has occurred relative to Pgt. These regions are further diverged by differing insertion loci of LTR-retrotransposon, Gypsy, Copia, Mutator, and Harbinger mobile elements. Uncharacterized Pt open reading frames were also found; these proteins are high in lysine and similar to multiple proteins in Pgt.

Conclusions

The three Pt loci are conserved in gene order, with a range of gene sequence divergence. Conservation of predicted haustoria expressed secreted protein genes between Pt and Pgt is extended to the more distant poplar rust, Melampsora larici-populina. The loci also reveal that genome expansion in Pt is in part due to higher occurrence of repeat-elements in this species.

Transcriptome analysis of enriched Golovinomyces orontii haustoria by deep 454 pyrosequencing

Powdery mildews are phytopathogenic ascomycetes that have an obligate biotrophic lifestyle and establish intimate relationships with their plant hosts. A crucial aspect of this plant-fungus interaction is the formation of specialized fungal infection structures termed haustoria. Although located within the cell boundaries of plant epidermal cells, haustoria remain separated from the plant cytoplasm by a host plasma membrane derivative, the extrahaustorial membrane. Haustoria are thought to represent pivotal sites of nutrient uptake and effector protein delivery. We enriched haustorial complexes from Arabidopsis thaliana plants infected with the powdery mildew fungus Golovinomyces orontii and performed in-depth transcriptome analysis by 454-based pyrosequencing of haustorial cDNAs. We assembled 7077 expressed sequence tag (EST) contigs with greater than 5-fold average coverage and analyzed these with regard to the respective predicted protein functions. We found that transcripts coding for gene products with roles in protein turnover, detoxification of reactive oxygen species and fungal pathogenesis are abundant in the haustorial EST contigs, while surprisingly transcripts encoding presumptive nutrient transporters were not highly represented in the haustorial cDNA library. A substantial proportion (∼38%) of transcripts coding for predicted secreted proteins comprises effector candidates. Our data provide valuable insights into the transcriptome of the key infection structure of a model obligate biotrophic phytopathogen.

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.

Obligate biotrophy features unraveled by the genomic analysis of rust fungi

Rust fungi are some of the most devastating pathogens of crop plants. They are obligate biotrophs, which extract nutrients only from living plant tissues and cannot grow apart from their hosts. Their lifestyle has slowed the dissection of molecular mechanisms underlying host invasion and avoidance or suppression of plant innate immunity. We sequenced the 101-Mb genome of Melampsora larici-populina, the causal agent of poplar leaf rust, and the 89-Mb genome of Puccinia graminis f. sp. tritici, the causal agent of wheat and barley stem rust. We then compared the 16,399 predicted proteins of M. larici-populina with the 17,773 predicted proteins of P. graminis f. sp tritici. Genomic features related to their obligate biotrophic lifestyle include expanded lineage-specific gene families, a large repertoire of effector-like small secreted proteins, impaired nitrogen and sulfur assimilation pathways, and expanded families of amino acid and oligopeptide membrane transporters. The dramatic up-regulation of transcripts coding for small secreted proteins, secreted hydrolytic enzymes, and transporters in planta suggests that they play a role in host infection and nutrient acquisition. Some of these genomic hallmarks are mirrored in the genomes of other microbial eukaryotes that have independently evolved to infect plants, indicating convergent adaptation to a biotrophic existence inside plant cells.