First Draft Genome Assemblies of Pleochaeta shiraiana and Phyllactinia moricola, Two Tree-Parasitic Powdery Mildew Fungi with Hemiendophytic Mycelia.

Powdery mildew fungi (Erysiphaceae) are widespread obligate biotrophic plant pathogens. Thus, applying genetic and omics approaches to study these fungi remains a major challenge, particularly for species with hemiendophytic mycelium. These belong to a distinct phylogenetic lineage within the family Erysiphaceae. To date, only a single draft genome assembly is available for this clade, obtained for Leveillula taurica. Here, we generated the first draft genome assemblies of Pleochaeta shiraiana and Phyllactinia moricola, two tree-parasitic powdery mildew species with hemiendophytic mycelium, representing two genera that have not yet been investigated with genomics tools. The Pleochaeta shiraiana assembly was 96,769,103 bp in length and consisted of 14,447 scaffolds, and the Phyllactinia moricola assembly was 180,382,532 bp in length on 45,569 scaffolds. Together with the draft genome of L. taurica, these resources will be pivotal for understanding the molecular basis of the lifestyle of these fungi, which is unique within the family Erysiphaceae.

The leucine-rich repeats in allelic barley MLA immune receptors define specificity towards sequence-unrelated powdery mildew avirulence effectors with a predicted common RNase-like fold.

Nucleotide-binding domain leucine-rich repeat-containing receptors (NLRs) in plants can detect avirulence (AVR) effectors of pathogenic microbes. The Mildew locus a (Mla) NLR gene has been shown to confer resistance against diverse fungal pathogens in cereal crops. In barley, Mla has undergone allelic diversification in the host population and confers isolate-specific immunity against the powdery mildew-causing fungal pathogen Blumeria graminis forma specialis hordei (Bgh). We previously isolated the Bgh effectors AVRA1, AVRA7, AVRA9, AVRA13, and allelic AVRA10/AVRA22, which are recognized by matching MLA1, MLA7, MLA9, MLA13, MLA10 and MLA22, respectively. Here, we extend our knowledge of the Bgh effector repertoire by isolating the AVRA6 effector, which belongs to the family of catalytically inactive RNase-Like Proteins expressed in Haustoria (RALPHs). Using structural prediction, we also identified RNase-like folds in AVRA1, AVRA7, AVRA10/AVRA22, and AVRA13, suggesting that allelic MLA recognition specificities could detect structurally related avirulence effectors. To better understand the mechanism underlying the recognition of effectors by MLAs, we deployed chimeric MLA1 and MLA6, as well as chimeric MLA10 and MLA22 receptors in plant co-expression assays, which showed that the recognition specificity for AVRA1 and AVRA6 as well as allelic AVRA10 and AVRA22 is largely determined by the receptors' C-terminal leucine-rich repeats (LRRs). The design of avirulence effector hybrids allowed us to identify four specific AVRA10 and five specific AVRA22 aa residues that are necessary to confer MLA10- and MLA22-specific recognition, respectively. This suggests that the MLA LRR mediates isolate-specific recognition of structurally related AVRA effectors. Thus, functional diversification of multi-allelic MLA receptors may be driven by a common structural effector scaffold, which could be facilitated by proliferation of the RALPH effector family in the pathogen genome.

Ultraviolet Mutagenesis Coupled with Next-Generation Sequencing as a Method for Functional Interrogation of Powdery Mildew Genomes.

Powdery mildews are obligate biotrophic fungal pathogens causing important diseases of plants worldwide. Very little is known about the requirements for their pathogenicity at the molecular level. This is largely due to the inability to culture these organisms in vitro or to modify them genetically. Here, we describe a mutagenesis procedure based on ultraviolet (UV) irradiation to accumulate mutations in the haploid genome of the barley powdery mildew pathogen Blumeria graminis f. sp. hordei. Exposure of B. graminis f. sp. hordei conidia to different durations of UV-C radiation (10 s to 12 min) resulted in a reduced number of macroscopically visible fungal colonies. B. graminis f. sp. hordei colony number was negatively correlated with exposure time and the total number of consecutive cycles of UV irradiation. Dark incubation following UV exposure further reduced fungal viability, implying that photoreactivation is an important component of DNA repair in B. graminis f. sp. hordei. After several rounds of UV mutagenesis, we selected two mutant isolates in addition to the parental B. graminis f. sp. hordei K1 isolate for whole-genome resequencing. By combining automated prediction of sequence variants and their manual validation, we identified unique UV-induced mutations in the genomes of the two isolates. Most of these mutations were in the up- or downstream regions of genes or in the intergenic space. Some of the variants detected in genes led to predicted missense mutations. As an additional insight, our bioinformatic analyses revealed a complex population structure within supposedly clonal B. graminis f. sp. hordei isolates.

Rapid evolution in plant-microbe interactions - a molecular genomics perspective.

Rapid (co-)evolution at multiple timescales is a hallmark of plant-microbe interactions. The mechanistic basis for the rapid evolution largely rests on the features of the genomes of the interacting partners involved. Here, we review recent insights into genomic characteristics and mechanisms that enable rapid evolution of both plants and phytopathogens. These comprise fresh insights in allelic series of matching pairs of resistance and avirulence genes, the generation of novel pathogen effectors, the recently recognised small RNA warfare, and genomic aspects of secondary metabolite biosynthesis. In addition, we discuss the putative contributions of permissive host environments, transcriptional plasticity and the role of ploidy on the interactions. We conclude that the means underlying the rapid evolution of plant-microbe interactions are multifaceted and depend on the particular nature of each interaction.

The Parauncinula polyspora Draft Genome Provides Insights into Patterns of Gene Erosion and Genome Expansion in Powdery Mildew Fungi.

Due to their comparatively small genome size and short generation time, fungi are exquisite model systems to study eukaryotic genome evolution. Powdery mildew fungi present an exceptional case because of their strict host dependency (termed obligate biotrophy) and the atypical size of their genomes (>100 Mb). This size expansion is largely due to the pervasiveness of transposable elements on 70% of the genome and is associated with the loss of multiple conserved ascomycete genes required for a free-living lifestyle. To date, little is known about the mechanisms that drove these changes, and information on ancestral powdery mildew genomes is lacking. We report genome analysis of the early-diverged and exclusively sexually reproducing powdery mildew fungus Parauncinula polyspora, which we performed on the basis of a natural leaf epiphytic metapopulation sample. In contrast to other sequenced species of this taxonomic group, the assembled P. polyspora draft genome is surprisingly small (<30 Mb), has a higher content of conserved ascomycete genes, and is sparsely equipped with transposons (<10%), despite the conserved absence of a common defense mechanism involved in constraining repetitive elements. We speculate that transposable element spread might have been limited by this pathogen's unique reproduction strategy and host features and further hypothesize that the loss of conserved ascomycete genes may promote the evolutionary isolation and host niche specialization of powdery mildew fungi. Limitations associated with this evolutionary trajectory might have been in part counteracted by the evolution of plastic, transposon-rich genomes and/or the expansion of gene families encoding secreted virulence proteins.IMPORTANCE Powdery mildew fungi are widespread and agronomically relevant phytopathogens causing major yield losses. Their genomes have disproportionately large numbers of mobile genetic elements, and they have experienced a significant loss of highly conserved fungal genes. In order to learn more about the evolutionary history of this fungal group, we explored the genome of an Asian oak tree pathogen, Parauncinula polyspora, a species that diverged early during evolution from the remaining powdery mildew fungi. We found that the P. polyspora draft genome is comparatively compact, has a low number of protein-coding genes, and, despite the absence of a dedicated genome defense system, lacks the massive proliferation of repetitive sequences. Based on these findings, we infer an evolutionary trajectory that shaped the genomes of powdery mildew fungi.

Smut infection of perennial hosts: the genome and the transcriptome of the Brassicaceae smut fungus Thecaphora thlaspeos reveal functionally conserved and novel effectors.

Biotrophic fungal plant pathogens can balance their virulence and form intricate relationships with their hosts. Sometimes, this leads to systemic host colonization over long time scales without macroscopic symptoms. However, how plant-pathogenic endophytes manage to establish their sustained systemic infection remains largely unknown. Here, we present a genomic and transcriptomic analysis of Thecaphora thlaspeos. This relative of the well studied grass smut Ustilago maydis is the only smut fungus adapted to Brassicaceae hosts. Its ability to overwinter with perennial hosts and its systemic plant infection including roots are unique characteristics among smut fungi. The T. thlaspeos genome was assembled to the chromosome level. It is a typical smut genome in terms of size and genome characteristics. In silico prediction of candidate effector genes revealed common smut effector proteins and unique members. For three candidates, we have functionally demonstrated effector activity. One of these, TtTue1, suggests a potential link to cold acclimation. On the plant side, we found evidence for a typical immune response as it is present in other infection systems, despite the absence of any macroscopic symptoms during infection. Our findings suggest that T. thlaspeos distinctly balances its virulence during biotrophic growth ultimately allowing for long-lived infection of its perennial hosts.

Small RNAs from cereal powdery mildew pathogens may target host plant genes.

Small RNAs (sRNAs) play a key role in eukaryotic gene regulation, for example by gene silencing via RNA interference (RNAi). The biogenesis of sRNAs depends on proteins that are generally conserved in all eukaryotic lineages, yet some species that lack part or all the components of the mechanism exist. Here we explored the presence of the RNAi machinery and its expression as well as the occurrence of sRNA candidates and their putative endogenous as well as host targets in phytopathogenic powdery mildew fungi. We focused on the species Blumeria graminis, which occurs in various specialized forms (formae speciales) that each have a strictly limited host range. B. graminis f. sp. hordei and B. graminis f. sp. tritici, colonizing barley and wheat, respectively, have genomes that are characterized by extensive gene loss. Nonetheless, we find that the RNAi machinery appears to be largely complete and expressed during infection. sRNA sequencing data enabled the identification of putative sRNAs in both pathogens. While a considerable part of the sRNA candidates have predicted target sites in endogenous genes and transposable elements, a small proportion appears to have targets in planta, suggesting potential cross-kingdom RNA transfer between powdery mildew fungi and their respective plant hosts.

Signatures of host specialization and a recent transposable element burst in the dynamic one-speed genome of the fungal barley powdery mildew pathogen.

BACKGROUND: Powdery mildews are biotrophic pathogenic fungi infecting a number of economically important plants. The grass powdery mildew, Blumeria graminis, has become a model organism to study host specialization of obligate biotrophic fungal pathogens. We resolved the large-scale genomic architecture of B. graminis forma specialis hordei (Bgh) to explore the potential influence of its genome organization on the co-evolutionary process with its host plant, barley (Hordeum vulgare). RESULTS: The near-chromosome level assemblies of the Bgh reference isolate DH14 and one of the most diversified isolates, RACE1, enabled a comparative analysis of these haploid genomes, which are highly enriched with transposable elements (TEs). We found largely retained genome synteny and gene repertoires, yet detected copy number variation (CNV) of secretion signal peptide-containing protein-coding genes (SPs) and locally disrupted synteny blocks. Genes coding for sequence-related SPs are often locally clustered, but neither the SPs nor the TEs reside preferentially in genomic regions with unique features. Extended comparative analysis with different host-specific B. graminis formae speciales revealed the existence of a core suite of SPs, but also isolate-specific SP sets as well as congruence of SP CNV and phylogenetic relationship. We further detected evidence for a recent, lineage-specific expansion of TEs in the Bgh genome. CONCLUSIONS: The characteristics of the Bgh genome (largely retained synteny, CNV of SP genes, recently proliferated TEs and a lack of significant compartmentalization) are consistent with a "one-speed" genome that differs in its architecture and (co-)evolutionary pattern from the "two-speed" genomes reported for several other filamentous phytopathogens.

The Plant-Dependent Life Cycle of Thecaphora thlaspeos: A Smut Fungus Adapted to Brassicaceae.

Smut fungi are globally distributed plant pathogens that infect agriculturally important crop plants such as maize or potato. To date, molecular studies on plant responses to smut fungi are challenging due to the genetic complexity of their host plants. Therefore, we set out to investigate the known smut fungus of Brassicaceae hosts, Thecaphora thlaspeos. T. thlaspeos infects different Brassicaceae plant species throughout Europe, including the perennial model plant Arabis alpina. In contrast to characterized smut fungi, mature and dry T. thlaspeos teliospores germinated only in the presence of a plant signal. An infectious filament emerges from the teliospore, which can proliferate as haploid filamentous cultures. Haploid filaments from opposite mating types mate, similar to sporidia of the model smut fungus Ustilago maydis. Consistently, the a and b mating locus genes are conserved. Infectious filaments can penetrate roots and aerial tissues of host plants, causing systemic colonization along the vasculature. Notably, we could show that T. thlaspeos also infects Arabidopsis thaliana. Exploiting the genetic resources of A. thaliana and Arabis alpina will allow us to characterize plant responses to smut infection in a comparative manner and, thereby, characterize factors for endophytic growth as well as smut fungi virulence in dicot plants.

Genetic Manipulation of the Plant Pathogen Ustilago maydis to Study Fungal Biology and Plant Microbe Interactions.

Gene deletion plays an important role in the analysis of gene function. One of the most efficient methods to disrupt genes in a targeted manner is the replacement of the entire gene with a selectable marker via homologous recombination. During homologous recombination, exchange of DNA takes place between sequences with high similarity. Therefore, linear genomic sequences flanking a target gene can be used to specifically direct a selectable marker to the desired integration site. Blunt ends of the deletion construct activate the cell's DNA repair systems and thereby promote integration of the construct either via homologous recombination or by non-homologous-end-joining. In organisms with efficient homologous recombination, the rate of successful gene deletion can reach more than 50% making this strategy a valuable gene disruption system. The smut fungus Ustilago maydis is a eukaryotic model microorganism showing such efficient homologous recombination. Out of its about 6,900 genes, many have been functionally characterized with the help of deletion mutants, and repeated failure of gene replacement attempts points at essential function of the gene. Subsequent characterization of the gene function by tagging with fluorescent markers or mutations of predicted domains also relies on DNA exchange via homologous recombination. Here, we present the U. maydis strain generation strategy in detail using the simplest example, the gene deletion.

PCR amplification of repetitive DNA: a limitation to genome editing technologies and many other applications.

Designer transcription-activator like effectors (TALEs) is a promising technology and made it possible to edit genomes with higher specificity. Such specific engineering and gene regulation technologies are also being developed using RNA-binding proteins like PUFs and PPRs. The main feature of TALEs, PUFs and PPRs is their repetitive DNA/RNA-binding domains which have single nucleotide binding specificity. Available kits today allow researchers to assemble these repetitive domains in any combination they desire when generating TALEs for gene targeting and editing. However, PCR amplifications of such repetitive DNAs are highly problematic as these mostly fail, generating undesired artifact products or deletions. Here we describe the molecular mechanisms leading to these artifacts. We tested our models also in plasmid templates containing one copy versus two copies of GFP-coding sequence arranged as either direct or inverted repeats. Some limited solutions in amplifying repetitive DNA regions are described.

Reevaluation of the reliability and usefulness of the somatic homologous recombination reporter lines.

A widely used approach for assessing genome instability in plants makes use of somatic homologous recombination (SHR) reporter lines. Here, we review the published characteristics and uses of SHR lines. We found a lack of detailed information on these lines and a lack of sufficient evidence that they report only homologous recombination. We postulate that instead of SHR, these lines might be reporting a number of alternative stress-induced stochastic events known to occur at transcriptional, posttranscriptional, and posttranslational levels. We conclude that the reliability and usefulness of the somatic homologous recombination reporter lines requires revision. Thus, more detailed information about these reporter lines is needed before they can be used with confidence to measure genome instability, including the complete sequences of SHR constructs, the genomic location of reporter genes and, importantly, molecular evidence that reconstituted gene expression in these lines is indeed a result of somatic recombination.