Reference-free identification and pangenome analysis of accessory chromosomes in a major fungal plant pathogen

Reference-free identification and pangenome analysis of accessory chromosomes in a major fungal plant pathogen | bioRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-M677548'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search New Results Reference-free identification and pangenome analysis of accessory chromosomes in a major fungal plant pathogen Anouk C. van Westerhoven , Like Fokkens , Kyran Wissink , View ORCID Profile Gert Kema , View ORCID Profile Martijn Rep , Michael F. Seidl doi: https://doi.org/10.1101/2024.12.12.627383 Anouk C. van Westerhoven 1 Theoretical Biology and Bioinformatics, Utrecht University , Padualaan 8, 3583CH, Utrecht, the Netherland 2 Laboratory of Phytopathology, Wageningen University & Research , Droevendaalsesteeg 1, 6708, Wageningen, the Netherland Find this author on Google Scholar Find this author on PubMed Search for this author on this site Like Fokkens 2 Laboratory of Phytopathology, Wageningen University & Research , Droevendaalsesteeg 1, 6708, Wageningen, the Netherland Find this author on Google Scholar Find this author on PubMed Search for this author on this site Kyran Wissink 1 Theoretical Biology and Bioinformatics, Utrecht University , Padualaan 8, 3583CH, Utrecht, the Netherland Find this author on Google Scholar Find this author on PubMed Search for this author on this site Gert Kema 2 Laboratory of Phytopathology, Wageningen University & Research , Droevendaalsesteeg 1, 6708, Wageningen, the Netherland Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Gert Kema Martijn Rep 3 Molecular Plant Pathology, Swammerdam Institute of Life Sciences, University of Amsterdam Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Martijn Rep Michael F. Seidl 1 Theoretical Biology and Bioinformatics, Utrecht University , Padualaan 8, 3583CH, Utrecht, the Netherland Find this author on Google Scholar Find this author on PubMed Search for this author on this site For correspondence: m.f.seidl{at}uu.nl Abstract Full Text Info/History Metrics Supplementary material Preview PDF Abstract Accessory chromosomes, found in some but not all individuals of a species, play an important role in pathogenicity and host specificity in fungal plant pathogens. However, their variability complicates reference-based analysis, especially when chromosomes are missing from reference genomes. Pangenome variation graphs offer a reference-free alternative for studying these chromosomes. Here, we constructed a pangenome variation graph for Fusarium oxysporum , a major fungal plant pathogen with a compartmentalized genome. To study accessory chromosomes, we constructed a chromosome similarity network and identified eleven conserved core chromosomes and many highly variable accessory chromosomes. Some of these are host-specific and are likely involved in determining host range, which we corroborate by analyzing nearly 600 F. oxysporum assemblies. By a reconstruction of pangenome variation graph per homologous chromosomes, we show that these evolve due to extensive structural variation as well as the exchange of genetic material between accessory chromosomes giving rise to these mosaic accessory chromosomes. Furthermore, we show that accessory chromosomes are horizontally transferred in natural populations. We demonstrate that pangenome variation graphs are a powerful approach to elucidate the evolutionary dynamics of accessory chromosomes in F. oxysporum and provides a computational framework for similar analyses in other species that encode accessory chromosomes. Introduction The number of chromosomes of individuals in a single species is generally conserved. However, in various plants, animals, oomycetes, and fungi, a variable number of chromosomes is identified within a species. ( D’Ambrosio et al ., 2017 ). These accessory chromosomes often encode very few genes, with low transcriptional activity, and little effect on the phenotype ( D’Ambrosio et al ., 2017 ; Ma et al ., 2017 ). In contrast, accessory chromosomes in various fungal plant pathogens encode genes with important roles in host infection ( Ma et al ., 2010 ; Habig et al ., 2017 ; Torres et al ., 2020 ; van Westerhoven et al ., 2024a ; Wei et al ., 2024 ). These accessory chromosomes are often highly variable across individuals and show extensive presence / absence variation ( Fouché et al ., 2018 ; Stalder et al ., 2022 ; van Westerhoven et al ., 2024a ). The presence of pathogenicity genes on these variable accessory chromosomes separates pathogenicity genes from house-keeping genes and has been thought to facilitate rapid adaptation to changing environments as well as the host immune system ( Croll & McDonald, 2012 ; Torres et al ., 2020 ). The increasing availability of continuous, chromosome-level genome assemblies enabled the identification of accessory chromosomes in numerous fungi ( Coleman et al ., 2009 ; Haas et al ., 2009 ; Ma et al ., 2010 ; Rouxel et al ., 2011 ; Goodwin et al ., 2011 ; Dong et al ., 2015 ; Henry et al ., 2021 ; Wacker et al ., 2023 ; van Westerhoven et al ., 2024a ; Barragan et al ., 2024 ), and allows for comparative genomic approaches that help to elucidate the evolution of accessory chromosomes ( Plissonneau et al ., 2016 ; Henry et al ., 2021 ; van Westerhoven et al ., 2024a ; Barragan et al ., 2024 ). However, the analysis of accessory chromosomes is limited by reference-based analyses, where the variation of accessory chromosomes can only be determined by comparisons to a single reference genome. Because of the variability of accessory regions, the reference genome often lacks specific accessory regions that will not be analyzed. To capture the complete variation in a collection of genome sequences, including the accessory chromosomes, reference-free analysis strategies are needed. Recent advances in pangenome analysis methods can now be used to capture all genetic variation in a species ( Tettelin et al ., 2005 ; The Computational Pan-Genomics Consortium, 2018). One approach is the pangenome variation graph ( Eizenga et al ., 2020 ; Hickey et al ., 2023 ), where nodes represent sequences and edges connect adjacent nodes; each path through the graph represents a single genome. This graph can be used as a framework to call variants in a population, annotate genes, and analyze chromosome structure ( Hickey et al ., 2020 ; Guarracino et al ., 2023 ; Skiadas et al ., 2024 ). Thus far, pangenome variation graphs have been mainly constructed to analyze the human genome as well as some plant genomes ( Gao et al ., 2019 ; Wang et al ., 2022 ; Guarracino et al ., 2023 ). Only recently a pangenome variation graph for the oomycete plant pathogen Peronospora effusa has been constructed to study processes driving genome evolution ( Skiadas et al ., 2024 ). In this study the construction of the pangenome graph was guided by a reference genome ( Hickey et al ., 2020 , 2023 ), which can offer important insights in organisms where the chromosome structure is largely conserved. However, pangenome variation graphs are more challenging to construct for fungal genomes that typically contain more chromosomal variations, especially in accessory chromosomes ( Croll & McDonald, 2012 ; Fouché et al ., 2018 ; Dijkstra et al ., 2024 ), and fungal pangenome variation graphs have not been published. Due to the lack of co-linearity, the analysis of accessory chromosomes using pangenome variation graphs is not straightforward. Grouping the chromosomes into homologous groups ( Garrison et al ., 2023 ) prior to construction of a variation graph is a promising strategy to analyze chromosomes in dynamic fungal species. However, its unknown how homologous groups should be constructed for species carrying accessory chromosomes. The fungal plant pathogen Fusarium oxysporum can infect various important crops ( van Dam et al ., 2016 ; Edel-Hermann & Lecomte, 2019 ; Armer et al ., 2024 ) and is known to have a compartmentalized genome where accessory regions can span entire chromosomes or can be attached to conserved core chromosomes ( Ma et al ., 2010 ; van Dam et al ., 2017 ; Armitage et al ., 2018 ; Li et al ., 2020 ; Henry et al ., 2021 ; van Westerhoven et al ., 2024a ; Dijkstra et al ., 2024 ). The core chromosomes are largely co- linear with limited large-scale variation such as deletions, insertions, inversions, or translocations ( Ma et al ., 2010 ; Fayyaz et al ., 2023 ; van Westerhoven et al ., 2024a ). In contrast, accessory chromosomes are diverse even between strains infecting the same host ( Henry et al ., 2021 ; van Westerhoven et al ., 2024a ). Importantly, horizontal transfer of accessory chromosomes between strains can transfer pathogenicity between strains in the lab ( Ma et al ., 2010 ; van Dam et al ., 2017 ; Li et al ., 2020 ). Thus far, the genome structure and diversity of accessory chromosomes in Fusarium oxysporum has mainly been analyzed for strains infecting a specific host in comparison to the widely studied reference genome F. oxysporum f.sp . lycopersici strain 4287 (Fol 4287 ) infecting tomato, which is known to have five accessory chromosomes ( Ma et al ., 2010 ). Little is known about occurrence, the frequency, and the variation of accessory chromosomes across the Fusarium oxysporum species complex. The presence of a compartmentalized genome together with the availability of various continuous genome assemblies make F. oxysporum a suitable model to analyze the dynamics of accessory chromosomes using a pangenome variation graph. Here, we constructed a pangenome variation graph to analyze the evolution of accessory chromosomes in F. oxysporum. Some of these accessory chromosomes were specific to strains infecting the same host, suggesting that these play a role in determining host specificity. Furthermore, accessory chromosomes were composed of mixed stretches of accessory regions and, likely, evolved through recombination of accessory chromosomes. These findings offer insights into the evolution of accessory chromosomes and show how pangenome variation graphs can be applied to analyze chromosome dynamics. Results Fusarium oxysporum contains eleven core chromosomes and various accessory chromosomes To study the diversity of core and accessory chromosomes in a large collection of diverse F. oxysporum strains, we obtained 73 highly continuous whole-genome assemblies (contig N50 > 2 Mb, <100 contigs) (Table S1). Together, these strains contain 1,261 contigs. Most (1,131 out of 1,261) contigs encode telomeric repeats on at least one end, and 446 represent complete chromosomes containing telomeric repeats on both sites, we therefore refer to all contigs as chromosomes in this manuscript. The dataset contains F. oxysporum strains assigned to 19 different formae speciales based on their capacity to infect specific host plants, three endophytes, and eight strains that cause plant disease but have not been assigned to a formae speciales ( Fig. 1a,b ). The strains have an average nucleotide similarity of 97.7%, and span all three known phylogenetic clades within the F. oxysporum species complex ( Fig. 1a ; O’Donnell et al ., 1998 ; Maryani et al ., 2019 ). Download figure Open in new tab Figure 1 – Eleven core-chromosomes as well as a high number of accessory chromosomes are present in Fusarium oxysporum . a ) Phylogenetic tree of 73 Fusarium oxysporum strains, based on 4,443 single-copy BUSCO genes. Different grey boxes highlight strains belonging to the three taxonomic lineages in the F. oxysporum species complex. The outer ring shows the different hosts that can be infected by a strain, colored according to (b). b ) Distribution of plant host that can be infected by the 73 F. oxysporum strains. Endophytic strains as well as plant pathogens not assigned to a formae speciales (Other) are separated. c ) The cartoon visualizes how sequences mappings between chromosomes are translated to a chromosome network. d ) Chromosome network displays all pair-wise mappings between the chromosomes of the 73 F. oxysporum strains. Nodes represent chromosomes and edges represent mappings between chromosomes. Colors indicate the different communities that are detected. Chromosome names correspond to the chromosomes found in the Fol 4287 reference genome, and these chromosomes are depicted by a large node in the graph. Orange circles highlight accessory chromosomes. All accessory chromosomes found in reference genome Fol 4287 occur in cluster 11 (dark blue cluster). e ) Mapping of communities to homologous chromosome in reference Fol 4287, colored according to (d) f ) Communities can be associated to specific chromosomes in Fol4 287. The squares display the fraction of mappings of the Fol 4287 chromosome (y-axis) with the corresponding community (x-axis). Core chromosomes map largely to a single community, while all accessory chromosomes are assigned to community 11. g ) Network statistics of the different communities in the chromosome network. Colors correspond to the communities in the network, see (d). To study the variation of chromosome content in these different F. oxysporum strains and to facilitate pangenome analysis, we first sought to detect groups of co-linear chromosomes based on all-vs-all homology mapping between all 1,261 chromosomes ( Fig. 1c ). We tried various thresholds for mapping, i.e., the percent identity and the length, to determine the optimal mapping strategy; the mappings need to be specific enough to map homologous chromosomes and to prevent spurious matches, yet sensitive enough to capture similarity between more distantly related chromosomes. These mappings between chromosomes were then translated into a chromosome network, and we used network properties to detect communities within the network that represent groups of homologous chromosomes ( Fig. 1d ). Since core chromosomes are expected to be conserved and co-linear, we expected to find communities of core chromosomes present in all strains, as well as several accessory chromosomes with a variable presence and absence pattern. To analyze the distribution of the chromosomes over the communities, we determined the presence of chromosomes of the F. oxysporum f. sp. lycopersici reference genome Fol 4287 in the communities ( Fig. 1e-f ). The core-chromosomes of Fol 4287 are grouped into eleven separate communities, whereas the five accessory chromosomes (chromosomes 3,6,14,16, and 17) all group together into a single large community, which contained 342 chromosomes in total. Interestingly, two out of 73 strains do not encode any accessory chromosome (36102; infecting banana and Fo5; endophytic, table S2). The core- chromosomes on the other hand are all present in all 73 strains. Of the 11 core-chromosome communities, six communities contain a chromosome of all 73 strains and five communities miss one or two strains (Table S2). These missing chromosomes are the result of chromosomal rearrangements, either due to a mis-assembly or due to interchromosomal translocation. These fused chromosomes have similarity to multiple communities but are only assigned to the community with most matches. This also means that none of the strains lacks one of the core chromosomes (Table S2). Irrespective of the mapping threshold that determines the minimal length of mappings between chromosomes, we identified eleven core chromosomes, with one copy in all strains, as well as one community that contained all accessory chromosomes ( Fig. 1d ). Pangenome graph captures variation of Fusarium oxysporum and highlights differences between core and accessory chromosomes All accessory chromosomes cluster in a single community (community 11; Fig. 1d ), suggesting that many accessory chromosomes in F. oxysporum share genetic material. However, community 11 has a lower average node degree (89) and a larger diameter (5) compared with the core communities (average node degree 136; average diameter 2.2). Moreover, the nodes in this cluster are generally less well connected (transitivity 0.74 vs avg. transitivity of 0.97, figure 1g ), collectively suggesting that some of the chromosomes in the community have limited similarity to each other but are connected through similarity with other chromosomes. To further analyze the chromosome structure and the similarities between chromosomes, we constructed a pangenome variation graph using the PanGenome Graph Builder (PGGB). This graph consists of eleven subgraphs, one graph per community, and contains 500 Mb of genetic material split over 33,651,077 nodes connected through 48,127,154 edges. 2,559,224 nodes (18.6 Mb) are core (present in all 73 genomes), 13,404,553 nodes (41.1 Mb) are accessory, and 17,687,283 nodes (440.6 Mb) are unique to a single genome ( Fig. 2a ). The pangenome graph constructed for these 73 F. oxysporum genomes is open (Heaps Alpha = 0.327, Fig. 2b ), meaning that every addition of a newly sequenced genome is expected to add new genetic material to the graph. Core, accessory, and unique material is found in different quantities in the different communities. The eleven core-chromosome communities contain genetic material present in all genomes. However, in five of these core-chromosome communities more accessory material is present ( Fig. 2c ), including communities containing the core- chromosomes 12, 13 and 15 of Fol 4287. Interestingly, these chromosomes have been described as more variable than the other core chromosomes and have been previously referred to as ‘fast-core’ chromosomes ( Fokkens et al. 2018 ). However, a similar pattern is also observed in chromosome 7 and 10, which are not considered part of the fast-core chromosomes, suggesting that these chromosomes similarly contain high amount of genomic variation. The amount of unique material is consistently low in the core chromosomes, and higher in the accessory community. As may be expected, the accessory community consists solely of accessory material ( Fig. 2c ), and not a single node is present in all accessory chromosomes. Download figure Open in new tab Figure 2 – Pangenome graph of Fusarium oxysporum is open and contains extensive accessory material. a) Histogram shows the amount of genetic material (in base pairs) that is present in an increasing number of samples (x-axis). Grey bar indicates the amount of material present in only one sample (unique material), the orange bars indicate the amount of material present in one to 65 strains, dark orange indicate the amount of material found in 65-72 strains, and dark grey shows the amount of material present in all strains (core). b) Pangenome growth graph of the complete F. oxysporum pangenome, as determined by Panacus ( Parmigiani et al ., 2024 ). Heaps Alpha < 1 indicates an open pangenome. c) Proportion (in base pairs) of core, accessory, and unique genetic material found in the different communities. Each dot represents a 500 bp window in the pangenome graph. As expected, community 11, containing all accessory chromosomes, lacks core material. Stars indicate the communities with more accessory than core material. d) Size of chromosomes in accessory community 11 (orange) compared to the size of chromosomes in the other core communities (grey). e) Inversion rate per 500 bp window in the pangenome graph in the core communities (grey) and accessory community 11 (orange). The inversion rate is determined by the number of inversions per genome and per node in the pangenome graphs. To further investigate the expected co-linearity and variation between chromosomes in the pangenome graph, we analyzed the orientation and location of the nodes. We observed a higher number of inversions and nodes occur in varying order in the pangenome graph of the accessory community ( Figure 2e , Fig S1). This shows that, while core chromosomes are largely co-linear and conserved, the sequence order between homologous regions in accessory chromosomes is not conserved. To determine the amount of similarity between different strains, we determined the node-similarity in the pangenome graph. We observed that based on the core communities, most genomes are highly similar (average sharedness 80%). In the accessory community the average sharedness between genomes is lower (11%), in line with the observed variation and lack of core nodes. The similarity of the accessory regions between strains infecting the same host is higher for some formae speciales, such as the tomato infecting strains (median similarity of 46%, across 64 strains) (Fig S2), but not for all host specific strains. For example, we only observed 14% similarity across 71 banana infecting strains (Fig. S1). The high similarity in tomato infecting strains supports earlier reports showing that F. oxysporum strains infecting tomato contain a host-specific accessory chromosome that is associated with pathogenicity ( Ma et al ., 2010 ). The lack of similarity between accessory chromosomes in strains infecting other hosts suggests these might not have host-specific pathogenicity chromosomes, as we have recently proposed for banana-infecting strains ( van Westerhoven et al ., 2024a ; Dijkstra et al ., 2024 ). Accessory chromosomes can be grouped into separate co-linear communities To further delineate subgroups of accessory chromosomes that share more extensive co-linearity, we performed an additional analysis where we only considered mappings between chromosomes that cover most of the corresponding chromosome, i.e., we exclude partial mappings between chromosomes. As a result, chromosomes are now grouped into 36 communities, composed of the same eleven core communities as well as 25 additional accessory communities, demonstrating that accessory chromosomes can be further separated into groups based on co-linearity between chromosomes. The retrieval of the same eleven core chromosome communities, even when excluding splitting, further supports that core chromosomes are colinear. Nevertheless, we observed 20 cases where a chromosome maps to two different communities. To construct the final pangenome graph, these chromosomes were manually split and added into the corresponding homologous community. Not all chromosomes present in the 73 F. oxysporum strains have a homologous chromosome in another isolate, we find 107 unique chromosomes, with an average size of 1.1 Mb, without any mappings to a homologous chromosome, suggesting that these are strains specific. The separation of accessory chromosomes into smaller co-linear communities enabled us to limit our analysis to highly similar accessory chromosomes. Interestingly, we observed that some of the co-linear accessory chromosomes are unique to, or at least enriched for one of the formae speciales ( Fig. 3c ). For example, chromosome 14 of Fol 4287 is well known to play a role in pathogenicity towards tomato ( Ma et al ., 2010 ), and has been found to be present in genetically diverse tomato infecting F. oxysporum strains ( Fokkens et al ., 2018 ). Similarly, we observed that community 13, containing Fol4 287 chromosome 14, solely consists of tomato infecting F. oxysporum strains. The other five accessory chromosomes of Fol 4287 (chromosomes 2, 6, 11, 16, and 17) are found in accessory community 12. The grouping of these five accessory chromosomes supports previous reports that suggested that the accessory chromosomes in F. oxysporum strain Fol 4287 are highly similar, due to extensive segmental duplications ( Ma et al ., 2010 ; van Westerhoven et al ., 2024a ). Download figure Open in new tab Figure 3 – Accessory chromosomes can be separated into several co-linear communities, revealing host-specific accessory chromosomes. a) The chromosome network displays all pair-wise mappings spanning the full length of chromosomes. Nodes represent chromosomes and edges represent mappings between chromosomes. Colors indicate the different communities that are detected in the network. Chromosome names correspond to the chromosomes found in the Fol 4287 reference genome assembly, and these chromosomes are depicted by larger nodes in the graph. Orange circles highlight accessory chromosomes. b) Network statistics of the different communities in the chromosome network. Colors correspond to the communities in the network, see (a). c) Distribution of the different strains and their host specificity over the different communities. All strains are found in the core communities (grey line, x-axis). The accessory communities (orange line, x-axis), on the other hand, contain a subset of strains and fifteen accessory communities are host specific; the number on top indicate the total number of genomes in the community. d) Distribution of different F. oxysporum strains and their host specificity over the different communities. In total, we observed 15 clusters that are specific to a formae speciales ( Fig. 3c ), suggesting that these might also represent pathogenicity chromosomes. However, most formae speciales were only represented by a limited number of genome assemblies ( Fig. 1b,c ). To include a broader set of strains, we sought to detect presence of accessory communities in a set of genome assemblies based on short- read sequencing data. We used a collection of 588 genome sequences that had been assembled from public data available at the Sequence Read Archive (Table S3). To rapidly determine the chromosomes, present in one of these genome assemblies, we compared the k-mer profiles of these genome assemblies against the k-mer profiles of the chromosome communities. As to be expected, core communities are present in all strains and the accessory communities show a variable presence/absence pattern ( Fig. 3d ). Using the formae speciales known for 276 strains, we could distinguish 14 host specific chromosomes (present in more than four strains; Fig 3d ), including the tomato-specific community 13 that contains Fol 4287 chromosome 14. These analyses show that some accessory communities are indeed host specific in a wide variety of F. oxysporum strains, yet many other accessory communities are not restricted to a single host. Accessory chromosomes are a mosaic of accessory regions To get a more detailed insight into the variation within homologous chromosomes, we constructed a pangenome variation graph per community. This resulted in a combined pangenome, consisting of 443 Mb of genetic material that is separated over 31,289,751 nodes and connected by 44,436,449 edges. This pangenome graph consists of 18.8 Mb core, 14 Mb soft-core (>90% of strains), 161 Mb accessory, and 250 Mb of unique material, and is open (Heaps Alpha = 0.41, Fig S3). As to be expected, the fraction of core material is considerably larger than in the previous pangenome graph ( Fig. 2 ), as accessory communities now consist of more similar material (Fig. S3). We observed that the accessory chromosomes contain a high number of inversions ( Fig 2c ) and are generally not co-linear ( Fig. 3 ). We distinguished three different types of accessory communities. Most accessory communities (22) are small with only five chromosomes. The chromosomes in these clusters are overall well connected (betweenness: 30.5, diameter: 1.15) and share a large amount of genetic material (49 Mb core, 21 Mb accessory, and 27 Mb unique). In addition, we find two large accessory communities (community 11 and community 12 with 69 and 31 members, respectively). These accessory communities have a larger diameter (10.0) and a higher betweenness (45.9) than the small accessory communities, indicating that the nodes in the communities are loosely connected; not all chromosomes share genetic material, analogous to accessory community 11 in our previous analysis ( Fig. 1 ). Moreover, the pangenome graph of community 11 and 12 does not contain any core nodes (Fig. S4), further supporting the lack of shared genetic material between the chromosomes in these communities. To determine what links the chromosomes in these large and loosely connected accessory communities, we visualized the presence and absence of nodes in the pangenome graph ( Fig. 4a ). This shows that not all chromosomes share regions with each other but stretches of similar regions are found in different combinations ( Fig. 4b ). For example, different parts of an accessory chromosome ( Fol 029, chromosome 12) are present in other accessory chromosomes ( Fig. 4c ), indicating that accessory chromosomes share different accessory region. Together, these findings supports the idea that accessory chromosomes evolve through chromosomal rearrangements, resulting in a mosaic of genetic material from different sources ( Fokkens et al ., 2018 ). Download figure Open in new tab Figure 4 – Accessory chromosomes are a mosaic of accessory regions from different genomes. a) Translation from a pair-wise chromosome alignment to a pangenome graph presence/absence matrix. The visualization places the nodes in the pangenome graph on the x-axis and the genomes on the y- axis. Colored blocks indicate a node that is present in a genome, while white indicates a node that is absent in the respective genome. b) Presence/absence of nodes in the pangenome graph based on the 69 members of community 12. Note, none of the nodes are present in all chromosomes. c) Alignment of eight chromosomes found in community 12. All chromosomes are aligned to Fol 0029 chromosome 12 (pink) and share different but similar genomic regions with the other chromosomes. The connecting lines indicated regions of shared similarity, colored according to query chromosome. Grey lines indicate similarity between two neighboring chromosomes, chromosome 2 and chromosome 6 of Fol 4287 are homologous. d) Community 8 shows a collection of core chromosomes surrounded by various accessory chromosomes (highlighted in orange). e) Node presence/absence matrix of a subset of chromosomes in community 8. The node presence and absence clearly demonstrate that the core part of the chromosomes is similar across all 73 genomes (grey line). Importantly, some core chromosomes contain an accessory region, while in other cases the accessory regions are separate from the core chromosome. Next to these two types of accessory communities, we also identified accessory chromosomes present in two core chromosome communities (community 8 and community 9 with 63 accessory and 14 accessory chromosomes, respectively; Fig. 4d ). This pattern can indicate that accessory chromosomes originated from core chromosomes or that accessory chromosomes fused to core chromosomes. We observed that similar core chromosomes are present in community 8, as well as additional separate accessory chromosomes that do not share homology to the core regions ( Fig. 4d,e ). However, this accessory region is attached to the core chromosome in some strains, whilst in others the accessory region is separate ( Fig. 4e ). This demonstrates that the link between core and accessory chromosomes in this community is caused by the fusion of an accessory region to a core chromosome, an arrangement that has been previously observed in various F. oxysporum genome assemblies ( Ma et al ., 2010 ; van Westerhoven et al ., 2024a ; Dijkstra et al ., 2024 ). Accessory chromosomes can be transferred between genetically different F. oxysporum strains Accessory chromosomes can be transferred between F. oxysporum strains as has been experimentally shown for a few chromosomes in some strains ( Ma et al ., 2010 ; van Dam et al ., 2017 ; Li et al ., 2020 ) Horizontal transfer can also occur in natural populations and has been identified based on the presence of similar accessory chromosomes between distantly related strains ( van Dam et al ., 2016 ; Fokkens et al ., 2018 ; Henry et al ., 2021 ). To determine horizontal transfer events in our dataset, we analyzed the similarity of accessory chromosomes between strains from different taxonomic clades. We observed that ten communities contain members from different taxonomic clades of the F. oxysporum species complex ( Fig. 5a ), suggesting that horizontal transfer might have taken place. For instance, community 11 contains strains from the taxonomic clade 1, 2, and 3. One strain from clade 1 ( Foc 013) shares 58% similarity with clade 2 genome Fom 0021, but only 12% similarity with closest clade 1 neighbor ( Fo 4, Fig. S5), indicating that horizontal transfer has occurred between strains from clade 1 and clade 2. Download figure Open in new tab Figure 5 – Similar accessory chromosomes can be found in different Fusarium oxysporum clades. a) Distribution of F. oxysporum clades over the different chromosome communities (x-axis). Ten accessory communities contain strains from different taxonomic clades, indicating that similar chromosomes can be shared by F. oxysporum strains from with large phylogenetic distances, even from different taxonomic clades. b ) Similarity of six accessory chromosomes found in the melon infecting F. oxysporum strain Fom 0021 (x-axis) across all 73 F. oxysporum strains (y-axis), ordered based on the phylogenic relationship of the F. oxysporum strains (left). Only alignments that span at least 5kb are shown. The color indicates the percent identity. Accessory regions from Fom 0021 are found in different strains in clade 2 as well as in one strain ( Foc 013) in clade 1, highlighted by an arrow. To get a more detailed insight into the traces of horizontal chromosome transfer between Fom 0021 and Foc 013, we performed whole-genome alignments between the accessory chromosomes found in Fom 021 to all 73 genomes. Interestingly, the accessory chromosomes are present in some but not all of the most closely related species, and similar genomic regions could even be found in a strain from different taxonomic clades. For example, chromosome 12, 13, 14 and 16 from Fom 0021 (formae speciales melonis ) are also found in Foc 013 (formae speciales cucumerinum ) from clade 1 and a region of accessory chromosome 14 of Fom 0021 can be found in four banana infecting genomes (formae speciales cubense ) in clade 1. This further corroborates that accessory chromosomes are shared across genomes and further highlights that horizontal transfer occurs between clades of F. oxysporum . Interestingly, the shared genomic regions are present in strains infecting different hosts, suggesting that the transfer is not always associated with the same host range. Discussion Pangenome variation graphs offer a reference-free genome representation of different individuals of the same species. This can provide insights into the genome organization, and can improve downstream genome-based analysis such as variant calling, and gene-annotation ( Li et al ., 2022 ; Wang et al ., 2022 ; Guarracino et al ., 2023 ; Skiadas et al ., 2024 ). Such reference free analysis would be especially useful for analyzing accessory chromosomes in filamentous plant pathogens, as these are not necessarily present in the reference genome assembly. Moreover, accessory chromosomes are often highly diverse hampering reference-based analyses ( Sherman & Salzberg, 2020 ; Barragan et al ., 2024 ). By combining a homologous chromosome mapping strategy with the construction of a pangenome variation graph, we captured all accessory chromosomes in a collection of 73 F. oxysporum strains and obtained insights into their evolutionary dynamics. We found eleven conserved core chromosomes and a large number of different accessory chromosomes, which evolve through rearrangements that result in a mosaic of different accessory genomic regions. Moreover, we found evidence for horizontal transfer of accessory chromosomes between naturally occurring F. oxysporum strains from different phylogenetic clades. These results highlight that pangenome variation graphs can be successfully reconstructed for species carrying accessory chromosomes and can guide the exploration of genetic variation that underlies specific phenotypes and provides insights into their genome organization to obtain insights into genome evolution. F. oxysporum is a highly diverse species complex that infects a wide variety of hosts ( Edel-Hermann & Lecomte, 2019 ; Fayyaz et al ., 2023 ; Armer et al ., 2024 ). We identified an open pangenome, meaning that every new strain adds yet unseen additional genomic regions to the pangenome, which is in line with previous gene-based pangenomes in F. oxysporum ( Fayyaz et al ., 2023 ; van Westerhoven et al ., 2024a ) and is expected for a highly diverse species ( Tettelin et al ., 2005 ). The rapid diversification of accessory chromosomes can underly the observed diversity ( Dong et al ., 2015 ; Yang et al ., 2020 ). Using the pangenome variation graph, we identified a large set of diverse accessory chromosomes with a mosaic composition, encoding a unique combination of accessory regions, similar to what has been observed in the accessory genome of other fungi ( Fokkens et al ., 2018 ; Langner et al ., 2021 ; Barragan et al ., 2024 ). Moreover, similar accessory chromosomes are found between strains from different clades, highlighting that accessory chromosomes can be horizontally transferred between distantly related strains. This supports previous studies suggesting that clades are not genetically isolated in F. oxysporum ( van Westerhoven et al ., 2024b ), and that genetic material can be exchanged not only within but also between taxonomic clades. The horizontal transfer of accessory chromosomes, together with rearrangement of accessory material, can provide genetic diversification that is especially important in asexual fungal species ( Seidl & Thomma, 2014 ; Fokkens & Rep, 2023 ). In F. oxysporum , many accessory chromosomes carry important pathogenicity genes and play a role in determining the host range ( Ma et al ., 2010 ; van Dam et al ., 2017 ; Armitage et al ., 2018 ). Similar to previous studies, we found highly similar and host-specific accessory chromosomes within several formae speciales. However, host-specific accessory chromosomes are not identified for all formae speciales. It remains unclear why some formae speciales have a clear pathogenicity related accessory chromosomes, like tomato ( Fokkens et al ., 2018 ), whereas others have variable accessory chromosomes ( Henry et al ., 2021 ; van Westerhoven et al ., 2024a ). This variability suggests that some host- specificities arose multiple times independently, or the genes required for this host specificity are reshuffled between accessory chromosomes or between accessory chromosomes and the core. Additionally, limited phenotyping data makes it difficult to study host specificity. In general, formae speciales are assigned based on the isolation source, but a strain may infect multiple hosts or cause different symptoms in the same host ( Henry et al ., 2021 ; Batson et al ., 2021 ). The pangenome constructed in this study can serve as a valuable resource for exploring chromosome diversity within and between formae speciales of F. oxysporum . This can enable the identification of conserved regions and essential pathogenicity genes, providing insights into the host range of strains. Additionally, the pangenome can place newly sequenced and assembled F. oxysporum strains into a larger context, this can reveal the presence of accessory chromosomes, show what accessory regions are shared and thereby provide an important framework to clarify the role of accessory regions in pathogenicity. Different forms of pangenomes can be constructed, based on genes, whole-genomes, or k-mer content ( Andreace et al ., 2023 ). The optimal approach depends on the biology of a species, as well as the exact biological question to be addressed. Research on pangenomes in fungi thus far mainly focused on gene- based pangenomes, representing all genes within a species ( Badet et al ., 2020 ; Barber et al ., 2021 ). This provides important insights into the evolution of gene content, especially of pathogenicity genes, but typically offers limited information on the genome organization and the presence and shape of accessory chromosomes. Insights into accessory chromosome dynamics have been previously obtained by pairwise or reference-based comparisons ( Feurtey et al ., 2023 ; van Westerhoven et al ., 2024a ; Barragan et al ., 2024 ). The here applied reference-free pangenome variation graph enables large-scale comparisons of accessory chromosomes and can be applied in different fungi to obtain insights into chromosome dynamics. Moreover, a pangenome variation graph is valuable data by improving variant calling ( Sirén et al ., 2021 ) as well as annotation of genes and transposable elements ( Fiddes et al ., 2018 ; Skiadas et al ., 2024 ; Groza et al ., 2024 ), thereby enhancing our capabilities to study a species’ biology. Future analyses of pangenome variation graph can thus not only help to reveal insights into the evolution of the genome organization, but also can help to reveal new insights into the population structure, pathogenicity genes, and other drivers of phenotypic variations. Material and Methods DNA isolation and whole-genome sequencing and assembly For strains Fol005, Fol007, Fom011, Fom014, Fom021, Fom024 and Fom025, strains were taken from glycerol stock, grown on PDA for three days at 25°C, after which 100 ml of NO3 was inoculated with an overgrown agar plug. Samples were grown in this liquid culture for 5-7 days at 25°C and 150 RPM. Mycelium was harvested by filtering through miracloth, stored in liquid N2 and freeze-dried overnight. Per sample DNA was isolated from ∼250 mg of ground mycelium using multiple rounds of phenol- chloroform extractions, two rounds of chloroform extractions, and was precipitated twice. Quality and quantity of DNA was checked with Nanodrop, Qubit and agarose gels. DNA samples were sequenced at KeyGene on three PromethION FLO_PRO002 cells where basecalling was performed with the high- accuracy model in MinKNOW 4.2.6 and reads with quality q >= 7 were selected. Adapters were trimmed with poreChop, and reads were filtered based on their quality and length with Filtlong (version 0.2.0, --min_length 1000, --min_mean_q 80 en --min_window_q 70). Per strains, the 80% longest reads were selected and assembled with Flye (v2.8). Strains Fol029 and FolMN25 were sequenced and assembled according to van Dam et al ., 2017 . Genome selection Fusarium oxysporum genome assemblies have been downloaded from NCBI (July 2024), retaining genome assemblies with less than a 100 contigs and an N50 larger than 2 Mb. In addition to the genomes from NCBI, we included twelve previously assembled whole-genome assemblies (Table S1). Not all genomes have been assembled and processed in the same way. To account for this heterogeneity, we applied a homogenous filtering strategy and excludes fragmented chromosomes. We filtered contigs smaller than 50 kb, removing between 0 and 80 contigs per genome assembly (Table S1). We then accessed the genome assembly quality using Quast version 5.0.2 ( Gurevich et al ., 2013 ) and determined the presence of telomeric repeats using tidk version 0.2.41 ( Brown et al ., 2023 ). Furthermore, we constructed a phylogenetic tree of the strains using conserved single-copy BUSCO genes with the BUSCO phylogenomic pipeline, available at https://github.com/jamiemcg/BUSCO_phylogenomics . Network construction To detect groups of homologous chromosomes, we first mapped all contigs to each other using wfmash version 0.10.3 ( Marco-Sola et al ., 2021 ). To determine the optimal settings for F. oxysporum , we tried different parameters for the length of mapped segments (-s, 2.5kb, 5kb, and 10kb) and the block length the chained segments (-l, from 5kb, 10kb, 25kb,50kb, and 100kb); increasing segment and block length also increases the length of co-linear blocks between chromosomes. The resulting all-vs-all homology mappings between chromosomes were used to construct a chromosome network. In addition to this network, we also created a network with splitting of chromosomes disabled (wfmash flag -N), removing partial chromosome mappings. These two all-vs-all homology mappings were treated similar for downstream comparisons. Chromosomes communities were detected in the all-vs-all homology mappings using the Leiden Algorithm ( Traag et al ., 2019 ) as described by ( Guarracino et al ., 2023 ; Garrison et al ., 2024 ). The chromosome homology networks were visualized using Gephi version 0.10 ( Bastian et al ., 2009 ), and the network layout was calculated using the YiFan Hu algorithm. Community statistics were determined using iGraph python module version 0.10.8 ( Csárdi & Nepusz, 2006 ). To understand the chromosome network in relation to the F oxysporum reference genome assembly Fol 4287, we plotted the communities assigned to chromosomes in the Fol 4287 using pafR version 0.0.2 ( https://dwinter.github.io/pafr/ ). Splitting of chromosomes Before we could construct the pangenome graphs per community, we needed to ensure that each community represent homologous chromosomes. However, in some cases one chromosome maps to two different communities, for example due to chromosome fusion or errors during genome assembly. To include chromosome that map to two different communities, we split these chromosomes at the breakpoint and assigned chromosomal regions to the corresponding communities. Chromosomes were split into two different communities when a consecutive stretch of 1 Mb mapped to a different community and when this part occurred at the flank of the chromosome, and thus we only split large co-linear fractions of the chromosome, preventing us to artificially splitting chromosomes into too small sections. Some windows in the chromosome did not map to any community. To maintain the chromosome structure, these unmapped regions were ignored when we determine the flanks of a chromosome and, these unmapped regions were assigned to the new community when they flank the split region. Pangenome variation graph construction We constructed a pangenome variation graph per community detected in the all-vs-all chromosome network. The pangenome graph was constructed with the PanGenomeGraphBuilder version 0.3.0 (PGGB, Garrison et al., 2024 ) with the following parameters: a segment length of 2,500 bp (-s), a block length of 50,000 bp (-l), a percent identity threshold of 95% (-p), and 72 haplotypes (n-1). The generated pangenome variation graph was then sorted and visualized using the optimized dynamic genome/graph implementation version 0.8.6.0 (odgi, Guarracino et al., 2022 ). The general graph statistics were calculated with odgi stats (-S), the pairwise genome similarity was assessed with odgi similarity, and the percentages of core and accessory nodes were determined using odgi stats (-a). To assess structural variation, the mean inversion rate was calculated using odgi bin -w 50,000 bp, and node presence/absence patterns were visualized with odgi viz. Subgraphs for different communities were combined using odgi squeeze to form the final F. oxysporum pangenome variation graph. Finally, this joined F. oxyporum variation graph was used to evaluate the pangenome growth curve with Panancus (version 0.2.3, Parmigiani et al., 2024 ). Short read genome assemblies and matching to the pangenome graph To include a larger variety of F. oxysporum strains, we compared the k-mer profile of 588 fragmented F. oxysporum assemblies (Table S3) to the k-mer profile of the co-linear accessory communities, based on our unsplit mapping approach. The strains were assembled from short read sequencing data downloaded from NCBI using spades version 3.13.0 ( Bankevich et al ., 2012 ). To compare the k-mer profiles, we used sourmash gather (--bp-threshold 0), version 4.8.9 ( Pierce et al ., 2019 ), and determined the presence of a chromosome community in the F. oxysporum strains based on the percentage of k-mer from the chromosome community that matched to the genome assemblies. The k-mer profile of a community is based on all chromosomes in the community, including similar chromosomes, and thus the short read k-mers will always match only a small percentage of the community. We considered that the maximal percentage of k-mers from a genome assembly to match a community to represent the optimal mapping. This maximum mapping percentage is used to normalize the matching percentage across all genome assemblies, and a community is considered present in a genome assembly when at least 20% of the maximum percentage is successfully mapped. Analysis of accessory communities To visualize similarity between accessory chromosomes we used PyGenomeViz version 1.0 ( https://moshi4.github.io/pyGenomeViz/ ), in addition to the pangenome based visualizations from odgi viz ( Guarracino et al ., 2022 ). To visualize similarities to the Fom0021 reference genome, we used nucmer –maxmatch to align all assemblies to the Fom021 assembly. We filtered the alignments to only retrieve alignments with length larger than 5 kb and a precent identity higher than 90% using custom Python scripts ( https://github.com/LikeFokkens/FOSC_multi-speed-genome ). Data availability Genome sequencing data and assemblies used in this study can be found in Table S1 and S3. Custom code generated for this study can be found https://github.com/Anouk-vw/FOSC_pangenome and https://github.com/LikeFokkens/FOSC_multi-speed-genome Acknowledgements A.C.W, and G.H.J.K were supported by a grant from the Bill and Melinda Gates Foundation to the International Institute of Tropical Agriculture, Agreement No - AG-5797. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Footnotes ↵ $ Current address: Biosystems Data Analysis, Swammerdam Institute of Life Sciences, University of Amsterdam References 1. ↵ Andreace F , Lechat P , Dufresne Y , Chikhi R . 2023 . Comparing methods for constructing and representing human pangenome graphs . Genome Biology 24 : 1 – 19 . OpenUrl CrossRef PubMed 2. ↵ Armer VJ , Kroll E , Darino M , Smith D , Urban M , Hammond-Kosack KE . 2024 . Navigating the Fusarium species complex: Host-Range Plasticity and Genome Variations . Fungal Biology . 3. ↵ Armitage AD , Taylor A , Sobczyk MK , Baxter L , Greenfield BPJ , Bates HJ , Wilson F , Jackson AC , Ott S , Harrison RJ , et al. 2018 . Characterisation of pathogen-specific regions and novel effector candidates in Fusarium oxysporum f. sp. cepae . Scientific Reports 8 : 13530 . OpenUrl CrossRef PubMed 4. ↵ Badet T , Oggenfuss U , Abraham L , McDonald BA , Croll D . 2020 . A 19-isolate reference-quality global pangenome for the fungal wheat pathogen Zymoseptoria tritici . BMC Biology 18 : 1 – 18 . OpenUrl CrossRef PubMed 5. ↵ Bankevich A , Nurk S , Antipov D , Gurevich AA , Dvorkin M , Kulikov AS , Lesin VM , Nikolenko SI , Pham S , Prjibelski AD , et al. 2012 . SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing . Journal of Computational Biology 19 : 455 – 477 . OpenUrl CrossRef PubMed 6. ↵ Barber AE , Sae-Ong T , Kang K , Seelbinder B , Li J , Walther G , Panagiotou G , Kurzai O . 2021 . Aspergillus fumigatus pan-genome analysis identifies genetic variants associated with human infection . Nature Microbiology 6 : 1526 – 1536 . OpenUrl CrossRef PubMed 7. ↵ Barragan AC , Latorre SM , Malmgren A , Harant A , Win J , Sugihara Y , Burbano HA , Kamoun S , Langner T . 2024 . Multiple horizontal mini-chromosome transfers drive genome evolution of clonal blast fungus lineages . Molecular Biology and Evolution: msae 164 . 8. ↵ Bastian M , Heymann S , Jacomy M . 2009 . Gephi: An Open Source Software for Exploring and Manipulating Networks . Proceedings of the International AAAI Conference on Web and Social Media 3 : 361 – 362 . OpenUrl CrossRef 9. ↵ Batson AM , Fokkens L , Rep M , du Toit LJ. 2021 . Putative Effector Genes Distinguish Two Pathogenicity Groups of Fusarium oxysporum f. sp. spinaciae . Molecular Plant-Microbe Interactions® 34 : 141 – 156 . OpenUrl CrossRef 10. ↵ Brown M , González De la Rosa PM , Mark B. 2023 . A Telomere Identification Toolkit. 11. ↵ Coleman JJ , Rounsley SD , Rodriguez-Carres M , Kuo A , Wasmann CC , Grimwood J , Schmutz J , Taga M , White GJ , Zhou S , et al. 2009 . The genome of Nectria haematococca : contribution of supernumerary chromosomes to gene expansion . PLOS Genetics 5 : e1000618 . OpenUrl CrossRef 12. ↵ Croll D , McDonald BA . 2012 . The Accessory Genome as a Cradle for Adaptive Evolution in Pathogens . PLOS Pathogens 8 : e1002608 . OpenUrl CrossRef PubMed 13. ↵ Csárdi G , Nepusz T. 2006 . The igraph software package for complex network research. 14. ↵ van Dam P , Fokkens L , Ayukawa Y , van der Gragt M , ter Horst A , Brankovics B , Houterman PM , Arie T , Rep M. 2017 . A mobile pathogenicity chromosome in Fusarium oxysporum for infection of multiple cucurbit species. Scientific Reports 7 : 9042 . OpenUrl CrossRef PubMed 15. ↵ van Dam P , Fokkens L , Schmidt SM , Linmans JHJ , Kistler HC , Ma L-J , Rep M. 2016 . Effector profiles distinguish formae speciales of Fusarium oxysporum . Environmental Microbiology 18 : 4087 – 4102 . OpenUrl CrossRef PubMed 16. ↵ D’Ambrosio U , Alonso-Lifante MP , Barros K , Kovařík A , Mas de Xaxars G , Garcia S. 2017 . B- chrom: a database on B-chromosomes of plants, animals and fungi . New Phytologist 216 : 635 – 642 . OpenUrl CrossRef PubMed 17. ↵ Dijkstra J , van Westerhoven AC , Gómez-Gil L , Aguilera-Galvez C , Nakasato-Tagami G , Garnier SD , Yamazaki M , Arie T , Kamakura T , Arazoe T , et al. 2024 . Extensive intrachromosomal duplications in a virulence-associated fungal accessory chromosome. : 2024.09.16.611982. 18. ↵ Dong S , Raffaele S , Kamoun S . 2015 . The two-speed genomes of filamentous pathogens: waltz with plants . Current Opinion in Genetics & Development 35 : 57 – 65 . OpenUrl CrossRef PubMed 19. ↵ Edel-Hermann V , Lecomte C . 2019 . Current status of Fusarium oxysporum formae speciales and races . Phytopathology 109 : 512 – 530 . OpenUrl CrossRef PubMed 20. ↵ Eizenga JM , Novak AM , Sibbesen JA , Heumos S , Ghaffaari A , Hickey G , Chang X , Seaman JD , Rounthwaite R , Ebler J , et al. 2020 . Pangenome graphs . Annual Review of Genomics and Human Genetics 21 : 139 – 162 . OpenUrl CrossRef PubMed 21. ↵ Fayyaz A , Robinson G , Chang PL , Bekele D , Yimer S , Carrasquilla-Garcia N , Negash K , Surendrarao A , von Wettberg EJB , Kemal S-A , et al. 2023 . Hiding in plain sight: genome-wide recombination and a dynamic accessory genome drive diversity in Fusarium oxysporum f.sp. ciceris . Proceedings of the National Academy of Sciences 120 : e2220570120 . OpenUrl CrossRef PubMed 22. ↵ Feurtey A , Lorrain C , McDonald MC , Milgate A , Solomon PS , Warren R , Puccetti G , Scalliet G , Torriani SFF , Gout L , et al. 2023 . A thousand-genome panel retraces the global spread and adaptation of a major fungal crop pathogen . Nature Communications 14 : 1059 . OpenUrl CrossRef PubMed 23. ↵ Fiddes IT , Armstrong J , Diekhans M , Nachtweide S , Kronenberg ZN , Underwood JG , Gordon D , Earl D , Keane T , Eichler EE , et al. 2018 . Comparative Annotation Toolkit (CAT)—simultaneous clade and personal genome annotation . Genome Research 28 : 1029 – 1038 . OpenUrl Abstract / FREE Full Text 24. ↵ Fokkens L , Rep M . 2023 . Population genomics reveals meiotic recombination in Fusarium oxysporum . Proceedings of the National Academy of Sciences 120 : e2309677120 . OpenUrl CrossRef PubMed 25. ↵ Fokkens L , Shahi S , Connolly LR , Stam R , Schmidt SM , Smith KM , Freitag M , Rep M. 2018 . The multi-speed genome of Fusarium oxysporum reveals association of histone modifications with sequence divergence and footprints of past horizontal chromosome transfer events . : 465070 . 26. ↵ Fouché S , Plissonneau C , McDonald BA , Croll D . 2018 . Meiosis Leads to Pervasive Copy-Number Variation and Distorted Inheritance of Accessory Chromosomes of the Wheat Pathogen Zymoseptoria tritici . Genome Biology and Evolution 10 : 1416 – 1429 . OpenUrl CrossRef PubMed 27. ↵ Gao L , Gonda I , Sun H , Ma Q , Bao K , Tieman DM , Burzynski-Chang EA , Fish TL , Stromberg KA , Sacks GL , et al. 2019 . The tomato pan-genome uncovers new genes and a rare allele regulating fruit flavor . Nature Genetics 51 : 1044 – 1051 . OpenUrl CrossRef PubMed 28. ↵ Garrison E , Guarracino A , Heumos S , Villani F , Bao Z , Tattini L , Hagmann J , Vorbrugg S , Marco-Sola S , Kubica C , et al. 2023 . Building pangenome graphs. : 2023.04.05.535718. 29. ↵ Garrison E , Guarracino A , Heumos S , Villani F , Bao Z , Tattini L , Hagmann J , Vorbrugg S , Marco-Sola S , Kubica C , et al. 2024 . Building pangenome graphs . Nature Methods : 1 – 5 . 30. ↵ Goodwin SB , M’Barek SB , Dhillon B , Wittenberg AHJ , Crane CF , Hane JK , Foster AJ , der Lee TAJV , Grimwood J , Aerts A , et al. 2011 . Finished Genome of the Fungal Wheat Pathogen Mycosphaerella graminicola Reveals Dispensome Structure, Chromosome Plasticity, and Stealth Pathogenesis . PLOS Genetics 7 : e1002070 . OpenUrl CrossRef PubMed 31. ↵ Groza C , Chen X , Wheeler TJ , Bourque G , Goubert C . 2024 . A unified framework to analyze transposable element insertion polymorphisms using graph genomes . Nature Communications 15 : 8915 . OpenUrl CrossRef PubMed 32. ↵ Guarracino A , Buonaiuto S , de Lima LG , Potapova T , Rhie A , Koren S , Rubinstein B , Fischer C , Abel HJ , Antonacci-Fulton LL , et al. 2023 . Recombination between heterologous human acrocentric chromosomes . Nature 617 : 335 – 343 . OpenUrl CrossRef PubMed 33. ↵ Guarracino A , Heumos S , Nahnsen S , Prins P , Garrison E . 2022 . ODGI: understanding pangenome graphs . Bioinformatics 38 : 3319 – 3326 . OpenUrl CrossRef PubMed 34. ↵ Gurevich A , Saveliev V , Vyahhi N , Tesler G . 2013 . QUAST: quality assessment tool for genome assemblies. Bioinformatics (Oxford , England ) 29 : 1072 – 1075 . OpenUrl 35. ↵ Haas BJ , Kamoun S , Zody MC , Jiang RHY , Handsaker RE , Cano LM , Grabherr M , Kodira CD , Raffaele S , Torto-Alalibo T , et al. 2009 . Genome sequence and analysis of the Irish potato famine pathogen Phytophthora infestans . Nature 461 : 393 – 398 . OpenUrl CrossRef PubMed Web of Science 36. ↵ Habig M , Quade J , Stukenbrock EH . 2017 . Forward Genetics Approach Reveals Host Genotype- Dependent Importance of Accessory Chromosomes in the Fungal Wheat Pathogen Zymoseptoria tritici . mBio 8 : doi: 10.1128/mbio.01919-17 . OpenUrl CrossRef 37. ↵ Henry PM , Pincot DDA , Jenner BN , Borrero C , Aviles M , Nam M-H , Epstein L , Knapp SJ , Gordon TR . 2021 . Horizontal chromosome transfer and independent evolution drive diversification in Fusarium oxysporum f. sp. fragariae . New Phytologist 230 : 327 – 340 . OpenUrl CrossRef PubMed 38. ↵ Hickey G , Heller D , Monlong J , Sibbesen JA , Sirén J , Eizenga J , Dawson ET , Garrison E , Novak AM , Paten B . 2020 . Genotyping structural variants in pangenome graphs using the vg toolkit . Genome Biology 21 : 1 – 17 . OpenUrl CrossRef PubMed 39. ↵ Hickey G , Monlong J , Ebler J , Novak AM , Eizenga JM , Gao Y , Marschall T , Li H , Paten B . 2023 . Pangenome graph construction from genome alignments with Minigraph-Cactus . Nature Biotechnology : 1 – 11 . 40. ↵ Langner T , Harant A , Gomez-Luciano LB , Shrestha RK , Malmgren A , Latorre SM , Burbano HA , Win J , Kamoun S . 2021 . Genomic rearrangements generate hypervariable mini-chromosomes in host- specific isolates of the blast fungus . PLOS Genetics 17 : e1009386 . OpenUrl CrossRef PubMed 41. ↵ Li J , Fokkens L , Conneely LJ , Rep M . 2020 . Partial pathogenicity chromosomes in Fusarium oxysporum are sufficient to cause disease and can be horizontally transferred . Environmental Microbiology 22 : 4985 – 5004 . OpenUrl CrossRef 42. ↵ Li H , Wang S , Chai S , Yang Z , Zhang Q , Xin H , Xu Y , Lin S , Chen X , Yao Z , et al. 2022 . Graph- based pan-genome reveals structural and sequence variations related to agronomic traits and domestication in cucumber . Nature Communications 13 : 682 . OpenUrl CrossRef PubMed 43. ↵ Ma L-J , van der Does HC , Borkovich KA , Coleman JJ , Daboussi M-J , Di Pietro A , Dufresne M , Freitag M , Grabherr M , Henrissat B , et al. 2010 . Comparative genomics reveals mobile pathogenicity chromosomes in Fusarium . Nature 464 : 367 – 373 . OpenUrl CrossRef PubMed Web of Science 44. ↵ Ma W , Gabriel TS , Martis MM , Gursinsky T , Schubert V , Vrána J , Doležel J , Grundlach H , Altschmied L , Scholz U , et al. 2017 . Rye B chromosomes encode a functional Argonaute-like protein with in vitro slicer activities similar to its A chromosome paralog . New Phytologist 213 : 916 – 928 . OpenUrl CrossRef PubMed 45. ↵ Marco-Sola S , Moure JC , Moreto M , Espinosa A . 2021 . Fast gap-affine pairwise alignment using the wavefront algorithm . Bioinformatics 37 : 456 – 463 . OpenUrl CrossRef PubMed 46. ↵ Maryani N , Lombard L , Poerba YS , Subandiyah S , Crous PW , Kema GHJ . 2019 . Phylogeny and genetic diversity of the banana fusarium wilt pathogen Fusarium oxysporum f. sp. cubense in the Indonesian centre of origin . Studies in Mycology 92 : 155 – 194 . OpenUrl CrossRef PubMed 47. ↵ O’Donnell K , Kistler HC , Cigelnik E , Ploetz RC . 1998 . Multiple evolutionary origins of the fungus causing panama disease of banana: concordant evidence from nuclear and mitochondrial gene genealogies . Proceedings of the National Academy of Sciences 95 : 2044 – 2049 . OpenUrl Abstract / FREE Full Text 48. ↵ Parmigiani L , Garrison E , Stoye J , Marschall T , Doerr D . 2024 . Panacus: fast and exact pangenome growth and core size estimation . : 2024 . 06 .11.598418. 49. ↵ Pierce NT , Irber L , Reiter T , Brooks P , Brown CT . 2019 . Large-scale sequence comparisons with sourmash . 50. ↵ Plissonneau C , Stürchler A , Croll D . 2016 . The Evolution of Orphan Regions in Genomes of a Fungal Pathogen of Wheat . mBio 7 : e01231 – 16 . OpenUrl CrossRef PubMed 51. ↵ Rouxel T , Grandaubert J , Hane JK , Hoede C , van de Wouw AP , Couloux A , Dominguez V , Anthouard V , Bally P , Bourras S , et al. 2011 . Effector diversification within compartments of the Leptosphaeria maculans genome affected by Repeat-Induced Point mutations . Nature Communications 2 : 202 . OpenUrl CrossRef PubMed 52. ↵ Seidl MF , Thomma BPHJ . 2014 . Sex or no sex: evolutionary adaptation occurs regardless. BioEssays: News and Reviews in Molecular , Cellular and Developmental Biology 36 : 335 – 345 . OpenUrl 53. ↵ Sherman RM , Salzberg SL . 2020 . Pan-genomics in the human genome era . Nature Reviews Genetics 21 : 243 – 254 . OpenUrl PubMed 54. ↵ Sirén J , Monlong J , Chang X , Novak AM , Eizenga JM , Markello C , Sibbesen JA , Hickey G , Chang P-C , Carroll A , et al. 2021 . Pangenomics enables genotyping of known structural variants in 5202 diverse genomes . Science 374 : abg8871. 55. ↵ Skiadas P , Vidal SR , Dommisse J , Mendel MN , Elberse J , den Ackerveken GV , de Jonge R , Seidl MF . 2024 . Pangenome graph analysis reveals extensive effector copy-number variation in spinach downy mildew . PLoS Genetics : 2024 . 05 .30.596583. OpenUrl 56. ↵ Stalder L , Oggenfuss U , Mohd-Assaad N , Croll D . 2022 . The population genetics of adaptation through copy number variation in a fungal plant pathogen . Molecular Ecology . 57. ↵ Tettelin H , Masignani V , Cieslewicz MJ , Donati C , Medini D , Ward NL , Angiuoli SV , Crabtree J , Jones AL , Durkin AS , et al. 2005 . Genome analysis of multiple pathogenic isolates of Streptococcus agalactiae : Implications for the microbial “pan-genome” . Proceedings of the National Academy of Sciences 102 : 13950 – 13955 . OpenUrl Abstract / FREE Full Text 58. The Computational Pan-Genomics Consortium. 2018 . Computational pan-genomics: status, promises and challenges. Briefings in Bioinformatics 19 : 118 – 135 . OpenUrl CrossRef PubMed 59. ↵ Torres DE , Oggenfuss U , Croll D , Seidl MF . 2020 . Genome evolution in fungal plant pathogens: looking beyond the two-speed genome model . Fungal Biology Reviews 34 : 136 – 143 . OpenUrl CrossRef 60. ↵ Traag VA , Waltman L , van Eck NJ. 2019 . From Louvain to Leiden: guaranteeing well-connected communities . Scientific Reports 9 : 5233 . OpenUrl CrossRef PubMed 61. ↵ van Westerhoven , Aguilera-Galvez C , Nakasato-Tagami G , Shi-Kunne X , Martinez de la Parte E , Chavarro-Carrero E , Meijer HJG , Feurtey A , Maryani N , Ordóñez N , et al. 2024a . Segmental duplications drive the evolution of accessory regions in a major crop pathogen . New Phytologist 242 : 610 – 625 . OpenUrl CrossRef PubMed 62. ↵ Wacker T , Helmstetter N , Wilson D , Fisher MC , Studholme DJ , Farrer RA . 2023 . Two-speed genome evolution drives pathogenicity in fungal pathogens of animals . Proceedings of the National Academy of Sciences 120 : e2212633120 . OpenUrl CrossRef PubMed 63. ↵ Wang T , Antonacci-Fulton L , Howe K , Lawson HA , Lucas JK , Phillippy AM , Popejoy AB , Asri M , Carson C , Chaisson MJP , et al. 2022 . The Human Pangenome Project: a global resource to map genomic diversity . Nature 604 : 437 – 446 . OpenUrl CrossRef PubMed 64. ↵ Wei H , Zhong Z , Li Z , Zhang Y , Stukenbrock EH , Tang B , Yang N , Baroncelli R , Peng L , Liu Z , et al. 2024 . Loss of the accessory chromosome converts a pathogenic tree-root fungus into a mutualistic endophyte. Plant Communications 5 : 100672 . OpenUrl CrossRef PubMed 65. ↵ van Westerhoven AC , Dijkstra J , Palop JLA , Wissink K , Bell J , Kema GHJ , Seidl MF. 2024b . Frequent genetic exchange shapes the pan-mitogenome of the fungal plant pathogen Fusarium oxysporum . mbio : 2024 . 06 .19.599757. OpenUrl 66. ↵ Yang H , Yu H , Ma L-J . 2020 . Accessory chromosomes in Fusarium oxysporum . Phytopathology® 110 : 1488 – 1496 . OpenUrl CrossRef View the discussion thread. Back to top Previous Next Posted December 13, 2024. Download PDF Supplementary Material Email Thank you for your interest in spreading the word about bioRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following Reference-free identification and pangenome analysis of accessory chromosomes in a major fungal plant pathogen Message Subject (Your Name) has forwarded a page to you from bioRxiv Message Body (Your Name) thought you would like to see this page from the bioRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share Reference-free identification and pangenome analysis of accessory chromosomes in a major fungal plant pathogen Anouk C. van Westerhoven , Like Fokkens , Kyran Wissink , Gert Kema , Martijn Rep , Michael F. Seidl bioRxiv 2024.12.12.627383; doi: https://doi.org/10.1101/2024.12.12.627383 Share This Article: Copy Citation Tools Reference-free identification and pangenome analysis of accessory chromosomes in a major fungal plant pathogen Anouk C. van Westerhoven , Like Fokkens , Kyran Wissink , Gert Kema , Martijn Rep , Michael F. Seidl bioRxiv 2024.12.12.627383; doi: https://doi.org/10.1101/2024.12.12.627383 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Genomics Subject Areas All Articles Animal Behavior and Cognition (7651) Biochemistry (17746) Bioengineering (13928) Bioinformatics (42066) Biophysics (21499) Cancer Biology (18650) Cell Biology (25579) Clinical Trials (138) Developmental Biology (13409) Ecology (19947) Epidemiology (2067) Evolutionary Biology (24374) Genetics (15633) Genomics (22557) Immunology (17775) Microbiology (40505) Molecular Biology (17217) Neuroscience (88796) Paleontology (667) Pathology (2845) Pharmacology and Toxicology (4836) Physiology (7664) Plant Biology (15179) Scientific Communication and Education (2047) Synthetic Biology (4304) Systems Biology (9839) Zoology (2272)