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.

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.

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.

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.