Chickpea production is constrained worldwide by the necrotrophic fungal pathogen Ascochyta rabiei, the causal agent of Ascochyta blight (AB). To reduce the impact of this disease, novel sources of resistance are required in chickpea cultivars. Here, we screened a new collection of wild Cicer accessions for AB resistance and identified accessions resistant to multiple, highly pathogenic isolates. In addition to this, analyses demonstrated that some collection sites of C. echinospermum harbor predominantly resistant accessions, knowledge that can inform future collection missions. Furthermore, a genome-wide association study identified regions of the C. reticulatum genome associated with AB resistance and investigation of these regions identified candidate resistance genes. Taken together, these results can be utilized to enhance the resistance of chickpea cultivars to this globally yield-limiting disease.

Ascochyta blight disease, caused by the necrotrophic fungusAscochyta 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 ofA. rabieiisolates 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 AustralianA. rabieipopulation 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 AustralianA. rabieipopulation. 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 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 an unknown function. Comparison of A. lentis genome features with the recently published reference assembly for closely related A. rabiei show 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.

Chocolate spot is a major fungal disease of faba bean caused by the ascomycete fungus, Botrytis fabae. B. fabae is also implicated in botrytis gray mold disease in lentils, along with B. cinerea. Here we have isolated and characterized two B. fabae isolates from chocolate spot lesions on faba bean leaves. In plant disease assays on faba bean and lentil, B. fabae was more aggressive than B. cinerea and we observed variation in susceptibility among a small set of cultivars for both plant hosts. Using light microscopy, we observed a spreading, generalized necrosis response in faba bean toward B. fabae. In contrast, the plant response to B. cinerea was localized to epidermal cells underlying germinated spores and appressoria. In addition to the species characterization of B. fabae, we produced genome assemblies for both B. fabae isolates using Illumina sequencing. Genome sequencing coverage and assembly size for B. fabae isolates, were 27x and 45x, and 43.2 and 44.5 Mb, respectively. Following genome assembly and annotation, carbohydrate-active enzyme (CAZymes) and effector genes were predicted. There were no major differences in the numbers of each of the major classes of CAZymes. We predicted 29 effector genes for B. fabae, and using the same selection criteria for B. cinerea, we predicted 34 putative effector genes. For five of the predicted effector genes, the pairwise dN/dS ratio between orthologs from B. fabae and B. cinerea was greater than 1.0, suggesting positive selection and the potential evolution of molecular mechanisms for host specificity in B. fabae. Furthermore, a homology search of secondary metabolite clusters revealed the absence of the B. cinerea phytotoxin botrydial and several other uncharacterized secondary metabolite biosynthesis genes from B. fabae. Although there were no obvious differences in the number or proportional representation of different transposable element classes, the overall proportion of AT-rich DNA sequence in B. fabae was double that of B. cinerea.

Ascochyta blight is a damaging disease that affects the stems, leaves, and pods of field pea ( Pisum sativum) and impacts yield and grain quality. In Australia, field pea Ascochyta blight is primarily caused by the necrotrophic fungal species Peyronellaea pinodes and Ascochyta koolunga. In this study, we screened 1,276 Pisum spp. germplasm accessions in seedling disease assays with a mix of three isolates of P. pinodes and 641 accessions with three mixed isolates of A. koolunga (513 accessions were screened with both species). A selection of three P. sativum accessions with low disease scores for either pathogen, or in some cases both, were crossed with Australian field pea varieties PBA Gunyah and PBA Oura, and recombinant inbred line populations were made. Populations at the F_3:4 and F_4:5 generation were phenotyped for their disease response to P. pinodes and A. koolunga, and genotypes were determined using the diversity arrays technology genotyping method. Marker-trait associations were identified using a genome-wide association study approach. Trait-associated loci were mapped to the published P. sativum genome assembly, and candidate resistance gene analogues were identified in the corresponding genomic regions. One locus on chromosome 2 (LG1) was associated with resistance to P. pinodes, and the 8 Mb genomic region contains 156 genes, two of which are serine/threonine protein kinases, putatively contributing to the resistance trait. A second locus on chromosome 5 (LG3) was associated with resistance to A. koolunga, and the 35 Mb region contains 488 genes, of which five are potential candidate resistance genes, including protein kinases, a mitogen-activated protein kinase, and an ethylene-responsive protein kinase homolog.