Genome structure and reproductive behaviour influence the evolutionary potential of a fungal phytopathogen

Transposable Element Genomic Fissuring in Pyrenophora teres Is Associated With Genome Expansion and Dynamics of Host-Pathogen Genetic Interactions

Pyrenophora teres, P. teres f. teres (PTT) and P. teres f. maculata (PTM) cause significant diseases in barley, but little is known about the large-scale genomic differences that may distinguish the two forms. Comprehensive genome assemblies were constructed from long DNA reads, optical and genetic maps. As repeat masking in fungal genomes influences the final gene annotations, an accurate and reproducible pipeline was developed to ensure comparability between isolates. The genomes of the two forms are highly collinear, each composed of 12 chromosomes. Genome evolution in P. teres is characterized by genome fissuring through the insertion and expansion of transposable elements (TEs), a process that isolates blocks of genic sequence. The phenomenon is particularly pronounced in PTT, which has a larger, more repetitive genome than PTM and more recent transposon activity measured by the frequency and size of genome fissures. PTT has a longer cultivated host association and, notably, a greater range of host-pathogen genetic interactions compared to other Pyrenophora spp., a property which associates better with genome size than pathogen lifestyle. The two forms possess similar complements of TE families with Tc1/Mariner and LINE-like Tad-1 elements more abundant in PTT. Tad-1 was only detectable as vestigial fragments in PTM and, within the forms, differences in genome sizes and the presence and absence of several TE families indicated recent lineage invasions. Gene differences between P. teres forms are mainly associated with gene-sparse regions near or within TE-rich regions, with many genes possessing characteristics of fungal effectors. Instances of gene interruption by transposons resulting in pseudogenization were detected in PTT. In addition, both forms have a large complement of secondary metabolite gene clusters indicating significant capacity to produce an array of different molecules. This study provides genomic resources for functional genetics to help dissect factors underlying the host-pathogen interactions.

A Six-Year Investigation of the Dynamics of Avirulence Allele Profiles, Blackleg Incidence, and Mating Type Alleles of Leptosphaeria maculans Populations Associated with Canola Crops in Manitoba, Canada

Blackleg, caused by the fungal pathogen Leptosphaeria maculans, is one of the most economically important diseases of canola (Brassica napus, oilseed rape) worldwide. This study assessed incidence of blackleg, the avirulence allele, and mating type distributions of L. maculans isolates collected in commercial canola fields in Manitoba, Canada, from 2010 to 2015. A total of 956 L. maculans isolates were collected from 2010 to 2015 to determine the presence of 12 avirulence alleles using differential canola cultivars and/or PCR assays specific for each avirulence allele. AvrLm2, AvrLm4, AvrLm5, AvrLm6, AvrLm7, AvrLm11, and AvrLmS were detected at frequencies ranging from 97 to 33%, where the AvrLm1, AvrLm3, AvrLm9, AvrLepR1, and AvrLepR2 alleles were the least abundant. When the race structure was examined, a total of 170 races were identified among the 956 isolates, with three major races, AvrLm-2-4-5-6-7-11, AvrLm-2-4-5-6-7-11-S, and Avr-1-4-5-6-7-11-(S) accounting for 15, 10, and 6% of the total fungal population, respectively. The distribution of the mating type alleles (MAT1-1 and MAT1-2) indicated that sexual reproduction was not inhibited in any of the nine Manitoba regions in any of the years L. maculans isolates were collected.

Distinct domains of the AVRPM3A2/F2 avirulence protein from wheat powdery mildew are involved in immune receptor recognition and putative effector function

Recognition of the AVRPM3A2/F2 avirulence protein from powdery mildew by the wheat PM3A/F immune receptor induces a hypersensitive response after co-expression in Nicotiana benthamiana. The molecular determinants of this interaction and how they shape natural AvrPm3a2/f2 allelic diversity are unknown. We sequenced the AvrPm3a2/f2 gene in a worldwide collection of 272 mildew isolates. Using the natural polymorphisms of AvrPm3a2/f2 as well as sequence information from related gene family members, we tested 85 single-residue-altered AVRPM3A2/F2 variants with PM3A, PM3F and PM3FL456P/Y458H (modified for improved signaling) in Nicotiana benthamiana for effects on recognition. An intact AvrPm3a2/f2 gene was found in all analyzed isolates and the protein variant recognized by PM3A/F occurred globally at high frequencies. Single-residue alterations in AVRPM3A2/F2 mostly disrupted, but occasionally enhanced, the recognition response by PM3A, PM3F and PM3FL456P/Y458H . Residues enhancing hypersensitive responses constituted a protein domain separate from both naturally occurring polymorphisms and positively selected residues of the gene family. These results demonstrate the utility of using gene family sequence diversity to screen residues for their role in recognition. This approach identified a putative interaction surface in AVRPM3A2/F2 not polymorphic in natural alleles. We conclude that molecular mechanisms besides recognition drive AvrPm3a2/f2 diversification.

Fine-Scale Recombination Maps of Fungal Plant Pathogens Reveal Dynamic Recombination Landscapes and Intragenic Hotspots

Meiotic recombination is an important driver of evolution. Variability in the intensity of recombination across chromosomes can affect sequence composition, nucleotide variation, and rates of adaptation. In many organisms, recombination events are concentrated within short segments termed recombination hotspots. The variation in recombination rate and positions of recombination hotspot can be studied using population genomics data and statistical methods. In this study, we conducted population genomics analyses to address the evolution of recombination in two closely related fungal plant pathogens: the prominent wheat pathogen Zymoseptoria tritici and a sister species infecting wild grasses Z. ardabiliae We specifically addressed whether recombination landscapes, including hotspot positions, are conserved in the two recently diverged species and if recombination contributes to rapid evolution of pathogenicity traits. We conducted a detailed simulation analysis to assess the performance of methods of recombination rate estimation based on patterns of linkage disequilibrium, in particular in the context of high nucleotide diversity. Our analyses reveal overall high recombination rates, a lack of suppressed recombination in centromeres, and significantly lower recombination rates on chromosomes that are known to be accessory. The comparison of the recombination landscapes of the two species reveals a strong correlation of recombination rate at the megabase scale, but little correlation at smaller scales. The recombination landscapes in both pathogen species are dominated by frequent recombination hotspots across the genome including coding regions, suggesting a strong impact of recombination on gene evolution. A significant but small fraction of these hotspots colocalize between the two species, suggesting that hotspot dynamics contribute to the overall pattern of fast evolving recombination in these species.

The birth and death of effectors in rapidly evolving filamentous pathogen genomes

Plant pathogenic fungi and oomycetes are major risks to food security due to their evolutionary success in overcoming plant defences. Pathogens produce effectors to interfere with host defences and metabolism. These effectors are often encoded in rapidly evolving compartments of the genome. We review how effector genes emerged and were lost in pathogen genomes drawing on the links between effector evolution and chromosomal rearrangements. Some new effectors entered pathogen genomes via horizontal transfer or introgression. However, new effector functions also arose through gene duplication or from previously non-coding sequences. The evolutionary success of an effector is tightly linked to its transcriptional regulation during host colonization. Some effectors converged on an epigenetic control of expression imposed by genomic defences against transposable elements. Transposable elements were also drivers of effector diversification and loss that led to mosaics in effector presence-absence variation. Such effector mosaics within species was the foundation for rapid pathogen adaptation.

Current Status and Challenges in Identifying Disease Resistance Genes in Brassica napus

Brassica napus is an economically important crop across different continents including temperate and subtropical regions in Europe, Canada, South Asia, China and Australia. Its widespread cultivation also brings setbacks as it plays host to fungal, oomycete and chytrid pathogens that can lead to serious yield loss. For sustainable crop production, identification of resistance (R) genes in B. napus has become of critical importance. In this review, we discuss four key pathogens affecting Brassica crops: Clubroot (Plasmodiophora brassicae), Blackleg (Leptosphaeria maculans and L. biglobosa), Sclerotinia Stem Rot (Sclerotinia sclerotiorum), and Downy Mildew (Hyaloperonospora parasitica). We first review current studies covering prevalence of these pathogens on Brassica crops and highlight the R genes and QTL that have been identified from Brassica species against these pathogens. Insights into the relationships between the pathogen and its Brassica host, the unique host resistance mechanisms and how these affect resistance outcomes is also presented. We discuss challenges in identification and deployment of R genes in B. napus in relation to highly specific genetic interactions between host subpopulations and pathogen pathotypes and emphasize the need for common or shared techniques and research materials or tighter collaboration between researchers to reconcile the inconsistencies in the research outcomes. Using current genomics tools, we provide examples of how characterization and cloning of R genes in B. napus can be carried out more effectively. Lastly, we put forward strategies to breed resistant cultivars through introgressions supported by genomic approaches and suggest prospects that can be implemented in the future for a better, pathogen-resistant B. napus.

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.

Sources of Fungal Genetic Variation and Associating It with Phenotypic Diversity

The first eukaryotic genome to be sequenced was fungal, and there continue to be more sequenced genomes in the kingdom Fungi than in any other eukaryotic kingdom. Comparison of these genomes reveals many sources of genetic variation, from single nucleotide polymorphisms to horizontal gene transfer and on to changes in the arrangement and number of chromosomes, not to mention endofungal bacteria and viruses. Population genomics shows that all sources generate variation all the time and implicate natural selection as the force maintaining genome stability. Variation in wild populations is a rich resource for associating genetic variation with phenotypic variation, whether through quantitative trait locus mapping, genome-wide association studies, or reverse ecology. Subjects of studies associating genetic and phenotypic variation include model fungi, e.g., Saccharomyces and Neurospora, but pioneering studies have also been made with fungi pathogenic to plants, e.g., Pyricularia (= Magnaporthe), Zymoseptoria, and Fusarium, and to humans, e.g., Coccidioides, Cryptococcus, and Candida.

Transposable Elements Direct The Coevolution between Plants and Microbes

Transposable elements are powerful drivers of genome evolution in many eukaryotes. Although they are mostly considered as 'selfish' genetic elements, increasing evidence suggests that they contribute to genetic variability; particularly under stress conditions. Over the past few years, the role of transposable elements during host-microbe interactions has been recognised. It has been proposed that many pathogenic microbes have evolved a 'two-speed' genome with regions that show increased variability and that are enriched in transposable elements and pathogenicity-related genes. Plants similarly display structured genomes with transposable-element-rich regions that mediate accelerated evolution. Immune receptor genes typically reside in such regions. Various mechanisms have recently been identified through which transposable elements contribute to the coevolution between plants and their associated microbes.

Complex Interactions between Fungal Avirulence Genes and Their Corresponding Plant Resistance Genes and Consequences for Disease Resistance Management

During infection, pathogens secrete an arsenal of molecules, collectively called effectors, key elements of pathogenesis which modulate innate immunity of the plant and facilitate infection. Some of these effectors can be recognized directly or indirectly by resistance (R) proteins from the plant and are then called avirulence (AVR) proteins. This recognition usually triggers defense responses including the hypersensitive response and results in resistance of the plant. R-AVR gene interactions are frequently exploited in the field to control diseases. Recently, the availability of fungal genomes has accelerated the identification of AVR genes in plant pathogenic fungi, including in fungi infecting agronomically important crops. While single AVR genes recognized by their corresponding R gene were identified, more and more complex interactions between AVR and R genes are reported (e.g., AVR genes recognized by several R genes, R genes recognizing several AVR genes in distinct organisms, one AVR gene suppressing recognition of another AVR gene by its corresponding R gene, two cooperating R genes both necessary to recognize an AVR gene). These complex interactions were particularly reported in pathosystems showing a long co-evolution with their host plant but could also result from the way agronomic crops were obtained and improved (e.g., through interspecific hybridization or introgression of resistance genes from wild related species into cultivated crops). In this review, we describe some complex R-AVR interactions between plants and fungi that were recently reported and discuss their implications for AVR gene evolution and R gene management.

Life, death and rebirth of avirulence effectors in a fungal pathogen of Brassica crops, Leptosphaeria maculans

Contents 526 I. 526 II. 527 III. 527 IV. 529 V. 529 VI. 530 VII. 530 531 References 531 SUMMARY: In agricultural systems, major (R) genes for resistance in plants exert strong selection pressure on cognate/corresponding avirulence effector genes of phytopathogens. However, a complex interplay often exists between trade-offs linked to effector function and the need to escape R gene recognition. Here, using the Leptosphaeria maculans-oilseed rape pathosystem we review evolution of effectors submitted to multiple resistance gene selection. Characteristics of this pathosystem include a crop in which resistance genes have been deployed intensively resulting in 'boom and bust' cycles; a fungal pathogen with a high adaptive potential in which seven avirulence genes are cloned and for which population surveys have been coupled with molecular analysis of events responsible for virulence. The mode of evolution of avirulence genes, all located in dispensable parts of the 'two-speed' genome, is a highly dynamic gene-specific process. In some instances, avirulence genes are readily deleted under selection. However, others, even when located in the most plastic genome regions, undergo only limited point mutations or their avirulence phenotype is 'camouflaged' by another avirulence gene. Thus, while hundreds of effector genes are present, some effectors are likely to have an important and nonredundant function, suggesting functional redundancy and dispensability of effectors might not be the rule.

Unusual evolutionary mechanisms to escape effector-triggered immunity in the fungal phytopathogen Leptosphaeria maculans

Leptosphaeria maculans is the fungus responsible for the stem canker disease of oilseed rape (Brassica napus). AvrLm3 and AvrLm4-7, two avirulence effector genes of L. maculans, are involved in an unusual relationship: AvrLm4-7 suppresses the Rlm3-mediated resistance. Here, we assessed AvrLm3 polymorphism in a collection of 235 L. maculans isolates. No field isolates exhibited deletion or inactivating mutations in AvrLm3, as observed for other L. maculans avirulence genes. Eleven isoforms of the AvrLm3 protein were found. In isolates virulent towards both Rlm3 and Rlm7 (a3a7), the loss of the Rlm3-mediated resistance response was due to two distinct mechanisms. First, when AvrLm4-7 was inactivated (deletion or inactivating mutations), amino acid substitutions in AvrLm3 generated virulent isoforms of the protein. Second, when only point mutations were observed in AvrLm4-7, a3a7 isolates still contained an avirulent allele of AvrLm3. Directed mutagenesis confirmed that some point mutations in AvrLm4-7 were sufficient for the fungus to escape Rlm7-mediated resistance while maintaining the suppression of the AvrLm3 phenotype. Signatures of positive selection were also identified in AvrLm3. The complex evolutionary mechanisms enabling L. maculans to escape Rlm3-mediated resistance while preserving AvrLm3 integrity, along with observed reduced aggressiveness of isolates silenced for AvrLm3, serves to emphasize the importance of this effector in pathogenicity towards B. napus. While the common response to resistance gene pressure is local selection of isolates depleted in the cognate avirulence gene, this example contributes to complexify the gene-for-gene concept of plant-pathogen evolution with a 'camouflaged' model allowing retention of nondispensable avirulence effectors.

New Multilocus Variable-Number Tandem-Repeat Analysis (MLVA) Scheme for Fine-Scale Monitoring and Microevolution-Related Study of Ralstonia pseudosolanacearum Phylotype I Populations

Bacterial wilt caused by the Ralstonia solanacearum species complex (RSSC) is considered one of the most harmful plant diseases in the world. Special attention should be paid to R. pseudosolanacearum phylotype I due to its large host range, its worldwide distribution, and its high evolutionary potential. So far, the molecular epidemiology and population genetics of this bacterium are poorly understood. Until now, the genetic structure of the RSSC has been analyzed on the worldwide and regional scales. Emerging questions regarding evolutionary forces in RSSC adaptation to hosts now require genetic markers that are able to monitor RSSC field populations. In this study, we aimed to evaluate the multilocus variable-number tandem-repeat analysis (MLVA) approach for its ability to discriminate genetically close phylotype I strains and for population genetics studies. We developed a new MLVA scheme (MLVA-7) allowing us to genotype 580 R. pseudosolanacearum phylotype I strains extracted from susceptible and resistant hosts and from different habitats (stem, soil, and rhizosphere). Based on specificity, polymorphism, and the amplification success rate, we selected seven fast-evolving variable-number tandem-repeat (VNTR) markers. The newly developed MLVA-7 scheme showed higher discriminatory power than the previously published MLVA-13 scheme when applied to collections sampled from the same location on different dates and to collections from different locations on very small scales. Our study provides a valuable tool for fine-scale monitoring and microevolution-related study of R. pseudosolanacearum phylotype I populations.IMPORTANCE Understanding the evolutionary dynamics of adaptation of plant pathogens to new hosts or ecological niches has become a key point for the development of innovative disease management strategies, including durable resistance. Whereas the molecular mechanisms underlying virulence or pathogenicity changes have been studied thoroughly, the population genetics of plant pathogen adaptation remains an open, unexplored field, especially for plant-pathogenic bacteria. MLVA has become increasingly popular for epidemiosurveillance and molecular epidemiology studies of plant pathogens. However, this method has been used mostly for genotyping and identification on a regional or global scale. In this study, we developed a new MLVA scheme, targeting phylotype I of the soilborne Ralstonia solanacearum species complex (RSSC), specifically to address the bacterial population genetics on the field scale. Such a MLVA scheme, based on fast-evolving loci, may be a tool of choice for field experimental evolution and spatial genetics studies.

The Evolution of Orphan Regions in Genomes of a Fungal Pathogen of Wheat

Fungal plant pathogens rapidly evolve virulence on resistant hosts through mutations in genes encoding proteins that modulate the host immune responses. The mutational spectrum likely includes chromosomal rearrangements responsible for gains or losses of entire genes. However, the mechanisms creating adaptive structural variation in fungal pathogen populations are poorly understood. We used complete genome assemblies to quantify structural variants segregating in the highly polymorphic fungal wheat pathogen Zymoseptoria tritici The genetic basis of virulence in Z. tritici is complex, and populations harbor significant genetic variation for virulence; hence, we aimed to identify whether structural variation led to functional differences. We combined single-molecule real-time sequencing, genetic maps, and transcriptomics data to generate a fully assembled and annotated genome of the highly virulent field isolate 3D7. Comparative genomics analyses against the complete reference genome IPO323 identified large chromosomal inversions and the complete gain or loss of transposable-element clusters, explaining the extensive chromosomal-length polymorphisms found in this species. Both the 3D7 and IPO323 genomes harbored long tracts of sequences exclusive to one of the two genomes. These orphan regions contained 296 genes unique to the 3D7 genome and not previously known for this species. These orphan genes tended to be organized in clusters and showed evidence of mutational decay. Moreover, the orphan genes were enriched in genes encoding putative effectors and included a gene that is one of the most upregulated putative effector genes during wheat infection. Our study showed that this pathogen species harbored extensive chromosomal structure polymorphism that may drive the evolution of virulence.

IMPORTANCE

Pathogen outbreak populations often harbor previously unknown genes conferring virulence. Hence, a key puzzle of rapid pathogen evolution is the origin of such evolutionary novelty in genomes. Chromosomal rearrangements and structural variation in pathogen populations likely play a key role. However, identifying such polymorphism is challenging, as most genome-sequencing approaches only yield information about point mutations. We combined long-read technology and genetic maps to assemble the complete genome of a strain of a highly polymorphic fungal pathogen of wheat. Comparisons against the reference genome of the species showed substantial variation in the chromosome structure and revealed large regions unique to each assembled genome. These regions were enriched in genes encoding likely effector proteins, which are important components of pathogenicity. Our study showed that pathogen populations harbor extensive polymorphism at the chromosome level and that this polymorphism can be a source of adaptive genetic variation in pathogen evolution.

Multi-environment QTL studies suggest a role for cysteine-rich protein kinase genes in quantitative resistance to blackleg disease in Brassica napus

BACKGROUND

Resistance to the blackleg disease of Brassica napus (canola/oilseed rape), caused by the hemibiotrophic fungal pathogen Leptosphaeria maculans, is determined by both race-specific resistance (R) genes and quantitative resistance loci (QTL), or adult-plant resistance (APR). While the introgression of R genes into breeding material is relatively simple, QTL are often detected sporadically, making them harder to capture in breeding programs. For the effective deployment of APR in crop varieties, resistance QTL need to have a reliable influence on phenotype in multiple environments and be well defined genetically to enable marker-assisted selection (MAS).

RESULTS

Doubled-haploid populations produced from the susceptible B. napus variety Topas and APR varieties AG-Castle and AV-Sapphire were analysed for resistance to blackleg in two locations over 3 and 4 years, respectively. Three stable QTL were detected in each population, with two loci appearing to be common to both APR varieties. Physical delineation of three QTL regions was sufficient to identify candidate defense-related genes, including a cluster of cysteine-rich receptor-like kinases contained within a 49 gene QTL interval on chromosome A01. Individual L. maculans isolates were used to define the physical intervals for the race-specific R genes Rlm3 and Rlm4 and to identify QTL common to both field studies and the cotyledon resistance response.

CONCLUSION

Through multi-environment QTL analysis we have identified and delineated four significant and stable QTL suitable for MAS of quantitative blackleg resistance in B. napus, and identified candidate genes which potentially play a role in quantitative defense responses to L. maculans.

Transposons passively and actively contribute to evolution of the two-speed genome of a fungal pathogen

Genomic plasticity enables adaptation to changing environments, which is especially relevant for pathogens that engage in “arms races” with their hosts. In many pathogens, genes mediating virulence cluster in highly variable, transposon-rich, physically distinct genomic compartments. However, understanding of the evolution of these compartments, and the role of transposons therein, remains limited. Here, we show that transposons are the major driving force for adaptive genome evolution in the fungal plant pathogen Verticillium dahliae. We show that highly variable lineage-specific (LS) regions evolved by genomic rearrangements that are mediated by erroneous double-strand repair, often utilizing transposons. We furthermore show that recent genetic duplications are enhanced in LS regions, against an older episode of duplication events. Finally, LS regions are enriched in active transposons, which contribute to local genome plasticity. Thus, we provide evidence for genome shaping by transposons, both in an active and passive manner, which impacts the evolution of pathogen virulence.

A game of hide and seek between avirulence genes AvrLm4-7 and AvrLm3 in Leptosphaeria maculans

Extending the durability of plant resistance genes towards fungal pathogens is a major challenge. We identified and investigated the relationship between two avirulence genes of Leptosphaeria maculans, AvrLm3 and AvrLm4-7. When an isolate possesses both genes, the Rlm3-mediated resistance of oilseed rape (Brassica napus) is not expressed due to the presence of AvrLm4-7 but virulent isolates toward Rlm7 recover the AvrLm3 phenotype. Combining genetic and genomic approaches (genetic mapping, RNA-seq, BAC (bacterial artificial chromosome) clone sequencing and de novo assembly) we cloned AvrLm3, a telomeric avirulence gene of L. maculans. AvrLm3 is located in a gap of the L. maculans reference genome assembly, is surrounded by repeated elements, encodes for a small secreted cysteine-rich protein and is highly expressed at early infection stages. Complementation and silencing assays validated the masking effect of AvrLm4-7 on AvrLm3 recognition by Rlm3 and we showed that the presence of AvrLm4-7 does not impede AvrLm3 expression in planta. Y2H assays suggest the absence of physical interaction between the two avirulence proteins. This unusual interaction is the basis for field experiments aiming to evaluate strategies that increase Rlm7 durability.

A Large Family of AvrLm6-like Genes in the Apple and Pear Scab Pathogens, Venturia inaequalis and Venturia pirina

Venturia inaequalis and V. pirina are Dothideomycete fungi that cause apple scab and pear scab disease, respectively. Whole genome sequencing of V. inaequalis and V. pirina isolates has revealed predicted proteins with sequence similarity to AvrLm6, a Leptosphaeria maculans effector that triggers a resistance response in Brassica napus and B. juncea carrying the resistance gene, Rlm6. AvrLm6-like genes are present as large families (>15 members) in all sequenced strains of V. inaequalis and V. pirina, while in L. maculans, only AvrLm6 and a single paralog have been identified. The Venturia AvrLm6-like genes are located in gene-poor regions of the genomes, and mostly in close proximity to transposable elements, which may explain the expansion of these gene families. An AvrLm6-like gene from V. inaequalis with the highest sequence identity to AvrLm6 was unable to trigger a resistance response in Rlm6-carrying B. juncea. RNA-seq and qRT-PCR gene expression analyses, of in planta- and in vitro-grown V. inaequalis, has revealed that many of the AvrLm6-like genes are expressed during infection. An AvrLm6 homolog from V. inaequalis that is up-regulated during infection was shown (using an eYFP-fusion protein construct) to be localized to the sub-cuticular stroma during biotrophic infection of apple hypocotyls.

The two-speed genomes of filamentous pathogens: waltz with plants

Fungi and oomycetes include deep and diverse lineages of eukaryotic plant pathogens. The last 10 years have seen the sequencing of the genomes of a multitude of species of these so-called filamentous plant pathogens. Already, fundamental concepts have emerged. Filamentous plant pathogen genomes tend to harbor large repertoires of genes encoding virulence effectors that modulate host plant processes. Effector genes are not randomly distributed across the genomes but tend to be associated with compartments enriched in repetitive sequences and transposable elements. These findings have led to the 'two-speed genome' model in which filamentous pathogen genomes have a bipartite architecture with gene sparse, repeat rich compartments serving as a cradle for adaptive evolution. Here, we review this concept and discuss how plant pathogens are great model systems to study evolutionary adaptations at multiple time scales. We will also introduce the next phase of research on this topic.

Rapid identification of the Leptosphaeria maculans avirulence gene AvrLm2 using an intraspecific comparative genomics approach

Five avirulence genes from Leptosphaeria maculans, the causal agent of blackleg of canola (Brassica napus), have been identified previously through map-based cloning. In this study, a comparative genomic approach was used to clone the previously mapped AvrLm2. Given the lack of a presence-absence gene polymorphism coincident with the AvrLm2 phenotype, 36 L. maculans isolates were resequenced and analysed for single-nucleotide polymorphisms (SNPs) in predicted small secreted protein-encoding genes present within the map interval. Three SNPs coincident with the AvrLm2 phenotype were identified within LmCys1, previously identified as a putative effector-coding gene. Complementation of a virulent isolate with LmCys1, as the candidate AvrLm2 allele, restored the avirulent phenotype on Rlm2-containing B. napus lines. AvrLm2 encodes a small cysteine-rich protein with low similarity to other proteins in the public databases. Unlike other avirulence genes, AvrLm2 resides in a small GC island within an AT-rich isochore of the genome, and was never found to be deleted completely in virulent isolates.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-891) contains supplementary material, which is available to authorized users.

An avirulence gene, AvrLmJ1, from the blackleg fungus, Leptosphaeria maculans, confers avirulence to Brassica juncea cultivars

The fungus Leptosphaeria maculans causes blackleg of Brassica species. Here, we report the mapping and subsequent cloning of an avirulence gene from L. maculans. This gene, termed AvrLmJ1, confers avirulence towards all three Brassica juncea cultivars tested. Analysis of RNA-seq data showed that AvrLmJ1 is housed in a region of the L. maculans genome which contains only one gene that is highly expressed in planta. The closest genes are 57 and 33 kb away and, like other avirulence genes of L. maculans, AvrLmJ1 is located within an AT-rich, gene-poor region of the genome. The encoded protein is 141 amino acids, has a predicted signal peptide and is cysteine rich. Two virulent isolates contain a premature stop codon in AvrLmJ1. Complementation of an isolate that forms cotyledonary lesions on B. juncea with the wild-type allele of AvrLmJ1 confers avirulence towards all three B. juncea cultivars tested, suggesting that the gene may confer species-specific avirulence activity.

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.

[Transposable elements reshaping genomes and favouring the evolutionary and adaptive potential of fungal phytopathogens]

Phytopathogenic fungi are a major threat for global food security and show an extreme plasticity in pathogenicity behaviours. They often have a high adaptive potential allowing them to rapidly counteract the control methods used by men in agrosystems. In this paper, we evaluate the link between genome plasticity and adaptive potential using genomics and comparative genomics approaches. Our model is the evolutionary series Leptosphaeria maculans-Leptosphaeria biglobosa, encompassing five distinct entities, whose conspecificity or heterospecificity status is unclear, and which all are pathogens of cruciferous plants. They however differ by their host range and pathogenicity. Compared to other species of the species complex, the species best adapted to oilseed rape, L. maculans "brassicae", causing important losses in the crop, has a genome that was submitted to a recent and massive burst of transposition by a few families of transposable elements (TEs). Whether the genome invasion contributed to speciation is still unclear to-date but there is a coincidence between this burst of TEs and divergence between two species. This TE burst contributed to diversification of effector proteins and thus to generation of novel pathogenic specificities. In addition, the location of effector genes within genome regions enriched in TEs has direct consequences on adaptation to plant resistance and favours a multiplicity of mutation events allowing "breakdown" of resistance. These data are substantiated by other examples in the literature showing that fungi tend to have a "two-speed" genome, in which a plastic compartment enriched in TE host genes is involved in pathogenicity and adaptation to host.

Quantitative resistance affects the speed of frequency increase but not the diversity of the virulence alleles overcoming a major resistance gene to Leptosphaeria maculans in oilseed rape

Quantitative resistance mediated by multiple genetic factors has been shown to increase the potential for durability of major resistance genes. This was demonstrated in the Leptosphaeria maculans/Brassica napus pathosystem in a 5year recurrent selection field experiment on lines harboring the qualitative resistance gene Rlm6 combined or not with quantitative resistance. The quantitative resistance limited the size of the virulent isolate population. In this study we continued this recurrent selection experiment in the same way to examine whether the pathogen population could adapt and render the major gene ineffective in the longer term. The cultivars Eurol, with a susceptible background, and Darmor, with quantitative resistance, were used. We confirmed that the combination of qualitative and quantitative resistance is an effective approach for controlling the pathogen epidemics over time. This combination did not prevent isolates virulent against the major gene from amplifying in the long term but the quantitative resistance significantly delayed for 5years the loss of effectiveness of the qualitative resistance and disease severity was maintained at a low level on the genotype with both types of resistance after the fungus population had adapted to the major gene. We also showed that diversity of AvrLm6 virulence alleles was comparable in isolates recovered after the recurrent selection on lines carrying either the major gene alone or in combination with quantitative resistance: a single repeat-induced point mutation and deletion events were observed in both situations. Breeding varieties which combine qualitative and quantitative resistance can effectively contribute to disease control by increasing the potential for durability of major resistance genes.

Identifying genetic diversity of avirulence genes in Leptosphaeria maculans using whole genome sequencing

Next generation sequencing technology allows rapid re-sequencing of individuals, as well as the discovery of single nucleotide polymorphisms (SNPs), for genomic diversity and evolutionary analyses. By sequencing two isolates of the fungal plant pathogen Leptosphaeria maculans, the causal agent of blackleg disease in Brassica crops, we have generated a resource of over 76 million sequence reads aligned to the reference genome. We identified over 21,000 SNPs with an overall SNP frequency of one SNP every 2,065 bp. Sequence validation of a selection of these SNPs in additional isolates collected throughout Australia indicates a high degree of polymorphism in the Australian population. In preliminary phylogenetic analysis, isolates from Western Australia clustered together and those collected from Brassica juncea stubble were identical. These SNPs provide a novel marker resource to study the genetic diversity of this pathogen. We demonstrate that re-sequencing provides a method of validating previously characterised SNPs and analysing differences in important genes, such as the disease related avirulence genes of L. maculans. Understanding the genetic characteristics of this devastating pathogen is vital in developing long-term solutions to managing blackleg disease in Brassica crops.

The dispensable chromosome of Leptosphaeria maculans shelters an effector gene conferring avirulence towards Brassica rapa

Phytopathogenic fungi frequently contain dispensable chromosomes, some of which contribute to host range or pathogenicity. In Leptosphaeria maculans, the stem canker agent of oilseed rape (Brassica napus), the minichromosome was previously suggested to be dispensable, without evidence for any role in pathogenicity. Using genetic and genomic approaches, we investigated the inheritance and molecular determinant of an L. maculans-Brassica rapa incompatible interaction. Single gene control of the resistance was found, while all markers located on the L. maculans minichromosome, absent in the virulent parental isolate, co-segregated with the avirulent phenotype. Only one candidate avirulence gene was identified on the minichromosome, validated by complementation experiments and termed AvrLm11. The minichromosome was frequently lost following meiosis, but the frequency of isolates lacking it remained stable in field populations sampled at a 10-yr time interval, despite a yearly sexual stage in the L. maculans life cycle. This work led to the cloning of a new 'lost in the middle of nowhere' avirulence gene of L. maculans, interacting with a B. rapa resistance gene termed Rlm11 and introgressed into B. napus. It demonstrated the dispensability of the L. maculans minichromosome and suggested that its loss generates a fitness deficit.

Microsatellite and single-nucleotide polymorphisms indicate recurrent transitions to asexuality in a microsporidian parasite

Assessing the mode of reproduction of microparasites remains a difficult task because direct evidence for sexual processes is often absent and the biological covariates of sex and asex are poorly known. Species with geographically divergent modes of reproduction offer the possibility to explore some of these covariates, for example, the influence of life-history traits, mode of transmission and life-cycle complexity. Here, we present a phylogeographical study of a microsporidian parasite, which allows us to relate population genetic structure and mode of reproduction to its geographically diverged life histories. We show that in microsporidians from the genus Hamiltosporidium, that use the cladoceran Daphnia as host, an epidemic population structure has evolved, most probably since the last Ice Age. We partially sequenced three housekeeping genes (alpha tubulin, beta tubulin and hsp70) and genotyped seven microsatellite loci in 51 Hamiltosporidium isolates sampled within Europe and the Middle East. We found two phylogenetically related asexual parasite lines, one each from Fennoscandia and Israel, which share the unique ability of being transmitted both vertically and horizontally from Daphnia to Daphnia. The sexual forms cannot transmit horizontally among Daphnia, but presumably have a complex life cycle with a second host species. In spite of the similarities between the two asexual lineages, a clustering analysis based on microsatellite polymorphisms shows that asexual Fennoscandian parasites do not share ancestry with any other Hamiltosporidium that we have sampled. Moreover, allele sequence divergence at the hsp70 locus is twice as large in Fennoscandian than in Israeli parasites. Our results indicate that asexual reproduction evolved twice independently, first in Fennoscandian and more recently in the Israeli parasites. We conclude that the independent origin of asexuality in these two populations is associated with the altered parasite mode of transmission and the underlying dynamics of host populations.

Experimental measures of pathogen competition and relative fitness

Competition among pathogen strains for limited host resources can have a profound effect on pathogen evolution. A better understanding of the principles and consequences of competition can be useful in designing more sustainable disease management strategies. The competitive ability and relative fitness of a pathogen strain are determined by its intrinsic biological properties, the resistance and heterogeneity of the corresponding host population, the population density and genetic relatedness of the competing strains, and the physical environment. Competitive ability can be inferred indirectly from fitness components, such as basic reproduction rate or transmission rate. However, pathogen strains that exhibit higher fitness components when they infect a host alone may not exhibit a competitive advantage when they co-infect the same host. The most comprehensive measures of competitive ability and relative fitness come from calculating selection coefficients in a mixed infection in a field setting. Mark-release-recapture experiments can be used to estimate fitness costs associated with unnecessary virulence and fungicide resistance.