Transposons Passively and Actively Contribute to Evolution of the Two-Speed Genome of a Fungal Pathogen
Total Page:16
File Type:pdf, Size:1020Kb
Downloaded from genome.cshlp.org on September 24, 2021 - Published by Cold Spring Harbor Laboratory Press Research Transposons passively and actively contribute to evolution of the two-speed genome of a fungal pathogen Luigi Faino,1,3 Michael F. Seidl,1,3 Xiaoqian Shi-Kunne,1 Marc Pauper,1 Grardy C.M. van den Berg,1 Alexander H.J. Wittenberg,2 and Bart P.H.J. Thomma1 1Laboratory of Phytopathology, Wageningen University, 6708 PB Wageningen, The Netherlands; 2Keygene N.V., 6708 PW Wageningen, The Netherlands 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 ge- nome 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 path- ogen virulence. [Supplemental material is available for this article.] Genomic plasticity enables organisms to adapt to environmental identified genomic rearrangements in V. dahliae that lead to exten- changes and occupy novel niches. Although such adaptation oc- sive chromosomal length polymorphisms (de Jonge et al. 2013). curs in any organism, this is particularly relevant for pathogens Moreover, the genomes of V. dahliae strains were found to contain that engage in coevolutionary “arms races” with their hosts highly dynamic, repeat-rich, lineage-specific (LS) regions (Raffaele and Kamoun 2012; Seidl and Thomma 2014; Dong (Klosterman et al. 2011; de Jonge et al. 2013). Intriguingly, these et al. 2015). In these interactions, hosts utilize their surveillance LS regions are enriched for in planta-induced effector genes that system to detect invaders and mount appropriate defenses, involv- contribute to fungal virulence (de Jonge et al. 2013). Similarly, ing detection of invasion patterns by immune receptors, whereas many filamentous pathogens are considered to have evolved so- pathogens secrete so-called effector molecules to support host col- called “two-speed” genomes with gene-rich, repeat-poor genomic onization and counteract immune responses (Rovenich et al. compartments that contain core genes that mediate general phys- 2014; Cook et al. 2015). This tight interaction exerts strong selec- iology and evolve slowly, whereas plastic, gene-poor and repeat- tion pressure on both partners and incites rapid genomic diversifi- rich compartments are enriched in effector genes that mediate cation (McDonald and Linde 2002). virulence in interactions with host plants and evolve relatively Although sexual reproduction drives genotypic diversity, not quickly (Raffaele and Kamoun 2012). all eukaryotes regularly reproduce sexually, including many asexu- Plastic, fast-evolving genomic compartments in plant and an- al fungal phyla (McDonald and Linde 2002; Heitman et al. 2007; imal pathogen genomes concern particular regions that are either Flot et al. 2013). However, in such asexual organisms adaptive ge- embedded within the core chromosomes or reside on conditional- nome evolution also occurs, mediated by various mechanisms ly dispensable chromosomes (Thon et al. 2006; Fedorova et al. ranging from single-nucleotide polymorphisms to large-scale 2008; Haas et al. 2009; Ma et al. 2010; Goodwin et al. 2011; structural variations that can affect chromosomal shape, organiza- Klosterman et al. 2011; Rouxel et al. 2011; de Jonge et al. 2013). tion, and gene content (Seidl and Thomma 2014). Verticillium dah- Irrespective of the type of organization of the two-speed genome, liae is a soil-borne fungal pathogen that infects susceptible hosts the fast-evolving compartment is generally enriched for transpos- through their roots and colonizes the water-conducting xylem ves- able elements (TEs) that are thought to actively promote genomic sels, leading to vascular wilt disease (Fradin and Thomma 2006). changes by causing DNA breaks during excision or by acting as a Despite its presumed asexual nature, V. dahliae is a highly success- substrate for rearrangement (Seidl and Thomma 2014). Neverthe- ful pathogen that causes disease on hundreds of plant hosts less, the exact role of TEs in the evolution of effector genes current- (Fradin and Thomma 2006; Klosterman et al. 2009; Inderbitzin ly remains unknown (Dong et al. 2015). We recently established and Subbarao 2014). Using comparative genomics, we recently gapless V. dahliae whole-genome assemblies using a combination of long-read sequencing and optical mapping (Faino et al. 2015; Thomma et al. 2015). Here, we exploit these novel assemblies for 3These authors contributed equally to this work. Corresponding author: [email protected] Article published online before print. Article, supplemental material, and publi- © 2016 Faino et al. This article, published in Genome Research, is available un- cation date are at http://www.genome.org/cgi/doi/10.1101/gr.204974.116. der a Creative Commons License (Attribution 4.0 International), as described at Freely available online through the Genome Research Open Access option. http://creativecommons.org/licenses/by/4.0/. 26:1091–1100 Published by Cold Spring Harbor Laboratory Press; ISSN 1088-9051/16; www.genome.org Genome Research 1091 www.genome.org Downloaded from genome.cshlp.org on September 24, 2021 - Published by Cold Spring Harbor Laboratory Press Faino et al. in-depth investigations into the molecular mechanisms responsi- ble for genomic variability, which is instrumental for adaptive ge- nome evolution in V. dahliae. Results Genomic rearrangements in Verticillium dahliae are associated with sequence similarity Whole-genome alignments between V. dahliae strains JR2 and VdLs17 revealed 24 synteny disruptions (Fig. 1A; Supplemental Fig. S1). To further increase their resolution, we aligned long (aver- age ∼9 kb) sequencing reads derived from strain VdLs17 (Faino et al. 2015) to the strain JR2 assembly, and reads of which the two sides aligned to two distinct genomic locations were used to determine breakpoints at high-resolution (Supplemental Fig. S2). After manual refinement, this procedure yielded 19 large-scale chromosomal rearrangements; 13 inter-chromosomal transloca- tions, and six intra-chromosomal inversions (Fig. 1; Table 1). Although seven rearrangements were confined to regions smaller than 100 bp, the remaining 12 concern larger genomic regions that could not be further refined due to the presence of strain-spe- cific or repeat-rich regions (Fig. 1A). We subsequently assessed the occurrence of the 13 inter- chromosomal rearrangements in nine other V. dahliae strains by querying paired-end reads derived from these genomes for pairs of discordantly mapped reads (i.e., both reads fail to map at the ex- pected distance or location) when mapped onto the strain JR2 as- sembly, revealing distinct rearrangement patterns (Supplemental Fig. S3). Although some synteny breakpoints identified in V. dah- liae strain JR2 are either specific to V. dahliae strain VdLs17 (Chr 1: 2,843,014–2,843,020) (Supplemental Fig. S4A) or common to all other V. dahliae strains (Chr 1: 8,013,865–8,013,868) (Supplemen- tal Fig. S4B), some are observed in only a subset of V. dahliae strains (Chr 3: 3,130,776–3,130,781) (Supplemental Fig. S4C). V. dahliae strains JR2 and VdLs17 are highly similar with only 8622 SNPs difference (0.024% nucleotide divergence). When fo- cusing on the regions (10 kb) upstream of and downstream from the 19 identified synteny breakpoints, we observe ∼200 SNPs (∼0.05% divergence), indicating a moderate increase in proximity of the breakpoints. Similarly, when considering all V. dahliae strains, we observe a genome-wide average of 4.1 SNPs/kb; whereas around breakpoints (10 kb upstream and downstream), 6.3 SNPs/ kb are found. Figure 1. Extensive rearrangements in Verticillium dahliae genomes are Genomic rearrangements are generally caused by double- mediated by repetitive elements. (A) Syntenic regions, indicated by rib- strand DNA breakages followed by unfaithful repair by mecha- bons, between chromosomes of the two highly similar V. dahliae strains JR2 (chromosomes displayed in white) and VdLs17 (chromosomes dis- nisms that utilize homologous sequences to repair such breaks played in gray) reveal multiple synteny breakpoints caused by inter-chro- (Hedges and Deininger 2007; Krejci et al. 2012; Seidl and mosomal rearrangements, highlighted by red arrows for the JR2 Thomma 2014). Repetitive genomic elements, such as transpos- genome. Red bars on the chromosomes indicate lineage-specific genomic able elements (TEs), may give rise to genomic rearrangements by regions (LS) that lack synteny in the other strain. To facilitate visibility, some chromosomes of V. dahliae strain VdLs17 have been reversed and providing an ectopic substrate that interferes with faithful repair complemented (indicated