Ascochyta rabiei is the causal organism of ascochyta blight of chickpea and is present in chickpea crops worldwide. Here we report the release of a high-quality PacBio genome assembly for the Australian A. rabiei isolate ArME14. We compare the ArME14 genome assembly with an Illumina assembly for Indian A. rabiei isolate, ArD2. The ArME14 assembly has gapless sequences for nine chromosomes with telomere sequences at both ends and 13 large contig sequences that extend to one telomere. The total length of the ArME14 assembly was 40,927,385 bp, which was 6.26 Mb longer than the ArD2 assembly. Division of the genome by OcculterCut into GC-balanced and AT-dominant segments reveals 21% of the genome contains gene-sparse, AT-rich isochores. Transposable elements and repetitive DNA sequences in the ArME14 assembly made up 15% of the genome. A total of 11,257 protein-coding genes were predicted compared with 10,596 for ArD2. Many of the predicted genes missing from the ArD2 assembly were in genomic regions adjacent to AT-rich sequence. We compared the complement of predicted transcription factors and secreted proteins for the two A. rabiei genome assemblies and found that the isolates contain almost the same set of proteins. The small number of differences could represent real differences in the gene complement between isolates or possibly result from the different sequencing methods used. Prediction pipelines were applied for carbohydrate-active enzymes, secondary metabolite clusters and putative protein effectors. We predict that ArME14 contains between 450 and 650 CAZymes, 39 putative protein effectors and 26 secondary metabolite clusters.

Ascochyta lentis causes ascochyta blight in lentil (Lens culinaris Medik.) and yield loss can be as high as 50%. With careful agronomic management practices, fungicide use and advances in breeding resistant lentil varieties, disease severity and impact to farmers have been largely controlled. However, evidence from major lentil producing countries, Canada and Australia, suggests that A. lentis isolates can change their virulence profile and level of aggressiveness over time and under different selection pressures. In this paper, we describe the first genome assembly for A. lentis for the Australian isolate Al4, through the integration of data from Illumina and PacBio SMRT sequencing. The Al4 reference genome assembly is almost 42 Mb in size and encodes 11,638 predicted genes. The Al4 genome comprises 21 full-length and gapless chromosomal contigs and two partial chromosome contigs each with one telomere. We predicted 31 secondary metabolite clusters, and 38 putative protein effectors, many of which were classified as having unknown function. Comparison of A. lentis genome features with the recently published reference assembly for closely related Ascochyta rabiei, shows that genome synteny between these species is highly conserved. However, there are several translocations and inversions of genome sequence. The location of secondary metabolite clusters near transposable element and repeat-rich genomic regions was common for A. lentis as has been reported for other fungal plant pathogens.

The plant immune system is made up of a complex response network that involves several lines of defense to fight invading pathogens. Fungal plant pathogens on the other hand, have evolved a range of ways to infect their host. The interaction between Ascochyta lentis and two lentil genotypes was explored to investigate the progression of ascochyta blight (AB) in lentils. In this study, we developed an Agrobacterium tumefaciens-mediated transformation system for A. lentis by constructing a new binary vector, pATMT-GpdGFP, for the constitutive expression of green fluorescent protein (EGFP). Green fluorescence was used as a highly efficient vital marker to study the developmental changes in A. lentis during AB disease progression on the susceptible and resistant lentil accessions, ILL6002 and ILL7537, respectively. The initial infection stages were similar in both the resistant and susceptible accessions where A. lentis uses infection structures such as germ tubes and appressoria to gain entry into the host while the host uses defense mechanisms to prevent pathogen entry. Penetration was observed at the junctions between neighbouring epidermal cells and occasionally, through the stomata. The pathogen attempted to penetrate and colonize ILL7537, but further fungal advancement appeared to be halted, and A. lentis did not enter the mesophyll. Successful entry and colonization of ILL6002 coincided with structural changes in A. lentis and the onset of necrotic lesions 5–7 days post inoculation. Once inside the leaf, A. lentis continued to grow, colonizing all parts of the leaf followed by plant cell collapse. Pycnidia-bearing spores appeared 14 days post inoculation, which marks the completion of the infection cycle. The use of fluorescent proteins in plant pathogenic fungi together with confocal laser scanning microscopy, provide a valuable tool to study the intracellular dynamics, colonization strategy and infection mechanisms during plant-pathogen interaction.

Ascochyta blight disease, caused by the necrotrophic fungus Ascochyta rabiei, is a major biotic constraint to chickpea production in Australia and worldwide. Detailed knowledge of the structure of the pathogen population and its potential to adapt to our farming practices is key to informing optimal management of the disease. This includes understanding the molecular diversity among isolates and the frequency and distribution of the isolates that have adapted to overcome host resistance across agroecologically distinct regions. Thanks to continuous monitoring efforts over the past 6 years, a comprehensive collection of A. rabiei isolates was collated from the major Australian chickpea production regions. To determine the molecular structure of the entire population, representative isolates from each collection year and growing region have been genetically characterized using a DArTseq genotyping-by-sequencing approach. The genotyped isolates were further phenotyped to determine their pathogenicity levels against a differential set of chickpea cultivars and genotype-phenotype associations were inferred. Overall, the Australian A. rabiei population displayed a far lower genetic diversity (average Nei’s gene diversity of 0.047) than detected in other populations worldwide. This may be explained by the presence of a single mating-type in Australia, MAT1-2, limiting its reproduction to a clonal mode. Despite the low detected molecular diversity, clonal selection appears to have given rise to a subset of adapted isolates that are highly pathogenic on commonly employed resistance sources, and that are occurring at an increasing frequency. Among these, a cluster of genetically similar isolates was identified, with a higher proportion of highly aggressive isolates than in the general population. The discovery of distinct genetic clusters associated with high and low isolate pathogenicity forms the foundation for the development of a molecular pathotyping tool for the Australian A. rabiei population. Application of such a tool, along with continuous monitoring of the genetic structure of the population will provide crucial information for the screening of breeding material and integrated disease management packages.

Ascochyta fabae Speg. is a serious foliar fungal disease of faba bean and a constraint to production worldwide. This study investigated the phenotypic and genotypic diversity of the A. fabae pathogen population in southern Australia and the pathogenic variability of the population was examined on a differential set of faba bean cultivars. The host set was inoculated with 154 A. fabae isolates collected from 2015 to 2018 and a range of disease reactions from high to low aggressiveness was observed. Eighty percent of isolates collected from 2015 to 2018 were categorized as pathogenicity group (PG) PG-2 (pathogenic on Farah) and were detected in every region in each year of collection. Four percent of isolates were non-pathogenic on Farah and designated as PG-1. A small group of isolates (16%) were pathogenic on the most resistant differential cultivars, PBA Samira or Nura, and these isolates were designated PG-3. Mating types of 311 isolates collected between 1991 and 2018 were determined and showed an equal ratio of MAT1–1 and MAT1–2 in the southern Australian population. The genetic diversity and population structure of 305 isolates were examined using DArTseq genotyping, and results suggest no association of genotype with any of the population descriptors viz.: collection year, region, host cultivar, mating type, or PG. A Genome-Wide Association Study (GWAS) was performed to assess genetic association with pathogenicity traits and a significant trait-associated genomic locus for disease in Farah AR and PBA Zahra, and PG was revealed. The high frequency of mating of A. fabae indicated by the wide distribution of the two mating types means changes to virulence genes would be quickly distributed to other genotypes. Continued monitoring of the A. fabae pathogen population through pathogenicity testing will be important to identify any increases in aggressiveness or emergence of novel PGs. GWAS and future genetic studies using biparental mating populations could be useful for identifying virulence genes responsible for the observed changes in pathogenicity.