Molecular Breeding for Ascochyta Blight Resistance in Lentil: Current Progress and Future Directions

Genomic selection for target traits in the Australian lentil breeding program

Genomic selection (GS) uses associations between markers and phenotypes to predict the breeding values of individuals. It can be applied early in the breeding cycle to reduce the cross-to-cross generation interval and thereby increase genetic gain per unit of time. The development of cost-effective, high-throughput genotyping platforms has revolutionized plant breeding programs by enabling the implementation of GS at the scale required to achieve impact. As a result, GS is becoming routine in plant breeding, even in minor crops such as pulses. Here we examined 2,081 breeding lines from Agriculture Victoria's national lentil breeding program for a range of target traits including grain yield, ascochyta blight resistance, botrytis grey mould resistance, salinity and boron stress tolerance, 100-grain weight, seed size index and protein content. A broad range of narrow-sense heritabilities was observed across these traits (0.24-0.66). Genomic prediction models were developed based on 64,781 genome-wide SNPs using Bayesian methodology and genomic estimated breeding values (GEBVs) were calculated. Forward cross-validation was applied to examine the prediction accuracy of GS for these targeted traits. The accuracy of GEBVs was consistently higher (0.34-0.83) than BLUP estimated breeding values (EBVs) (0.22-0.54), indicating a higher expected rate of genetic gain with GS. GS-led parental selection using early generation breeding materials also resulted in higher genetic gain compared to BLUP-based selection performed using later generation breeding lines. Our results show that implementing GS in lentil breeding will fast track the development of high-yielding cultivars with increased resistance to biotic and abiotic stresses, as well as improved seed quality traits.

Virulence Profiles and Genome-Wide Association Study for Ascochyta lentis Isolates Collected from Australian Lentil-Growing Regions

Ascochyta lentis, the causal organism of Ascochyta blight (AB) of lentil (Lens culinaris), has been shown to produce an avirulence effector protein that mediates AB resistance in certain lentil cultivars. The two known forms of the effector protein were identified from a biparental mapping population between isolates that have reciprocal virulence on 'PBA Hurricane XT' and 'Nipper'. The effector AlAvr1-1 was described for the PBA Hurricane XT-avirulent isolate P94-24 and AlAvr1-2 characterized in the PBA Hurricane XT-virulent isolate AlKewell. Here, we performed a genome-wide association study to identify other loci associated with AB for a differential set of lentil cultivars from a diverse panel of isolates collected in the Australian lentil-growing regions from 2013 to 2020. The chromosome 3 AlAvr1 locus was strongly associated with the PBA Hurricane XT, 'Indianhead', and Nipper disease responses, but one other genomic region on chromosome 11 was also associated with the Nipper disease trait. Our results corroborate earlier work that identified the AlAvr1 locus for field-collected isolates that span the period before release and after widespread adoption of PBA Hurricane XT. A multiplex PCR assay was developed to differentiate the genes AlAvr1-1 and AlAvr1-2 to predict PBA Hurricane XT avirulence and pathotype designation in the diversity panel. Increasing numbers of the PBA Hurricane XT-virulent pathotype 2 isolates across that time indicate strong selection for isolates with the AlAvr1-2 allele. Furthermore, one other region of the A. lentis genome may contribute to the pathogen-host interaction for lentil AB.

The Prospects of gene introgression from crop wild relatives into cultivated lentil for climate change mitigation

Crop wild relatives (CWRs), landraces and exotic germplasm are important sources of genetic variability, alien alleles, and useful crop traits that can help mitigate a plethora of abiotic and biotic stresses and crop yield reduction arising due to global climatic changes. In the pulse crop genus Lens, the cultivated varieties have a narrow genetic base due to recurrent selections, genetic bottleneck and linkage drag. The collection and characterization of wild Lens germplasm resources have offered new avenues for the genetic improvement and development of stress-tolerant, climate-resilient lentil varieties with sustainable yield gains to meet future food and nutritional requirements. Most of the lentil breeding traits such as high-yield, adaptation to abiotic stresses and resistance to diseases are quantitative and require the identification of quantitative trait loci (QTLs) for marker assisted selection and breeding. Advances in genetic diversity studies, genome mapping and advanced high-throughput sequencing technologies have helped identify many stress-responsive adaptive genes, quantitative trait loci (QTLs) and other useful crop traits in the CWRs. The recent integration of genomics technologies with plant breeding has resulted in the generation of dense genomic linkage maps, massive global genotyping, large transcriptomic datasets, single nucleotide polymorphisms (SNPs), expressed sequence tags (ESTs) that have advanced lentil genomic research substantially and allowed for the identification of QTLs for marker-assisted selection (MAS) and breeding. Assembly of lentil and its wild species genomes (~4Gbp) opens up newer possibilities for understanding genomic architecture and evolution of this important legume crop. This review highlights the recent strides in the characterization of wild genetic resources for useful alleles, development of high-density genetic maps, high-resolution QTL mapping, genome-wide studies, MAS, genomic selections, new databases and genome assemblies in traditionally bred genus Lens for future crop improvement amidst the impending global climate change.

De Novo Long-Read Whole-Genome Assemblies and the Comparative Pan-Genome Analysis of Ascochyta Blight Pathogens Affecting Field Pea

Ascochyta Blight (AB) is a major disease of many cool-season legumes globally. In field pea, three fungal pathogens have been identified to be responsible for this disease in Australia, namely Peyronellaea pinodes, Peyronellaea pinodella and Phoma koolunga. Limited genomic resources for these pathogens have been generated, which has hampered the implementation of effective management strategies and breeding for resistant cultivars. Using Oxford Nanopore long-read sequencing, we report the first high-quality, fully annotated, near-chromosome-level nuclear and mitochondrial genome assemblies for 18 isolates from the Australian AB complex. Comparative genome analysis was performed to elucidate the differences and similarities between species and isolates using phylogenetic relationships and functional diversity. Our data indicated that P. pinodella and P. koolunga are heterothallic, while P. pinodes is homothallic. More homology and orthologous gene clusters are shared between P. pinodes and P. pinodella compared to P. koolunga. The analysis of the repetitive DNA content showed differences in the transposable repeat composition in the genomes and their expression in the transcriptomes. Significant repeat expansion in P. koolunga's genome was seen, with strong repeat-induced point mutation (RIP) activity being evident. Phylogenetic analysis revealed that genetic diversity can be exploited for species marker development. This study provided the much-needed genetic resources and characterization of the AB species to further drive research in key areas such as disease epidemiology and host-pathogen interactions.

Genetic Augmentation of Legume Crops Using Genomic Resources and Genotyping Platforms for Nutritional Food Security

Recent advances in next generation sequencing (NGS) technologies have led the surge of genomic resources for the improvement legume crops. Advances in high throughput genotyping (HTG) and high throughput phenotyping (HTP) enable legume breeders to improve legume crops more precisely and efficiently. Now, the legume breeder can reshuffle the natural gene combinations of their choice to enhance the genetic potential of crops. These genomic resources are efficiently deployed through molecular breeding approaches for genetic augmentation of important legume crops, such as chickpea, cowpea, pigeonpea, groundnut, common bean, lentil, pea, as well as other underutilized legume crops. In the future, advances in NGS, HTG, and HTP technologies will help in the identification and assembly of superior haplotypes to tailor the legume crop varieties through haplotype-based breeding. This review article focuses on the recent development of genomic resource databases and their deployment in legume molecular breeding programmes to secure global food security.

Development and application of the Faba_bean_130K targeted next-generation sequencing SNP genotyping platform based on transcriptome sequencing

KEY MESSAGE

Large-scale faba bean transcriptome data are available, and the first genotyping platform based on liquid-phase probe targeted capture technology was developed for genetic and molecular breeding studies. Faba bean (Vicia faba L., 2n = 12) is an important food legume crop that is widely grown for multiple uses worldwide. However, no reference genome is currently available due to its very large genome size (approximately 13 Gb) and limited single nucleotide polymorphism (SNP) markers as well as highly efficient genotyping tools have been reported for faba bean. In this study, 16.7 billion clean reads were obtained from transcriptome libraries of flowers and leaves of 102 global faba bean accessions. A total of 243,120 unigenes were de novo assembled and functionally annotated. Moreover, a total of 1,579,411 SNPs were identified and further filtered according to a selection pipeline to develop a high-throughput, flexible, low-cost Faba_bean_130K targeted next-generation sequencing (TNGS) genotyping platform. A set of 69 Chinese faba bean accessions were genotyped with the TNGS genotyping platform, and the average mapping rate of captured reads to reference transcripts was 93.14%, of which 53.23% were located in the targeted regions. The TNGS genotyping results were validated by Sanger sequencing and the average consistency rate reached 93.6%. Comprehensive population genetic analysis was performed on the 69 Chinese faba bean accessions and identified four genetic subgroups correlated with the geographic distribution. This study provides valuable genomic resources and a reliable genotyping tool that could be implemented in genetic and molecular breeding studies to accelerate new cultivar development and improvement in faba bean.

Lens orientalis Contributes Quantitative Trait Loci and Candidate Genes Associated With Ascochyta Blight Resistance in Lentil

Australian lentil production is affected by several major biotic constraints including Ascochyta blight (AB), caused by Ascochyta lentis, a devastating fungal disease. Cultivation of AB resistant cultivars, alongside agronomic management including fungicide application, is the current most economically viable control strategy. However, the breakdown of AB resistance in cultivars, such as Northfield and Nipper, suggests the need for introgression of new and diverse resistance genes. Successful introgression entails an understanding of the genetic basis of resistance. In this context, a biparental mapping population derived from a cross between a recently identified AB resistant accession ILWL 180 (Lens orientalis) and a susceptible cultivar ILL 6002 was produced. A genetic linkage map was constructed from single-nucleotide polymorphism markers generated using a genotyping-by-sequencing transcript approach. Genetic dissection of the mapping population revealed a major quantitative trait loci (QTL) region nested with three QTLs on linkage group 5 and explained 9.5-11.5 percent (%) of phenotypic variance for AB resistance. Another QTL was identified on LG2 with phenotypic variance of 9.6%. The identified QTL regions harbored putative candidate genes potentially associated with defense responses to A. lentis infection. The QTL analysis and the candidate gene information are expected to contribute to the development of diagnostic markers and enable marker-assisted resistance selection in lentil breeding programmes.

Reference genome assembly for Australian Ascochyta lentis isolate Al4

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.

Intelligent Characterization of Lentil Genetic Resources: Evolutionary History, Genetic Diversity of Germplasm, and the Need for Well-Represented Collections

The genetic and phenotypic characterization of crops allows us to elucidate their evolutionary and domestication history, the genetic basis of important traits, and the use of variation present in landraces and wild relatives to enhance resilience. In this context, we aim to provide an overview of the main genetic resources developed for lentil and their main outcomes, and to suggest protocols for continued work on this important crop. Lens culinaris is the third-most-important cool-season grain and its use is increasing as a quick-cooking, nutritious, plant-based source of protein. L. culinaris was domesticated in the Fertile Crescent, and six additional wild taxa (L. orientalis, L. tomentosus, L. odemensis, L. lamottei, L. ervoides, and L. nigricans) are recognized. Numerous genetic diversity studies have shown that wild relatives present high levels of genetic variation and provide a reservoir of alleles that can be used for breeding programs. Furthermore, the integration of genetics/genomics and breeding techniques has resulted in identification of quantitative trait loci and genes related to attributes of interest. Genetic maps, massive genotyping, marker-assisted selection, and genomic selection are some of the genetic resources generated and applied in lentil. In addition, despite its size (∼4 Gbp) and complexity, the L. culinaris genome has been assembled, allowing a deeper understanding of its architecture. Still, major knowledge gaps exist in lentil, and a deeper understanding and characterization of germplasm resources, including wild relatives, is critical to lentil breeding and improvement. © 2021 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Recording of lentil seed descriptors Basic Protocol 2: Lentil seed imaging Basic Protocol 3: Lentil seed increase Basic Protocol 4: Recording of primary lentil seed INCREASE descriptors.

Targeted Disruption of Scytalone Dehydratase Gene Using Agrobacterium tumefaciens-Mediated Transformation Leads to Altered Melanin Production in Ascochyta lentis

Sustainable crop production is constantly challenged by the rapid evolution of fungal pathogens equipped with an array of host infection strategies and survival mechanisms. One of the devastating fungal pathogens that infect lentil is the ascomycete Ascochyta lentis which causes black spot or ascochyta blight (AB) on all above ground parts of the plant. In order to explore the mechanisms involved in the pathogenicity of A. lentis, we developed a targeted gene replacement method using Agrobacterium tumefaciens mediated transformation (ATMT) to study and characterize gene function. In this study, we investigated the role of scytalone dehydratase (SCD) in the synthesis of 1,8-dihydroxynaphthalene (DHN)-melanin in AlKewell. Two SCD genes have been identified in AlKewell, AlSCD1 and AlSCD2. Phylogenetic analysis revealed that AlSCD1 clustered with the previously characterized fungal SCDs; thus, AlSCD1 was disrupted using the targeted gene replacement vector, pTAR-hyg-SCD1. The vector was constructed in a single step process using Gibson Assembly, which facilitated an easy and seamless assembly of multiple inserts. The resulting AlKewell scd1::hyg transformants appeared light brown/brownish-pink in contrast to the dark brown pycnidia of the WT strain and ectopic transformant, indicating an altered DHN-melanin production. Disruption of AlSCD1 gene did not result in a change in the virulence profile of AlKewell towards susceptible and resistant lentil varieties. This is the first report of a targeted gene manipulation in A. lentis which serves as a foundation for the functional gene characterization to provide a better understanding of molecular mechanisms involved in pathogen diversity and host specificity.

Characterization of Genetic and Allelic Diversity Amongst Cultivated and Wild Lentil Accessions for Germplasm Enhancement

Intensive breeding of cultivated lentil has resulted in a relatively narrow genetic base, which limits the options to increase crop productivity through selection. Assessment of genetic diversity in the wild gene pool of lentil, as well as characterization of useful and novel alleles/genes that can be introgressed into elite germplasm, presents new opportunities and pathways for germplasm enhancement, followed by successful crop improvement. In the current study, a lentil collection consisting of 467 wild and cultivated accessions that originated from 10 diverse geographical regions was assessed, to understand genetic relationships among different lentil species/subspecies. A total of 422,101 high-confidence SNP markers were identified against the reference lentil genome (cv. CDC Redberry). Phylogenetic analysis clustered the germplasm collection into four groups, namely, Lens culinaris/Lens orientalis, Lens lamottei/Lens odemensis, Lens ervoides, and Lens nigricans. A weak correlation was observed between geographical origin and genetic relationship, except for some accessions of L. culinaris and L. ervoides. Genetic distance matrices revealed a comparable level of variation within the gene pools of L. culinaris (Nei’s coefficient 0.01468–0.71163), L. ervoides (Nei’s coefficient 0.01807–0.71877), and L. nigricans (Nei’s coefficient 0.02188–1.2219). In order to understand any genic differences at species/subspecies level, allele frequencies were calculated from a subset of 263 lentil accessions. Among all cultivated and wild lentil species, L. nigricans exhibited the greatest allelic differentiation across the genome compared to all other species/subspecies. Major differences were observed on six genomic regions with the largest being on Chromosome 1 (c. 1 Mbp). These results indicate that L. nigricans is the most distantly related to L. culinaris and additional structural variations are likely to be identified from genome sequencing studies. This would provide further insights into evolutionary relationships between cultivated and wild lentil germplasm, for germplasm improvement and introgression.

Agrobacterium tumefaciens-mediated transformation and expression of GFP in Ascochyta lentis to characterize ascochyta blight disease progression in lentil

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.

Advanced Imaging for Quantitative Evaluation of Aphanomyces Root Rot Resistance in Lentil

Aphanomyces root rot (ARR) is a soil-borne disease that results in severe yield losses in lentil. The development of resistant cultivars is one of the key strategies to control this pathogen. However, the evaluation of disease severity is limited to visual scores that can be subjective. This study utilized image-based phenotyping approaches to evaluate Aphanomyces euteiches resistance in lentil genotypes in greenhouse (351 genotypes from lentil single plant/LSP derived collection and 191 genotypes from recombinant inbred lines/RIL using digital Red-Green-Blue/RGB and hyperspectral imaging) and field (173 RIL genotypes using unmanned aerial system-based multispectral imaging) conditions. Moderate to strong correlations were observed between RGB, multispectral, and hyperspectral derived features extracted from lentil shoots/roots and visual scores. In general, root features extracted from RGB imaging were found to be strongly associated with disease severity. With only three root traits, elastic net regression model was able to predict disease severity across and within multiple datasets (R 2 = 0.45-0.73 and RMSE = 0.66-1.00). The selected features could represent visual disease scores. Moreover, we developed twelve normalized difference spectral indices (NDSIs) that were significantly correlated with disease scores: two NDSIs for lentil shoot section - computed from wavelengths of 1170, 1160, 1270, and 1280 nm (0.12 ≤ |r| ≤ 0.24, P < 0.05) and ten NDSIs for lentil root sections - computed from wavelengths in the range of 630-670, 700-840, and 1320-1530 nm (0.10 ≤ |r| ≤ 0.50, P < 0.05). Root-derived NDSIs were more accurate in predicting disease scores with an R 2 of 0.54 (RMSE = 0.86), especially when the model was trained and tested on LSP accessions, compared to R 2 of 0.25 (RMSE = 1.64) when LSP and RIL genotypes were used as train and test datasets, respectively. Importantly, NDSIs - computed from wavelengths of 700, 710, 730, and 790 nm - had strong positive correlations with disease scores (0.35 ≤r ≤ 0.50, P < 0.0001), which was confirmed in field phenotyping with similar correlations using vegetation index with red edge wavelength (normalized difference red edge, 0.36 ≤ |r| ≤ 0.57, P < 0.0001). The adopted image-based phenotyping approaches can help plant breeders to objectively quantify ARR resistance and reduce the subjectivity in selecting potential genotypes.

Construction of a high-density interspecific (Lens culinaris x L. odemensis) genetic map based on functional markers for mapping morphological and agronomical traits, and QTLs affecting resistance to Ascochyta in lentil

Usage of high-throughput sequencing approaches allow for the generation and characterization of reference transcriptome datasets that support gene-based marker discovery, which in turn can be used to build genetic maps among other purposes. We have obtained a transcriptome assembly including 49,453 genes for the lentil (Lens culinaris Medik.) cultivar Alpo using RNAseq methodology. This transcriptome was used as reference to obtain 6,306 quality polymorphic markers (SNPs and short indels) analyzing genotype data from a RIL population at F7 generation derived from the interspecific cross between L. culinaris cv. Alpo and L. odemensis accession ILWL235. L. odemensis is a wild species included in the secondary gene pool and can be used as a source for gene introgression in lentil breeding programs. Marker data were used to construct the first genetic interspecific map between these two species. This linkage map has been used to precisely identify regions of the CDC-Redberry lentil draft genome in which the candidate genes for some qualitative traits (seed coat spotting pattern, flower color, and stem pigmentation) could be located. The genome regions corresponding to a significant single quantitative trait locus (QTL) controlling "time to flowering" located in chromosome 6 and three QTLs regulating seed size and positioned in chromosomes 1 and 5 (two QTLs) were also identified. Significant QTLs for Ascochyta blight resistance in lentil were mapped to chromosome 6 in the genome region or close to it where QTLs for Ascochyta blight resistance have previously been reported.