Pathotype-specific genetic factors in chickpea (Cicer arietinum L.) for quantitative resistance to ascochyta blight

Identification of Pathogen-Specific Novel Sources of Genetic Resistance Against Ascochyta Blight and Identification of Their Underlying Genetic Control

Global chickpea production is restricted by Ascochyta blight caused by the necrotrophic fungi Ascochyta rabiei. Developing locally adapted disease-resistant cultivars is an economically and environmentally sustainable approach to combat this disease. However, the lack of genetic variability in cultivated chickpeas and breeder-friendly markers poses a significant challenge to Ascochyta blight-resistant breeding efforts in chickpeas. In this study, we screened the mini-core germplasm of Cicer reticulatum against a local pathotype of A. rabiei. A modified mini-dome screening approach resulted in the identification of five accessions showing a high level of resistance. The mean disease score of resistant accessions ranged between 1.75 ± 0.3 and 2.88 ± 0.4 compared to susceptible accessions, where the mean disease score ranged between 3.59 ± 0.62 and 8.86 ± 0.14. Genome-wide association study revealed a strong association on chromosome 5, explaining ∼58% of the phenotypic variance. The underlying region contained two candidate genes (Cr_14190.1_v2 and Cr_14189.1_v2), the characterization of which showed the presence of a DNA-binding domain (cl28899 and cd18793) in Cr_14190.1_v2 and its orthologs in C. arietinum, whereas Cr_14190.1_v2 carried an additional N-terminal domain (cl31759). qPCR expression analysis in resistant and susceptible accessions revealed ∼3- and ∼110-fold higher transcript abundance for Cr_14189.1 and Cr_14190.1, respectively.

Genome wide association study of genes controlling resistance to Didymella rabiei Pathotype IV through genotyping by sequencing in chickpeas (Cicer arietinum)

Ascochyta blight (AB) is a major disease in chickpeas (Cicer arietinum L.) that can cause a yield loss of up to 100%. Chickpea germplasm collections at the center of origin offer great potential to discover novel sources of resistance to pests and diseases. Herein, 189 Cicer arietinum samples were genotyped via genotyping by sequencing. This chickpea collection was phenotyped for resistance to an aggressive Turkish Didymella rabiei Pathotype IV isolate. Genome-wide association studies based on different models revealed 19 single nucleotide polymorphism (SNP) associations on chromosomes 1, 2, 3, 4, 7, and 8. Although eight of these SNPs have been previously reported, to the best of our knowledge, the remaining ten were associated with AB resistance for the first time. The regions identified in this study can be addressed in future studies to reveal the genetic mechanism underlying AB resistance and can also be utilized in chickpea breeding programs to improve AB resistance in new chickpea varieties.

Genomics-assisted genetics of complex region from chickpea chromosome 4 reveals two candidate genes for Ascochyta blight resistance

Ascochyta blight (AB) disease caused by the fungus Ascochyta rabiei is a major threat to global chickpea production. Molecular breeding for improved AB resistance requires the identification of robust fine-mapped QTLs/candidate genes and associated markers. Earlier, we identified three QTLs (qABR4.1, qABR4.2, and qABR4.3) for AB resistance on chickpea chromosome 4 by employing multiple quantitative trait loci sequencing strategy on an intra-specific (FLIP84-92C x PI359075) and an inter-specific (FLIP84-92C x PI599072) crosses derived recombinant inbred lines. Here, we report the identification of AB resistance providing candidate genes under the fine mapped qABR4.2 and qABR4.3 genomic region by combining genetic mapping, haplotype block inheritance, and expression analysis. The qABR4.2 region was narrowed down from 5.94Mb to ~800kb. Among 34 predicted gene models, a secreted class III peroxidase encoding gene showed higher expression in AB resistant parent after A. rabiei conidia inoculation. Under qABR4.3, we identified a frame-shift mutation in a cyclic nucleotide-gated channel CaCNGC1 gene leading to the truncated N-terminal domain in resistant accession of chickpea. This N-terminal domain of CaCNGC1 interacts with chickpea calmodulin. Thus, our analysis has revealed narrowed genomic regions and their associated polymorphic markers, CaNIP43 and CaCNGCPD1. These co-dominant markers significantly associate with AB resistance on qABR4.2 and qABR4.3 regions. Our genetic analysis revealed that the presence of resistant alleles for two major QTLs (qABR4.1 and qABR4.2) together provide AB resistance in the field while minor QTL qABR4.3 determines the degree of resistance. The identified candidate genes and their diagnostic markers will help in biotechnological and AB resistance introgression into farmers adapted local chickpea varieties.

The R2R3-MYB-SG7 transcription factor CaMYB39 orchestrates surface phenylpropanoid metabolism and pathogen resistance in chickpea

Flavonoids are important plant pigments and defense compounds; understanding the transcriptional regulation of flavonoid biosynthesis may enable engineering crops with improved nutrition and stress tolerance. Here, we characterize R2R3-MYB domain subgroup 7 transcription factor CaMYB39, which regulates flavonol biosynthesis primarily in chickpea trichomes. CaMYB39 overexpression in chickpea was accompanied by a change in flux availability for the phenylpropanoid pathway, particularly flavonol biosynthesis. Lines overexpressing CaMYB39 showed higher isoflavonoid levels, suggesting its role in regulating isoflavonoid pathway. CaMYB39 transactivates the transcription of early flavonoid biosynthetic genes (EBG). FLAVONOL SYNTHASE2, an EBG, encodes an enzyme with higher substrate specificity for dihydrokaempferol than other dihydroflavonols explaining the preferential accumulation of kaempferol derivatives as prominent flavonols in chickpea. Interestingly, CaMYB39 overexpression increased trichome density and enhanced the accumulation of diverse flavonol derivatives in trichome-rich tissues. Moreover, CaMYB39 overexpression reduced reactive oxygen species levels and induced defense gene expression which aids in partially blocking the penetration efficiency of the fungal pathogen, Ascochyta rabiei, resulting in lesser symptoms, thus establishing its role against deadly Ascochyta blight (AB) disease. Overall, our study reports an instance where R2R3-MYB-SG7 member, CaMYB39, besides regulating flavonol biosynthesis, modulates diverse pathways like general phenylpropanoid, isoflavonoid, trichome density, and defense against necrotrophic fungal infection in chickpea.

The nuclear effector ArPEC25 from the necrotrophic fungus Ascochyta rabiei targets the chickpea transcription factor CaβLIM1a and negatively modulates lignin biosynthesis, increasing host susceptibility

Fungal pathogens deploy a barrage of secreted effectors to subvert host immunity, often by evading, disrupting, or altering key components of transcription, defense signaling, and metabolic pathways. However, the underlying mechanisms of effectors and their host targets are largely unexplored in necrotrophic fungal pathogens. Here, we describe the effector protein Ascochyta rabiei PEXEL-like Effector Candidate 25 (ArPEC25), which is secreted by the necrotroph A. rabiei, the causal agent of Ascochyta blight disease in chickpea (Cicer arietinum), and is indispensable for virulence. After entering host cells, ArPEC25 localizes to the nucleus and targets the host LIM transcription factor CaβLIM1a. CaβLIM1a is a transcriptional regulator of CaPAL1, which encodes phenylalanine ammonia lyase (PAL), the regulatory, gatekeeping enzyme of the phenylpropanoid pathway. ArPEC25 inhibits the transactivation of CaβLIM1a by interfering with its DNA-binding ability, resulting in negative regulation of the phenylpropanoid pathway and decreased levels of intermediates of lignin biosynthesis, thereby suppressing lignin production. Our findings illustrate the role of fungal effectors in enhancing virulence by targeting a key defense pathway that leads to the biosynthesis of various secondary metabolites and antifungal compounds. This study provides a template for the study of less explored necrotrophic effectors and their host target functions.

Chickpea (Cicer arietinum L.) Biology and Biotechnology: From Domestication to Biofortification and Biopharming

Chickpea (Cicer arietinum L.), the world's second most consumed legume crop, is cultivated in more than 50 countries around the world. It is a boon for diabetics and is an excellent source of important nutrients such as vitamins A, C, E, K, B1-B3, B5, B6, B9 and minerals (Fe, Zn, Mg and Ca) which all have beneficial effects on human health. By 2050, the world population can cross 9 billion, and in order to feed the teaming millions, chickpea production should also be increased, as it is a healthy alternative to wheat flour and a boon for diabetics. Moreover, it is an important legume that is crucial for food, nutrition, and health security and the livelihood of the small-scale farmers with poor resources, in developing countries. Although marvelous improvement has been made in the development of biotic and abiotic stress-resistant varieties, still there are many lacunae, and to fulfill that, the incorporation of genomic technologies in chickpea breeding (genomics-assisted breeding, high-throughput and precise-phenotyping and implementation of novel breeding strategies) will facilitate the researchers in developing high yielding, climate resilient, water use efficient, salt-tolerant, insect/pathogen resistant varieties, acceptable to farmers, consumers, and industries. This review focuses on the origin and distribution, nutritional profile, genomic studies, and recent updates on crop improvement strategies for combating abiotic and biotic stresses in chickpea.

Ascochyta rabiei: A threat to global chickpea production

The necrotrophic fungus Ascochyta rabiei causes Ascochyta blight (AB) disease in chickpea. A. rabiei infects all aerial parts of the plant, which results in severe yield loss. At present, AB disease occurs in most chickpea-growing countries. Globally increased incidences of A. rabiei infection and the emergence of new aggressive isolates directed the interest of researchers toward understanding the evolution of pathogenic determinants in this fungus. In this review, we summarize the molecular and genetic studies of the pathogen along with approaches that are helping in combating the disease. Possible areas of future research are also suggested.

TAXONOMY

kingdom Mycota, phylum Ascomycota, class Dothideomycetes, subclass Coelomycetes, order Pleosporales, family Didymellaceae, genus Ascochyta, species rabiei.

PRIMARY HOST

A. rabiei survives primarily on Cicer species.

DISEASE SYMPTOMS

A. rabiei infects aboveground parts of the plant including leaves, petioles, stems, pods, and seeds. The disease symptoms first appear as watersoaked lesions on the leaves and stems, which turn brown or dark brown. Early symptoms include small circular necrotic lesions visible on the leaves and oval brown lesions on the stem. At later stages of infection, the lesions may girdle the stem and the region above the girdle falls off. The disease severity increases at the reproductive stage and rounded lesions with concentric rings, due to asexual structures called pycnidia, appear on leaves, stems, and pods. The infected pod becomes blighted and often results in shrivelled and infected seeds.

DISEASE MANAGEMENT STRATEGIES

Crop failures may be avoided by judicious practices of integrated disease management based on the use of resistant or tolerant cultivars and growing chickpea in areas where conditions are least favourable for AB disease development. Use of healthy seeds free of A. rabiei, seed treatments with fungicides, and proper destruction of diseased stubbles can also reduce the fungal inoculum load. Crop rotation with nonhost crops is critical for controlling the disease. Planting moderately resistant cultivars and prudent application of fungicides is also a way to combat AB disease. However, the scarcity of AB-resistant accessions and the continuous evolution of the pathogen challenges the disease management process.

USEFUL WEBSITES

https://www.ndsu.edu/pubweb/pulse-info/resourcespdf/Ascochyta%20blight%20of%20chickpea.pdf https://saskpulse.com/files/newsletters/180531_ascochyta_in_chickpeas-compressed.pdf http://www.pulseaus.com.au/growing-pulses/bmp/chickpea/ascochyta-blight http://agriculture.vic.gov.au/agriculture/pests-diseases-and-weeds/plant-diseases/grains-pulses-and-cereals/ascochyta-blight-of-chickpea http://www.croppro.com.au/crop_disease_manual/ch05s02.php https://www.northernpulse.com/uploads/resources/722/handout-chickpeaascochyta-nov13-2011.pdf http://oar.icrisat.org/184/1/24_2010_IB_no_82_Host_Plant https://www.crop.bayer.com.au/find-crop-solutions/by-pest/diseases/ascochyta-blight.

Conventional and molecular breeding for disease resistance in chickpea: status and strategies

Chickpea (Cicer arietinum L.) is an important grain legume at the global level. Among different biotic stresses, diseases are the most important factor limiting its production, causing yield losses up to 100% in severe condition. The major diseases that adversely affect yield of chickpea include Fusarium wilt, Ascochyta blight and Botrytis gray mold. However, dry root rot, collar rot, Sclerotinia stem rot, rust, stunt disease and phyllody have been noted as emerging biotic threats to chickpea production in many production regions. Identification and incorporation of different morphological and biochemical traits are required through breeding to enhance genetic gain for disease resistance. In recent years, remarkable progress has been made in the development of trait-specific breeding lines, genetic and genomic resources in chickpea. Advances in genomics technologies have opened up new avenues to introgress genes from secondary and tertiary gene pools for improving disease resistance in chickpea. In this review, we have discussed important diseases, constraints and improvement strategies for enhancing disease resistance in chickpea.

Novel Genomic Regions Linked to Ascochyta Blight Resistance in Two Differentially Resistant Cultivars of Chickpea

Ascochyta blight (AB), caused by the fungal pathogen Ascochyta rabiei, is a devastating foliar disease of chickpea (Cicer arietinum L.). The genotyping-by-sequencing (GBS)-based approach was deployed for mapping QTLs associated with AB resistance in chickpea in two recombinant inbred line populations derived from two crosses (AB₃₂₇₉ derived from ILC 1929 × ILC 3279 and AB₄₈₂ derived from ILC 1929 × ILC 482) and tested in six different environments. Twenty-one different genomic regions linked to AB resistance were identified in regions CalG02 and CalG04 in both populations AB₃₂₇₉ and AB₄₈₂. These regions contain 1,118 SNPs significantly associated with AB resistance (p ≤ 0.001), which explained 11.2-39.3% of the phenotypic variation (PVE). Nine of the AB resistance-associated genomic regions were newly detected in this study, while twelve regions were known from previous AB studies. The proposed physical map narrows down AB resistance to consistent genomic regions identified across different environments. Gene ontology (GO) assigned these QTLs to 319 genes, many of which were associated with stress and disease resistance, and with most important genes belonging to resistance gene families such as leucine-rich repeat (LRR) and transcription factor families. Our results indicate that the flowering-associated gene GIGANTEA is a possible key factor in AB resistance in chickpea. The results have identified AB resistance-associated regions on the physical genetic map of chickpea and allowed for the identification of associated markers that will help in breeding of AB-resistant varieties.

Breeding and Genomics Interventions for Developing Ascochyta Blight Resistant Grain Legumes

Grain legumes are a key food source for ensuring global food security and sustaining agriculture. However, grain legume production is challenged by growing disease incidence due to global climate change. Ascochyta blight (AB) is a major disease, causing substantial yield losses in grain legumes worldwide. Harnessing the untapped reserve of global grain legume germplasm, landraces, and crop wild relatives (CWRs) could help minimize yield losses caused by AB infection in grain legumes. Several genetic determinants controlling AB resistance in various grain legumes have been identified following classical genetic and conventional breeding approaches. However, the advent of molecular markers, biparental quantitative trait loci (QTL) mapping, genome-wide association studies, genomic resources developed from various genome sequence assemblies, and whole-genome resequencing of global germplasm has revealed AB-resistant gene(s)/QTL/genomic regions/haplotypes on various linkage groups. These genomics resources allow plant breeders to embrace genomics-assisted selection for developing/transferring AB-resistant genomic regions to elite cultivars with great precision. Likewise, advances in functional genomics, especially transcriptomics and proteomics, have assisted in discovering possible candidate gene(s) and proteins and the underlying molecular mechanisms of AB resistance in various grain legumes. We discuss how emerging cutting-edge next-generation breeding tools, such as rapid generation advancement, field-based high-throughput phenotyping tools, genomic selection, and CRISPR/Cas9, could be used for fast-tracking AB-resistant grain legumes to meet the increasing demand for grain legume-based protein diets and thus ensuring global food security.

A kabuli chickpea ideotype

The concept of 'crop ideotype' is coined as a desirable plant model expected to better perform for seed yield, oils and other useful characteristics when developed as a cultivar, and it consists of two major approaches, namely, (i) 'defect elimination', that is, integration of disease resistance to a susceptible genotype from a resistant genotype and (ii) 'selection for yield' by improving yield after crosses between desirable parents. For consideration of these approaches, here we introduced an ideotype in kabuli chickpea (Cicer arietinum L.) which is high-yielding, extra-large-seeded, and double- or multi-podded, has high plant height and imparipinnate-leafed traits, and is heat tolerant and resistant to ascochyta blight [Ascochyta rabiei (Pass.) Labr.], which causes considerable yield losses, via marker-assisted selection. F3 and F4 lines were evaluated for agro-morphological traits divided into six classes, namely, (i) imparipinnate-leafed and single-podded progeny, (ii) imparipinnate-leafed and double-podded progeny, (iii) imparipinnate-leafed and multi-podded progeny, (iv) unifoliolate-leafed and single-podded progeny, (v) unifoliolate-leafed and double-podded progeny, (vi) unifoliolate-leafed and multi-podded progeny. F3:4 lines having 100-seed weight ≥ 45 g and double- or multi-podded traits were additionally assessed for resistance to ascochyta blight using molecular markers including SCY17590 and CaETR-1. Superior lines having higher values than their best parents were determined for all studied traits indicating that economic and important traits including yield and seed size in chickpea could be improved by crossing suitable parents. Imparipinnate-leafed and multi-podded plants had not only the highest number of pods and seeds per plant but also the highest yield. On the other hand, imparipinnate-leafed and single podded progeny had the largest seed size, followed by imparipinnate-leafed and double-podded progeny. Multi-podded plants produced 23% more seed yield than that of single-podded plants, while multi-podded plants attained 7.6% more seed yield than that of double-podded plants. SCY17590 and CaETR-1 markers located on LG4 related to QTLAR2 and QTLAR1 were found in 14 lines among 152 F3:4 lines. Six superior lines were selected for being double- or multi-podded, imparipinnate-leafed, suitable for combine harvest, heat-tolerant, and resistant to ascochyta blight, and having both of two resistance markers and extra-large seeds as high as 50-60 g per 100-seed weight. Resistance alleles from two different backgrounds for resistance to ascochyta blight were integrated with double- or multi-podded kabuli chickpea lines having high yield, extra-large seeds, high plant height, imparipinnate-leaves and high heat tolerance, playing a crucial role for future demands of population and food security. These approaches seem to be applicable in ideotype breeding for other important crop plants.

Genome-Wide Association Analyses Track Genomic Regions for Resistance to Ascochyta rabiei in Australian Chickpea Breeding Germplasm

Ascochyta blight (AB), caused by a necrotrophic fungus, Ascochyta rabiei (syn. Phoma rabiei) has the potential to destroy the chickpea industry worldwide, due to limited sources of genetic resistance in the cultivated gene pool, high evolutionary potential of the pathogen and challenges with integrated disease management. Therefore, the deployment of stable genetic resistance in new cultivars could provide an effective disease control strategy. To investigate the genetic basis of AB resistance, genotyping-by-sequencing based DArTseq-single nucleotide polymorphism (SNP) marker data along with phenotypic data of 251 advanced breeding lines and chickpea cultivars were used to perform genome-wide association (GWAS) analysis. Host resistance was evaluated seven weeks after sowing using two highly aggressive single spore isolates (F17191-1 and TR9571) of A. rabiei. GWAS analyses based on single-locus and multi-locus mixed models and haplotyping trend regression identified twenty-six genomic regions on Ca1, Ca4, and Ca6 that showed significant association with resistance to AB. Two haplotype blocks (HB) on chromosome Ca1; HB5 (992178-1108145 bp), and HB8 (1886221-1976301 bp) were associated with resistance against both isolates. Nine HB on the chromosome, Ca4, spanning a large genomic region (14.9-56.6 Mbp) were also associated with resistance, confirming the role of this chromosome in providing resistance to AB. Furthermore, trait-marker associations in two F₃ derived populations for resistance to TR9571 isolate at the seedling stage under glasshouse conditions were also validated. Eighty-nine significantly associated SNPs were located within candidate genes, including genes encoding for serine/threonine-protein kinase, Myb protein, quinone oxidoreductase, and calmodulin-binding protein all of which are implicated in disease resistance. Taken together, this study identifies valuable sources of genetic resistance, SNP markers and candidate genes underlying genomic regions associated with AB resistance which may enable chickpea breeding programs to make genetic gains via marker-assisted/genomic selection strategies.

Effects of Host Plant Resistance and Fungicide Applications on Ascochyta Blight Symptomology and Yield of Chickpea

Ascochyta blight (AB), caused by the pathogen Ascochyta rabiei, is a major threat to chickpea production worldwide, causing major yield losses and decreasing quality. Control of AB requires integrating pest management options including resistant cultivars and fungicide applications. To address this, fungicides with different modes of action were evaluated on three chickpea cultivars with differing levels of susceptibility to AB under irrigated and dryland conditions in 2015 to 2017. The fungicides were applied once or twice and compared with a no-fungicide application control on AB score and yield. The mean grain yields across locations and years were 1,753, 1,283, and 981 kg/ha, with a corresponding AB mean score of 2.6, 3.2, and 3.3 on 0 to 7 scale (where 0 is no disease and 7 is completely dead) for the moderately resistant, moderately susceptible, and susceptible chickpea cultivars, respectively. Fungicide application was not enough to control disease throughout the season. The use of AB-resistant cultivars had the most significant impact on minimizing the disease and maximizing yield, irrespective of year and location. This study supports previous research indicating that planting AB-resistant chickpea cultivars is essential for disease control, regardless of the fungicides applied.

Molecular mapping of quantitative trait loci for ascochyta blight and botrytis grey mould resistance in an inter-specific cross in chickpea (Cicer arietinum L.) using genotyping by sequencing

Ascochyta blight (AB) and botrytis grey mould (BGM) are the most devastating fungal diseases of chickpea worldwide. The wild relative of chickpea, C. reticulatum acc. ILWC 292 was found resistant to BGM whereas, GPF2 (Cicer arietinum L.) is resistant to AB. A total of 187 F₈ Recombinant Inbred Lines (RILs) developed from an inter-specific cross of GPF2 × C. reticulatum acc. ILWC 292 were used to identify quantitative trait loci (QTLs) responsible for resistance to AB and BGM. RILs along with parents were evaluated under artificial epiphytotic field/laboratory conditions for two years. Highly significant differences (P < 0.001) were observed for reaction to both pathogens in both years. Parents and RILs were genotyped-by-sequencing to identify genome wide single nucleotide polymorphism (SNPs). A total of 1365 filtered and parental polymorphic SNPs were used for linkage map construction, of which, 673 SNPs were arranged on eight linkage groups. Composite interval mapping revealed three QTLs for AB and four QTLs for BGM resistance. Out of which, two QTLs for AB and three QTLs for BGM were consistent in both years. These QTLs can be targeted for further fine mapping for deployment of resistance to AB and BGM in elite chickpea cultivars using marker-assisted-selection.

Identification of Novel Sources of Resistance to Ascochyta Blight in a Collection of Wild Cicer Accessions

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.

Genome-wide SNP discovery for development of high-density genetic map and QTL mapping of ascochyta blight resistance in chickpea (Cicer arietinum L.)

A high-density linkage map of chickpea using 3430 SNPs was constructed and used to identify QTLs and candidate genes for ascochyta blight resistance in chickpea. Chickpea cultivation in temperate conditions is highly vulnerable to ascochyta blight infection. Cultivation of resistant cultivars in combination with fungicide application within an informed disease management package is the most effective method to control ascochyta blight in chickpeas. Identifying new sources of resistance is critical for continued improvement in ascochyta blight resistance in chickpea. The objective of this study was to identify genetic loci and candidate genes controlling the resistance to ascochyta blight in recombinant inbred lines derived from crossing cultivars Amit and ICCV 96029. The RILs were genotyped using the genotyping-by-sequencing procedure and Illumina^® GoldenGate array. The RILs were evaluated in the field over three site-years and in three independent greenhouse experiments. A genetic map with eight linkage groups was constructed using 3430 SNPs. Eight QTLs for resistance were identified on chromosomes 2, 3, 4, 5 and 6. The QTLs individually explained 7-40% of the phenotypic variations. The QTLs on chromosomes 2 and 6 were associated with the resistance at vegetative stage only. The QTLs on chromosomes 2 and 4 that were previously reported to be conserved across diverse genetic backgrounds and against different isolates of Ascochyta rabiei were confirmed in this study. Candidate genes were identified within the QTL regions. Their co-localization with the underlying QTLs was confirmed by genetic mapping. The candidate gene-based SNP markers would lead to more efficient marker-assisted selection for ascochyta blight resistance and would provide a framework for fine mapping and subsequent cloning of the genes associated with the resistance.

QTL sequencing strategy to map genomic regions associated with resistance to ascochyta blight in chickpea

Whole-genome sequencing-based bulked segregant analysis (BSA) for mapping quantitative trait loci (QTL) provides an efficient alternative approach to conventional QTL analysis as it significantly reduces the scale and cost of analysis with comparable power to QTL detection using full mapping population. We tested the application of next-generation sequencing (NGS)-based BSA approach for mapping QTLs for ascochyta blight resistance in chickpea using two recombinant inbred line populations CPR-01 and CPR-02. Eleven QTLs in CPR-01 and six QTLs in CPR-02 populations were mapped on chromosomes Ca1, Ca2, Ca4, Ca6 and Ca7. The QTLs identified in CPR-01 using conventional biparental mapping approach were used to compare the efficiency of NGS-based BSA in detecting QTLs for ascochyta blight resistance. The QTLs on chromosomes Ca1, Ca4, Ca6 and Ca7 overlapped with the QTLs previously detected in CPR-01 using conventional QTL mapping method. The QTLs on chromosome Ca4 were detected in both populations and overlapped with the previously reported QTLs indicating conserved region for ascochyta blight resistance across different chickpea genotypes. Six candidate genes in the QTL regions identified using NGS-based BSA on chromosomes Ca2 and Ca4 were validated for their association with ascochyta blight resistance in the CPR-02 population. This study demonstrated the efficiency of NGS-based BSA as a rapid and cost-effective method to identify QTLs associated with ascochyta blight in chickpea.

mQTL-seq and classical mapping implicates the role of an AT-HOOK MOTIF CONTAINING NUCLEAR LOCALIZED (AHL) family gene in Ascochyta blight resistance of chickpea

Ascochyta blight (AB) caused by the fungal pathogen Ascochyta rabiei is a serious foliar disease of chickpea (Cicer arietinum L.). Despite many genetic studies on chickpea-Ascochyta interaction, genome-wide scan of chickpea for the identification of AB-associated quantitative trait loci (QTLs) and their gene(s) has not been accomplished. To elucidate narrow QTLs for AB resistance, here, we report the use of multiple QTL-sequencing approach on 2 sets of extreme AB phenotype bulks derived from Cicer intraspecific and interspecific crosses. Two major QTLs, qABR4.1 and qABR4.2, and a minor QTL, qABR4.3, were identified on assembled chickpea pseudomolecule 4. We narrowed qABR4.1 to a "robust region" at 4.568-4.618 Mb through mapping on a larger intraspecific cross-derived population and comparative analysis. Among 4 genes, the CaAHL18 gene showed higher expression under Ascochyta stress in AB resistant parent suggesting that it is the candidate gene under "robust qABR4.1." Dual-luciferase assay with CaAHL18 polymorphic cis-regulatory sequences showed that allelic variation is associated with higher expression. Thus, our findings on chickpea-Ascochyta interaction have narrowed down AB resistance associated QTLs on chickpea physical map. The narrowed QTLs and gene-associated markers will help in biotechnological and breeding programs for chickpea improvement.

Current advances in chickpea genomics: applications and future perspectives

Chickpea genomics promises to illuminate our understanding of genome organization, structural variations, evolutionary and domestication-related insights and fundamental biology of legume crops. Unprecedented advancements of next generation sequencing (NGS) technologies have enabled in decoding of multiple chickpea genome sequences and generating huge genomic resources in chickpea both at functional and structural level. This review is aimed to update the current progress of chickpea genomics ranging from high density linkage map development, genome-wide association studies (GWAS), functional genomics resources for various traits, emerging role of abiotic stress responsive coding and non-coding RNAs after the completion of draft chickpea genome sequences. Additionally, the current efforts of whole genome re-sequencing (WGRS) approach of global chickpea germplasm to capture the global genetic diversity existing in the historically released varieties across the world and increasing the resolution of the previously identified candidate gene(s) of breeding importance have been discussed. Thus, the outcomes of these genomics resources will assist in genomics-assisted selection and facilitate breeding of climate-resilient chickpea cultivars for sustainable agriculture.

Evidence and Consequence of a Highly Adapted Clonal Haplotype within the Australian Ascochyta rabiei Population

The Australian Ascochyta rabiei (Pass.) Labr. (syn. Phoma rabiei) population has low genotypic diversity with only one mating type detected to date, potentially precluding substantial evolution through recombination. However, a large diversity in aggressiveness exists. In an effort to better understand the risk from selective adaptation to currently used resistance sources and chemical control strategies, the population was examined in detail. For this, a total of 598 isolates were quasi-hierarchically sampled between 2013 and 2015 across all major Australian chickpea growing regions and commonly grown host genotypes. Although a large number of haplotypes were identified (66) through short sequence repeat (SSR) genotyping, overall low gene diversity (Hexp = 0.066) and genotypic diversity (D = 0.57) was detected. Almost 70% of the isolates assessed were of a single dominant haplotype (ARH01). Disease screening on a differential host set, including three commonly deployed resistance sources, revealed distinct aggressiveness among the isolates, with 17% of all isolates identified as highly aggressive. Almost 75% of these were of the ARH01 haplotype. A similar pattern was observed at the host level, with 46% of all isolates collected from the commonly grown host genotype Genesis090 (classified as "resistant" during the term of collection) identified as highly aggressive. Of these, 63% belonged to the ARH01 haplotype. In conclusion, the ARH01 haplotype represents a significant risk to the Australian chickpea industry, being not only widely adapted to the diverse agro-geographical environments of the Australian chickpea growing regions, but also containing a disproportionately large number of aggressive isolates, indicating fitness to survive and replicate on the best resistance sources in the Australian germplasm.

Genome Analysis Identified Novel Candidate Genes for Ascochyta Blight Resistance in Chickpea Using Whole Genome Re-sequencing Data

Ascochyta blight (AB) is a fungal disease that can significantly reduce chickpea production in Australia and other regions of the world. In this study, 69 chickpea genotypes were sequenced using whole genome re-sequencing (WGRS) methods. They included 48 Australian varieties differing in their resistance ranking to AB, 16 advanced breeding lines from the Australian chickpea breeding program, four landraces, and one accession representing the wild chickpea species Cicer reticulatum. More than 800,000 single nucleotide polymorphisms (SNPs) were identified. Population structure analysis revealed relatively narrow genetic diversity amongst recently released Australian varieties and two groups of varieties separated by the level of AB resistance. Several regions of the chickpea genome were under positive selection based on Tajima's D test. Both Fst genome- scan and genome-wide association studies (GWAS) identified a 100 kb region (AB4.1) on chromosome 4 that was significantly associated with AB resistance. The AB4.1 region co-located to a large QTL interval of 7 Mb∼30 Mb identified previously in three different mapping populations which were genotyped at relatively low density with SSR or SNP markers. The AB4.1 region was validated by GWAS in an additional collection of 132 advanced breeding lines from the Australian chickpea breeding program, genotyped with approximately 144,000 SNPs. The reduced level of nucleotide diversity and long extent of linkage disequilibrium also suggested the AB4.1 region may have gone through selective sweeps probably caused by selection of the AB resistance trait in breeding. In total, 12 predicted genes were located in the AB4.1 QTL region, including those annotated as: NBS-LRR receptor-like kinase, wall-associated kinase, zinc finger protein, and serine/threonine protein kinases. One significant SNP located in the conserved catalytic domain of a NBS-LRR receptor-like kinase led to amino acid substitution. Transcriptional analysis using qPCR showed that some predicted genes were significantly induced in resistant lines after inoculation compared to non-inoculated plants. This study demonstrates the power of combining WGRS data with relatively simple traits to rapidly develop "functional makers" for marker-assisted selection and genomic selection.

Genetic Analysis of NBS-LRR Gene Family in Chickpea and Their Expression Profiles in Response to Ascochyta Blight Infection

Ascochyta blight is one of the major diseases of chickpea worldwide. The genetic resistance to ascochyta blight in chickpea is complex and governed by multiple QTLs. However, the molecular mechanism of quantitative disease resistance to ascochyta blight and the genes underlying these QTLs are still unknown. Most often disease resistance is determined by resistance (R) genes. The most predominant R-genes contain nucleotide binding site and leucine rich repeat (NBS-LRR) domains. A total of 121 NBS-LRR genes were identified in the chickpea genome. Ninety-eight of these genes contained all essential conserved domains while 23 genes were truncated. The NBS-LRR genes were grouped into eight distinct classes based on their domain architecture. Phylogenetic analysis grouped these genes into two major clusters based on their structural variation, the first cluster with toll or interleukin-1 like receptor (TIR) domain and the second cluster either with or without a coiled-coil domain. The NBS-LRR genes are distributed unevenly across the eight chickpea chromosomes and nearly 50% of the genes are present in clusters. Thirty of the NBS-LRR genes were co-localized with nine of the previously reported ascochyta blight QTLs and were tested as potential candidate genes for ascochyta blight resistance. Expression pattern of these genes was studied in two resistant (CDC Corinne and CDC Luna) and one susceptible (ICCV 96029) genotypes at different time points after ascochyta blight infection using real-time quantitative PCR. Twenty-seven NBS-LRR genes showed differential expression in response to ascochyta blight infection in at least one genotype at one time point. Among these 27 genes, the majority of the NBS-LRR genes showed differential expression after inoculation in both resistant and susceptible genotypes which indicates the involvement of these genes in response to ascochyta blight infection. Five NBS-LRR genes showed genotype specific expression. Our study provides a new insight of NBS-LRR gene family in chickpea and the potential involvement of NBS-LRR genes in response to ascochyta blight infection.

A Genome-wide Combinatorial Strategy Dissects Complex Genetic Architecture of Seed Coat Color in Chickpea

The study identified 9045 high-quality SNPs employing both genome-wide GBS- and candidate gene-based SNP genotyping assays in 172, including 93 cultivated (desi and kabuli) and 79 wild chickpea accessions. The GWAS in a structured population of 93 sequenced accessions detected 15 major genomic loci exhibiting significant association with seed coat color. Five seed color-associated major genomic loci underlying robust QTLs mapped on a high-density intra-specific genetic linkage map were validated by QTL mapping. The integration of association and QTL mapping with gene haplotype-specific LD mapping and transcript profiling identified novel allelic variants (non-synonymous SNPs) and haplotypes in a MATE secondary transporter gene regulating light/yellow brown and beige seed coat color differentiation in chickpea. The down-regulation and decreased transcript expression of beige seed coat color-associated MATE gene haplotype was correlated with reduced proanthocyanidins accumulation in the mature seed coats of beige than light/yellow brown seed colored desi and kabuli accessions for their coloration/pigmentation. This seed color-regulating MATE gene revealed strong purifying selection pressure primarily in LB/YB seed colored desi and wild Cicer reticulatum accessions compared with the BE seed colored kabuli accessions. The functionally relevant molecular tags identified have potential to decipher the complex transcriptional regulatory gene function of seed coat coloration and for understanding the selective sweep-based seed color trait evolutionary pattern in cultivated and wild accessions during chickpea domestication. The genome-wide integrated approach employed will expedite marker-assisted genetic enhancement for developing cultivars with desirable seed coat color types in chickpea.

Genetic diversity and association mapping of iron and zinc concentrations in chickpea (Cicer arietinum L.)

Chickpea (Cicer arietinum L.) is the world's second most important pulse crop after common bean. Chickpea has historically been an important daily staple in the diet of millions of people, especially in the developing countries. Current chickpea breeding programs have mainly been directed toward high yield, biotic and abiotic stress resilience that has increased global production, but less attention has been directed toward improving micronutrient concentrations in seeds. In an effort to develop micronutrient-dense chickpea lines, a study to examine the variability and to identify SNP alleles associated with seed iron and zinc concentrations was conducted using 94 diverse accessions of chickpea. The results indicated that there is substantial variability present in chickpea germplasm for seed iron and zinc concentrations. In the current set of germplasm, zinc is negatively correlated with grain yield across all locations and years; whereas the negative correlation between iron and grain yield was only significant at the Elrose locality. Eight SNP loci associated with iron and (or) zinc concentrations in chickpea seeds were identified. One SNP located on chromosome 1 (chr1) is associated with both iron and zinc concentrations. On chr4, three SNPs associated with zinc concentration and two SNPs for iron concentration were identified. Two additional SNP loci, one on chr6 and the other on chr7, were also found to be associated with iron and zinc concentrations, respectively. The results show potential opportunity for molecular breeding for improvement of seed iron and zinc concentrations in chickpea.

A combinatorial approach of comprehensive QTL-based comparative genome mapping and transcript profiling identified a seed weight-regulating candidate gene in chickpea

High experimental validation/genotyping success rate (94-96%) and intra-specific polymorphic potential (82-96%) of 1536 SNP and 472 SSR markers showing in silico polymorphism between desi ICC 4958 and kabuli ICC 12968 chickpea was obtained in a 190 mapping population (ICC 4958 × ICC 12968) and 92 diverse desi and kabuli genotypes. A high-density 2001 marker-based intra-specific genetic linkage map comprising of eight LGs constructed is comparatively much saturated (mean map-density: 0.94 cM) in contrast to existing intra-specific genetic maps in chickpea. Fifteen robust QTLs (PVE: 8.8-25.8% with LOD: 7.0-13.8) associated with pod and seed number/plant (PN and SN) and 100 seed weight (SW) were identified and mapped on 10 major genomic regions of eight LGs. One of 126.8 kb major genomic region harbouring a strong SW-associated robust QTL (Caq'SW1.1: 169.1-171.3 cM) has been delineated by integrating high-resolution QTL mapping with comprehensive marker-based comparative genome mapping and differential expression profiling. This identified one potential regulatory SNP (G/A) in the cis-acting element of candidate ERF (ethylene responsive factor) TF (transcription factor) gene governing seed weight in chickpea. The functionally relevant molecular tags identified have potential to be utilized for marker-assisted genetic improvement of chickpea.

Fast track genetic improvement of ascochyta blight resistance and double podding in chickpea by marker-assisted backcrossing

Ascochyta blight (AB) caused by the fungus Ascochyta rabiei Pass. Lab. is one of the major diseases of chickpea worldwide and a constraint to production in western Canada. The use of varieties with high levels of resistance is considered the most economical solution for long-term ascochyta blight management in chickpea. QTL for resistance to ascochyta blight have been identified in chickpea. The availability of molecular markers associated with QTL for ascochyta blight resistant and double podding provides an opportunity to apply marker-assisted backcrossing to introgress the traits into adapted chickpea cultivars. In the present study, molecular markers that were linked to the QTL for ascochyta blight resistance and the double podding trait, and those unlinked to the resistance were used in foreground and background selection, respectively, in backcrosses between moderately resistant donors (CDC Frontier and CDC 425-14) and the adapted varieties (CDC Xena, CDC Leader and FLIP98-135C). The strategy included two backcrosses and selection for two QTL for ascochyta blight resistance and a locus associated with double podding. The fixation of the elite genetic background was monitored with 16-22 SSR markers to accelerate restoration of the genetic background at each backcross. By the BC2F1 generation, plants with improved ascochyta blight resistance and double podding were identified. The selected plants possessed the majority of elite parental type SSR alleles on all fragments analyzed except the segment of LG 4, LG 6 and LG 8 that possessed the target QTL. The results showed that the adapted variety could be efficiently converted into a variety with improved resistance in two backcross generations.

Multi-environment QTL analyses for drought-related traits in a recombinant inbred population of chickpea (Cicer arietinum L.)

A recombinant inbred line (RIL) population, comprising 181 lines derived from ILC588 × ILC3279, was evaluated in 10 environments across three locations with different moisture gradients. A drought resistance score (DRS) and three phenology traits-plant height (PLHT), days to flowering (DFLR), and days to maturity (MAT)-were recorded along with seven yield-related traits-grain yield (GY), biological yield (BY), harvest index (HI), the number of pods/3 plants (Pod), percentage of empty pods (%Epod), 100 seed weight (100 sw), and seed number/3 plants (SN). Two RILs (152, 162) showed the best GYs and DRSs under stressed and non-stressed environments. The quantitative trait loci (QTLs) analyses detected 93 significant QTLs (LOD ≥ 2.0) across the genome × environment interactions. The highest phenotypic variation (>24 %) was explained by the QTLDFLR in Terbol-11. Four common possible pleiotropic QTLs on LG3 and LG4 were identified as associated with DFLR, DRS, GY, MAT, HI, SN, and Pod. No significant epistatic interactions were found between these QTLs and the other markers. However, the QTL for DRS was detected as a conserved QTL in three late planting environments. The markers H6C-07 (on LG3) and H5G01 (on LG4) were associated with QTLs for many traits in all environments studied except two. The allele 'A' of marker H6C07 (from the tolerant parent ILC588) explained 80 % of the yield increase under late planting and 29.8 % of that under dry environments. Concentrating on LG3 and LG4 in molecular breeding programs for drought could speed up improvement for these traits.

Characterization and genetic analysis of an EIN4-like sequence (CaETR-1) located in QTL(AR1) implicated in ascochyta blight resistance in chickpea

Two alleles of a chickpea (Cicer arietinum L.) ethylene receptor-like sequence (CaETR-1) were sequence-characterized using synteny analysis with genome sequences of Medicago truncatula L. The full length of the sequence obtained in the accession FLIP84-92C resistant to ascochyta blight (CaETR-1a) span 4,428 bp, including the polyadenylation signal in the 3'-untranslated region (UTR), whereas it has a 730 bp deletion in the 3'-UTR region in the susceptible accession PI359075 (CaETR-1b). The deduced protein belongs to subfamily II of the ethylene receptors and contains all the domains that define EIN4 homologs in Arabidopsis. The EIN4-like sequence (CaETR-1) has been mapped using a recombinant inbred line (RIL) population derived from an intraspecific cross between ILC3279 and WR315, resistant and susceptible to blight, respectively. The locus was located in LGIVa of the genetic map, flanked by markers NCPGR91 and GAA47 (at distances of 11.3 and 17.9 cM, respectively). This is the first potentially functional sequence identified under a QTL peak for ascochyta blight resistance in chickpea (QTL(AR1)). This EIN4-like (CaETR-1) sequence explained up to 33.8% of the total phenotypic variation. This sequence could be directly related to blight resistance, together with other QTLs that have been found to be involved in resistance to this major chickpea disease.

EST-derived genic molecular markers: development and utilization for generating an advanced transcript map of chickpea

Well-saturated linkage maps especially those based on expressed sequence tag (EST)-derived genic molecular markers (GMMs) are a pre-requisite for molecular breeding. This is especially true in important legumes such as chickpea where few simple sequence repeats (SSR) and even fewer GMM-based maps have been developed. Therefore, in this study, 2,496 ESTs were generated from chickpea seeds and utilized for the development of 487 novel EST-derived functional markers which included 125 EST-SSRs, 151 intron targeted primers (ITPs), 109 expressed sequence tag polymorphisms (ESTPs), and 102 single nucleotide polymorphisms (SNPs). Whereas ESTSSRs, ITPs, and ESTPs were developed by in silico analysis of the developed EST sequences, SNPs were identified by allele resequencing and their genotyping was performedusing the Illumina GoldenGate Assay. Parental polymorphism was analyzed between C. arietinum ICC4958 and C. reticulatum PI489777, parents of the reference chickpea mapping population, using a total of 872 markers: 487 new gene-based markers developed in this study along with 385 previously published markers, of which 318 (36.5%) were found to be polymorphic and were used for genotyping. The genotypic data were integrated with the previously published data of 108 markers and an advanced linkage map was generated that contained 406 loci distributed on eight linkage groups that spanned 1,497.7 cM. The average marker density was 3.68 cM and the average number of markers per LG was 50.8. Among the mapped markers, 303 new genomic locations were defined that included 177 gene-based and 126 gSSRs (genomic SSRs) thereby producing the most advanced gene-rich map of chickpea solely based on co-dominant markers.

Future prospects for ascochyta blight resistance breeding in cool season food legumes

Legume cultivation is strongly hampered by the occurrence of ascochyta blights. Strategies of control have been developed but only marginal successes achieved. Breeding for disease resistance is regarded the most cost efficient method of control. Significant genetic variation for disease resistance exists in most legume crops with numerous germplasm lines maintained, providing an excellent resource for plant breeders. Fast and reliable screening methods have been adjusted to fulfill breeding program needs. However, the complex inheritance controlled quantitatively by multiple genes, has been difficult to manipulate. Successful application of biotechnology to ascochyta blight resistance breeding in legume crops will facilitate a good biological knowledge both of the crops-pathogen interaction and of the mechanisms underlying resistance. The current focus in applied breeding is leveraging biotechnological tools to develop more and better markers to speed up the delivery of improved cultivars to the farmer. To date, however, progress in marker development and delivery of useful markers has been slow in most legumes. The limited saturation of the genomic regions bearing putative QTLs in legume crops makes difficult to identify the most tightly linked markers and to determine the accurate position of QTLs. The application of next generation sequencing technologies will contribute to the development of new markers and the identification of candidate genes for ascochyta blight resistance.

Bacterial artificial chromosome libraries of pulse crops: characteristics and applications

Pulse crops are considered minor on a global scale despite their nutritional value for human consumption. Therefore, they are relatively less extensively studied in comparison with the major crops. The need to improve pulse crop production and quality will increase with the increasing global demand for food security and people's awareness of nutritious food. The improvement of pulse crops will require fully utilizing all their genetic resources. Bacterial artificial chromosome (BAC) libraries of pulse crops are essential genomic resources that have the potential to accelerate gene discovery and enhance molecular breeding in these crops. Here, we review the availability, characteristics, applications, and potential applications of the BAC libraries of pulse crops.

Advancing the STMS genomic resources for defining new locations on the intraspecific genetic linkage map of chickpea (Cicer arietinum L.)

Background

Chickpea (Cicer arietinum L.) is an economically important cool season grain legume crop that is valued for its nutritive seeds having high protein content. However, several biotic and abiotic stresses and the low genetic variability in the chickpea genome have continuously hindered the chickpea molecular breeding programs. STMS (Sequence Tagged Microsatellite Sites) markers which are preferred for the construction of saturated linkage maps in several crop species, have also emerged as the most efficient and reliable source for detecting allelic diversity in chickpea. However, the number of STMS markers reported in chickpea is still limited and moreover exhibit low rates of both inter and intraspecific polymorphism, thereby limiting the positions of the SSR markers especially on the intraspecific linkage maps of chickpea. Hence, this study was undertaken with the aim of developing additional STMS markers and utilizing them for advancing the genetic linkage map of chickpea which would have applications in QTL identification, MAS and for de novo assembly of high throughput whole genome sequence data.

Results

A microsatellite enriched library of chickpea (enriched for (GT/CA)_n and (GA/CT)_n repeats) was constructed from which 387 putative microsatellite containing clones were identified. From these, 254 STMS primers were designed of which 181 were developed as functional markers. An intraspecific mapping population of chickpea, [ICCV-2 (single podded) × JG-62 (double podded)] and comprising of 126 RILs, was genotyped for mapping. Of the 522 chickpea STMS markers (including the double-podding trait, screened for parental polymorphism, 226 (43.3%) were polymorphic in the parents and were used to genotype the RILs. At a LOD score of 3.5, eight linkage groups defining the position of 138 markers were obtained that spanned 630.9 cM with an average marker density of 4.57 cM. Further, based on the common loci present between the current map and the previously published chickpea intraspecific map, integration of maps was performed which revealed improvement of marker density and saturation of the region in the vicinity of sfl (double-podding) gene thereby bringing about an advancement of the current map.

Conclusion

An arsenal of 181 new chickpea STMS markers was reported. The developed intraspecific linkage map defined map positions of 138 markers which included 101 new locations.Map integration with a previously published map was carried out which revealed an advanced map with improved density. This study is a major contribution towards providing advanced genomic resources which will facilitate chickpea geneticists and molecular breeders in developing superior genotypes with improved traits.

QTL mapping of internal heat necrosis in tetraploid potato

Internal heat necrosis (IHN) is a physiological disorder of potato tubers. We developed a linkage map of tetraploid potato using AFLP and SSR markers, and mapped QTL for mean severity and percent incidence of IHN. Phenotypic data indicated that the distribution of IHN is skewed toward resistance. Late foliage maturity was slightly but significantly correlated with increased IHN symptoms. The linkage map for 'Atlantic', the IHN-susceptible parent, covered 1034.4 cM and included 13 linkage groups, and the map for B1829-5, the IHN-resistant parent, covered 940.2 cM and contained 14 linkage groups. QTL for increased resistance to IHN were located on chromosomes IV, V, and groups VII and X of 'Atlantic', and on group VII of B1829-5 in at least 2 of 3 years. The QTL explained between 4.5 and 29.4% of the variation for mean severity, and from 3.7 to 14.5% of the variation for percent incidence. Most QTL detected were dominant, and associated with decreased IHN symptoms. One SSR and 13 AFLP markers that were linked to IHN were tested in a second population. One AFLP marker was associated with decreased symptoms in both populations. The SSR marker was not associated with IHN in the second population, but was closely linked in repulsion to another marker that was associated with IHN, and had the same (negative) effect on the trait as the SSR marker did in the first population. The correlation between maturity and IHN may be partially explained by the presence of markers on chromosome V that are linked to both traits. This research represents the first molecular genetic research of IHN in potato.

Mapping of QTLs governing agronomic and yield traits in chickpea

Chickpea is one of the most important leguminous cool season food crops, cultivated prevalently in South Asia and Middle East. The main objective of this study was to identify quantitative trait loci (QTLs) associated with seven agronomic and yield traits in two recombinant inbred line populations of chickpea derived from the crosses JG62 × Vijay (JV population) and Vijay × ICC4958 (VI population) from at least three environments. Single locus QTL analysis involved composite interval mapping (CIM) for individual traits and multiple-trait composite interval mapping (MCIM) for correlated traits to detect pleiotropic QTLs. Two-locus analysis was conducted to identify the main effect QTLs (M-QTLs), epistatic QTLs (E-QTLs) and QTL × environment interactions. Through CIM analysis, a total of 106 significant QTLs (41 in JV and 65 in VI populations) were identified for the seven traits, of which one QTL each for plant height and days to maturity was common in both the populations. Six pleiotropic QTLs that were consistent over the environments were also identified. LG2 in JV and LG1a in VI contained at least one QTL for each trait. Hence, concentrating on these LGs in molecular breeding programs is most likely to bring simultaneous improvement in these traits.

A BAC/BIBAC-based physical map of chickpea, Cicer arietinum L

Background

Chickpea (Cicer arietinum L.) is the third most important pulse crop worldwide. Despite its importance, relatively little is known about its genome. The availability of a genome-wide physical map allows rapid fine mapping of QTL, development of high-density genome maps, and sequencing of the entire genome. However, no such a physical map has been developed in chickpea.

Results

We present a genome-wide, BAC/BIBAC-based physical map of chickpea developed by fingerprint analysis. Four chickpea BAC and BIBAC libraries, two of which were constructed in this study, were used. A total of 67,584 clones were fingerprinted, and 64,211 (~11.7 ×) of the fingerprints validated and used in the physical map assembly. The physical map consists of 1,945 BAC/BIBAC contigs, with each containing an average of 28.3 clones and having an average physical length of 559 kb. The contigs collectively span approximately 1,088 Mb. By using the physical map, we identified the BAC/BIBAC contigs containing or closely linked to QTL4.1 for resistance to Didymella rabiei (RDR) and QTL8 for days to first flower (DTF), thus further verifying the physical map and confirming its utility in fine mapping and cloning of QTL.

Conclusion

The physical map represents the first genome-wide, BAC/BIBAC-based physical map of chickpea. This map, along with other genomic resources previously developed in the species and the genome sequences of related species (soybean, Medicago and Lotus), will provide a foundation necessary for many areas of advanced genomics research in chickpea and other legume species. The inclusion of transformation-ready BIBACs in the map greatly facilitates its utility in functional analysis of the legume genomes.

Genetic analyses and conservation of QTL for ascochyta blight resistance in chickpea (Cicer arietinum L.)

Ascochyta blight (AB) caused by Ascochyta rabiei (teleomorph, Didymella rabiei) Pass. Lab. is an important fungal disease of chickpea worldwide. Only moderate sources of resistance are available within the cultivated species and we hypothesized that the available sources may carry different genes for resistance, which could be pyramided to improve field resistance to AB. Four divergent moderately resistant cultivars CDC Frontier, CDC Luna, CDC Corinne, and Amit were each crossed to a highly susceptible germplasm ICCV 96029. Parents, F(1) and F(2) generations were evaluated under controlled conditions for their reactions to AB. A total of 144 simple sequence repeat (SSR) markers were first mapped to eight linkage groups (LG) for the CDC Frontier x ICCV 96029 population. Then based on the evidence from this population, 76, 61, and 42 SSR markers were systematically chosen and mapped in CDC Luna, CDC Corinne, and Amit populations, respectively. Frequency distributions of the AB rating in the F(2) generation varied among the four populations. Composite interval mapping revealed five QTLs (QTL1-5), one on each of LG 2, 3, 4, 6, and 8, respectively, distributed across different sources, controlling resistance to AB. CDC Frontier contained QTL2, 3, and 4 that simultaneously accounted for 56% of phenotypic variations. CDC Luna contained QTL 1 and 3. CDC Corinne contained QTL 3 and 5, while only QTL 2 was identified in Amit. Altogether these QTL explained 48, 38, and 14% of the estimated phenotypic variations in CDC Luna, CDC Corinne, and Amit populations, respectively. The results suggested that these QTLs could be combined into a single genotype to enhance field resistance to AB.

The marker SCK13(603) associated with resistance to ascochyta blight in chickpea is located in a region of a putative retrotransposon

The sequence characterized amplified region (SCAR) marker SCK13(603), associated with ascochyta blight resistance in a chickpea recombinant inbred line (RIL) population, was used as anchored sequence for genome walking. The PCRs performed in the walking steps to walk in the same direction produced eight bands in 5' direction and five bands in 3' direction with a length ranking from 530 to 2,871 bp. The assembly of the bands sequences along with the sequence of SCK13(603) resulted in 7,815 bp contig. Blastn analyses showed stretches of DNA sequence mainly distributed from the nucleotides 1,500 to 4,500 significantly similar to Medicago truncatula genomic DNA. Three open reading frames (ORFs) were identified and blastp analysis of predicted amino acids sequences revealed that ORF1, ORF2 and ORF3 had significant similarity to a CCHC zinc finger protein, to an integrase, and to a precursor of the glucoamylase s1/s2, respectively, from M. truncatula. The high homology of the putative proteins derived from ORF1 and ORF2 with retrotransposon proteins and the prediction of the existence of conserved domains usually present in retrotransposon proteins indicate that the marker SCK13(603) is located in a region of a putative retrotransposon. The information generated in this study has contributed to increase the knowledge of this important region for blight resistance in chickpea.

Genetic mapping of ascochyta blight resistance in chickpea (Cicer arietinum L.) using a simple sequence repeat linkage map

Ascochyta blight, caused by the fungus Ascochyta rabiei (Pass.) Lab., is one of the most devastating diseases of chickpea (Cicer arietinum L.) worldwide. Research was conducted to map genetic factors for resistance to ascochyta blight using a linkage map constructed with 144 simple sequence repeat markers and 1 morphological marker (fc, flower colour). Stem cutting was used to vegetatively propagate 186 F2 plants derived from a cross between Cicer arietinum L. 'ICCV96029' and 'CDC Frontier'. A total of 556 cutting-derived plants were evaluated for their reaction to ascochyta blight under controlled conditions. Disease reaction of the F1 and F2 plants demonstrated that the resistance was dominantly inherited. A Fain's test based on the means and variances of the ascochyta blight reaction of the F3 families showed that a few genes were segregating in the population. Composite interval mapping identified 3 genomic regions that were associated with the reaction to ascochyta blight. One quantitative trait locus (QTL) on each of LG3, LG4, and LG6 accounted for 13%, 29%, and 12%, respectively, of the total estimated phenotypic variation for the reaction to ascochyta blight. Together, these loci controlled 56% of the total estimated phenotypic variation. The QTL on LG4 and LG6 were in common with the previously reported QTL for ascochyta blight resistance, whereas the QTL on LG3 was unique to the current population.

Development of an integrated intraspecific map of chickpea (Cicer arietinum L.) using two recombinant inbred line populations

A composite intraspecific linkage map of chickpea was developed by integrating individual maps developed from two F(8:9) RIL populations with one common parent. Different molecular markers viz. RAPD, ISSR, RGA, SSR and ASAP were analyzed along with three yield related traits: double podding, seeds per pod and seed weight. A total of 273 markers and 186 RILs were used to generate the map with eight linkage groups at a LOD score of >/=3.0 and maximum recombination fraction of 0.4. The map spanned 739.6 cM with 230 markers at an average distance of 3.2 cM between markers. The predominantly used SSR markers facilitated identification of homologous linkage groups from the previously published interspecific linkage map of chickpea and confirmed conservation of the SSR markers across the two maps as well as the variation in terms of marker distance and order. The double podding gene was tagged by the markers NCPGR33 and UBC249z at 2.0 and 1.1 cM, respectively. Whereas, seeds per pod, was tagged by the markers TA2x and UBC465 at 0.1 and 1.8 cM, respectively. Eight QTLs were identified that influence seed weight. The joint map approach allowed mapping a large number of markers with a moderate coverage of the chickpea genome and few linkage gaps.

Expression profiling of chickpea genes differentially regulated during a resistance response to Ascochyta rabiei

Using microarray technology and a set of chickpea (Cicer arietinum L.) unigenes, grasspea (Lathyrus sativus L.) expressed sequence tags (ESTs) and lentil (Lens culinaris Med.) resistance gene analogues, the ascochyta blight (Ascochyta rabiei (Pass.) L.) resistance response was studied in four chickpea genotypes, including resistant, moderately resistant, susceptible and wild relative (Cicer echinospermum L.) genotypes. The experimental system minimized environmental effects and was conducted in reference design, in which samples from mock-inoculated controls acted as reference against post-inoculation samples. Robust data quality was achieved through the use of three biological replicates (including a dye swap), the inclusion of negative controls and strict selection criteria for differentially expressed genes, including a fold change cut-off determined by self-self hybridizations, Student's t-test and multiple testing correction (P < 0.05). Microarray observations were also validated by quantitative reverse transcriptase-polymerase chain reaction (RT-PCR). The time course expression patterns of 756 microarray features resulted in the differential expression of 97 genes in at least one genotype at one time point. k-means clustering grouped the genes into clusters of similar observations for each genotype, and comparisons between A. rabiei-resistant and A. rabiei-susceptible genotypes revealed potential gene 'signatures' predictive of effective A. rabiei resistance. These genes included several pathogenesis-related proteins, SNAKIN2 antimicrobial peptide, proline-rich protein, disease resistance response protein DRRG49-C, environmental stress-inducible protein, leucine-zipper protein, polymorphic antigen membrane protein, Ca-binding protein and several unknown proteins. The potential involvement of these genes and their pathways of induction are discussed. This study represents the first large-scale gene expression profiling in chickpea, and future work will focus on the functional validation of the genes of interest.

Mapping quantitative trait loci in chickpea associated with time to flowering and resistance to Didymella rabiei the causal agent of Ascochyta blight

Drought is the major constraint to chickpea (Cicer arietinum L.) productivity worldwide. Utilizing early-flowering genotypes and advancing sowing from spring to autumn have been suggested as strategies for drought avoidance. However, Ascochyta blight (causal agent: Didymella rabiei (Kov.) v. Arx.) is a major limitation for chickpea winter cultivation. Most efforts to introgress resistance to the pathogen into Kabuli germplasm resulted in relatively late flowering germplasm. With the aim to explore the feasibility of combining earliness and resistance, RILs derived from a cross between a Kabuli cultivar and a Desi accession were evaluated under field conditions and genotyped with SSR markers. Three quantitative trait loci (QTLs) with significant effects on resistance were identified: two linked loci located on LG4 in epistatic interaction and a third locus on LG8. Two QTLs were detected for time to flowering: one in LG1 and another on LG2. When resistance and time to flowering were analyzed together, the significance of the resistance estimates obtained for the LG8 locus increased and the locus effect on days to flowering, previously undetected, was significantly different from zero. The identification of a locus linked both to resistance and time to flowering may account for the correlation observed between these traits in this and other breeding attempts.

Detection of two quantitative trait loci for resistance to ascochyta blight in an intra-specific cross of chickpea (Cicer arietinum L.): development of SCAR markers associated with resistance

Two quantitative trait loci (QTLs), (QTL(AR1) and QTL(AR2)) associated with resistance to ascochyta blight, caused by Ascochyta rabiei, have been identified in a recombinant inbred line population derived from a cross of kabulixdesi chickpea. The population was evaluated in two cropping seasons under field conditions and the QTLs were found to be located in two different linkage groups (LG4a and LG4b). LG4b was saturated with RAPD markers and four of them associated with resistance were sequenced to give sequence characterized amplified regions (SCARs) that segregated with QTL(AR2). This QTL explained 21% of the total phenotypic variation. However, QTL(AR1), located in LG4a, explained around 34% of the total phenotypic variation in reaction to ascochyta blight when scored in the second cropping season. This LG4a region only includes a few markers, the flower colour locus (B/b), STMS GAA47, a RAPD marker and an inter-simple-sequence-repeat and corresponds with a previously reported QTL. From the four SCARs tagging QTL(AR2), SCAR (SCY17(590)) was co-dominant, and the other three were dominant. All SCARs segregated in a 1:1 (presence:absence) ratio and the scoring co-segregated with their respective RAPD markers. QTL(AR2) on LG4b was mapped in a highly saturated genomic region covering a genetic distance of 0.8 cM with a cluster of nine markers (three SCARs, two sequence-tagged microsatellite sites (STMS) and four RAPDs). Two of the four SCARs showed significant alignment with genes or proteins related to disease resistance in other species and one of them (SCK13(603)) was sited in the highly saturated region linked to QTL(AR2). STMS TA72 and TA146 located in LG4b were described in previous maps where QTL for blight resistance were also localized in both inter and intraspecific crosses. These findings may improve the precision of molecular breeding for QTL(AR2) as they will allow the choice of as much polymorphism as possible in any population and could be the starting point for finding a candidate resistant gene for ascochyta blight resistance in chickpea.

A linkage map of chickpea (Cicer arietinum L.) based on populations from Kabuli x Desi crosses: location of genes for resistance to fusarium wilt race 0

Two recombinant inbred line (RIL) populations derived from intraspecific crosses with a common parental line (JG62) were employed to develop a chickpea genetic map. Molecular markers, flower colour, double podding, seed coat thickness and resistance to fusarium wilt race 0 (FOC-0) were included in the study. Joint segregation analysis involved a total of 160 markers and 159 RILs. Ten linkage groups (LGs) were obtained that included morphological markers and 134 molecular markers (3 ISSRs, 13 STMSs and 118 RAPDs). Flower colour (B/b) and seed coat thickness (Tt/tt) appeared to be linked to STMS (GAA47). The single-/double-podding locus was located on LG9 jointly with two RAPD markers and STMS TA80. LG3 included a gene for resistance to FOC-0 (Foc0(1)/foc0(1)) flanked by RAPD marker OPJ20(600) and STMS marker TR59. The association of this LG with FOC-0 resistance was confirmed by QTL analysis in the CA2139 x JG62 RIL population where two genes were involved in the resistance reaction. The STMS markers enabled comparison of LGs with preceding maps.