Search

Proteases are key components of the hydrolytic enzyme arsenal employed by fungal pathogens to invade their host plants. The recent advances in -omics era have facilitated identification of functional proteases involved in plant-fungus interactions. By comparison of the publically available sequences of fungal genomes we found that the number of protease genes present in the genome of Cladosporium fulvum, a biotrophic tomato pathogen, is comparable with that of hemibiotrophs. To identify host plant inducible protease genes and their products, we performed transcriptome and proteome analyses of C. fulvumin vitro and in planta by means of RNA-Seq/RT-qPCR and mass spectrometry. Transcriptome data showed that 14 out of the 59 predicted proteases are expressed during in vitro and in planta growth of C. fulvum, of which nine belong to serine proteases S8 and S10 and the rest belong to metallo- and aspartic proteases. Mass spectrometry confirmed the presence of six proteases at proteome level during plant infection. Expression of limited number of proteases by C. fulvum might sustain biotrophic growth and benefits its stealth pathogenesis.

Gasparite-(La), La(AsO4), is a new mineral (IMA 2018-079) from Mn ores of the Ushkatyn-III deposit, Central Kazakhstan (type locality) and from alpine fissures in metamorphic rocks of the Wanni glacier, Binn Valley, Switzerland (co-type locality). Gasparite-(La) is named for its dominant lanthanide, according to current nomenclature of rare-earth minerals. The occurrences and parageneses in both localities are distinct: minute isometric grains up to 15 μm in size, associated with friedelite, jacobsite, pennantite, manganhumite series minerals (alleghanyite, sonolite), sarkinite, tilasite, and retzian-(La) are typically embedded into calcite-rhodochrosite veinlets (Ushkatyn-III deposit) vs. elongated crystals up to 2 mm in size in classical alpine fissures in two-mica gneiss without indicative associated minerals (Wanni glacier). Their chemical compositions have been studied by EDX and WDX; crystal-chemical formulas of gasparite-(La) from the Ushkatyn-III deposit (holotype specimen) and Wanni glacier (co-type specimen) are (La0.65Ce0.17Nd0.07Ca0.06Mn0.05Pr0.02)1.02[(As0.70V0.28P0.02)1.00O4] and (La0.59Ce0.37Nd0.02 Ca0.02Th0.01)1.01[(As0.81P0.16Si0.02S0.02)1.01O4], respectively. In polished sections, crystals are yellow and translucent with bright submetallic luster. Selected reflectance values R1/R2 (λ, nm) for the holotype specimen in air are: 11.19/9.05 (400), 11.45/9.44 (500), 10.85/8.81 (600), 11.23/9.08 (700). The structural characteristics of gasparite-(La) were studied by means of EBSD (holotype specimen), XRD, and SREF (co-type specimen). Gasparite-(La) has a monoclinic structure with the space group P21/n. Our studies revealed that gasparite-(La) from the Ushkatyn-III deposit and Wanni glacier have different origins. La/Ce and As/P/V ratios in gasparite-(La) may be used as an indicator of formation conditions.

Genetic gain in potato is hampered by the heterozygous tetraploid genome of cultivated potato. Converting potato into a diploid inbred-line based F1-hybrid crop provides a promising route towards increased genetic gain. The introduction of a dominant S-locus inhibitor (Sli) gene into diploid potato germplasm allows efficient generation of self-fertilized seeds and thus the development of potato inbred lines. Little is known about the structure and function of the Sli locus. Here we describe the mapping of Sli to a 12.6 kb interval on chromosome 12 using a recombinant screen approach. One of two candidate genes present in this interval shows a unique sequence that is exclusively present in self-compatible lines. We describe an expression vector that converts self-incompatible genotypes into self-compatible and a CRISPR-Cas9 vector that converts SC genotypes into SI. The Sli gene encodes an F-box protein that is specifically expressed in pollen from self-compatible plants. A 533 bp insertion in the promotor of that gene leads to a gain of function mutation, which overcomes self-pollen rejection., The S-locus inhibitor (Sli) gene could allow potato breeding by facilitating production of diploid inbred lines. Here the authors show that Sli encodes an F-box protein with a promoter insertion enhancing expression in pollen can overcome pollen rejection in the styles of diploid potato.

Use of genomic prediction (GP) in tetraploid is becoming more common. Therefore, we think it is the right time for a comparison of GP models for tetraploid potato. GP models were compared that contrasted shrinkage with variable selection, parametric vs. non-parametric models and different ways of accounting for non-additive genetic effects. As a complement to GP, association studies were carried out in an attempt to understand the differences in prediction accuracy. We compared our GP models on a data set consisting of 147 cultivars, representing worldwide diversity, with over 39 k GBS markers and measurements on four tuber traits collected in six trials at three locations during 2 years. GP accuracies ranged from 0.32 for tuber count to 0.77 for dry matter content. For all traits, differences between GP models that utilised shrinkage penalties and those that performed variable selection were negligible. This was surprising for dry matter, as only a few additive markers explained over 50% of phenotypic variation. Accuracy for tuber count increased from 0.35 to 0.41, when dominance was included in the model. This result is supported by Genome Wide Association Study (GWAS) that found additive and dominance effects accounted for 37% of phenotypic variation, while significant additive effects alone accounted for 14%. For tuber weight, the Reproducing Kernel Hilbert Space (RKHS) model gave a larger improvement in prediction accuracy than explicitly modelling epistatic effects. This is an indication that capturing the between locus epistatic effects of tuber weight can be done more effectively using the semi-parametric RKHS model. Our results show good opportunities for GP in 4x potato.

With the rapid expansion of the application of genomics and sequencing in plant breeding, there is a constant drive for better reference genomes. In potato (Solanum tuberosum), the third largest food crop in the world, the related speciesS. phureja, designated “DM”, has been used as the most popular reference genome for the last 10 years. Here, we introduce thede novosequenced genome of Solyntus as the next standard reference in potato genome studies. A trueSolanum tuberosummade up of 116 contigs that is also highly homozygous, diploid, vigorous and self-compatible, Solyntus provides a more direct and contiguous reference then ever before available. It was constructed by sequencing with state-of-the-art long and short read technology and assembled with Canu. The 116 contigs were assembled into scaffolds to form each pseudochromosome, with three contigs to 17 contigs per chromosome. This assembly contains 93.7% of the single-copy gene orthologs from the Solanaceae set and has an N50 of 63.7 Mbp. The genome and related files can be found athttps://www.plantbreeding.wur.nl/Solyntus/. With the release of this research line and its draft genome we anticipate many exciting developments in (diploid) potato research.

Tomato leaf mould disease is caused by the biotrophic fungusCladosporium fulvum. During infection,C. fulvumproduces extracellular small secreted protein (SSP) effectors that function to promote colonization of the leaf apoplast. Resistance to the disease is governed byCfimmune receptor genes that encode receptor-like proteins (RLPs). These RLPs recognize specific SSP effectors to initiate a hypersensitive response (HR) that renders the pathogen avirulent.C. fulvumstrains capable of overcoming one or more of all clonedCfgenes have now emerged. To combat these strains, newCfgenes are required. An effectoromics approach was employed to identify wild tomato accessions carrying newCfgenes. Proteomics and transcriptome sequencing were first used to identify 70 apoplasticin planta-inducedC. fulvumSSPs. Based on sequence homology, 61 of these SSPs were novel or lacked known functional domains. Seven, however, had predicted structural homology to antimicrobial proteins, suggesting a possible role in mediating antagonistic microbe−microbe interactionsin planta. Wild tomato accessions were then screened for HR-associated recognition of 41 SSPs using thePotato virus X-based transient expression system. Nine SSPs were recognized by one or more accessions, suggesting that these plants carry newCfgenes available for incorporation into cultivated tomato.

The Cf-5 gene of tomato confers resistance to strains of the fungal pathogen Cladosporium fulvum carrying the avirulence gene Avr5. Although Cf-5 has been cloned, Avr5 has remained elusive. We report the cloning of Avr5 using a combined bioinformatic and transcriptome sequencing approach. RNA-Seq was performed on the sequenced race 0 strain (0WU; carrying Avr5), as well as a race 5 strain (IPO 1979; lacking a functional Avr5 gene) during infection of susceptible tomato. Forty-four in planta–induced C. fulvum candidate effector (CfCE) genes of 0WU were identified that putatively encode a secreted, small cysteine-rich protein. An expressed transcript sequence comparison between strains revealed two polymorphic CfCE genes in IPO 1979. One of these conferred avirulence to IPO 1979 on Cf-5 tomato following complementation with the corresponding 0WU allele, confirming identification of Avr5. Complementation also led to increased fungal biomass during infection of susceptible tomato, signifying a role for Avr5 in virulence. Seven of eight race 5 strains investigated escape Cf-5-mediated resistance through deletion of the Avr5 gene. Avr5 is heavily flanked by repetitive elements, suggesting that repeat instability, in combination with Cf-5-mediated selection pressure, has led to the emergence of race 5 strains deleted for the Avr5 gene.

Fungal Wor1-like proteins are conserved transcriptional regulators that are reported to regulate the virulence of several plant pathogenic fungi by affecting the expression of virulence genes. Here, we report the functional analysis of CfWor1, the homologue of Wor1 in Cladosporium fulvum. ?cfwor1 mutants produce sclerotium-like structures and rough hyphae, which are covered with a black extracellular matrix. These mutants do not sporulate and are no longer virulent on tomato. A CE.CfWor1 transformant that constitutively expresses CfWor1 produces fewer spores with altered morphology and is also reduced in virulence. RNA-seq and RT-qrtPCR analyses suggest that reduced virulence of ?cfwor1 mutants is due to global downregulation of transcription, translation and mitochondrial respiratory chain. The reduced virulence of the CE.CfWor1 transformant is likely due to downregulation of effector genes. Complementation of a non-virulent ?fosge1 (Wor1-homologue) mutant of Fusarium oxysporum f. sp. lycopersici with CfWor1 restored expression of the SIX effector genes in this fungus, but not its virulence. Chimeric proteins of CfWor1/FoSge1 also only partially restored defects of the ?fosge1 mutant, suggesting that these transcriptional regulators have functionally diverged. Altogether, our results suggest that CfWor1 primarily regulates development of C.?fulvum, which indirectly affects the expression of a subset of virulence genes.

Pseudogenes are genes with significant homology to functional genes but contain disruptive mutations (DMs) leading to production of non- or partially functional proteins. Little is known about pseudogenization in pathogenic fungi with different lifestyles. Here we report on identification of DMs causing pseudogenes in the genomes of the fungal plant pathogens Botrytis cinerea, Cladosporium fulvum, Dothistroma septosporum, Mycosphaerella fijiensis, Verticillium dahliae and Zymoseptoria tritici. In these fungi we have identified 1740 gene models containing 2795 DMs obtained by an alignment-based gene prediction method. The contribution of sequencing errors to DMs was minimized by analyses of resequenced genomes to obtain a refined data set of 924 gene models containing 1666 true DMs. The frequency of pseudogenes varied from 1 to 5% in the gene catalogues of these fungi, being the highest in the asexually reproducing fungi C. fulvum (4.9%), followed by D. septosporum (2.4%) and V. dahliae (2.1%). The majority of pseudogenes does not represent recent gene duplications, but members of multi-gene families and unitary genes. In general there was no bias for pseudogenization of specific genes in the six fungi. Single exceptions are those encoding secreted proteins including proteases which appeared more frequently pseudogenized in C. fulvum than in D. septosporum. Most pseudogenes present in these two phylogenically closely related fungi are not shared suggesting that they are related to adaptation to a different host (tomato versus pine) and lifestyle (biotroph versus hemi-biotroph)

Spliceosomal introns are noncoding sequences that separate exons in eukaryotic genes and are removed from pre-messenger RNAs by the splicing machinery. Their origin has remained a mystery in biology since their discovery [ [1] and [2]] because intron gains seem to be infrequent in many eukaryotic lineages [ [3] and [4]]. Although a few recent intron gains have been reported [ [5] and [6]], none of the proposed gain mechanisms [7] can convincingly explain the high number of introns in present-day eukaryotic genomes. Here we report on particular spliceosomal introns that share high sequence similarity and are reminiscent of introner elements [8]. These elements multiplied in unrelated genes of six fungal genomes and account for the vast majority of intron gains in these fungal species. Such introner-like elements (ILEs) contain all typical characteristics of regular spliceosomal introns (RSIs) [ [9] and [10]] but are longer and predicted to harbor more stable secondary structures. However, dating of multiplication events showed that they degenerate in sequence and length within 100,000 years to eventually become indistinguishable from RSIs. We suggest that ILEs not only account for intron gains in six fungi but also in ancestral eukaryotes to give rise to most RSIs by a yet unknown multiplication mechanism

International audience; Spliceosomal introns are key components of the eukaryotic gene structure. Although they contributed to the emergence of eukaryotes, their origin remains elusive. In fungi, they might originate from the multiplication of invasive introns named Introner-Like Elements (ILEs). However, so far ILEs have been observed in six fungal species only, including Fulvia fulva and Dothistroma septosporum (Dothideomycetes), arguing against ILE insertion as a general mechanism for intron gain. Here, we identified novel ILEs in eight additional fungal species that are phylogenetically related to F. fulva and D. septosporum using PCR amplification with primers derived from previously identified ILEs. The ILE content appeared unique to each species, suggesting independent multiplication events. Interestingly, we identified four genes each containing two gained ILEs. By analysing intron positions in orthologues of these four genes in Ascomycota, we found that three ILEs had inserted within a 15 bp window that contains regular spliceosomal introns in other fungal species. These three positions are not the result of intron sliding because ILEs are newly gained introns. Furthermore, the alternative hypothesis of an inferred ancestral gain followed by independent losses contradicts the observed degeneration of ILEs. These observations clearly indicate three parallel intron gains in four genes that were randomly identified. Our findings suggest that parallel intron gain is a phenomenon that has been highly underestimated in ILE-containing fungi, and likely in the whole fungal kingdom.

Fungal Wor1-like proteins are conserved transcriptional regulators that are reported to regulate the virulence of several plant pathogenic fungi by affecting the expression of virulence genes. Here, we report the functional analysis of CfWor1, the homologue of Wor1 in Cladosporium fulvum. Δcfwor1 mutants produce sclerotium-like structures and rough hyphae, which are covered with a black extracellular matrix. These mutants do not sporulate and are no longer virulent on tomato. A CE.CfWor1 transformant that constitutively expresses CfWor1 produces fewer spores with altered morphology and is also reduced in virulence. RNA-seq and RT-qrtPCR analyses suggest that reduced virulence of Δcfwor1 mutants is due to global downregulation of transcription, translation and mitochondrial respiratory chain. The reduced virulence of the CE.CfWor1 transformant is likely due to downregulation of effector genes. Complementation of a non-virulent Δfosge1 (Wor1-homologue) mutant of Fusarium oxysporum f. sp. lycopersici with CfWor1 restored expression of the SIX effector genes in this fungus, but not its virulence. Chimeric proteins of CfWor1/FoSge1 also only partially restored defects of the Δfosge1 mutant, suggesting that these transcriptional regulators have functionally diverged. Altogether, our results suggest that CfWor1 primarily regulates development of C. fulvum, which indirectly affects the expression of a subset of virulence genes.

Pseudogenes are genes with significant homology to functional genes, but contain disruptive mutations (DMs) leading to the production of non‐ or partially functional proteins. Little is known about pseudogenization in pathogenic fungi with different lifestyles. Here, we report the identification of DMs causing pseudogenes in the genomes of the fungal plant pathogens Botrytis cinerea, Cladosporium fulvum, Dothistroma septosporum, Mycosphaerella fijiensis, Verticillium dahliae and Zymoseptoria tritici. In these fungi, we identified 1740 gene models containing 2795 DMs obtained by an alignment‐based gene prediction method. The contribution of sequencing errors to DMs was minimized by analyses of resequenced genomes to obtain a refined dataset of 924 gene models containing 1666 true DMs. The frequency of pseudogenes varied from 1% to 5% in the gene catalogues of these fungi, being the highest in the asexually reproducing fungus C. fulvum (4.9%), followed by D. septosporum (2.4%) and V. dahliae (2.1%). The majority of pseudogenes do not represent recent gene duplications, but members of multi‐gene families and unitary genes. In general, there was no bias for pseudogenization of specific genes in the six fungi. Single exceptions were those encoding secreted proteins, including proteases, which appeared more frequently pseudogenized in C. fulvum than in D. septosporum. Most pseudogenes present in these two phylogenetically closely related fungi are not shared, suggesting that they are related to adaptation to a different host (tomato versus pine) and lifestyle (biotroph versus hemibiotroph).

Background Automated gene-calling is still an error-prone process, particularly for the highly plastic genomes of fungal species. Improvement through quality control and manual curation of gene models is a time-consuming process that requires skilled biologists and is only marginally performed. The wealth of available fungal genomes has not yet been exploited by an automated method that applies quality control of gene models in order to obtain more accurate genome annotations. Results We provide a novel method named alignment-based fungal gene prediction (ABFGP) that is particularly suitable for plastic genomes like those of fungi. It can assess gene models on a gene-by-gene basis making use of informant gene loci. Its performance was benchmarked on 6,965 gene models confirmed by full-length unigenes from ten different fungi. 79.4% of all gene models were correctly predicted by ABFGP. It improves the output of ab initio gene prediction software due to a higher sensitivity and precision for all gene model components. Applicability of the method was shown by revisiting the annotations of six different fungi, using gene loci from up to 29 fungal genomes as informants. Between 7,231 and 8,337 genes were assessed by ABFGP and for each genome between 1,724 and 3,505 gene model revisions were proposed. The reliability of the proposed gene models is assessed by an a posteriori introspection procedure of each intron and exon in the multiple gene model alignment. The total number and type of proposed gene model revisions in the six fungal genomes is correlated to the quality of the genome assembly, and to sequencing strategies used in the sequencing centre, highlighting different types of errors in different annotation pipelines. The ABFGP method is particularly successful in discovering sequence errors and/or disruptive mutations causing truncated and erroneous gene models. Conclusions The ABFGP method is an accurate and fully automated quality control method for fungal gene catalogues that can be easily implemented into existing annotation pipelines. With the exponential release of new genomes, the ABFGP method will help decreasing the number of gene models that require additional manual curation.

We sequenced and compared the genomes of the Dothideomycete fungal plant pathogens Cladosporium fulvum (Cfu) (syn. Passalora fulva) and Dothistroma septosporum (Dse) that are closely related phylogenetically, but have different lifestyles and hosts. Although both fungi grow extracellularly in close contact with host mesophyll cells, Cfu is a biotroph infecting tomato, while Dse is a hemibiotroph infecting pine. The genomes of these fungi have a similar set of genes (70% of gene content in both genomes are homologs), but differ significantly in size (Cfu >61.1-Mb; Dse 31.2-Mb), which is mainly due to the difference in repeat content (47.2% in Cfu versus 3.2% in Dse). Recent adaptation to different lifestyles and hosts is suggested by diverged sets of genes. Cfu contains an α-tomatinase gene that we predict might be required for detoxification of tomatine, while this gene is absent in Dse. Many genes encoding secreted proteins are unique to each species and the repeat-rich areas in Cfu are enriched for these species-specific genes. In contrast, conserved genes suggest common host ancestry. Homologs of Cfu effector genes, including Ecp2 and Avr4, are present in Dse and induce a Cf-Ecp2- and Cf-4-mediated hypersensitive response, respectively. Strikingly, genes involved in production of the toxin dothistromin, a likely virulence factor for Dse, are conserved in Cfu, but their expression differs markedly with essentially no expression by Cfu in planta. Likewise, Cfu has a carbohydrate-degrading enzyme catalog that is more similar to that of necrotrophs or hemibiotrophs and a larger pectinolytic gene arsenal than Dse, but many of these genes are not expressed in planta or are pseudogenized. Overall, comparison of their genomes suggests that these closely related plant pathogens had a common ancestral host but since adapted to different hosts and lifestyles by a combination of differentiated gene content, pseudogenization, and gene regulation., Author Summary We compared the genomes of two closely related pathogens with very different lifestyles and hosts: C. fulvum (Cfu), a biotroph of tomato, and D. septosporum (Dse), a hemibiotroph of pine. Some differences in gene content were identified that can be directly related to their different hosts, such as the presence of a gene involved in degradation of a tomato saponin only in Cfu. However, in general the two species share a surprisingly large proportion of genes. Dse has functional homologs of Cfu effector genes, while Cfu has genes for biosynthesis of dothistromin, a toxin probably associated with virulence in Dse. Cfu also has an unexpectedly large content of genes for biosynthesis of other secondary metabolites and degradation of plant cell walls compared to Dse, contrasting with its host preference and lifestyle. However, many of these genes were not expressed in planta or were pseudogenized. These results suggest that evolving species may retain genetic signatures of the host and lifestyle preferences of their ancestor and that evolution of new genes, gene regulation, and pseudogenization are important factors in adaptation.

Background MicroRNAs (miRNAs), short ~21-nucleotide RNA molecules, play an important role in post-transcriptional regulation of gene expression. The number of known miRNA hairpins registered in the miRBase database is rapidly increasing, but recent reports suggest that many miRNAs with restricted temporal or tissue-specific expression remain undiscovered. Various strategies for in silico miRNA identification have been proposed to facilitate miRNA discovery. Notably support vector machine (SVM) methods have recently gained popularity. However, a drawback of these methods is that they do not provide insight into the biological properties of miRNA sequences. Results We here propose a new strategy for miRNA hairpin prediction in which the likelihood that a genomic hairpin is a true miRNA hairpin is evaluated based on statistical distributions of observed biological variation of properties (descriptors) of known miRNA hairpins. These distributions are transformed into a single and continuous outcome classifier called the L score. Using a dataset of known miRNA hairpins from the miRBase database and an exhaustive set of genomic hairpins identified in the genome of Caenorhabditis elegans, a subset of 18 most informative descriptors was selected after detailed analysis of correlation among and discriminative power of individual descriptors. We show that the majority of previously identified miRNA hairpins have high L scores, that the method outperforms miRNA prediction by threshold filtering and that it is more transparent than SVM classifiers. Conclusion The L score is applicable as a prediction classifier with high sensitivity for novel miRNA hairpins. The L- score approach can be used to rank and select interesting miRNA hairpin candidates for downstream experimental analysis when coupled to a genome-wide set of in silico-identified hairpins or to facilitate the analysis of large sets of putative miRNA hairpin loci obtained in deep-sequencing efforts of small RNAs. Moreover, the in-depth analyses of miRNA hairpins descriptors preceding and determining the L score outcome could be used as an extension to miRBase entries to help increase the reliability and biological relevance of the miRNA registry.

The Gene Ontology (GO) is a widely used controlled vocabulary for the description of gene function. In this study we quantify the usage of multiple and hierarchically independent GO terms in the curated genome annotations of seven well-studied species. In most genomes, significant proportions (6-60%) of genes have been annotated with multiple and hierarchically independent terms. This may be necessary to attain adequate specificity of description. One noticeable exception is Arabidopsis thaliana, in which genes are much less frequently annotated with multiple terms (6-14%). In contrast, an analysis of the occurrence of InterPro hits in the proteomes of the seven species, followed by a mapping of the hits to GO terms, did not reveal an aberrant pattern for the A. thaliana genome. This study shows the widespread usage of multiple hierarchically independent GO terms in the functional annotation of genes. By consequence, probabilistic methods that aim to predict gene function automatically through integration of diverse genomic datasets, and that employ the GO, must be able to predict such multiple terms. We attribute the low frequency with which multiple GO terms are used in Arabidopsis to deviating practices in the genome annotation and curation process between communities of annotators. This may bias genome-scale comparisons of gene function between different species. GO term assignment should therefore be performed according to strictly similar rules and standards.

Background Modern omics research involves the application of high-throughput technologies that generate vast volumes of data. These data need to be pre-processed, analyzed and integrated with existing knowledge through the use of diverse sets of software tools, models and databases. The analyses are often interdependent and chained together to form complex workflows or pipelines. Given the volume of the data used and the multitude of computational resources available, specialized pipeline software is required to make high-throughput analysis of large-scale omics datasets feasible. Results We have developed a generic pipeline system called Cyrille2. The system is modular in design and consists of three functionally distinct parts: 1) a web based, graphical user interface (GUI) that enables a pipeline operator to manage the system; 2) the Scheduler, which forms the functional core of the system and which tracks what data enters the system and determines what jobs must be scheduled for execution, and; 3) the Executor, which searches for scheduled jobs and executes these on a compute cluster. Conclusion The Cyrille2 system is an extensible, modular system, implementing the stated requirements. Cyrille2 enables easy creation and execution of high throughput, flexible bioinformatics pipelines.

The recent discovery of introner-like elements (ILEs) in six fungal species shed new light on the origin of regular spliceosomal introns (RSIs) and the mechanism of intron gains. These novel spliceosomal introns are found in hundreds of copies, are longer than RSIs and harbor stable predicted secondary structures. Yet, they are prone to degeneration in sequence and length to become undistinguishable from RSIs, suggesting that ILEs are predecessors of most RSIs. In most fungi, other near-identical introns were found duplicated in lower numbers in the same gene or in unrelated genes, indicating that intron duplication is a widespread phenomenon. However, ILEs are associated with the majority of intron gains, suggesting that the other types of duplication are of minor importance to the overall gains of introns. Our data support the hypothesis that ILEs’ multiplication corresponds to the main mechanism of intron gain in fungi.