Phylotranscriptomics of the Pentapetalae Reveals Frequent Regulatory Variation in Plant Local Responses to the Fungal Pathogen Sclerotinia sclerotiorum

Time-course dual RNA-seq analyses and gene identification during early stages of plant-Phytophthora infestans interactions

Late blight caused by Phytophthora infestans is a major threat to global potato and tomato production. Sustainable management of late blight requires the development of resistant crop cultivars. This process can be facilitated by high-throughput identification of functional genes involved in late blight pathogenesis. In this study, we generated a high-quality transcriptomic time-course dataset focusing on the initial 24 h of contact between P. infestans and 3 solanaceous plant species, tobacco (Nicotiana benthamiana), tomato (Solanum lycopersicum), and potato (Solanum tuberosum). Our results demonstrate species-specific transcriptional regulation in early stages of the infection. Transient silencing of putative RIBOSE-5-PHOSPHATE ISOMERASE and HMG-CoA REDUCTASE genes in N. benthamiana lowered plant resistance against P. infestans. Furthermore, heterologous expression of a putative tomato Golgi-localized nucleosugar transporter-encoding gene exacerbated P. infestans infection of N. benthamiana. In comparison, bioassays using transgenic tomato lines showed that the quantitative disease resistance genes were required but insufficient for late blight resistance; genetic knock-out of the susceptibility gene enhanced resistance. The same RNA-seq dataset was exploited to examine the transcriptional landscape of P. infestans and revealed host-specific gene expression patterns in the pathogen. This temporal transcriptomic diversity, in combination with genomic distribution features, identified the P. infestans IPI-B family GLYCINE-RICH PROTEINs as putative virulence factors that promoted disease severity or induced plant tissue necrosis when transiently expressed in N. benthamiana. These functional genes underline the effectiveness of functional gene-mining through a time-course dual RNA-seq approach and provide insight into the molecular interactions between solanaceous plants and P. infestans.

The Marchantia polymorpha pangenome reveals ancient mechanisms of plant adaptation to the environment

Plant adaptation to terrestrial life started 450 million years ago and has played a major role in the evolution of life on Earth. The genetic mechanisms allowing this adaptation to a diversity of terrestrial constraints have been mostly studied by focusing on flowering plants. Here, we gathered a collection of 133 accessions of the model bryophyte Marchantia polymorpha and studied its intraspecific diversity using selection signature analyses, a genome-environment association study and a pangenome. We identified adaptive features, such as peroxidases or nucleotide-binding and leucine-rich repeats (NLRs), also observed in flowering plants, likely inherited from the first land plants. The M. polymorpha pangenome also harbors lineage-specific accessory genes absent from seed plants. We conclude that different land plant lineages still share many elements from the genetic toolkit evolved by their most recent common ancestor to adapt to the terrestrial habitat, refined by lineage-specific polymorphisms and gene family evolution.

Phylotranscriptomics provides a treasure trove of flood-tolerance mechanisms in the Cardamineae tribe

Flooding events are highly detrimental to most terrestrial plant species. However, there is an impressive diversity of plant species that thrive in flood-prone regions and represent a treasure trove of unexplored flood-resilience mechanisms. Here we surveyed a panel of four species from the Cardamineae tribe representing a broad tolerance range. This included the flood-tolerant Cardamine pratensis, Rorippa sylvestris and Rorippa palustris and the flood-sensitive species Cardamine hirsuta. All four species displayed a quiescent strategy, evidenced by the repression of shoot growth underwater. Comparative transcriptomics analyses between the four species and the sensitive model species Arabidopsis thaliana were facilitated via de novo transcriptome assembly and identification of 16 902 universal orthogroups at a high resolution. Our results suggest that tolerance likely evolved separately in the Cardamine and Rorippa species. While the Rorippa response was marked by a strong downregulation of cell-cycle genes, Cardamine minimized overall transcriptional regulation. However, a weak starvation response was a universal trait of tolerant species, potentially achieved in multiple ways. It could result from a strong decline in cell-cycle activity, but is also intertwined with autophagy, senescence, day-time photosynthesis and night-time fermentation capacity. Our data set provides a rich source to study adaptational mechanisms of flooding tolerance.

High-Resolution Disease Phenotyping Reveals Distinct Resistance Mechanisms of Tomato Crop Wild Relatives against Sclerotinia sclerotiorum

Besides the well-understood qualitative disease resistance, plants possess a more complex quantitative form of resistance: quantitative disease resistance (QDR). QDR is commonly defined as a partial but more durable form of resistance and, therefore, might display a valuable target for resistance breeding. The characterization of QDR phenotypes, especially of wild crop relatives, displays a bottleneck in deciphering QDR's genomic and regulatory background. Moreover, the relationship between QDR parameters, such as infection frequency, lag-phase duration, and lesion growth rate, remains elusive. High hurdles for applying modern phenotyping technology, such as the low availability of phenotyping facilities or complex data analysis, further dampen progress in understanding QDR. Here, we applied a low-cost (<1.000 €) phenotyping system to measure lesion growth dynamics of wild tomato species (e.g., Solanum pennellii or Solanum pimpinellifolium). We provide insight into QDR diversity of wild populations and derive specific QDR mechanisms and their cross-talk. We show how temporally continuous observations are required to dissect end-point severity into functional resistance mechanisms. The results of our study show how QDR can be maintained by facilitating different defense mechanisms during host-parasite interaction and that the capacity of the QDR toolbox highly depends on the host's genetic context. We anticipate that the present findings display a valuable resource for more targeted functional characterization of the processes involved in QDR. Moreover, we show how modest phenotyping technology can be leveraged to help answer highly relevant biological questions.

Till death do us pair: Co-evolution of plant-necrotroph interactions

Plants use programmed cell death as a potent defense response against biotrophic pathogens that require living host cells to thrive. However, cell death can promote infection by necrotrophic pathogens. This discrepancy creates specific co-evolutionary dynamics in the interaction between plants and necrotrophs. Necrotrophic pathogens produce diverse cell death-inducing effectors that act redundantly on several plant targets and sometimes suppress plant immune responses as an additional function. Plants use surface receptors that recognize necrotrophic effectors to increase quantitative disease resistance, some of which evolved independently in several plant lineages. Co-evolution has shaped molecular mechanisms involved in plant-necrotroph interactions into robust systems, relying on degenerate and multifunctional modules, general-purpose components, and compartmentalized functioning.

Comparative transcriptome profiling reveals the importance of GmSWEET15 in soybean susceptibility to Sclerotinia sclerotiorum

Soybean sclerotinia stem rot (SSR) is a disease caused by Sclerotinia sclerotiorum that causes incalculable losses in soybean yield each year. Considering the lack of effective resistance resources and the elusive resistance mechanisms, we are urged to develop resistance genes and explore their molecular mechanisms. Here, we found that loss of GmSWEET15 enhanced the resistance to S. sclerotiorum, and we explored the molecular mechanisms by which gmsweet15 mutant exhibit enhanced resistance to S. sclerotiorum by comparing transcriptome. At the early stage of inoculation, the wild type (WT) showed moderate defense response, whereas gmsweet15 mutant exhibited more extensive and intense transcription reprogramming. The gmsweet15 mutant enriched more biological processes, including the secretory pathway and tetrapyrrole metabolism, and it showed stronger changes in defense response, protein ubiquitination, MAPK signaling pathway-plant, plant-pathogen interaction, phenylpropanoid biosynthesis, and photosynthesis. The more intense and abundant transcriptional reprogramming of gmsweet15 mutant may explain how it effectively delayed colonization by S. sclerotiorum. In addition, we identified common and specific differentially expressed genes between WT and gmsweet15 mutant after inoculation with S. sclerotiorum, and gene sets and genes related to gmsweet15_24 h were identified through Gene Set Enrichment Analysis. Moreover, we constructed the protein-protein interaction network and gene co-expression networks and identified several groups of regulatory networks of gmsweet15 mutant in response to S. sclerotiorum, which will be helpful for the discovery of candidate functional genes. Taken together, our results elucidate molecular mechanisms of delayed colonization by S. sclerotiorum after loss of GmSWEET15 in soybean, and we propose novel resources for improving resistance to SSR.

Pathogen-derived mechanical cues potentiate the spatio-temporal implementation of plant defense

BACKGROUND

The ongoing adaptation of plants to their environment is the basis for their survival. In this adaptation, mechanoperception of gravity and local curvature plays a role of prime importance in finely regulating growth and ensuring a dynamic balance preventing buckling. However, the abiotic environment is not the exclusive cause of mechanical stimuli. Biotic interactions between plants and microorganisms also involve physical forces and potentially mechanoperception. Whether pathogens trigger mechanoperception in plants and the impact of mechanotransduction on the regulation of plant defense remains however elusive.

RESULTS

Here, we found that the perception of pathogen-derived mechanical cues by microtubules potentiates the spatio-temporal implementation of plant immunity to fungus. By combining biomechanics modeling and image analysis of the post-invasion stage, we reveal that fungal colonization releases plant cell wall-born tension locally, causing fluctuations of tensile stress in walls of healthy cells distant from the infection site. In healthy cells, the pathogen-derived mechanical cues guide the reorganization of mechanosensing cortical microtubules (CMT). The anisotropic patterning of CMTs is required for the regulation of immunity-related genes in distal cells. The CMT-mediated mechanotransduction of pathogen-derived cues increases Arabidopsis disease resistance by 40% when challenged with the fungus Sclerotinia sclerotiorum.

CONCLUSIONS

CMT anisotropic patterning triggered by pathogen-derived mechanical cues activates the implementation of early plant defense in cells distant from the infection site. We propose that the mechano-signaling triggered immunity (MTI) complements the molecular signals involved in pattern and effector-triggered immunity.

Genetic co-option into plant-filamentous pathogen interactions

Plants are engaged in a coevolutionary arms race with their pathogens that drives rapid diversification and specialization of genes involved in resistance and virulence. However, some major innovations in plant-pathogen interactions, such as molecular decoys, trans-kingdom RNA interference, two-speed genomes, and receptor networks, evolved through the expansion of the functional landscape of genes. This is a typical outcome of genetic co-option, the evolutionary process by which available genes are recruited into new biological functions. Co-option into plant-pathogen interactions emerges generally from (i) cis-regulatory variation, (ii) horizontal gene transfer (HGT), (iii) mutations altering molecular promiscuity, and (iv) rewiring of gene networks and protein complexes. Understanding these molecular mechanisms is key for the functional and predictive biology of plant-pathogen interactions.

Phylogenetic analysis and expression profiles of jasmonate ZIM-domain gene family provide insight into abiotic stress resistance in sunflower

Jasmonate ZIM-domain (JAZ) proteins act as inhibitory factors of the jasmonic acid (JA) pathway, which is involved in regulating plant development and defense responses. However, there are no extensive studies available on JAZ genes in sunflower (Helianthus annuus L.). In this study, the phylogenetic analysis of 139 putative JAZ genes from eight plants demonstrated that these JAZs could be divided into five groups (Groups I-V), and the 27 sunflower JAZs (HaJAZs) were classified into these five groups. All groups contained genes from both monocotyledons and dicotyledons, indicating that the emergence of JAZ genes predates the differentiation of monocotyledons and dicotyledons. Both segmental and tandem duplications contributed greatly to this gene family's expansion in sunflower, especially in Group II. Moreover, the expression profiles of HaJAZ genes under normal conditions, hormone treatments or abiotic stresses were analyzed based on RNA-seq data. HaJAZ2 may be undergoing pseudogenization as a nonfunctional gene because it was not expressed in any tissue. Many HaJAZ genes in roots upregulated their expression when involved in responding to exogenous hormones, especially methyl-jasmonate. The abiotic stress treatments of sunflower showed that HaJAZ5, HaJAZ15, HaJAZ17, HaJAZ20, and HaJAZ21 tend to be sensitive to certain abiotic stresses. HaJAZs from different groups may share similar functions but also exercise their unique functions when responding to abiotic stresses. We speculated that this gene family was conserved in sequence but varied in its expression among duplicated HaJAZ genes, which implies that they may confer neofunctionalization in the adaptation to abiotic stresses; this work provides insight into the resistance of sunflowers and their adaptation to diverse environmental conditions.

Genome-Wide Identification, Evolution and Expression Analysis of the Glutathione S-Transferase Supergene Family in Euphorbiaceae

Euphorbiaceae, a family of plants mainly grown in the tropics and subtropics, is also widely distributed all over the world and is well known for being rich in rubber, oil, medicinal materials, starch, wood and other economically important plant products. Glutathione S-transferases (GSTs) constitute a family of proteins encoded by a large supergene family and are widely expressed in animals, bacteria, fungi and plants, but with few reports of them in Euphorbiaceae plants. These proteins participate in and regulate the detoxification and oxidative stress response of heterogeneous organisms, resistance to stress, growth and development, signal transduction and other related processes. In this study, we identified and analyzed the whole genomes of four species of Euphorbiaceae, namely Ricinus communis, Jatropha curcas, Hevea brasiliensis, and Manihot esculenta, which have high economic and practical value. A total of 244 GST genes were identified. Based on their sequence characteristics and conserved domain types, the GST supergene family in Euphorbiaceae was classified into 10 subfamilies. The GST supergene families of Euphorbiaceae and Arabidopsis have been found to be highly conserved in evolution, and tandem repeats and translocations in these genes have made the greatest contributions to gene amplification here and have experienced strong purification selection. An evolutionary analysis showed that Euphorbiaceae GST genes have also evolved into new subtribes (GSTO, EF1BG, MAPEG), which may play a specific role in Euphorbiaceae. An analysis of expression patterns of the GST supergene family in Euphorbiaceae revealed the functions of these GSTs in different tissues, including resistance to stress and participation in herbicide detoxification. In addition, an interaction analysis was performed to determine the GST gene regulatory mechanism. The results of this study have laid a foundation for further analysis of the functions of the GST supergene family in Euphorbiaceae, especially in stress and herbicide detoxification. The results have also provided new ideas for the study of the regulatory mechanism of the GST supergene family, and have provided a reference for follow-up genetics and breeding work.

The genome of the Paleogene relic tree Bretschneidera sinensis: insights into trade-offs in gene family evolution, demographic history, and adaptive SNPs

Among relic species, genomic information may provide the key to inferring their long-term survival. Therefore, in this study, we investigated the genome of the Paleogene relic tree species, Bretschneidera sinensis, which is a rare endemic species within southeastern Asia. Specifically, we assembled a high-quality genome for B. sinensis using PacBio high-fidelity and high-throughput chromosome conformation capture reads and annotated it with long and short RNA sequencing reads. Using the genome, we then detected a trade-off between active and passive disease defences among the gene families. Gene families involved in salicylic acid and MAPK signalling pathways expanded as active defence mechanisms against disease, but families involved in terpene synthase activity as passive defences contracted. When inferring the long evolutionary history of B. sinensis, we detected population declines corresponding to historical climate change around the Eocene–Oligocene transition and to climatic fluctuations in the Quaternary. Additionally, based on this genome, we identified 388 single nucleotide polymorphisms (SNPs) that were likely under selection, and showed diverse functions in growth and stress responses. Among them, we further found 41 climate-associated SNPs. The genome of B. sinensis and the SNP dataset will be important resources for understanding extinction/diversification processes using comparative genomics in different lineages.

SWEET Gene Family in Sugar Beet (Beta vulgaris): Genome-Wide Survey, Phylogeny and Expression Analysis

&lt;b&gt;Background and Objective:&lt;/b&gt; The SWEET (Sugars Will Eventually be Exported Transporter) proteins play important roles in modulating the growth and development processes in plants. However, little information is available on the SWEET family in sugar beet (&lt;i&gt;Beta vulgaris&lt;/i&gt;). The objectives of this present study were to genome-wide identify and characterize the BvSWEET family in sugar beet. &lt;b&gt;Materials and Methods:&lt;/b&gt; Based on the available genome, proteome and transcriptome databases of sugar beet, various computational tools have been used to analyze the nucleotide and full-length protein sequences of members of the BvSWEET family. &lt;b&gt;Results:&lt;/b&gt; A total of 16 members of the BvSWEET family has been identified in sugar beet at the genome-wide scale. Structural analysis indicated that the BvSWEET family exhibited variable characteristics. Furthermore, the BvSWEET family in sugar beet could be categorized into four distinct groups like in other plant species. Of our interest, we found that some &lt;i&gt;BvSWEET&lt;/i&gt; genes exhibited strongly preferential expression in major organs/tissues under adverse environmental stimuli. &lt;b&gt;Conclusion:&lt;/b&gt; The results provided a comprehensive foundation for further functional characterization of the &lt;i&gt;BvSWEET &lt;/i&gt;gene family.

Transcriptional response to host chemical cues underpins the expansion of host range in a fungal plant pathogen lineage

The host range of parasites is an important factor in assessing the dynamics of disease epidemics. The evolution of pathogens to accommodate new hosts may lead to host range expansion, a process the molecular bases of which are largely enigmatic. The fungus Sclerotinia sclerotiorum has been reported to parasitize more than 400 plant species from diverse eudicot families while its close relative, S. trifoliorum, is restricted to plants from the Fabaceae family. We analyzed S. sclerotiorum global transcriptome reprogramming on hosts from six botanical families and reveal a flexible, host-specific transcriptional program. We generated a chromosome-level genome assembly for S. trifoliorum and found near-complete gene space conservation in two representative strains of broad and narrow host range Sclerotinia species. However, S. trifoliorum showed increased sensitivity to the Brassicaceae defense compound camalexin. Comparative analyses revealed a lack of transcriptional response to camalexin in the S. trifoliorum strain and suggest that regulatory variation in detoxification and effector genes at the population level may associate with the genetic accommodation of Brassicaceae in the Sclerotinia host range. Our work proposes transcriptional plasticity and the co-existence of signatures for generalist and polyspecialist adaptive strategies in the genome of a plant pathogen.

PRGdb 4.0: an updated database dedicated to genes involved in plant disease resistance process

The Plant Resistance Genes database (PRGdb; http://prgdb.org/prgdb4/) has been greatly expanded, keeping pace with the increasing amount of available knowledge and data (sequenced proteomes, cloned genes, public analysis data, etc.). The easy-to-use style of the database website has been maintained, while an updated prediction tool, more data and a new section have been added. This new section will contain plant resistance transcriptomic experiments, providing additional easy-to-access experimental information. DRAGO3, the tool for automatic annotation and prediction of plant resistance genes behind PRGdb, has been improved in both accuracy and sensitivity, leading to more reliable predictions. PRGdb offers 199 reference resistance genes and 586.652 putative resistance genes from 182 sequenced proteomes. Compared to the previous release, PRGdb 4.0 has increased the number of reference resistance genes from 153 to 199, the number of putative resistance genes from 177K from 76 proteomes to 586K from 182 sequenced proteomes. A new section has been created that collects plant-pathogen transcriptomic data for five species of agricultural interest. Thereby, with these improvements and data expansions, PRGdb 4.0 aims to serve as a reference to the plant scientific community and breeders worldwide, helping to further study plant resistance mechanisms that contribute to fighting pathogens.

Status and advances in mining for blackleg (Leptosphaeria maculans) quantitative resistance (QR) in oilseed rape (Brassica napus)

KEY MESSAGE

Quantitative resistance (QR) loci discovered through genetic and genomic analyses are abundant in the Brassica napus genome, providing an opportunity for their utilization in enhancing blackleg resistance. Quantitative resistance (QR) has long been utilized to manage blackleg in Brassica napus (canola, oilseed rape), even before major resistance genes (R-genes) were extensively explored in breeding programmes. In contrast to R-gene-mediated qualitative resistance, QR reduces blackleg symptoms rather than completely eliminating the disease. As a polygenic trait, QR is controlled by numerous genes with modest effects, which exerts less pressure on the pathogen to evolve; hence, its effectiveness is more durable compared to R-gene-mediated resistance. Furthermore, combining QR with major R-genes has been shown to enhance resistance against diseases in important crops, including oilseed rape. For these reasons, there has been a renewed interest among breeders in utilizing QR in crop improvement. However, the mechanisms governing QR are largely unknown, limiting its deployment. Advances in genomics are facilitating the dissection of the genetic and molecular underpinnings of QR, resulting in the discovery of several loci and genes that can be potentially deployed to enhance blackleg resistance. Here, we summarize the efforts undertaken to identify blackleg QR loci in oilseed rape using linkage and association analysis. We update the knowledge on the possible mechanisms governing QR and the advances in searching for the underlying genes. Lastly, we lay out strategies to accelerate the genetic improvement of blackleg QR in oilseed rape using improved phenotyping approaches and genomic prediction tools.

Cross-species analysis between the maize smut fungi Ustilago maydis and Sporisorium reilianum highlights the role of transcriptional change of effector orthologs for virulence and disease

The constitution and regulation of effector repertoires shape host-microbe interactions. Ustilago maydis and Sporisorium reilianum are two closely related smut fungi, which both infect maize but cause distinct disease symptoms. Understanding how effector orthologs are regulated in these two pathogens can therefore provide insights into the evolution of different infection strategies. We tracked the infection progress of U. maydis and S. reilianum in maize leaves and used two distinct infection stages for cross-species RNA-sequencing analyses. We identified 207 of 335 one-to-one effector orthologs as differentially regulated during host colonization, which might reflect the distinct disease development strategies. Using CRISPR-Cas9-mediated gene conversion, we identified two differentially expressed effector orthologs with conserved function between two pathogens. Thus, differential expression of functionally conserved genes might contribute to species-specific adaptation and symptom development. Interestingly, another differentially expressed orthogroup (UMAG_05318/Sr10075) showed divergent protein function, providing a possible case for neofunctionalization. Collectively, we demonstrated that the diversification of effector genes in related pathogens can be caused both by alteration on the transcriptional level and through functional diversification of the encoded effector proteins.

Differential Regulation of Maize and Sorghum Orthologs in Response to the Fungal Pathogen Exserohilum turcicum

Pathogens that infect more than one host offer an opportunity to study how resistance mechanisms have evolved across different species. Exserohilum turcicum infects both maize and sorghum and the isolates are host-specific, offering a unique system to examine both compatible and incompatible interactions. We conducted transcriptional analysis of maize and sorghum in response to maize-specific and sorghum-specific E. turcicum isolates and identified functionally related co-expressed modules. Maize had a more robust transcriptional response than sorghum. E. turcicum responsive genes were enriched in core orthologs in both crops, but only up to 16% of core orthologs showed conserved expression patterns. Most changes in gene expression for the core orthologs, including hub genes, were lineage specific, suggesting a role for regulatory divergent evolution. We identified several defense-related shared differentially expressed (DE) orthologs with conserved expression patterns between the two crops, suggesting a role for parallel evolution of those genes in both crops. Many of the differentially expressed genes (DEGs) during the incompatible interaction were related to quantitative disease resistance (QDR). This work offers insights into how different hosts with relatively recent divergence interact with a common pathogen. Our results are important for developing resistance to this critical pathogen and understanding the evolution of host-pathogen interactions.

Silent control: microbial plant pathogens evade host immunity without coding sequence changes

Both animals and plants have evolved a robust immune system to surveil and defeat invading pathogenic microbes. Evasion of host immune surveillance is the key for pathogens to initiate successful infection. To evade the host immunity, plant pathogens evolved a variety of strategies such as masking themselves from host immune recognitions, blocking immune signaling transductions, reprogramming immune responses and adapting to immune microenvironmental changes. Gain of new virulence genes, sequence and structural variations enables plant pathogens to evade host immunity through changes in the genetic code. However, recent discoveries demonstrated that variations at the transcriptional, post-transcriptional, post-translational and glycome level enable pathogens to cope with the host immune system without coding sequence changes. The biochemical modification of pathogen associated molecular patterns and silencing of effector genes emerged as potent ways for pathogens to hide from host recognition. Altered processing in mRNA activities provide pathogens with resilience to microenvironment changes. Importantly, these hiding variants are directly or indirectly modulated by catalytic enzymes or enzymatic complexes and cannot be revealed by classical genomics alone. Unveiling these novel host evasion mechanisms in plant pathogens enables us to better understand the nature of plant disease and pinpoints strategies for rational diseases management in global food protection.

Genome-wide alternative splicing profiling in the fungal plant pathogen Sclerotinia sclerotiorum during the colonization of diverse host families

Sclerotinia sclerotiorum is a notorious generalist plant pathogen that threatens more than 600 host plants, including wild and cultivated species. The molecular bases underlying the broad compatibility of S. sclerotiorum with its hosts is not fully elucidated. In contrast to higher plants and animals, alternative splicing (AS) is not well studied in plant-pathogenic fungi. AS is a common regulated cellular process that increases cell protein and RNA diversity. In this study, we annotated spliceosome genes in the genome of S. sclerotiorum and characterized their expression in vitro and during the colonization of six host species. Several spliceosome genes were differentially expressed in planta, suggesting that AS was altered during infection. Using stringent parameters, we identified 1,487 S. sclerotiorum genes differentially expressed in planta and exhibiting alternative transcripts. The most common AS events during the colonization of all plants were retained introns and the alternative 3' receiver site. We identified S. sclerotiorum genes expressed in planta for which (a) the relative accumulation of alternative transcripts varies according to the host being colonized and (b) alternative transcripts harbour distinct protein domains. This notably included 42 genes encoding predicted secreted proteins showing high-confidence AS events. This study indicates that AS events are taking place in the plant pathogenic fungus S. sclerotiorum during the colonization of host plants and could generate functional diversity in the repertoire of proteins secreted by S. sclerotiorum during infection.