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.

Purification of fungal haustoria from infected plant tissue by flow cytometry.

A hallmark of biotrophy in many fungal plant pathogens is the formation of a specialized pathogenic structure called the haustorium from infectious hyphae. This is the major parasitic structure, where nutrients are taken up from the host and pathogenicity factors are exported to the host tissue. Obligate biotrophic fungi can typically be cultured in vivo only to a limited extent and do not produce haustoria under these conditions. This has hampered the application of classic molecular biology techniques to haustoria-forming pathogens. The lectin Concanavalin A (Con A), which binds specifically to sugars present on the exterior of rust haustoria, was first used in a column-based affinity purification procedure in 1992 (Hahn and Mendgen, Protoplasma 170:95-103, 1992). Here we describe a new technique where we combine initial gradient purification of haustoria with flow-sorting based on labeling of haustoria with fluorescent Con A. Our method allows haustorial isolation with purity above 98 % and yields ten times more isolated haustoria in a single experiment than the previous procedure.

Comparative genomics of Australian isolates of the wheat stem rust pathogen Puccinia graminis f. sp. tritici reveals extensive polymorphism in candidate effector genes.

The wheat stem rust fungus Puccinia graminis f. sp. tritici (Pgt) is one of the most destructive pathogens of wheat. In this study, a draft genome was built for a founder Australian Pgt isolate of pathotype (pt.) 21-0 (collected in 1954) by next generation DNA sequencing. A combination of reference-based assembly using the genome of the previously sequenced American Pgt isolate CDL 75-36-700-3 (p7a) and de novo assembly were performed resulting in a 92 Mbp reference genome for Pgt isolate 21-0. Approximately 13 Mbp of de novo assembled sequence in this genome is not present in the p7a reference assembly. This novel sequence is not specific to 21-0 as it is also present in three other Pgt rust isolates of independent origin. The new reference genome was subsequently used to build a pan-genome based on five Australian Pgt isolates. Transcriptomes from germinated urediniospores and haustoria were separately assembled for pt. 21-0 and comparison of gene expression profiles showed differential expression in ∼10% of the genes each in germinated spores and haustoria. A total of 1,924 secreted proteins were predicted from the 21-0 transcriptome, of which 520 were classified as haustorial secreted proteins (HSPs). Comparison of 21-0 with two presumed clonal field derivatives of this lineage (collected in 1982 and 1984) that had evolved virulence on four additional resistance genes (Sr5, Sr11, Sr27, SrSatu) identified mutations in 25 HSP effector candidates. Some of these mutations could explain their novel virulence phenotypes.

Strategies for Wheat Stripe Rust Pathogenicity Identified by Transcriptome Sequencing.

Stripe rust caused by the fungus Puccinia striiformis f.sp. tritici (Pst) is a major constraint to wheat production worldwide. The molecular events that underlie Pst pathogenicity are largely unknown. Like all rusts, Pst creates a specialized cellular structure within host cells called the haustorium to obtain nutrients from wheat, and to secrete pathogenicity factors called effector proteins. We purified Pst haustoria and used next-generation sequencing platforms to assemble the haustorial transcriptome as well as the transcriptome of germinated spores. 12,282 transcripts were assembled from 454-pyrosequencing data and used as reference for digital gene expression analysis to compare the germinated uredinospores and haustoria transcriptomes based on Illumina RNAseq data. More than 400 genes encoding secreted proteins which constitute candidate effectors were identified from the haustorial transcriptome, with two thirds of these up-regulated in this tissue compared to germinated spores. RT-PCR analysis confirmed the expression patterns of 94 effector candidates. The analysis also revealed that spores rely mainly on stored energy reserves for growth and development, while haustoria take up host nutrients for massive energy production for biosynthetic pathways and the ultimate production of spores. Together, these studies substantially increase our knowledge of potential Pst effectors and provide new insights into the pathogenic strategies of this important organism.

Survey and analysis of microsatellites from transcript sequences in Phytophthora species: frequency, distribution, and potential as markers for the genus.

BACKGROUND: Members of the genus Phytophthora are notorious pathogens with world-wide distribution. The most devastating species include P. infestans, P. ramorum and P. sojae. In order to develop molecular methods for routinely characterizing their populations and to gain a better insight into the organization and evolution of their genomes, we used an in silico approach to survey and compare simple sequence repeats (SSRs) in transcript sequences from these three species. We compared the occurrence, relative abundance, relative density and cross-species transferability of the SSRs in these oomycetes. RESULTS: The number of SSRs in oomycetes transcribed sequences is low and long SSRs are rare. The in silico transferability of SSRs among the Phytophthora species was analyzed for all sets generated, and primers were selected on the basis of similarity as possible candidates for transferability to other Phytophthora species. Sequences encoding putative pathogenicity factors from all three Phytophthora species were also surveyed for presence of SSRs. However, no correlation between gene function and SSR abundance was observed. The SSR survey results, and the primer pairs designed for all SSRs from the three species, were deposited in a public database. CONCLUSION: In all cases the most common SSRs were trinucleotide repeat units with low repeat numbers. A proportion (7.5%) of primers could be transferred with 90% similarity between at least two species of Phytophthora. This information represents a valuable source of molecular markers for use in population genetics, genetic mapping and strain fingerprinting studies of oomycetes, and illustrates how genomic databases can be exploited to generate data-mining filters for SSRs before experimental validation.