Large-scale transcriptomics to dissect 2 years of the life of a fungal phytopathogen interacting with its host plant.

BACKGROUND: The fungus Leptosphaeria maculans has an exceptionally long and complex relationship with its host plant, Brassica napus, during which it switches between different lifestyles, including asymptomatic, biotrophic, necrotrophic, and saprotrophic stages. The fungus is also exemplary of "two-speed" genome organisms in the genome of which gene-rich and repeat-rich regions alternate. Except for a few stages of plant infection under controlled conditions, nothing is known about the genes mobilized by the fungus throughout its life cycle, which may last several years in the field. RESULTS: We performed RNA-seq on samples corresponding to all stages of the interaction of L. maculans with its host plant, either alive or dead (stem residues after harvest) in controlled conditions or in field experiments under natural inoculum pressure, over periods of time ranging from a few days to months or years. A total of 102 biological samples corresponding to 37 sets of conditions were analyzed. We show here that about 9% of the genes of this fungus are highly expressed during its interactions with its host plant. These genes are distributed into eight well-defined expression clusters, corresponding to specific infection lifestyles or to tissue-specific genes. All expression clusters are enriched in effector genes, and one cluster is specific to the saprophytic lifestyle on plant residues. One cluster, including genes known to be involved in the first phase of asymptomatic fungal growth in leaves, is re-used at each asymptomatic growth stage, regardless of the type of organ infected. The expression of the genes of this cluster is repeatedly turned on and off during infection. Whatever their expression profile, the genes of these clusters are enriched in heterochromatin regions associated with H3K9me3 or H3K27me3 repressive marks. These findings provide support for the hypothesis that part of the fungal genes involved in niche adaptation is located in heterochromatic regions of the genome, conferring an extreme plasticity of expression. CONCLUSION: This work opens up new avenues for plant disease control, by identifying stage-specific effectors that could be used as targets for the identification of novel durable disease resistance genes, or for the in-depth analysis of chromatin remodeling during plant infection, which could be manipulated to interfere with the global expression of effector genes at crucial stages of plant infection.

Hatching of Globodera pallida Induced by Root Exudates Is Not Influenced by Soil Microbiota Composition.

Plant-parasitic nematodes are among the most harmful pests of cultivated crops causing important economic losses. The ban of chemical nematicides requires the development of alternative agroecological approaches to protect crops against nematodes. For cyst nematodes, egg hatching is stimulated by host plant root exudates. Inducing "suicide hatching" of nematode second-stage juveniles (J2), using root exudates in the absence of the host plant, may constitute an effective and innovative biocontrol method to control cyst nematodes. However, before considering the development of this approach, understanding the effect of soil biotic component on cyst nematode hatching by root exudates is a major issue. The effectiveness of this approach could be modulated by other soil organisms consuming root exudates for growth as soil microbiota, and this must be evaluated. To do that, four different native agricultural soils were selected based on their physicochemical properties and their microbiota composition were characterized by rDNA metabarcoding. To disentangle the effect of microbiota from that of soil on hatching, four recolonized artificial soils were obtained by inoculating a common sterile soil matrix with the microbiota proceeding from each agricultural soil. Each soil was then inoculated with cysts of the potato cyst nematode, Globodera pallida, and low or high doses of potato root exudates (PREs) were applied. After 40 days, viable J2 remaining in cysts were counted to determine the efficiency of root exudates to stimulate hatching in different soils. Results showed that (i) when physicochemical and microbiota compositions varied among native soils, the hatching rates remained very high albeit small differences were measured and no dose effect was detected and (ii) when only microbiota composition varied among recolonized soils, the hatching rates were also high at the highest dose of PREs, but a strong dose effect was highlighted. This study shows that abiotic and biotic factors may not compromise the development of methods based on suicide hatching of cyst nematodes, using root exudates, molecules inducing J2 hatch, or trap crops.

Soil microbiota influences clubroot disease by modulating Plasmodiophora brassicae and Brassica napus transcriptomes.

The contribution of surrounding plant microbiota to disease development has led to the 'pathobiome' concept, which represents the interaction between the pathogen, the host plant and the associated biotic microbial community, resulting or not in plant disease. The aim herein is to understand how the soil microbial environment may influence the functions of a pathogen and its pathogenesis, and the molecular response of the plant to the infection, with a dual-RNAseq transcriptomics approach. We address this question using Brassica napus and Plasmodiophora brassicae, the pathogen responsible for clubroot. A time-course experiment was conducted to study interactions between P. brassicae, two B. napus genotypes and three soils harbouring high, medium or low microbiota diversities and levels of richness. The soil microbial diversity levels had an impact on disease development (symptom levels and pathogen quantity). The P. brassicae and B. napus transcriptional patterns were modulated by these microbial diversities, these modulations being dependent on the host genotype plant and the kinetic time. The functional analysis of gene expressions allowed the identification of pathogen and plant host functions potentially involved in the change of plant disease level, such as pathogenicity-related genes (NUDIX effector) in P. brassicae and plant defence-related genes (glucosinolate metabolism) in B. napus.

Temporal dynamics of bacterial and fungal communities during the infection of Brassica rapa roots by the protist Plasmodiophora brassicae.

The temporal dynamics of rhizosphere and root microbiota composition was compared between healthy and infected Chinese cabbage plants by the pathogen Plasmodiophora brassicae. When inoculated with P. brassicae, disease was measured at five sampling dates from early root hair infection to late gall development. The first symptoms of clubroot disease appeared 14 days after inoculation (DAI) and increased drastically between 14 and 35 DAI. The structure of microbial communities associated to rhizosphere soil and root from healthy and inoculated plants was characterized through high-throughput DNA sequencing of bacterial (16S) and fungal (18S) molecular markers and compared at each sampling date. In healthy plants, Proteobacteria and Bacteroidetes bacterial phyla dominated the rhizosphere and root microbiota of Chinese cabbage. Rhizosphere bacterial communities contained higher abundances of Actinobacteria and Firmicutes compared to the roots. Moreover, a drastic shift of fungal communities of healthy plants occurred between the two last sampling dates, especially in plant roots, where most of Ascomycota fungi dominated until they were replaced by a fungus assigned to the Chytridiomycota phylum. Parasitic invasion by P. brassicae disrupted the rhizosphere and root-associated community assembly at a late step during the root secondary cortical infection stage of clubroot disease. At this stage, Flavisolibacter and Streptomyces in the rhizosphere, and Bacillus in the roots, were drastically less abundant upon parasite invasion. Rhizosphere of plants colonized by P. brassicae was significantly more invaded by the Chytridiomycota fungus, which could reflect a mutualistic relationship in this compartment between these two microorganisms.

De novo assembly and annotation of three Leptosphaeria genomes using Oxford Nanopore MinION sequencing.

Leptosphaeria maculans and Leptosphaeria biglobosa are ascomycete phytopathogens of Brassica napus (oilseed rape, canola). Here we report the complete sequence of three Leptosphaeria genomes (L. maculans JN3, L. maculans Nz-T4 and L. biglobosa G12-14). Nz-T4 and G12-14 genome assemblies were generated de novo and the reference JN3 genome assembly was improved using Oxford Nanopore MinION reads. The new assembly of L. biglobosa showed the existence of AT rich regions and pointed to a genome compartmentalization previously unsuspected following Illumina sequencing. Moreover nanopore sequencing allowed us to generate a chromosome-level assembly for the L. maculans reference isolate, JN3. The genome annotation was supported by integrating conserved proteins and RNA sequencing from Leptosphaeria-infected samples. The newly produced high-quality assemblies and annotations of those three Leptosphaeria genomes will allow further studies, notably focused on the tripartite interaction between L. maculans, L. biglobosa and oilseed rape. The discovery of as yet unknown effectors will notably allow progress in B. napus breeding towards L. maculans resistance.

Can soil microbial diversity influence plant metabolites and life history traits of a rhizophagous insect? A demonstration in oilseed rape.

Interactions between plants and phytophagous insects play an important part in shaping the biochemical composition of plants. Reciprocally plant metabolites can influence major life history traits in these insects and largely contribute to their fitness. Plant rhizospheric microorganisms are an important biotic factor modulating plant metabolites and adaptation to stress. While plant-insects or plant-microorganisms interactions and their consequences on the plant metabolite signature are well-documented, the impact of soil microbial communities on plant defenses against phytophagous insects remains poorly known. In this study, we used oilseed rape (Brassica napus) and the cabbage root fly (Delia radicum) as biological models to tackle this question. Even though D. radicum is a belowground herbivore as a larva, its adult life history traits depend on aboveground signals. We therefore tested whether soil microbial diversity influenced emergence rate and fitness but also fly oviposition behavior, and tried to link possible effects to modifications in leaf and root metabolites. Through a removal-recolonization experiment, 3 soil microbial modalities ("high," "medium," "low") were established and assessed through amplicon sequencing of 16S and 18S ribosomal RNA genes. The "medium" modality in the rhizosphere significantly improved insect development traits. Plant-microorganism interactions were marginally associated to modulations of root metabolites profiles, which could partly explain these results. We highlighted the potential role of plant-microbial interaction in plant defenses against Delia radicum. Rhizospheric microbial communities must be taken into account when analyzing plant defenses against herbivores, being either below or aboveground.

Different waves of effector genes with contrasted genomic location are expressed by Leptosphaeria maculans during cotyledon and stem colonization of oilseed rape.

Leptosphaeria maculans, the causal agent of stem canker disease, colonizes oilseed rape (Brassica napus) in two stages: a short and early colonization stage corresponding to cotyledon or leaf colonization, and a late colonization stage during which the fungus colonizes systemically and symptomlessly the plant during several months before stem canker appears. To date, the determinants of the late colonization stage are poorly understood; L. maculans may either successfully escape plant defences, leading to stem canker development, or the plant may develop an 'adult-stage' resistance reducing canker incidence. To obtain an insight into these determinants, we performed an RNA-sequencing (RNA-seq) pilot project comparing fungal gene expression in infected cotyledons and in symptomless or necrotic stems. Despite the low fraction of fungal material in infected stems, sufficient fungal transcripts were detected and a large number of fungal genes were expressed, thus validating the feasibility of the approach. Our analysis showed that all avirulence genes previously identified are under-expressed during stem colonization compared with cotyledon colonization. A validation RNA-seq experiment was then performed to investigate the expression of candidate effector genes during systemic colonization. Three hundred and seven 'late' effector candidates, under-expressed in the early colonization stage and over-expressed in the infected stems, were identified. Finally, our analysis revealed a link between the regulation of expression of effectors and their genomic location: the 'late' effector candidates, putatively involved in systemic colonization, are located in gene-rich genomic regions, whereas the 'early' effector genes, over-expressed in the early colonization stage, are located in gene-poor regions of the genome.

Crystal structure of the effector AvrLm4-7 of Leptosphaeria maculans reveals insights into its translocation into plant cells and recognition by resistance proteins.

The avirulence gene AvrLm4-7 of Leptosphaeria maculans, the causal agent of stem canker in Brassica napus (oilseed rape), confers a dual specificity of recognition by two resistance genes (Rlm4 and Rlm7) and is strongly involved in fungal fitness. In order to elucidate the biological function of AvrLm4-7 and understand the specificity of recognition by Rlm4 and Rlm7, the AvrLm4-7 protein was produced in Pichia pastoris and its crystal structure was determined. It revealed the presence of four disulfide bridges, but no close structural analogs could be identified. A short stretch of amino acids in the C terminus of the protein, (R/N)(Y/F)(R/S)E(F/W), was well-conserved among AvrLm4-7 homologs. Loss of recognition of AvrLm4-7 by Rlm4 is caused by the mutation of a single glycine to an arginine residue located in a loop of the protein. Loss of recognition by Rlm7 is governed by more complex mutational patterns, including gene loss or drastic modifications of the protein structure. Three point mutations altered residues in the well-conserved C-terminal motif or close to the glycine involved in Rlm4-mediated recognition, resulting in the loss of Rlm7-mediated recognition. Transient expression in Nicotiana benthamiana (tobacco) and particle bombardment experiments on leaves from oilseed rape suggested that AvrLm4-7 interacts with its cognate R proteins inside the plant cell, and can be translocated into plant cells in the absence of the pathogen. Translocation of AvrLm4-7 into oilseed rape leaves is likely to require the (R/N)(Y/F)(R/S)E(F/W) motif as well as an RAWG motif located in a nearby loop that together form a positively charged region.

Transposable element-assisted evolution and adaptation to host plant within the Leptosphaeria maculans-Leptosphaeria biglobosa species complex of fungal pathogens.

BACKGROUND: Many plant-pathogenic fungi have a tendency towards genome size expansion, mostly driven by increasing content of transposable elements (TEs). Through comparative and evolutionary genomics, five members of the Leptosphaeria maculans-Leptosphaeria biglobosa species complex (class Dothideomycetes, order Pleosporales), having different host ranges and pathogenic abilities towards cruciferous plants, were studied to infer the role of TEs on genome shaping, speciation, and on the rise of better adapted pathogens. RESULTS: L. maculans 'brassicae', the most damaging species on oilseed rape, is the only member of the species complex to have a TE-invaded genome (32.5%) compared to the other members genomes (<4%). These TEs had an impact at the structural level by creating large TE-rich regions and are suspected to have been instrumental in chromosomal rearrangements possibly leading to speciation. TEs, associated with species-specific genes involved in disease process, also possibly had an incidence on evolution of pathogenicity by promoting translocations of effector genes to highly dynamic regions and thus tuning the regulation of effector gene expression in planta. CONCLUSIONS: Invasion of L. maculans 'brassicae' genome by TEs followed by bursts of TE activity allowed this species to evolve and to better adapt to its host, making this genome species a peculiarity within its own species complex as well as in the Pleosporales lineage.

Epigenetic control of effector gene expression in the plant pathogenic fungus Leptosphaeria maculans.

Plant pathogens secrete an arsenal of small secreted proteins (SSPs) acting as effectors that modulate host immunity to facilitate infection. SSP-encoding genes are often located in particular genomic environments and show waves of concerted expression at diverse stages of plant infection. To date, little is known about the regulation of their expression. The genome of the Ascomycete Leptosphaeria maculans comprises alternating gene-rich GC-isochores and gene-poor AT-isochores. The AT-isochores harbor mosaics of transposable elements, encompassing one-third of the genome, and are enriched in putative effector genes that present similar expression patterns, namely no expression or low-level expression during axenic cultures compared to strong induction of expression during primary infection of oilseed rape (Brassica napus). Here, we investigated the involvement of one specific histone modification, histone H3 lysine 9 methylation (H3K9me3), in epigenetic regulation of concerted effector gene expression in L. maculans. For this purpose, we silenced the expression of two key players in heterochromatin assembly and maintenance, HP1 and DIM-5 by RNAi. By using HP1-GFP as a heterochromatin marker, we observed that almost no chromatin condensation is visible in strains in which LmDIM5 was silenced by RNAi. By whole genome oligoarrays we observed overexpression of 369 or 390 genes, respectively, in the silenced-LmHP1 and -LmDIM5 transformants during growth in axenic culture, clearly favouring expression of SSP-encoding genes within AT-isochores. The ectopic integration of four effector genes in GC-isochores led to their overexpression during growth in axenic culture. These data strongly suggest that epigenetic control, mediated by HP1 and DIM-5, represses the expression of at least part of the effector genes located in AT-isochores during growth in axenic culture. Our hypothesis is that changes of lifestyle and a switch toward pathogenesis lift chromatin-mediated repression, allowing a rapid response to new environmental conditions.

Incidence of genome structure, DNA asymmetry, and cell physiology on T-DNA integration in chromosomes of the phytopathogenic fungus Leptosphaeria maculans.

The ever-increasing generation of sequence data is accompanied by unsatisfactory functional annotation, and complex genomes, such as those of plants and filamentous fungi, show a large number of genes with no predicted or known function. For functional annotation of unknown or hypothetical genes, the production of collections of mutants using Agrobacterium tumefaciens-mediated transformation (ATMT) associated with genotyping and phenotyping has gained wide acceptance. ATMT is also widely used to identify pathogenicity determinants in pathogenic fungi. A systematic analysis of T-DNA borders was performed in an ATMT-mutagenized collection of the phytopathogenic fungus Leptosphaeria maculans to evaluate the features of T-DNA integration in its particular transposable element-rich compartmentalized genome. A total of 318 T-DNA tags were recovered and analyzed for biases in chromosome and genic compartments, existence of CG/AT skews at the insertion site, and occurrence of microhomologies between the T-DNA left border (LB) and the target sequence. Functional annotation of targeted genes was done using the Gene Ontology annotation. The T-DNA integration mainly targeted gene-rich, transcriptionally active regions, and it favored biological processes consistent with the physiological status of a germinating spore. T-DNA integration was strongly biased toward regulatory regions, and mainly promoters. Consistent with the T-DNA intranuclear-targeting model, the density of T-DNA insertion correlated with CG skew near the transcription initiation site. The existence of microhomologies between promoter sequences and the T-DNA LB flanking sequence was also consistent with T-DNA integration to host DNA mediated by homologous recombination based on the microhomology-mediated end-joining pathway.

Effector diversification within compartments of the Leptosphaeria maculans genome affected by Repeat-Induced Point mutations.

Fungi are of primary ecological, biotechnological and economic importance. Many fundamental biological processes that are shared by animals and fungi are studied in fungi due to their experimental tractability. Many fungi are pathogens or mutualists and are model systems to analyse effector genes and their mechanisms of diversification. In this study, we report the genome sequence of the phytopathogenic ascomycete Leptosphaeria maculans and characterize its repertoire of protein effectors. The L. maculans genome has an unusual bipartite structure with alternating distinct guanine and cytosine-equilibrated and adenine and thymine (AT)-rich blocks of homogenous nucleotide composition. The AT-rich blocks comprise one-third of the genome and contain effector genes and families of transposable elements, both of which are affected by repeat-induced point mutation, a fungal-specific genome defence mechanism. This genomic environment for effectors promotes rapid sequence diversification and underpins the evolutionary potential of the fungus to adapt rapidly to novel host-derived constraints.