Segmental duplications drive the evolution of accessory regions in a major crop pathogen.

Many pathogens evolved compartmentalized genomes with conserved core and variable accessory regions (ARs) that carry effector genes mediating virulence. The fungal plant pathogen Fusarium oxysporum has such ARs, often spanning entire chromosomes. The presence of specific ARs influences the host range, and horizontal transfer of ARs can modify the pathogenicity of the receiving strain. However, how these ARs evolve in strains that infect the same host remains largely unknown. We defined the pan-genome of 69 diverse F. oxysporum strains that cause Fusarium wilt of banana, a significant constraint to global banana production, and analyzed the diversity and evolution of the ARs. Accessory regions in F. oxysporum strains infecting the same banana cultivar are highly diverse, and we could not identify any shared genomic regions and in planta-induced effectors. We demonstrate that segmental duplications drive the evolution of ARs. Furthermore, we show that recent segmental duplications specifically in accessory chromosomes cause the expansion of ARs in F. oxysporum. Taken together, we conclude that extensive recent duplications drive the evolution of ARs in F. oxysporum, which contribute to the evolution of virulence.

Molecular characterization reveals no functional evidence for naturally occurring cross-kingdom RNA interference in the early stages of Botrytis cinerea-tomato interaction.

Plant immune responses are triggered during the interaction with pathogens. The fungus Botrytis cinerea has previously been reported to use small RNAs (sRNAs) as effector molecules capable of interfering with the host immune response. Conversely, a host plant produces sRNAs that may interfere with the infection mechanism of an intruder. We used high-throughput sequencing to identify sRNAs produced by B. cinerea and Solanum lycopersicum (tomato) during early phases of interaction and to examine the expression of their predicted mRNA targets in the other organism. A total of 7042 B. cinerea sRNAs were predicted to target 3185 mRNAs in tomato. Of the predicted tomato target genes, 163 were indeed transcriptionally down-regulated during the early phase of infection. Several experiments were performed to study a causal relation between the production of B. cinerea sRNAs and the down-regulation of predicted target genes in tomato. We generated B. cinerea mutants in which a transposon region was deleted that is the source of c.10% of the fungal sRNAs. Furthermore, mutants were generated in which both Dicer-like genes (Bcdcl1 and Bcdcl2) were deleted and these displayed a >99% reduction of transposon-derived sRNA production. Neither of these mutants was significantly reduced in virulence on any plant species tested. Our results reveal no evidence for any detectable role of B. cinerea sRNAs in the virulence of the fungus.

Analysis of plant cell death-inducing proteins of the necrotrophic fungal pathogens Botrytis squamosa and Botrytis elliptica.

Fungal plant pathogens secrete proteins that manipulate the host in order to facilitate colonization. Necrotrophs have evolved specialized proteins that actively induce plant cell death by co-opting the programmed cell death machinery of the host. Besides the broad host range pathogen Botrytis cinerea, most other species within the genus Botrytis are restricted to a single host species or a group of closely related hosts. Here, we focused on Botrytis squamosa and B. elliptica, host specific pathogens of onion (Allium cepa) and lily (Lilium spp.), respectively. Despite their occurrence on different hosts, the two fungal species are each other's closest relatives. Therefore, we hypothesize that they share a considerable number of proteins to induce cell death on their respective hosts. In this study, we first confirmed the host-specificity of B. squamosa and B. elliptica. Then we sequenced and assembled high quality genomes. The alignment of these two genomes revealed a high level of synteny with few balanced structural chromosomal arrangements. To assess the cell death-inducing capacity of their secreted proteins, we produced culture filtrates of B. squamosa and B. elliptica that induced cell death responses upon infiltration in host leaves. Protein composition of the culture filtrate was analysed by mass spectrometry, and we identified orthologous proteins that were present in both samples. Subsequently, the expression of the corresponding genes during host infection was compared. RNAseq analysis showed that the majority of the orthogroups of the two sister species display similar expression patterns during infection of their respective host. The analysis of cell death-inducing proteins of B. squamosa and B. elliptica provides insights in the mechanisms used by these two Botrytis species to infect their respective hosts.

Dynamics in Secondary Metabolite Gene Clusters in Otherwise Highly Syntenic and Stable Genomes in the Fungal Genus Botrytis.

Fungi of the genus Botrytis infect >1,400 plant species and cause losses in many crops. Besides the broad host range pathogen Botrytis cinerea, most other species are restricted to a single host. Long-read technology was used to sequence genomes of eight Botrytis species, mostly pathogenic on Allium species, and the related onion white rot fungus, Sclerotium cepivorum. Most assemblies contained <100 contigs, with the Botrytis aclada genome assembled in 16 gapless chromosomes. The core genome and pan-genome of 16 Botrytis species were defined and the secretome, effector, and secondary metabolite repertoires analyzed. Among those genes, none is shared among all Allium pathogens and absent from non-Allium pathogens. The genome of each of the Allium pathogens contains 8-39 predicted effector genes that are unique for that single species, none stood out as potential determinant for host specificity. Chromosome configurations of common ancestors of the genus Botrytis and family Sclerotiniaceae were reconstructed. The genomes of B. cinerea and B. aclada were highly syntenic with only 19 rearrangements between them. Genomes of Allium pathogens were compared with ten other Botrytis species (nonpathogenic on Allium) and with 25 Leotiomycetes for their repertoire of secondary metabolite gene clusters. The pattern was complex, with several clusters displaying patchy distribution. Two clusters involved in the synthesis of phytotoxic metabolites are at distinct genomic locations in different Botrytis species. We provide evidence that the clusters for botcinic acid production in B. cinerea and Botrytis sinoallii were acquired by horizontal transfer from taxa within the same genus.

Dynamic virulence-related regions of the plant pathogenic fungus Verticillium dahliae display enhanced sequence conservation.

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.

In silico prediction and characterisation of secondary metabolite clusters in the plant pathogenic fungus Verticillium dahliae.

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.

The Genome of the Fungal Pathogen Verticillium dahliae Reveals Extensive Bacterial to Fungal Gene Transfer.

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 interkingdom 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 interkingdom HGT events in the genome of V. dahliae and provide evidence for extensive horizontal gene acquisition from bacterial origin.

Evolution within the fungal genus Verticillium is characterized by chromosomal rearrangement and gene loss.

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.

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.

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.

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.

Draft Genome Sequence of a Strain of Cosmopolitan Fungus Trichoderma atroviride.

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 Genome of the Saprophytic Fungus Verticillium tricorpus Reveals a Complex Effector Repertoire Resembling That of Its Pathogenic Relatives.

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.