Next Article in Journal
Silicon Application Induced Alleviation of Aluminum Toxicity in Xaraés Palisadegrass
Previous Article in Journal
Nitrous Oxide Emission from Forage Plantain and Perennial Ryegrass Swards Is Affected by Belowground Resource Allocation Dynamics
Previous Article in Special Issue
Potential Application of Genomic Technologies in Breeding for Fungal and Oomycete Disease Resistance in Pea
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of Genomics Approaches for the Improvement in Ascochyta Blight Resistance in Chickpea

1
AgriBio, Centre for AgriBioscience, Agriculture Victoria, Bundoora, VIC 3083, Australia
2
School of Applied Systems Biology, La Trobe University, Bundoora, VIC 3086, Australia
3
NSW Department of Primary Industries, 4 Marsden Park Rd, Calala, NSW 2340, Australia
*
Author to whom correspondence should be addressed.
Agronomy 2021, 11(10), 1937; https://doi.org/10.3390/agronomy11101937
Submission received: 14 July 2021 / Revised: 17 September 2021 / Accepted: 20 September 2021 / Published: 27 September 2021

Abstract

:
Advancements in high-throughput genotyping and sequencing technologies are enabling the development of a vast range of genomic tools and resources for a new revolution in plant breeding. Several genotyping-by-sequencing (GBS) methods including capture-based, genome complexity reduction and sequencing of cDNA (GBS-t) are available for application in trait dissection, association mapping, and genomic selection (GS) in crop plants. The aims of this study were to identify genomic regions conferring resistance to Ascochyta blight (AB) introgressed from the wild Cicer echinospernum into the domesticated C. arietinum, through a conventional recombinant inbred population genotyped using a variety of GBS methods. Evaluation of GBS methods revealed that capture-based approaches are robust and reproducible while GBS-t is rapid and flexible. A genetic linkage map consisting of 5886 polymorphic loci spanning 717.26 cM was generated. Using field phenotyping data from two years, a single genomic region on LG4 was identified with quantitative trait loci (QTL) mapping. Both GBS methods reported in this study are well suited for applications in genomics assisted plant breeding. Linked markers for AB resistance, identified in the current study, provide an important resource for the deployment into chickpea breeding programs for marker-assisted selection (MAS).

1. Introduction

Chickpea (Cicer arietinum L.) is the third most important grain legume crop after soybean and pea, grown worldwide on c. 17.8 Mha with an annual production of 17.2 million tonnes [1]. Due to its high protein content, chickpea serves as an important protein source in the vegetarian diet of many developing nations. Under optimal growth conditions, chickpea productivity can reach up to 6 t/ha; however, exposure to biotic and abiotic stresses can reduce the average yield to <1 t/ha [2]. Among these stresses, Ascochyta blight (AB), caused by the fungus Ascochyta rabiei (Pass.) Labr. [teleomorph: Didymella rabiei (Kovacheski) von Arx (= Mycosphaerella rabiei Kovacheski)], is one of the major diseases that can cause a yield loss of 10–100%. Different strategies can be applied for the management of AB, such as short-acting seed treatments, strategic foliar fungicide application, and growing resistant varieties [3]. Development of resistant varieties is the most effective and viable option for ongoing sustainability and has been one of the major aims for chickpea breeding programs across the world. Conventional breeding approaches have been successful, to some extent, in improving the levels of AB resistance in breeder’s germplasm especially for cultivated chickpea. Previous studies have reported several major resistance genes along with multiple minor genes controlling AB [4,5,6,7]. However, major resistance sources being incorporated and widely used in breeding programs are always at a risk of boom-and-bust cycles due to constant evolution of pathogen structure [6]. Introgression of minor resistance sources along with major sources through genomics assisted breeding is anticipated to provide durable AB resistance. To maintain AB resistance in chickpea varieties, it is necessary to explore new sources of resistance not only within the cultivated gene pool but also from wild gene pools. Wild chickpea relatives adapted to different conditions and geographical locations offer tremendous genetic diversity [8]. Previous studies have identified new sources of AB resistance in C. bijugum K. H. Rech., C. echinospermum P. H. Davis, C. pinnatifidum Jaub. & Sp., and C. reticulatum Ladiz. [6,9]. Due to the feasibility of introgression and generation of fertile hybrids, C. echinospermum and C. reticulatum are known to be the common wild relatives of chickpea used in germplasm enhancement through conventional breeding [10].
Advances in next-generation sequencing technologies have provided opportunities to enrich genomic and genetic resources for chickpea. Currently, large-scale genomic resources are focused on a collection of reference genomes (kabuli type [11], desi type [12], wild chickpea [13]) as well as collection of SNP markers [3,14] and genetic maps [2]. Discovery of genome wide SNPs using genotyping-by-sequencing (GBS) approaches is a vital tool for marker trait association studies leading to genetic enhancement. High-throughput SNP genotyping approaches include complexity reduction-based methods [15,16,17], skim genome sequencing, and target capture assays [18]. Recently, high-throughput SNP genotyping methods, including restriction site-associated DNA sequencing (RAD-Seq) as well as skim genome sequencing, were used to genotype chickpea mapping populations [19,20]. In this study, a combination of approaches including capture-based and GBS-transcriptomics (GBS-t) sequencing were utilised to genotype an interspecific RIL mapping population and identify regions associated with AB resistance. Target capture-based systems are more reproducible and specific and they provide enrichment of targeted regions [21]; however, they require prior knowledge of the genome sequence and structure in order to select specific regions for amplification. GBS-t is a type of target capture assay that also serves as complexity reduction method by targeting the whole exome of any given species. This strategy is scalable, economical, and allows for deeper sequencing coverage compared to whole genome approaches. It is also highly suitable for species with limited genomic resources available and has been exemplified across a broad range of species and ploidy levels [22,23,24].
Molecular breeding offers many advantages over traditional breeding in terms of time and cost savings as well as breeding efficiency and selection accuracy. Genetic linkage map construction for use in trait dissection and marker-assisted breeding has been a part of chickpea breeding programs. As a result, a large number of maps based on different types of markers such as SSRs and SNPs have been generated [3,25,26]. Both intraspecific as well as interspecific maps were published; however, more emphasis has been given to the interspecific maps, particularly those between C. arietinum and C. reticulatum [27]. The rationale for the use of interspecific mapping populations has been the low levels of genetic diversity that are known to be present within the gene pool of cultivated chickpea [28]. Multiple quantitative trait loci (QTLs) for resistance to ascochyta blight have been identified, from both inter and intra-specific crosses, for resistance at different growth stages (including seedling, flowering, and pod filling) under greenhouse and field conditions [3,5,7,20]. These QTL analyses have historically accelerated the dissection of genetic control of AB, thus allowing marker-assisted selection (MAS) to be implemented in breeding processes [3,29]. A limitation of MAS is that it focuses on the selection of major resistance genes that can be dramatically overcome by rapid adaptation of A. rabiei pathogen. Moreover, there is limited research on genes underlaying these QTLs and the resistance mechanism against AB infection.
As genetic diversity in chickpea suffers from a domestication bottleneck, it is imperative to utilize wild relatives to introgress agronomically important traits and expand the gene pool [30]. Introgression from wild relatives will likely deliver a single genomic region that would provide enhanced major gene or gene set resistance. Major gene resistance, when used in combination with other sources to improve the quantitative resistance level within the germplasm, is likely to be the most effective mean of control of AB pathogens. In this study, we used a recombinant inbred mapping population derived from a cross between C. arietinum (genotype Sonali) and C. echinospermum (genotype 04067-81-2-1-1) that was generated to introgress AB resistance into commonly cultivated genetics. This study evaluated capture-based and GBS-t sequencing approaches utilising Illumina sequencing. SNPs identified using both approaches were used for genetic linkage map construction and QTL analysis for identification of genomic region for AB resistance.

2. Materials and Methods

2.1. Plant Materials and Nucleic Acid Extraction

Seeds from a total of seven chickpea varieties/advanced breeding lines (Genesis114, Howzat, ICC3996, Lasseter, Sonali, Yorker, 04067-81-2-1-1) were obtained from NSW-DPI, Tamworth, Australia. These seven lines were extensively used in chickpea breeding programs and were the parents of RIL mapping populations segregating for key traits of interest. Multiple seeds of each genotype were germinated and maintained for 2–3 weeks in glasshouse at 22  ±  2 °C under a 16/8 h (light/dark) photoperiod, in individual pots filled with potting mix (AGH Mix, Bio Gro, South Australia, Australia) at the premises of Agriculture Victoria, Bundoora, Victoria, and Australia.
An interspecific RIL mapping population was developed by crossing a susceptible cultivar, Sonali (C. arietinum), and a highly resistant line, 04067-81-2-1-1 (C. echinospermum). The RIL mapping population consisted of 134 individuals based on single seed descent from F2 to F6 progeny in the glasshouse (NSW-DPI, Tamworth, Australia).
For nucleic acid extractions, 3–4 young leaflets were harvested, snap frozen in liquid nitrogen and stored at −80 °C until required. Genomic DNA extraction was performed using DNeasy® 96 Plant Kit (QIAGEN, Hilden, Germany) and RNA extractions using RNeasy® 96 Kit (QIAGEN) following manufacturer’s instructions, with the exception that Dithiothreitol (DTT) (Sigma-Aldrich, St. Louis, MO, USA) was used instead of β-mercaptoethanol. Tissue disruption was performed in a 96-well plate (QIAGEN Microtubes) using a Mixer Mill 300 (Retsch®, Haan, Germany) and was eluted in 40 µL of elution buffer and stored at −80 °C. Concentration and quality being confirmed by a Thermo Scientific™ NanoDropTM UV-Visible spectrophotometer (Thermo Fisher Scientific Inc., Waltham, MA, USA) and integrity evaluated using a TapeStation 2200 platform, using either gDNA ScreenTape or RNA ScreenTape (Agilent Technologies, Santa Clara, CA, USA).

2.2. GBS-t Sequencing Library Preparation, SNP Identification and Target Enrichment Array Design

Following RNA extraction from above mentioned seven chickpea varieties/advanced breeding lines, paired-end sequencing libraries (SureSelect Strand Specific RNA Library Preparation Kit [Agilent Technologies]) were constructed with insert sizes of 350 bp. The resulting libraries were evaluated using the TapeStation 2200 platform (Agilent Technologies). Each library was generated with a unique barcode and pooled in equal mass. Pooled libraries were quantified using the KAPA library quantification kit (KAPA Biosystems, Wilmington, MA, USA). Sequencing was performed using an Illumina HiSeq 2000 Sequencing platform producing 2× 150 bp reads, and each sample was sequenced to a depth of 5–6 million paired reads per sample (Illumina Inc., San Diego, CA, USA).
Following fastq data generation, the raw sequence reads were filtered using a custom perl script to remove adaptor sequences along with reads and bases of low quality (Q ≤ 30). Reads with 3 consecutive unassigned nucleotides (N) were also trimmed and, finally, any reads shorter than 50 bp in length were removed from the final set. The remaining high-quality trimmed sequence reads were aligned to the chickpea reference genome sequence from cultivar CDC Frontier [11] using TopHat2 (version 2.1.1) [31]. Variant calling was performed using SAMtools (version-1.5) [32] and BCFtools [33]. The final VCF output was then filtered based on various parameters: depth (DP → 10), minor allelic frequency (MAF 0.05), maximum missing data (20%), and base quality (Q30) using vcftools [33].
For target enrichment array design, a final set of SNPs was selected based on uniform distribution across the genome and flanking sequences 100 bp both upstream and downstream from each target SNP were sent to Agilent Technologies for target capture assay design (Table S1).

2.3. Genotyping the RIL Mapping Population

GBS-t sequencing libraries (SureSelect Strand Specific RNA Library Preparation Kit [Agilent Technologies]) were used, and SNP identification of RIL mapping population was performed, as described above, with slight modifications on the number of paired-end reads generated (around 2–3 million paired-end reads for each RILs).
Target enrichment was performed according to the manufacturer’s instructions (SureSelect Target Enrichment [Agilent Technologies]), except for a reduction in reaction volume (1/4 of the reaction). High-quality DNA (200 ng) was randomly fragmented using MspJI restriction enzyme (New England Biolabs, Ipswich, MA, USA), as described previously [34]. All reads were pair-end sequenced using the MiSeq platform (Illumina Inc.). The generated sequence data were trimmed as before and then aligned to the chickpea reference genome sequence, i.e., the CDC Frontier [11] using BWA MEM [35]. SNP calling was performed using SAMtools (version-1.5) [32] and BCFtools [33] with a defined SNP list from target enrichment assay. VCF output was processed with the same parameters used for GBS-t sequencing data filtering using vcftools [33].
Additionally, 12 markers (2 SSRs and 10 SNP markers) flanking the AB resistance QTLs from a previous study [3] were screened on Sonali X 04067-81-2-1-1 RIL population. SSR genotyping was performed, as described previously [3], as well as SNP genotyping using KASPTM genotyping chemistry (LGC Genomics, Middlesex, UK) with some modifications [36].

2.4. Genetic Linkage Mapping

SNP genotyping data showing contrasting alleles between the parental genotypes were selected for linkage mapping. A chi-squared test was performed on the filtered SNPs and any markers showing significant deviation from 1:1 ratio (p < 0.01) were excluded from further analysis. Linkage maps were constructed from the final set of SNPs using the statistical software QTL IciMapping V4.1 [37] with a logarithm of odds (LOD) score of 14. To finalize the linkage map, marker ordering and rippling were performed using RECORD and COUNT algorithms, respectively. The Mapchart 2.2 software [38] was used to visualise the final genetic map. Visual comparisons of the genetic linkage maps with the physical map was performed using the Strudel software package [39].

2.5. AB Inoculation, Disease Rating and Statistical Analysis of Phenotyping Data

The phenotypic evaluation of AB resistance in Sonali X 04067-81-2-1-1 RIL population was performed, as described previously [3]. Briefly, the RIL population and parents were sown in a 0.6 m single row plot (25 seeds/plot) in field in Tamworth, NSW, Australia in 2014 and 2015 in a randomized complete block design with three replicates. Once plants were established, inoculation was performed with the spore suspension (500,000 conidia/mL at a rate of 100 L per hectare). Disease symptoms were scored using a scale of 1–9 based on whole plot severity.
The means of the phenotypic scores were estimated after the adjustment for any spatial patterning within the trial. Models were fitted using residual maximum likelihood (REML), as implemented in GenStat (GenStat Committee, 2002) [40]. Phenotypic scores from both trials (2014 and 2015) were plotted and the score distributions were compared in R [41]. Best Linear Unbiased Estimates (BLUEs) analysis was used to calculate the broad sense heritability (H2) using ASReml v4.1 [42].

2.6. Identification of Genomic Regions for AB Resistance

To identify genomic regions associated with AB resistance, QTL mapping was performed. Both simple interval mapping (SIM) and composite interval mapping (CIM) were performed in Windows QTL Cartographer version 2.5 [43]. An arbitrary LOD threshold of 2.5 was used to determine significance for SIM. For CIM, significance levels for LOD thresholds were determined using 1000 permutations. Based on the sequences of SNP loci flanking QTL intervals, genomic regions containing putative candidate genes were retrieved (based on annotation).

3. Results

3.1. Target Enrichment Array Design from Advanced Varieties/Advanced Breeding Lines

A total of 789,525,381 high-quality trimmed paired end transcriptome reads obtained from seven chickpea varieties/advanced breeding lines were aligned to the reference genome of chickpea (v1.0). Table S2 provides details of paired end reads that were mapped to the reference from each genotype. A total of 368,449 base variants were identified that were filtered based on depth, quality score, missing data, and minor allele frequency to generate a final set of 26,563 SNPs. Finally, 11,222 SNPs were selected based on uniform distribution across the chickpea genome and used for target capture design. Details of these selected SNPs, genomic location, and probe sequences are provided in Table S1.

3.2. Genotyping the RIL Mapping Population and Genetic Linkage Map Construction

To identify SNPs in the RIL mapping population, capture-based and GBS-t sequencing approaches were performed. In the capture-based approach, sequencing generated a total of 48,503,448 paired end reads with an average of 350,000 reads per line. Most of the probes (10,732; 95%) were identified as being on-target sequences with reads distributed evenly across all targeted regions. Among them, 55% of the probes (5854) were polymorphic between the two RIL parents. A final set of 2479 SNP markers were selected for linkage mapping after additional filtering based on missing data (<20%), heterozygosity (<10%), and χ2 analysis (p < 0.05) (Table S3).
From GBS-t sequencing of the RIL mapping population, a total of 329 million high-quality reads were used for analysis. After the initial filtering based on read quality, read depth and missing data, 3856 polymorphic SNP markers were identified. After performing a χ2 analysis (p < 0.05), a final set containing 3591 segregating SNP markers was used for linkage mapping (Table S3). From previously published AB resistance QTL flanking markers, two SSR and three SNP markers were also included in linkage mapping (Table S3).
Polymorphic marker data from all these genotyping approaches were combined (6075); only 5919 unique markers (based on the genomic position) were used for linkage map construction. Mapping resulted in final set of 5886 markers that were assigned to eight major linkage groups (that equaled to the eight chromosome) and three satellites (Figure 1, Table 1 and Table S4). The proportion of loci assigned to LGs was 99%, with a small proportion marker unlinked. Total map length was calculated to be 717 cM with an average inter-marker distance of 0.12 cM. The marker distribution across the LG’s were not uniform; LG4 had the highest number of markers (1934) and LG5 with the lowest (105). Marker order on LG’s were compared to the chickpea genome, which identified significant commonality (Figure S1). However, inconsistencies were observed for a handful of markers which were assigned to different LG’s and some marker order was not co-linear (especially on LG7).

3.3. Phenotypic Analysis of RIL Populations and the Identification of Genomic Regions for AB Resistance

The AB disease score showed greater variation in the 2014 trial (range 1.8–5.2, mean 2.74 ± 0.79) than in 2015 (range 3–5, mean 3.68 ± 0.57), depicted in histograms provided in Figure S2. Correlation (Pearson) between the mean AB disease score from the two trials was 0.52 (Figure S3). The estimated heritability of the combined data set was 0.53.
SIM and CIM analyses detected a single genomic region on LG4 associated with resistance to AB (Figure 1). The amount of phenotypic variance (Vp) explained was 34% and 41% from 2014 and 2015 phenotyping trials, respectively (Table 2). Resistance determinant was derived from the C. echinospermum parent (04067-81-2-1-1). AB_echino QTL interval was found to be located on the Ca4 between 14–16 Mbp (Figure S3). From the physical interval, potential candidate genes were identified based on annotation, and they include serine/threonine protein phosphatase, receptor protein kinase-like protein, LRR receptor-like serine/threonine-protein kinase, and cysteine proteinase, which have been demonstrated to be involved in disease resistance (Table S5).

4. Discussion

High-throughput and low-cost sequencing technologies have extensive applications in marker-assisted breeding through linkage mapping, QTL analysis, and MAS [44]. Different GBS approaches vary greatly in the types/amount of data generated, their costs per marker data point, and errors/potential biases. Hence, selecting the best suited method depends on the genome size of the crop, population structure, availability of prior genomic information, reproducibility, ease of application, and the relative scale of the industry. In this study, two high-throughput approaches i.e., capture-based and GBS-t sequencing, were evaluated on a RIL mapping population. Genotyping by capture-based approach had a high proportion (95%) of the targeted SNPs recovered, with ~50% (5854 markers) being polymorphic in the RIL population, whilst the GBS-t approach identified 3856 polymorphic SNP markers, 3591 of which were used for linkage mapping. These results indicate that capture-based approach is more robust, reproducible, and scalable as compared to GBS-t. However, the main caveat/drawback of this method is the lack of flexibility of SNP panels, especially to apply on different germplasm collections such as exotic lines and wild relatives. To address this limitation, custom SNP panels need to be designed and integrated into existing assays, which will incur inevitable additional costs. The alternative is to apply a more flexible approach, such as GBS-t. Even though GBS-t is more flexible as compared to capture-based approach and can deliver functionally associated gene-based markers for the identifcaiton of expression QTLs [45], it still has its own disadvantages including the requirement of high-quality plant material for extraction of RNA and comparatively higher proportion of missing data compared to capture-based approaches. Due to relatively small genome and simple ploidy in chickpea, skim genome sequencing is another suitable approach for detection of genome wide polymorphisms. As the cost of sequencing is declining rapidly with the advent of new experimental protocols, equipment, and applications of sequencing, the limitation of missing data can be alleviated by higher coverage/depth [46]. In brief, two high-throughput genotyping approaches performed in this study have their own pros and cons. GBS-t is an effective genotyping approach for species with limited genomic resources to implement high-throughput genotyping quickly or to evaluate crosses with wild relatives. The development of capture-based assays will then be an obvious extension that incorporates the necessary variants for application in a breeding scenario or germplasm enhancement. Generic costs for both approaches are similar with obvious cost savings for operating at scale; this can potentially deliver greater savings to capture assays. With the availability of efficient computational methods, low-coverage sequencing data can be efficiently imputed to obtain dense marker coverage. There are several computational packages available for imputation, such as Minimac [47], Beagle [48], and LinkImpute [49]. Performance of these packages depends on several factors including the size and structure of population, effective population size, SNP density and sequencing coverage [50].
In this study, SNP markers from GBS data have successfully been used to build a high-density genetic map (5886 markers). The majority of the markers displayed conserved assignments with physical positions on the genome. Uneven marker distribution on certain regions of LG’s was observed, which could be due to lack of polymorphism in genes (suppressed for meiotic crossovers hence reduced recombination rates), simple repeats, retroelements, as well as gene-poor or heterochromatic regions interspersed in the genome [51]. The high-density Sonali x 04067-81-2-1-1 linkage map provided the basis for identification of a single genomic region on LG4 conferring AB resistance. The slight differences observed in the Vp proportions (34% and 41%) in 2014 and 2015 could be due to variability between field conditions of the two screening experiments. Currently, in Australia, chickpea cultivars only possess partial resistance to A. rabiei and this level of resistance can easily breakdown as the pathogen is highly variable. Hence, there is an urgent need to introduce novel sources of disease resistance from global germplasm collections, including landraces and wild relatives, into chickpea breeding programs to develop superior AB resistant varieties. A previous study has successfully incorporated high levels of AB and botrytis grey mould (BGM) resistance in cultivated chickpea from wild relatives including C. judicium and C. pinnatifidum. Populations developed from these crosses were later evaluated, and a proportion of derivatives were resistant to both AB and BGM [52]. Multiple studies have previously been conducted to identify superior sources of AB resistance, identifying both single gene (either dominant or recessive) and multiple genes controlling AB resistance [3,5,20,53]. A recent study which performed the physical mapping of previously reported AB resistance QTLs using flanking SSR and SNP marker sequences the placed majority of QTLs in clusters on Ca2 and Ca4, and only few QTLs on Ca3, Ca5, and Ca6. Previously identified QTLs on Ca4, qAB4.1, and qAB4.2 were mapped at 6–8 Mbp and 21–26 Mbp, respectively [20]. The genomic region identified from this study revealed a physical position at 14–16 Mbp. The fact that the resistance source comes from C. echinospermum (genotype 04067-81-2-1-1) suggests it is highly likely to be a novel allele (or set of alleles) that can add value and diversity to the overall resistance breeding strategy. An efficient combination of genotyping approaches along with high-throughput phenotyping can help to determine the precise location of the causal gene responsible for resistant/tolerant phenotypes. GBS approaches for mapping/fine-mapping and gene analysis has been exploited in different species, including pigeonpea for Fusarium wilt and sterility mosaic disease resistance [54], in groundnut for bacterial wilt [55], and in cotton for drought tolerance [56], allowing the identification of associated pathways and candidate genes/SNPs useful for genomics-assisted breeding. Examining the QTL interval from this study helped to identify potential candidate genes, serine/threonine-protein phosphatase, receptor protein kinase-like protein, LRR receptor-like serine/threonine-protein kinase, and cysteine proteinase, all of which play key roles in disease resistance in plants [57]. These genes could be plausible candidates, although further work should be performed for validation. As the resistance source identified in this study originates from a wild relative, there could be additional novel genes present within the QTL region that couldn’t be identified using the cultivated genome reference. In the future, long-read sequencing platforms can be used for the comprehensive fine-mapping of this region thereby absolutely defining the gene contents within QTL interval.
QTL mapping has been the method of choice over decades to identify trait-linked markers that can be deployed in breeding for MAS. Genomic regions and markers identified in this study will be highly useful for immediate deployment into chickpea breeding using MAS. However, QTL mapping poses some limitations for traits controlled by multiple genes such as AB [58]. Previous studies on AB in chickpea have focused on the identification of QTLs of moderate to major effects that tend to be highly effective but vulnerable to breakdown with rapid changes in pathogen races. Furthermore, small population sizes as well as the usage of biparental mapping populations in individual studies had limited statistical power to detect quantitative genes conferring partial resistance, which is known to also be an important component for achieving durable resistance to AB. By collating pre-existing datasets and performing meta-analysis, it should be possible to achieve a population size that provides sufficient statistical power to identify novel quantitative genes conferring partial resistance to AB. Several studies have shown that GS outperforms MAS for complex trait in terms of gain per unit time and cost [59,60]. It also enables the simultaneous introgression of multiple loci from cultivated, as well as wild, backgrounds into elite breeding germplasm for durable resistance. Studies have been published comparing different GS prediction models, the effect of marker density, population structure, and population size on prediction accuracy [61,62]. The effects of prediction accuracy varied for each trait under investigation, especially as accuracy decreases with the increase in the trait complexity [2]. Regardless of population type or trait, the pre-requisite for GS is to develop a training population comprising diverse lines with complete phenotypic and genotypic data. Advances in high-throughput phenotyping and genotyping platforms will play a major role in the application of GS for rapid advancement of genetic gain in breeding programs. A foreground selection based on genotyping/diagnostic markers, followed by GS, could be the way forward for chickpea breeding programs.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/agronomy11101937/s1. Figure S1: Visual representation of the comparison of linkage groups of Sonali X 04067-81-2-1-1 RIL population map with chickpea chromosomes. Corresponding lines indicate the common region. Figure S2: Frequency distribution histograms generated from AB phenotyping using disease scores for 2014 and 2015 trials in Sonali X 04067-81-2-1-1 RIL population. Figure S3: Comparison of the AB disease scores for 2014 and 2015 trials in Sonali X 04067-81-2-1-1 RIL population. Figure S4: Comparison of AB_echino QTL interval from Sonali X 04067-81-2-1-1 RIL population with previously published QTLs on Ca4. Table S1: Name and sequence information of the SNP markers from the target enrichment assay. Table S2: Details on the reads aligned to the chickpea genome from advanced varieties/advanced breeding lines. Table S3: Details on the type and number of markers used for linkage mapping. Table S4. Linkage map statistics from Sonali X 04067-81-2-1-1 RIL population. Table S5: Details of plausible candidate genes identified in the QTL interval and their annotation.

Author Contributions

S.S. performed transcriptome sequencing, SNP genotyping, map construction, QTL analysis, candidate gene identification and contributed to drafting the manuscript. H.V.K. phenotypic data analysis and assisted in drafting the manuscript. K.H. developed the mapping population, performed the phenotyping experiments and contributed to data interpretation. S.B. performed DNA extraction and assisted in genotyping. N.D. assisted in the phenotyping experiments. N.O.I.C. co-conceptualised project and assisted in drafting the manuscript. S.K. co-conceptualised and coordinated the project and assisted in drafting the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by funding from the Victorian Department of Jobs, Precincts and Regions, Australia and the Grains Research and Development Council, Australia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

The authors would like to thank Ben G. Cocks for helpful critical comments on the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Food Agriculture Organization of the United Nations. FAOSTAT Statistical Database; FAO: Rome, Italy, 2020. [Google Scholar]
  2. Roorkiwal, M.; Rathore, A.; Das, R.R.; Singh, M.K.; Jain, A.; Srinivasan, S.; Gaur, P.; Chellapilla, B.; Tripathi, S.; Li, Y.; et al. Genome-Enabled Prediction Models for Yield Related Traits in Chickpea. Front. Plant Sci. 2016, 7, 1666. [Google Scholar] [CrossRef] [Green Version]
  3. Stephens, A.; Lombardi, M.; Cogan, N.O.I.; Forster, J.W.; Hobson, K.; Materne, M.; Kaur, S. Genetic marker discovery, intraspecific linkage map construction and quantitative trait locus analysis of ascochyta blight resistance in chickpea (Cicer arietinum L.). Mol. Breed. 2013, 33, 297–313. [Google Scholar] [CrossRef]
  4. Li, Y.; Ruperao, P.; Batley, J.; Edwards, D.; Davidson, J.; Hobson, K.; Sutton, T. Genome Analysis Identified Novel Candidate Genes for Ascochyta Blight Resistance in Chickpea Using Whole Genome Re-sequencing Data. Front. Plant Sci. 2017, 8, 359. [Google Scholar] [CrossRef] [Green Version]
  5. Dey, S.K.; Singh, G. Resistance to Ascochyta blight in chickpea—Genetic basis. Euphytica 1993, 68, 147–153. [Google Scholar] [CrossRef]
  6. Sharma, M.; Ghosh, R. An Update on Genetic Resistance of Chickpea to Ascochyta Blight. Agronomy 2016, 6, 18. [Google Scholar] [CrossRef] [Green Version]
  7. Udupa, S.M.; Baum, M. Genetic dissection of pathotype-specific resistance to ascochyta blight disease in chickpea (Cicer arietinum L.) using microsatellite markers. Theor. Appl. Genet. 2002, 106, 1196–1202. [Google Scholar] [CrossRef]
  8. Coyne, C.J.; Kumar, S.; Von Wettberg, E.J.; Marques, E.; Berger, J.D.; Redden, R.J.; Ellis, T.N.; Brus, J.; Zablatzká, L.; Smýkal, P. Potential and limits of exploitation of crop wild relatives for pea, lentil, and chickpea improvement. Legum. Sci. 2020, 2, e36. [Google Scholar] [CrossRef] [Green Version]
  9. Collard, B.C.Y.; Ades, P.K.; Pang, E.C.K.; Brouwer, J.B.; Taylor, P.W.J. Prospecting for sources of resistance to ascochyta blight in wild Cicer species. Australas. Plant Pathol. 2001, 30, 271–276. [Google Scholar] [CrossRef]
  10. Pundir, R.P.S.; Mengesha, M.H. Cross compatibility between chickpea and its wild relative, Cicer echinospermum Davis. Euphytica 1995, 83, 241–245. [Google Scholar] [CrossRef] [Green Version]
  11. Varshney, R.K.; Song, C.; Saxena, R.K.; Azam, S.; Yu, S.; Sharpe, A.G.; Cannon, S.; Baek, J.; Rosen, B.D.; Tar’An, B.; et al. Draft genome sequence of chickpea (Cicer arietinum) provides a resource for trait improvement. Nat. Biotechnol. 2013, 31, 240–246. [Google Scholar] [CrossRef] [Green Version]
  12. Jain, M.; Misra, G.; Patel, R.; Priya, P.; Jhanwar, S.; Khan, A.W.; Shah, N.; Singh, V.K.; Garg, R.; Jeena, G.; et al. A draft genome sequence of the pulse crop chickpea (Cicer arietinum L.). Plant J. 2013, 74, 715–729. [Google Scholar] [CrossRef]
  13. Gupta, S.; Nawaz, K.; Parween, S.; Roy, R.; Sahu, K.; Pole, A.K.; Khandal, H.; Srivastava, R.; Parida, S.K.; Chattopadhyay, D. Draft genome sequence of Cicer reticulatum L., the wild progenitor of chickpea provides a resource for agronomic trait improvement. DNA Res. 2016, 24, 1–10. [Google Scholar] [CrossRef] [Green Version]
  14. Roorkiwal, M.; Jain, A.; Kale, S.M.; Doddamani, D.; Chitikineni, A.; Thudi, M.; Varshney, R.K. Development and evaluation of high-density Axiom® CicerSNP Array for high-resolution genetic mapping and breeding applications in chickpea. Plant Biotechnol. J. 2017, 16, 890–901. [Google Scholar] [CrossRef] [Green Version]
  15. Elshire, R.; Glaubitz, J.C.; Sun, Q.; Poland, J.; Kawamoto, K.; Buckler, E.; Mitchell, S.E. A Robust, Simple Genotyping-by-Sequencing (GBS) Approach for High Diversity Species. PLoS ONE 2011, 6, e19379. [Google Scholar] [CrossRef] [Green Version]
  16. Baird, N.A.; Etter, P.D.; Atwood, T.S.; Currey, M.C.; Shiver, A.L.; Lewis, Z.; Selker, E.U.; Cresko, W.; Johnson, E.A. Rapid SNP Discovery and Genetic Mapping Using Sequenced RAD Markers. PLoS ONE 2008, 3, e3376. [Google Scholar] [CrossRef]
  17. Miller, M.R.; Dunham, J.P.; Amores, A.; Cresko, W.; Johnson, E.A. Rapid and cost-effective polymorphism identification and genotyping using restriction site associated DNA (RAD) markers. Genome Res. 2007, 17, 240–248. [Google Scholar] [CrossRef] [Green Version]
  18. Dapprich, J.; Ferriola, D.; Mackiewicz, K.; Clark, P.M.; Rappaport, E.; D’Arcy, M.; Sasson, A.; Gai, X.; Schug, J.; Kaestner, K.H.; et al. The next generation of target capture technologies—Large DNA fragment enrichment and sequencing determines regional genomic variation of high complexity. BMC Genom. 2016, 17, 486. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Kale, S.; Jaganathan, D.; Ruperao, P.; Chen, C.; Punna, R.; Kudapa, H.; Thudi, M.; Roorkiwal, M.; Katta, M.A.; Doddamani, D.; et al. Prioritization of candidate genes in “QTL-hotspot” region for drought tolerance in chickpea (Cicer arietinum L.). Sci. Rep. 2015, 5, 15296. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Deokar, A.; Sagi, M.; Tar’An, B. Genome-wide SNP discovery for development of high-density genetic map and QTL mapping of ascochyta blight resistance in chickpea (Cicer arietinum L.). Theor. Appl. Genet. 2019, 132, 1861–1872. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Giolai, M.; Paajanen, P.; Verweij, W.; Percival-Alwyn, L.; Baker, D.; Witek, K.; Jupe, F.; Bryan, G.; Hein, I.; Jones, J.D.G.; et al. Targeted capture and sequencing of gene-sized DNA molecules. BioTechniques 2016, 61, 315–322. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Malmberg, M.M.; Pembleton, L.W.; Baillie, R.C.; Drayton, M.C.; Sudheesh, S.; Kaur, S.; Shinozuka, H.; Verma, P.; Spangenberg, G.C.; Daetwyler, H.D.; et al. Genotyping-by-sequencing through transcriptomics: Implementation in a range of crop species with varying reproductive habits and ploidy levels. Plant Biotechnol. J. 2017, 16, 877–889. [Google Scholar] [CrossRef]
  23. Puritz, J.B.; Lotterhos, K.E. Expressed exome capture sequencing: A method for cost-effective exome sequencing for all organisms. Mol. Ecol. Resour 2018, 18, 1209–1222. [Google Scholar] [CrossRef] [Green Version]
  24. Warr, A.; Robert, C.; Hume, D.; Archibald, A.; Deeb, N.; Watson, M. Exome Sequencing: Current and Future Perspectives. G3: Genes Genomes Genet. 2015, 5, 1543–1550. [Google Scholar] [CrossRef] [Green Version]
  25. Gaur, R.; Azam, S.; Jeena, G.; Khan, A.W.; Choudhary, S.; Jain, M.; Yadav, G.; Tyagi, A.K.; Chattopadhyay, D.; Bhatia, S. High-Throughput SNP Discovery and Genotyping for Constructing a Saturated Linkage Map of Chickpea (Cicer arietinum L.). DNA Res. 2012, 19, 357–373. [Google Scholar] [CrossRef]
  26. Kujur, A.; Bajaj, D.; Upadhyaya, H.D.; Das, S.; Ranjan, R.; Shree, T.; Saxena, M.S.; Badoni, S.; Kumar, V.; Tripathi, S.; et al. A genome-wide SNP scan accelerates trait-regulatory genomic loci identification in chickpea. Sci. Rep. 2015, 5, 11166. [Google Scholar] [CrossRef] [Green Version]
  27. Khajuria, Y.P.; Saxena, M.S.; Gaur, R.; Chattopadhyay, D.; Jain, M.; Parida, S.K.; Bhatia, S. Development and Integration of Genome-Wide Polymorphic Microsatellite Markers onto a Reference Linkage Map for Constructing a High-Density Genetic Map of Chickpea. PLoS ONE 2015, 10, e0125583. [Google Scholar] [CrossRef] [Green Version]
  28. Jhanwar, S.; Priya, P.; Garg, R.; Parida, S.K.; Tyagi, A.K.; Jain, M. Transcriptome sequencing of wild chickpea as a rich resource for marker development. Plant Biotechnol. J. 2012, 10, 690–702. [Google Scholar] [CrossRef] [PubMed]
  29. Deokar, A.; Sagi, M.; Daba, K.; Tar’An, B. QTL sequencing strategy to map genomic regions associated with resistance to ascochyta blight in chickpea. Plant Biotechnol. J. 2018, 17, 275–288. [Google Scholar] [CrossRef] [PubMed]
  30. Abbo, S.; Berger, J.; Turner, N. Viewpoint: Evolution of cultivated chickpea: Four bottlenecks limit diversity andconstrain adaptation. Funct. Plant Biol. 2003, 30, 1081–1087. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Trapnell, C.; Pachter, L.; Salzberg, S. TopHat: Discovering splice junctions with RNA-Seq. Bioinformatics 2009, 25, 1105–1111. [Google Scholar] [CrossRef] [PubMed]
  32. Li, H.; Handsaker, R.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R. The Sequence Alignment/Map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Danecek, P.; Auton, A.; Abecasis, G.; Albers, C.A.; Banks, E.; DePristo, M.A.; Handsaker, R.; Lunter, G.; Marth, G.T.; Sherry, S.T.; et al. The variant call format and VCFtools. Bioinformatics 2011, 27, 2156–2158. [Google Scholar] [CrossRef] [PubMed]
  34. Shinozuka, H.; Cogan, N.O.; Shinozuka, M.; Marshall, A.; Kay, P.; Lin, Y.-H.; Spangenberg, G.C.; Forster, J.W. A simple method for semi-random DNA amplicon fragmentation using the methylation-dependent restriction enzyme MspJI. BMC Biotechnol. 2015, 15, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Li, H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv 2013, arXiv:1303.3997 2013. [Google Scholar]
  36. Javid, M.; Rosewarne, G.M.; Sudheesh, S.; Kant, P.; Leonforte, A.; Lombardi, M.; Kennedy, P.R.; Cogan, N.O.I.; Slater, A.T.; Kaur, S. Validation of molecular markers associated with boron tolerance, powdery mildew resistance and salinity tolerance in field peas. Front. Plant Sci. 2015, 6, 917. [Google Scholar] [CrossRef] [Green Version]
  37. Meng, L.; Li, H.; Zhang, L.; Wang, J. QTL IciMapping: Integrated software for genetic linkage map construction and quantitative trait locus mapping in biparental populations. Crop. J. 2015, 3, 269–283. [Google Scholar] [CrossRef] [Green Version]
  38. Voorrips, R.E. MapChart: Software for the Graphical Presentation of Linkage Maps and QTLs. J. Hered. 2002, 93, 77–78. [Google Scholar] [CrossRef] [Green Version]
  39. Bayer, M.; Milne, I.; Stephen, G.; Shaw, P.; Cardle, L.; Wright, F.; Marshall, D. Comparative visualization of genetic and physical maps with Strudel. Bioinformatics 2011, 27, 1307–1308. [Google Scholar] [CrossRef] [Green Version]
  40. Payne, R.; Baird, D.; Cherry, M.; Gilmour, A.; Harding, S.; Lane, P.; Morgan, G.; Murray, D.; Soutar, D.; Thompson, R. GenStat Release 6.1 Reference Manual. Part Directives; VSN International: Hemel Hempstead, UK, 2002. [Google Scholar]
  41. R Core Team. R: A Language and Environment for Statistical Computing; R Core Team: Vienna, Austria, 2020. [Google Scholar]
  42. Gilmour, A.R.; Gogel, B.J.; Cullis, B.R.; Welham, S.J.; Thompson, R. ASReml User Guide Release 4.1; VSN International Ltd.: Hemel Hempstead, UK, 2015. [Google Scholar]
  43. Wang, S.; Basten, C.; Zeng, Z. Windows QTL Cartographer 2.Department of Statistics; North Carolina State University: Raleigh, NC, USA, 2012. [Google Scholar]
  44. You, Q.; Yang, X.; Peng, Z.; Xu, L.; Wang, J. Development and Applications of a High Throughput Genotyping Tool for Polyploid Crops: Single Nucleotide Polymorphism (SNP) Array. Front. Plant Sci. 2018, 9, 104. [Google Scholar] [CrossRef] [Green Version]
  45. Lutz, S.; Brion, C.; Kliebhan, M.; Albert, F.W. DNA variants affecting the expression of numerous genes in trans have diverse mechanisms of action and evolutionary histories. PLoS Genet. 2019, 15, e1008375. [Google Scholar] [CrossRef] [Green Version]
  46. Reuter, J.A.; Spacek, D.V.; Snyder, M.P. High-Throughput Sequencing Technologies. Mol. Cell 2015, 58, 586–597. [Google Scholar] [CrossRef] [Green Version]
  47. Fuchsberger, C.; Abecasis, G.R.; Hinds, D. minimac2: Faster genotype imputation. Bioinformatics 2014, 31, 782–784. [Google Scholar] [CrossRef] [Green Version]
  48. Browning, B.L.; Zhou, Y.; Browning, S. A One-Penny Imputed Genome from Next-Generation Reference Panels. Am. J. Hum. Genet. 2018, 103, 338–348. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Money, D.; Gardner, K.; Migicovsky, Z.; Schwaninger, H.; Zhong, G.-Y.; Myles, S. LinkImpute: Fast and Accurate Genotype Imputation for Nonmodel Organisms. G3: Genes Genomes Genet. 2015, 5, 2383–2390. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Shi, S.; Yuan, N.; Yang, M.; Du, Z.; Wang, J.; Sheng, X.; Wu, J.; Xiao, J. Comprehensive Assessment of Genotype Imputation Performance. Hum. Hered. 2018, 83, 107–116. [Google Scholar] [CrossRef] [PubMed]
  51. Pootakham, W.; Shearman, J.R.; Ruang-Areerate, P.; Sonthirod, C.; Sangsrakru, D.; Jomchai, N.; Yoocha, T.; Triwitayakorn, K.; Tragoonrung, S.; Tangphatsornruang, S. Large-Scale SNP Discovery through RNA Sequencing and SNP Genotyping by Targeted Enrichment Sequencing in Cassava (Manihot esculenta Crantz). PLoS ONE 2014, 9, e116028. [Google Scholar] [CrossRef] [Green Version]
  52. Kaur, L.; Sirari, A.; Kumar, D.; Singh Sandhu, J.; Singh, S.; Kapoor, K.; Singh, I.; Gowda, C.L.L.; Pande, S.; Gaur, P. Combining Ascochyta blight and Botrytis grey mould resistance in chickpea through interspecific hybridization. Phytopathol. Mediterr. 2013, 52, 157–165. [Google Scholar]
  53. Tewari, S.K.; Pandey, M.P. Genetics of resistance to ascochyta blight in chickpea (Cicer arietinum L.). Euphytica 1986, 35, 211–215. [Google Scholar] [CrossRef]
  54. Singh, V.; Khan, A.W.; Saxena, R.K.; Kumar, V.; Kale, S.; Sinha, P.; Chitikineni, A.; Pazhamala, L.T.; Garg, V.; Sharma, M.; et al. Next-generation sequencing for identification of candidate genes for Fusarium wilt and sterility mosaic disease in pigeonpea (Cajanus cajan). Plant Biotechnol. J. 2015, 14, 1183–1194. [Google Scholar] [CrossRef] [Green Version]
  55. Luo, H.; Pandey, M.K.; Khan, A.W.; Wu, B.; Guo, J.; Ren, X.; Zhou, X.; Chen, Y.; Chen, W.; Huang, L.; et al. Next-generation sequencing identified genomic region and diagnostic markers for resistance to bacterial wilt on chromosome B02 in peanut (Arachis hypogaea L.). Plant Biotechnol. J. 2019, 17, 2356–2369. [Google Scholar] [CrossRef] [Green Version]
  56. Magwanga, R.O.; Lu, P.; Kirungu, J.N.; Diouf, L.; Dong, Q.; Hu, Y.; Cai, X.; Xu, Y.; Hou, Y.; Zhou, Z.; et al. GBS Mapping and Analysis of Genes Conserved between Gossypium tomentosum and Gossypium hirsutum Cotton Cultivars that Respond to Drought Stress at the Seedling Stage of the BC2F2 Generation. Int. J. Mol. Sci. 2018, 19, 1614. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Sari, E.; Bhadauria, V.; Ramsay, L.; Borhan, M.H.; Lichtenzveig, J.; Bett, K.E.; Vandenberg, A.; Banniza, S. Defense responses of lentil (Lens culinaris) genotypes carrying non-allelic ascochyta blight resistance genes to Ascochyta lentis infection. PLoS ONE 2018, 13, e0204124. [Google Scholar] [CrossRef] [Green Version]
  58. Fikere, M.; Barbulescu, D.M.; Malmberg, M.M.; Shi, F.; Koh, J.C.; Slater, A.T.; MacLeod, I.M.; Bowman, P.J.; Salisbury, P.A.; Spangenberg, G.C.; et al. Genomic Prediction Using Prior Quantitative Trait Loci Information Reveals a Large Reservoir of Underutilised Blackleg Resistance in Diverse Canola (Brassica napus L.) Lines. Plant Genome 2018, 11, 170100. [Google Scholar] [CrossRef] [Green Version]
  59. Olatoye, M.O.; Clark, L.V.; Wang, J.; Yang, X.; Yamada, T.; Sacks, E.J.; Lipka, A.E. Evaluation of genomic selection and marker-assisted selection in Miscanthus and energycane. Mol. Breed. 2019, 39, 171. [Google Scholar] [CrossRef] [Green Version]
  60. Jannink, J.-L. Dynamics of long-term genomic selection. Genet. Sel. Evol. 2010, 42, 35. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Etayeh, N.; Eklein, A.; Paslier, M.-C.E.; Ejacquin, F.; Ehoutin, H.; Erond, C.; Echabert-Martinello, M.; Emagnin-Robert, J.-B.; Emarget, P.; Aubert, G.; et al. Genomic Prediction in Pea: Effect of Marker Density and Training Population Size and Composition on Prediction Accuracy. Front. Plant Sci. 2015, 6, 941. [Google Scholar] [CrossRef] [Green Version]
  62. Spindel, J.; Begum, H.; Akdemir, D.; Virk, P.; Collard, B.; Redoña, E.; Atlin, G.; Jannink, J.-L.; McCouch, S.R. Genomic Selection and Association Mapping in Rice (Oryza sativa): Effect of Trait Genetic Architecture, Training Population Composition, Marker Number and Statistical Model on Accuracy of Rice Genomic Selection in Elite, Tropical Rice Breeding Lines. PLoS Genet. 2015, 11, e1004982. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Genetic linkage map of the Sonali X 04067-81-2-1-1 cross, showing the location of genomic region for AB resistance. The name is provided at the top of each LG. Distances of the loci (cM) are shown to the left and names of loci are shown to the right side of LGs. For presentation purposes, only selected markers are shown on the map.
Figure 1. Genetic linkage map of the Sonali X 04067-81-2-1-1 cross, showing the location of genomic region for AB resistance. The name is provided at the top of each LG. Distances of the loci (cM) are shown to the left and names of loci are shown to the right side of LGs. For presentation purposes, only selected markers are shown on the map.
Agronomy 11 01937 g001
Table 1. Marker distribution over the LGs of Sonali X 04067-81-2-1-1 genetic linkage map.
Table 1. Marker distribution over the LGs of Sonali X 04067-81-2-1-1 genetic linkage map.
Linkage GroupNumber of Mapped MarkersCumulative Map Length (cM)Average Marker Density
LG193662.230.07
LG220676.050.37
LG2.22116.080.77
LG33132.561.05
LG3.250136.540.07
LG41934105.760.05
LG5105110.831.06
LG6103156.810.06
LG6.2337.680.23
LG7676148.820.22
LG841263.900.16
Total5886717.260.12
Table 2. Identification of QTLs for AB resistance on Sonali X 04067-81-2-1-1 RIL population genetic maps based on CIM.
Table 2. Identification of QTLs for AB resistance on Sonali X 04067-81-2-1-1 RIL population genetic maps based on CIM.
QTL NameChromosomeFlanking MarkersPosition (cM)LOD ThresholdMaximum LOD Threshold% Phenotypic Variance Additive Effect
AB_echino_2014 LG4Ca_Ce_1844548.154.513.7340.45
Ca_Ce_18594, a_Ce_1857748.68
Ca_Ce_1865649.14
AB_echino_2015LG4Ca_Ce_1844548.15416.0410.40
Ca_Ce_18594, a_Ce_1857748.68
Ca_Ce_1865649.14
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sudheesh, S.; Kahrood, H.V.; Braich, S.; Dron, N.; Hobson, K.; Cogan, N.O.I.; Kaur, S. Application of Genomics Approaches for the Improvement in Ascochyta Blight Resistance in Chickpea. Agronomy 2021, 11, 1937. https://doi.org/10.3390/agronomy11101937

AMA Style

Sudheesh S, Kahrood HV, Braich S, Dron N, Hobson K, Cogan NOI, Kaur S. Application of Genomics Approaches for the Improvement in Ascochyta Blight Resistance in Chickpea. Agronomy. 2021; 11(10):1937. https://doi.org/10.3390/agronomy11101937

Chicago/Turabian Style

Sudheesh, Shimna, Hossein V. Kahrood, Shivraj Braich, Nicole Dron, Kristy Hobson, Noel O. I. Cogan, and Sukhjiwan Kaur. 2021. "Application of Genomics Approaches for the Improvement in Ascochyta Blight Resistance in Chickpea" Agronomy 11, no. 10: 1937. https://doi.org/10.3390/agronomy11101937

APA Style

Sudheesh, S., Kahrood, H. V., Braich, S., Dron, N., Hobson, K., Cogan, N. O. I., & Kaur, S. (2021). Application of Genomics Approaches for the Improvement in Ascochyta Blight Resistance in Chickpea. Agronomy, 11(10), 1937. https://doi.org/10.3390/agronomy11101937

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop