Comparative genome analyses reveal sequence features reflecting distinct modes of host-adaptation between dicot and monocot powdery mildew

The massive 340 megabase genome of Anisogramma anomala, a biotrophic ascomycete that causes eastern filbert blight of hazelnut

BACKGROUND

The ascomycete fungus Anisogramma anomala causes Eastern Filbert Blight (EFB) on hazelnut (Corylus spp.) trees. It is a minor disease on its native host, the American hazelnut (C. americana), but is highly destructive on the commercially important European hazelnut (C. avellana). In North America, EFB has historically limited commercial production of hazelnut to west of the Rocky Mountains. A. anomala is an obligately biotrophic fungus that has not been grown in continuous culture, rendering its study challenging. There is a 15-month latency before symptoms appear on infected hazelnut trees, and only a sexual reproductive stage has been observed. Here we report the sequencing, annotation, and characterization of its genome.

RESULTS

The genome of A. anomala was assembled into 108 scaffolds totaling 342,498,352 nt with a GC content of 34.46%. Scaffold N50 was 33.3 Mb and L50 was 5. Nineteen scaffolds with lengths over 1 Mb constituted 99% of the assembly. Telomere sequences were identified on both ends of two scaffolds and on one end of another 10 scaffolds. Flow cytometry estimated the genome size of A. anomala at 370 Mb. The genome exhibits two-speed evolution, with 93% of the assembly as AT-rich regions (32.9% GC) and the other 7% as GC-rich (57.1% GC). The AT-rich regions consist predominantly of repeats with low gene content, while 90% of predicted protein coding genes were identified in GC-rich regions. Copia-like retrotransposons accounted for more than half of the genome. Evidence of repeat-induced point mutation (RIP) was identified throughout the AT-rich regions, and two copies of the rid gene and one of dim-2, the key genes in the RIP mutation pathway, were identified in the genome. Consistent with its homothallic sexual reproduction cycle, both MAT1-1 and MAT1-2 idiomorphs were found. We identified a large suite of genes likely involved in pathogenicity, including 614 carbohydrate active enzymes, 762 secreted proteins and 165 effectors.

CONCLUSIONS

This study reveals the genomic structure, composition, and putative gene function of the important pathogen A. anomala. It provides insight into the molecular basis of the pathogen's life cycle and a solid foundation for studying EFB.

Leaf abaxial immunity to powdery mildew in Arabidopsis is conferred by multiple defense mechanisms

Powdery mildew fungi are obligate biotrophic pathogens that only invade plant epidermal cells. There are two epidermal surfaces in every plant leaf: the adaxial (upper) side and the abaxial (lower) side. While both leaf surfaces can be susceptible to adapted powdery mildew fungi in many plant species, there have been observations of leaf abaxial immunity in some plant species including Arabidopsis. The genetic basis of such leaf abaxial immunity remains unknown. In this study, we tested a series of Arabidopsis mutants defective in one or more known defense pathways with the adapted powdery mildew isolate Golovinomyces cichoracearum UCSC1. We found that leaf abaxial immunity was significantly compromised in mutants impaired for both the EDS1/PAD4- and PEN2/PEN3-dependent defenses. Consistently, expression of EDS1-yellow fluorescent protein and PEN2-green fluorescent protein fusions from their respective native promoters in the respective eds1-2 and pen2-1 mutant backgrounds was higher in the abaxial epidermal cells than in the adaxial epidermal cells. Altogether, our results indicate that leaf abaxial immunity against powdery mildew in Arabidopsis is at least partially due to enhanced EDS1/PAD4- and PEN2/PEN3-dependent defenses. Such transcriptionally pre-programmed defense mechanisms may underlie leaf abaxial immunity in other plant species such as hemp and may be exploited for engineering adaxial immunity against powdery mildew fungi in crop plants.

Spray-induced gene silencing to identify powdery mildew gene targets and processes for powdery mildew control

Spray-induced gene silencing (SIGS) is an emerging tool for crop pest protection. It utilizes exogenously applied double-stranded RNA to specifically reduce pest target gene expression using endogenous RNA interference machinery. In this study, SIGS methods were developed and optimized for powdery mildew fungi, which are widespread obligate biotrophic fungi that infect agricultural crops, using the known azole-fungicide target cytochrome P450 51 (CYP51) in the Golovinomyces orontii-Arabidopsis thaliana pathosystem. Additional screening resulted in the identification of conserved gene targets and processes important to powdery mildew proliferation: apoptosis-antagonizing transcription factor in essential cellular metabolism and stress response; lipid catabolism genes lipase a, lipase 1, and acetyl-CoA oxidase in energy production; and genes involved in manipulation of the plant host via abscisic acid metabolism (9-cis-epoxycarotenoid dioxygenase, xanthoxin dehydrogenase, and a putative abscisic acid G-protein coupled receptor) and secretion of the effector protein, effector candidate 2. Powdery mildew is the dominant disease impacting grapes and extensive powdery mildew resistance to applied fungicides has been reported. We therefore developed SIGS for the Erysiphe necator-Vitis vinifera system and tested six successful targets identified using the G. orontii-A. thaliana system. For all targets tested, a similar reduction in powdery mildew disease was observed between systems. This indicates screening of broadly conserved targets in the G. orontii-A. thaliana pathosystem identifies targets and processes for the successful control of other powdery mildew fungi. The efficacy of SIGS on powdery mildew fungi makes SIGS an exciting prospect for commercial powdery mildew control.

A chromosome-scale genome assembly of the grape powdery mildew pathogen Erysiphe necator reveals its genomic architecture and previously unknown features of its biology

Erysiphe necator is an obligate fungal pathogen that causes grape powdery mildew, globally the most important disease on grapevines. Previous attempts to obtain a quality genome assembly for this pathogen were hindered by its high repetitive DNA content. Here, chromatin conformation capture (Hi-C) with long-read PacBio sequencing was combined to obtain a chromosome-scale assembly and a high-quality annotation for E. necator isolate EnFRAME01. The resulting 81.1 Mb genome assembly is 98% complete and consists of 34 scaffolds, 11 of which represent complete chromosomes. All chromosomes contain large centromeric-like regions and lack synteny to the 11 chromosomes of the cereal PM pathogen Blumeria graminis. Further analysis of their composition showed that repeats and transposable elements (TEs) occupy 62.7% of their content. TEs were almost evenly interspersed outside centromeric and telomeric regions and massively overlapped with regions of annotated genes, suggesting that they could have a significant functional impact. Abundant gene duplicates were observed as well, particularly in genes encoding candidate secreted effector proteins. Moreover, younger in age gene duplicates exhibited more relaxed selection pressure and were more likely to be located physically close in the genome than older duplicates. A total of 122 genes with copy number variations among six isolates of E. necator were also identified and were enriched in genes that were duplicated in EnFRAME01, indicating they may reflect an adaptive variation. Taken together, our study illuminates higher-order genomic architectural features of E. necator and provides a valuable resource for studying genomic structural variations in this pathogen. IMPORTANCE Grape powdery mildew caused by the ascomycete fungus Erysiphe necator is economically the most important and recurrent disease in vineyards across the world. The obligate biotrophic nature of E. necator hinders the use of typical genetic methods to elucidate its pathogenicity and adaptation to adverse conditions, and thus comparative genomics has been a major method to study its genome biology. However, the current reference genome of E. necator isolate C-strain is highly fragmented with many non-coding regions left unassembled. This incompleteness prohibits in-depth comparative genomic analyses and the study of genomic structural variations (SVs) that are known to affect several aspects of microbial life, including fitness, virulence, and host adaptation. By obtaining a chromosome-scale genome assembly and a high-quality gene annotation for E. necator, we reveal the organization of its chromosomal content, unearth previously unknown features of its biology, and provide a reference for studying genomic SVs in this pathogen.

Annual (2023) taxonomic update of RNA-directed RNA polymerase-encoding negative-sense RNA viruses (realm Riboviria: kingdom Orthornavirae: phylum Negarnaviricota)

In April 2023, following the annual International Committee on Taxonomy of Viruses (ICTV) ratification vote on newly proposed taxa, the phylum Negarnaviricota was amended and emended. The phylum was expanded by one new family, 14 new genera, and 140 new species. Two genera and 538 species were renamed. One species was moved, and four were abolished. This article presents the updated taxonomy of Negarnaviricota as now accepted by the ICTV.

Key Advances in the New Era of Genomics-Assisted Disease Resistance Improvement of Brassica Species

Disease resistance improvement remains a major focus in breeding programs as diseases continue to devastate Brassica production systems due to intensive cultivation and climate change. Genomics has paved the way to understand the complex genomes of Brassicas, which has been pivotal in the dissection of the genetic underpinnings of agronomic traits driving the development of superior cultivars. The new era of genomics-assisted disease resistance breeding has been marked by the development of high-quality genome references, accelerating the identification of disease resistance genes controlling both qualitative (major) gene and quantitative resistance. This facilitates the development of molecular markers for marker assisted selection and enables genome editing approaches for targeted gene manipulation to enhance the genetic value of disease resistance traits. This review summarizes the key advances in the development of genomic resources for Brassica species, focusing on improved genome references, based on long-read sequencing technologies and pangenome assemblies. This is further supported by the advances in pathogen genomics, which have resulted in the discovery of pathogenicity factors, complementing the mining of disease resistance genes in the host. Recognizing the co-evolutionary arms race between the host and pathogen, it is critical to identify novel resistance genes using crop wild relatives and synthetic cultivars or through genetic manipulation via genome-editing to sustain the development of superior cultivars. Integrating these key advances with new breeding techniques and improved phenotyping using advanced data analysis platforms will make disease resistance improvement in Brassica species more efficient and responsive to current and future demands.

Beyond Nuclear Ribosomal DNA Sequences: Evolution, Taxonomy, and Closest Known Saprobic Relatives of Powdery Mildew Fungi (Erysiphaceae) Inferred From Their First Comprehensive Genome-Scale Phylogenetic Analyses

Powdery mildew fungi (Erysiphaceae), common obligate biotrophic pathogens of many plants, including important agricultural and horticultural crops, represent a monophyletic lineage within the Ascomycota. Within the Erysiphaceae, molecular phylogenetic relationships and DNA-based species and genera delimitations were up to now mostly based on nuclear ribosomal DNA (nrDNA) phylogenies. This is the first comprehensive genome-scale phylogenetic analysis of this group using 751 single-copy orthologous sequences extracted from 24 selected powdery mildew genomes and 14 additional genomes from Helotiales, the fungal order that includes the Erysiphaceae. Representative genomes of all powdery mildew species with publicly available whole-genome sequencing (WGS) data that were of sufficient quality were included in the analyses. The 24 powdery mildew genomes included in the analysis represented 17 species belonging to eight out of 19 genera recognized within the Erysiphaceae. The epiphytic genera, all but one represented by multiple genomes, belonged each to distinct, well-supported lineages. Three hemiendophytic genera, each represented by a single genome, together formed the hemiendophytic lineage. Out of the 14 other taxa from the Helotiales, Arachnopeziza araneosa, a saprobic species, was the only taxon that grouped together with the 24 genome-sequenced powdery mildew fungi in a monophyletic clade. The close phylogenetic relationship between the Erysiphaceae and Arachnopeziza was revealed earlier by a phylogenomic study of the Leotiomycetes. Further analyses of powdery mildew and Arachnopeziza genomes may discover signatures of the evolutionary processes that have led to obligate biotrophy from a saprobic way of life. A separate phylogeny was produced using the 18S, 5.8S, and 28S nrDNA sequences of the same set of powdery mildew specimens and compared to the genome-scale phylogeny. The nrDNA phylogeny was largely congruent to the phylogeny produced using 751 orthologs. This part of the study has revealed multiple contamination and other quality issues in some powdery mildew genomes. We recommend that the presence of 28S, internal transcribed spacer (ITS), and 18S nrDNA sequences in powdery mildew WGS datasets that are identical to those determined by Sanger sequencing should be used to assess the quality of assemblies, in addition to the commonly used Benchmarking Universal Single-Copy Orthologs (BUSCO) values.

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.

New Records of Powdery Mildews from Taiwan: Erysiphe ipomoeae comb. nov., E. aff. betae on Buckwheat, and E. neolycopersici comb. nov. on Cardiospermum halicacabum

Erysiphe is the largest genus of powdery mildews (PMs), a group of obligate plant pathogenic fungi. Exploration of biodiversity generally relies on regional surveys and our aim is to investigate the PMs in Taiwan. Collections of the fungi on five plant species around agricultural environments were subjected to morphological and molecular characterization, using both internal transcribed spacer (ITS) and β-tubulin gene (TUB2) regions for the phylogenetic analyses. Erysipheipomoeae comb. nov., a species able to infect Ipomoea obscura and I. aquatica demonstrated by pathogenicity tests, has been neotypified. The two buckwheat species, Fagopyrum esculentum and F. tataricum, are found to be hosts of E. aff. betae. These results suggest that hosts in some plant families can be infected by more than one Erysiphe pathogen, e.g., Convolvulaceae by E. ipomoeae and E. convolvuli and Polygonaceae by E. polygoni and E. aff. betae, respectively. In addition, phylogenetic analyses of PMs on Cardiospermum halicacabum and tomato belonging to the E. aquilegiae complex are allocated under E. neolycopersici comb. nov. This extends the potential host range of E. aquilegiae complex to the plant family Sapindaceae. We conclude that awareness of the host associations of PMs can potentially benefit crop disease management.

Endoparasitic plants and fungi show evolutionary convergence across phylogenetic divisions

Endoparasitic plants are the most reduced flowering plants, spending most of their lives as a network of filaments within the tissues of their hosts. Despite their extraordinary life form, we know little about their biology. Research into a few species has revealed unexpected insights, such as the total loss of plastome, the reduction of the vegetative phase to a proembryonic stage, and elevated information exchange between host and parasite. To consolidate our understanding, we review life history, anatomy, and molecular genetics across the four independent lineages of endoparasitic plants. We highlight convergence across these clades and a striking trans-kingdom convergence in life history among endoparasitic plants and disparate lineages of fungi at the molecular and physiological levels. We hypothesize that parasitism of woody plants preselected for the endoparasitic life history, providing parasites a stable host environment and the necessary hydraulics to enable floral gigantism and/or high reproductive output. Finally, we propose a broader view of endoparasitic plants that connects research across disciplines, for example, pollen-pistil and graft incompatibility interactions and plant associations with various fungi. We shine a light on endoparasitic plants and their hosts as under-explored ecological microcosms ripe for identifying unexpected biological processes, interactions and evolutionary convergence.

Virulence Traits and Population Genomics of the Black Yeast Aureobasidium melanogenum

The black yeast-like fungus Aureobasidium melanogenum is an opportunistic human pathogen frequently found indoors. Its traits, potentially linked to pathogenesis, have never been systematically studied. Here, we examine 49 A. melanogenum strains for growth at 37 °C, siderophore production, hemolytic activity, and assimilation of hydrocarbons and human neurotransmitters and report within-species variability. All but one strain grew at 37 °C. All strains produced siderophores and showed some hemolytic activity. The largest differences between strains were observed in the assimilation of hydrocarbons and human neurotransmitters. We show for the first time that fungi from the order Dothideales can assimilate aromatic hydrocarbons. To explain the background, we sequenced the genomes of all 49 strains and identified genes putatively involved in siderophore production and hemolysis. Genomic analysis revealed a fairly structured population of A.melanogenum, raising the possibility that some phylogenetic lineages have higher virulence potential than others. Population genomics indicated that the species is strictly clonal, although more than half of the genomes were diploid. The existence of relatively heterozygous diploids in an otherwise clonal species is described for only the second time in fungi. The genomic and phenotypic data from this study should help to resolve the non-trivial taxonomy of the genus Aureobasidium and reduce the medical hazards of exploiting the biotechnological potential of other, non-pathogenic species of this genus.

Genome-Wide Identification of Powdery Mildew Resistance in Common Bean (Phaseolus vulgaris L.)

Genome-wide association studies (GWAS) have been utilized to detect genetic variations related to several agronomic traits and disease resistance in common bean. However, its application in the powdery mildew (PM) disease to identify candidate genes and their location in the common bean genome has not been fully addressed. Single-nucleotide polymorphism (SNP) genotyping with a BeadChip containing 5398 SNPs was used to detect genetic variations related to PM disease resistance in a panel of 211 genotypes grown under two field conditions for two consecutive years. Significant SNPs identified on chromosomes Pv04 and Pv10 were repeatable, ensuring the phenotypic data’s reliability and the causal relationship. A cluster of resistance genes was revealed on the Pv04 of the common bean genome, coiled-coil-nucleotide-binding site–leucine-rich repeat (CC-NBS-LRR, CNL), and Toll/interleukin-1 receptor-nucleotide-binding site–leucine-rich repeat type (TIR-NBS-LRR, TNL)-like resistance genes were identified. Furthermore, two resistance genes, Phavu_010G1320001g and Phavu_010G136800g, were also identified on Pv10. Further sequence analysis showed that these genes were homologs to the disease-resistance protein (RLM1A-like) and the putative disease-resistance protein (At4g11170.1) in Arabidopsis. Significant SNPs related to two LRR receptor-like kinases (RLK) were only identified on Pv11 in 2018. Many genes encoding the auxin-responsive protein, TIFY10A protein, growth-regulating factor five-like, ubiquitin-like protein, and cell wall RBR3-like protein related to PM disease resistance were identified nearby significant SNPs. These results suggested that the resistance to PM pathogen involves a network of many genes constitutively co-expressed.

AtSTP8, an endoplasmic reticulum-localised monosaccharide transporter from Arabidopsis, is recruited to the extrahaustorial membrane during powdery mildew infection

Biotrophic pathogens are believed to strategically manipulate sugar transport in host cells to enhance their access to carbohydrates. However, mechanisms of sugar translocation from host cells to biotrophic fungi such as powdery mildew across the plant-haustorium interface remain poorly understood. To investigate this question, systematic subcellular localisation analysis was performed for all the 14 members of the monosaccharide sugar transporter protein (STP) family in Arabidopsis thaliana. The best candidate AtSTP8 was further characterised for its transport properties in Saccharomyces cerevisiae and potential role in powdery mildew infection by gene ablation and overexpression in Arabidopsis. Our results showed that AtSTP8 was mainly localised to the endoplasmic reticulum (ER) and appeared to be recruited to the host-derived extrahaustorial membrane (EHM) induced by powdery mildew. Functional complementation assays in S. cerevisiae suggested that AtSTP8 can transport a broad spectrum of hexose substrates. Moreover, transgenic Arabidopsis plants overexpressing AtSTP8 showed increased hexose concentration in leaf tissues and enhanced susceptibility to powdery mildew. Our data suggested that the ER-localised sugar transporter AtSTP8 may be recruited to the EHM where it may be involved in sugar acquisition by haustoria of powdery mildew from host cells in Arabidopsis.

Application of Gamma Ray-Responsive Genes for Transcriptome-Based Phytodosimetry in Rice

Transcriptome-based dose–response curves were recently applied to the phytodosimetry of gamma radiation in a dicot plant, Arabidopsis thaliana, as an alternative biological assessment of genotoxicity using DNA damage response (DDR) genes. In the present study, we characterized gamma ray-responsive marker genes for transcriptome-based phytodosimetry in a monocot plant, rice (Oryza sativa L.), and compared different phytodosimetry models between rice and Arabidopsis using gamma-H2AX, comet, and quantitative transcriptomic assays. The transcriptome-based dose–response curves of four marker genes (OsGRG, OsMutS, OsRAD51, and OsRPA1) were reliably fitted to quadratic or exponential decay equations (r² > 0.99). However, the single or integrated dose–response curves of these genes were distinctive from the conventional models obtained by the gamma-H2AX or comet assays. In comparison, rice displayed a higher dose-dependency in the comet signal and OsRAD51 transcription, while the gamma-H2AX induction was more dose-dependent in Arabidopsis. The dose-dependent transcriptions of the selected gamma-ray-inducible marker genes, including OsGRG, OsMutS, OsRAD51, and OsRPA1 in rice and AtGRG, AtPARP1, AtRAD51, and AtRPA1E in Arabidopsis, were maintained similarly at different vegetative stages. These results suggested that the transcriptome-based phytodosimetry model should be further corrected with conventional genotoxicity- or DDR-based models despite the high reliability or dose-dependency of the model. In addition, the relative weighting of each gene in the integrated transcriptome-based dose–response model using multiple genes needs to be considered based on the trend and amplitude of the transcriptional change.

A haustorial-expressed lytic polysaccharide monooxygenase from the cucurbit powdery mildew pathogen Podosphaera xanthii contributes to the suppression of chitin-triggered immunity

Podosphaera xanthii is the main causal agent of cucurbit powdery mildew and a limiting factor of crop productivity. The lifestyle of this fungus is determined by the development of specialized parasitic structures inside epidermal cells, termed haustoria, that are responsible for the acquisition of nutrients and the release of effectors. A typical function of fungal effectors is the manipulation of host immunity, for example the suppression of pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI). Chitin is a major component of fungal cell walls, and chitin oligosaccharides are well-known PAMP elicitors. In this work, we examined the role of PHEC27213, the most highly expressed, haustorium-specific effector candidate of P. xanthii. According to different computational predictions, the protein folding of PHEC27213 was similar to that of lytic polysaccharide monooxygenases (LPMOs) and included a conserved histidine brace; however, PHEC27213 had low sequence similarity with LPMO proteins and displayed a putative chitin-binding domain that was different from the canonical carbohydrate-binding module. Binding and enzymatic assays demonstrated that PHEC27213 was able to bind and catalyse colloidal chitin, as well as chitooligosaccharides, acting as an LPMO. Furthermore, RNAi silencing experiments showed the potential of this protein to prevent the activation of chitin-triggered immunity. Moreover, proteins with similar features were found in other haustorium-forming fungal pathogens. Our results suggest that this protein is a new fungal LPMO that catalyses chitooligosaccharides, thus contributing to the suppression of plant immunity during haustorium development. To our knowledge, this is the first mechanism identified in the haustorium to suppress chitin signalling.

Categorization of Orthologous Gene Clusters in 92 Ascomycota Genomes Reveals Functions Important for Phytopathogenicity

Phytopathogenic Ascomycota are responsible for substantial economic losses each year, destroying valuable crops. The present study aims to provide new insights into phytopathogenicity in Ascomycota from a comparative genomic perspective. This has been achieved by categorizing orthologous gene groups (orthogroups) from 68 phytopathogenic and 24 non-phytopathogenic Ascomycota genomes into three classes: Core, (pathogen or non-pathogen) group-specific, and genome-specific accessory orthogroups. We found that (i) ~20% orthogroups are group-specific and accessory in the 92 Ascomycota genomes, (ii) phytopathogenicity is not phylogenetically determined, (iii) group-specific orthogroups have more enriched functional terms than accessory orthogroups and this trend is particularly evident in phytopathogenic fungi, (iv) secreted proteins with signal peptides and horizontal gene transfers (HGTs) are the two functional terms that show the highest occurrence and significance in group-specific orthogroups, (v) a number of other functional terms are also identified to have higher significance and occurrence in group-specific orthogroups. Overall, our comparative genomics analysis determined positive enrichment existing between orthogroup classes and revealed a prediction of what genomic characteristics make an Ascomycete phytopathogenic. We conclude that genes shared by multiple phytopathogenic genomes are more important for phytopathogenicity than those that are unique in each genome.

Genome Resource of Podosphaera xanthii, the Host-Specific Fungal Pathogen That Causes Cucurbit Powdery Mildew

Approximately 33 types of commonly consumed fruits and vegetables are members of the family Cucurbitaceae, making it an important crop family worldwide. However, pathogen resistance to pesticides and fungicides has become a growing problem in cultivation practices. The identification of the effector proteins in each unique fungus-host pair would help toward the development of strategies for preventing the infection of important crops. In this study, we characterized the genome of Podosphaera xanthii, the fungal pathogen that causes powdery mildew disease in cucurbitaceous plants. A first-draft genome of 209.08 MB was assembled and compared with those of 25 other fungal pathogens, particularly for identifying candidate secreted effector proteins. This draft genome can serve as a valuable resource for future genomic and proteomic studies of P. xanthii and its host-specific pathogenesis.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

A Hybrid Genome Assembly Resource for Podosphaera xanthii, the Main Causal Agent of Powdery Mildew Disease in Cucurbits

Podosphaera xanthii is the main causal agent of powdery mildew in cucurbits and, arguably, the most important fungal pathogen of cucurbit crops. Here, we present the first reference genome assembly for P. xanthii. We performed a hybrid genome assembly, using reads from Illumina NextSeq550 and PacBio Sequel S3. The short and long reads were assembled into 1,727 scaffolds with an N₅₀ size of 163,173 bp, resulting in a 142-Mb genome size. The combination of homology-based and ab initio predictions allowed the prediction of 14,911 complete genes. Repetitive sequences comprised 76.2% of the genome. Our P. xanthii genome assembly improves considerably the molecular resources for research on P. xanthii-cucurbit interactions and provides new opportunities for further genomics, transcriptomics, and evolutionary studies in powdery mildew fungi.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

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.

Effector mining from the Erysiphe pisi haustorial transcriptome identifies novel candidates involved in pea powdery mildew pathogenesis

Pea powdery mildew (PM) is an important fungal disease caused by an obligate biotroph, Erysiphe pisi (Ep), which significantly impacts pea production worldwide. The phytopathogen secretes a plethora of effectors, primarily through specialized infection structures termed haustoria, to establish a dynamic relationship with its host. To identify Ep effector candidates, a cDNA library of enriched haustoria from Ep-infected pea leaves was sequenced. The Ep transcriptome encodes 622 Ep candidate secreted proteins (CSPs), of which 167 were predicted to be candidate secreted effector proteins (CSEPs). Phylogenetic analysis indicates that Ep CSEPs are highly diverse, but, unlike cereal PM CSEPs, exhibit extensive sequence similarity with effectors from other PMs. Quantitative real-time PCR of a subset of EpCSEP/CSPs revealed that the majority are preferentially expressed in haustoria and exhibit infection stage-specific expression patterns. The functional roles of EpCSEP001, EpCSEP009 and EpCSP083 were probed by host-induced gene silencing (HIGS) via a double-stranded (ds) RNA-mediated RNAi approach. Foliar application of individual EpCSEP/CSP dsRNAs resulted in a marked reduction in PM disease symptoms. These findings were consistent with microscopic and molecular studies, suggesting that these Ep CSEP/CSPs play important roles in pea PM pathogenesis. Homology modelling revealed that EpCSEP001 and EpCSEP009 are analogous to fungal ribonucleases and belong to the RALPH family of effectors. This is the first study to identify and functionally validate candidate effectors from the agriculturally relevant pea PM, and highlights the utility of transcriptomics and HIGS to elucidate the key proteins associated with Ep pathogenesis.

The cytological basis of powdery mildew resistance in wild Chinese Vitis species

The wild Chinese grapevines (Vitis spp.) show varying levels of resistance to powdery mildew caused by Erysiphe necator that is an economically important disease of cultivated grapevines (Vitis vinifera). However, little information is available regarding the cytological mechanisms of powdery mildew resistance in these wild relatives. Here, we studied the cytological responses of three wild Chinese grapevine accessions after they were infected with E. necator (En) NAFU1 in comparison to the susceptible V. vinifera cv. 'Thompson Seedless' grape. The hyphal growth and sporulation of En NAFU1 were significantly restricted in wild species compared to 'Thompson Seedless', which appears to be associated with early cell wall deposition at the attempt sites, encasement of haustoria, and hypersensitive response-like cell death of penetrated epidermal cells. Moreover, endogenous free salicylic acid (SA) was more abundant in wild Chinese Vitis species than in 'Thompson Seedless' under pathogen-free condition. During En NAFU1 colonization, SA conjugates accumulated higher in wild grapevines than in 'Thompson Seedless'. In addition, the species-specific expression patterns of defense-associated genes during En NAFU1 colonization indicated that mechanisms underlying powdery mildew resistance are divergent among different wild Chinese Vitis species. These results contribute to understanding of mechanisms underlying defense responses of wild Chinese Vitis species against powdery mildew.

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

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.

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.

Whole-Genome Sequence of the Phlox Powdery Mildew Pathogen Golovinomyces magnicellulatus Strain FPH2017-1

Powdery mildew (PM) fungi are obligate biotrophs capable of infecting diverse plant hosts, ranging from monocotyledonous agricultural crops to dicotyledonous ornamental crops. The PM lifestyle poses significant challenges for studying these pathogens in isolation from their host. We present a draft genome of Golovinomyces magnicellulatus, a host-specific PM on Phlox species.

Powdery mildew (PM) fungi are obligate biotrophs capable of infecting diverse plant hosts, ranging from monocotyledonous agricultural crops to dicotyledonous ornamental crops. The PM lifestyle poses significant challenges for studying these pathogens in isolation from their host. We present a draft genome of Golovinomyces magnicellulatus, a host-specific PM on Phlox species.

The haustorial transcriptome of the cucurbit pathogen Podosphaera xanthii reveals new insights into the biotrophy and pathogenesis of powdery mildew fungi

Background

Podosphaera xanthii is the main causal agent of powdery mildew disease in cucurbits and is responsible for important yield losses in these crops worldwide. Powdery mildew fungi are obligate biotrophs. In these parasites, biotrophy is determined by the presence of haustoria, which are specialized structures of parasitism developed by these fungi for the acquisition of nutrients and the delivery of effectors. Detailed molecular studies of powdery mildew haustoria are scarce due mainly to difficulties in their isolation. Therefore, their analysis is considered an important challenge for powdery mildew research. The aim of this work was to gain insights into powdery mildew biology by analysing the haustorial transcriptome of P. xanthii.

Results

Prior to RNA isolation and massive-scale mRNA sequencing, a flow cytometric approach was developed to isolate P. xanthii haustoria free of visible contaminants. Next, several commercial kits were used to isolate total RNA and to construct the cDNA and Illumina libraries that were finally sequenced by the Illumina NextSeq system. Using this approach, the maximum amount of information from low-quality RNA that could be obtained was used to accomplish the de novo assembly of the P. xanthii haustorial transcriptome. The subsequent analysis of this transcriptome and comparison with the epiphytic transcriptome allowed us to identify the importance of several biological processes for haustorial cells such as protection against reactive oxygen species, the acquisition of different nutrients and genetic regulation mediated by non-coding RNAs. In addition, we could also identify several secreted proteins expressed exclusively in haustoria such as cell adhesion proteins that have not been related to powdery mildew biology to date.

Conclusions

This work provides a novel approach to study the molecular aspects of powdery mildew haustoria. In addition, the results of this study have also allowed us to identify certain previously unknown processes and proteins involved in the biology of powdery mildews that could be essential for their biotrophy and pathogenesis.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5938-0) contains supplementary material, which is available to authorized users.

Two novel fungal negative-strand RNA viruses related to mymonaviruses and phenuiviruses in the shiitake mushroom (Lentinula edodes)

There is still limited information on the diversity of (-)ssRNA viruses that infect fungi. Here, we have discovered two novel (-)ssRNA mycoviruses in the shiitake mushroom (Lentinula edodes). The first virus has a monopartite RNA genome and relates to that of mymonaviruses (Mononegavirales), especially to Hubei rhabdo-like virus 4 from arthropods and thus designated as Lentinula edodes negative-strand RNA virus 1. The second virus has a putative bipartite RNA genome and is related to the recently discovered bipartite or tripartite phenui-like viruses (Bunyavirales) associated with plants and ticks, and designated as Lentinula edodes negative-strand RNA virus 2 (LeNSRV2). LeNSRV2 is likely the first segmented (-)ssRNA virus known to infect fungi. Its smaller RNA segment encodes a putative nucleocapsid and a plant MP-like protein using a potential ambisense coding strategy. These findings enhance our understanding of the diversity, evolution and spread of (-)ssRNA viruses in fungi.

Powdery Mildews Are Characterized by Contracted Carbohydrate Metabolism and Diverse Effectors to Adapt to Obligate Biotrophic Lifestyle

Powdery mildew is a widespread plant disease caused by obligate biotrophic fungal pathogens involving species-specific interactions between host and parasite. To gain genomic insights into the underlying obligate biotrophic mechanisms, we analyzed 15 microbial genomes covering powdery and downy mildews and rusts. We observed a genome-wide, massive contraction of multiple gene families in powdery mildews, such as enzymes in the carbohydrate metabolism pathway, when compared with ascomycete phytopathogens, while the fatty acid metabolism pathway maintained its integrity. We also observed significant differences in candidate secreted effector protein (CSEP) families between monocot and dicot powdery mildews, perhaps due to different selection forces. While CSEPs in monocot mildews are likely subject to positive selection causing rapid expansion, CSEP families in dicot mildews are shrinking under strong purifying selection. Our results not only illustrate obligate biotrophic mechanisms of powdery mildews driven by gene family evolution in nutrient metabolism, but also demonstrate how the divergence of CSEPs between monocot and dicot lineages might contribute to species-specific adaption.