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

A perspective on varied fungal virulence factors causing infection in host plants

Pathogenic fungi and their spores are ubiquitously present and invade the tissues of higher living plants causing pathogenesis and inevitably death or retarded growth. A group of fungi kills its hosts and consume the dead tissues (necrotrophs), while others feed on living tissue (biotrophs) or combination of two (hemibiotrophs). A number of virulent factors is used by fungal pathogens to inhabit new hosts and cause illness. Fungal pathogens develop specialized structures for complete invasion into plant organs to regulate pathogenic growth. Virulence factors like effectors, mycotoxins, cell wall degrading enzymes and organic acids have varied roles depending on the infection strategy and assist the pathogens to possess control on living tissues of the plants. Infection strategies employed by fungi generally masks the plant defense mechanism, however necrotrophs are best known to harm plant tissues with their poisonous secretion. Interestingly, the effector chemicals released by Biotrophs reduce plant cell growth and regulate plant metabolism in their advantage causing no direct death. All these virulence tools cause huge loss to the agricultural product of pre- harvest crops and post-harvest yields causing low output leading to huge economic losses. This review focusses on comprehensive study of range of virulence factors of the pathogenic fungi responsible for their invasion inside the healthy tissues of plants. The compiled information would influence researchers to design antidote against all virulence factors of fungi relevant to their area of research which could pave way for protection against plant pathogenesis.

Genome biology and evolution of mating-type loci in four cereal rust fungi

Permanent heterozygous loci, such as sex- or mating-compatibility regions, often display suppression of recombination and signals of genomic degeneration. In Basidiomycota, two distinct loci confer mating compatibility. These loci encode homeodomain (HD) transcription factors and pheromone receptor (Pra)-ligand allele pairs. To date, an analysis of genome level mating-type (MAT) loci is lacking for obligate biotrophic basidiomycetes in the Pucciniales, an order containing serious agricultural plant pathogens. Here, we focus on four species of Puccinia that infect oat and wheat, including P. coronata f. sp. avenae, P. graminis f. sp. tritici, P. triticina and P. striiformis f. sp. tritici. MAT loci are located on two separate chromosomes supporting previous hypotheses of a tetrapolar mating compatibility system in the Pucciniales. The HD genes are multiallelic in all four species while the PR locus appears biallelic, except for P. graminis f. sp. tritici, which potentially has multiple alleles. HD loci are largely conserved in their macrosynteny, both within and between species, without strong signals of recombination suppression. Regions proximal to the PR locus, however, displayed signs of recombination suppression and genomic degeneration in the three species with a biallelic PR locus. Our observations support a link between recombination suppression, genomic degeneration, and allele diversity of MAT loci that is consistent with recent mathematical modelling and simulations. Finally, we confirm that MAT genes are expressed during the asexual infection cycle, and we propose that this may support regulating nuclear maintenance and pairing during infection and spore formation. Our study provides insights into the evolution of MAT loci of key pathogenic Puccinia species. Understanding mating compatibility can help predict possible combinations of nuclear pairs, generated by sexual reproduction or somatic recombination, and the potential evolution of new virulent isolates of these important plant pathogens.

Overview and Management of the Most Common Eukaryotic Diseases of Flax (Linum usitatissimum)

Flax is an important crop cultivated for its seeds and fibers. It is widely grown in temperate regions, with an increase in cultivation areas for seed production (linseed) in the past 50 years and for fiber production (fiber flax) in the last decade. Among fiber-producing crops, fiber flax is the most valuable species. Linseed is the highest omega-3 oleaginous crop, and its consumption provides several benefits for animal and human health. However, flax production is impacted by various abiotic and biotic factors that affect yield and quality. Among biotic factors, eukaryotic diseases pose a significant threat to both seed production and fiber quality, which highlights the economic importance of controlling these diseases. This review focuses on the major eukaryotic diseases that affect flax in the field, describing the pathogens, their transmission modes and the associated plant symptoms. Moreover, this article aims to identify the challenges in disease management and provide future perspectives to overcome these biotic stresses in flax cultivation. By emphasizing the key diseases and their management, this review can aid in promoting sustainable and profitable flax production.

Molecular Advances to Combat Different Biotic and Abiotic Stresses in Linseed (Linum usitatissimum L.): A Comprehensive Review

Flax, or linseed, is considered a "superfood", which means that it is a food with diverse health benefits and potentially useful bioactive ingredients. It is a multi-purpose crop that is prized for its seed oil, fibre, nutraceutical, and probiotic qualities. It is suited to various habitats and agro-ecological conditions. Numerous abiotic and biotic stressors that can either have a direct or indirect impact on plant health are experienced by flax plants as a result of changing environmental circumstances. Research on the impact of various stresses and their possible ameliorators is prompted by such expectations. By inducing the loss of specific alleles and using a limited number of selected varieties, modern breeding techniques have decreased the overall genetic variability required for climate-smart agriculture. However, gene banks have well-managed collectionns of landraces, wild linseed accessions, and auxiliary Linum species that serve as an important source of novel alleles. In the past, flax-breeding techniques were prioritised, preserving high yield with other essential traits. Applications of molecular markers in modern breeding have made it easy to identify quantitative trait loci (QTLs) for various agronomic characteristics. The genetic diversity of linseed species and the evaluation of their tolerance to abiotic stresses, including drought, salinity, heavy metal tolerance, and temperature, as well as resistance to biotic stress factors, viz., rust, wilt, powdery mildew, and alternaria blight, despite addressing various morphotypes and the value of linseed as a supplement, are the primary topics of this review.

Prediction of effector protein structures from fungal phytopathogens enables evolutionary analyses

Elucidating the similarity and diversity of pathogen effectors is critical to understand their evolution across fungal phytopathogens. However, rapid divergence that diminishes sequence similarities between putatively homologous effectors has largely concealed the roots of effector evolution. Here we modelled the structures of 26,653 secreted proteins from 14 agriculturally important fungal phytopathogens, six non-pathogenic fungi and one oomycete with AlphaFold 2. With 18,000 successfully predicted folds, we performed structure-guided comparative analyses on two aspects of effector evolution: uniquely expanded sequence-unrelated structurally similar (SUSS) effector families and common folds present across the fungal species. Extreme expansion of lineage-specific SUSS effector families was found only in several obligate biotrophs, Blumeria graminis and Puccinia graminis. The highly expanded effector families were the source of conserved sequence motifs, such as the Y/F/WxC motif. We identified new classes of SUSS effector families that include known virulence factors, such as AvrSr35, AvrSr50 and Tin2. Structural comparisons revealed that the expanded structural folds further diversify through domain duplications and fusion with disordered stretches. Putatively sub- and neo-functionalized SUSS effectors could reconverge on regulation, expanding the functional pools of effectors in the pathogen infection cycle. We also found evidence that many effector families could have originated from ancestral folds conserved across fungi. Collectively, our study highlights diverse effector evolution mechanisms and supports divergent evolution as a major force in driving SUSS effector evolution from ancestral proteins.

Biotrophic Fungal Pathogens: a Critical Overview

Biotrophic fungi are one group of heterogeneous organisms and these fungi differ in their traits like mode of nutrition, types of reproduction, and dispersal systems. Generally, based on the nutritional mode, fungi are classified into three broad categories, viz. biotrophs, necrotrophs, and hemi-biotrophs. Biotrophs derive their nutrients and energy from living plant cells and survive within the interstitial space of the cells. Biotrophic fungi cause serious crop diseases but are highly challenging to investigate and develop a treatment strategy. Blumeria (Erysiphe) graminis, Uromyces fabae, Ustilago maydis, Cladosporium fulvum, Puccinia graminis, and Phytophthora infestans are some of the significant biotrophic fungi that affect mainly plants. One among the biotrophic fungus, Pneumocystis jirovecii (Taphrinomycotina subphylum of the Ascomycota) exclusively a human pathogen, can cause lung diseases such as "pneumocystis." Biotrophic fungus widely parasitizing Solanaceae family crops (Tomato and potato) has done massive damage to the crops and has led to economic impact worldwide. During infection and for nutrient absorption, biotrophs develops external appendages such as appressoria or haustoria. The hyphae or appressorium adheres to the plant cell wall and collapses the layers for their nutrient absorption. The pathogen also secretes effector molecules to escape from the plant defense mechanism. Later, plants activate their primary and secondary defense mechanisms; however, the pathogen induces virulence genes to escape the host immune responses. Obligate biotrophic fungi pathogenicity has not been fully understood at the molecular level because of the complex interaction, recognition, and signaling with the host. This review summarizes the mechanism of infection in the host, and immune response to emphasize the understanding of the biotrophic fungal biology and pathogenesis in crops. Thus, the detailed review will pave the way to design methods to overcome the resistance of biotrophic fungi and develop disease-free crops.

A decade after the first Pucciniales genomes: A bibliometric snapshot of (post) genomics studies in three model rust fungi

Pucciniales (rust fungi) are one of the largest fungal order of plant pathogens. They collectively infect key crops such as wheat and soybean, and threaten global food security. In the early 2010s, the genome sequences of three rust fungi were released: Melampsora larici-populina (the poplar leaf rust fungus), Puccinia graminis f. sp. tritici (the wheat stem rust fungus), and Puccinia striiformis f. sp. triciti (the wheat stripe rust or wheat yellow rust fungus). The availability of those genomes has forwarded rust biology into the post-genomic era, sparking a series of genomics, transcriptomics, in silico, and functional studies. Here, we snapshot the last 10 years of post-genomics studies addressing M. larici-populina, P. graminis f. sp. tritici, and/or P. striiformis f. sp. tritici. This mini-review notably reveals the model species-centered structure of the research community, and highlights the drastic increase of the number of functional studies focused on effectors since 2014, which notably revealed chloroplasts as a central host compartment targeted by rust fungi. This mini-review also discusses genomics-facilitated studies in other rust species, and emerging post-genomic research trends related to fully-phased rust genomes.

Current Status and Future Perspectives of Genomics Research in the Rust Fungi

Rust fungi in Pucciniales have caused destructive plant epidemics, have become more aggressive with new virulence, rapidly adapt to new environments, and continually threaten global agriculture. With the rapid advancement of genome sequencing technologies and data analysis tools, genomics research on many of the devastating rust fungi has generated unprecedented insights into various aspects of rust biology. In this review, we first present a summary of the main findings in the genomics of rust fungi related to variations in genome size and gene composition between and within species. Then we show how the genomics of rust fungi has promoted our understanding of the pathogen virulence and population dynamics. Even with great progress, many questions still need to be answered. Therefore, we introduce important perspectives with emphasis on the genome evolution and host adaptation of rust fungi. We believe that the comparative genomics and population genomics of rust fungi will provide a further understanding of the rapid evolution of virulence and will contribute to monitoring the population dynamics for disease management.

Integrated omics approaches for flax improvement under abiotic and biotic stress: Current status and future prospects

Flax (Linum usitatissimum L.) or linseed is one of the important industrial crops grown all over the world for seed oil and fiber. Besides oil and fiber, flax offers a wide range of nutritional and therapeutic applications as a feed and food source owing to high amount of α-linolenic acid (omega-3 fatty acid), lignans, protein, minerals, and vitamins. Periodic losses caused by unpredictable environmental stresses such as drought, heat, salinity-alkalinity, and diseases pose a threat to meet the rising market demand. Furthermore, these abiotic and biotic stressors have a negative impact on biological diversity and quality of oil/fiber. Therefore, understanding the interaction of genetic and environmental factors in stress tolerance mechanism and identification of underlying genes for economically important traits is critical for flax improvement and sustainability. In recent technological era, numerous omics techniques such as genomics, transcriptomics, metabolomics, proteomics, phenomics, and ionomics have evolved. The advancements in sequencing technologies accelerated development of genomic resources which facilitated finer genetic mapping, quantitative trait loci (QTL) mapping, genome-wide association studies (GWAS), and genomic selection in major cereal and oilseed crops including flax. Extensive studies in the area of genomics and transcriptomics have been conducted post flax genome sequencing. Interestingly, research has been focused more for abiotic stresses tolerance compared to disease resistance in flax through transcriptomics, while the other areas of omics such as metabolomics, proteomics, ionomics, and phenomics are in the initial stages in flax and several key questions remain unanswered. Little has been explored in the integration of omic-scale data to explain complex genetic, physiological and biochemical basis of stress tolerance in flax. In this review, the current status of various omics approaches for elucidation of molecular pathways underlying abiotic and biotic stress tolerance in flax have been presented and the importance of integrated omics technologies in future research and breeding have been emphasized to ensure sustainable yield in challenging environments.

Transcriptome Analysis of Fusarium–Tomato Interaction Based on an Updated Genome Annotation of Fusarium oxysporum f. sp. lycopersici Identifies Novel Effector Candidates That Suppress or Induce Cell Death in Nicotiana benthamiana

Fusarium oxysporum f. sp. lycopersici (Fol) causes vascular wilt disease in tomato. Upon colonization of the host, Fol secretes many small effector proteins into the xylem sap to facilitate infection. Besides known SIX (secreted in xylem) proteins, the identity of additional effectors that contribute to Fol pathogenicity remains largely unexplored. We performed a deep RNA-sequencing analysis of Fol race 2-infected tomato, used the sequence data to annotate a published genome assembly generated via PacBio SMRT sequencing of the Fol race 2 reference strain Fol4287, and analysed the resulting transcriptome to identify Fol effector candidates among the newly annotated genes. We examined the Fol-infection expression profiles of all 13 SIX genes present in Fol race 2 and identified 27 new candidate effector genes that were likewise significantly upregulated upon Fol infection. Using Agrobacterium-mediated transformation, we tested the ability of 22 of the new candidate effector genes to suppress or induce cell death in leaves of Nicotiana benthamiana. One effector candidate designated Fol-EC19, encoding a secreted guanyl-specific ribonuclease, was found to trigger cell death and two effector candidates designated Fol-EC14 and Fol-EC20, encoding a glucanase and a secreted trypsin, respectively, were identified that can suppress Bax-mediated cell death. Remarkably, Fol-EC14 and Fol-EC20 were also found to suppress I-2/Avr2- and I/Avr1-mediated cell death. Using the yeast secretion trap screening system, we showed that these three biologically-active effector candidates each contain a functional signal peptide for protein secretion. Our findings provide a basis for further understanding the virulence functions of Fol effectors.

Infection Strategies and Pathogenicity of Biotrophic Plant Fungal Pathogens

Biotrophic plant pathogenic fungi are widely distributed and are among the most damaging pathogenic organisms of agriculturally important crops responsible for significant losses in quality and yield. However, the pathogenesis of obligate parasitic pathogenic microorganisms is still under investigation because they cannot reproduce and complete their life cycle on an artificial medium. The successful lifestyle of biotrophic fungal pathogens depends on their ability to secrete effector proteins to manipulate or evade plant defense response. By integrating genomics, transcriptomics, and effectoromics, insights into how the adaptation of biotrophic plant fungal pathogens adapt to their host populations can be gained. Efficient tools to decipher the precise molecular mechanisms of rust–plant interactions, and standardized routines in genomics and functional pipelines have been established and will pave the way for comparative studies. Deciphering fungal pathogenesis not only allows us to better understand how fungal pathogens infect host plants but also provides valuable information for plant diseases control, including new strategies to prevent, delay, or inhibit fungal development. Our review provides a comprehensive overview of the efforts that have been made to decipher the effector proteins of biotrophic fungal pathogens and demonstrates how rapidly research in the field of obligate biotrophy has progressed.

The Melampsora americana Population on Salix purpurea in the Great Lakes Region Is Highly Diverse with a Contributory Influence of Clonality

Shrub willows (Salix spp.) are emerging as a viable lignocellulosic, second-generation bioenergy crop with many growth characteristics favorable for marginal lands in New York State and surrounding areas. Willow rust, caused by members of the genus Melampsora, is the most limiting disease of shrub willow in this region and remains extremely understudied. In this study, genetic diversity, genetic structure, and pathogen clonality were examined in Melampsora americana over two growing seasons via genotyping-by-sequencing to identify single-nucleotide polymorphism markers. In conjunction with this project, a reference genome of rust isolate R15-033-03 was generated to aid in variant discovery. Sampling between years allowed regional and site-specific investigation into population dynamics, in the context of both wild and cultivated hosts within high-density plantings. This work revealed that this pathogen is largely panmictic over the sampled areas, with few sites showing moderate genetic differentiation. These data support the hypothesis of sexual recombination between growing seasons because no genotype persisted across the two years of sampling. Additionally, clonality was determined as a driver of pathogen populations within cultivated fields and single shrubs; however, there is also evidence of high genetic diversity of rust isolates in all settings. This work provides a framework for M. americana population structure in the Great Lakes region, providing crucial information that can aid in future resistance breeding efforts.

The haustorium: The root of biotrophic fungal pathogens

Biotrophic plant pathogenic fungi are among the dreadful pathogens that continuously threaten the production of economically important crops. The interaction of biotrophic fungal pathogens with their hosts necessitates the development of unique infection mechanisms and involvement of various virulence-associated components. Biotrophic plant pathogenic fungi have an exceptional lifestyle that supports nutrient acquisition from cells of a living host and are fully dependent on the host for successful completion of their life cycle. The haustorium, a specialized infection structure, is the key organ for biotrophic fungal pathogens. The haustorium is not only essential in the uptake of nutrients without killing the host, but also in the secretion and delivery of effectors into the host cells to manipulate host immune system and defense responses and reprogram the metabolic flow of the host. Although there is a number of unanswered questions in this area yet, results from various studies indicate that the haustorium is the root of biotrophic fungal pathogens. This review provides an overview of current knowledge of the haustorium, its structure, composition, and functions, which includes the most recent haustorial transcriptome studies.

An automated and combinative method for the predictive ranking of candidate effector proteins of fungal plant pathogens

Fungal plant-pathogens promote infection of their hosts through the release of 'effectors'-a broad class of cytotoxic or virulence-promoting molecules. Effectors may be recognised by resistance or sensitivity receptors in the host, which can determine disease outcomes. Accurate prediction of effectors remains a major challenge in plant pathology, but if achieved will facilitate rapid improvements to host disease resistance. This study presents a novel tool and pipeline for the ranking of predicted effector candidates-Predector-which interfaces with multiple software tools and methods, aggregates disparate features that are relevant to fungal effector proteins, and applies a pairwise learning to rank approach. Predector outperformed a typical combination of secretion and effector prediction methods in terms of ranking performance when applied to a curated set of confirmed effectors derived from multiple species. We present Predector ( https://github.com/ccdmb/predector ) as a useful tool for the ranking of predicted effector candidates, which also aggregates and reports additional supporting information relevant to effector and secretome prediction in a simple, efficient, and reproducible manner.

Fatty acid dioxygenase-cytochrome P450 fusion enzymes of filamentous fungal pathogens

Oxylipins designate oxygenated unsaturated C₁₈ fatty acids. Many filamentous fungi pathogens contain dioxygenases (DOX) in oxylipin biosynthesis with homology to human cyclooxygenases. They contain a DOX domain, which is often fused to a functional cytochrome P450 at the C-terminal end. A Tyr radical in the DOX domain initiates dioxygenation of linoleic acid by hydrogen abstraction with formation of 8-, 9-, or 10-hydroperoxy metabolites. The P450 domains can catalyze heterolytic cleavage of 8- and 10-hydroperoxides with oxidation of the heme thiolate iron for hydroxylation at C-5, C-7, C-9, or C-11 and for epoxidation of the 12Z double bond; thus displaying linoleate diol synthase (LDS) and epoxy alcohol synthase (EAS) activities. LSD activities are present in the rice blast pathogen Magnaporthe oryzae, Botrytis cinerea causing grey mold and the black scurf pathogen Rhizoctonia solani. 10R-DOX-EAS has been found in M. oryzae and Fusarium oxysporum. The P450 domains may also catalyze homolytic cleavage of 8- and 9-hydroperoxy fatty acids and dehydration to produce epoxides with an adjacent double bond, i.e., allene oxides, thus displaying 8- and 9-DOX-allene oxide synthases (AOS). F. oxysporum, F. graminearum, and R. solani express 9S-DOX-AOS and Zymoseptoria tritici 8S-and 9R-DOX-AOS. Homologues are present in endemic human-pathogenic fungi with extensive studies in Aspergillus fumigatus, A. flavus (also a plant pathogen) as well as the genetic model A. nidulans. 8R-and 10R-DOX appear to bind fatty acids "headfirst" in the active site, whereas 9S-DOX binds them "tail first" in analogy with cyclooxygenases. The biological relevance of 8R-DOX-5,8-LDS (also designated PpoA) was first discovered in relation to sporulation of A. nidulans and recently for development and programmed hyphal branching of A. fumigatus. Gene deletion DOX-AOS homologues in F. verticillioides, A. flavus, and A. nidulans alters, inter alia, mycotoxin production, sporulation, and gene expression.

De novo transcriptome assembly, polymorphic SSR markers development and population genetics analyses for southern corn rust (Puccinia polysora)

Southern corn rust is a destructive maize disease caused by Puccinia polysora Underw that can lead to severe yield losses. However, genomic information and microsatellite markers are currently unavailable for this disease. In this study, we generated a total of 27,295,216 high-quality cDNA sequence reads using Illumina sequencing technology. These reads were assembled into 17,496 unigenes with an average length of 1015 bp. The functional annotation indicated that 8113 (46.37%), 1933 (11.04%) and 5516 (31.52%) unigenes showed significant similarity to known proteins in the NCBI Nr, Nt and Swiss-Prot databases, respectively. In addition, 2921 (16.70%) unigenes were assigned to KEGG database categories; 4218 (24.11%), to KOG database categories; and 6,603 (37.74%), to GO database categories. Furthermore, we identified 8,798 potential SSRs among 6653 unigenes. A total of 9 polymorphic SSR markers were developed to evaluate the genetic diversity and population structure of 96 isolates collected from Guangdong Province in China. Clonal reproduction of P. polysora in Guangdong was dominant. The YJ (Yangjiang) population had the highest genotypic diversity and the greatest number of the multilocus genotypes, followed by the HY (Heyuan), HZ (Huizhou) and XY (Xinyi) populations. These results provide valuable information for the molecular genetic analysis of P. polysora and related species.

WITHDRAWN: Fatty acid dioxygenase-cytochrome P450 fusion enzymes of the top 10 fungal pathogens in molecular plant pathology and human-pathogenic fungi

This article has been withdrawn at the request of the author(s) and/or editor. The Publisher apologizes for any inconvenience this may cause. The full Elsevier Policy on Article Withdrawal can be found at https://www.elsevier.com/about/our-business/policies/article-withdrawal

Distribution of methionine sulfoxide reductases in fungi and conservation of the free-methionine-R-sulfoxide reductase in multicellular eukaryotes

Methionine, either as a free amino acid or included in proteins, can be oxidized into methionine sulfoxide (MetO), which exists as R and S diastereomers. Almost all characterized organisms possess thiol-oxidoreductases named methionine sulfoxide reductase (Msr) enzymes to reduce MetO back to Met. MsrA and MsrB reduce the S and R diastereomers of MetO, respectively, with strict stereospecificity and are found in almost all organisms. Another type of thiol-oxidoreductase, the free-methionine-R-sulfoxide reductase (fRMsr), identified so far in prokaryotes and a few unicellular eukaryotes, reduces the R MetO diastereomer of the free amino acid. Moreover, some bacteria possess molybdenum-containing enzymes that reduce MetO, either in the free or protein-bound forms. All these Msrs play important roles in the protection of organisms against oxidative stress. Fungi are heterotrophic eukaryotes that colonize all niches on Earth and play fundamental functions, in organic matter recycling, as symbionts, or as pathogens of numerous organisms. However, our knowledge on fungal Msrs is still limited. Here, we performed a survey of msr genes in almost 700 genomes across the fungal kingdom. We show that most fungi possess one gene coding for each type of methionine sulfoxide reductase: MsrA, MsrB, and fRMsr. However, several fungi living in anaerobic environments or as obligate intracellular parasites were devoid of msr genes. Sequence inspection and phylogenetic analyses allowed us to identify non-canonical sequences with potentially novel enzymatic properties. Finaly, we identified several ocurences of msr horizontal gene transfer from bacteria to fungi.

Identification of two compounds able to improve flax resistance towards Fusarium oxysporum infection

Plants are surrounded by a diverse range of microorganisms that causes serious crop losses and requires the use of pesticides. Flax is a major crop in Normandy used for its fibres and is regularly challenged by the pathogenic fungus Fusarium oxysporum (Fo) f. sp. lini. To protect themselves, plants use "innate immunity" as a first line of defense level against pathogens. Activation of plant defense with elicitors could be an alternative for crop plant protection. A previous work was conducted by screening a chemical library and led to the identification of compounds able to activate defense responses in Arabidopsis thaliana. Four compounds were tested for their abilities to improve resistance of two flax varieties against Fo. Two of them, one natural (holaphyllamine or HPA) and one synthetic (M4), neither affected flax nor Fo growth. HPA and M4 induced oxidative burst and callose deposition. Furthermore, HPA and M4 caused changes in the expression patterns of defense-related genes coding a glucanase and a chitinase-like. Finally, plants pre-treated with HPA or M4 exhibited a significant decrease in the disease symptoms. Together, these findings demonstrate that HPA and M4 are able to activate defense responses in flax and improve its resistance against Fo infection.

Evolution of isoprenyl diphosphate synthase-like terpene synthases in fungi

Terpene synthases (TPSs) and trans-isoprenyl diphosphate synthases (IDSs) are among the core enzymes for creating the enormous diversity of terpenoids. Despite having no sequence homology, TPSs and IDSs share a conserved “α terpenoid synthase fold” and a trinuclear metal cluster for catalysis, implying a common ancestry with TPSs hypothesized to evolve from IDSs anciently. Here we report on the identification and functional characterization of novel IDS-like TPSs (ILTPSs) in fungi that evolved from IDS relatively recently, indicating recurrent evolution of TPSs from IDSs. Through large-scale bioinformatic analyses of fungal IDSs, putative ILTPSs that belong to the geranylgeranyl diphosphate synthase (GGDPS) family of IDSs were identified in three species of Melampsora. Among the GGDPS family of the two Melampsora species experimentally characterized, one enzyme was verified to be bona fide GGDPS and all others were demonstrated to function as TPSs. Melampsora ILTPSs displayed kinetic parameters similar to those of classic TPSs. Key residues underlying the determination of GGDPS versus ILTPS activity and functional divergence of ILTPSs were identified. Phylogenetic analysis implies a recent origination of these ILTPSs from a GGDPS progenitor in fungi, after the split of Melampsora from other genera within the class of Pucciniomycetes. For the poplar leaf rust fungus Melampsora larici-populina, the transcripts of its ILTPS genes were detected in infected poplar leaves, suggesting possible involvement of these recently evolved ILTPS genes in the infection process. This study reveals the recurrent evolution of TPSs from IDSs since their ancient occurrence and points to the possibility of a wide distribution of ILTPS genes in three domains of life.

High-Quality Genome Assembly of Fusarium oxysporum f. sp. lini

In the present work, a highly pathogenic isolate of Fusarium oxysporum f. sp. lini, which is the most harmful pathogen affecting flax (Linum usitatissimum L.), was sequenced for the first time. To achieve a high-quality genome assembly, we used the combination of two sequencing platforms - Oxford Nanopore Technologies (MinION system) with long noisy reads and Illumina (HiSeq 2500 instrument) with short accurate reads. Given the quality of DNA is crucial for Nanopore sequencing, we developed the protocol for extraction of pure high-molecular-weight DNA from fungi. Sequencing of DNA extracted using this protocol allowed us to obtain about 85x genome coverage with long (N50 = 29 kb) MinION reads and 30x coverage with 2 × 250 bp HiSeq reads. Several tools were developed for genome assembly; however, they provide different results depending on genome complexity, sequencing data volume, read length and quality. We benchmarked the most requested assemblers (Canu, Flye, Shasta, wtdbg2, and MaSuRCA), Nanopore polishers (Medaka and Racon), and Illumina polishers (Pilon and POLCA) on our sequencing data. The assembly performed with Canu and polished with Medaka and POLCA was considered the most full and accurate. After further elimination of redundant contigs using Purge Haplotigs, we achieved a high-quality genome of F. oxysporum f. sp. lini with a total length of 59 Mb, N50 of 3.3 Mb, and 99.5% completeness according to BUSCO. We also obtained a complete circular mitochondrial genome with a length of 38.7 kb. The achieved assembly expands studies on F. oxysporum and plant-pathogen interaction in flax.

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.

Evolution of virulence in rust fungi - multiple solutions to one problem

Rust fungi are major pathogens that negatively affect crops and ecosystems. Recent rust disease epidemics driven by the emergence of strains with novel virulence profiles demand a better understanding of the evolutionary mechanisms of these organisms. Here, we review research advances in genome-scale analysis coupled with functional validation of effector candidate genes that have been instrumental to elucidate processes that contribute to changes in virulence phenotypes. We highlight how haplotype-phased genome references have paved the road to link these processes to the reproductive phases of rust fungi and have provided evidence for somatic exchange between strains as an important mechanism for generating diversity in asexual populations. With increasing data availability, we envision the future development of molecular virulence diagnostic tools.

Flax rust infection transcriptomics reveals a transcriptional profile that may be indicative for rust Avr genes

Secreted effectors of fungal pathogens are essential elements for disease development. However, lack of sequence conservation among identified effectors has long been a problem for predicting effector complements in fungi. Here we have explored the expression characteristics of avirulence (Avr) genes and candidate effectors of the flax rust fungus, Melampsora lini. We performed transcriptome sequencing and real-time quantitative PCR (qPCR) on RNA extracted from ungerminated spores, germinated spores, isolated haustoria and flax seedlings inoculated with M. lini isolate CH5 during plant infection. Genes encoding two categories of M. lini proteins, namely Avr proteins and plant cell wall degrading enzymes (CWDEs), were investigated in detail. Analysis of the expression profiles of 623 genes encoding predicted secreted proteins in the M. lini transcriptome shows that the six known Avr genes (i.e. AvrM (avrM), AvrM14, AvrL2, AvrL567, AvrP123 (AvrP) and AvrP4) fall within a group of 64 similarly expressed genes that are induced in planta and show a peak of expression early in infection with a subsequent decline towards sporulation. Other genes within this group include two paralogues of AvrL2, an AvrL567 virulence allele, and a number of genes encoding putative effector proteins. By contrast, M. lini genes encoding CWDEs fall into different expression clusters with their distribution often unrelated to their catalytic activity or substrate targets. These results suggest that synthesis of M. lini Avr proteins may be regulated in a coordinated fashion and that the expression profiling-based analysis has significant predictive power for the identification of candidate Avr genes.

Structural genomics applied to the rust fungus Melampsora larici-populina reveals two candidate effector proteins adopting cystine knot and NTF2-like protein folds

Rust fungi are plant pathogens that secrete an arsenal of effector proteins interfering with plant functions and promoting parasitic infection. Effectors are often species-specific, evolve rapidly, and display low sequence similarities with known proteins. How rust fungal effectors function in host cells remains elusive, and biochemical and structural approaches have been scarcely used to tackle this question. In this study, we produced recombinant proteins of eleven candidate effectors of the leaf rust fungus Melampsora larici-populina in Escherichia coli. We successfully purified and solved the three-dimensional structure of two proteins, MLP124266 and MLP124017, using NMR spectroscopy. Although both MLP124266 and MLP124017 show no sequence similarity with known proteins, they exhibit structural similarities to knottins, which are disulfide-rich small proteins characterized by intricate disulfide bridges, and to nuclear transport factor 2-like proteins, which are molecular containers involved in a wide range of functions, respectively. Interestingly, such structural folds have not been reported so far in pathogen effectors, indicating that MLP124266 and MLP124017 may bear novel functions related to pathogenicity. Our findings show that sequence-unrelated effectors can adopt folds similar to known proteins, and encourage the use of biochemical and structural approaches to functionally characterize effector candidates.

De Novo Genome Assembly and Comparative Genomics of the Barley Leaf Rust Pathogen Puccinia hordei Identifies Candidates for Three Avirulence Genes

Puccinia hordei (Ph) is a damaging pathogen of barley throughout the world. Despite its importance, almost nothing is known about the genomics of this pathogen, and a reference genome is lacking. In this study, the first reference genome was assembled for an Australian isolate of Ph ("Ph560") using long-read SMRT sequencing. A total of 838 contigs were assembled, with a total size of 207 Mbp. This included both haplotype collapsed and separated regions, consistent with an estimated haploid genome size of about 150Mbp. An annotation pipeline that combined RNA-Seq of Ph-infected host tissues and homology to proteins from four other Puccinia species predicted 25,543 gene models of which 1,450 genes were classified as encoding secreted proteins based on the prediction of a signal peptide and no transmembrane domain. Genome resequencing using short-read technology was conducted for four additional Australian strains, Ph612, Ph626, Ph608 and Ph584, which are considered to be simple mutational derivatives of Ph560 with added virulence to one or two of three barley leaf rust resistance genes (viz. Rph3, Rph13 and Rph19). To identify candidate genes for the corresponding avirulence genes AvrRph3, AvrRph13 and AvrRph19, genetic variation in predicted secreted protein genes between the strains was correlated to the virulence profiles of each, identifying 35, 29 and 46 candidates for AvrRph13, AvrRph3 and AvrRph19, respectively. The identification of these candidate genes provides a strong foundation for future efforts to isolate these three avirulence genes, investigate their biological properties, and develop new diagnostic tests for monitoring pathogen virulence.

Chromium sequencing: the doors open for genomics of obligate plant pathogens

It is challenging to sequence and assemble genomes of obligate plant pathogens and microorganisms because of limited amounts of DNA, comparatively large genomes and high numbers of repeat regions. We sequenced the 1.2 gigabase genome of an obligate rust fungus, Austropuccinia psidii, the cause of rust on Myrtaceae, with a Chromium 10X library. This technology has mostly been applied for single-cell sequencing in immunological studies of mammals. We compared scaffolds of a genome assembled from the Chromium library with one assembled from combined paired-end and mate-pair libraries, sequenced with Illumina HiSeq. Chromium 10X provided a superior assembly, in terms of number of scaffolds, N50 and number of genes recovered. It required less DNA than other methods and was sequenced and assembled at a lower cost. Chromium sequencing could provide a solution to sequence and assemble genomes of obligate plant pathogens where the amount of available DNA is a limiting factor.

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.

Genome sequencing and transcript analysis of Hemileia vastatrix reveal expression dynamics of candidate effectors dependent on host compatibility

Coffee leaf rust caused by the fungus Hemileia vastatrix is one of the most important leaf diseases of coffee plantations worldwide. Current knowledge of the H. vastatrix genome is limited and only a small fraction of the total fungal secretome has been identified. In order to obtain a more comprehensive understanding of its secretome, we aimed to sequence and assemble the entire H. vastatrix genome using two next-generation sequencing platforms and a hybrid assembly strategy. This resulted in a 547 Mb genome of H. vastatrix race XXXIII (Hv33), with 13,364 predicted genes that encode 13,034 putative proteins with transcriptomic support. Based on this proteome, 615 proteins contain putative secretion peptides, and lack transmembrane domains. From this putative secretome, 111 proteins were identified as candidate effectors (EHv33) unique to H. vastatrix, and a subset consisting of 17 EHv33 genes was selected for a temporal gene expression analysis during infection. Five genes were significantly induced early during an incompatible interaction, indicating their potential role as pre-haustorial effectors possibly recognized by the resistant coffee genotype. Another nine genes were significantly induced after haustorium formation in the compatible interaction. Overall, we suggest that this fungus is able to selectively mount its survival strategy with effectors that depend on the host genotype involved in the infection process.

Apollo: Democratizing genome annotation

Genome annotation is the process of identifying the location and function of a genome's encoded features. Improving the biological accuracy of annotation is a complex and iterative process requiring researchers to review and incorporate multiple sources of information such as transcriptome alignments, predictive models based on sequence profiles, and comparisons to features found in related organisms. Because rapidly decreasing costs are enabling an ever-growing number of scientists to incorporate sequencing as a routine laboratory technique, there is widespread demand for tools that can assist in the deliberative analytical review of genomic information. To this end, we present Apollo, an open source software package that enables researchers to efficiently inspect and refine the precise structure and role of genomic features in a graphical browser-based platform. Some of Apollo’s newer user interface features include support for real-time collaboration, allowing distributed users to simultaneously edit the same encoded features while also instantly seeing the updates made by other researchers on the same region in a manner similar to Google Docs. Its technical architecture enables Apollo to be integrated into multiple existing genomic analysis pipelines and heterogeneous laboratory workflow platforms. Finally, we consider the implications that Apollo and related applications may have on how the results of genome research are published and made accessible.

Structural and functional insights into the modulation of the activity of a flax cytokinin oxidase by flax rust effector AvrL567-A

During infection, plant pathogens secrete effector proteins to facilitate colonization. In comparison with our knowledge of bacterial effectors, the current understanding of how fungal effectors function is limited. In this study, we show that the effector AvrL567-A from the flax rust fungus Melampsora lini interacts with a flax cytosolic cytokinin oxidase, LuCKX1.1, using both yeast two-hybrid and in planta bimolecular fluorescence assays. Purified LuCKX1.1 protein shows catalytic activity against both N6-(Δ2-isopentenyl)-adenine (2iP) and trans-zeatin (tZ) substrates. Incubation of LuCKX1.1 with AvrL567-A results in increased catalytic activity against both substrates. The crystal structure of LuCKX1.1 and docking studies with AvrL567-A indicate that the AvrL567 binding site involves a flexible surface-exposed region that surrounds the cytokinin substrate access site, which may explain its effect in modulating LuCKX1.1 activity. Expression of AvrL567-A in transgenic flax plants gave rise to an epinastic leaf phenotype consistent with hormonal effects, although no difference in overall cytokinin levels was observed. We propose that, during infection, plant pathogens may differentially modify the levels of extracellular and intracellular cytokinins.

Candidate Effectors From Uromyces appendiculatus, the Causal Agent of Rust on Common Bean, Can Be Discriminated Based on Suppression of Immune Responses

Rust fungi are devastating pathogens for several important crop plants. The biotrophic lifestyle of rust fungi requires that they influence their host plants to create a favorable environment for growth and reproduction. Rust fungi secrete a variety of effector proteins that manipulate host target proteins to alter plant metabolism and suppress defense responses. Because of the obligate biotrophic lifestyle of rust fungi, direct evidence for effector function is difficult to obtain, and so suites of experiments utilizing expression in heterologous systems are necessary. Here, we present results from a yeast cell death suppression assay and assays for suppression of PAMP-triggered immunity (PTI) and effector triggered immunity (ETI) based on delivery of effectors through the bacterial type III secretion system. In addition, subcellular localization was tested using transient expression of GFP fusion proteins in Nicotiana benthamiana through Agrobacterium infiltration. We tested 31 representative effector candidates from the devastating common bean rust pathogen Uromyces appendiculatus. These effector candidates were selected based on features of their gene families, most important lineage specificity. We show that several of our effector candidates suppress plant defense. Some of them also belong to families of effector candidates that are present in multiple rust species where their homologs probably also have effector functions. In our analysis of candidate effector mRNA expression, some of those effector candidates that gave positive results in the other assays were not up-regulated during plant infection, indicating that either these proteins have functions at multiple life stages or that strong up-regulation of RNA level in planta may not be as important a criterion for identifying effectors as previously thought. Overall, our pipeline for selecting effector candidates based on sequence features followed by screening assays using heterologous expression systems was successful in discriminating effector candidates. This work lays the foundation for functional characterization of U. appendiculatus effectors, the identification of effector targets, and identification of novel sources for resistance in common bean.

Arms race: diverse effector proteins with conserved motifs

Effector proteins play important roles in the infection by pathogenic oomycetes and fungi or the colonization by endophytic and mycorrhizal fungi. They are either translocated into the host plant cells via specific translocation mechanisms and function in the host's cytoplasm or nucleus, or they reside in the apoplast of the plant cells and act at the extracellular host-microbe interface. Many effector proteins possess conserved motifs (such as the RXLR, CRN, LysM, RGD, DELD, EAR, RYWT, Y/F/WXC or CFEM motifs) localized in their N- or C-terminal regions. Analysis of the functions of effector proteins, especially so-called "core effectors", is crucial for the understanding of pathogenicity/symbiosis mechanisms and plant defense strategies, and helps to develop breeding strategies for pathogen-resistant cultivars, and to increase crop yield and quality as well as abiotic stress resistance. This review summarizes current knowledge about these effector proteins with the conversed motifs and their involvement in pathogenic or mutualistic plant/fungal interactions.

Detection of Race-Specific Resistance Against Puccinia coronata f. sp. avenae in Brachypodium Species

Oat crown rust caused by Puccinia coronata f. sp. avenae is the most destructive foliar disease of cultivated oat. Characterization of genetic factors controlling resistance responses to Puccinia coronata f. sp. avenae in nonhost species could provide new resources for developing disease protection strategies in oat. We examined symptom development and fungal colonization levels of a collection of Brachypodium distachyon and B. hybridum accessions infected with three North American P. coronata f. sp. avenae isolates. Our results demonstrated that colonization phenotypes are dependent on both host and pathogen genotypes, indicating a role for race-specific responses in these interactions. These responses were independent of the accumulation of reactive oxygen species. Expression analysis of several defense-related genes suggested that salicylic acid and ethylene-mediated signaling but not jasmonic acid are components of resistance reaction to P. coronata f. sp. avenae. Our findings provide the basis to conduct a genetic inheritance study to examine whether effector-triggered immunity contributes to nonhost resistance to P. coronata f. sp. avenae in Brachypodium spp.

New insights into Phakopsora pachyrhizi infection based on transcriptome analysis in planta

Asian soybean rust (ASR) is one of the most destructive diseases affecting soybeans. The causative agent of ASR, the fungus Phakopsora pachyrhizi, presents characteristics that make it difficult to study in vitro, limiting our knowledge of plant-pathogen dynamics. Therefore, this work used leaf lesion laser microdissection associated with deep sequencing to determine the pathogen transcriptome during compatible and incompatible interactions with soybean. The 36,350 generated unisequences provided an overview of the main genes and biological pathways that were active in the fungus during the infection cycle. We also identified the most expressed transcripts, including sequences similar to other fungal virulence and signaling proteins. Enriched P. pachyrhizi transcripts in the resistant (PI561356) soybean genotype were related to extracellular matrix organization and metabolic signaling pathways and, among infection structures, in amino acid metabolism and intracellular transport. Unisequences were further grouped into gene families along predicted sequences from 15 other fungi and oomycetes, including rust fungi, allowing the identification of conserved multigenic families, as well as being specific to P. pachyrhizi. The results revealed important biological processes observed in P. pachyrhizi, contributing with information related to fungal biology and, consequently, a better understanding of ASR.

Improved prediction of fungal effector proteins from secretomes with EffectorP 2.0

Plant-pathogenic fungi secrete effector proteins to facilitate infection. We describe extensive improvements to EffectorP, the first machine learning classifier for fungal effector prediction. EffectorP 2.0 is now trained on a larger set of effectors and utilizes a different approach based on an ensemble of classifiers trained on different subsets of negative data, offering different views on classification. EffectorP 2.0 achieves an accuracy of 89%, compared with 82% for EffectorP 1.0 and 59.8% for a small size classifier. Important features for effector prediction appear to be protein size, protein net charge as well as the amino acids serine and cysteine. EffectorP 2.0 decreases the number of predicted effectors in secretomes of fungal plant symbionts and saprophytes by 40% when compared with EffectorP 1.0. However, EffectorP 1.0 retains value, and combining EffectorP 1.0 and 2.0 results in a stringent classifier with a low false positive rate of 9%. EffectorP 2.0 predicts significant enrichments of effectors in 12 of 13 sets of infection-induced proteins from diverse fungal pathogens, whereas a small cysteine-rich classifier detects enrichment in only seven of 13. EffectorP 2.0 will fast track the prioritization of high-confidence effector candidates for functional validation and aid in improving our understanding of effector biology. EffectorP 2.0 is available at http://effectorp.csiro.au.

Annotation resource of tandem repeat-containing secretory proteins in sixty fungi

Fungal secretory proteins that interact with host plants are regarded as effectors. Because fungal effectors rarely contain conserved sequence features, identification and annotation of fungal effectors from predicted secretory proteins are difficult using outward comparison methods such as BLAST or hidden Markov model. In desire of more sequence features to prioritize research interests of fungal secretory proteins, this study developed a pipeline to identify tandem repeat (TR) domain within putative secretory proteins and tested a hypothesis that at least one type of TR domain in non-orthologous secretory proteins has emerged from convergent evolution for plant pathogenicity. There were 2804 types of TR domains and a total of 2925 TR-containing secretory proteins found from 60 fungi. There was no conserved type of TR domain shared only by plant pathogens, indicating functional divergence for different types of TR domain and TR-containing secretory proteins. The annotation resource of putative fungal TR-containing secretory proteins provides new sequence features that will be useful for the community interested in fungal effector biology.

A Near-Complete Haplotype-Phased Genome of the Dikaryotic Wheat Stripe Rust Fungus Puccinia striiformis f. sp. tritici Reveals High Interhaplotype Diversity

A long-standing biological question is how evolution has shaped the genomic architecture of dikaryotic fungi. To answer this, high-quality genomic resources that enable haplotype comparisons are essential. Short-read genome assemblies for dikaryotic fungi are highly fragmented and lack haplotype-specific information due to the high heterozygosity and repeat content of these genomes. Here, we present a diploid-aware assembly of the wheat stripe rust fungus Puccinia striiformis f. sp. tritici based on long reads using the FALCON-Unzip assembler. Transcriptome sequencing data sets were used to infer high-quality gene models and identify virulence genes involved in plant infection referred to as effectors. This represents the most complete Puccinia striiformis f. sp. tritici genome assembly to date (83 Mb, 156 contigs, N50 of 1.5 Mb) and provides phased haplotype information for over 92% of the genome. Comparisons of the phase blocks revealed high interhaplotype diversity of over 6%. More than 25% of all genes lack a clear allelic counterpart. When we investigated genome features that potentially promote the rapid evolution of virulence, we found that candidate effector genes are spatially associated with conserved genes commonly found in basidiomycetes. Yet, candidate effectors that lack an allelic counterpart are more distant from conserved genes than allelic candidate effectors and are less likely to be evolutionarily conserved within the P. striiformis species complex and Pucciniales In summary, this haplotype-phased assembly enabled us to discover novel genome features of a dikaryotic plant-pathogenic fungus previously hidden in collapsed and fragmented genome assemblies.IMPORTANCE Current representations of eukaryotic microbial genomes are haploid, hiding the genomic diversity intrinsic to diploid and polyploid life forms. This hidden diversity contributes to the organism's evolutionary potential and ability to adapt to stress conditions. Yet, it is challenging to provide haplotype-specific information at a whole-genome level. Here, we take advantage of long-read DNA sequencing technology and a tailored-assembly algorithm to disentangle the two haploid genomes of a dikaryotic pathogenic wheat rust fungus. The two genomes display high levels of nucleotide and structural variations, which lead to allelic variation and the presence of genes lacking allelic counterparts. Nonallelic candidate effector genes, which likely encode important pathogenicity factors, display distinct genome localization patterns and are less likely to be evolutionary conserved than those which are present as allelic pairs. This genomic diversity may promote rapid host adaptation and/or be related to the age of the sequenced isolate since last meiosis.

De Novo Assembly and Phasing of Dikaryotic Genomes from Two Isolates of Puccinia coronata f. sp. avenae, the Causal Agent of Oat Crown Rust

Oat crown rust, caused by the fungus Pucinnia coronata f. sp. avenae, is a devastating disease that impacts worldwide oat production. For much of its life cycle, P. coronata f. sp. avenae is dikaryotic, with two separate haploid nuclei that may vary in virulence genotype, highlighting the importance of understanding haplotype diversity in this species. We generated highly contiguous de novo genome assemblies of two P. coronata f. sp. avenae isolates, 12SD80 and 12NC29, from long-read sequences. In total, we assembled 603 primary contigs for 12SD80, for a total assembly length of 99.16 Mbp, and 777 primary contigs for 12NC29, for a total length of 105.25 Mbp; approximately 52% of each genome was assembled into alternate haplotypes. This revealed structural variation between haplotypes in each isolate equivalent to more than 2% of the genome size, in addition to about 260,000 and 380,000 heterozygous single-nucleotide polymorphisms in 12SD80 and 12NC29, respectively. Transcript-based annotation identified 26,796 and 28,801 coding sequences for isolates 12SD80 and 12NC29, respectively, including about 7,000 allele pairs in haplotype-phased regions. Furthermore, expression profiling revealed clusters of coexpressed secreted effector candidates, and the majority of orthologous effectors between isolates showed conservation of expression patterns. However, a small subset of orthologs showed divergence in expression, which may contribute to differences in virulence between 12SD80 and 12NC29. This study provides the first haplotype-phased reference genome for a dikaryotic rust fungus as a foundation for future studies into virulence mechanisms in P. coronata f. sp. avenaeIMPORTANCE Disease management strategies for oat crown rust are challenged by the rapid evolution of Puccinia coronata f. sp. avenae, which renders resistance genes in oat varieties ineffective. Despite the economic importance of understanding P. coronata f. sp. avenae, resources to study the molecular mechanisms underpinning pathogenicity and the emergence of new virulence traits are lacking. Such limitations are partly due to the obligate biotrophic lifestyle of P. coronata f. sp. avenae as well as the dikaryotic nature of the genome, features that are also shared with other important rust pathogens. This study reports the first release of a haplotype-phased genome assembly for a dikaryotic fungal species and demonstrates the amenability of using emerging technologies to investigate genetic diversity in populations of P. coronata f. sp. avenae.

Suppression or Activation of Immune Responses by Predicted Secreted Proteins of the Soybean Rust Pathogen Phakopsora pachyrhizi

Rust fungi, such as the soybean rust pathogen Phakopsora pachyrhizi, are major threats to crop production. They form specialized haustoria that are hyphal structures intimately associated with host-plant cell membranes. These haustoria have roles in acquiring nutrients and secreting effector proteins that manipulate host immune systems. Functional characterization of effector proteins of rust fungi is important for understanding mechanisms that underlie their virulence and pathogenicity. Hundreds of candidate effector proteins have been predicted for rust pathogens, but it is not clear how to prioritize these effector candidates for further characterization. There is a need for high-throughput approaches for screening effector candidates to obtain experimental evidence for effector-like functions, such as the manipulation of host immune systems. We have focused on identifying effector candidates with immune-related functions in the soybean rust fungus P. pachyrhizi. To facilitate the screening of many P. pachyrhizi effector candidates (named PpECs), we used heterologous expression systems, including the bacterial type III secretion system, Agrobacterium infiltration, a plant virus, and a yeast strain, to establish an experimental pipeline for identifying PpECs with immune-related functions and establishing their subcellular localizations. Several PpECs were identified that could suppress or activate immune responses in nonhost Nicotiana benthamiana, N. tabacum, Arabidopsis, tomato, or pepper plants.

Identification of Plasmodiophora brassicae effectors - A challenging goal

Clubroot is an economically important disease affecting Brassica plants worldwide. Plasmodiophora brassicae is the protist pathogen associated with the disease, and its soil-borne obligate parasitic nature has impeded studies related to its biology and the mechanisms involved in its infection of the plant host. The identification of effector proteins is key to understanding how the pathogen manipulates the plant's immune response and the genes involved in resistance. After more than 140 years studying clubroot and P. brassicae, very little is known about the effectors playing key roles in the infection process and subsequent disease progression. Here we analyze the information available for identified effectors and suggest several features of effector genes that can be used in the search for others. Based on the information presented in this review, we propose a comprehensive bioinformatics pipeline for effector identification and provide a list of the bioinformatics tools available for such.

Secretome Characterization and Correlation Analysis Reveal Putative Pathogenicity Mechanisms and Identify Candidate Avirulence Genes in the Wheat Stripe Rust Fungus Puccinia striiformis f. sp. tritici

Stripe (yellow) rust, caused by Puccinia striiformis f. sp. tritici (Pst), is one of the most destructive diseases of wheat worldwide. Planting resistant cultivars is an effective way to control this disease, but race-specific resistance can be overcome quickly due to the rapid evolving Pst population. Studying the pathogenicity mechanisms is critical for understanding how Pst virulence changes and how to develop wheat cultivars with durable resistance to stripe rust. We re-sequenced 7 Pst isolates and included additional 7 previously sequenced isolates to represent balanced virulence/avirulence profiles for several avirulence loci in seretome analyses. We observed an uneven distribution of heterozygosity among the isolates. Secretome comparison of Pst with other rust fungi identified a large portion of species-specific secreted proteins, suggesting that they may have specific roles when interacting with the wheat host. Thirty-two effectors of Pst were identified from its secretome. We identified candidates for Avr genes corresponding to six Yr genes by correlating polymorphisms for effector genes to the virulence/avirulence profiles of the 14 Pst isolates. The putative AvYr76 was present in the avirulent isolates, but absent in the virulent isolates, suggesting that deleting the coding region of the candidate avirulence gene has produced races virulent to resistance gene Yr76. We conclude that incorporating avirulence/virulence phenotypes into correlation analysis with variations in genomic structure and secretome, particularly presence/absence polymorphisms of effectors, is an efficient way to identify candidate Avr genes in Pst. The candidate effector genes provide a rich resource for further studies to determine the evolutionary history of Pst populations and the co-evolutionary arms race between Pst and wheat. The Avr candidates identified in this study will lead to cloning avirulence genes in Pst, which will enable us to understand molecular mechanisms underlying Pst-wheat interactions, to determine the effectiveness of resistance genes and further to develop durable resistance to stripe rust.

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.

Fungal phytopathogens encode functional homologues of plant rapid alkalinization factor (RALF) peptides

In this article, we describe the presence of genes encoding close homologues of an endogenous plant peptide, rapid alkalinization factor (RALF), within the genomes of 26 species of phytopathogenic fungi. Members of the RALF family are key growth factors in plants, and the sequence of the RALF active region is well conserved between plant and fungal proteins. RALF1-like sequences were observed in most cases; however, RALF27-like sequences were present in the Sphaerulina musiva and Septoria populicola genomes. These two species are pathogens of poplar and, interestingly, the closest relative to their respective RALF genes is a poplar RALF27-like sequence. RALF peptides control cellular expansion during plant development, but were originally defined on the basis of their ability to induce rapid alkalinization in tobacco cell cultures. To test whether the fungal RALF peptides were biologically active in plants, we synthesized RALF peptides corresponding to those encoded by two sequenced genomes of the tomato pathogen Fusarium oxysporum f. sp. lycopersici. One of these peptides inhibited the growth of tomato seedlings and elicited responses in tomato and Nicotiana benthamiana typical of endogenous plant RALF peptides (reactive oxygen species burst, induced alkalinization and mitogen-activated protein kinase activation). Gene expression analysis confirmed that a RALF-encoding gene in F. oxysporum f. sp. lycopersici was expressed during infection on tomato. However, a subsequent reverse genetics approach revealed that the RALF peptide was not required by F. oxysporum f. sp. lycopersici for infection on tomato roots. This study has demonstrated the presence of functionally active RALF peptides encoded within phytopathogens that harbour an as yet undetermined role in plant-pathogen interactions.

Effectors of Filamentous Plant Pathogens: Commonalities amid Diversity

Fungi and oomycetes are filamentous microorganisms that include a diversity of highly developed pathogens of plants. These are sophisticated modulators of plant processes that secrete an arsenal of effector proteins to target multiple host cell compartments and enable parasitic infection. Genome sequencing revealed complex catalogues of effectors of filamentous pathogens, with some species harboring hundreds of effector genes. Although a large fraction of these effector genes encode secreted proteins with weak or no sequence similarity to known proteins, structural studies have revealed unexpected similarities amid the diversity. This article reviews progress in our understanding of effector structure and function in light of these new insights. We conclude that there is emerging evidence for multiple pathways of evolution of effectors of filamentous plant pathogens but that some families have probably expanded from a common ancestor by duplication and diversification. Conserved folds, such as the oomycete WY and the fungal MAX domains, are not predictive of the precise function of the effectors but serve as a chassis to support protein structural integrity while providing enough plasticity for the effectors to bind different host proteins and evolve unrelated activities inside host cells. Further effector evolution and diversification arise via short linear motifs, domain integration and duplications, and oligomerization.

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.

The complete genome sequence of the phytopathogenic fungus Sclerotinia sclerotiorum reveals insights into the genome architecture of broad host range pathogens

Sclerotinia sclerotiorum is a phytopathogenic fungus with over 400 hosts including numerous economically important cultivated species. This contrasts many economically destructive pathogens that only exhibit a single or very few hosts. Many plant pathogens exhibit a “two-speed” genome. So described because their genomes contain alternating gene rich, repeat sparse and gene poor, repeat-rich regions. In fungi, the repeat-rich regions may be subjected to a process termed repeat-induced point mutation (RIP). Both repeat activity and RIP are thought to play a significant role in evolution of secreted virulence proteins, termed effectors. We present a complete genome sequence of S. sclerotiorum generated using Single Molecule Real-Time Sequencing technology with highly accurate annotations produced using an extensive RNA sequencing data set. We identified 70 effector candidates and have highlighted their in planta expression profiles. Furthermore, we characterized the genome architecture of S. sclerotiorum in comparison to plant pathogens that exhibit “two-speed” genomes. We show that there is a significant association between positions of secreted proteins and regions with a high RIP index in S. sclerotiorum but we did not detect a correlation between secreted protein proportion and GC content. Neither did we detect a negative correlation between CDS content and secreted protein proportion across the S. sclerotiorum genome. We conclude that S. sclerotiorum exhibits subtle signatures of enhanced mutation of secreted proteins in specific genomic compartments as a result of transposition and RIP activity. However, these signatures are not observable at the whole-genome scale.

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.

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.

Biology, dynamics, and applications of transposable elements in basidiomycete fungi

The phylum Basidiomycota includes filamentous fungi and yeast species with different ecological and genomic characteristics. Transposable elements (TEs) are abundant components of most eukaryotic genomes, and their transition from being genomic parasites to key drivers of genomic architecture, functionality, and evolution is a subject receiving much attention. In light of the abundant genomic information released during the last decade, the aims of this mini-review are to discuss the dynamics and impact of TEs in basidiomycete fungi. To do this, we surveyed and explored data from 75 genomes, which encompass the phylogenetic diversity of the phylum Basidiomycota. We describe annotation approaches and analyze TE distribution in the context of species phylogeny and genome size. Further, we review the most relevant literature about the role of TEs in species lifestyle, their impact on genome architecture and functionality, and the defense mechanisms evolved to control their proliferation. Finally, we discuss potential applications of TEs that can drive future innovations in fungal research.