{"paper_id":"1e39ddb7-d016-45f3-8d63-496ddd5bc410","body_text":"Transposable elements create distinct genomic niches for effector evolution among Magnaporthe oryzae lineages | 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 Transposable elements create distinct genomic niches for effector evolution among Magnaporthe oryzae lineages View ORCID Profile Ana Margarida Sampaio , View ORCID Profile Daniel Croll doi: https://doi.org/10.1101/2025.05.23.655716 Ana Margarida Sampaio 1 Laboratory of Evolutionary Genetics, Institute of Biology, University of Neuchâtel , 2000 Neuchâtel, Switzerland Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Ana Margarida Sampaio Daniel Croll 1 Laboratory of Evolutionary Genetics, Institute of Biology, University of Neuchâtel , 2000 Neuchâtel, Switzerland Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Daniel Croll For correspondence: daniel.croll{at}unine.ch Abstract Full Text Info/History Metrics Supplementary material Preview PDF Abstract Plant-pathogen interactions are characterized by evolutionary arms races. At the molecular level, fungal effectors can target important plant functions, while plants evolve to improve effector recognition. Rapid evolution in genes encoding effectors can be facilitated by transposable elements (TEs). In Magnaporthe oryzae , the causal agent of blast disease in several cereals and grasses, TEs play important roles in chromosomal evolution as well as the gain or loss of effector genes in host specialized lineages. However, a global understanding of TE dynamics driving effector evolution at population scale and across lineages is lacking. Here, we focus on 16 AVR effector loci assessed across a global sampling of 11 reference genomes and 447 newly generated draft genome assemblies across all major M. oryzae lineages and outgroups. We classified each effector based on evidence for duplication, deletion and translocation processes among lineages. Next, we determined AVR gain and loss dynamics across lineages allowing for a broad categorization of effector dynamics. Each AVR was integrated in a distinct genomic niche determined by the TE activity profile contributing to the diversification at the locus. We quantified TE contributions to effector niches and found that TE identity helped diversify AVR loci. We used the large genomic dataset to recapitulate the evolution of the rice blast AVR1-CO39 locus. Taken together, our work demonstrates how TE dynamics are an integral component of M. oryzae effector evolution, likely facilitating escape from host recognition. In-depth tracking of effector loci is a valuable tool to predict the durability of host resistance. Introduction Plant-pathogen interactions are characterized by continuous evolutionary arms races, where hosts adapt to resist infection while pathogens evolve to overcome host defenses ( Sironi et al., 2015 ). Effector-receptor dynamics are crucial for this process, where pathogens deploy effectors ( i.e. avirulence factors, AVR ) to manipulate host immunity ( Lo Presti et al., 2015 ; Remick et al., 2023 ), while host receptors are responsible to detect them and trigger defense responses ( Kanyuka & Rudd, 2019 ; Stukenbrock & McDonald, 2009 ). Genes encoding such effectors are among the most rapidly evolving genes in pathogen genomes, driven by high mutation rates and strong selection pressure ( Sánchez-Vallet et al., 2018 ), enabling pathogens to evade host resistance and posing a challenge for effective disease control. Effector genes are often encoded in subtelomeric chromosomal regions, which are repeat-rich regions and evolve more rapidly than regions encoding housekeeping genes ( Dong et al., 2015 ; Faino et al., 2016 ). These genomic regions are also frequently associated with transposable element (TE) activity, which can disrupt coding sequences or promote regulatory changes ( Chuma et al., 2011 ; Fouché et al., 2018 ; Kang et al., 2001 ; Sampaio et al., 2025 ; Seidl & Thomma, 2017 ; Whisson et al., 2012 ). Moreover, these chromosomal regions are also exhibiting a higher propensity for stochastic epigenetic regulation ( Aparicio et al., 1991 ; De Las Peñas et al., 2003 ; Fan et al., 2008 ), and chromosomal rearrangements resulting in rapid effector gene evolution ( De Jonge et al., 2013 ; Shi-Kunne et al., 2018 ; Wang et al., 2020 ). Plant pathogens capable of infecting both wild host plants and cultivated crops are of particular concern given their propensity to switch to new hosts ( Anderson et al., 2004 ). Magnaporthe oryzae can infect over 50 wild and cultivated grass species including major cereal crops such as rice ( Oryza sativa ) and wheat ( Triticum aestivum ) ( Wilson & Talbot, 2009 ). Despite the wide host range, M. oryzae genotypes are grouping into host-specialized forms and recognized as different pathotypes ( Choi et al., 2013 ; Kim et al., 2019 ; Valent, 1990 ). The most studied pathotypes include M. oryzae Oryza (MoO), M. oryzae Triticum (MoT), and M. oryzae Lolium (MoL), causing blast disease on rice, wheat and ryegrass, respectively. The origin of new pathotypes was triggered by host jumps including the emergence of wheat and rice blast disease ( Couch et al., 2005 ; Inoue et al., 2017 ). Wheat blast emerged from a host jump from ryegrass ( Lolium spp.) to Brazilian wheat cultivars lacking the RWT3 gene, which is the resistance gene expressing the receptor recognizing the PWT3 effector. While both PWT3 and PWT4 effector genes are found in the Lolium pathotype, cultivation of rwt3 wheat cultivars (lacking the ability to recognize PWT3) allowed the emergence of a new Triticum pathotype isolates carrying the PWT3 effector but losing PWT4, an effector that would be recognized due to the presence of its complementary resistance gene ( RWT4 ) in wheat cultivars. Subsequent loss-of-function mutations arose due to the nearby cultivation of wheat cultivars carrying RWT3 , along with the spread of pathogens to common wheat varieties ( Inoue et al., 2017 ). These findings support the hypothesis that host specialization was mainly driven by genetic changes at effector gene loci including gain or loss of effector functions ( Yoshida et al., 2016 ). Numerous M. oryzae AVRs were cloned and characterized for their interaction with plant resistance factors ( De Wit et al., 2009 ). Among these, AVR-Pita interacts with the cognate Pi-ta resistance protein and shows homology to fungal zinc-dependent metalloproteases ( Jia et al., 2000 ; Orbach et al., 2000 ). AVR-ACE1 is involved in the production of secondary metabolites and activates the rice resistance factor Pi33 ( Böhnert et al., 2004 ; Collemare et al., 2008 ). AVR-Piz-t is recognized by Piz-t and suppresses pathogen-associated molecular pattern (PAMP)-triggered immunity ( Li et al., 2009 ; Park et al., 2012 ). PWL effectors are primarily encoded by MoO and are rapidly evolving, small, glycine-rich secreted proteins ( Kang et al., 1995 ; Sweigard et al., 1995 ). Beyond functional differences in their encoded proteins, AVR genes exhibit a high degree of genetic instability and are often localized in telomeric regions ( Farman, 2007 ; Rehmeyer et al., 2006 ). Such high rates of sequence changes are likely increasing their adaptative potential to evade host recognition. Sequence diversification occurred mostly through simple point mutations as reported for AVR-Pita and AVR-Pik ( Dai et al., 2010 ; Kanzaki et al., 2012 ; Longya et al., 2019 ), or sequence rearrangements causing segmental deletion of coding sequences as for AVR-Pita and AVR-Pib ( Orbach et al., 2000 ; Zhang et al., 2015 ) associated with gains in virulence. There is also evidence for horizontal transfer of PWT4 from M. pennisetigena to an Avena isolate from Brazil ( Inoue et al., 2021 ). There is strong evidence that TEs impacted AVR loci and facilitate rearrangements. TEs facilitated virulence gains of MoO through loss-of-function mutations such as the insertion of a Mg-SINE into the AvrPi9 coding sequence ( Wu et al., 2015 ) or gains of virulence linked to AVR-Pita and AVR-Pib due to a Pot3 TE insertion ( Hu et al., 2022 ; Kang et al., 2001 ). TEs were also likely facilitating the translocations observed for several AVR genes including AVR-Pita ( Chuma et al., 2011 ). The loss of telomeric ends resulting in the elimination of AVR-Pita and AVR-Pii ( Chuma et al., 2011 ; Khang et al., 2008 ) were likely also favored by the repetitive nature of subtelomeric regions. Overall, the rapid evolution of M. oryzae AVRs to evade host recognition is likely facilitated by TE dynamics. TEs have specifically expanded in MoO compared to MoT and MoL ( Lin et al., 2024 ; Nakamoto et al., 2023 ). Furthermore, TE insertions also mediated the divergence of M. oryzae populations infecting different rice subspecies ( Lin et al., 2024 ; Nakamoto et al., 2023 ). However, a population genomics perspective on TE impacts on effectors leveraging the vast available genomic datasets on M. oryzae lineages is lacking. Here, we used extensive genomic datasets covering all major M. oryzae host-associated pathotypes to recapitulate insertion dynamics near 16 AVRs and included M. grisea and M. pennisetigena as outgroups. We recapitulated chromosomal rearrangements affecting AVR loci in the M. oryzae pathotypes MoO, MoT, and MoL using reference-quality genomes. We combined synteny analyses with draft assemblies from available sequencing datasets to assess AVR gain/loss patterns in conjunction with TE insertion dynamics surrounding the AVRs . Material and Methods Magnaporthe genomic datasets and genome assemblies We performed analyses on a global collection of genomic datasets comprising 458 Magnaporthe isolates (Supplementary Table S1 and S2). All genomes were accessed from public databases reported by previous studies ( Gladieux et al., 2018 ; Latorre et al., 2023 ; Thierry et al., 2022 ; Zhong et al., 2018 ). The geographic origin of isolates and the host organism information was retrieved from metadata attributes and cross-checked with the associated literature. The isolates were from Asia (n=256), South America (n=96), Africa (n=69), Europe (n=11), North America (n=8), and a small number without reported origin (n=18). Out of the 458 isolates, 448 belonged to M. oryzae collected from 18 different host species, including cereals and grasses. We also included as outgroups two isolates of M. pennisetigena (collected on Pennisetum sp.) and eight isolates of M. grisea isolated from Cenchrus and Digitaria sp. Overall, 447 isolates were sequenced using Illumina paired-end whole-genome sequencing (WGS) (Supplementary Table S1). Illumina sequencing data was initially filtered using fastp v0.23.4 with default settings ( Chen et al., 2018 ) to remove adapter sequence and low-quality reads. De novo draft assemblies were generated using the software SPAdes v3.15.5 ( Bankevich et al., 2012 ) with the “careful” method and automated k-mer selection. All genomes were verified to have more than 95% completeness using BUSCO version 5.8.2 ( Simão et al., 2015 ) searching the Ascomycota orthology database. We used QUAST to calculate assembly metrics ( Gurevich et al., 2013 ). We retained assemblies with N50 (length of the shortest contig for which longer length contigs cover at least 50% of the assembly) above 17,034 bp and a total assembly size above 36.72 Mbp. Eleven Magnaporthe reference-quality genome assemblies were also included in this study ( Dean et al., 2005 ; Gomez-Luciano et al., 2019 ; Liu et al., 2024 ; Peng et al., 2019 ). These genomes included seven M. oryzae genomes (3 MoO, 3 MoT, 1 MoL) and four outgroup genomes (2 M. grisea and 2 M. pennisetigena ) (Supplementary Table S2). Phylogenetic analyses Phylogenetic relationships were assessed separately for the set of Magnaporthe reference-quality genomes and for 405 Magnaporthe draft genomes associated with the most frequently sampled host species. For the phylogenomic tree, we used single-copy genes predicted by AUGUSTUS v3.5.0 ( Stanke et al., 2004 ) using the pretrained gene prediction database available for the M. grisea genome. Predicted protein sequences were used for orthology analyses performed with Orthofinder v2.5.5 ( Emms & Kelly, 2019 ). Single-copy orthologs present in at least 90% of the total number of isolates were kept. For the phylogenetic reconstruction of the large Magnaporthe worldwide collection, we retained 300 randomly selected single-copy orthologs to reduce computational load. Selected ortholog protein sequences were aligned as a supermatrix using the AMAS tool v1.0 ( Borowiec, 2016 ). Phylogenetic trees were built using RAxML v8 ( Stamatakis, 2014 ) to construct a maximum-likelihood phylogenetic tree with the parameters -m PROTGAMMAAUTO for protein sequences with 1000 bootstrap replicates. Identification of effector homologues We analyzed the 447 draft genomes produced by SPAdes and the 11 reference-quality genomes (including outgroups) to search for homologs of 16 M. oryzae effectors. We focused on cloned and well-characterized effectors identified in Oryza or Triticum pathotype isolates ( Table 1 ). For effectors present in more than one pathotype, effector sequences with best hits in the reference-quality genomes were used as query for BLASTn analyses ( Camacho et al., 2009 ). Hits were filtered for a maximum e -value of 10 -5 , followed by individual minimum length filtering for each effector based on visual inspection of alignment length distributions. The ORF finder tool ( https://www.ncbi.nlm.nih.gov/orffinder/ ) was used to refine effector open reading frames detected by BLASTn hits. View this table: View inline View popup Download powerpoint Table 1: Overview of M. oryzae AVR effectors analyzed in this study. The AVR pathotype identifies the M. oryzae pathotype in which the effector was first characterized. Transposable element annotation and analysis To detect TE insertions near effectors, all draft assemblies, reference-quality and outgroup genomes were annotated with the TE consensus sequences reported for Magnaporthe ( Lin et al., 2024 ) (available from https://github.com/S-t-ing/mBio-data-availablility/blob/main/Mo.TE_Consensus.fasta ). For this, we used RepeatMasker v4.0.9_p2 with “-no_is” and “-nolow” parameters to skip simple repeat and low complexity region annotations. Annotated TEs shorter than 50 bp were filtered out. The search for TEs near effector loci was restricted to a 1000 bp window up- and downstream the effector gene. Only MoO and MoT isolates clearly identified in large clades based on the phylogenetic tree were retained these analyses. Percentage of TE superfamilies per AVR in MoO and MoT was calculated per effector based on the total TE number per pathotype. Percentage of TE superfamilies per insertion site (bp) has been calculated per effector based on the total number of isolates containing at least one TE for this AVR . We used the ANOVA function in R to test for the impact of TE family, AVR identity, and the AVR pathotype (MoO and MoT) on variation in TE counts near effectors. The proportion of variance explained by each factor and interactions was calculated based on the sum of squares. Analyses of the AVR1-CO39 locus To assess the AVR1-CO39 locus organization in genomes of different pathotypes, we first analyzed reference-quality genomes. Synteny was plotted with genoplotR v0.8.11 ( Guy et al., 2010 ) for the Magnaporthe reference genomes using gene and TE annotation ( Lin et al., 2024 ). For M. oryzae isolates, we mapped coding sequences from the MoO 70-15 (GCF_000002495.2) gene annotation ( Dean et al., 2005 ), and for outgroup M. pennisetigena we mapped coding sequences from the Br36 (GCA_004337985.1) gene annotation ( Gomez-Luciano et al., 2019 ). The two annotations of 70-15 and r36 provided complementary coverage of major haplotypes found at the locus. AUGUSTUS gene annotations as described above were used to infer coding sequences in non- MoO isolates. Magnaporthe genomes assembly having AVR1-CO39 were inspected for TE presence surrounding the effector gene. Results Samples distribution, genome assembling and reference genomes To unravel M. oryzae AVR evolution, we assembled a collection of genomic datasets for 458 Magnaporthe spp. isolates collected across continents and diverse hosts ( Figure 1A ) ( Gladieux et al., 2018 ; Latorre et al., 2023 ; Thierry et al., 2022 ; Zhong et al., 2018 ). Most isolates were collected in Asia (59%) and predominantly infecting O. sativa (65%) ( Figure 1A ), reflecting the high incidence of rice blast disease in this region. Isolates infecting cereals such as Triticum sp. as well as grasses such as Lolium sp. and Eleusine sp. across continents were also included ( Figure 1A ). From the 458 analyzed isolates, 448 belonged to M. oryzae (Supplementary Table S1). Genomic data included 441 short read sequencing datasets and 7 reference-quality genomes from MoO, MoL and MoT (Supplementary Table S2) ( Figure 1B ). Ten outgroup genomes from two distinct Magnaporthe species ( M. grisea and M. pennisetigena ) were also included in the analysis of which four genomes were of reference quality (Supplementary Table S2). The selected outgroup Magnaporthe species M. grisea and M. pennisetigena clustered separately from each other and separated well from M. oryzae in the phylogenomic tree analysis ( Figure 1B ). To complement the available reference genomes, we either accessed or assembled draft genomes for 447 additional isolates. Assembly genome sizes ranged from 39.1 Mbp in isolates infecting Leersia spp. and 42.9 Mbp in isolates infecting Cenchrus sp. ( Figure 1C ). Assembled genomes for this study showed acceptable contiguity for the purpose of analyzing coding regions with the N50 (length of the shortest contig for which longer length contigs cover at least 50% of the assembly) averaging between 19,433 bp in isolates infecting Pennisetum sp. and 103,607 bp in isolates infecting Setaria sp. ( Figure 1C ). Assembly genome sizes were not meaningfully correlated with the assembly contiguity ( i.e. N50), which supports the notion that the totality of the genome is reasonably well covered despite variation in assembly quality (Supplementary Figure S1). The reference-quality genomes showed BUSCO completeness scores ranging from 96.5% in PM1 ( M. pennisetigena ) to 97.9% in M. oryzae isolates infecting Oryza spp. ( Figure 1D ). All draft genome assemblies similarly exhibited >95% BUSCO completeness and are hence comparable to the reference genomes in terms of gene content. A small number of Oryza -infecting isolates showed slightly lower BUSCO completeness scores compared to the reference genomes, yet the assembly genome sizes were similar. Download figure Open in new tab Figure 1: Global panel of analyzed Magnaporthe spp. genomes. A) Geographical distribution of the Magnaporthe spp. isolates. The color identifies the reported host genus, and the size defines the number of samples from the same location. B) Phylogenomic tree of reference-quality genomes. The host genus is reported in parentheses. Bootstrap confidence values >80% are displayed. C) De novo draft assembly genome size for genomes assembled from short-read sequencing data. Genomes are grouped by reported host genus. D) Analyses of genome completeness based on BUSCO completeness percentages for reference-quality genomes (orange on top) and de novo draft assembled genomes. Rearrangements of AVRs among M. oryzae reference genomes Some M. oryzae AVR underwent chromosomal translocations. To comprehensively track AVR localization among isolates, we focused on 16 cloned and characterized AVR effectors ( Table 1 ) in seven reference M. oryzae genomes from MoO, MoL and MoT. We expected MoO effectors to be shared among most MoO reference genomes. However, from the 12 MoO AVR effectors ( Table 1 ), only seven were shared among all MoO reference genomes and only four were found in all M. oryzae reference genomes ( Figure 2A ). These four AVRs include AVR-Pi9 and AVR-Pi54 being at conserved chromosomal locations within pathotypes, while ACE1 and AVR-Pik showed translocations among MoO isolates ( Figure 2C ). AVR-Pita1 , PWL1 and PWL2 were present in all MoO reference genomes but have undergone duplication events ( Figure 2D ). We also identified deletions of two MoO AVRs ( AVR-Piz-t and AVR-Pib ) in at least one MoO reference genome ( Figure 2E ). AVR-Pii was absent in two MoO reference genomes ( Figure 2E ) but exhibited a duplication in two MoT reference genomes ( Figure 2D ). AVR-Pii was the only MoO effector that does not exhibit a duplication in a reference genome for its associated pathotype but in other M. oryzae pathotype reference genomes. AVR1-CO39 and AVR-Pia suffered the most dramatic loss, being absent in all MoO reference genomes ( Figure 2E ). AVR1-CO39 is shared among all MoL and MoT reference genomes though ( Figure 2E ). MoT effectors showed similar loss patterns to MoO effectors. Two MoT effectors were shared among the MoT reference genomes ( Figure 2B ). PWT3 is shared among all M. oryzae reference genomes and located at a conserved position on chromosome 5 ( Figure 2B ). PWT4 is absent in most M. oryzae reference genomes except MoT. PWT6 exhibited the most pronounced pattern of effector loss with the AVR being retained in only two MoO reference genomes ( Figure 2D ). Overall, AVR localizations are highly dynamic among the seven reference genomes with contributions by translocations, deletions and duplication events. Download figure Open in new tab Figure 2: Overview of AVR effector synteny among M. oryzae reference genomes from MoO, MoL and MoT. Chromosomal localizations of AVR s are reported for (A) MoO and (B) MoT. Overview of AVRs showing evidence for (C) translocation among chromosomes, as well as (D) duplication and (E) deletion events. Dashed lines connect AVR localizations between reference genomes. If an effector was missing from a genome, a dashed line without a circle symbolizes absence of the effector gene. Effector gain and loss dynamics across global M. oryzae lineages To comprehensively map the evolution trajectory of AVRs across M. oryzae lineages, we searched for effector homologs across a dataset of 458 M. oryzae genomes (Supplementary Table S3). Isolates infecting Cenchrus and Pennisetum sp. were used as outgroups to clarify gains and loss patterns. Relationships among isolates were assessed based on a phylogenomic tree including 300 single-copy orthologs. The phylogenetic grouping is consistent with previous studies of M. oryzae lineage diversification ( Figure 3 ). Isolates belonging to the Avena , Lolium and Triticum pathotypes clustered together. Among those, isolates collected from Triticum sp. were the most dispersed across the tree, corroborating the high genetic diversity reported for MoT pathotype isolates ( Rahnama et al., 2023 ). We verified consistency of phylogenetic placements for 24 isolates overlapping with a previous phylogenomic analysis ( Gladieux et al., 2018 ). Three isolates reported as collected on Oryza sp. leaves (A-PHL-64, ARG-60 and ARG-61) ( Onaga et al., 2020 ) did not cluster with the remaining MoO isolates. Furthermore, four MoT isolates were not clustering as closely with MoL and MoA as expected but rather with MoS ( Figure 3 ). Download figure Open in new tab Figure 3: Effectors gain and loss dynamics across Magnaporthe spp. Maximum likelihood phylogenomic tree of 405 Magnaporthe genomes based on 300 protein sequence alignments of single-copy genes. Bootstrap confidence values >80% are displayed. Numbers in parentheses indicate the number of genomes. Pie charts indicate the proportion of isolates carrying specific effectors per clade. The presence/absence category refers to a summary classification of AVR gain/loss patterns observed across the phylogeny: 1) stable AVR presence; 2) AVR frequency increase in pathotype originally described for the AVR function; 3) general decrease in AVR frequency across M. oryzae pathotypes; 4) AVR frequency decrease in pathotype originally described for the AVR function. To categorize effector gain and loss dynamics, we assessed the frequency of effectors according to their reported host ( Table 1 ). Most AVRs were present in at least one of the outgroups infecting Cenchrus and Pennisetum sp., respectively, suggesting that most M. oryzae effectors were present in the common ancestor to all extant M. oryzae lineages. Exceptions include AVR-Pi9 , AVR-Piz-t , AVR1-CO39 and AVR-Rmg8 being absent in the outgroups. We categorized AVRs into four types according to their frequency across M. oryzae lineages: 1) AVRs with largely stable frequencies among lineages; 2) AVRs increasing in frequency in the pathotype reportedly linked to the AVR function; 3) AVR largely lost in most M. oryzae pathotypes; 4) AVR largely lost in the pathotype reportedly linked to the AVR function. Among MoO isolates, all four gain/loss categories were observed for MoO effectors. AVR-Pi9 , AVR-Pi54 , ACE1 and AVR-Piz-t are well conserved across the M. oryzae phylogeny (category 1; Figure 3 ). The near fixation of the AVRs suggests that recognition by the host is not widely distributed among host varieties or that the AVR serves an additional function. Transposable element colonization near M. oryzae effectors The high frequency of rearrangements at AVR loci is consistent with high repetitive DNA content nearby. TE insertions are known to facilitate the creation of structural variation. Hence, we investigated patterns of TE dynamics surrounding AVRs and potential links to effector presence/absence variation. We annotated assembled contigs encoding the different AVRs for the presence of TE sequences considering a window of +/- 1000 bp of the effector coding sequence (Supplementary Table S4). The most frequent and conserved effectors ( AVR-Pi9 , AVR-Pi54 , ACE1 and AVR-Piz-t ) loci were devoid or nearly devoid of TEs in proximity in both MoO and MoT isolates ( Figure 4A and B ). This suggests that the conservation of AVR effectors is facilitated by suppressed TE activity nearby. On the contrary, AVR1-CO39 , an effector lost nearly entirely in MoO but remaining at high frequency in MoT isolates, showed the highest percentage of TEs among MoT isolates ( Figure 4C ). Hence, the high rate of TE insertions nearby could have facilitated the loss of the effector. Interestingly, the frequency of TE families is associated with the AVR identity. The two DNA transposons POT2 and POT3 are at high frequency near PWT2 , AVR-Pib and AVR-Pik loci ( Figure 4C ). Retrotransposons were also abundant near specific AVR s. LTR-RETRO5 and LTR-RETRO7 were commonly found near the MoO effectors AVR-Pita1 , AVR-Pii , and AVR1-CO39 , while LTR-RETRO6 was more frequent near the MoT effector PWT6 . The LTR-MGL3 retrotransposon showed a high frequency near AVR1-CO39 and PWT4 . Notably, LTR-Pyret was exclusively detected near PWT3 in both MoO and MoT isolates. Additionally, the MGL-LTR retrotransposon was most frequent near effectors first characterized in MoT ( Figure 4C ). The distance between TEs and the nearby effector gene varied. Effectors in MoO had typically a higher distance to TEs than effectors in MoT isolates suggesting that sequence rearrangements could have affected the spacing between TEs and effectors ( Figure 4A and B ). Download figure Open in new tab Figure 4: Transposable element (TE) dynamics and diversification near M. oryzae AVR . A) Presence of TEs in a +/- 1000 bp window surrounding AVR loci in Oryza pathotype (MoO) isolates. The circle color reflects TE identity and classification. The circle size indicates the percentage of isolates carrying the respective AVR having a specific TE present. 100% refers to all isolates with the AVR sharing a specific TE at a specific position. Counts ( n ) indicate the number of isolates exhibiting at least one TE for a given AVR . Negative bp distance values represent AVR upstream regions. B) Triticum pathotype (MoT) isolates. C) Relative abundance of TEs among AVR loci shown separately for MoO and MoT isolates. D) Analyses of variance (ANOVA) of factors potentially explaining variation in TE abundance among AVRs . The percentage of variance explained by the different models based on factorial combinations for TE, AVR and pathotype (for which the AVR was first described) and their interactions. Given the heterogeneity in TE occupancy near AVRs across pathotypes, we sought to formally assess what factors explain best the variability in TE content. Using a multi-factorial ANOVA, we found that the TE identity ( i.e. classification), the identity of the AVR and the pathotype explained a significant portion of the variance in TE occupancy ( Figure 4D , Supplementary Table S5). The largest proportion of variance (51%) in TE occupancy was explained by the interaction of AVR and TE identity. This suggests that AVR loci may have co-evolved with distinct TE families, reinforcing the idea that certain TEs have affinity for specific chromosomal regions. AVR1-CO39 locus dynamics AVR1-CO39 was the only effector at low frequency in isolates of the pathotype (MoO) in which the effector function was originally described and at high frequency in all the other pathotypes ( Figure 3 ). AVR1-CO39 was previously characterized for a sequence rearrangement at the origin of the host switch to rice ( Farman et al., 2002 ; Peyyala & Farman, 2006 ; Tosa et al., 2005 ). The two variants were described as the G- and J-type in MoO lacking a functional AVR and an alternative W-type associated with an intact AVR1-CO39 in isolates infecting weeping love grass ( Farman et al., 2002 ). The most frequent type in MoO (G-type) is characterized by a complete loss of the coding sequence through a deletion and TE replacement. The J-type consists of a loss-of-function version caused by a repetitive element called REP1 ( Farman et al., 2002 ). REP1 corresponds to LTR-RETRO6 in more recent TE annotations ( Figure 5A ). We analyzed the evolution of the AVR1-CO39 locus organization across pathotypes by inspecting first the reference genomes. We used gene models annotated in the MoO genome 70-15 to obtain gene annotations in the additional M. oryzae reference genomes included in the comparison ( Figure 5B ). The W-type was represented by the MoL and MoT reference genomes ( Figure 5B ). The W-type carries genes adjacent to AVR1-CO39 , which are absent in 70-15 (empty chromosomal region of the W-type; Figure 5B ). We identified no locus synteny with the outgroup reference genome (Br36) ( Figure 5B ). As previously reported for the MoO isolate Guy11 ( Farman et al., 2002 ), the 70-15 MoO reference genome shows a substantial contraction of the region adjacent to AVR1-CO39 compared to the MoL and MoT genomes ( Figure 5B ). Download figure Open in new tab Figure 5: Chromosomal synteny analyses of the AVR1-CO39 locus. A) Sketch of AVR1-CO39 locus types as reported in the literature ( Farman et al., 2002 ). B) Synteny plot of the AVR1-CO39 region on chromosome 2 of MoO, MoL and MoT reference genomes and M. pennisetigena . Grey blocks between chromosomes indicate homologous regions. C) Comparison of TE types upstream of AVR1-CO39 among Magnaporthe lineages. The phylogenomic tree was based on 300 single-copy orthologs. Bootstrap confidence values >80% are displayed. The contracted region in 70-15 exhibits different TEs compared to the W - type, with LTR-RETRO5 and LTR-RETRO6 occupying the contracted region. The LTR-RETRO5 has been previously identified in the G-type and the adjacent element was called REP1 , which we identified here as LTR-RETRO6 . The MoT isolate B71 exhibits yet distinct pattern with the AVR1-CO39 upstream region fully occupied by TEs ( Figure 5B ). Two of these TEs were intact LTR-RETRO7 , suggesting a recent insertion. One of the LTR-RETRO7 insertion resulted in the truncation of AVR1-CO39 in B71, differing in size from the J-type rearrangement. However, both retained ORF3, which encodes the AVR1-CO39 domain interacting with the corresponding resistance gene ( Ribot et al., 2013 ). The Zambia wheat blast outbreak isolate ZM1-2 (MoT) was reported to share a similar AVR1-CO39 locus structure as B71 ( Latorre et al., 2023 ). This highlights that the AVR1-CO39 loss-of-function is not restricted to MoO isolates but is also occurring independently in MoT isolates through a different TE insertion. Since the different AVR1-CO39 types are characterized by distinct TE content upstream of the coding sequence, we expanded our investigations to all isolates of the collection carrying the effector. We find that the LTR-MGL3 , characteristic of the W-type rearrangement, was the most frequent among Magnaporthe lineages ( Figure 5C ). In contrast, the J-type was rare and found in only 10 MoO isolates ( Figure 5C ). The LTR-RETRO7 characteristic of the J-type was found in the reference genome B71 and a small number of additional MoT isolates ( Figure 5C ). We identified two further TEs, LTR-Pyret and POT2 , which were also found directly upstream of AVR1-CO39 in our dataset. Our investigations of the locus indicate that the W-type is likely ancestral and has undergone multiple rearrangements followed by different TE insertions leading ultimately to the loss of the effector ( i.e. G-type). Discussion AVR effectors of M. oryzae have undergone rapid evolution to circumvent the matching plant receptors and increase pathogen virulence ( Huang et al., 2014 ). Point mutations, insertions and deletions are widely described mechanisms responsible for the loss of avirulence function of AVR genes ( Chuma et al., 2011 ; Hu et al., 2022 ; Orbach et al., 2000 ; Zhang et al., 2015 ). Abundance of TE sequences near effector regions can act as a mutagen or epigenetic regulator for accelerated effector evolution ( Abraham et al., 2024 ; Bao et al., 2017 ; Jeon et al., 2015 ). In M. oryzae , TE content correlates with host identity and contributes to genetic differentiation among lineages ( Lin et al., 2024 ; Nakamoto et al., 2023 ). Here, we show that effectors evolved well differentiated presence/absence patterns across Magnaporthe pathotypes, which at least in part reflects the diversity in recognition mechanisms by the different hosts. Across pathotypes of M. oryzae , we show that TE insertion dynamics most likely underpin the observed effector rearrangements. In contrast, conserved AVRs show no recent TE activity in the surrounding regions. Variable AVR frequencies among M. oryzae host-specific pathotypes may be driven by adaptation to new hosts or the deployment of a new resistance gene. Here, we show that multiple MoO effectors were lost in different M. oryzae pathotypes including rice blast. This is consistent with the high rates of gene losses reported for rice-infecting M. oryzae compared to Triticum or Avena sp. infecting isolates ( Yoshida et al., 2016 ). How closely the AVR frequencies among hosts reflects pathogenicity and recognition capabilities remains to be determined. Effector loss may be sufficient to escape host recognition and therefore contributing to host range expansion. This has been reported for instance for PWT3 , with its loss coinciding with the widespread deployment of the complementary R gene RWT3 ( Inoue et al., 2017 ). In contrast, the more recently described MAX effectors do not necessarily reflect host specificity ( Naour-Vernet et al., 2023 ). This is consistent with the absence of certain effectors, which may not reflect avoidance of host recognition. Well-conserved effectors across multiple pathotypes and an origin outside of M. oryzae may reflect a conserved function in virulence maintained by purifying selection. However, contributions to virulence may still vary among pathotypes, as observed for the ACE1 effector ( Collemare et al., 2008 ; Vy et al., 2024 ). TEs represent approximately 10% of the M. oryzae genome ( Bao et al., 2017 ; Dean et al., 2005 ) and can be inserted in or around AVR genes altering their virulence spectrum through transcriptional silencing, loss-of-function or loss of avirulence ( Hu et al., 2022 ; Kang et al., 2001 ; Li et al., 2023 ; Li et al., 2009 ; Miki et al., 2009 ; Zhang et al., 2015 ). TE-mediated disruptions can result in the permanent loss of AVR genes, being a possible explanation for the observed deletion patterns. In AVR-Pib , a POT3 transposon insertion mediates effector loss-of-function in Philippine MoO isolates ( Olukayode et al., 2019 ). We found a striking association of the TE POT3 in the neighboring regions of AVR-Pib in MoO, suggesting that the TE plays a functional role. TE insertions can also promote the emergence of new virulent effector variants, as observed for the MGL retrotransposon insertion into the ACE1 gene ( Fudal et al., 2005 ) and POT3 insertion in AVR-Pita1 ( Kang et al., 2001 ; Zhou et al., 2007 ) and AVR-Piz-t ( Li et al., 2009 ) ensuring the maintenance of such effectors. In our collection, all the four MoO isolates exhibiting TEs in the ACE1 surrounding regions exhibited this same TE overlapping with the ACE1 coding sequence. The remaining 270 MoO isolates were devoid of TEs near the effector. TEs can also facilitate the translocation of AVR genes such as AVR-Pita ( Chuma et al., 2011 ). Here, we focused on MoO reference-quality genomes to investigate translocations of ACE1 and AVR-Pik . However, we found no association between these translocations and specific TEs. M. oryzae shows lineage-specific TE activity with LTR-retrotransposons expanded in certain lineages, rather than undergoing a single expansion followed by selective deletion ( Nakamoto et al., 2023 ). TEs were often found near genes with presence/absence variation including effectors ( Joubert & Krasileva, 2024 ). Here, we observed that well maintained AVRs had almost no TE activity in the surrounding area. This is possibly explained by their chromosomal location. In contrast to most of the M. oryzae AVR genes, usually located in telomeric or subtelomeric regions ( Chen et al., 2007 ; Chuma et al., 2010 , 2011 ; Orbach et al., 2000 ; Yoshida et al., 2009 ), highly conserved AVRs could have been selected to occupy euchromatic and repeat-poor regions. AVR-Pi9 , a widely present effector among M. oryzae pathotypes, is located in a genomic region close to the chromosome 7 centromere, which is one of the most stable regions of the genome ( Wu et al., 2015 ). Similarly, AVR-Piz-t is found at high frequency among pathotypes and located 230 kb from AVR-Pi9 ( Li et al., 2009 ). Effectors with highly variable presence/absence frequencies among pathotypes include AVR-Pia , AVR-Pii , AVR-Pita and PWL and are well-known to be located in unstable chromosomal regions such as subtelomeres prone to rearrangement ( Kang et al., 1995 ; Orbach et al., 2000 ; Sweigard et al., 1995 ; Yasuda et al., 2006 ; Yoshida et al., 2009 ). Despite the broad spectrum in chromosomal locations, most of the effectors are expressed during the early stages of infection. AVR-Pi9 accumulates in the biotrophic interfacial complex structure and is translocated in the early stages of infection to the host cell (Wu 2015). A similar expression and translocation pattern was also observed for less conserved effectors such as AVR-Pia and AVR-Pita ( Han et al., 2021 ; Khang et al., 2012 ; Sornkom et al., 2017 ). TE activity is typically high in telomeric and subtelomeric regions ( Croll & McDonald, 2012 ; Faino et al., 2016 ). Here, we show that AVRs with highly variable frequencies among pathotypes are indeed surrounded by TEs. The presence of TE can also affect epigenetic regulation through changes in heterochromatic structure ( Aparicio et al., 1991 ; De Las Peñas et al., 2003 ; Fan et al., 2008 ). These results reinforce the idea that in M. oryzae different TE environments impact effector evolution in distinct ways. Our findings also indicate that certain TE superfamilies may have more affinity for specific AVR s loci, reflected also in the fact that the TE presence is largely conserved among pathotypes. This suggests that despite differences in TE activity among fungal plant pathogens ( Gourlie et al., 2022 ; Oggenfuss et al., 2021 ; Shirke et al., 2016 ), including in M. oryzae ( Lin et al., 2024 ; Nakamoto et al., 2023 ), the effector regions may have been selected for a variety of genomic niche features. The clearest evidence for TE insertions creating transitional states at effector loci was found for AVR1-CO39 , with TEs likely contributing to the loss in MoO. Overall, we show that M. oryzae AVR locus evolution was characterized by parallel and well differentiation dynamics in TEs. Spanning the spectrum of effectors exhibiting rapid changes in frequencies across M. oryzae pathotypes to conserved effector loci retained at high frequency and devoid of TEs. This highlights that host-mediated selection plays not only a role in AVR frequencies but that the genomic niche of the effectors displays likely similar associated dynamics. The dynamics of TE insertions within plant pathogens as a response to host selection remains understudied. Our work demonstrates though that TE dynamics can be an integral component of genomic niche evolution. In-depth tracking of effector niches will likely augment our ability to predict the durability of host resistance. Declarations Competing interests The authors declare that no competing interests exist. Data availability The data analyzed in the frame of this study were retrieved from NCBI repositories as indicated in the Supplementary Tables S1 and S2. Draft genome assembly data is available from Zenodo ( https://doi.org/10.5281/zenodo.15496085 ). Funding This study was supported by a Swiss National Science Foundation grant to DC (201149). Acknowledgements We are grateful to group members for critical discussions. Funder Information Declared Swiss National Science Foundation, https://ror.org/00yjd3n13 , 201149 References ↵ Abraham , L. N. , Oggenfuss , U. , & Croll , D . ( 2024 ). Population-level transposable element expression dynamics influence trait evolution in a fungal crop pathogen . MBio , 15 ( 3 ). doi: 10.1128/MBIO.02840-23 OpenUrl CrossRef ↵ Anderson , P. K. , Cunningham , A. A. , Patel , N. G. , Morales , F. J. , Epstein , P. R. , & Daszak , P . ( 2004 ). Emerging infectious diseases of plants: pathogen pollution, climate change and agrotechnology drivers . Trends in Ecology & Evolution , 19 ( 10 ), 535 – 544 . doi: 10.1016/J.TREE.2004.07.021 OpenUrl CrossRef PubMed Anh , V. L. , Anh , N. T. , Tagle , A. G. , Vy , T. T. P. , Inoue , Y. , Takumi , S. , Chuma , I. , & Tosa , Y . ( 2015 ). Rmg8 , a new gene for resistance to Triticum isolates of Pyricularia oryzae in hexaploid wheat . Phytopathology , 105 ( 12 ), 1568 – 1572 . doi: 10.1094/PHYTO-02-15-0034-R/ASSET/IMAGES/LARGE/PHYTO-02-15-0034-R_F3-1450427894492.JPEG OpenUrl CrossRef PubMed Anh , V. L. , Inoue , Y. , Asuke , S. , Vy , T. T. P. , Anh , N. T. , Wang , S. , Chuma , I. , & Tosa , Y . ( 2018 ). Rmg8 and Rmg7 , wheat genes for resistance to the wheat blast fungus, recognize the same avirulence gene AVR-Rmg8 . Molecular Plant Pathology , 19 ( 5 ), 1252 – 1256 . doi: 10.1111/MPP.12609 OpenUrl CrossRef PubMed ↵ Aparicio , O. M. , Billington , B. L. , & Gottschling , D. E . ( 1991 ). Modifiers of position effect are shared between telomeric and silent mating-type loci in S. cerevisiae . Cell , 66 ( 6 ), 1279 – 1287 . doi: 10.1016/0092-8674(91)90049-5 OpenUrl CrossRef PubMed Web of Science Asuke , S. , Umehara , Y. , Inoue , Y. , Vy , T. T. P. , Iwakawa , M. , Matsuoka , Y. , Kato , K. , & Tosa , Y . ( 2021 ). Origin and Dynamics of Rwt6 , a Wheat Gene for Resistance to Nonadapted Pathotypes of Pyricularia oryzae . Phytopathology , 111 ( 11 ), 2023 – 2029 . doi: 10.1094/PHYTO-02-21-0080-R/ASSET/IMAGES/LARGE/PHYTO-02-21-0080-RF5-1639396102537.JPEG OpenUrl CrossRef PubMed ↵ Bankevich , A. , Nurk , S. , Antipov , D. , Gurevich , A. A. , Dvorkin , M. , Kulikov , A. S. , Lesin , V. M. , Nikolenko , S. I. , Pham , S. , Prjibelski , A. D. , Pyshkin , A. V. , Sirotkin , A. V. , Vyahhi , N. , Tesler , G. , Alekseyev , M. A. , & Pevzner , P. A . ( 2012 ). SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing . Journal of Computational Biology: A Journal of Computational Molecular Cell Biology , 19 ( 5 ), 455 – 477 . doi: 10.1089/CMB.2012.0021 OpenUrl CrossRef PubMed ↵ Bao , J. , Chen , M. , Zhong , Z. , Tang , W. , Lin , L. , Zhang , X. , Jiang , H. , Zhang , D. , Miao , C. , Tang , H. , Zhang , J. , Lu , G. , Ming , R. , Norvienyeku , J. , Wang , B. , & Wang , Z . ( 2017 ). PacBio Sequencing Reveals Transposable Elements as a Key Contributor to Genomic Plasticity and Virulence Variation in Magnaporthe oryzae . Molecular Plant , 10 ( 11 ), 1465 – 1468 . doi: 10.1016/J.MOLP.2017.08.008 OpenUrl CrossRef PubMed ↵ Böhnert , H. U. , Fudal , I. , Dioh , W. , Tharreau , D. , Notteghem , J. L. , & Lebrun , M. H . ( 2004 ). A putative polyketide synthase/peptide synthetase from Magnaporthe grisea signals pathogen attack to resistant rice . The Plant Cell , 16 ( 9 ), 2499 – 2513 . doi: 10.1105/TPC.104.022715 OpenUrl Abstract / FREE Full Text ↵ Borowiec , M. L . ( 2016 ). AMAS: a fast tool for alignment manipulation and computing of summary statistics . PeerJ , 4 ( 1 ). doi: 10.7717/PEERJ.1660 OpenUrl CrossRef ↵ Camacho , C. , Coulouris , G. , Avagyan , V. , Ma , N. , Papadopoulos , J. , Bealer , K. , & Madden , T. L . ( 2009 ). BLAST+: Architecture and applications . BMC Bioinformatics , 10 ( 1 ), 1 – 9 . doi: 10.1186/1471-2105-10-421/FIGURES/4 OpenUrl CrossRef PubMed ↵ Chen , Q. H. , Wang , Y. C. , Li , A. N. , Zhang , Z. G. , & Zheng , X. B . ( 2007 ). Molecular mapping of two cultivar-specific avirulence genes in the rice blast fungus Magnaporthe grisea . Molecular Genetics and Genomics: MGG , 277 ( 2 ), 139 – 148 . doi: 10.1007/S00438-006-0179-8 OpenUrl CrossRef PubMed ↵ Chen , S. , Zhou , Y. , Chen , Y. , & Gu , J . ( 2018 ). fastp: an ultra-fast all-in-one FASTQ preprocessor . Bioinformatics , 34 ( 17 ), i884 – i890 . doi: 10.1093/BIOINFORMATICS/BTY560 OpenUrl CrossRef ↵ Choi , J. , Park , S. Y. , Kim , B. R. , Roh , J. H. , Oh , I. S. , Han , S. S. , & Lee , Y. H . ( 2013 ). Comparative analysis of pathogenicity and phylogenetic relationship in Magnaporthe grisea species complex . PLOS ONE , 8 ( 2 ), e57196 . doi: 10.1371/JOURNAL.PONE.0057196 OpenUrl CrossRef PubMed ↵ Chuma , I. , Isobe , C. , Hotta , Y. , Ibaragi , K. , Futamata , N. , Kusaba , M. , Yoshida , K. , Terauchi , R. , Fujita , Y. , Nakayashiki , H. , Valent , B. , & Tosa , Y . ( 2011 ). Multiple Translocation of the AVR-Pita Effector Gene among Chromosomes of the Rice Blast Fungus Magnaporthe oryzae and Related Species . PLoS Pathogens , 7 ( 7 ). doi: 10.1371/JOURNAL.PPAT.1002147 OpenUrl CrossRef ↵ Chuma , I. , Zhan , S. W. , Asano , S. , Nga , N. T. T. , Vy , T. T. P. , Shirai , M. , Ibaragi , K. , & Tosa , Y . ( 2010 ). PWT1 , an avirulence gene of Magnaporthe oryzae tightly linked to the rDNA Locus, is recognized by two staple crops, common wheat and barley . Phytopathology , 100 ( 5 ), 436 – 443 . doi: 10.1094/PHYTO-100-5-0436 OpenUrl CrossRef PubMed ↵ Collemare , J. , Pianfetti , M. , Houlle , A. E. , Morin , D. , Camborde , L. , Gagey , M. J. , Barbisan , C. , Fudal , I. , Lebrun , M. H. , & Böhnert , H. U . ( 2008 ). Magnaporthe grisea avirulence gene ACE1 belongs to an infection-specific gene cluster involved in secondary metabolism . The New Phytologist , 179 ( 1 ), 196 – 208 . doi: 10.1111/J.1469-8137.2008.02459.X OpenUrl CrossRef PubMed ↵ Couch , B. C. , Fudal , I. , Lebrun , M. H. , Tharreau , D. , Valent , B. , Van Kim , P. , Nottéghem , J. L. , & Kohn , L. M. ( 2005 ). Origins of Host-Specific Populations of the Blast Pathogen Magnaporthe oryzae in Crop Domestication With Subsequent Expansion of Pandemic Clones on Rice and Weeds of Rice . Genetics , 170 ( 2 ), 613 . doi: 10.1534/GENETICS.105.041780 OpenUrl Abstract / FREE Full Text ↵ Croll , D. , & McDonald , B. A . ( 2012 ). The Accessory Genome as a Cradle for Adaptive Evolution in Pathogens . PLOS Pathogens , 8 ( 4 ), e1002608 . doi: 10.1371/JOURNAL.PPAT.1002608 OpenUrl CrossRef PubMed ↵ Dai , Y. , Jia , Y. , Correll , J. , Wang , X. , & Wang , Y . ( 2010 ). Diversification and evolution of the avirulence gene AVR-Pita1 in field isolates of Magnaporthe oryzae . Fungal Genetics and Biology: FG & B , 47 ( 12 ), 973 – 980 . doi: 10.1016/J.FGB.2010.08.003 OpenUrl CrossRef ↵ De Jonge , R. , Bolton , M. D. , Kombrink , A. , Van Den Berg , G. C. M. , Yadeta , K. A. , & Thomma , B. P. H. J. ( 2013 ). Extensive chromosomal reshuffling drives evolution of virulence in an asexual pathogen . Genome Research , 23 ( 8 ), 1271 – 1282 . doi: 10.1101/GR.152660.112 OpenUrl Abstract / FREE Full Text ↵ De Las Peñas , A. , Pan , S. J. , Castaño , I. , Alder , J. , Cregg , R. , & Cormack , B. P . ( 2003 ). Virulence-related surface glycoproteins in the yeast pathogen Candida glabrata are encoded in subtelomeric clusters and subject to RAP1 - and SIR -dependent transcriptional silencing . Genes & Development , 17 ( 18 ), 2245 . doi: 10.1101/GAD.1121003 OpenUrl Abstract / FREE Full Text ↵ De Wit , P. J. G. M. , Mehrabi , R. , Van Den Burg , H. A. , & Stergiopoulos , I. ( 2009 ). Fungal effector proteins: past, present and future . Molecular Plant Pathology , 10 ( 6 ), 735 – 747 . doi: 10.1111/J.1364-3703.2009.00591.X OpenUrl CrossRef PubMed ↵ Dean , R. A. , Talbot , N. J. , Ebbole , D. J. , Farman , M. L. , Mitchell , T. K. , Orbach , M. J. , Thon , M. , Kulkarni , R. , Xu , J. R. , Pan , H. , Read , N. D. , Lee , Y. I. , Carbone , I. , Brown , D. , Yeon , Y. O. , Donofrio , N. , Jun , S. J. , Soanes , D. M. , Djonovic , S. , … Dirren , B. W . ( 2005 ). The genome sequence of the rice blast fungus Magnaporthe grisea . Nature , 434 ( 7036 ), 980 – 986 . doi: 10.1038/NATURE03449 OpenUrl CrossRef PubMed Web of Science ↵ 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 . doi: 10.1016/J.GDE.2015.09.001 OpenUrl CrossRef PubMed ↵ Emms , D. M. , & Kelly , S . ( 2019 ). OrthoFinder: Phylogenetic orthology inference for comparative genomics . Genome Biology , 20 ( 1 ), 1 – 14 . doi: 10.1186/S13059-019-1832-Y/FIGURES/5 OpenUrl CrossRef ↵ Faino , L. , Seidl , M. F. , Shi-Kunne , X. , Pauper , M. , Van Den Berg , G. C. M. , Wittenberg , A. H. J. , & Thomma , B. P. H. J. ( 2016 ). Transposons passively and actively contribute to evolution of the two-speed genome of a fungal pathogen . Genome Research , 26 ( 8 ), 1091 – 1100 . doi: 10.1101/GR.204974.116/-/DC1 OpenUrl Abstract / FREE Full Text ↵ Fan , C. , Zhang , Y. , Yu , Y. , Rounsley , S. , Long , M. , & Wing , R. A . ( 2008 ). The subtelomere of Oryza sativa chromosome 3 short arm as a hot bed of new gene origination in rice . Molecular Plant , 1 ( 5 ), 839 . doi: 10.1093/MP/SSN050 OpenUrl CrossRef PubMed Web of Science ↵ Farman , M. L . ( 2007 ). Telomeres in the rice blast fungus Magnaporthe oryzae : the world of the end as we know it . FEMS Microbiology Letters , 273 ( 2 ), 125 – 132 . doi: 10.1111/J.1574-6968.2007.00812.X OpenUrl CrossRef PubMed ↵ Farman , M. L. , Eto , Y. , Nakao , T. , Tosa , Y. , Nakayashiki , H. , Mayama , S. , & Leong , S. A . ( 2002 ). Analysis of the structure of the AVR1-CO39 avirulence locus in virulent rice-infecting isolates of Magnaporthe grisea . Molecular Plant-Microbe Interactions , 15 ( 1 ), 6 – 16 . doi: 10.1094/MPMI.2002.15.1.6 OpenUrl CrossRef PubMed Web of Science Farman , M. L. , & Leong , S. A . ( 1998 ). Chromosome Walking to the AVR1-CO39 Avirulence Gene of Magnaporthe grisea : Discrepancy Between the Physical and Genetic Maps . Genetics , 150 ( 3 ), 1049 – 1058 . doi: 10.1093/GENETICS/150.3.1049 OpenUrl Abstract / FREE Full Text ↵ Fouché , S. , Plissonneau , C. , & Croll , D . ( 2018 ). The birth and death of effectors in rapidly evolving filamentous pathogen genomes . Current Opinion in Microbiology , 46 , 34 – 42 . doi: 10.1016/J.MIB.2018.01.020 OpenUrl CrossRef PubMed ↵ Fudal , I. , Böhnert , H. U. , Tharreau , D. , & Lebrun , M. H . ( 2005 ). Transposition of MINE, a composite retrotransposon, in the avirulence gene ACE1 of the rice blast fungus Magnaporthe grisea . Fungal Genetics and Biology: FG & B , 42 ( 9 ), 761 – 772 . doi: 10.1016/J.FGB.2005.05.001 OpenUrl CrossRef ↵ Gladieux , P. , Condon , B. , Ravel , S. , Soanes , D. , Maciel , J. L. N. , Nhani , A. , Chen , L. , Terauchi , R. , Lebrun , M. H. , Tharreau , D. , Mitchell , T. , Pedley , K. F. , Valent , B. , Talbot , N. J. , Farman , M. , & Fournier , E . ( 2018 ). Gene flow between divergent cereal- and grass-specific lineages of the rice blast fungus Magnaporthe oryzae . MBio , 9 ( 1 ). doi: 10.1128/MBIO.01219-17/SUPPL_FILE/MBO001183743ST3.PDF OpenUrl CrossRef ↵ Gomez-Luciano , L. B. , Tsai , I. J. , Chuma , I. , Tosa , Y. , Chen , Y. H. , Li , J. Y. , Li , M. Y. , Lu , M. Y. J. , Nakayashiki , H. , & Li , W. H . ( 2019 ). Blast Fungal Genomes Show Frequent Chromosomal Changes, Gene Gains and Losses, and Effector Gene Turnover . Molecular Biology and Evolution , 36 ( 6 ), 1148 – 1161 . doi: 10.1093/MOLBEV/MSZ045 OpenUrl CrossRef PubMed ↵ Gourlie , R. , McDonald , M. , Hafez , M. , Ortega-Polo , R. , Low , K. E. , Abbott , D. W. , Strelkov , S. E. , Daayf , F. , & Aboukhaddour , R . ( 2022 ). The pangenome of the wheat pathogen Pyrenophora tritici-repentis reveals novel transposons associated with necrotrophic effectors ToxA and ToxB . BMC Biology , 20 ( 1 ), 239 . doi: 10.1186/S12915-022-01433-W OpenUrl CrossRef PubMed ↵ Gurevich , A. , Saveliev , V. , Vyahhi , N. , & Tesler , G . ( 2013 ). QUAST: quality assessment tool for genome assemblies. Bioinformatics (Oxford , England ) , 29 ( 8 ), 1072 – 1075 . doi: 10.1093/BIOINFORMATICS/BTT086 OpenUrl CrossRef ↵ Guy , L. , Kultima , J. R. , Andersson , S. G. E. , & Quackenbush , J . ( 2010 ). genoPlotR: comparative gene and genome visualization in R . Bioinformatics (Oxford, England) , 26 ( 18 ), 2334 – 2335 . doi: 10.1093/BIOINFORMATICS/BTQ413 OpenUrl CrossRef PubMed Web of Science ↵ Han , J. , Wang , X. , Wang , F. , Zhao , Z. , Li , G. , Zhu , X. , Su , J. , & Chen , L . ( 2021 ). The Fungal Effector Avr-Pita Suppresses Innate Immunity by Increasing COX Activity in Rice Mitochondria . Rice , 14 ( 1 ), 1 – 11 . doi: 10.1186/S12284-021-00453-4/FIGURES/7 OpenUrl CrossRef PubMed ↵ Hu , Z. J. , Huang , Y. Y. , Lin , X. Y. , Feng , H. , Zhou , S. X. , Xie , Y. , Liu , X. X. , Liu , C. , Zhao , R. M. , Zhao , W. S. , Feng , C. H. , Pu , M. , Ji , Y. P. , Hu , X. H. , Li , G. B. , Zhao , J. H. , Zhao , Z. X. , Wang , H. , Zhang , J. W. , … Wang , W. M . ( 2022 ). Loss and Natural Variations of Blast Fungal Avirulence Genes Breakdown Rice Resistance Genes in the Sichuan Basin of China . Frontiers in Plant Science , 13 , 788876 . doi: 10.3389/FPLS.2022.788876/BIBTEX OpenUrl CrossRef PubMed ↵ Huang , J. , Si , W. , Deng , Q. , Li , P. , & Yang , S . ( 2014 ). Rapid evolution of avirulence genes in rice blast fungus Magnaporthe oryzae . BMC Genetics , 15 , 45 . doi: 10.1186/1471-2156-15-45 OpenUrl CrossRef PubMed ↵ Inoue , Y. , Vy , T. T. P. , Tani , D. , & Tosa , Y . ( 2021 ). Suppression of wheat blast resistance by an effector of Pyricularia oryzae is counteracted by a host specificity resistance gene in wheat . New Phytologist , 229 ( 1 ), 488 – 500 . doi: 10.1111/NPH.16894 OpenUrl CrossRef PubMed ↵ Inoue , Y. , Vy , T. T. P. , Yoshida , K. , Asano , H. , Mitsuoka , C. , Asuke , S. , Anh , V. L. , Cumagun , C. J. R. , Chuma , I. , Terauchi , R. , Kato , K. , Mitchell , T. , Valent , B. , Farman , M. , & Tosa , Y . ( 2017 ). Evolution of the wheat blast fungus through functional losses in a host specificity determinant . Science (New York, N.Y.) , 357 ( 6346 ), 80 – 83 . doi: 10.1126/SCIENCE.AAM9654 OpenUrl CrossRef ↵ Jeon , J. , Choi , J. , Lee , G. W. , Park , S. Y. , Huh , A. , Dean , R. A. , & Lee , Y. H . ( 2015 ). Genome-wide profiling of DNA methylation provides insights into epigenetic regulation of fungal development in a plant pathogenic fungus, Magnaporthe oryzae . Scientific Reports , 5 , 8567 . doi: 10.1038/SREP08567 OpenUrl CrossRef PubMed ↵ Jia , Y. , McAdams , S. A. , Bryan , G. T. , Hershey , H. P. , & Valent , B . ( 2000 ). Direct interaction of resistance gene and avirulence gene products confers rice blast resistance . The EMBO Journal , 19 ( 15 ), 4004 – 4014 . doi: 10.1093/EMBOJ/19.15.4004 OpenUrl Abstract / FREE Full Text ↵ Joubert , P. M. , & Krasileva , K. V . ( 2024 ). Distinct genomic contexts predict gene presence–absence variation in different pathotypes of Magnaporthe oryzae . Genetics , 226 ( 4 ), iyae012 . doi: 10.1093/GENETICS/IYAE012 OpenUrl CrossRef PubMed ↵ Kang , S. , Lebrun , M. H. , Farrall , L. , & Valent , B . ( 2001 ). Gain of virulence caused by insertion of a Pot3 transposon in a Magnaporthe grisea avirulence gene . Molecular Plant-Microbe Interactions: MPMI , 14 ( 5 ), 671 – 674 . doi: 10.1094/MPMI.2001.14.5.671 OpenUrl CrossRef PubMed Web of Science ↵ Kang , S. , Sweigard , J. A. , & Valent , B . ( 1995 ). The PWL host specificity gene family in the blast fungus Magnaporthe grisea . Molecular Plant-Microbe Interactions: MPMI , 8 ( 6 ), 939 – 948 . doi: 10.1094/MPMI-8-0939 OpenUrl CrossRef PubMed Web of Science ↵ Kanyuka , K. , & Rudd , J. J . ( 2019 ). Cell surface immune receptors: the guardians of the plant’s extracellular spaces . Current Opinion in Plant Biology , 50 , 1 – 8 . doi: 10.1016/J.PBI.2019.02.005 OpenUrl CrossRef PubMed ↵ Kanzaki , H. , Yoshida , K. , Saitoh , H. , Fujisaki , K. , Hirabuchi , A. , Alaux , L. , Fournier , E. , Tharreau , D. , & Terauchi , R . ( 2012 ). Arms race co-evolution of Magnaporthe oryzae AVR-Pik and rice Pik genes driven by their physical interactions . The Plant Journal , 72 ( 6 ), 894 – 907 . doi: 10.1111/J.1365-313X.2012.05110.X OpenUrl CrossRef PubMed ↵ Khang , C. H. , Berruyer , R. , Giraldo , M. C. , Kankanala , P. , Park , S. Y. , Czymmek , K. , Kang , S. , & Valent , B . ( 2012 ). Translocation of Magnaporthe oryzae Effectors into Rice Cells and Their Subsequent Cell-to-Cell Movement . The Plant Cell , 22 ( 4 ), 1388 – 1403 . doi: 10.1105/TPC.109.069666 OpenUrl CrossRef ↵ Khang , C. H. , Park , S. Y. , Lee , Y. H. , Valent , B. , & Kang , S . ( 2008 ). Genome organization and evolution of the AVR-Pita avirulence gene family in the Magnaporthe grisea species complex . Molecular Plant-Microbe Interactions: MPMI , 21 ( 5 ), 658 – 670 . doi: 10.1094/MPMI-21-5-0658 OpenUrl CrossRef PubMed Web of Science ↵ Kim , K. T. , Ko , J. , Song , H. , Choi , G. , Kim , H. , Jeon , J. , Cheong , K. , Kang , S. , & Lee , Y. H . ( 2019 ). Evolution of the Genes Encoding Effector Candidates Within Multiple Pathotypes of Magnaporthe oryzae . Frontiers in Microbiology , 10 , 491941 . doi: 10.3389/FMICB.2019.02575/BIBTEX OpenUrl CrossRef ↵ Latorre , S. M. , Were , V. M. , Foster , A. J. , Langner , T. , Malmgren , A. , Harant , A. , Asuke , S. , Reyes-Avila , S. , Gupta , D. R. , Jensen , C. , Ma , W. , Mahmud , N. U. , Mehebub , M. S. , Mulenga , R. M. , Muzahid , A. N. M. , Paul , S. K. , Rabby , S. M. F. , Al Mahbub Rahat , A. , Ryder , L. , … Kamoun , S. ( 2023 ). Genomic surveillance uncovers a pandemic clonal lineage of the wheat blast fungus . PLOS Biology , 21 ( 4 ), e3002052 . doi: 10.1371/JOURNAL.PBIO.3002052 OpenUrl CrossRef PubMed ↵ Li , J. , Lu , L. , Li , C. , Wang , Q. , & Shi , Z . ( 2023 ). Insertion of Transposable Elements in AVR-Pib of Magnaporthe oryzae Leading to LOSS of the Avirulent Function . International Journal of Molecular Sciences , 24 ( 21 ), 15542 . doi: 10.3390/IJMS242115542/S1 OpenUrl CrossRef PubMed ↵ Li , W. , Wang , B. , Wu , J. , Lu , G. , Hu , Y. , Zhang , X. , Zhang , Z. , Zhao , Q. , Feng , Q. , Zhang , H. , Wang , Z. , Wang , G. L. , Han , B. , Wang , Z. , & Zhou , B . ( 2009 ). The Magnaporthe oryzae Avirulence Gene The AvrPiz-t Encodes a Predicted Secreted Protein That Triggers the Immunity in Rice Mediated by the Blast Resistance Gene Piz-t . doi: 10.1094/MPMI-22-4-0411 , 22 ( 4 ), 411 – 420 . 10.1094/MPMI-22-4-0411 OpenUrl CrossRef PubMed Web of Science ↵ Lin , L. , Sun , T. , Guo , J. , Lin , L. , Chen , M. , Wang , Z. , Bao , J. , Norvienyeku , J. , Zhang , D. , Han , Y. , Lu , G. , Rensing , C. , Zheng , H. , Zhong , Z. , & Wang , Z . ( 2024 ). Transposable elements impact the population divergence of rice blast fungus Magnaporthe oryzae . MBio , 15 ( 5 ). doi: 10.1128/MBIO.00086-24 OpenUrl CrossRef ↵ Liu , S. , Lin , G. , Ramachandran , S. R. , Daza , L. C. , Cruppe , G. , Tembo , B. , Singh , P. K. , Cook , D. , Pedley , K. F. , & Valent , B . ( 2024 ). Rapid mini-chromosome divergence among fungal isolates causing wheat blast outbreaks in Bangladesh and Zambia . New Phytologist , 241 ( 3 ), 1266 – 1276 . doi: 10.1111/NPH.19402 OpenUrl CrossRef PubMed ↵ Lo Presti , L. , Lanver , D. , Schweizer , G. , Tanaka , S. , Liang , L. , Tollot , M. , Zuccaro , A. , Reissmann , S. , & Kahmann , R. ( 2015 ). Fungal Effectors and Plant Susceptibility . Annual Review of Plant Biology , 66 , 513 – 545 . doi: 10.1146/ANNUREV-ARPLANT-043014-114623 OpenUrl CrossRef PubMed ↵ Longya , A. , Chaipanya , C. , Franceschetti , M. , Maidment , J. H. R. , Banfield , M. J. , & Jantasuriyarat , C . ( 2019 ). Gene duplication and mutation in the emergence of a novel aggressive allele of the AVR-PIK effector in the rice blast fungus . Molecular Plant-Microbe Interactions , 32 ( 6 ), 740 – 749 . doi: 10.1094/MPMI-09-18-0245-R/ASSET/IMAGES/LARGE/MPMI-09-18-0245-R_F8-1557354838425.JPEG OpenUrl CrossRef PubMed ↵ Miki , S. , Matsui , K. , Kito , H. , Otsuka , K. , Ashizawa , T. , Yasuda , N. , Fukiya , S. , Sato , J. , Hirayae , K. , Fujita , Y. , Nakajima , T. , Tomita , F. , & Sone , T . ( 2009 ). Molecular cloning and characterization of the AVR-Pia locus from a Japanese field isolate of Magnaporthe oryzae . Molecular Plant Pathology , 10 ( 3 ), 361 – 374 . doi: 10.1111/J.1364-3703.2009.00534.X OpenUrl CrossRef PubMed ↵ Nakamoto , A. A. , Joubert , P. M. , & Krasileva , K. V . ( 2023 ). Intraspecific Variation of Transposable Elements Reveals Differences in the Evolutionary History of Fungal Phytopathogen Pathotypes . Genome Biology and Evolution , 15 ( 12 ). doi: 10.1093/GBE/EVAD206 OpenUrl CrossRef ↵ Naour-Vernet , M. Le , Charriat , F. , Gracy , J. , Cros-Arteil , S. , Ravel , S. , Veillet , F. , Meusnier , I. , Padilla , A. , Kroj , T. , Cesari , S. , & Gladieux , P. ( 2023 ). Adaptive evolution in virulence effectors of the rice blast fungus Pyricularia oryzae . PLOS Pathogens , 19 ( 9 ), e1011294 . doi: 10.1371/JOURNAL.PPAT.1011294 OpenUrl CrossRef PubMed ↵ Oggenfuss , U. , Badet , T. , Wicker , T. , Hartmann , F. E. , Singh , N. K. , Abraham , L. , Karisto , P. , Vonlanthen , T. , Mundt , C. , McDonald , B. A. , & Croll , D . ( 2021 ). A population-level invasion by transposable elements triggers genome expansion in a fungal pathogen . ELife , 10 . doi: 10.7554/ELIFE.69249 OpenUrl CrossRef ↵ Olukayode , T. , Quime , B. , Shen , Y. C. , Yanoria , M. J. , Zhang , S. , Yang , J. , Zhu , X. , Shen , W. C. , von Tiedemann , A. , & Zhou , B. ( 2019 ). Dynamic insertion of POT3 in AvrPib prevailing in a field rice blast population in the Philippines led to the high virulence frequency against the resistance gene Pib in rice . Phytopathology , 109 ( 5 ), 870 – 877 . doi: 10.1094/PHYTO-06-18-0198-R/ASSET/IMAGES/LARGE/PHYTO-06-18-0198-R_F5-1555764579507.JPEG OpenUrl CrossRef PubMed ↵ Onaga , G. , Suktrakul , W. , Wanjiku , M. , Lorenzo Quibod , I. , Baka , J .-4, Entfellner , D. , Bigirimana , J. , Habarugira , G. , Murori , R. , Asea , G. , Ismail , A. M. , Jantasuriyarat , C. , & Oliva , R. ( 2020 ). Magnaporthe oryzae populations in Sub-Saharan Africa are diverse and show signs of local adaptation . BioRxiv , 2020.11.17.377325 . doi: 10.1101/2020.11.17.377325 OpenUrl Abstract / FREE Full Text ↵ Orbach , M. J. , Farrall , L. , Sweigard , J. A. , Chumley , F. G. , & Valent , B . ( 2000 ). A Telomeric Avirulence Gene Determines Efficacy for the Rice Blast Resistance Gene Pi-ta . The Plant Cell , 12 ( 11 ), 2019 – 2032 . doi: 10.1105/TPC.12.11.2019 OpenUrl Abstract / FREE Full Text ↵ Park , C. H. , Chen , S. , Shirsekar , G. , Zhou , B. , Khang , C. H. , Songkumarn , P. , Afzal , A. J. , Ning , Y. , Wang , R. , Bellizzi , M. , Valent , B. , & Wang , G. L . ( 2012 ). The Magnaporthe oryzae Effector AvrPiz-t Targets the RING E3 Ubiquitin Ligase APIP6 to Suppress Pathogen-Associated Molecular Pattern–Triggered Immunity in Rice . The Plant Cell , 24 ( 11 ), 4748 – 4762 . doi: 10.1105/TPC.112.105429 OpenUrl Abstract / FREE Full Text ↵ Peng , Z. , Oliveira-Garcia , E. , Lin , G. , Hu , Y. , Dalby , M. , Migeon , P. , Tang , H. , Farman , M. , Cook , D. , White , F. F. , Valent , B. , & Liu , S . ( 2019 ). Effector gene reshuffling involves dispensable mini-chromosomes in the wheat blast fungus . PLoS Genetics , 15 ( 9 ). doi: 10.1371/JOURNAL.PGEN.1008272 OpenUrl CrossRef ↵ Peyyala , R. , & Farman , M. L . ( 2006 ). Magnaporthe oryzae isolates causing gray leaf spot of perennial ryegrass possess a functional copy of the AVR1-CO39 avirulence gene . Molecular Plant Pathology , 7 ( 3 ), 157 – 165 . doi: 10.1111/J.1364-3703.2006.00325.X OpenUrl CrossRef PubMed ↵ Rahnama , M. , Condon , B. , Ascari , J. P. , Dupuis , J. R. , Del Ponte , E. M. , Pedley , K. F. , Martinez , S. , Valent , B. , & Farman , M. L . ( 2023 ). Recent co-evolution of two pandemic plant diseases in a multi-hybrid swarm . Nature Ecology & Evolution 2023 7:12 , 7 ( 12 ), 2055 – 2066 . doi: 10.1038/s41559-023-02237-z OpenUrl CrossRef PubMed Ray , S. , Singh , P. K. , Gupta , D. K. , Mahato , A. K. , Sarkar , C. , Rathour , R. , Singh , N. K. , & Sharma , T. R . ( 2016 ). Analysis of Magnaporthe oryzae genome reveals a fungal effector, which is able to induce resistance response in transgenic rice line containing resistance gene, Pi54 . Frontiers in Plant Science , 7 ( AUG2016 ), 210426 . doi: 10.3389/FPLS.2016.01140/BIBTEX OpenUrl CrossRef ↵ Rehmeyer , C. , Li , W. , Kusaba , M. , Kim , Y. S. , Brown , D. , Staben , C. , Dean , R. , & Farman , M . ( 2006 ). Organization of chromosome ends in the rice blast fungus, Magnaporthe oryzae . Nucleic Acids Research , 34 ( 17 ), 4685 . doi: 10.1093/NAR/GKL588 OpenUrl CrossRef PubMed Web of Science ↵ Remick , B. C. , Gaidt , M. M. , & Vance , R. E . ( 2023 ). Effector-Triggered Immunity . Annual Review of Immunology , 41 , 453 – 481 . doi: 10.1146/ANNUREV-IMMUNOL-101721-031732 OpenUrl CrossRef PubMed ↵ Ribot , C. , Césari , S. , Abidi , I. , Chalvon , V. , Bournaud , C. , Vallet , J. , Lebrun , M. H. , Morel , J. B. , & Kroj , T . ( 2013 ). The Magnaporthe oryzae effector AVR1-CO39 is translocated into rice cells independently of a fungal-derived machinery . The Plant Journal: For Cell and Molecular Biology , 74 ( 1 ), 1 – 12 . doi: 10.1111/TPJ.12099 OpenUrl CrossRef ↵ Sampaio , A. M. , Tralamazza , S. M. , Mohamadi , F. , De Oliveira , Y. , Enjalbert , J. , Saintenac , C. , & Croll , D. ( 2025 ). Diversification, loss, and virulence gains of the major effector AvrStb6 during continental spread of the wheat pathogen Zymoseptoria tritici . PLOS Pathogens , 21 ( 3 ), e1012983 . doi: 10.1371/JOURNAL.PPAT.1012983 OpenUrl CrossRef PubMed ↵ Sánchez-Vallet , A. , Fouché , S. , Fudal , I. , Hartmann , F. E. , Soyer , J. L. , Tellier , A. , & Croll , D . ( 2018 ). The Genome Biology of Effector Gene Evolution in Filamentous Plant Pathogens . Annual Review of Phytopathology , 56 , 21 – 40 . doi: 10.1146/ANNUREV-PHYTO-080516-035303 OpenUrl CrossRef PubMed ↵ Seidl , M. F. , & Thomma , B. P. H. J . ( 2017 ). Transposable Elements Direct The Coevolution between Plants and Microbes . Trends in Genetics , 33 ( 11 ), 842 – 851 . doi: 10.1016/J.TIG.2017.07.003 OpenUrl CrossRef PubMed ↵ Shi-Kunne , X. , Faino , L. , van den Berg , G. C. M. , Thomma , B. P. H. J. , & Seidl , M. F. ( 2018 ). Evolution within the fungal genus Verticillium is characterized by chromosomal rearrangement and gene loss . Environmental Microbiology , 20 ( 4 ), 1362 – 1373 . doi: 10.1111/1462-2920.14037 OpenUrl CrossRef ↵ Shirke , M. D. , Mahesh , H. B. , & Gowda , M . ( 2016 ). Genome-Wide Comparison of Magnaporthe Species Reveals a Host-Specific Pattern of Secretory Proteins and Transposable Elements . PLoS ONE , 11 ( 9 ), e0162458 . doi: 10.1371/JOURNAL.PONE.0162458 OpenUrl CrossRef PubMed ↵ Simão , F. A. , Waterhouse , R. M. , Ioannidis , P. , Kriventseva , E. V. , & Zdobnov , E. M . ( 2015 ). BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs . Bioinformatics (Oxford, England) , 31 ( 19 ), 3210 – 3212 . doi: 10.1093/BIOINFORMATICS/BTV351 OpenUrl CrossRef PubMed ↵ Sironi , M. , Cagliani , R. , Forni , D. , & Clerici , M . ( 2015 ). Evolutionary insights into host–pathogen interactions from mammalian sequence data . Nature Reviews Genetics 2015 16:4 , 16 ( 4 ), 224 – 236 . doi: 10.1038/nrg3905 OpenUrl CrossRef PubMed ↵ Sornkom , W. , Miki , S. , Takeuchi , S. , Abe , A. , Asano , K. , & Sone , T . ( 2017 ). Fluorescent reporter analysis revealed the timing and localization of AVR-Pia expression, an avirulence effector of Magnaporthe oryzae . Molecular Plant Pathology , 18 ( 8 ), 1138 – 1149 . doi: 10.1111/MPP.12468 OpenUrl CrossRef PubMed ↵ Stamatakis , A . ( 2014 ). RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies . Bioinformatics , 30 ( 9 ), 1312 – 1313 . doi: 10.1093/BIOINFORMATICS/BTU033 OpenUrl CrossRef PubMed Web of Science ↵ Stanke , M. , Steinkamp , R. , Waack , S. , & Morgenstern , B . ( 2004 ). AUGUSTUS: a web server for gene finding in eukaryotes . Nucleic Acids Research , 32 ( Web Server issue ), W309 . doi: 10.1093/NAR/GKH379 OpenUrl CrossRef PubMed Web of Science ↵ Stukenbrock , E. H. , & McDonald , B. A . ( 2009 ). Population genetics of fungal and oomycete effectors involved in gene-for-gene interactions . Molecular Plant-Microbe Interactions , 22 ( 4 ), 371 – 380 . doi: 10.1094/MPMI-22-4-0371 OpenUrl CrossRef PubMed Web of Science ↵ Sweigard , J. A. , Carroll , A. M. , Kang , S. , Farrall , L. , Chumley , F. G. , & Valent , B . ( 1995 ). Identification, cloning, and characterization of PWL2, a gene for host species specificity in the rice blast fungus . The Plant Cell , 7 ( 8 ), 1221 – 1233 . doi: 10.1105/TPC.7.8.1221 OpenUrl Abstract / FREE Full Text ↵ Thierry , M. , Charriat , F. , Milazzo , J. , Adreit , H. , Ravel , S. , Cros-Arteil , S. , Borron , S. , Sella , V. , Kroj , T. , Ioos , R. , Fournier , E. , Tharreau , D. , & Gladieux , P . ( 2022 ). Maintenance of divergent lineages of the Rice Blast Fungus Pyricularia oryzae through niche separation, loss of sex and post-mating genetic incompatibilities . PLoS Pathogens , 18 ( 7 ). doi: 10.1371/JOURNAL.PPAT.1010687 OpenUrl CrossRef ↵ Tosa , Y. , Osue , J. , Eto , Y. , Oh , H. S. , Nakayashiki , H. , Mayama , S. , & Leong , S. A . ( 2005 ). Evolution of an avirulence gene, AVR1 -CO39, concomitant with the evolution and differentiation of Magnaporthe oryzae . Molecular Plant-Microbe Interactions: MPMI , 18 ( 11 ), 1148 – 1160 . doi: 10.1094/MPMI-18-1148 OpenUrl CrossRef PubMed Web of Science ↵ Valent , B . ( 1990 ). Rice Blast as a Model System for Plant Pathology . Phytopathology , 80 ( 1 ), 33 . doi: 10.1094/PHYTO-80-33 OpenUrl CrossRef Web of Science ↵ Vy , T. T. P. , Inoue , Y. , Asuke , S. , Chuma , I. , Nakayashiki , H. , & Tosa , Y . ( 2024 ). The ACE1 secondary metabolite gene cluster is a pathogenicity factor of wheat blast fungus . Communications Biology 2024 7:1 , 7 ( 1 ), 1 – 10 . doi: 10.1038/s42003-024-06517-7 OpenUrl CrossRef PubMed ↵ Wang , Q. , Sun , M. , Zhang , Y. , Song , Z. , Zhang , S. , Zhang , Q. , Xu , J. R. , & Liu , H . ( 2020 ). Extensive chromosomal rearrangements and rapid evolution of novel effector superfamilies contribute to host adaptation and speciation in the basal ascomycetous fungi . Molecular Plant Pathology , 21 ( 3 ), 330 – 348 . doi: 10.1111/MPP.12899 OpenUrl CrossRef PubMed ↵ Whisson , S. C. , Vetukuri , R. R. , Avrova , A. O. , & Dixelius , C . ( 2012 ). Can silencing of transposons contribute to variation in effector gene expression in Phytophthora infestans ? Mobile Genetic Elements , 2 ( 2 ), 110 . doi: 10.4161/MGE.20265 OpenUrl CrossRef PubMed ↵ Wilson , R. A. , & Talbot , N. J . ( 2009 ). Under pressure: investigating the biology of plant infection by Magnaporthe oryzae . Nature Reviews. Microbiology , 7 ( 3 ), 185 – 195 . doi: 10.1038/NRMICRO2032 OpenUrl CrossRef PubMed Web of Science ↵ Wu , J. , Kou , Y. , Bao , J. , Li , Y. , Tang , M. , Zhu , X. , Ponaya , A. , Xiao , G. , Li , J. , Li , C. , Song , M. Y. , Cumagun , C. J. R. , Deng , Q. , Lu , G. , Jeon , J. S. , Naqvi , N. I. , & Zhou , B . ( 2015 ). Comparative genomics identifies the Magnaporthe oryzae avirulence effector AvrPi9 that triggers Pi9 -mediated blast resistance in rice . New Phytologist , 206 ( 4 ), 1463 – 1475 . doi: 10.1111/NPH.13310 OpenUrl CrossRef PubMed ↵ Yasuda , N. , Tsujimoto Noguchi , M. , & Fujita , Y . ( 2006 ). Partial mapping of avirulence genes AVR-Pii and AVR-Pia in the rice blast fungus Magnaporthe oryzae . Canadian Journal of Plant Pathology , 28 ( 3 ), 494 – 498 . doi: 10.1080/07060660609507325 OpenUrl CrossRef ↵ Yoshida , K. , Saitoh , H. , Fujisawa , S. , Kanzaki , H. , Matsumura , H. , Yoshida , K. , Tosa , Y. , Chuma , I. , Takano , Y. , Win , J. , Kamoun , S. , & Terauchi , R . ( 2009 ). Association Genetics Reveals Three Novel Avirulence Genes from the Rice Blast Fungal Pathogen Magnaporthe oryzae . The Plant Cell , 21 ( 5 ), 1573 – 1591 . doi: 10.1105/TPC.109.066324 OpenUrl Abstract / FREE Full Text ↵ Yoshida , K. , Saunders , D. G. O. , Mitsuoka , C. , Natsume , S. , Kosugi , S. , Saitoh , H. , Inoue , Y. , Chuma , I. , Tosa , Y. , Cano , L. M. , Kamoun , S. , & Terauchi , R . ( 2016 ). Host specialization of the blast fungus Magnaporthe oryzae is associated with dynamic gain and loss of genes linked to transposable elements . BMC Genomics , 17 ( 1 ), 1 – 18 . doi: 10.1186/S12864-016-2690-6/FIGURES/7 OpenUrl CrossRef PubMed ↵ Zhang , S. , Wang , L. , Wu , W. , He , L. , Yang , X. , & Pan , Q. ( 2015 ). Function and evolution of Magnaporthe oryzae avirulence gene AvrPib responding to the rice blast resistance gene Pib . Scientific Reports , 5 . doi: 10.1038/SREP11642 OpenUrl CrossRef ↵ Zhong , Z. , Chen , M. , Lin , L. , Han , Y. , Bao , J. , Tang , W. , Lin , L. , Lin , Y. , Somai , R. , Lu , L. , Zhang , W. , Chen , J. , Hong , Y. , Chen , X. , Wang , B. , Shen , W. C. , Lu , G. , Norvienyeku , J. , Ebbole , D. J. , & Wang , Z . ( 2018 ). Population genomic analysis of the rice blast fungus reveals specific events associated with expansion of three main clades . The ISME Journal , 12 ( 8 ), 1867 – 1878 . doi: 10.1038/S41396-018-0100-6 OpenUrl CrossRef PubMed ↵ Zhou , E. , Jia , Y. , Singh , P. , Correll , J. C. , & Lee , F. N . ( 2007 ). Instability of the Magnaporthe oryzae avirulence gene AVR-Pita alters virulence . Fungal Genetics and Biology: FG & B , 44 ( 10 ), 1024 – 1034 . doi: 10.1016/J.FGB.2007.02.003 OpenUrl CrossRef View the discussion thread. Back to top Previous Next Posted May 23, 2025. 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 Transposable elements create distinct genomic niches for effector evolution among Magnaporthe oryzae lineages 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 Transposable elements create distinct genomic niches for effector evolution among Magnaporthe oryzae lineages Ana Margarida Sampaio , Daniel Croll bioRxiv 2025.05.23.655716; doi: https://doi.org/10.1101/2025.05.23.655716 Share This Article: Copy Citation Tools Transposable elements create distinct genomic niches for effector evolution among Magnaporthe oryzae lineages Ana Margarida Sampaio , Daniel Croll bioRxiv 2025.05.23.655716; doi: https://doi.org/10.1101/2025.05.23.655716 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 Microbiology Subject Areas All Articles Animal Behavior and Cognition (7622) Biochemistry (17648) Bioengineering (13870) Bioinformatics (41880) Biophysics (21423) Cancer Biology (18553) Cell Biology (25458) Clinical Trials (138) Developmental Biology (13364) Ecology (19866) Epidemiology (2067) Evolutionary Biology (24290) Genetics (15589) Genomics (22475) Immunology (17711) Microbiology (40327) Molecular Biology (17145) Neuroscience (88472) Paleontology (666) Pathology (2826) Pharmacology and Toxicology (4815) Physiology (7635) Plant Biology (15114) Scientific Communication and Education (2044) Synthetic Biology (4286) Systems Biology (9815) Zoology (2268)","source_license":"CC-BY-4.0","license_restricted":false}