Genome structure impacts molecular evolution at the AvrLm1 avirulence locus of the plant pathogen Leptosphaeria maculans

Leptosphaeria maculans, a dothideomycete fungus causing stem canker on oilseed rape, develops gene-for-gene interactions with its host plants. It has the ability to rapidly adapt to selection pressure exerted by cultivars harbouring novel resistance genes as exemplified recently by the 3-year evolution towards virulence at the AvrLm1 locus in French populations. The AvrLm1 avirulence gene was recently cloned and shown to be a solo gene within a 269 kb non-coding, heterochromatin-like region. Here we describe the sequencing of the AvrLm1 genomic region in one avirulent and two virulent isolates to investigate the molecular basis of evolution towards virulence at the AvrLm1 locus. For these virulent isolates, the gain of virulence was linked to a 260 kb deletion of a chromosomal segment spanning AvrLm1 and deletion breakpoints were identical or similar. Among the 460 isolates analysed from France, Australia and Mexico, a similar large deletion was apparent in > 90% of the virulent isolates. Deletion breakpoints were also strongly conserved in most of the virulent isolates, which led to the hypothesis that a unique deletion event leading to the avrLm1 virulence has diffused in pathogen populations. These data finally suggest that retrotransposons are key drivers in genome evolution and adaptation to novel selection pressure in L. maculans.     

Lost in the middle of nowhere: the AvrLm1 avirulence gene of the Dothideomycete Leptosphaeria maculans

Leptosphaeria maculans, a Dothideomycete causing stem canker on oilseed rape (Brassica napus), develops gene-for-gene interactions with its host plants. To date, nine resistance genes (Rlm1-9) have been identified in Brassica spp. The corresponding nine avirulence genes (AvrLm1-9) in L. maculans have been mapped at four independent loci, thereby revealing two clusters of three and four linked avirulence genes. Here, we report the completion of map-based cloning of AvrLm1. AvrLm1 was genetically delineated within a 7.3 centimorgan interval corresponding to a 439 kb BAC contig. AvrLm1 is a single copy gene isolated within a 269 kb non-coding, heterochromatin-like region. The region comprised a number of degenerated, nested copies of four long-terminal repeat (LTR) retrotransposons, including Pholy and three novel Gypsy-like retrotransposons. AvrLm1 restored the avirulent phenotype on Rlm1 cultivars following functional complementation of virulent isolates. AvrLm1 homologues were not detected in other Leptosphaeria species or in known fungal genomes including the closely related species Stagonospora nodorum. The predicted AvrLm1 protein is composed of 205 amino acids, of which only one is a cysteine residue. It contains a peptide signal suggesting extracellular localization. Unlike most other fungal avirulence genes, AvrLm1 is constitutively expressed, with a probable increased level of expression upon plant infection, suggesting the absence of tight regulation of AvrLm1 expression.     

Genetic mapping of the <Emphasis Type="Italic">Leptosphaeria maculans</Emphasis> avirulence gene corresponding to the <Emphasis Type="Italic">LepR1</Emphasis> resistance gene of <Emphasis Type="Italic">Brassica napus</Emphasis>

AvrLepR1 of the fungal pathogen Leptosphaeria maculans is the avirulence gene that corresponds to Brassica LepR1, a plant gene controlling dominant, race-specific resistance to this pathogen. An in vitro cross between the virulent L. maculans isolate, 87-41, and the avirulent isolate, 99-56, was performed in order to map the AvrLepR1 gene. The disease reactions of the 94 of the resulting F₁ progenies were tested on the canola line ddm-12-6s-1, which carries LepR1. There were 44 avirulent progenies and 50 virulent progenies suggesting a 1:1 segregation ratio and that the avirulence of 99-56 on ddm-12-6s-1 is controlled by a single gene. Tetrad analysis also indicated a 1:1 segregation ratio. The AvrLepR1 gene was positioned on a genetic map of L. maculans relative to 259 sequence-related amplified polymorphism (SRAP) markers, two cloned avirulence genes (AvrLm1 and AvrLm4-7) and the mating type locus (MAT1). The genetic map consisted of 36 linkage groups, ranging in size from 13.1 to 163.7 cM, and spanned a total of 2,076.4 cM. The AvrLepR1 locus was mapped to linkage group 4, in the 13.1 cM interval flanked by the SRAP markers SBG49-110 and FT161-223. The AvrLm4-7 locus was also positioned on linkage group 4, close to but distinct from the AvrLepR1 locus, in the 5.4 cM interval flanked by FT161-223 and P1314-300. This work will make possible the further characterization and map-based cloning of AvrLepR1. A combination of genetic mapping and pathogenicity tests demonstrated that AvrLepR1 is different from each of the L. maculans avirulence genes that have been characterized previously.     

Genetic mapping of Pyrenophora teres f. teres genes conferring avirulence on barley

A Pyrenophora teres f. teres cross between isolates 0-1 and 15A was used to evaluate the genetics of avirulence associated with barley lines Canadian Lake Shore (CLS), Tifang, and Prato. 15A is avirulent on Tifang and CLS, but virulent on Prato. Conversely, 0-1 is avirulent on Prato, but virulent on Tifang and CLS. Avirulence:virulence on Tifang and CLS segregated 1:1, whereas avirulence:virulence on Prato segregated 3:1. An AFLP-based linkage map was constructed and used to identify a single locus derived from 15A (AvrHar) conferring avirulence to Tifang and CLS. Virulence on Prato was conferred by two epistatic genes (AvrPra1 and AvrPra2). AvrPra2 co-segregated with AvrHar, but the two genes from opposite parents conferred opposite reactions. This work provides the foundation for the isolation of these avirulence genes.     

Avirulence in the wheat septoria tritici leaf blotch fungus Mycosphaerella graminicola is controlled by a single locus

Segregation of avirulence in Mycosphaerella graminicola, a heterothallic ascomycete that causes wheat septoria tritici leaf blotch, was studied in F1, BC1, and F2 populations by inoculation assays on five wheat cultivars in the seedling stage and by amplified fragment length polymorphism and random amplified polymorphic DNA analyses. F1 was generated by crossing isolates IPO323 (avirulent) and IPO94269 (virulent). All F1, BC1, and F2 progeny isolates were virulent on the susceptible check cultivar Taichung 29 and were avirulent on the resistant check cultivar Kavkav-K4500. Avirulence segregation was observed in F1 and in several BC1 and F2 generations on the differential cultivars Shafir, Kavkaz, and Veranopolis at a 1:1 ratio. Avirulence for the three differential cultivars always cosegregated. We conclude that avirulence in isolate IPO323 is controlled by a single, seemingly complex locus.     

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.     

Dominance and recessiveness at loci for virulence against potato and tomato in Phytophthora infestans

Summary In this study we investigated the genetic control of virulence in the diploid fungal pathogen, Phytophthora infestans, against host resistance genes R1, R2, R3, and R4 (potato) and Ph1 (tomato). For four of these virulence traits, the presence or absence of segregation indicated conclusively which phenotype was dominant. We observed a 31 (virulentavirulent) segregation on R2 in the progeny of parents which were both virulent, suggesting that virulence is dominant and both parents are heterozygous. In a cross in which one parent was virulent and the other avirulent on potato gene R3, all progeny tested were avirulent, so avirulence against R3 is dominant. The same virulent parent crossed with a different avirulent parent produced virulent and avirulent progeny in a 13 ratio, indicating that a second locus may be involved. The progeny of two parents virulent on R4 segregated for virulence and avirulence, so virulence against R4 is dominant. For Ph1, a 13 segregation in the progeny of two avirulent parents showed that the avirulent phenotype is dominant, and a 31 ration in a second cross suggested the involvement of a second locus. The segregations for virulence against R1 did not indicate which phenotype was dominant, but did suggest singlelocus control.     

A genomic map enriched for markers linked to Avr1 in Cronartium quercuum f.sp. fusiforme

A novel approach is presented to map avirulence gene Avr1 in the basidiomycete Cronartium quercuum f.sp. fusiforme, the causal agent of fusiform rust disease in pines. DNA markers tightly linked to resistance gene Fr1 in loblolly pine tree 10-5 were used to classify 10-5 seedling progeny as either resistant or susceptible. A single dikaryotic isolate (P2) heterozygous at the corresponding Avr1 gene was developed by crossing Fr1 avirulent isolate SC20-21 with Fr1 virulent isolate NC2-40. Bulk basidiospore inoculum derived from isolate P2 was used to challenge the pine progeny. The ability to unambiguously marker classify 10-5 progeny as resistant (selecting for virulence) or susceptible (non-selecting) permitted the genetic mapping of the corresponding Avr1 gene by bulked segregant analysis. Using this approach, 14 genetic markers significantly linked to Avr1 were identified and placed within the context of a genome-wide linkage map produced for isolate P2 using samples from susceptible seedlings.     

Use of high resolution melting analysis to genotype the avirulence AvrLm4‐7 gene of Leptosphaeria maculans,a fungal pathogen of Brassica napus

The oilseed rape (Brassica napus) stem canker disease, due to the fungal pathogen Leptosphaeria maculans, is mainly controlled by host genetic resistance. Since 2004, the specific resistance gene Rlm7 is widely used in France. Specific resistance is effective when fungal populations are mainly composed of avirulent isolates. The development of molecular tools for the identification of virulent isolates towards Rlm7 was needed to undertake large‐scale surveys and to monitor the emergence of virulent populations in fields. Previous studies have described a large diversity of molecular events leading to virulence towards Rlm7, rendering conventional polymerase chain reaction (PCR) methods inapplicable to identify virulent isolates. Interestingly, a very limited nucleotide polymorphism was observed for avirulent, AvrLm7, alleles. Such characteristics were exploited here to develop a diagnostic method based on high resolution melting (HRM) analysis of the AvrLm4‐7 gene. High resolution melting analysis of a collection of 206 reference isolates revealed only four different profiles within 100 avirulent isolates and 87% of virulent isolates showed either no amplification or HRM curves distinct from those of avirulent isolates. The reliability of the method was confirmed using a second set of 119 unknown isolates, for which biological phenotyping and HRM genotyping were in agreement for 93% of the isolates. HRM combined with the PCR amplification of a larger fragment encompassing AvrLm4‐7 led to a correct diagnostic for 97.5% of the isolates.     

Chromosomal deletion in isolates of Phytophthora infestans correlates with virulence on R3, R10, and R11 potato lines.

In Phytophthora infestans, a cluster of three dominant avirulence genes is located on the distal part of linkage group VIII. In a mapping population from a cross between two Dutch field isolates, probe M5.1, derived from an amplified fragment length polymorphism (AFLP) marker linked to the Avr3-Avr10-Avr11 cluster, hybridized only to DNA from the parent and F1 progeny that is avirulent on potato lines carrying the R3, R10, and R11 resistance gene. In the virulent parent and the virulent progeny, no M5.1 homologue was detected, demonstrating a deletion on that part of linkage group VIII. P. infestans is diploid, so the avirulent strains must be hemizygous for the region concerned. A similar situation was found in another mapping population from two Mexican strains. The deletion was also found to occur in many field isolates. In a large set of unique isolates collected in The Netherlands from 1980 to 1991, 37% had no M5.1 homologue and the deletion correlated strongly with gain of virulence on potato lines carrying R3, R10, and R11. Also, in some old isolates that belong to a single clonal lineage (US-1) and are thus highly homogenous, deletions at the M5.1 locus were detected, indicating that this region is unstable.     

Magnaporthe Grisea Genes for Pathogenicity and Virulence Identified through a Series of Backcrosses

We have identified genes for pathogenicity toward rice (Oryza sativa) and genes for virulence toward specific rice cultivars in the plant pathogenic fungus Magnaporthe grisea. A genetic cross was conducted between the weeping lovegrass (Eragrostis curvula) pathogen 4091-5-8, a highly fertile, hermaphroditic laboratory strain, and the rice pathogen O-135, a poorly fertile, female-sterile field isolate that infects weeping lovegrass as well as rice. A six-generation backcrossing scheme was then undertaken with the rice pathogen as the recurrent parent. One goal of these crosses was to generate rice pathogenic progeny with the high fertility characteristic of strain 4091-5-8, which would permit rigorous genetic analysis of rice pathogens. Therefore, progeny strains to be used as parents for backcross generations were chosen only on the basis of fertility. The ratios of pathogenic to nonpathogenic (and virulent to avirulent) progeny through the backcross generations suggested that the starting parent strains differ in two types of genes that control the ability to infect rice. First, they differ by polygenic factors that determine the extent of lesion development achieved by those progeny that infect rice. These genes do not appear to play a role in infection of weeping lovegrass because both parents and all progeny infect weeping lovegrass. Second, the parents differ by simple Mendelian determinants, ``avirulence genes,' that govern virulence toward specific rice cultivars in all-or-none fashion. Several crosses confirm the segregation of three unlinked avirulence genes, Avr1-CO39, Avr1-M201 and Avr1-YAMO, alleles of which determine avirulence on rice cultivars CO39, M201, and Yashiro-mochi, respectively. Interestingly, avirulence alleles of Avr1-CO39, Avr1-M201 and Avr1-YAMO were inherited from the parent strain 4091-5-8, which is a nonpathogen of rice. Middle repetitive DNA sequences (``MGR sequences'), present in approximately 40-50 copies in the genome of the rice pathogen parent, and in very low copy number in the genome of the nonpathogen of rice, were used as physical markers to monitor restoration of the rice pathogen genetic background during introgression of fertility. The introgression of highest levels of fertility into the most successful rice pathogen progeny was incomplete by the sixth generation, perhaps a consequence of genetic linkage between genes for fertility and genes for rice pathogenicity. One chromosomal DNA segment with MGR sequence homology appeared to be linked to the gene Avr1-CO39. Finally, many of the crosses described in this paper exhibited a characteristic common to many crosses involving M. grisea rice pathogen field isolates. That is, pigment-defective mutants frequently appeared among the progeny.     

One gene‐one name: the AvrLmJ1 avirulence gene of Leptosphaeria maculans is AvrLm5

Leptosphaeria maculans, the causal agent of blackleg disease, interacts with Brassica napus (oilseed rape, canola) and other Brassica hosts in a gene‐for‐gene manner. The avirulence gene AvrLmJ1 has been cloned previously and shown to interact with an unidentified Brassica juncea resistance gene. In this study, we show that the AvrLmJ1 gene maps to the same position as the AvrLm5 locus. Furthermore, isolates complemented with the AvrLmJ1 locus confer avirulence towards B. juncea genotypes harbouring Rlm5. These findings demonstrate that AvrLmJ1 is AvrLm5 and highlight the need for shared resources to characterize accurately avirulence and/or resistance genes.     

Identification and characterization of DNA markers associated with a locus conferring virulence on barley in the plant pathogenic fungus Cochliobolus sativus

Cochliobolus sativus is a plant pathogenic fungus that causes spot blotch on barley and wheat. Virulence of a pathotype-2 isolate (ND90Pr) on barley cultivar Bowman was previously determined to be controlled by a single locus. To identify DNA markers associated with this virulence locus, amplified fragment length polymorphism (AFLP) analysis was conducted on 104 progeny isolates derived from a cross between isolates ND90Pr (exhibiting high virulence on Bowman) and ND93-1 (exhibiting low virulence on Bowman). Among 115 AFLP markers identified, 14 were linked to the virulence locus VHv1 in isolate ND90Pr, six of which co-segregated with VHv1. Two (E-AG/M-CA-207 and E-AG/M-CG-121) of the six co-segregating AFLP markers were cloned and used to probe genomic DNAs from the fungal parents and progeny. Both markers hybridized only with DNAs from ND90Pr and the virulent progeny. These two cloned markers were also used as probes to survey field isolates of C. sativus collected from different regions of the world and again only hybridized to DNAs from isolates that had the same virulence phenotype as ND90Pr. The results of this study indicate that E-AG/M-CA-207 and E-AG/M-CG-121 are closely linked to VHv1 and are unique to isolates carrying the virulence locus. Development of a linkage group, coupled with the identification of closely linked molecular markers, will facilitate the cloning of the virulence gene VHv1 in C. sativus by map-based cloning.     

A genetic map of the lettuce downy mildew pathogen,Bremia lactucae,constructed from molecular markers and avirulence genes

The genetic map of Bremia lactucae was expanded utilizing 97 F(1) progeny derived from a cross between Finnish and Californian isolates (SF5xC82P24). Genetic maps were constructed for each parent utilizing 7 avirulence genes, 83 RFLP markers, and 347 AFLP markers, and a consensus map was constructed from the complete data set. The framework map for SF5 contained 24 linkage groups distributed over 835cM; the map for C82P24 contained 21 linkage groups distributed over 606cM. The consensus map contained 12 linkage groups with markers from both parents and 24 parent-specific groups. Six avirulence genes mapped to different linkage groups; four were located at the ends of linkage groups. The closest linkages between molecular markers and avirulence genes were 3cM to Avr4 and 1cM to Avr7. Mating type seemed to be determined by a single locus, where the heterozygote determined the B(2) type and the homozygous recessive genotype determined the B(1) type.     

Marker-based cloning of the region containing the UhAvr1 avirulence gene from the basidiomycete barley pathogen Ustilago hordei

Race-cultivar specialization during the interaction of the basidiomycete smut pathogen Ustilago hordei with its barley host was described in the 1940s. Subsequent genetic analyses revealed the presence of dominant avirulence genes in the pathogen that conform to the gene-for-gene theory. This pathosystem therefore presents an opportunity for the molecular genetic characterization of fungal genes controlling avirulence. We performed a cross between U. hordei strains to obtain 54 progeny segregating for three dominant avirulence genes on three differential barley cultivars. Bulked segregant analysis was used to identify RAPD and AFLP markers tightly linked to the avirulence gene UhAvr1. The UhAvr1 gene is located in an area containing repetitive DNA and this region is undetectable in cosmid libraries prepared from the avirulent parental strain. PCR and hybridization probes developed from the linked markers were therefore used to identify cosmid clones from the virulent (Uhavr1) parent. By walking on Uhavr1-linked cosmid clones, a nonrepetitive, nearby probe was found that recognized five overlapping BAC clones spanning 170 kb from the UhAvr1 parent. A contig of the clones in the UhAvr1 region was constructed and selected probes were used for RFLP analysis of the segregating population. This approach genetically defined an approximately 80-kb region that carries the UhAvr1 gene and provided cloned sequences for subsequent genetic analysis. UhAvr1 represents the first avirulence gene cloned from a basidiomycete plant pathogen.     

A gene from Xanthomonas campestris pv. vesicatoria that determines avirulence in tomato is related to avrBs3.

Strains of Xanthomonas campestris pv. vesicatoria that were avirulent in tomato leaves but virulent in pepper leaves were identified. A cloned gene, avrBsP, from one of the strains, Xv 87-7, converted a virulent strain in tomato to avirulent in tomato. A 1.7-kb subclone containing the avirulence gene cross-hybridized with the avirulence gene, which determines race 1 within the pepper group of strains (avrBs3). However, the two avirulence genes differ in their biological activity. The base sequences of the two avirulence genes were almost identical through the 1.7-kb segment of avrBsP, with significant differences only in some bases in the repeat region.     

Obviation of wheat resistance to the Hessian fly through systemic induced susceptibility

Unlike most documented plant-insect interactions, Hessian fly-resistance [Mayetiola destructor (Say)] in wheat (Triticum aestivum L.) is initiated by a gene-for-gene recognition event in which plants carrying a specific R gene recognize salivary effectors encoded by a corresponding larval avirulence gene. However, dual infestation resulting from oviposition by virulent insects from 5 d before to 3 d after oviposition by avirulent insects on the same host plant, lead to systemic induced susceptibility, obviation of resistance, and ultimately the survival of both virulent and genetically avirulent progeny to adulthood. Simultaneous oviposition allowed greater survival of avirulent progeny than ovipositions separated by larger intervals. Because of the induction of plant resistance, hatch of avirulent larvae before virulent was more detrimental to rate of development than hatch of virulent before avirulent larvae. Obviation of resistance was not localized to the leaf being attacked by the virulent larvae, but also functioned across spatial distance into younger leaves. This research suggests that virulent Hessian fly larvae directly suppress the defense response of wheat, thus providing a refuge for avirulent genotypes, preserving diversity in field populations and increasing durability of deployed resistance genes.     

Genome Structure and Reproductive Behaviour Influence the Evolutionary Potential of a Fungal Phytopathogen

Modern agriculture favours the selection and spread of novel plant diseases. Furthermore, crop genetic resistance against pathogens is often rendered ineffective within a few years of its commercial deployment. Leptosphaeria maculans, the cause of phoma stem canker of oilseed rape, develops gene-for-gene interactions with its host plant, and has a high evolutionary potential to render ineffective novel sources of resistance in crops. Here, we established a four-year field experiment to monitor the evolution of populations confronted with the newly released Rlm7 resistance and to investigate the nature of the mutations responsible for virulence against Rlm7. A total of 2551 fungal isolates were collected from experimental crops of a Rlm7 cultivar or a cultivar without Rlm7. All isolates were phenotyped for virulence and a subset was genotyped with neutral genetic markers. Virulent isolates were investigated for molecular events at the AvrLm4-7 locus. Whilst virulent isolates were not found in neighbouring crops, their frequency had reached 36% in the experimental field after four years. An extreme diversity of independent molecular events leading to virulence was identified in populations, with large-scale Repeat Induced Point mutations or complete deletion of AvrLm4-7 being the most frequent. Our data suggest that increased mutability of fungal genes involved in the interactions with plants is directly related to their genomic environment and reproductive system. Thus, rapid allelic diversification of avirulence genes can be generated in L. maculans populations in a single field provided that large population sizes and sexual reproduction are favoured by agricultural practices.     

Evolution of linked avirulence effectors in Leptosphaeria maculans is affected by genomic environment and exposure to resistance genes in host plants

Brassica napus (canola) cultivars and isolates of the blackleg fungus, Leptosphaeria maculans interact in a 'gene for gene' manner whereby plant resistance (R) genes are complementary to pathogen avirulence (Avr) genes. Avirulence genes encode proteins that belong to a class of pathogen molecules known as effectors, which includes small secreted proteins that play a role in disease. In Australia in 2003 canola cultivars with the Rlm1 resistance gene suffered a breakdown of disease resistance, resulting in severe yield losses. This was associated with a large increase in the frequency of virulence alleles of the complementary avirulence gene, AvrLm1, in fungal populations. Surprisingly, the frequency of virulence alleles of AvrLm6 (complementary to Rlm6) also increased dramatically, even though the cultivars did not contain Rlm6. In the L. maculans genome, AvrLm1 and AvrLm6 are linked along with five other genes in a region interspersed with transposable elements that have been degenerated by Repeat-Induced Point (RIP) mutations. Analyses of 295 Australian isolates showed deletions, RIP mutations and/or non-RIP derived amino acid substitutions in the predicted proteins encoded by these seven genes. The degree of RIP mutations within single copy sequences in this region was proportional to their proximity to the degenerated transposable elements. The RIP alleles were monophyletic and were present only in isolates collected after resistance conferred by Rlm1 broke down, whereas deletion alleles belonged to several polyphyletic lineages and were present before and after the resistance breakdown. Thus, genomic environment and exposure to resistance genes in B. napus has affected the evolution of these linked avirulence genes in L. maculans.