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

Large-scale population survey of Leptosphaeria maculans in France highlights both on-going breakdowns and potentially effective resistance genes in oilseed rape

BACKGROUND

Leptosphaeria maculans, the cause of stem canker of oilseed rape, develops gene-for-gene interactions with its host and shows a high evolutionary potential to 'break down' novel resistance genes (R, Rlm) deployed in cultivars over large areas. For optimal management of R genes, updated knowledge of the population structure of the pathogen is needed. In France, large-scale surveys have been done at 10-year intervals since 2000. Here we report the characterization of a large L. maculans population collected in France in 2019-2020.

RESULTS

A total of 844 isolates were collected from 11 sites in ten French departments and were phenotyped for their virulence against nine Brassica napus R genes. All isolates were virulent toward Rlm2 and Rlm9. Very few isolates were avirulent on Rlm1 (1.8%) and Rlm4 (0.6%). Avirulent isolates toward Rlm7 ('AvrLm7') varied from 67% to 11.3%, depending on the site sampled, illustrating the ongoing breakdown of Rlm7. The decrease of AvrLm7 isolates (29.2% at the national level) compared to the 2010 survey (96.5%) was accompanied by an increase of avirulent isolates on Rlm3 (0% in 2010; 54% in 2019-2020). However, virulent isolates on both Rlm3 and Rlm7, previously rarely detected, were found in all sites with a frequency of 17.3%. Finally, most or all isolates were avirulent on Rlm11 (96.1%), LepR2 (RlmS, 99.8%), and Rlm6 (100%), suggesting these three genes still effectively control the disease.

CONCLUSION

These data will help guide strategies for breeding and deploying resistant oilseed rape varieties against L. maculans in France. © 2023 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Regulation of effector gene expression as concerted waves in Leptosphaeria maculans: a two-player game

Effector genes, encoding molecules involved in disease establishment, are concertedly expressed throughout the lifecycle of plant-pathogenic fungi. However, little is known about how effector gene expression is regulated. Since many effector genes are located in repeat-rich regions, the role of chromatin remodeling in their regulation was recently investigated, notably establishing that the repressive histone modification H3K9me3, deposited by KMT1, was involved in several fungal species including Leptosphaeria maculans. Nevertheless, previous data suggest that a second regulatory layer, probably involving a specific transcription factor (TF), might be required. In L. maculans, a Dothideomycete causing stem canker of oilseed rape, we identified the ortholog of Pf2, a TF belonging to the Zn2Cys6 fungal-specific family, and described as essential for pathogenicity and effector gene expression. We investigated its role together with KMT1, by inactivating and over-expressing LmPf2 in a wild-type strain and a ∆kmt1 mutant. Functional analyses of the corresponding transformants highlighted an essential role of LmPf2 in the establishment of pathogenesis and we found a major effect of LmPf2 on the induction of effector gene expression once KMT1 repression is lifted. Our results show, for the first time, a dual control of effector gene expression.

Population genomic analyses suggest recent dispersal events of the pathogen Cercospora zeina into East and Southern African maize cropping systems

A serious factor hampering global maize production is gray leaf spot disease. Cercospora zeina is one of the causative pathogens, but population genomics analysis of C. zeina is lacking. We conducted whole-genome Illumina sequencing of a representative set of 30 C. zeina isolates from Kenya and Uganda (East Africa) and Zambia, Zimbabwe, and South Africa (Southern Africa). Selection of the diverse set was based on microsatellite data from a larger collection of the pathogen. Pangenome analysis of the C. zeina isolates was done by (1) de novo assembly of the reads with SPAdes, (2) annotation with BRAKER, and (3) protein clustering with OrthoFinder. A published long-read assembly of C. zeina (CMW25467) from Zambia was included and annotated using the same pipeline. This analysis revealed 790 non-shared accessory and 10,677 shared core orthogroups (genes) between the 31 isolates. Accessory gene content was largely shared between isolates from all countries, with a few genes unique to populations from Southern Africa (32) or East Africa (6). There was a significantly higher proportion of effector genes in the accessory secretome (44%) compared to the core secretome (24%). PCA, ADMIXTURE, and phylogenetic analysis using a neighbor-net network indicated a population structure with a geographical subdivision between the East African isolates and the Southern African isolates, although gene flow was also evident. The small pangenome and partial population differentiation indicated recent dispersal of C. zeina into Africa, possibly from 2 regional founder populations, followed by recurrent gene flow owing to widespread maize production across sub-Saharan Africa.

The neighbouring genes AvrLm10A and AvrLm10B are part of a large multigene family of cooperating effector genes conserved in Dothideomycetes and Sordariomycetes

Fungal effectors (small-secreted proteins) have long been considered as species or even subpopulation-specific. The increasing availability of high-quality fungal genomes and annotations has allowed the identification of trans-species or trans-genera families of effectors. Two avirulence effectors, AvrLm10A and AvrLm10B, of Leptosphaeria maculans, the fungus causing stem canker of oilseed rape, are members of such a large family of effectors. AvrLm10A and AvrLm10B are neighbouring genes, organized in divergent transcriptional orientation. Sequence searches within the L. maculans genome showed that AvrLm10A/AvrLm10B belong to a multigene family comprising five pairs of genes with a similar tail-to-tail organization. The two genes, in a pair, always had the same expression pattern and two expression profiles were distinguished, associated with the biotrophic colonization of cotyledons and/or petioles and stems. Of the two protein pairs further investigated, AvrLm10A_like1/AvrLm10B_like1 and AvrLm10A_like2/AvrLm10B_like2, the second one had the ability to physically interact, similarly to what was previously described for the AvrLm10A/AvrLm10B pair, and cross-interactions were also detected for two pairs. AvrLm10A homologues were identified in more than 30 Dothideomycete and Sordariomycete plant-pathogenic fungi. One of them, SIX5, is an effector from Fusarium oxysporum f. sp. lycopersici physically interacting with the avirulence effector Avr2. We found that AvrLm10A/SIX5 homologues were associated with at least eight distinct putative effector families, suggesting that AvrLm10A/SIX5 is able to cooperate with different effectors. These results point to a general role of the AvrLm10A/SIX5 proteins as "cooperating proteins", able to interact with diverse families of effectors whose encoding gene is co-regulated with the neighbouring AvrLm10A homologue.

Polymorphism of Avirulence Genes and Adaptation to Brassica Resistance Genes Is Gene-Dependent in the Phytopathogenic Fungus Leptosphaeria maculans

The fungal phytopathogen Leptosphaeria maculans, which causes stem canker (blackleg) of rapeseed (Brassica napus), is mainly controlled worldwide by genetic resistance, which includes major resistance genes (Rlm). This model is one of those for which the highest number of avirulence genes (AvrLm) has been cloned. In many systems, including the L. maculans-B. napus interaction, intense use of resistance genes exerts strong selection pressure on the corresponding avirulent isolates, and the fungi may rapidly escape resistance through various molecular events which modify the avirulence genes. In the literature, the study of polymorphism at avirulence loci is often focused on single genes under selection pressure. In this study, we investigate allelic polymorphism at 11 avirulence loci in a French population of 89 L. maculans isolates collected on a trap cultivar in four geographic locations in the 2017-2018 cropping season. The corresponding Rlm genes have been (i) used for a long time, (ii) recently used, or (iii) unused in agricultural practice. The sequence data generated indicate an extreme diversity of situations. For example, genes submitted to an ancient selection may have either been deleted in populations (AvrLm1) or replaced by a single-nucleotide mutated virulent version (AvrLm2, AvrLm5-9). Genes that have never been under selection may either be nearly invariant (AvrLm6, AvrLm10A, AvrLm10B), exhibit rare deletions (AvrLm11, AvrLm14), or display a high diversity of alleles and isoforms (AvrLmS-Lep2). These data suggest that the evolutionary trajectory of avirulence/virulence alleles is gene-dependent and independent of selection pressure in L. maculans. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Comprehending the dynamism of B chromosomes in their journey towards becoming unselfish

Investigated for more than a century now, B chromosomes (Bs) research has come a long way from Bs being considered parasitic or neutral to becoming unselfish and bringing benefits to their hosts. B chromosomes exist as accessory chromosomes along with the standard A chromosomes (As) across eukaryotic taxa. Represented singly or in multiple copies, B chromosomes are largely heterochromatic but also contain euchromatic and organellar segments. Although B chromosomes are derived entities, they follow their species-specific evolutionary pattern. B chromosomes fail to pair with the standard chromosomes during meiosis and vary in their number, size, composition and structure across taxa and ensure their successful transmission through non-mendelian mechanisms like mitotic, pre-meiotic, meiotic or post-meiotic drives, unique non-disjunction, self-pairing or even imparting benefits to the host when they lack drive. B chromosomes have been associated with cellular processes like sex determination, pathogenicity, resistance to pathogens, phenotypic effects, and differential gene expression. With the advancements in B-omics research, novel insights have been gleaned on their functions, some of which have been associated with the regulation of gene expression of A chromosomes through increased expression of miRNAs or differential expression of transposable elements located on them. The next-generation sequencing and emerging technologies will further likely unravel the cellular, molecular and functional behaviour of these enigmatic entities. Amidst the extensive fluidity shown by B chromosomes in their structural and functional attributes, we perceive that the existence and survival of B chromosomes in the populations most likely seem to be a trade-off between the drive efficiency and adaptive significance versus their adverse effects on reproduction.

Chromosome-level analysis of the Colletotrichum graminicola genome reveals the unique characteristics of core and minichromosomes

The fungal pathogen Colletotrichum graminicola causes the anthracnose of maize (Zea mays) and is responsible for significant yield losses worldwide. The genome of C. graminicola was sequenced in 2012 using Sanger sequencing, 454 pyrosequencing, and an optical map to obtain an assembly of 13 pseudochromosomes. We re-sequenced the genome using a combination of short-read (Illumina) and long-read (PacBio) technologies to obtain a chromosome-level assembly. The new version of the genome sequence has 13 chromosomes with a total length of 57.43 Mb. We detected 66 (23.62 Mb) structural rearrangements in the new assembly with respect to the previous version, consisting of 61 (21.98 Mb) translocations, 1 (1.41 Mb) inversion, and 4 (221 Kb) duplications. We annotated the genome and obtained 15,118 predicted genes and 3,614 new gene models compared to the previous version of the assembly. We show that 25.88% of the new assembly is composed of repetitive DNA elements (13.68% more than the previous assembly version), which are mostly found in gene-sparse regions. We describe genomic compartmentalization consisting of repeat-rich and gene-poor regions vs. repeat-poor and gene-rich regions. A total of 1,140 secreted proteins were found mainly in repeat-rich regions. We also found that ~75% of the three smallest chromosomes (minichromosomes, between 730 and 551 Kb) are strongly affected by repeat-induced point mutation (RIP) compared with 28% of the larger chromosomes. The gene content of the minichromosomes (MCs) comprises 121 genes, of which 83.6% are hypothetical proteins with no predicted function, while the mean percentage of Chr1-Chr10 is 36.5%. No predicted secreted proteins are present in the MCs. Interestingly, only 2% of the genes in Chr11 have homologs in other strains of C. graminicola, while Chr12 and 13 have 58 and 57%, respectively, raising the question as to whether Chrs12 and 13 are dispensable. The core chromosomes (Chr1-Chr10) are very different with respect to the MCs (Chr11-Chr13) in terms of the content and sequence features. We hypothesize that the higher density of repetitive elements and RIPs in the MCs may be linked to the adaptation and/or host co-evolution of this pathogenic fungus.

Molecular Interactions Between Leptosphaeria maculans and Brassica Species

Canola is an important oilseed crop, providing food, feed, and fuel around the world. However, blackleg disease, caused by the ascomycete Leptosphaeria maculans, causes significant yield losses annually. With the recent advances in genomic technologies, the understanding of the Brassica napus-L. maculans interaction has rapidly increased, with numerous Avr and R genes cloned, setting this system up as a model organism for studying plant-pathogen associations. Although the B. napus-L. maculans interaction follows Flor's gene-for-gene hypothesis for qualitative resistance, it also puts some unique spins on the interaction. This review discusses the current status of the host-pathogen interaction and highlights some of the future gaps that need addressing moving forward.

Non-Mendelian transmission of accessory chromosomes in fungi

Non-Mendelian transmission has been reported for various genetic elements, ranging from small transposons to entire chromosomes. One prime example of such a transmission pattern are B chromosomes in plants and animals. Accessory chromosomes in fungi are similar to B chromosomes in showing presence/absence polymorphism and being non-essential. How these chromosomes are transmitted during meiosis is however poorly understood—despite their often high impact on the fitness of the host. For several fungal organisms, a non-Mendelian transmission or a mechanistically unique meiotic drive of accessory chromosomes have been reported. In this review, we provide an overview of the possible mechanisms that can cause the non-Mendelian transmission or meiotic drives of fungal accessory chromosomes. We compare processes responsible for the non-Mendelian transmission of accessory chromosomes for different fungal eukaryotes and discuss the structural traits of fungal accessory chromosomes affecting their meiotic transmission. We conclude that research on fungal accessory chromosomes, due to their small size, ease of sequencing, and epigenetic profiling, can complement the study of B chromosomes in deciphering factors that influence and regulate the non-Mendelian transmission of entire chromosomes.

Independent breakdown events of the Brassica napus Rlm7 resistance gene including via the off-target impact of a dual-specificity avirulence interaction

Protection of many crops is achieved through the use of genetic resistance. Leptosphaeria maculans, the causal agent of blackleg disease of Brassica napus, has emerged as a model for understanding gene-for-gene interactions that occur between plants and pathogens. Whilst many of the characterized avirulence effector genes interact with a single resistance gene in the host, the AvrLm4-7 avirulence gene is recognized by two resistance genes, Rlm4 and Rlm7. Here, we report the "breakdown" of the Rlm7 resistance gene in Australia, under two different field conditions. The first, and more typical, breakdown probably resulted from widescale use of Rlm7-containing cultivars whereby selection has led to an increase of individuals in the L. maculans population that have undergone repeat-induced point (RIP) mutations at the AvrLm4-7 locus. This has rendered the AvrLm4-7 gene ineffective and therefore these isolates have become virulent towards both Rlm4 and Rlm7. The second, more atypical, situation was the widescale use of Rlm4 cultivars. Whilst a single-nucleotide polymorphism is the more common mechanism of virulence towards Rlm4, in this field situation, RIP mutations have been selected leading to the breakdown of resistance for both Rlm4 and Rlm7. This is an example of a resistance gene being rendered ineffective without having grown cultivars with the corresponding resistance gene due to the dual specificity of the avirulence gene. These findings highlight the value of pathogen surveillance in the context of expanded knowledge about potential complexities for Avr-R interactions for the deployment of appropriate resistance gene strategies.

Two independent approaches converge to the cloning of a new Leptosphaeria maculans avirulence effector gene, AvrLmS-Lep2

Brassica napus (oilseed rape, canola) seedling resistance to Leptosphaeria maculans, the causal agent of blackleg (stem canker) disease, follows a gene-for-gene relationship. The avirulence genes AvrLmS and AvrLep2 were described to be perceived by the resistance genes RlmS and LepR2, respectively, present in B. napus 'Surpass 400'. Here we report cloning of AvrLmS and AvrLep2 using two independent methods. AvrLmS was cloned using combined in vitro crossing between avirulent and virulent isolates with sequencing of DNA bulks from avirulent or virulent progeny (bulked segregant sequencing). AvrLep2 was cloned using a biparental cross of avirulent and virulent L. maculans isolates and a classical map-based cloning approach. Taking these two approaches independently, we found that AvrLmS and AvrLep2 are the same gene. Complementation of virulent isolates with this gene confirmed its role in inducing resistance on Surpass 400, Topas-LepR2, and an RlmS-line. The gene, renamed AvrLmS-Lep2, encodes a small cysteine-rich protein of unknown function with an N-terminal secretory signal peptide, which is a common feature of the majority of effectors from extracellular fungal plant pathogens. The AvrLmS-Lep2/LepR2 interaction phenotype was found to vary from a typical hypersensitive response through intermediate resistance sometimes towards susceptibility, depending on the inoculation conditions. AvrLmS-Lep2 was nevertheless sufficient to significantly slow the systemic growth of the pathogen and reduce the stem lesion size on plant genotypes with LepR2, indicating the potential efficiency of this resistance to control the disease in the field.

Efficacy of Blackleg Major Resistance Genes in B. napus in Germany

Leptosphaeria maculans is one of the major pathogens of oilseed rape (B. napus). It causes blackleg disease, which accounts for significant yield losses worldwide. Using cultivars that harbor major resistance (R) genes is one of the most effective control methods. However, the efficacy of major R genes is related to the frequency of the corresponding avirulence (Avr) genes in a L. maculans population. In this paper, we report the Avr profiles of L. maculans populations and the ratio of its mating types in Northern and Central regions of Germany. Eleven Avr genes in five-hundred and seventy-four isolates were characterized either by applying cotyledon tests on a B. napus differential set or by amplifying avirulence gene-specific PCR markers. Fifty-two races were determined, among which the most dominant race was Avrlm6, -7, -11, AvrlepR1, -R2. Results showed that the resistance gene Rlm2 is 100% ineffective, some other major R genes such as Rlm1, Rlm3, Rlm4 and LepR3 are partially effective (with corresponding Avr frequencies ≤ 42%), while LepR1, LepR2, Rlm6, Rlm11 and Rlm7 can still provide relatively effective resistance in the German fields investigated (with corresponding Avr frequencies of 63–100%). Sexual reproduction is a factor that enhances the potential of L. maculans to evolve under selection pressure. Mating types of the L. maculans populations did not deviate from the ratio of 1:1 in the examined regions, indicating that sexual reproduction and ascospores play central roles in the L. maculans lifecycle. Overall, this study provides an important dataset for the establishment of a strategic plan to preserve the efficacies of major R genes in Germany by applying cultivar rotations of oilseed rape.

Leptosphaeria maculans-Brassica napus Battle: A Comparison of Incompatible vs. Compatible Interactions Using Dual RNASeq

Leptosphaeria maculans causes blackleg disease, which is one of the most destructive diseases of canola (Brassica napus L.). Due to the erosion of the current resistance in B. napus, it is pivotal to introduce new resistant genotypes to the growers. This study evaluated the potential of Rlm7 gene as resistance to its corresponding avirulence AvrLm7 gene is abundant. The Rlm7 line was inoculated with L. maculans isolate with AvrLm7; UMAvr7; and the CRISPR/Cas9 knockout AvrLm7 mutant, umavr7, of the same isolate to cause incompatible and compatible interactions, respectively. Dual RNA-seq showed differential gene expressions in both interactions. High expressions of virulence-related pathogen genes-CAZymes, merops, and effector proteins after 7-dpi in compatible interactions but not in incompatible interaction—confirmed that the pathogen was actively virulent only in compatible interactions. Salicyclic and jasmonic acid biosynthesis and signaling-related genes, defense-related PR1 gene (GSBRNA2T00150001001), and GSBRNA2T00068522001 in the NLR gene family were upregulated starting as early as 1- and 3-dpi in the incompatible interaction and the high upregulation of those genes after 7-dpi in compatible interactions confirmed the early recognition of the pathogen by the host and control it by early activation of host defense mechanisms in the incompatible interaction.

A chromosome-scale genome assembly of the tomato pathogen Cladosporium fulvum reveals a compartmentalized genome architecture and the presence of a dispensable chromosome

Cladosporium fulvum is a fungal pathogen that causes leaf mould of tomato. The reference genome of this pathogen was released in 2012 but its high repetitive DNA content prevented a contiguous assembly and further prohibited the analysis of its genome architecture. In this study, we combined third generation sequencing technology with the Hi-C chromatin conformation capture technique, to produce a high-quality and near complete genome assembly and gene annotation of a Race 5 isolate of C. fulvum. The resulting genome assembly contained 67.17 Mb organized into 14 chromosomes (Chr1-to-Chr14), all of which were assembled telomere-to-telomere. The smallest of the chromosomes, Chr14, is only 460 kb in size and contains 25 genes that all encode hypothetical proteins. Notably, PCR assays revealed that Chr14 was absent in 19 out of 24 isolates of a world-wide collection of C. fulvum, indicating that Chr14 is dispensable. Thus, C. fulvum is currently the second species of Capnodiales shown to harbour dispensable chromosomes. The genome of C. fulvum Race 5 is 49.7 % repetitive and contains 14 690 predicted genes with an estimated completeness of 98.9%, currently one of the highest among the Capnodiales. Genome structure analysis revealed a compartmentalized architecture composed of gene-dense and repeat-poor regions interspersed with gene-sparse and repeat-rich regions. Nearly 39.2 % of the C. fulvum Race 5 genome is affected by Repeat-Induced Point (RIP) mutations and evidence of RIP leakage toward non-repetitive regions was observed in all chromosomes, indicating the RIP plays an important role in the evolution of this pathogen. Finally, 345 genes encoding candidate effectors were identified in C. fulvum Race 5, with a significant enrichment of their location in gene-sparse regions, in accordance with the 'two-speed genome' model of evolution. Overall, the new reference genome of C. fulvum presents several notable features and is a valuable resource for studies in plant pathogens.

Evolution of B Chromosomes: From Dispensable Parasitic Chromosomes to Essential Genomic Players

B chromosomes represent additional chromosomes found in many eukaryotic organisms. Their origin is not completely understood but recent genomic studies suggest that they mostly arise through rearrangements and duplications from standard chromosomes. They can occur in single or multiple copies in a cell and are usually present only in a subset of individuals in the population. Because B chromosomes frequently show unstable inheritance, their maintenance in a population is often associated with meiotic drive or other mechanisms that increase the probability of their transmission to the next generation. For all these reasons, B chromosomes have been commonly considered to be nonessential, selfish, parasitic elements. Although it was originally believed that B chromosomes had little or no effect on an organism's biology and fitness, a growing number of studies have shown that B chromosomes can play a significant role in processes such as sex determination, pathogenicity and resistance to pathogens. In some cases, B chromosomes became an essential part of the genome, turning into new sex chromosomes or germline-restricted chromosomes with important roles in the organism's fertility. Here, we review such cases of "cellular domestication" of B chromosomes and show that B chromosomes can be important genomic players with significant evolutionary impact.

A new avirulence gene of Leptosphaeria maculans, AvrLm14, identifies a resistance source in American broccoli (Brassica oleracea) genotypes

In many cultivated crops, sources of resistance to diseases are sparse and rely on introgression from wild relatives. Agricultural crops often are allopolyploids resulting from interspecific crosses between related species, which are sources of diversity for resistance genes. This is the case for Brassica napus (oilseed rape, canola), an interspecific hybrid between Brassica rapa (turnip) and Brassica oleracea (cabbage). B. napus has a narrow genetic basis and few effective resistance genes against stem canker (blackleg) disease, caused by the fungus Leptosphaeria maculans, are currently available. B. rapa diversity has proven to be a valuable source of resistance (Rlm, LepR) genes, while B. oleracea genotypes were mostly considered susceptible. Here we identified a new resistance source in B. oleracea genotypes from America, potentially effective against French L. maculans isolates under both controlled and field conditions. Genetic analysis of fungal avirulence and subsequent cloning and validation identified a new avirulence gene termed AvrLm14 and suggested a typical gene-for-gene interaction between AvrLm14 and the postulated Rlm14 gene. AvrLm14 shares all the usual characteristics of L. maculans avirulence genes: it is hosted in a genomic region enriched in transposable elements and heterochromatin marks H3K9me3, its expression is repressed during vegetative growth but shows a strong overexpression 5-9 days following cotyledon infection, and it encodes a small secreted protein enriched in cysteine residues with few matches in databases. Similar to the previously cloned AvrLm10-A, AvrLm14 contributes to reduce lesion size on susceptible cotyledons, pointing to a complex interplay between effectors promoting or reducing lesion development.

Genetic diversity and population structure of Leptosphaeria maculans isolates in Western Canada

Leptosphaeria maculans is a serious concern for canola production worldwide. For effective disease management, knowledge of the pathogen's genetic variability and population structure is a prerequisite. In this study, whole-genome sequencing was performed for 162 of 1590 L. maculans isolates collected in the years 2007-2008 and 2012-2014 in Western Canada. DNA variants in genome-wide and specific regions including avirulence (Avr) genes were characterized. A total of 31,870 high-quality polymorphic DNA variants were used to study L. maculans genetic diversity and population structure. Cluster analysis showed that 150 isolates were clustered into 2 main groups and 4 subgroups by DNA variants located in either Avr or small secreted protein-encoding genes and into 2 main groups and 6 subgroups by genome-wide variants. The analysis of nucleotide diversity and differentiation also confirmed genetic variation within a population and among populations. Principal component analysis with genome-wide variants showed that the isolates collected in 2012-2014 were more genetically diverse than those collected in 2007-2008. Population structure analysis discovered three distinct sub-populations. Although isolates from Saskatchewan and Alberta were of similar genetic composition, Manitoba isolates were highly diverse. Genome-wide association study detected DNA variants in genes AvrLm4-7, Lema_T86300, and Lema_T86310 associated with the years of collection.

Epigenetic modifications affect the rate of spontaneous mutations in a pathogenic fungus

Mutations are the source of genetic variation and the substrate for evolution. Genome-wide mutation rates appear to be affected by selection and are probably adaptive. Mutation rates are also known to vary along genomes, possibly in response to epigenetic modifications, but causality is only assumed. In this study we determine the direct impact of epigenetic modifications and temperature stress on mitotic mutation rates in a fungal pathogen using a mutation accumulation approach. Deletion mutants lacking epigenetic modifications confirm that histone mark H3K27me3 increases whereas H3K9me3 decreases the mutation rate. Furthermore, cytosine methylation in transposable elements (TE) increases the mutation rate 15-fold resulting in significantly less TE mobilization. Also accessory chromosomes have significantly higher mutation rates. Finally, we find that temperature stress substantially elevates the mutation rate. Taken together, we find that epigenetic modifications and environmental conditions modify the rate and the location of spontaneous mutations in the genome and alter its evolutionary trajectory.

Remote homology clustering identifies lowly conserved families of effector proteins in plant-pathogenic fungi

Plant diseases caused by fungal pathogens are typically initiated by molecular interactions between 'effector' molecules released by a pathogen and receptor molecules on or within the plant host cell. In many cases these effector-receptor interactions directly determine host resistance or susceptibility. The search for fungal effector proteins is a developing area in fungal-plant pathology, with more than 165 distinct confirmed fungal effector proteins in the public domain. For a small number of these, novel effectors can be rapidly discovered across multiple fungal species through the identification of known effector homologues. However, many have no detectable homology by standard sequence-based search methods. This study employs a novel comparison method (RemEff) that is capable of identifying protein families with greater sensitivity than traditional homology-inference methods, leveraging a growing pool of confirmed fungal effector data to enable the prediction of novel fungal effector candidates by protein family association. Resources relating to the RemEff method and data used in this study are available from https://figshare.com/projects/Effector_protein_remote_homology/87965.

A Modified Meiotic Recombination in Brassica napus Largely Improves Its Breeding Efficiency

Simple Summary

The selection of varieties more resilient to disease and climate change requires generating new genetic diversity for breeding. The main mechanism for reshuffling genetic information is through the recombination of chromosomes during meiosis. We showed in oilseed rape (Brassica napus, AACC, 2n = 4x = 38), which is a natural hybrid formed from a cross between turnip (B. rapa, AA, 2n = 2x = 20) and cabbage (B. oleracea, CC, 2n = 2x = 18), that there is significantly more crossovers occurring along the entire A chromosomes in allotriploid AAC (crossbetween B. napus and B. rapa) than in diploid AA or allotetraploid AACC hybrids. We demonstrated that these allotriploid AAC hybrids are highly efficient to introduce new variability within oilseed rape varieties, notably by enabling the introduction of small genomic regions carrying genes controlling agronomically interesting traits.

Abstract

Meiotic recombination is the main tool used by breeders to generate biodiversity, allowing genetic reshuffling at each generation. It enables the accumulation of favorable alleles while purging deleterious mutations. However, this mechanism is highly regulated with the formation of one to rarely more than three crossovers, which are not randomly distributed. In this study, we showed that it is possible to modify these controls in oilseed rape (Brassica napus, AACC, 2n = 4x = 38) and that it is linked to AAC allotriploidy and not to polyploidy per se. To that purpose, we compared the frequency and the distribution of crossovers along A chromosomes from hybrids carrying exactly the same A nucleotide sequence, but presenting three different ploidy levels: AA, AAC and AACC. Genetic maps established with 202 SNPs anchored on reference genomes revealed that the crossover rate is 3.6-fold higher in the AAC allotriploid hybrids compared to AA and AACC hybrids. Using a higher SNP density, we demonstrated that smaller and numerous introgressions of B. rapa were present in AAC hybrids compared to AACC allotetraploid hybrids, with 7.6 Mb vs. 16.9 Mb on average and 21 B. rapa regions per plant vs. nine regions, respectively. Therefore, this boost of recombination is highly efficient to reduce the size of QTL carried in cold regions of the oilseed rape genome, as exemplified here for a QTL conferring blackleg resistance.

Validating the Strategic Deployment of Blackleg Resistance Gene Groups in Commercial Canola Fields on the Canadian Prairies

Blackleg, caused by the fungal pathogen Leptosphaeria maculans, is a serious threat to canola (Brassica napus L.) production in western Canada. Crop scouting and extended crop rotation, along with the use of effective genetic resistance, have been key management practices available to mitigate the impact of the disease. In recent years, new pathogen races have reduced the effectiveness of some of the resistant cultivars deployed. Strategic deployment and rotation of major resistance (R) genes in cultivars have been used in France and Australia to help increase the longevity of blackleg resistance. Canada also introduced a grouping system in 2017 to identify blackleg R genes in canola cultivars. The main objective of this study was to examine and validate the concept of R gene deployment through monitoring the avirulence (Avr) profile of L. maculans population and disease levels in commercial canola fields within the Canadian prairies. Blackleg disease incidence and severity was collected from 146 cultivars from 53 sites across Manitoba, Saskatchewan, and Alberta in 2018 and 2019, and the results varied significantly between gene groups, which is likely influenced by the pathogen population. Isolates collected from spring and fall stubble residues were examined for the presence of Avr alleles AvrLm1, AvrLm2, AvrLm3, AvrLm4, AvrLm5, AvrLm6, AvrLm7, AvrLm9, AvrLm10, AvrLm11, AvrLepR1, AvrLepR2, AvrLep3, and AvrLmS using a set of differential host genotypes carrying known resistance genes or PCR-based markers. The Simpson's evenness index was very low, due to two dominant L. maculans races (AvrLm2-4-5-6-7-10-11 and AvrLm2-5-6-7-10-11) representing 49% of the population, but diversity of the population was high from the 35 L. maculans races isolated in Manitoba. AvrLm6 and AvrLm11 were found in all 254 L. maculans isolates collected in Manitoba. Knowledge of the blackleg disease levels in relation to the R genes deployed, along with the L. maculans Avr profile, helps to measure the effectiveness of genetic resistance.

Genome-wide mapping of histone modifications during axenic growth in two species of Leptosphaeria maculans showing contrasting genomic organization

Leptosphaeria maculans 'brassicae' (Lmb) and Leptosphaeria maculans 'lepidii' (Lml) are closely related phytopathogenic species that exhibit a large macrosynteny but contrasting genome structure. Lmb has more than 30% of repeats clustered in large repeat-rich regions, while the Lml genome has only a small amount of evenly distributed repeats. Repeat-rich regions of Lmb are enriched in effector genes, expressed during plant infection. The distinct genome structures of Lmb and Lml provide an excellent model for comparing the organization of pathogenicity genes in relation to the chromatin landscape in two closely related phytopathogenic fungi. Here, we performed chromatin immunoprecipitation (ChIP) during axenic culture, targeting histone modifications typical for heterochromatin or euchromatin, combined with transcriptomic analysis to analyze the influence of chromatin organization on gene expression. In both species, we found that facultative heterochromatin is enriched with genes lacking functional annotation, including numerous effector and species-specific genes. Notably, orthologous genes located in H3K27me3 domains are enriched with effector genes. Compared to other fungal species, including Lml, Lmb is distinct in having large H3K9me3 domains associated with repeat-rich regions that contain numerous species-specific effector genes. Discovery of these two distinctive heterochromatin landscapes now raises questions about their involvement in the regulation of pathogenicity, the dynamics of these domains during plant infection and the selective advantage to the fungus to host effector genes in H3K9me3 or H3K27me3 domains.

The Rlm13 Gene, a New Player of Brassica napus-Leptosphaeria maculans Interaction Maps on Chromosome C03 in Canola

Canola exhibits an extensive genetic variation for resistance to blackleg disease, caused by the fungal pathogen Leptosphaeria maculans. Despite the identification of several Avr effectors and R (race-specific) genes, specific interactions between Avr-R genes are not yet fully understood in the Brassica napus-L. maculans pathosystem. In this study, we investigated the genetic basis of resistance in an F_{2 : 3} population derived from Australian canola varieties CB-Telfer (Rlm4)/ATR-Cobbler (Rlm4) using a single-spore isolate of L. maculans, PHW1223. A genetic linkage map of the CB-Telfer/ATR-Cobbler population was constructed using 7,932 genotyping-by-sequencing-based DArTseq markers and subsequently utilized for linkage and haplotype analyses. Genetic linkage between DArTseq markers and resistance to PHW1223 isolate was also validated using the B. napus 60K Illumina Infinium array. Our results revealed that a major locus for resistance, designated as Rlm13, maps on chromosome C03. To date, no R gene for resistance to blackleg has been reported on the C subgenome in B. napus. Twenty-four candidate R genes were predicted to reside within the quantitative trait locus (QTL) region. We further resequenced both the parental lines of the mapping population (CB-Telfer and ATR-Cobbler, > 80 × coverage) and identified several structural sequence variants in the form of single-nucleotide polymorphisms (SNPs), insertions/deletions (InDels), and presence/absence variations (PAVs) near Rlm13. Comparative mapping revealed that Rlm13 is located within the homoeologous A03/C03 region in ancestral karyotype block "R" of Brassicaceae. Our results provide a "target" for further understanding the Avr-Rlm13 gene interaction as well as a valuable tool for increasing resistance to blackleg in canola germplasm.

Leptosphaeria maculans Isolates Reveal Their Allele Frequency in Western Canada

Blackleg, caused by Leptosphaeria maculans, is a major disease of canola in Canada, Australia, and Europe. For effective deployment of resistant varieties and disease management, it is crucial to understand the population structure of L. maculans. In this study, we analyzed L. maculans isolates from commercial fields in western Canada from 2014 to 2016 for the presence and frequency of avirulence (Avr) genes. A total of 1,584 isolates were examined for the presence of Avr genes AvrLm1, AvrLm2, AvrLm3, AvrLm4, AvrLm6, AvrLm7, AvrLm9, AvrLepR1, AvrLepR2, and AvrLmS via a set of differential host genotypes carrying known resistance genes and a PCR assay. Several Avr genes showed a higher frequency in the pathogen population, such as AvrLm6 and AvrLm7, which were present in >90% of isolates, whereas AvrLm3, AvrLm9, and AvrLepR2 showed frequencies of <10%. A total of 189 races (different combinations of Avr genes) were detected, with Avr-2-4-6-7-S, Avr-1-4-6-7, and Avr-2-4-6-7 as the three predominant races. When the effect of crop rotation was assessed, only a 3-year rotation showed a significantly higher frequency of AvrLm2 relative to shorter rotations. This study provides the information for producers to select effective canola varieties for blackleg management and for breeders to deploy new R genes in disease resistance breeding in western Canada.

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.

Genomic rearrangements generate hypervariable mini-chromosomes in host-specific isolates of the blast fungus

Supernumerary mini-chromosomes–a unique type of genomic structural variation–have been implicated in the emergence of virulence traits in plant pathogenic fungi. However, the mechanisms that facilitate the emergence and maintenance of mini-chromosomes across fungi remain poorly understood. In the blast fungus Magnaporthe oryzae (Syn. Pyricularia oryzae), mini-chromosomes have been first described in the early 1990s but, until very recently, have been overlooked in genomic studies. Here we investigated structural variation in four isolates of the blast fungus M. oryzae from different grass hosts and analyzed the sequences of mini-chromosomes in the rice, foxtail millet and goosegrass isolates. The mini-chromosomes of these isolates turned out to be highly diverse with distinct sequence composition. They are enriched in repetitive elements and have lower gene density than core-chromosomes. We identified several virulence-related genes in the mini-chromosome of the rice isolate, including the virulence-related polyketide synthase Ace1 and two variants of the effector gene AVR-Pik. Macrosynteny analyses around these loci revealed structural rearrangements, including inter-chromosomal translocations between core- and mini-chromosomes. Our findings provide evidence that mini-chromosomes emerge from structural rearrangements and segmental duplication of core-chromosomes and might contribute to adaptive evolution of the blast fungus.

The genomes of plant pathogens often exhibit an architecture that facilitates high rates of dynamic rearrangements and genetic diversification in virulence associated regions. These regions, which tend to be gene sparse and repeat rich, are thought to serve as a cradle for adaptive evolution. Supernumerary chromosomes, i.e. chromosomes that are only present in some but not all individuals of a species, are a special type of structural variation that have been observed in plants, animals, and fungi. Here we identified and studied supernumerary mini-chromosomes in the blast fungus Magnaporthe oryzae, a pathogen that causes some of the most destructive plant diseases. We found that rice, foxtail millet and goosegrass isolates of this pathogen contain mini-chromosomes with distinct sequence composition. All mini-chromosomes are rich in repetitive genetic elements and have lower gene densities than core-chromosomes. Further, we identified virulence-related genes on the mini-chromosome of the rice isolate. We observed large-scale genomic rearrangements around these loci, indicative of a role of mini-chromosomes in facilitating genome dynamics. Taken together, our results indicate that mini-chromosomes contribute to genome rearrangements and possibly adaptive evolution of the blast fungus.

Recent Findings Unravel Genes and Genetic Factors Underlying Leptosphaeria maculans Resistance in Brassica napus and Its Relatives

Among the Brassica oilseeds, canola (Brassica napus) is the most economically significant globally. However, its production can be limited by blackleg disease, caused by the fungal pathogen Lepstosphaeria maculans. The deployment of resistance genes has been implemented as one of the key strategies to manage the disease. Genetic resistance against blackleg comes in two forms: qualitative resistance, controlled by a single, major resistance gene (R gene), and quantitative resistance (QR), controlled by numerous, small effect loci. R-gene-mediated blackleg resistance has been extensively studied, wherein several genomic regions harbouring R genes against L. maculans have been identified and three of these genes were cloned. These studies advance our understanding of the mechanism of R gene and pathogen avirulence (Avr) gene interaction. Notably, these studies revealed a more complex interaction than originally thought. Advances in genomics help unravel these complexities, providing insights into the genes and genetic factors towards improving blackleg resistance. Here, we aim to discuss the existing R-gene-mediated resistance, make a summary of candidate R genes against the disease, and emphasise the role of players involved in the pathogenicity and resistance. The comprehensive result will allow breeders to improve resistance to L. maculans, thereby increasing yield.

Impact of a resistance gene against a fungal pathogen on the plant host residue microbiome: The case of the Leptosphaeria maculans-Brassica napus pathosystem

Oilseed rape residues are a crucial determinant of stem canker epidemiology as they support the sexual reproduction of the fungal pathogen Leptosphaeria maculans. The aim of this study was to characterize the impact of a resistance gene against L. maculans infection on residue microbial communities and to identify microorganisms interacting with this pathogen during residue degradation. We used near-isogenic lines to obtain healthy and infected host plants. The microbiome associated with the two types of plant residues was characterized by metabarcoding. A combination of linear discriminant analysis and ecological network analysis was used to compare the microbial communities and to identify microorganisms interacting with L. maculans. Fungal community structure differed between the two lines at harvest, but not subsequently, suggesting that the presence/absence of the resistance gene influences the microbiome at the base of the stem whilst the plant is alive, but that this does not necessarily lead to differential colonization of the residues by fungi. Direct interactions with other members of the community involved many fungal and bacterial amplicon sequence variants (ASVs). L. maculans appeared to play a minor role in networks, whereas one ASV affiliated to Plenodomus biglobosus (synonym Leptosphaeria biglobosa) from the Leptosphaeria species complex may be considered a keystone taxon in the networks at harvest. This approach could be used to identify and promote microorganisms with beneficial effects against residue-borne pathogens and, more broadly, to decipher the complex interactions between multispecies pathosystems and other microbial components in crop residues.

Development of a specific marker for detection of a functional AvrLm9 allele and validating the interaction between AvrLm7 and AvrLm9 in Leptosphaeria maculans

Blackleg, which is caused by the fungus Leptosphaeria maculans (L. maculans), is a major disease of canola in western Canada and worldwide. Long-term use of one source of resistance could cause the breakdown of its effectiveness. Therefore, appropriate use of R genes is very important, and knowledge about the distribution of avirulence genes is a prerequisite for effectively deploying resistance. Of the 14 avirulence genes identified in L. maculans, AvrLm5 and AvrLm9 were recognized as the two alleles of the same gene based on two single nucleotide polymorphisms, C⁸⁵T and G¹⁶⁴A/C. In this study, a specific marker was developed to identify AvrLm5 and AvrLm9 based on two single nucleotide polymorphisms, C⁸⁵T and G¹⁶⁴A/C, which are responsible for the function of AvrLm9. The specific marker can be used to discriminate the AvrLm9 from avrLm9 accurately in L. maculans isolates, which is consistent with inoculation tests in isolates without AvrLm4-7. This specific marker was used to screen 1229 isolates collected from fields in the years 2014 through 2016 in Manitoba. From 68 to 84% of the isolates were found to contain the AvrLm9 allele; while 4-7% of them were avirulent on the variety Goéland with Rlm9 loci. Furthermore, no isolates having both AvrLm9 and AvrLm7 were detected using a cotyledon test, while 67% to 84% of isolates contained both avirulence genes via PCR detection, implying suppression of AvrLm9 by AvrLm7. In addition, avirulence gene profiles of the other 10 avirulence alleles were examined with the 1229 isolates using cotyledon tests or PCR amplifications. Taken together, this research enables the fast identification of AvrLm5/9, provides the Avr genes' landscape of western Canada and elaborates the relationship between AvrLm9 and AvrLm7 using isolates from grower fields.

Genome compartmentalization predates species divergence in the plant pathogen genus Zymoseptoria

Background

Antagonistic co-evolution can drive rapid adaptation in pathogens and shape genome architecture. Comparative genome analyses of several fungal pathogens revealed highly variable genomes, for many species characterized by specific repeat-rich genome compartments with exceptionally high sequence variability. Dynamic genome structure may enable fast adaptation to host genetics. The wheat pathogen Zymoseptoria tritici with its highly variable genome, has emerged as a model organism to study genome evolution of plant pathogens. Here, we compared genomes of Z. tritici isolates and of sister species infecting wild grasses to address the evolution of genome composition and structure.

Results

Using long-read technology, we sequenced and assembled genomes of Z. ardabiliae, Z. brevis, Z. pseudotritici and Z. passerinii, together with two isolates of Z. tritici. We report a high extent of genome collinearity among Zymoseptoria species and high conservation of genomic, transcriptomic and epigenomic signatures of compartmentalization. We identify high gene content variability both within and between species. In addition, such variability is mainly limited to the accessory chromosomes and accessory compartments. Despite strong host specificity and non-overlapping host-range between species, predicted effectors are mainly shared among Zymoseptoria species, yet exhibiting a high level of presence-absence polymorphism within Z. tritici. Using in planta transcriptomic data from Z. tritici, we suggest different roles for the shared orthologs and for the accessory genes during infection of their hosts.

Conclusion

Despite previous reports of high genomic plasticity in Z. tritici, we describe here a high level of conservation in genomic, epigenomic and transcriptomic composition and structure across the genus Zymoseptoria. The compartmentalized genome allows the maintenance of a functional core genome co-occurring with a highly variable accessory genome.

Host-specificity factors in plant pathogenic fungi

Fortunately, no fungus can cause disease on all plant species, and although some plant-pathogenic fungi have quite a broad host range, most are highly limited in the range of plant species or even cultivars that they cause disease in. The mechanisms of host specificity have been extensively studied in many plant-pathogenic fungi, especially in fungal pathogens causing disease on economically important crops. Specifically, genes involved in host specificity have been identified during the last few decades. In this overview, we describe and discuss these host-specificity genes. These genes encode avirulence (Avr) proteins, proteinaceous host-specific toxins or secondary metabolites. We discuss the genomic context of these genes, their expression, polymorphism, horizontal transfer and involvement in pathogenesis.

Generation and Characterization of a Virulent Leptosphaeria maculans Isolate Carrying a Mutated AvrLm7 Gene Using the CRISPR/Cas9 System

Blackleg, caused by the fungal pathogen Leptosphaeria maculans, is the most important disease affecting canola (Brassica napus) crops worldwide. We employed the clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated (Cas) system to generate the mutant isolate umavr7 from a point mutation of the AvrLm7 coding region in a L. maculans isolate (UMAvr7). Reverse transcription PCR and transcriptome data confirmed that the AvrLm7 gene was knocked out in the mutant isolate. Pathogenicity tests indicated that umavr7 can cause large lesions on a set of Brassica differential genotypes that express different resistance (R) genes. Comparative pathogenicity tests between UMAvr7 (wild type) and umavr7 on the corresponding B. napus genotype 01-23-2-1 (with Rlm7) showed that umavr7 is a mutant isolate, producing large gray/green lesions on cotyledons. The pathogenicity of the mutant isolate was shifted from avirulent to virulent on the B. napus Rlm7 genotype. Therefore, this mutant is virulence on the identified resistant genes to blackleg disease in B. napus genotypes. Superoxide accumulated differently in cotyledons in response to infection with UMAvr7 and umavr7, especially in resistant B. napus genotype 01-23-2-1. Resistance/susceptibility was further evaluated on 123 B. napus genotypes with the mutant isolate, umavr7. Only 6 of the 123 genotypes showed resistance to umavr7. The identification of these six resistant B. napus genotypes will lead to further studies on the development of blackleg disease resistance through breeding and the identification of novel R genes.

Identification of resistance loci in Chinese and Canadian canola/rapeseed varieties against Leptosphaeria maculans based on genome-wide association studies

Background

The fungal pathogen Leptosphaeria maculans (Lm). causes blackleg disease on canola/rapeseed in many parts of the world. It is important to use resistant cultivars to manage the disease and minimize yield losses. In this study, twenty-two Lm isolates were used to identify resistance genes in a collection of 243 canola/rapeseed (Brassica napus L.) accessions from Canada and China. These Lm isolates carry different compliments of avirulence genes, and the investigation was based on a genome-wide association study (GWAS) and genotype-by-sequencing (GBS).

Results

Using the CROP-SNP pipeline, a total of 81,471 variants, including 78,632 SNPs and 2839 InDels, were identified. The GWAS was performed using TASSEL 5.0 with GLM + Q model. Thirty-two and 13 SNPs were identified from the Canadian and Chinese accessions, respectively, tightly associated with blackleg resistance with P values < 1 × 10^− 4. These SNP loci were distributed on chromosomes A03, A05, A08, A09, C01, C04, C05, and C07, with the majority of them on A08 followed by A09 and A03. The significant SNPs identified on A08 were all located in a 2010-kb region and associated with resistance to 12 of the 22 Lm isolates. Furthermore, 25 resistance gene analogues (RGAs) were identified in these regions, including two nucleotide binding site (NBS) domain proteins, fourteen RLKs, three RLPs and six TM-CCs. These RGAs can be the potential candidate genes for blackleg resistance.

Conclusion

This study provides insights into potentially new genomic regions for discovery of additional blackleg resistance genes. The identified regions associated with blackleg resistance in the germplasm collection may also contribute directly to the development of canola varieties with novel resistance genes against blackleg of canola.

An Avirulence Gene Cluster in the Wheat Stripe Rust Pathogen (Puccinia striiformis f. sp. tritici) Identified through Genetic Mapping and Whole-Genome Sequencing of a Sexual Population

Puccinia striiformis f. sp. tritici, the causal agent of wheat stripe (yellow) rust, is an obligate, biotrophic fungus. It was difficult to study the genetics of the pathogen due to the lack of sexual reproduction. The recent discovery of alternate hosts for P. striiformis f. sp. tritici makes it possible to study inheritance and map genes involved in its interaction with plant hosts. To identify avirulence (Avr) genes in P. striiformis f. sp. tritici, we developed a segregating population by selfing isolate 12-368 on barberry (Berberis vulgaris) plants under controlled conditions. The dikaryotic sexual population segregated for avirulent/virulent phenotypes on nine Yr single-gene lines. The parental and progeny isolates were whole-genome sequenced at >30× coverage using Illumina HiSeq PE150 technology. A total of 2,637 high-quality markers were discovered by mapping the whole-genome sequencing (WGS) reads to the reference genome of strain 93-210 and used to construct a genetic map, consisting of 41 linkage groups, spanning 7,715.0 centimorgans (cM) and covering 68 Mb of the reference genome. The recombination rate was estimated to be 1.81 ± 2.32 cM/10 kb. Quantitative trait locus analysis mapped six Avr gene loci to the genetic map, including an Avr cluster harboring four Avr genes, AvYr7, AvYr43, AvYr44, and AvYrExp2 Aligning the genetic map to the reference genome identified Avr candidates and narrowed them to a small genomic region (<200 kb). The discovery of the Avr gene cluster is useful for understanding pathogen evolution, and the identification of candidate genes is an important step toward cloning Avr genes for studying molecular mechanisms of pathogen-host interactions.IMPORTANCE Stripe rust is a destructive disease of wheat worldwide. Growing resistant cultivars is the most effective, easy-to-use, economical, and environmentally friendly strategy for the control of the disease. However, P. striiformis f. sp. tritici can produce new virulent races that may circumvent race-specific resistance. Therefore, understanding the genetic basis of the interactions between wheat genes for resistance and P. striiformis f. sp. tritici genes for avirulence is useful for improving cultivar resistance for more effective control of the disease. This study developed a high-quality map that facilitates genomic and genetic studies of important traits related to pathogen pathogenicity and adaptation to different environments and crop cultivars carrying different resistance genes. The information on avirulence/virulence genes identified in this study can be used for guiding breeding programs to select combinations of genes for developing new cultivars with effective resistance to mitigate this devastating disease.

Polycomb Repression without Bristles: Facultative Heterochromatin and Genome Stability in Fungi

Genome integrity is essential to maintain cellular function and viability. Consequently, genome instability is frequently associated with dysfunction in cells and associated with plant, animal, and human diseases. One consequence of relaxed genome maintenance that may be less appreciated is an increased potential for rapid adaptation to changing environments in all organisms. Here, we discuss evidence for the control and function of facultative heterochromatin, which is delineated by methylation of histone H3 lysine 27 (H3K27me) in many fungi. Aside from its relatively well understood role in transcriptional repression, accumulating evidence suggests that H3K27 methylation has an important role in controlling the balance between maintenance and generation of novelty in fungal genomes. We present a working model for a minimal repressive network mediated by H3K27 methylation in fungi and outline challenges for future research.

The Many Questions about Mini Chromosomes in Colletotrichum spp

Many fungal pathogens carry accessory regions in their genome, which are not required for vegetative fitness. Often, although not always, these regions occur as relatively small chromosomes in different species. Such mini chromosomes appear to be a typical feature of many filamentous plant pathogens. Since these regions often carry genes coding for effectors or toxin-producing enzymes, they may be directly related to virulence of the respective pathogen. In this review, we outline the situation of small accessory chromosomes in the genus Colletotrichum, which accounts for ecologically important plant diseases. We summarize which species carry accessory chromosomes, their gene content, and chromosomal makeup. We discuss the large variation in size and number even between different isolates of the same species, their potential roles in host range, and possible mechanisms for intra- and interspecies exchange of these interesting genetic elements.

A Comparative Genomic Analysis of the Barley Pathogen Pyrenophora teres f. teres Identifies Subtelomeric Regions as Drivers of Virulence

Pyrenophora teres f. teres causes net form net blotch of barley and is an economically important pathogen throughout the world. However, P. teres f. teres is lacking in the genomic resources necessary to characterize the mechanisms of virulence. Recently a high-quality reference genome was generated for P. teres f. teres isolate 0-1. Here, we present the reference quality sequence and annotation of four new isolates and we use the five available P. teres f. teres genomes for an in-depth comparison, resulting in the generation of hypotheses pertaining to the potential mechanisms and evolution of virulence. Comparative analyses were performed between all five P. teres f. teres genomes, examining genomic organization, structural variations, and core and accessory genomic content, specifically focusing on the genomic characterization of known virulence loci and the localization of genes predicted to encode secreted and effector proteins. We showed that 14 of 15 currently published virulence quantitative trait loci (QTL) span accessory genomic regions, consistent with these accessory regions being important drivers of host adaptation. Additionally, these accessory genomic regions were frequently found in subtelomeric regions of chromosomes, with 10 of the 14 accessory region QTL localizing to subtelomeric regions. Comparative analysis of the subtelomeric regions of P. teres f. teres chromosomes revealed translocation events in which homology was detected between nonhomologous chromosomes at a significantly higher rate than the rest of the genome. These results indicate that the subtelomeric accessory genomic compartments not only harbor most of the known virulence loci but, also, that these regions have the capacity to rapidly evolve.

LmCBP1, a secreted chitin-binding protein, is required for the pathogenicity of Leptosphaeria maculans on Brassica napus

Leptosphaeria maculans is the causal agent of blackleg disease on Brassica napus. Determining the underlying functions of genes required for pathogenesis is essential for understanding the infection process. A chitin-binding protein (LmCBP1) was discovered as a pathogenicity factor for the infection of B. napus by L. maculans through gene knockout using the CRISPR-Cas9 system. Chitin-binding activity was demonstrated through a chitin-protein binding assay. A secreted signal peptide was detected using a yeast secreted-signal peptide trap assay. An increased expression level during the infection stage was also observed, suggesting that LmCBP1 is a secreted protein. The knockout mutants showed decreased infection on B. napus, with reduced pathogenicity on ten cultivars with/without diverse R genes. The mutants were more sensitive to H₂O₂ compared to wild type L. maculans isolate JN3. This study provides evidence of the virulence of a novel chitin-binding protein LmCBP1 on B. napus through mutants created via the CRISPR-Cas9 system.

Repeat-Induced Point Mutations Drive Divergence between Fusarium circinatum and Its Close Relatives

The Repeat-Induced Point (RIP) mutation pathway is a fungal-specific genome defense mechanism that counteracts the deleterious effects of transposable elements. This pathway permanently mutates its target sequences by introducing cytosine to thymine transitions. We investigated the genome-wide occurrence of RIP in the pitch canker pathogen, Fusarium circinatum, and its close relatives in the Fusarium fujikuroi species complex (FFSC). Our results showed that the examined fungi all exhibited hallmarks of RIP, but that they differed in terms of the extent to which their genomes were affected by this pathway. RIP mutations constituted a large proportion of all the FFSC genomes, including both core and dispensable chromosomes, although the latter were generally more extensively affected by RIP. Large RIP-affected genomic regions were also much more gene sparse than the rest of the genome. Our data further showed that RIP-directed sequence diversification increased the variability between homologous regions of related species, and that RIP-affected regions can interfere with homologous recombination during meiosis, thereby contributing to post-mating segregation distortion. Taken together, these findings suggest that RIP can drive the independent divergence of chromosomes, alter chromosome architecture, and contribute to the divergence among F. circinatum and other members of this economically important group of fungi.

Effect of Elevated CO2 Concentration on the Disease Severity of Compatible and Incompatible Interactions of Brassica napus-Leptosphaeria maculans Pathosystem

Global warming by increased atmospheric CO2 concentration has been widely accepted. Yet, there has not been any consistent conclusion on the doubled CO2 concentration that in the future will affect plant disease incidence and severity. Blackleg disease, mainly caused by Leptosphaeria maculans, is a major disease on canola production globally. Brassica napus and L. maculans have a gene-for-gene interaction, which causes an incompatible reaction between canola plants carrying resistance genes and L. maculans isolates carrying corresponding avirulence genes. In this study, B. napus varieties and lines inoculated with different Leptosphaeria isolates were subjected to simulated growth conditions, namely, growth chambers with normal environments and with controlled CO2 concentrations of 400, 600, and 800 ppm. The results indicated that the elevated CO2 concentrations have no noticeable effect on the inferred phenotypes of the canola-blackleg interactions. However, the disease severity decreased in most of the B. napus-L. maculans interactions at extremely high CO2 concentration (800 ppm). The varied pathogenicity changes of the B. napus-L. maculans pathosystem under elevated CO2 concentrations at 400 or 600 ppm may be due to the genetic background or physiological differences in plants and pathogenicity differences in L. maculans isolates having different Avr gene profiles. The mechanisms by which elevated CO2 concentrations affect the B. napus-L. maculans pathosystem will help us understand how climate change will impact crops and diseases.

An evolutionary framework for host shifts - jumping ships for survival

Host jumping is a process by which pathogens settle in new host groups. It is a cornerstone in the evolution of pathogens, as it leads to pathogen diversification. It is unsurprising that host jumping is observed in facultative pathogens, as they can reproduce even if they kill their hosts. However, host jumps were thought to be rare in obligate biotrophic pathogens, but molecular phylogenetics has revealed that the opposite is true. Here, I review some concepts and recent findings and present several hypotheses on the matter. In short, pathogens evolve and diversify via host jumps, followed by radiation, specialisation and speciation. Host jumps are facilitated by, for example, effector innovations, stress, compatible pathogens and physiological similarities. Host jumping, subsequent establishment, and speciation takes place rapidly - within centuries and millennia rather than over millions of years. If pathogens are unable to evolve into neutral or mutualistic interactions with their hosts, they will eventually be removed from the host population, despite balancing trade-offs. Thus, generally, plant pathogens only survive in the course of evolution if they jump hosts. This is also reflected by the diversity patterns observed in many genera of plant pathogens, where it leads to a mosaic pattern of host groups over time, in which the original host group becomes increasingly obscure.

Effector gene reshuffling involves dispensable mini-chromosomes in the wheat blast fungus

Newly emerged wheat blast disease is a serious threat to global wheat production. Wheat blast is caused by a distinct, exceptionally diverse lineage of the fungus causing rice blast disease. Through sequencing a recent field isolate, we report a reference genome that includes seven core chromosomes and mini-chromosome sequences that harbor effector genes normally found on ends of core chromosomes in other strains. No mini-chromosomes were observed in an early field strain, and at least two from another isolate each contain different effector genes and core chromosome end sequences. The mini-chromosome is enriched in transposons occurring most frequently at core chromosome ends. Additionally, transposons in mini-chromosomes lack the characteristic signature for inactivation by repeat-induced point (RIP) mutation genome defenses. Our results, collectively, indicate that dispensable mini-chromosomes and core chromosomes undergo divergent evolutionary trajectories, and mini-chromosomes and core chromosome ends are coupled as a mobile, fast-evolving effector compartment in the wheat pathogen genome.

The emerging blast disease on wheat is proving even harder to control than the ancient, still-problematic rice blast disease. Potential wheat resistance identified using strains isolated soon after disease emergence are no longer effective in controlling recent aggressive field isolates from wheat in South America and South Asia. We construct a high-quality assembly of an aggressive, recently-isolated wheat blast fungal strain and the first assembled mini-chromosome genome sequence of wheat and rice blast pathogens. We report that recent wheat pathogens can contain one or two highly-variable dispensable mini-chromosomes, each with an amalgamation of fungal effector genes and other sequences that are duplicated or absent from indispensable core chromosome ends. Well-studied effectors found on different core chromosomes in rice pathogens appear side-by-side in wheat pathogen mini-chromosomes. The rice pathogen often overcomes deployed resistance genes by deleting triggering effector genes. We propose that the fast-evolving effector-rich compartment of the wheat blast fungus is a combination of core chromosome ends and mobile mini-chromosomes that are easily lost from individual strains. Localization of effectors on mini-chromosomes would therefore accelerate pathogen adaptation in the field.

Characterization of the Race Structure of Leptosphaeria maculans Causing Blackleg of Winter Canola in Oklahoma and Kansas

Blackleg, caused by the fungus Leptosphaeria maculans, is a widespread disease of winter canola (Brassica napus) in Oklahoma and Kansas. Deployment of genetic resistance is the primary strategy for managing blackleg. Resistance genes (Rlm) in canola interact with avirulence genes in the fungus (AvrLm) in a gene-for-gene manner. Little is known about the diversity and frequency of avirulence genes and the race structure in the region. Isolates of Leptosphaeria spp. were collected from diseased leaves in nine counties in Oklahoma and one county in Kansas from 2009 to 2013. Based on pathogenicity and PCR amplification of mating type and species-specific internal transcribed spacer loci, most isolates (n = 90) were L. maculans. The presence of avirulence genes was evaluated using phenotypic interactions on cotyledons of differential cultivars with Rlm1, Rlm2, Rlm3, and Rlm4 and amplification of AvrLm1, AvrLm4-7, and AvrLm6 by PCR. The avirulence alleles AvrLm6 and AvrLm7 were present in the entire L. maculans population. AvrLm1 was found in 34% of the population, AvrLm2 in 4%, and AvrLm4 in only 1%. A total of five races, defined as combinations of avirulence alleles, were identified that included AvrLm1-2-6-7, AvrLm2-6-7, AvrLm4-6-7, AvrLm1-6-7, and AvrLm6-7. Races virulent on the most Rlm genes, AvrLm1-6-7 at 32% and AvrLm6-7 at 62%, were predominant. Defining the avirulence allele frequency and race structure of L. maculans should be useful for the identification and development of resistant cultivars and hybrids for blackleg management in the region. The results suggest that Rlm6 and Rlm7 would be effective, although their deployment should be integrated with quantitative resistance and cultural practices, such as crop rotation, that limit selection pressure on Rlm genes.

Effector gene reshuffling involves dispensable mini-chromosomes in the wheat blast fungus

Newly emerged wheat blast disease is a serious threat to global wheat production. Wheat blast is caused by a distinct, exceptionally diverse lineage of the fungus causing rice blast disease. Through sequencing a recent field isolate, we report a reference genome that includes seven core chromosomes and mini-chromosome sequences that harbor effector genes normally found on ends of core chromosomes in other strains. No mini-chromosomes were observed in an early field strain, and at least two from another isolate each contain different effector genes and core chromosome end sequences. The mini-chromosome is enriched in transposons occurring most frequently at core chromosome ends. Additionally, transposons in mini-chromosomes lack the characteristic signature for inactivation by repeat-induced point (RIP) mutation genome defenses. Our results, collectively, indicate that dispensable mini-chromosomes and core chromosomes undergo divergent evolutionary trajectories, and mini-chromosomes and core chromosome ends are coupled as a mobile, fast-evolving effector compartment in the wheat pathogen genome.

Transcriptional Insight Into Brassica napus Resistance Genes LepR3 and Rlm2-Mediated Defense Response Against the Leptosphaeria maculans Infection

The phytopathogenic fungus Leptosphaeria maculans causes the blackleg disease on Brassica napus, resulting in severe loss of rapeseed production. Breeding of resistant cultivars containing race-specific resistance genes is provably effective to combat this disease. While two allelic resistance genes LepR3 and Rlm2 recognizing L. maculans avirulence genes AvrLm1 and AvrLm2 at plant apoplastic space have been cloned in B. napus, the downstream gene expression network underlying the resistance remains elusive. In this study, transgenic lines expressing LepR3 and Rlm2 were created in the susceptible "Westar" cultivar and inoculated with L. maculans isolates containing different sets of AvrLm1 and AvrLm2 for comparative transcriptomic analysis. Through grouping the RNA-seq data based on different levels of defense response, we find LepR3 and Rlm2 orchestrate a hierarchically regulated gene expression network, consisting of induced ABA acting independently of the disease reaction, activation of signal transduction pathways with gradually increasing intensity from compatible to incompatible interaction, and specifically induced enzymatic and chemical actions contributing to hypersensitive response with recognition of AvrLm1 and AvrLm2. This study provides an unconventional investigation into LepR3 and Rlm2-mediated plant defense machinery and adds novel insight into the interaction between surface-localized receptor-like proteins (RLPs) and apoplastic fungal pathogens.

A two genes - for - one gene interaction between Leptosphaeria maculans and Brassica napus

Interactions between Leptosphaeria maculans, causal agent of stem canker of oilseed rape, and its Brassica hosts are models of choice to explore the multiplicity of 'gene-for-gene' complementarities and how they diversified to increased complexity in the course of plant-pathogen co-evolution. Here, we support this postulate by investigating the AvrLm10 avirulence that induces a resistance response when recognized by the Brassica nigra resistance gene Rlm10. Using genome-assisted map-based cloning, we identified and cloned two AvrLm10 candidates as two genes in opposite transcriptional orientation located in a subtelomeric repeat-rich region of the genome. The AvrLm10 genes encode small secreted proteins and show expression profiles in planta similar to those of all L. maculans avirulence genes identified so far. Complementation and silencing assays indicated that both genes are necessary to trigger Rlm10 resistance. Three assays for protein-protein interactions showed that the two AvrLm10 proteins interact physically in vitro and in planta. Some avirulence genes are recognized by two distinct resistance genes and some avirulence genes hide the recognition specificities of another. Our L. maculans model illustrates an additional case where two genes located in opposite transcriptional orientation are necessary to induce resistance. Interestingly, orthologues exist for both L. maculans genes in other phytopathogenic species, with a similar genome organization, which may point to an important conserved effector function linked to heterodimerization of the two proteins.

Dissecting R gene and host genetic background effect on the Brassica napus defense response to Leptosphaeria maculans

While our understanding of the genetics underlying the Brassica-Leptosphaeria pathosystem has advanced greatly in the last decade, differences in molecular responses due to interaction between resistance genes and host genetic background has not been studied. We applied RNAseq technology to monitor the transcriptome profiles of Brassica napus (Bn) lines carrying one of four blackleg R genes (Rlm2, Rlm3, LepR1 & LepR2) in Topas or Westar background, during the early stages of infection by a Leptosphaeria maculans (Lm) isolate carrying the corresponding Avr genes. We observed upregulation of host genes involved in hormone signalling, cell wall thickening, response to chitin and glucosinolate production in all R gene lines at 3 day after inoculation (dai) albeit having higher level of expression in LepR1 and Rlm2 than in Rlm3 and LepR2 lines. Bn-SOBIR1 (Suppressor Of BIR1-1), a receptor like kinase (RLK) that forms complex receptor like proteins (RLPs) was highly expressed in LepR1 and Rlm2 at 3 dai. In contrast Bn-SOBIR1 induction was low in Rlm3 line, which could indicate that Rlm3 may function independent of SOBIR1. Expression of Salicylic acid (SA) related defense was enhanced in LepR1 and Rlm2 at 3 dai. In contrast to SA, expression of Bn genes with homology to PDF1.2, a jasmonic acid (JA) pathway marker, were increased in all Rlm and LepR lines at 6 and 9 dai. Effect of host genetic background on induction of defense, was determined by comparison of LepR1 and LepR2 in Topas vs Westar genotype (i.e. T-LepR1 vs W-LepR1 and T-LepR2 vs W-LepR2). In both cases (regardless of R gene) overall number of defense related genes at the earliest time point (3 dai) was higher in Tops compared to Westar. SA and JA markers genes such as PR1 and PDF1.2 were more induced in Topas compared to Westar introgression lines at this time point. Even in the absence of any R gene, effect of Topas genotype in enhanced defense, was also evident by the induction of PDF1.2 that started at a low level at 3 dai and peaked at 6 and 9 dai, while no induction in Westar genotype was observed at any of these time points. Overall, variation in time and intensity of expression of genes related to defense, was clearly dependent on both R gene and the host genotype.

The Modern View of B Chromosomes Under the Impact of High Scale Omics Analyses

Supernumerary B chromosomes (Bs) are extra karyotype units in addition to A chromosomes, and are found in some fungi and thousands of animals and plant species. Bs are uniquely characterized due to their non-Mendelian inheritance, and represent one of the best examples of genomic conflict. Over the last decades, their genetic composition, function and evolution have remained an unresolved query, although a few successful attempts have been made to address these phenomena. A classical concept based on cytogenetics and genetics is that Bs are selfish and abundant with DNA repeats and transposons, and in most cases, they do not carry any function. However, recently, the modern quantum development of high scale multi-omics techniques has shifted B research towards a new-born field that we call "B-omics". We review the recent literature and add novel perspectives to the B research, discussing the role of new technologies to understand the mechanistic perspectives of the molecular evolution and function of Bs. The modern view states that B chromosomes are enriched with genes for many significant biological functions, including but not limited to the interesting set of genes related to cell cycle and chromosome structure. Furthermore, the presence of B chromosomes could favor genomic rearrangements and influence the nuclear environment affecting the function of other chromatin regions. We hypothesize that B chromosomes might play a key function in driving their transmission and maintenance inside the cell, as well as offer an extra genomic compartment for evolution.

Characterization of Callose Deposition and Analysis of the Callose Synthase Gene Family of Brassica napus in Response to Leptosphaeria maculans

Callose plays a critical role in different biological processes including development as well as in the response to multiple biotic and abiotic stresses. In this study, we characterized the callose deposition in cotyledons of different Brassica napus varieties post-inoculated with different Leptosphaeria maculans isolates. Further, members of the callose synthase gene were identified from the whole genome of B. napus using the 12 Arabidopsis thaniana callose synthase protein sequences, and were then classified into three groups based on their phylogenetic relationships. Chromosomal location and duplication patterns indicated uneven distribution and segmental duplication patterns of BnCalS genes in the B. napus genome. Subsequently, gene structures, conserved domains analysis, and protein properties were analyzed for BnCalS genes. In addition, 12 B. napus orthologs of the AtCalS were selected for investigating the tissue expression pattern, indicating diverse expression patterns for these BnCalS genes. Responses of the selected 12 orthologs and all the BnCalS genes were characterized in the different types (AvrLm1-Rlm1, AvrLm4-Rlm4, AvrLepR1-LepR1) of B. napus–L. maculans interactions and B. napus-Leptosphaeria biglobosa interactions, implying their potential roles in response to Leptosphaeria infection.