Single-Molecule Real-Time Sequencing Combined with Optical Mapping Yields Completely Finished Fungal Genome

Verticillium dahliae LysM effectors differentially contribute to virulence on plant hosts

Chitin-binding lysin motif (LysM) effectors contribute to the virulence of various plant-pathogenic fungi that are causal agents of foliar diseases. Here, we report the LysM effectors of the soil-borne fungal vascular wilt pathogen Verticillium dahliae. Comparative genomics revealed three core LysM effectors that are conserved in a collection of V. dahliae strains. Remarkably, and in contrast with the previously studied LysM effectors of other plant pathogens, no expression of core LysM effectors was monitored in planta in a taxonomically diverse panel of host plants. Moreover, targeted deletion of the individual LysM effector genes in V. dahliae strain JR2 did not compromise virulence in infections on Arabidopsis, tomato or Nicotiana benthamiana. Interestingly, an additional lineage-specific LysM effector is encoded in the genome of V. dahliae strain VdLs17, but not in any other V. dahliae strain sequenced to date. Remarkably, this lineage-specific effector is expressed in planta and contributes to the virulence of V. dahliae strain VdLs17 on tomato, but not on Arabidopsis or N. benthamiana. Functional analysis revealed that this LysM effector binds chitin, is able to suppress chitin-induced immune responses and protects fungal hyphae against hydrolysis by plant hydrolytic enzymes. Thus, in contrast with the core LysM effectors of V. dahliae, this lineage-specific LysM effector of strain VdLs17 contributes to virulence in planta.

Complete genome sequence of Bacillus velezensis LM2303, a biocontrol strain isolated from the dung of wild yak inhabited Qinghai-Tibet plateau

Bacillus velezensis LM2303 is a biocontrol strain with a broad inhibitory spectrum against plant pathogens, isolated from the dung of wild yak inhabited Qinghai-Tibet plateau, China. Here we present its complete genome sequence, which consists of a single, circular chromosome of 3,989,393bp with a 46.68% G+C content. Genome analysis revealed genes encoding specialized functions for the biosynthesis of antifungal metabolites and antibacterial metabolites, the promotion of plant growth, the alleviation of oxidative stress and nutrient utilization. And the biosynthesis of antimicrobial metabolites in strain LM2303 was confirmed by biochemical analysis, while its plant growth promoting traits were confirmed by inoculation tests. Our results will establish a better foundation for further studies and biocontrol application of B. velezensis LM2303.

Genome sequence of Plasmopara viticola and insight into the pathogenic mechanism

Plasmopara viticola causes downy mildew disease of grapevine which is one of the most devastating diseases of viticulture worldwide. Here we report a 101.3 Mb whole genome sequence of P. viticola isolate 'JL-7-2' obtained by a combination of Illumina and PacBio sequencing technologies. The P. viticola genome contains 17,014 putative protein-coding genes and has ~26% repetitive sequences. A total of 1,301 putative secreted proteins, including 100 putative RXLR effectors and 90 CRN effectors were identified in this genome. In the secretome, 261 potential pathogenicity genes and 95 carbohydrate-active enzymes were predicted. Transcriptional analysis revealed that most of the RXLR effectors, pathogenicity genes and carbohydrate-active enzymes were significantly up-regulated during infection. Comparative genomic analysis revealed that P. viticola evolved independently from the Arabidopsis downy mildew pathogen Hyaloperonospora arabidopsidis. The availability of the P. viticola genome provides a valuable resource not only for comparative genomic analysis and evolutionary studies among oomycetes, but also enhance our knowledge on the mechanism of interactions between this biotrophic pathogen and its host.

Genome sequence of Plasmopara viticola and insight into the pathogenic mechanism

Plasmopara viticola causes downy mildew disease of grapevine which is one of the most devastating diseases of viticulture worldwide. Here we report a 101.3 Mb whole genome sequence of P. viticola isolate 'JL-7-2' obtained by a combination of Illumina and PacBio sequencing technologies. The P. viticola genome contains 17,014 putative protein-coding genes and has ~26% repetitive sequences. A total of 1,301 putative secreted proteins, including 100 putative RXLR effectors and 90 CRN effectors were identified in this genome. In the secretome, 261 potential pathogenicity genes and 95 carbohydrate-active enzymes were predicted. Transcriptional analysis revealed that most of the RXLR effectors, pathogenicity genes and carbohydrate-active enzymes were significantly up-regulated during infection. Comparative genomic analysis revealed that P. viticola evolved independently from the Arabidopsis downy mildew pathogen Hyaloperonospora arabidopsidis. The availability of the P. viticola genome provides a valuable resource not only for comparative genomic analysis and evolutionary studies among oomycetes, but also enhance our knowledge on the mechanism of interactions between this biotrophic pathogen and its host.

Rhamnose synthase activity is required for pathogenicity of the vascular wilt fungus Verticillium dahliae

The initial interaction of a pathogenic fungus with its host is complex and involves numerous metabolic pathways and regulatory proteins. Considerable attention has been devoted to proteins that play a crucial role in these interactions, with an emphasis on so-called effector molecules that are secreted by the invading microbe to establish the symbiosis. However, the contribution of other types of molecules, such as glycans, is less well appreciated. Here, we present a random genetic screen that enabled us to identify 58 novel candidate genes that are involved in the pathogenic potential of the fungal pathogen Verticillium dahliae, which causes vascular wilt diseases in over 200 dicotyledonous plant species, including economically important crops. One of the candidate genes that was identified concerns a putative biosynthetic gene involved in nucleotide sugar precursor formation, as it encodes a putative nucleotide-rhamnose synthase/epimerase-reductase (NRS/ER). This enzyme has homology to bacterial enzymes involved in the biosynthesis of the nucleotide sugar deoxy-thymidine diphosphate (dTDP)-rhamnose, a precursor of L-rhamnose, which has been shown to be required for virulence in several human pathogenic bacteria. Rhamnose is known to be a minor cell wall glycan in fungi and has therefore not been suspected as a crucial molecule in fungal-host interactions. Nevertheless, our study shows that deletion of the VdNRS/ER gene from the V. dahliae genome results in complete loss of pathogenicity on tomato and Nicotiana benthamiana plants, whereas vegetative growth and sporulation are not affected. We demonstrate that VdNRS/ER is a functional enzyme in the biosynthesis of uridine diphosphate (UDP)-rhamnose, and further analysis has revealed that VdNRS/ER deletion strains are impaired in the colonization of tomato roots. Collectively, our results demonstrate that rhamnose, although only a minor cell wall component, is essential for the pathogenicity of V. dahliae.

Purification of High Molecular Weight Genomic DNA from Powdery Mildew for Long-Read Sequencing

The powdery mildew fungi are a group of economically important fungal plant pathogens. Relatively little is known about the molecular biology and genetics of these pathogens, in part due to a lack of well-developed genetic and genomic resources. These organisms have large, repetitive genomes, which have made genome sequencing and assembly prohibitively difficult. Here, we describe methods for the collection, extraction, purification and quality control assessment of high molecular weight genomic DNA from one powdery mildew species, Golovinomyces cichoracearum. The protocol described includes mechanical disruption of spores followed by an optimized phenol/chloroform genomic DNA extraction. A typical yield was 7 µg DNA per 150 mg conidia. The genomic DNA that is isolated using this procedure is suitable for long-read sequencing (i.e., > 48.5 kbp). Quality control measures to ensure the size, yield, and purity of the genomic DNA are also described in this method. Sequencing of the genomic DNA of the quality described here will allow for the assembly and comparison of multiple powdery mildew genomes, which in turn will lead to a better understanding and improved control of this agricultural pathogen.

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.

Broad taxonomic characterization of Verticillium wilt resistance genes reveals an ancient origin of the tomato Ve1 immune receptor

Plant-pathogenic microbes secrete effector molecules to establish themselves on their hosts, whereas plants use immune receptors to try and intercept such effectors in order to prevent pathogen colonization. The tomato cell surface-localized receptor Ve1 confers race-specific resistance against race 1 strains of the soil-borne vascular wilt fungus Verticillium dahliae which secrete the Ave1 effector. Here, we describe the cloning and characterization of Ve1 homologues from tobacco (Nicotiana glutinosa), potato (Solanum tuberosum), wild eggplant (Solanum torvum) and hop (Humulus lupulus), and demonstrate that particular Ve1 homologues govern resistance against V. dahliae race 1 strains through the recognition of the Ave1 effector. Phylogenetic analysis shows that Ve1 homologues are widely distributed in land plants. Thus, our study suggests an ancient origin of the Ve1 immune receptor in the plant kingdom.

Comparative analyses of secreted proteins in plant pathogenic smut fungi and related basidiomycetes

In the ten years since the genome sequence of the basidiomycete corn smut fungus Ustilago maydis was published, additional genomes of smut species infecting different hosts became available. In addition, the genomes of related Malassezia species causing skin diseases and of Pseudozyma species not known to infect plants were determined. As secreted proteins are critical virulence determinants in U. maydis we compare here the secretomes of 12 basidiomycete species to gain information about their composition and conservation. For this we classify secreted proteins into those with and without domains using InterPro scans. Homology among proteins is inferred by building clusters based on pairwise similarities and cluster presence is then assessed in the different species. We detect in particular a strong correspondence between the secretomes of Pseudozyma species and plant infecting smuts. Furthermore, we identify a high proportion of secreted proteins to be part of gene families and present an advancement of the CRISPR-Cas9 technology for simultaneous disruption of multiple genes in U. maydis using five genes of the eff1 family as example.

Centrochromatin of Fungi

The centromere is an essential chromosomal locus that dictates the nucleation point for assembly of the kinetochore and subsequent attachment of spindle microtubules during chromosome segregation. Research over the last decades demonstrated that centromeres are defined by a combination of genetic and epigenetic factors. Recent work showed that centromeres are quite diverse and flexible and that many types of centromere sequences and centromeric chromatin ("centrochromatin") have evolved. The kingdom of the fungi serves as an outstanding example of centromere plasticity, including organisms with centromeres as diverse as 0.15-300 kb in length, and with different types of chromatin states for most species examined thus far. Some of the species in the less familiar taxa provide excellent opportunities to help us better understand centromere biology in all eukaryotes, which may improve treatment options against fungal infection, and biotechnologies based on fungi. This review summarizes the current knowledge of fungal centromeres and centrochromatin, including an outlook for future research.

A gapless genome sequence of the fungus Botrytis cinerea

Following earlier incomplete and fragmented versions of a genome sequence for the grey mould Botrytis cinerea, a gapless, near-finished genome sequence for B. cinerea strain B05.10 is reported. The assembly comprised 18 chromosomes and was confirmed by an optical map and a genetic map based on approximately 75 000 single nucleotide polymorphism (SNP) markers. All chromosomes contained fully assembled centromeric regions, and 10 chromosomes had telomeres on both ends. The genetic map consisted of 4153 cM and a comparison of the genetic distances with the physical distances identified 40 recombination hotspots. The linkage map also identified two mutations, located in the previously described genes Bos1 and BcsdhB, that conferred resistance to the fungicides boscalid and iprodione. The genome was predicted to encode 11 701 proteins. RNAseq data from >20 different samples were used to validate and improve gene models. Manual curation of chromosome 1 revealed interesting features, such as the occurrence of a dicistronic transcript and fully overlapping genes in opposite orientations, as well as many spliced antisense transcripts. Manual curation also revealed that the untranslated regions (UTRs) of genes can be complex and long, with many UTRs exceeding lengths of 1 kb and possessing multiple introns. Community annotation is in progress.

The Age of Effectors: Genome-Based Discovery and Applications

Microbial pathogens cause devastating diseases on economically and ecologically important plant species, threatening global food security, and causing billions of dollars of losses annually. During the infection process, pathogens secrete so-called effectors that support host colonization, often by deregulating host immune responses. Over the last decades, much of the research on molecular plant-microbe interactions has focused on the identification and functional characterization of such effectors. The increasing availability of sequenced plant pathogen genomes has enabled genomics-based discovery of effector candidates. Nevertheless, identification of full plant pathogen effector repertoires is often hampered by erroneous gene annotation and the localization effector genes in genomic regions that are notoriously difficult to assemble. Here, we argue that recent advances in genome sequencing technologies, genome assembly, gene annotation, as well as effector identification methods hold promise to disclose complete and correct effector repertoires. This allows to exploit complete effector repertoires, and knowledge of their diversity within pathogen populations, to develop durable and sustainable resistance breeding strategies, disease control, and management of plant pathogens.

Living apart together: crosstalk between the core and supernumerary genomes in a fungal plant pathogen

BACKGROUND

Eukaryotes display remarkable genome plasticity, which can include supernumerary chromosomes that differ markedly from the core chromosomes. Despite the widespread occurrence of supernumerary chromosomes in fungi, their origin, relation to the core genome and the reason for their divergent characteristics are still largely unknown. The complexity of genome assembly due to the presence of repetitive DNA partially accounts for this.

RESULTS

Here we use single-molecule real-time (SMRT) sequencing to assemble the genome of a prominent fungal wheat pathogen, Fusarium poae, including at least one supernumerary chromosome. The core genome contains limited transposable elements (TEs) and no gene duplications, while the supernumerary genome holds up to 25 % TEs and multiple gene duplications. The core genome shows all hallmarks of repeat-induced point mutation (RIP), a defense mechanism against TEs, specific for fungi. The absence of RIP on the supernumerary genome accounts for the differences between the two (sub)genomes, and results in a functional crosstalk between them. The supernumerary genome is a reservoir for TEs that migrate to the core genome, and even large blocks of supernumerary sequence (>200 kb) have recently translocated to the core. Vice versa, the supernumerary genome acts as a refuge for genes that are duplicated from the core genome.

CONCLUSIONS

For the first time, a mechanism was determined that explains the differences that exist between the core and supernumerary genome in fungi. Different biology rather than origin was shown to be responsible. A "living apart together" crosstalk exists between the core and supernumerary genome, accelerating chromosomal and organismal evolution.

Transposons passively and actively contribute to evolution of the two-speed genome of a fungal pathogen

Genomic plasticity enables adaptation to changing environments, which is especially relevant for pathogens that engage in “arms races” with their hosts. In many pathogens, genes mediating virulence cluster in highly variable, transposon-rich, physically distinct genomic compartments. However, understanding of the evolution of these compartments, and the role of transposons therein, remains limited. Here, we show that transposons are the major driving force for adaptive genome evolution in the fungal plant pathogen Verticillium dahliae. We show that highly variable lineage-specific (LS) regions evolved by genomic rearrangements that are mediated by erroneous double-strand repair, often utilizing transposons. We furthermore show that recent genetic duplications are enhanced in LS regions, against an older episode of duplication events. Finally, LS regions are enriched in active transposons, which contribute to local genome plasticity. Thus, we provide evidence for genome shaping by transposons, both in an active and passive manner, which impacts the evolution of pathogen virulence.

Extensive Mobilome-Driven Genome Diversification in Mouse Gut-Associated Bacteroides vulgatus mpk

Like many other Bacteroides species, Bacteroides vulgatus strain mpk, a mouse fecal isolate which was shown to promote intestinal homeostasis, utilizes a variety of mobile elements for genome evolution. Based on sequences collected by Pacific Biosciences SMRT sequencing technology, we discuss the challenges of assembling and studying a bacterial genome of high plasticity. Additionally, we conducted comparative genomics comparing this commensal strain with the B. vulgatus type strain ATCC 8482 as well as multiple other Bacteroides and Parabacteroides strains to reveal the most important differences and identify the unique features of B. vulgatus mpk. The genome of B. vulgatus mpk harbors a large and diverse set of mobile element proteins compared with other sequenced Bacteroides strains. We found evidence of a number of different horizontal gene transfer events and a genome landscape that has been extensively altered by different mobilization events. A CRISPR/Cas system could be identified that provides a possible mechanism for preventing the integration of invading external DNA. We propose that the high genome plasticity and the introduced genome instabilities of B. vulgatus mpk arising from the various mobilization events might play an important role not only in its adaptation to the challenging intestinal environment in general, but also in its ability to interact with the gut microbiota.

Mind the gap; seven reasons to close fragmented genome assemblies

Like other domains of life, research into the biology of filamentous microbes has greatly benefited from the advent of whole-genome sequencing. Next-generation sequencing (NGS) technologies have revolutionized sequencing, making genomic sciences accessible to many academic laboratories including those that study non-model organisms. Thus, hundreds of fungal genomes have been sequenced and are publically available today, although these initiatives have typically yielded considerably fragmented genome assemblies that often lack large contiguous genomic regions. Many important genomic features are contained in intergenic DNA that is often missing in current genome assemblies, and recent studies underscore the significance of non-coding regions and repetitive elements for the life style, adaptability and evolution of many organisms. The study of particular types of genetic elements, such as telomeres, centromeres, repetitive elements, effectors, and clusters of co-regulated genes, but also of phenomena such as structural rearrangements, genome compartmentalization and epigenetics, greatly benefits from having a contiguous and high-quality, preferably even complete and gapless, genome assembly. Here we discuss a number of important reasons to produce gapless, finished, genome assemblies to help answer important biological questions.