Search

Plant pathogens secrete effector molecules during host invasion to promote colonization. However, some of these effectors become recognized by host receptors to mount a defence response and establish immunity. Recently, a novel resistance was identified in wild tomato, mediated by the single dominant V2 locus, to control strains of the soil-borne vascular wilt fungus Verticillium dahliae that belong to race 2. With comparative genomics of race 2 strains and resistance-breaking race 3 strains, we identified the avirulence effector that activates V2 resistance, termed Av2. We identified 277 kb of race 2-specific sequence comprising only two genes encoding predicted secreted proteins that are expressed during tomato colonization. Subsequent functional analysis based on genetic complementation into race 3 isolates and targeted deletion from the race 1 isolate JR2 and race 2 isolate TO22 confirmed that one of the two candidates encodes the avirulence effector Av2 that is recognized in V2 tomato plants. Two Av2 allelic variants were identified that encode Av2 variants that differ by a single acid. Thus far, a role in virulence could not be demonstrated for either of the two variants.

Transposable elements (TEs) are a major source of genetic and regulatory variation in their host genome and are consequently thought to play important roles in evolution. Many fungal and oomycete plant pathogens have evolved dynamic and TE-rich genomic regions containing genes that are implicated in host colonization. TEs embedded in these regions have typically been thought to accelerate the evolution of these genomic compartments, but little is known about their dynamics in strains that harbor them. Here, we used whole-genome sequencing data of 42 strains of the fungal plant pathogen Verticillium dahliae to systematically identify polymorphic TEs that may be implicated in genomic as well as in gene expression variation. We identified 2,523 TE polymorphisms and characterize a subset of 8% of the TEs as dynamic elements that are evolutionary younger, less methylated, and more highly expressed when compared with the remaining 92% of the TE complement. As expected, the dynamic TEs are enriched in the dynamic genomic regions. Besides, we observed an association of dynamic TEs with pathogenicity-related genes that localize nearby and that display high expression levels. Collectively, our analyses demonstrate that TE dynamics in V. dahliae contributes to genomic variation, correlates with expression of pathogenicity-related genes, and potentially impacts the evolution of dynamic genomic regions.Significance statementTransposable elements (TEs) are ubiquitous components of genomes and are major sources of genetic and regulatory variation. Many plant pathogens have evolved TE-rich genomic regions containing genes with roles in host colonization, and TEs are thought to contribute to accelerated evolution of these dynamic regions. We analyzed the fungal plant pathogen Verticillium dahliae to identify TE variation between strains and to demonstrate that polymorphic TEs have specific characteristic that separates them from the majority of TEs. Polymorphic TEs are enriched in dynamic genomic regions and are associated with structural variants and highly expressed pathogenicity-related genes. Collectively, our results provide evidence for the hypothesis that dynamic TEs contribute to increased genomic diversity, functional variation, and the evolution of dynamic genomic regions.

The socio-economic impact of Fusarium odoratissimum, which is colloquially called tropical race 4 (TR4), is escalating as this fungal pathogen spreads to new banana-growing areas. Hence, the development of simple, reliable and rapid detection technologies is indispensable for implementing quarantine measures. Here, a versatile loop-mediated isothermal amplification (LAMP) assay has been developed that is applicable under field and laboratory conditions. DNA markers unique to TR4 isolates were obtained by diversity arrays technology sequencing (DArTseq), a genotyping by sequencing technology that was conducted on 27 genotypes, comprising 24 previously reported vegetative compatibility groups (VCGs) and three TR4 isolates. The developed LAMP TR4 assay was successfully tested using 22 TR4 isolates and 45 non-target fungal and bacterial isolates, as well as on infected plants under greenhouse and field conditions. The detection limit was 1 pg µL−1 pure TR4 DNA or 102 copies plasmid-localized TR4 unique sequence (SeqA) per reaction, which was not affected by background DNA in complex samples. The LAMP TR4 assay offers a powerful tool for the routine and unambiguous identification of banana plants infected with TR4, contributing to advanced diagnosis in field situations and monitoring of fusarium wilt.

Differential growth conditions typically trigger global transcriptional responses in filamentous fungi. Such fungal responses to environmental cues involve epigenetic regulation, including chemical histone modifications. It has been proposed that conditionally expressed genes, such as those that encode secondary metabolites but also effectors in pathogenic species, are often associated with a specific histone modification, lysine27 methylation of H3 (H3K27me3). However, thus far no analyses on the global H3K27me3 profiles have been reported under differential growth conditions in order to assess if H3K27me3 dynamics governs differential transcriptional. Using ChIP- and RNA-sequencing data from the plant pathogenic fungus Verticillium dahliae grown on three in vitro cultivation media, we now show that a substantial number of the identified H3K27me3 domains globally display stable profiles among these growth conditions. However, we do observe local quantitative differences in H3K27me3 ChIP-seq signal that associate with differentially transcribed genes between media. Comparing the in vitro results to expression during plant infection suggests that in planta induced genes require chromatin remodelling to achieve expression. Overall, our results demonstrate that although some loci display H3K27me3 dynamics that can contribute to transcriptional variation, other loci do not show such dependence. Thus, we conclude that while H3K27me3 is required for transcriptional repression, it is not a conditionally responsive global regulator of differential transcription. We propose that the H3K27me3 domains that do not undergo dynamic methylation may contribute to transcription through other mechanisms or may serve additional genomic regulatory functions. IMPORTANCE In many organisms, including filamentous fungi, epigenetic mechanisms that involve chemical and physical modifications of DNA without changing the genetic sequence have been implicated in transcriptional responses upon developmental or environmental cues. In fungi, facultative heterochromatin that can de-condense to allow transcription in response to developmental changes or to environmental stimuli is characterized by tri-methylation of lysine 27 on histone H3 (H3K27me3), and H3K27me3 has been implicated in transcriptional regulation although precise mechanisms and functions remain enigmatic. Based on ChIP- and RNA-sequencing data, we show for the soil-borne broad host range vascular wilt plant pathogenic fungus Verticillium dahliae that although some loci display H3K27me3 dynamics that can contribute to transcriptional variation, other loci do not show such dependence. Thus, although we recognize that H3K27me3 is required for transcriptional repression, we also conclude that this mark is not a conditionally responsive global regulator of differential transcription upon responses to environmental cues.

Background DNA methylation is an important epigenetic control mechanism that in many fungi is restricted to genomic regions containing transposable elements (TEs). Two DNA methyltransferases, Dim2 and Dnmt5, are known to perform methylation at cytosines in fungi. While most ascomycete fungi encode both Dim2 and Dnmt5, only few functional studies have been performed in species containing both. Methods In this study, we report functional analysis of both Dim2 and Dnmt5 in the plant pathogenic fungus Verticillium dahliae. Results Our results show that Dim2, but not Dnmt5 or the putative sexual-cycle-related DNA methyltransferase Rid, is responsible for the majority of DNA methylation under the tested conditions. Single or double DNA methyltransferase mutants did not show altered development, virulence, or transcription of genes or TEs. In contrast, Hp1 and Dim5 mutants that are impacted in chromatin-associated processes upstream of DNA methylation are severely affected in development and virulence and display transcriptional reprogramming in specific hypervariable genomic regions (so-called adaptive genomic regions) that contain genes associated with host colonization. As these adaptive genomic regions are largely devoid of DNA methylation and of Hp1- and Dim5-associated heterochromatin, the differential transcription is likely caused by pleiotropic effects rather than by differential DNA methylation. Conclusion Overall, our study suggests that Dim2 is the main DNA methyltransferase in V. dahliae and, in conjunction with work on other fungi, is likely the main active DNMT in ascomycetes, irrespective of Dnmt5 presence. We speculate that Dnmt5 and Rid act under specific, presently enigmatic, conditions or, alternatively, act in DNA-associated processes other than DNA methylation.

DNA methylation is an important epigenetic control mechanism that in many fungi is restricted to genomic regions containing transposons. Two DNA methyltransferases, Dim2 and Dnmt5, are known to perform methylation at cytosines in fungi. While most ascomycete fungi encode both Dim2 and Dnmt5, only few functional studies have been performed in species containing both. In this study, we report functional analysis of bothDim2andDnmt5in the plant pathogenic fungusVerticillium dahliae. Our results show that Dim2, but not Dnmt5 or the putative sexual-cycle related DNA methyltransferase Rid, is responsible for nearly all DNA methylation. Single or double DNA methyltransferase mutants did not show altered development, virulence, or transcription of genes or transposons. In contrast,Hp1andDim5mutants that are impacted in chromatin-associated processes upstream of DNA methylation are severely affected in development and virulence and display extensive transcriptional reprogramming in specific hypervariable genomic regions (so-called lineage-specific (LS) regions) that contain genes associated with host colonization. As these LS regions are largely devoid of DNA methylation and of Hp1- and Dim5-associated heterochromatin, the differential transcription is likely caused by pleiotropic effects rather than by differential DNA methylation. Overall, our study suggests that Dim2 is the main DNA methyltransferase inV. dahliaeand, in conjunction with work on other fungi, is likely the main active DNMT in ascomycetes, irrespective ofDnmt5presence. We speculate that Dnmt5 acts under specific, presently enigmatic, conditions or, alternatively, acts in DNA-associated processes other than DNA methylation.

SUMMARYPlant pathogens secrete effector molecules during host invasion to promote host colonization. However, some of these effectors become recognized by host receptors, encoded by resistance genes, to mount defense response and establish immunity. Recently, a novel resistance was identified in tomato, mediated by the single dominant V2 locus, to control strains of the soil-borne vascular wilt fungus Verticillium dahliae that belong to race 2. We performed comparative genomics between race 2 strains and resistance-breaking race 3 strains to identify the avirulence effector that activates V2 resistance, termed Av2. We identified 277 kb of race 2-specific sequence comprising only two genes that encode predicted secreted proteins, both of which are expressed by V. dahliae during tomato colonization. Subsequent functional analysis based on genetic complementation into race 3 isolates confirmed that one of the two candidates encodes the avirulence effector Av2 that is recognized in V2 tomato plants. The identification of Av2 will not only be helpful to select tomato cultivars that are resistant to race 2 strains of V. dahliae, as the corresponding V2 resistance gene has not yet been mapped, but also to monitor adaptations in the V. dahliae population to deployment of V2-containing tomato cultivars in agriculture.

The fungal genus Verticillium contains ten species, some of which are notorious plant pathogens causing vascular wilt diseases in host plants, while others are known as saprophytes and opportunistic plant pathogens. Whereas the genome of V. dahliae, the most notorious plant pathogen of the genus, has been well characterized, evolution and speciation of other members of the genus received little attention thus far. Here, we sequenced the genomes of the nine haploid Verticillium spp. to study evolutionary trajectories of their divergence from a last common ancestor. Frequent occurrence of chromosomal rearrangement and gene family loss was identified. In addition to ∼11 000 genes that are shared at least between two species, only 200-600 species-specific genes occur. Intriguingly, these species-specific genes show different features than the shared genes.

Transposable elements are powerful drivers of genome evolution in many eukaryotes. Although they are mostly considered as 'selfish' genetic elements, increasing evidence suggests that they contribute to genetic variability; particularly under stress conditions. Over the past few years, the role of transposable elements during host-microbe interactions has been recognised. It has been proposed that many pathogenic microbes have evolved a 'two-speed' genome with regions that show increased variability and that are enriched in transposable elements and pathogenicity-related genes. Plants similarly display structured genomes with transposable-element-rich regions that mediate accelerated evolution. Immune receptor genes typically reside in such regions. Various mechanisms have recently been identified through which transposable elements contribute to the coevolution between plants and their associated microbes. Transposable elements are powerful drivers of adaptive genome evolution in plants and in symbiotic microbes, and contribute to their coevolution.Analysis of large next-generation sequencing datasets fuels a better understanding of the precise role of transposable elements in genome evolution, revealing active as well as passive contributions.Transposable elements contribute to genome evolution through diverse mechanisms, by mediating structural variations, gene inactivation, gene copy variation, but also by affecting gene expression.Further mechanistic understanding of the role of transposable elements in genome evolution will come from detailed structural analysis of chromatin.

Plant pathogens continuously evolve to evade host immune responses. During host colonization, many fungal pathogens secrete effectors to perturb such responses, but these in turn may become recognized by host immune receptors. To facilitate the evolution of effector repertoires, such as the elimination of recognized effectors, effector genes often reside in genomic regions that display increased plasticity, a phenomenon that is captured in the two‐speed genome hypothesis. The genome of the vascular wilt fungus Verticillium dahliae displays regions with extensive presence/absence polymorphisms, so‐called lineage‐specific regions, that are enriched in in planta‐induced putative effector genes. As expected, comparative genomics reveals differential degrees of sequence divergence between lineage‐specific regions and the core genome. Unanticipated, lineage‐specific regions display markedly higher sequence conservation in coding as well as noncoding regions than the core genome. We provide evidence that disqualifies horizontal transfer to explain the observed sequence conservation and conclude that sequence divergence occurs at a slower pace in lineage‐specific regions of the V. dahliae genome. We hypothesize that differences in chromatin organisation may explain lower nucleotide substitution rates in the plastic, lineage‐specific regions of V. dahliae.

Along withPlasmopara destructor, Peronosopora belbahriihas arguably been the economically most important newly emerging downy mildew pathogen of the past two decades. Originating from Africa, it has started devastating basil production throughout the world, most likely due to the distribution of infested seed material. Here we present the genome of this pathogen and results from comparisons of its genomic features to other oomycetes. The assembly of the nuclear genome was ca. 35.4 Mbp in length, with an N50 scaffold length of ca. 248 kbp and an L50 scaffold count of 46. The circular mitochondrial genome consisted of ca. 40.1 kbp. From the repeat-masked genome 9049 protein-coding genes were predicted, out of which 335 were predicted to have extracellular functions, representing the smallest secretome so far found in peronosporalean oomycetes. About 16 % of the genome consists of repetitive sequences, and based on simple sequence repeat regions, we provide a set of microsatellites that could be used for population genetic studies ofPe. belbahrii. Peronospora belbahriihas undergone a high degree of convergent evolution, reflecting its obligate biotrophic lifestyle. Features of its secretome, signalling networks, and promoters are presented, and some patterns are hypothesised to reflect the high degree of host specificity inPeronosporaspecies. In addition, we suggest the presence of additional virulence factors apart from classical effector classes that are promising candidates for future functional studies.

The oomycete pathogen Phytophthora infestans causes potato and tomato late blight, a disease that is a serious threat to agriculture. P. infestans is a hemibiotrophic pathogen, and during infection, it scavenges nutrients from living host cells for its own proliferation. To date, the nutrient flux from host to pathogen during infection has hardly been studied, and the interlinked metabolisms of the pathogen and host remain poorly understood. Here, we reconstructed an integrated metabolic model of P. infestans and tomato (Solanum lycopersicum) by integrating two previously published models for both species. We used this integrated model to simulate metabolic fluxes from host to pathogen and explored the topology of the model to study the dependencies of the metabolism of P. infestans on that of tomato. This showed, for example, that P. infestans, a thiamine auxotroph, depends on certain metabolic reactions of the tomato thiamine biosynthesis. We also exploited dual-transcriptome data of a time course of a full late blight infection cycle on tomato leaves and integrated the expression of metabolic enzymes in the model. This revealed profound changes in pathogen-host metabolism during infection. As infection progresses, P. infestans performs less de novo synthesis of metabolites and scavenges more metabolites from tomato. This integrated metabolic model for the P. infestans-tomato interaction provides a framework to integrate data and generate hypotheses about in planta nutrition of P. infestans throughout its infection cycle. IMPORTANCE Late blight disease caused by the oomycete pathogen Phytophthora infestans leads to extensive yield losses in tomato and potato cultivation worldwide. To effectively control this pathogen, a thorough understanding of the mechanisms shaping the interaction with its hosts is paramount. While considerable work has focused on exploring host defense mechanisms and identifying P. infestans proteins contributing to virulence and pathogenicity, the nutritional strategies of the pathogen are mostly unresolved. Genome-scale metabolic models (GEMs) can be used to simulate metabolic fluxes and help in unravelling the complex nature of metabolism. We integrated a GEM of tomato with a GEM of P. infestans to simulate the metabolic fluxes that occur during infection. This yields insights into the nutrients that P. infestans obtains during different phases of the infection cycle and helps in generating hypotheses about nutrition in planta.

Horizontal gene transfer (HGT) involves the transmission of genetic material between distinct evolutionary lineages and can be an important source of biological innovation. Reports of inter-kingdom HGT to eukaryotic microbial pathogens have accumulated over recent years. Verticillium dahliae is a notorious plant pathogen that causes vascular wilt disease on hundreds of plant species, resulting in high economic losses every year. Previously, the effector gene Ave1 and a glucosyltransferase-encoding gene were identified as virulence factor-encoding genes that were proposed to be horizontally acquired from a plant and a bacterial donor, respectively. However, to what extent HGT contributed to the overall genome composition of V. dahliae remained elusive. Here, we systematically searched for evidence of inter-kingdom HGT events in the genome of V. dahliae and provide evidence for extensive horizontal gene acquisition from bacterial origin.

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.

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.

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.

Contains fulltext : 168245.pdf (Publisher’s version ) (Open Access) 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.

Fungi are renowned producers of natural compounds, also known as secondary metabolites (SMs) that display a wide array of biological activities. Typically, the genes that are involved in the biosynthesis of SMs are located in close proximity to each other in so-called secondary metabolite clusters. Many plant-pathogenic fungi secrete SMs during infection in order to promote disease establishment, for instance as cytocoxic compounds. Verticillium dahliae is a notorious plant pathogen that can infect over 200 host plants worldwide. However, the SM repertoire of this vascular pathogen remains mostly uncharted. To unravel the potential of V. dahliae to produce SMs, we performed in silico predictions and in-depth analyses of its secondary metabolite clusters. Using distinctive traits of gene clusters and the conserved signatures of core genes 25 potential SM gene clusters were identified. Subsequently, phylogenetic and comparative genomics analyses were performed, revealing that two putative siderophores, ferricrocin and TAFC, DHN-melanin and fujikurin may belong to the SM repertoire of V. dahliae., In total, 25 potential SMCs were identified in the genome of the plant pathogenic fungus Verticillium dahliae, including loci that can be implicated in DHN-melanin, ferricrocin, triacetyl fusarinine and fujikurin production.

Hybridization is an important evolutionary mechanism that can enable organisms to adapt to environmental challenges. It has previously been shown that the fungal allodiploid species Verticillium longisporum, causal agent of Verticillium stem striping in rape seed, has originated from at least three independent hybridization events between two haploid Verticillium species. To reveal the impact of genome duplication as a consequence of the hybridization, we studied the genome and transcriptome dynamics upon two independent V. longisporum hybridization events, represented by the hybrid lineages “A1/D1” and “A1/D3”. We show that the V. longisporum genomes are characterized by extensive chromosomal rearrangements, including between parental chromosomal sets. V. longisporum hybrids display signs of evolutionary dynamics that are typically associated with the aftermath of allodiploidization, such as haploidization and a more relaxed gene evolution. Expression patterns of the two sub-genomes within the two hybrid lineages are more similar than those of the shared A1 parent between the two lineages, showing that expression patterns of the parental genomes homogenized within a lineage. However, as genes that display differential parental expression in planta do not typically display the same pattern in vitro, we conclude that sub-genome-specific responses occur in both lineages. Overall, our study uncovers the genomic and transcriptomic plasticity during evolution of the filamentous fungal hybrid V. longisporum and illustrate its adaptive potential.ImportanceVerticillium is a genus of plant-associated fungi that include a handful of plant pathogens that collectively affect a wide range of hosts. On several occasions, haploid Verticillium species hybridized into the stable allodiploid species Verticillium longisporum, which is, in contrast to haploid Verticillium species, a Brassicaceae specialist. Here, we studied the evolutionary genome and transcriptome dynamics of V. longisporum and the impact of the hybridization. V. longisporum genomes display a mosaic structure due do genomic rearrangements between the parental chromosome sets. Similar to other allopolyploid hybrids, V. longisporum displays an ongoing loss of heterozygosity and a more relaxed gene evolution. Also, differential parental gene expression is observed, with an enrichment for genes that encode secreted proteins. Intriguingly, the majority of these genes displays sub-genome-specific responses under differential growth conditions. In conclusion, hybridization has incited the genomic and transcriptomic plasticity that enables adaptation to environmental changes in a parental allele-specific fashion.

Selection pressure impacts genomes unevenly, as different genes adapt with differential speed to establish an organism’s optimal fitness. Plant pathogens co-evolve with their hosts, which implies continuously adaption to evade host immunity. Effectors are secreted proteins that mediate immunity evasion, but may also typically become recognized by host immune receptors. To facilitate effector repertoire alterations, in many pathogens, effector genes reside in dynamic genomic regions that are thought to display accelerated evolution, a phenomenon that is captured by the two-speed genome hypothesis. The genome of the vascular wilt pathogen Verticillium dahliae has been proposed to obey to a similar two-speed regime with dynamic, lineage-specific regions that are characterized by genomic rearrangements, increased transposable element activity and enrichment in in planta-induced effector genes. However, little is known of the origin of, and sequence diversification within, these lineage-specific regions. Based on comparative genomics among Verticillium spp. we now show differential sequence divergence between core and lineage-specific genomic regions of V. dahliae. Surprisingly, we observed that lineage-specific regions display markedly increased sequence conservation. Since single nucleotide diversity is reduced in these regions, host adaptation seems to be merely achieved through presence/absence polymorphisms. Increased sequence conservation of genomic regions important for pathogenicity is an unprecedented finding for filamentous plant pathogens and signifies the diversity of genomic dynamics in host-pathogen co-evolution.

Allopolyploidization, genome duplication through interspecific hybridization, is an important evolutionary mechanism that can enable organisms to adapt to environmental changes or stresses. This increased adaptive potential of allopolyploids can be particularly relevant for plant pathogens in their quest for host immune response evasion. Allodiploidization likely caused the shift in host range of the fungal pathogen plant Verticillium longisporum, as V. longisporum mainly infects Brassicaceae plants in contrast to haploid Verticillium spp. In this study, we investigated the allodiploid genome structure of V. longisporum and its evolution in the hybridization aftermath. The nuclear genome of V. longisporum displays a mosaic structure, as numerous contigs consists of sections of both parental origins. V. longisporum encountered extensive genome rearrangements, whereas the contribution of gene conversion is negligible. Thus, the mosaic genome structure mainly resulted from genomic rearrangements between parental chromosome sets. Furthermore, a mosaic structure was also found in the mitochondrial genome, demonstrating its bi-parental inheritance. In conclusion, the nuclear and mitochondrial genomes of V. longisporum parents interacted dynamically in the hybridization aftermath. Conceivably, novel combinations of DNA sequence of different parental origin facilitated genome stability after hybridization and consecutive niche adaptation of V. longisporum.

The Alternaria genus consists of saprophytic fungi as well as plant-pathogenic species that have significant economic impact. To date, the genomes of multiple Alternaria species have been sequenced. These studies have yielded valuable data for molecular studies on Alternaria fungi. However, most of the current Alternaria genome assemblies are highly fragmented, thereby hampering the identification of genes that are involved in causing disease. Here, we report a gapless genome assembly of A. solani, the causal agent of early blight in tomato and potato. The genome assembly is a significant step toward a better understanding of pathogenicity of A. solani.

SUMMARYThe fungal genus Verticillium contains ten species, some of which are notorious plant pathogens causing vascular wilt diseases in host plants, while others are known as saprophytes and opportunistic plant pathogens. Whereas the genome of V. dahliae, the most notorious plan pathogen of the genus, has been well characterized, evolution and speciation of other members of the genus received little attention thus far. Here, we sequenced the genomes of the nine haploid Verticillium spp. to study evolutionary trajectories of their divergence from a last common ancestor. Frequent occurrence of chromosomal rearrangement and gene family loss was identified. In addition to ~11,000 core genes that are shared among all species, only 200-600 species-specific genes occur. Intriguingly, these species-specific genes show different features than core genes.

A DNA modification - methylation of cytosines and adenines - has important roles in diverse processes such as regulation of gene expression and genome stability, yet until recently adenine methylation had been considered to be only a hallmark of prokaryotes. A new study identifies abundant adenine methylation of transcriptionally active genes in early-diverging fungi that, together with recent other work, emphasizes the importance of adenine methylation in eukaryotes.

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.

All species continuously evolve to adapt to changing environments. The genetic variation that fosters such adaptation is caused by a plethora of mechanisms, including meiotic recombination that generates novel allelic combinations in the progeny of two parental lineages. However, a considerable number of eukaryotic species, including many fungi, do not have an apparent sexual cycle and are consequently thought to be limited in their evolutionary potential. As such organisms are expected to have reduced capability to eliminate deleterious mutations, they are often considered as evolutionary dead ends. However, inspired by recent reports we argue that such organisms can be as persistent as organisms with conventional sexual cycles through the use of other mechanisms, such as genomic rearrangements, to foster adaptation.

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 aggressiveness 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. We now show that transposons are the major driving force for adaptive genome evolution in the fungal plant pathogen Verticillium dahliae. Highly variable genomic regions evolved by frequent segmental duplications mediated by erroneous homologous recombination, often utilizing transposons, leading to genetic material that is free to diverge. Intriguingly, the duplicated regions are enriched in active transposons that further 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 aggressiveness.

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.

An organism’s adaptation to changing environments is fueled by its genetic variability, which is established by mechanisms ranging from single-nucleotide polymorphisms to large-scale structural variations, all of which affect chromosomal shape, organization, and gene content [1]. These processes are particularly relevant for pathogens that must respond to continual selection pressure arising from host immune systems that evolved to detect the presence or activity of potentialmicrobial pathogens through a variety of invasion patterns [2]. In their adaptive response, pathogens evolve strategies, often involving secreted effectormolecules, to overcome host immunity and support host colonization [3]. Thus, it can be anticipated that this coevolutionary arms race leads to highly specific interactions between adapted pathogens and their specific hosts. Paradoxically, particular pathogens successfully colonize a broad range of hosts, yet how such pathogens cope in arms races with such a diversity of hosts remains unknown.

Lettuce downy mildew (Bremia lactucae) is a rapidly adapting oomycete pathogen affecting commercial lettuce cultivation. Oomycetes are known to use a diverse arsenal of secreted proteins (effectors) to manipulate their hosts. Two classes of effector are known to be translocated by the host: the RXLRs and Crinklers. To gain insight into the repertoire of effectors used by B. lactucae to manipulate its host, we performed massively parallel sequencing of cDNA derived from B. lactucae spores and infected lettuce (Lactuca sativa) seedlings. From over 2.3 million 454 GS FLX reads, 59 618 contigs were assembled representing both plant and pathogen transcripts. Of these, 19 663 contigs were determined to be of B. lactucae origin as they matched pathogen genome sequences (SOLiD) that were obtained from >270 million reads of spore-derived genomic DNA. After correction of cDNA sequencing errors with SOLiD data, translation into protein models and filtering, 16 372 protein models remained, 1023 of which were predicted to be secreted. This secretome included elicitins, necrosis and ethylene-inducing peptide 1-like proteins, glucanase inhibitors and lectins, and was enriched in cysteine-rich proteins. Candidate host-translocated effectors included 78 protein models with RXLR effector features. In addition, we found indications for an unknown number of Crinkler-like sequences. Similarity clustering of secreted proteins revealed additional effector candidates. We provide a first look at the transcriptome of B. lactucae and its encoded effector arsenal.

This article deals with: Bananas: their origin and global rollout; genetic diversity of Fusarium oxysporum f.sp. cubense, the causal agent of Panama Disease; Panama Disease: history repeats itself; tropical race 4, a single pathogen clone, threatens global banana production; strategies for sustainable Panama Disease management.

An unknown fungus has been isolated as a contaminant of in vitro -grown fungal cultures. In an attempt to identify the contamination, we isolated the causal agent and performed whole-genome sequencing. BLAST analysis of the internal transcribed spacer (ITS) sequence against the NCBI database showed 100% identity to Trichoderma atroviride , and further alignment of the genome assembly confirmed the unknown fungus to be T. atroviride . Here, we report the draft genome sequence of a T. atroviride strain.

The cosmopolitan fungal genusAlternariaconsists of multiple saprophytic and pathogenic species. Based on phylogenetic and morphological studies, the genus is currently divided into 26 sections.Alternariasect.Alternariacontains most of the small-sporedAlternariaspecies with concatenated conidia, including important plant, human and postharvest pathogens. Species within sect.Alternariahave been mostly described based on morphology and / or host-specificity, yet molecular variation between them is minimal. To investigate whether the described morphospecies within sect.Alternariaare supported by molecular data, whole-genome sequencing of nineAlternariamorphospecies supplemented with transcriptome sequencing of 12Alternariamorphospecies as well as multi-gene sequencing of 168Alternariaisolates was performed. The assembled genomes ranged in size from 33.3–35.2 Mb within sect.Alternariaand from 32.0–39.1 Mb for allAlternariagenomes. The number of repetitive sequences differed significantly between the differentAlternariagenomes; ranging from 1.4–16.5 %. The repeat content within sect.Alternariawas relatively low with only 1.4–2.7 % of repeats. Whole-genome alignments revealed 96.7–98.2 % genome identity between sect.Alternariaisolates, compared to 85.1–89.3 % genome identity for isolates from other sections to theA. alternatareference genome. Similarly, 1.4–2.8 % and 0.8–1.8 % single nucleotide polymorphisms (SNPs) were observed in genomic and transcriptomic sequences, respectively, between isolates from sect.Alternaria, while the percentage of SNPs found in isolates from different sections compared to theA. alternatareference genome was considerably higher; 8.0–10.3 % and 6.1–8.5 %. The topology of a phylogenetic tree based on the whole-genome and transcriptome reads was congruent with multi-gene phylogenies based on commonly used gene regions. Based on the genome and transcriptome data, a set of core proteins was extracted, and primers were designed on two gene regions with a relatively low degree of conservation within sect.Alternaria(96.8 and 97.3 % conservation). Their potential discriminatory power within sect.Alternariawas tested next to nine commonly used gene regions in sect.Alternaria, namely the SSU, LSU, ITS,gapdh,rpb2,tef1,Alt a 1,endoPGand OPA10-2 gene regions. The phylogenies from the two gene regions with a relatively low conservation, KOG1058 and KOG1077, could not distinguish the described morphospecies within sect.Alternariamore effectively than the phylogenies based on the commonly used gene regions forAlternaria. Based on genome and transcriptome comparisons and molecular phylogenies,Alternariasect.Alternariaconsists of only 11 phylogenetic species and one species complex. Thirty-five morphospecies, which cannot be distinguished based on the multi-gene phylogeny, are synonymised underA. alternata. By providing guidelines for the naming and identification of phylogenetic species inAlternariasect.Alternaria, this manuscript provides a clear and stable species classification in this section.

Cytokinins represent a group of plant hormones that have been shown to be essential for plant growth and development. A recent large-scale phylogenetic analysis of components of the cytokinin signal transduction pathway revealed, among other findings, the existence of a second, previously unknown subfamily of cytokinin receptors. Here we report that the cytokinin binding domains of the members of the 2 subfamilies contain residues that are highly conserved in either or in both subfamilies. Experiments using fluorescence microscopy hint at an ER and a plasma membrane localization for 2 members of the newly identified subfamily. These data provide new insights in the conservation of sequence and localization properties among the 2 subfamilies.

Next-generation sequencing (NGS) technologies have increased the scalability, speed, and resolution of genomic sequencing and, thus, have revolutionized genomic studies. However, eukaryotic genome sequencing initiatives typically yield considerably fragmented genome assemblies. Here, we assessed various state-of-the-art sequencing and assembly strategies in order to produce a contiguous and complete eukaryotic genome assembly, focusing on the filamentous fungus Verticillium dahliae. Compared with Illumina-based assemblies of the V. dahliae genome, hybrid assemblies that also include PacBio-generated long reads establish superior contiguity. Intriguingly, provided that sufficient sequence depth is reached, assemblies solely based on PacBio reads outperform hybrid assemblies and even result in fully assembled chromosomes. Furthermore, the addition of optical map data allowed us to produce a gapless and complete V. dahliae genome assembly of the expected eight chromosomes from telomere to telomere. Consequently, we can now study genomic regions that were previously not assembled or poorly assembled, including regions that are populated by repetitive sequences, such as transposons, allowing us to fully appreciate an organism’s biological complexity. Our data show that a combination of PacBio-generated long reads and optical mapping can be used to generate complete and gapless assemblies of fungal genomes., IMPORTANCE Studying whole-genome sequences has become an important aspect of biological research. The advent of next-generation sequencing (NGS) technologies has nowadays brought genomic science within reach of most research laboratories, including those that study nonmodel organisms. However, most genome sequencing initiatives typically yield (highly) fragmented genome assemblies. Nevertheless, considerable relevant information related to genome structure and evolution is likely hidden in those nonassembled regions. Here, we investigated a diverse set of strategies to obtain gapless genome assemblies, using the genome of a typical ascomycete fungus as the template. Eventually, we were able to show that a combination of PacBio-generated long reads and optical mapping yields a gapless telomere-to-telomere genome assembly, allowing in-depth genome analyses to facilitate functional studies into an organism’s biology.

Vascular wilts caused by Verticillium spp. are destructive plant diseases affecting hundreds of hosts. Only a few Verticillium spp. are causal agents of vascular wilt diseases, of which V. dahliae is the most notorious pathogen, and several V. dahliae genomes are available. In contrast, V. tricorpus is mainly known as a saprophyte and causal agent of opportunistic infections. Based on a hybrid approach that combines second and third generation sequencing, a near-gapless V. tricorpus genome assembly was obtained. With comparative genomics, we sought to identify genomic features in V. dahliae that confer the ability to cause vascular wilt disease. Unexpectedly, both species encode similar effector repertoires and share a genomic structure with genes encoding secreted proteins clustered in genomic islands. Intriguingly, V. tricorpus contains significantly fewer repetitive elements and an extended spectrum of secreted carbohydrate- active enzymes when compared with V. dahliae. In conclusion, we highlight the technical advances of a hybrid sequencing and assembly approach and show that the saprophyte V. tricorpus shares many hallmark features with the pathogen V. dahliae.

The gene composition of present-day genomes has been shaped by a complicated evolutionary history, resulting in diverse distributions of genes across genomes. The pattern of presence and absence of a gene in different genomes is called its phylogenetic profile. It has been shown that proteins whose encoding genes have highly similar profiles tend to be functionally related: As these genes were gained and lost together, their encoded proteins can probably only perform their full function if both are present. However, a large proportion of genes encoding interacting proteins do not have matching profiles. In this study, we analysed one possible reason for this, namely that phylogenetic profiles can be affected by multi-functional proteins such as shared subunits of two or more protein complexes. We found that by considering triplets of proteins, of which one protein is multi-functional, a large fraction of disturbed co-occurrence patterns can be explained., Author Summary Every genome of current day species contains a very unique selection of genes. Why a specific genome is composed of exactly those genes is determined by many factors, but often not resolvable. It seems plausible that interacting genes would either occur together or be absent together, because if one of them is alone, it might not be able to perform its function properly, just as a bolt can only perform its function together with a nut and vice versa. However, it turns out that interacting genes very often do not nicely co-occur across a wide range of species, and frequently one gene can be found but the other not. In this study, we investigated the co-occurrences of multi-functional proteins and found that they are often maintained in a genome, even if one of their interaction partners is lost. This is because they can still perform some functions with other interaction partners that are still present. We can show that this has a noticeable effect on genome compositions and can explain otherwise surprisingly mismatching co-occurrence patterns of interacting genes.

Plant infection by oomycete pathogens is a complex process. It requires precise expression of a plethora of genes in the pathogen that contribute to a successful interaction with the host. Whereas much effort has been made to uncover the molecular systems underlying this infection process, mechanisms of transcriptional regulation of the genes involved remain largely unknown. We performed the first systematic de-novo DNA motif discovery analysis in Phytophthora. To this end, we utilized the genome sequence of the late blight pathogen Phytophthora infestans and two related Phytophthora species (P. ramorum and P. sojae), as well as genome-wide in planta gene expression data to systematically predict 19 conserved DNA motifs. This catalog describes common eukaryotic promoter elements whose functionality is supported by the presence of orthologs of known general transcription factors. Together with strong functional enrichment of the common promoter elements towards effector genes involved in pathogenicity, we obtained a new and expanded picture of the promoter structure in P. infestans. More intriguingly, we identified specific DNA motifs that are either highly abundant or whose presence is significantly correlated with gene expression levels during infection. Several of these motifs are observed upstream of genes encoding transporters, RXLR effectors, but also transcriptional regulators. Motifs that are observed upstream of known pathogenicity-related genes are potentially important binding sites for transcription factors. Our analyses add substantial knowledge to the as of yet virtually unexplored question regarding general and specific gene regulation in this important class of pathogens. We propose hypotheses on the effects of cis-regulatory motifs on the gene regulation of pathogenicity-related genes and pinpoint motifs that are prime targets for further experimental validation.

Biotrophic plant pathogens secrete effector proteins that are important for infection of the host. The aim of this study was to identify effectors of the downy mildew pathogen Hyaloperonospora arabidopsidis (Hpa) that are expressed during infection of its natural host Arabidopsis thaliana. Infection-related transcripts were identified from Expressed Sequence Tags (ESTs) derived from leaves of the susceptible Arabidopsis Ws eds1-1 mutant inoculated with the highly virulent Hpa isolate Waco9. Assembly of 6364 ESTs yielded 3729 unigenes, of which 2164 were Hpa-derived. From the translated Hpa unigenes, 198 predicted secreted proteins were identified. Of these, 75 were found to be Hpa-specific and six isolate Waco9-specific. Among 42 putative effectors identified there were three Elicitin-like proteins, 16 Cysteine-rich proteins and 18 host-translocated RXLR effectors. Sequencing of alleles in different Hpa isolates revealed that five RXLR genes show signatures of diversifying selection. Thus, EST analysis of Hpa-infected Arabidopsis is proving to be a powerful method for identifying pathogen effector candidates expressed during infection. Delivery of the Waco9-specific protein RXLR29 in planta revealed that this effector can suppress PAMP-triggered immunity and enhance disease susceptibility. We propose that differences in host colonization can be conditioned by isolate-specific effectors.

Oomycetes comprise a diverse group of organisms that morphologically resemble fungi but belong to the stramenopile lineage within the supergroup of chromalveolates. Recent studies have shown that plant pathogenic oomycetes have expanded gene families that are possibly linked to their pathogenic lifestyle. We analyzed the protein domain organization of 67 eukaryotic species including four oomycete and five fungal plant pathogens. We detected 246 expanded domains in fungal and oomycete plant pathogens. The analysis of genes differentially expressed during infection revealed a significant enrichment of genes encoding expanded domains as well as signal peptides linking a substantial part of these genes to pathogenicity. Overrepresentation and clustering of domain abundance profiles revealed domains that might have important roles in host-pathogen interactions but, as yet, have not been linked to pathogenicity. The number of distinct domain combinations (bigrams) in oomycetes was significantly higher than in fungi. We identified 773 oomycete-specific bigrams, with the majority composed of domains common to eukaryotes. The analyses enabled us to link domain content to biological processes such as host-pathogen interaction, nutrient uptake, or suppression and elicitation of plant immune responses. Taken together, this study represents a comprehensive overview of the domain repertoire of fungal and oomycete plant pathogens and points to novel features like domain expansion and species-specific bigram types that could, at least partially, explain why oomycetes are such remarkable plant pathogens.

Oomycetes in the class Saprolegniomycetidae of the Eukaryotic kingdom Stramenopila have evolved as severe pathogens of amphibians, crustaceans, fish and insects, resulting in major losses in aquaculture and damage to aquatic ecosystems. We have sequenced the 63 Mb genome of the fresh water fish pathogen, Saprolegnia parasitica. Approximately 1/3 of the assembled genome exhibits loss of heterozygosity, indicating an efficient mechanism for revealing new variation. Comparison of S. parasitica with plant pathogenic oomycetes suggests that during evolution the host cellular environment has driven distinct patterns of gene expansion and loss in the genomes of plant and animal pathogens. S. parasitica possesses one of the largest repertoires of proteases (270) among eukaryotes that are deployed in waves at different points during infection as determined from RNA-Seq data. In contrast, despite being capable of living saprotrophically, parasitism has led to loss of inorganic nitrogen and sulfur assimilation pathways, strikingly similar to losses in obligate plant pathogenic oomycetes and fungi. The large gene families that are hallmarks of plant pathogenic oomycetes such as Phytophthora appear to be lacking in S. parasitica, including those encoding RXLR effectors, Crinkler's, and Necrosis Inducing-Like Proteins (NLP). S. parasitica also has a very large kinome of 543 kinases, 10% of which is induced upon infection. Moreover, S. parasitica encodes several genes typical of animals or animal-pathogens and lacking from other oomycetes, including disintegrins and galactose-binding lectins, whose expression and evolutionary origins implicate horizontal gene transfer in the evolution of animal pathogenesis in S. parasitica., Author Summary Fish are an increasingly important source of animal protein globally, with aquaculture production rising dramatically over the past decade. Saprolegnia is a fungal-like oomycete and one of the most destructive fish pathogens, causing millions of dollars in losses to the aquaculture industry annually. Saprolegnia has also been linked to a worldwide decline in wild fish and amphibian populations. Here we describe the genome sequence of the first animal pathogenic oomycete and compare the genome content with the available plant pathogenic oomycetes. We found that Saprolegnia lacks the large effector families that are hallmarks of plant pathogenic oomycetes, showing evolutionary adaptation to the host. Moreover, Saprolegnia harbors pathogenesis-related genes that were derived by lateral gene transfer from the host and other animal pathogens. The retrotransposon LINE family also appears to be acquired from animal lineages. By transcriptome analysis we show a high rate of allelic variation, which reveals rapidly evolving genes and potentially adaptive evolutionary mechanisms coupled to selective pressures exerted by the animal host. The genome and transcriptome data, as well as subsequent biochemical analyses, provided us with insight in the disease process of Saprolegnia at a molecular and cellular level, providing us with targets for sustainable control of Saprolegnia.