Characterization and mapping of sequence-tagged microsatellite sites in the chickpea (Cicer arietinum L.) genome

Meta QTL analysis for dissecting abiotic stress tolerance in chickpea

BACKGROUND

Chickpea is prone to many abiotic stresses such as heat, drought, salinity, etc. which cause severe loss in yield. Tolerance towards these stresses is quantitative in nature and many studies have been done to map the loci influencing these traits in different populations using different markers. This study is an attempt to meta-analyse those reported loci projected over a high-density consensus map to provide a more accurate information on the regions influencing heat, drought, cold and salinity tolerance in chickpea.

RESULTS

A meta-analysis of QTL reported to be responsible for tolerance to drought, heat, cold and salinity stress tolerance in chickpeas was done. A total of 1512 QTL responsible for the concerned abiotic stress tolerance were collected from literature, of which 1189 were projected on a chickpea consensus genetic map. The QTL meta-analysis predicted 59 MQTL spread over all 8 chromosomes, responsible for these 4 kinds of abiotic stress tolerance in chickpea. The physical locations of 23 MQTL were validated by various marker-trait associations and genome-wide association studies. Out of these reported MQTL, CaMQAST1.1, CaMQAST4.1, CaMQAST4.4, CaMQAST7.8, and CaMQAST8.2 were suggested to be useful for different breeding approaches as they were responsible for high per cent variance explained (PVE), had small intervals and encompassed a large number of originally reported QTL. Many putative candidate genes that might be responsible for directly or indirectly conferring abiotic stress tolerance were identified in the region covered by 4 major MQTL- CaMQAST1.1, CaMQAST4.4, CaMQAST7.7, and CaMQAST6.4, such as heat shock proteins, auxin and gibberellin response factors, etc. CONCLUSION: The results of this study should be useful for the breeders and researchers to develop new chickpea varieties which are tolerant to drought, heat, cold, and salinity stresses.

Comparison of different screening methods for the selection of Ascochyta blight disease on chickpea (Cicer arietinum L.) genotypes

Chickpea (Cicer arietinum L.) is the second most important edible food grain legume, widely grown all over the world. However, the cultivation and production of chickpea are mainly affected by the Ascochyta blight (AB) disease, which causes losses of up to 100% in areas with high humidity and warm temperature conditions. Various screening methods are used in the selection of chickpea genotypes for resistance to AB disease. These methods are natural field condition (NFC), artificial epidemic field condition (AEC), marker-assisted selection (MAS), and real-time PCR (RT-PCR). The study was conducted with 88 chickpea test genotypes between the 2014 and 2016 growing seasons. The results of the screening were used to sort the genotypes into three categories: susceptible (S), moderately resistant (MR), and resistant (R). Using MAS screening, 13, 21, and 54 chickpea genotypes were identified as S, MR, and R, respectively. For RT-PCR screening, 39 genotypes were S, 31 genotypes were MR, and 18 genotypes were R. In the AEC method for NFC screening, 7, 17, and 64 genotypes were S, MR, and R, while 74 and 6 genotypes were S and MR, and 8 genotypes were R-AB disease. As a result of screening chickpea genotypes for AB disease, it was determined that the most effective method was artificial inoculation (AEC) under field conditions. In the study, Azkan, ICC3996, Tüb-19, and Tüb-82 were determined as resistant within all methods for Pathotype 1.

Functional allele of a MATE gene selected during domestication modulates seed color in chickpea

Seed color is one of the key target traits of domestication and artificial selection in chickpeas due to its implications on consumer preference and market value. The complex seed color trait has been well dissected in several crop species; however, the genetic mechanism underlying seed color variation in chickpea remains poorly understood. Here, we employed an integrated genomics strategy involving QTL mapping, high-density mapping, map-based cloning, association analysis, and molecular haplotyping in an inter-specific RIL mapping population, association panel, wild accessions, and introgression lines (ILs) of Cicer gene pool. This delineated a MATE gene, CaMATE23, encoding a Transparent Testa (TT) and its natural allele (8-bp insertion) and haplotype underlying a major QTL governing seed color on chickpea chromosome 4. Signatures of selective sweep and a strong purifying selection reflected that CaMATE23, especially its 8-bp insertion natural allelic variant, underwent selection during chickpea domestication. Functional investigations revealed that the 8-bp insertion containing the third cis-regulatory RY-motif element in the CaMATE23 promoter is critical for enhanced binding of CaFUSCA3 transcription factor, a key regulator of seed development and flavonoid biosynthesis, thereby affecting CaMATE23 expression and proanthocyanidin (PA) accumulation in the seed coat to impart varied seed color in chickpea. Consequently, overexpression of CaMATE23 in Arabidopsis tt12 mutant partially restored the seed color phenotype to brown pigmentation, ascertaining its functional role in PA accumulation in the seed coat. These findings shed new light on the seed color regulation and evolutionary history, and highlight the transcriptional regulation of CaMATE23 by CaFUSCA3 in modulating seed color in chickpea. The functionally relevant InDel variation, natural allele, and haplotype from CaMATE23 are vital for translational genomic research, including marker-assisted breeding, for developing chickpea cultivars with desirable seed color that appeal to consumers and meet global market demand.

SSR Genotyping and Marker-Trait Association with Yield Components in a Kazakh Germplasm Collection of Chickpea (Cicer arietinum L.)

Genetic diversity and marker-trait association with yield-related components were assessed in 39 chickpea accessions from a germplasm collection with either spring or autumn-sown seeds in South-Eastern Kazakhstan. Chickpea accessions originated from Azerbaijan, Germany, Kazakhstan, Moldova, Russia, Türkiye, Ukraine, Syria, and the International Center for Agricultural Research in the Dry Areas (ICARDA). Eleven SSR markers were used for molecular genotyping. Yield and yield components were evaluated in nine traits in experiments with spring and autumn seed sowing. The number of alleles of polymorphic markers varied from 2 to 11. The greatest polymorphism was found in the studied chickpea genotypes using SSR marker TA22 (11 alleles), while NCPGR6 and NCPGR12 markers were monomorphic. In the studied chickpea accessions, unique alleles of the SSR loci TA14, TA46, TA76s, and TA142 were found that were not previously described by other authors. An analysis of correlation relationships between yield-related traits in chickpea revealed the dependence of yield on plant height, branching, and the setting of a large number of beans. These traits showed maximal values in experiments with chickpea plants from autumn seed sowing. An analysis of the relationship between the SSR markers applied and morphological yield-related traits revealed several informative markers associated with important traits, such as plant height, height to first pod, number of branches, number of productive nodes, number of pods per plant, hundred seed weight, seed weight per plant, and seed yield.

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.

Genetic variability and population structure of Ethiopian chickpea (Cicer arietinum L.) germplasm

Evaluation of the genetic diversity and an understanding of the genetic structure and relationships of chickpea genotypes are valuable to design efficient germplasm conservation strategies and crop breeding programs. Information is limited, in these regards, for Ethiopian chickpea germplasms. Therefore, the present study was carried out to estimate the genetic diversity, population structure, and relationships of 152 chickpea genotypes using simple sequence repeats (SSR) markers. Twenty three SSR markers exhibited polymorphism producing a total of 133 alleles, with a mean of 5.8 alleles per locus. Analyses utilizing various genetic-based statistics included pairwise population Nei’s genetic distance, heterozygosity, Shannon’s information index, polymorphic information content, and percent polymorphism. These analyses exemplified the existence of high genetic variation within and among chickpea genotypes. The 152 genotypes were divided into two major clusters based on Nei’s genetic distances. The exotic genotypes were grouped in one cluster exclusively showing that these genotypes are distinct to Ethiopian genotypes, while the patterns of clustering of Ethiopian chickpea genotypes based on their geographic region were not consistent because of the seed exchange across regions. Model-based population structure clustering identified two discrete populations. These finding provides useful insight for chickpea collections and ex-situ conservation and national breeding programs for widening the genetic base of chickpea.

Unraveling genetics of semi-determinacy and identification of markers for indeterminate stem growth habit in chickpea (Cicer arietinum L.)

Chickpea (Cicer arietinum L.) is predominantly an indeterminate plant and tends to generate vegetative growth when the ambient is conducive for soil moisture, temperature and certain other environmental conditions. The semi-determinate (SDT) types are comparatively early, resistant to lodging and found to be similar in their yield potential to indeterminate (IDT) lines. Indeterminate and semi-determinate genotypes are found to be similar during early stage, which makes it difficult to distinguish between them. Thus, there is a need to identify molecular markers linked either to indeterminate or semi-determinate plant types. The present study was carried out to study the genetics of semi-determinacy and identify molecular markers linked to stem growth habit. The study was undertaken in the cross involving BG 362(IDT) × BG 3078-1(SDT). All F₁ plants were indeterminate, which indicates that indeterminate stem type is dominant over semi-determinate. In further advancement to F₂ generation, F₂ plants are segregated in the ratio of 3(Indeterminate): 1(Semi-determinate) that indicates that the IDT and SDT parents which are involved in the cross differed for a single gene. The segregation pattern observed in F₂ is confirmed in F₃ generation. The parental polymorphic survey was undertaken for molecular analysis using total of 245 SSR markers, out of which 41 polymorphic markers were found to distinguish the parents and were utilized for bulked segregant analysis (BSA). The segregation pattern in F₂ indicates that the IDT (Indeterminate) and SDT (Semi-determinate) parents which are involved in the cross differed for single gene. The segregation pattern of F₂ and F₃ derived from the cross BG 362 (IDT) × BG 3078-1 (SDT) confirmed the genotypic structure of the newly found SDT genotype BG 3078-1 as dt1dt1Dt2Dt2. Three SSR markers TA42, Ca_GPSSR00560 and H3DO5 were found to be putatively linked to Dt1 locus regulating IDT stem growth habit. Our results indicate that the SSR markers identified for Dt1 locus helps to differentiate stem growth habit of chickpea in its early growth stage itself and can be efficiently utilized in Marker Assisted Selection (MAS) for changed plant type in chickpea.

Cross-genera amplification of Cajanus spp. specific SSR markers in Clitoria ternatea (L.) and their application in genetic diversity studies

Clitoria ternatea (L.) is a medicinal leguminous plant and is cultivated to cater the need of herbal industries and asthetic purposes. The unavailability of steady molecular marker impedes the genetic improvement of C. ternatea. In the present study, transferability of 98 pairs of Cajanus spp. specific SSR primers were assessed among 14 genotypes of C. ternatea, varied for their flower color, floral architecture and bio-metabolite (taraxerol and delphinidin) content, and out of them 43 had successfully amplified the fragments. Among them, 36 pairs of primers showed 100% transferability, whereas rest seven varied from 42.86 to 92.85% transferability. The transferable 43 pairs of SSR primers generated 196 alleles across the 14 genotypes and the AMOVA analysis showed moderate genetic variation (55.1%) among the genotypes of C. ternatea, which was also reinforced by Nei's genetic distance and gene identity estimates derived haplotype matrix. Similarly, both the principal coordinate analysis and dendrogram grouped these 14 genotypes of C. ternatea into two major clusters based on SSR allele distribution and frequency, and the clustering pattern is in accordance with petal color but in contrast to floral architecture. MCheza based outlier analysis revealed 16 alleles for balancing selection, which are putatively involved in the maintenance of genetic polymorphism in C. ternatea. Moreover, the estimates of molecular diversity and bio-metabolite content revealed the possible use of these genotypes in future breeding programme of this species.

Population structure and association analysis of heat stress relevant traits in chickpea (Cicer arietinum L.)

Understanding genetic diversity and population structure is prerequisite to broaden the cultivated base of any crop. In the current investigation, we report discovery of a total of 319 alleles by assaying 81 SSRs on 71 chickpea genotypes. The cluster analysis based on Jaccard coefficient and unweighted neighbor joining algorithm categorized all genotypes into two major clusters. Cultivars grown within the same agro-climatic zones were clustered together, whereas the remaining genotypes particularly advanced breeding lines and accessions assigned to another cluster. Population structure analysis separated the entire collection into two subpopulations (K = 2) and the clustering pattern remained in close agreement with those of distance-based methods. Importantly, we also discovered marker trait association for membrane stability index (MSI) and leaf chlorophyll content measured as SPAD chlorophyll meter reading (SCMR), the two important physiological parameters indicative of heat stress (HS) tolerance in chickpea. Association analysis using both general linear and mixed linear models of the mean phenotypic data of traits recorded in 2016 and 2017 uncovered significant association of NCPGR206 and H2L102 with the MSI trait. Likewise, SSR markers GA9, TR31 and TA113 exhibited significant association with SCMR trait. The genomic regions putatively linked with two traits may be investigated in greater detail to further improve knowledge about the genetic architecture of HS tolerance in chickpea.

Population structure and association analysis of heat stress relevant traits in chickpea (Cicer arietinum L.)

Understanding genetic diversity and population structure is prerequisite to broaden the cultivated base of any crop. In the current investigation, we report discovery of a total of 319 alleles by assaying 81 SSRs on 71 chickpea genotypes. The cluster analysis based on Jaccard coefficient and unweighted neighbor joining algorithm categorized all genotypes into two major clusters. Cultivars grown within the same agro-climatic zones were clustered together, whereas the remaining genotypes particularly advanced breeding lines and accessions assigned to another cluster. Population structure analysis separated the entire collection into two subpopulations (K = 2) and the clustering pattern remained in close agreement with those of distance-based methods. Importantly, we also discovered marker trait association for membrane stability index (MSI) and leaf chlorophyll content measured as SPAD chlorophyll meter reading (SCMR), the two important physiological parameters indicative of heat stress (HS) tolerance in chickpea. Association analysis using both general linear and mixed linear models of the mean phenotypic data of traits recorded in 2016 and 2017 uncovered significant association of NCPGR206 and H2L102 with the MSI trait. Likewise, SSR markers GA9, TR31 and TA113 exhibited significant association with SCMR trait. The genomic regions putatively linked with two traits may be investigated in greater detail to further improve knowledge about the genetic architecture of HS tolerance in chickpea.

DRDB: An Online Date Palm Genomic Resource Database

Background: Date palm (Phoenix dactylifera L.) is a cultivated woody plant with agricultural and economic importance in many countries around the world. With the advantages of next generation sequencing technologies, genome sequences for many date palm cultivars have been released recently. Short sequence repeat (SSR) and single nucleotide polymorphism (SNP) can be identified from these genomic data, and have been proven to be very useful biomarkers in plant genome analysis and breeding. Results: Here, we first improved the date palm genome assembly using 130X of HiSeq data generated in our lab. Then 246,445 SSRs (214,901 SSRs and 31,544 compound SSRs) were annotated in this genome assembly; among the SSRs, mononucleotide SSRs (58.92%) were the most abundant, followed by di- (29.92%), tri- (8.14%), tetra- (2.47%), penta- (0.36%), and hexa-nucleotide SSRs (0.19%). The high-quality PCR primer pairs were designed for most (174,497; 70.81% out of total) SSRs. We also annotated 6,375,806 SNPs with raw read depth≥3 in 90% cultivars. To further reduce false positive SNPs, we only kept 5,572,650 (87.40% out of total) SNPs with at least 20% cultivars support for downstream analyses. The high-quality PCR primer pairs were also obtained for 4,177,778 (65.53%) SNPs. We reconstructed the phylogenetic relationships among the 62 cultivars using these variants and found that they can be divided into three clusters, namely North Africa, Egypt - Sudan, and Middle East - South Asian, with Egypt - Sudan being the admixture of North Africa and Middle East - South Asian cultivars; we further confirmed these clusters using principal component analysis. Moreover, 34,346 SSRs and 4,177,778 SNPs with PCR primers were assigned to shared cultivars for cultivar classification and diversity analysis. All these SSRs, SNPs and their classification are available in our database, and can be used for cultivar identification, comparison, and molecular breeding. Conclusion:DRDB is a comprehensive genomic resource database of date palm. It can serve as a bioinformatics platform for date palm genomics, genetics, and molecular breeding. DRDB is freely available at http://drdb.big.ac.cn/home.

Marker-assisted introgression of resistance to fusarium wilt race 2 in Pusa 256, an elite cultivar of desi chickpea

Fusarium wilt caused by F. oxysporum f. sp. ciceris causes extensive damage to chickpea (Cicer arietinum L.) in many parts of the world. In the central part of India, pathogen race 2 (Foc 2) causes severe yield losses. We initiated molecular marker-assisted backcrossing (MABC) using desi cultivar, Vijay, as a donor to introgress resistance to this race (Foc2) in Pusa 256, another elite desi cultivar of chickpea. To confirm introgression of resistance for this race, foreground selection was undertaken using two SSR markers (TA 37 and TA110), with background selection to observe the recovery of recurrent parent genome using 45 SSRs accommodated in 8 multiplexes. F₁ plants were confirmed with molecular markers and backcrossed with Pusa 256, followed by cycles of foreground and background selection at each stage to generate 161 plants in BC₃F₂ during the period 2009-2013. Similarly, 46 BC₃F₁ plants were also generated in another set during the same period. On the basis of foreground selection, 46 plants were found homozygotes in BC₃F₂. Among them, 17 plants recorded >91% background recovery with the highest recovery percentage of 96%. In BC₃F₁ also, 14 hybrid plants recorded a background recovery of >85% with the highest background recovery percentage of >94%. The identified plants were selfed to obtain 1341 BC₃F₃ and 2198 BC₃F₂ seeds which were screened phenotypically for resistance to fusarium wilt (race 2) besides doing marker analysis. Finally, 17 BC₃F₄ and 11 BC₃F₃ lines were obtained which led to identification of 5 highly resistant lines of Pusa 256 with Foc 2 gene introgressed in them. Development of these lines will help in horizontal as well as vertical expansion of chickpea in central part of India.

Draft genome sequence of Cicer reticulatum L., the wild progenitor of chickpea provides a resource for agronomic trait improvement

Cicer reticulatum L. is the wild progenitor of the fourth most important legume crop chickpea (C. arietinum L.). We assembled short-read sequences into 416 Mb draft genome of C. reticulatum and anchored 78% (327 Mb) of this assembly to eight linkage groups. Genome annotation predicted 25,680 protein-coding genes covering more than 90% of predicted gene space. The genome assembly shared a substantial synteny and conservation of gene orders with the genome of the model legume Medicago truncatula. Resistance gene homologs of wild and domesticated chickpeas showed high sequence homology and conserved synteny. Comparison of gene sequences and nucleotide diversity using 66 wild and domesticated chickpea accessions suggested that the desi type chickpea was genetically closer to the wild species than the kabuli type. Comparative analyses predicted gene flow between the wild and the cultivated species during domestication. Molecular diversity and population genetic structure determination using 15,096 genome-wide single nucleotide polymorphisms revealed an admixed domestication pattern among cultivated (desi and kabuli) and wild chickpea accessions belonging to three population groups reflecting significant influence of parentage or geographical origin for their cultivar-specific population classification. The assembly and the polymorphic sequence resources presented here would facilitate the study of chickpea domestication and targeted use of wild Cicer germplasms for agronomic trait improvement in chickpea.

Molecular Mapping of Flowering Time Major Genes and QTLs in Chickpea (Cicer arietinum L.)

Flowering time is an important trait for adaptation and productivity of chickpea in the arid and the semi-arid environments. This study was conducted for molecular mapping of genes/quantitative trait loci (QTLs) controlling flowering time in chickpea using F₂ populations derived from four crosses (ICCV 96029 × CDC Frontier, ICC 5810 × CDC Frontier, BGD 132 × CDC Frontier and ICC 16641 × CDC Frontier). Genetic studies revealed monogenic control of flowering time in the crosses ICCV 96029 × CDC Frontier, BGD 132 × CDC Frontier and ICC 16641 × CDC Frontier, while digenic control with complementary gene action in ICC 5810 × CDC Frontier. The intraspecific genetic maps developed from these crosses consisted 75, 75, 68 and 67 markers spanning 248.8 cM, 331.4 cM, 311.1 cM and 385.1 cM, respectively. A consensus map spanning 363.8 cM with 109 loci was constructed by integrating four genetic maps. Major QTLs corresponding to flowering time genes efl-1 from ICCV 96029, efl-3 from BGD 132 and efl-4 from ICC 16641 were mapped on CaLG04, CaLG08 and CaLG06, respectively. The QTLs and linked markers identified in this study can be used in marker-assisted breeding for developing early maturing chickpea.

Fine mapping for double podding gene in chickpea

For the first time, fine mapping for sfl locus was carried out using a battery of new STMS and SNP markers. The target region was delimited to 92.6 Kb where seven annotated genes were found that could be candidate genes for the simple/double podding trait in chickpea. Four recombinant inbred populations (RIP-1, RIP-7, RIP-11, and CPR-01) were used to map the double podding gene (sfl) in chickpea. In RIP-1, the gene was initially mapped on linkage group (LG) 6 between the two sequence-tagged microsatellite site (STMS) markers TA120 and TR1. Eight new STMS markers were added onto LG6 in the target region and sfl locus was finally located between CAGM27819 and CAGM27777 markers within an interval of 2 cM. Seven out of the eight markers were mapped in RIP-7 and its reciprocal RIP-11 confirming the location of the sfl locus to a 4.8 cM interval flanked by TR44 and CAGM27705. Furthermore, using a high-density single nucleotide polymorphism (SNP) map of CPR-01, sfl was mapped to the same genomic region in a 5.1 cM interval between TR44 and the SNP scaffold1646p97220. Five pairs of near isogenic lines (NILs) and eight recombinant inbred lines (RILs) were used to refine this region in the chickpea physical map. Combining data from linkage analysis in four RIPs, marker physical positions and recombination events obtained in both pairs of NILs and selected RILs, sfl could be placed within a genomic window of 92.6 Kb. Seven annotated genes were extracted from this region. The regulator of axillary meristem-predicted gene could be a candidate gene for the simple/double podding gene. This study provides additional set of markers flanking and tightly linked to sfl locus that are useful for marker-assisted selection.

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.

A genome-scale integrated approach aids in genetic dissection of complex flowering time trait in chickpea

A combinatorial approach of candidate gene-based association analysis and genome-wide association study (GWAS) integrated with QTL mapping, differential gene expression profiling and molecular haplotyping was deployed in the present study for quantitative dissection of complex flowering time trait in chickpea. Candidate gene-based association mapping in a flowering time association panel (92 diverse desi and kabuli accessions) was performed by employing the genotyping information of 5724 SNPs discovered from 82 known flowering chickpea gene orthologs of Arabidopsis and legumes as well as 832 gene-encoding transcripts that are differentially expressed during flower development in chickpea. GWAS using both genome-wide GBS- and candidate gene-based genotyping data of 30,129 SNPs in a structured population of 92 sequenced accessions (with 200-250 kb LD decay) detected eight maximum effect genomic SNP loci (genes) associated (34% combined PVE) with flowering time. Six flowering time-associated major genomic loci harbouring five robust QTLs mapped on a high-resolution intra-specific genetic linkage map were validated (11.6-27.3% PVE at 5.4-11.7 LOD) further by traditional QTL mapping. The flower-specific expression, including differential up- and down-regulation (>three folds) of eight flowering time-associated genes (including six genes validated by QTL mapping) especially in early flowering than late flowering contrasting chickpea accessions/mapping individuals during flower development was evident. The gene haplotype-based LD mapping discovered diverse novel natural allelic variants and haplotypes in eight genes with high trait association potential (41% combined PVE) for flowering time differentiation in cultivated and wild chickpea. Taken together, eight potential known/candidate flowering time-regulating genes [efl1 (early flowering 1), FLD (Flowering locus D), GI (GIGANTEA), Myb (Myeloblastosis), SFH3 (SEC14-like 3), bZIP (basic-leucine zipper), bHLH (basic helix-loop-helix) and SBP (SQUAMOSA promoter binding protein)], including novel markers, QTLs, alleles and haplotypes delineated by aforesaid genome-wide integrated approach have potential for marker-assisted genetic improvement and unravelling the domestication pattern of flowering time in chickpea.

Development of gene-based markers for use in construction of the chickpea (Cicer arietinum L.) genetic linkage map and identification of QTLs associated with seed weight and plant height

Seed weight and plant height are important agronomic traits and contribute to seed yield. The objective of this study was to identify QTLs underlying these traits using an intra-specific mapping population of chickpea. A F11 population of 177 recombinant inbred lines derived from a cross between SBD377 (100-seed weight--48 g and plant height--53 cm) and BGD112 (100-seed weight--15 g and plant height--65 cm) was used. A total of 367 novel EST-derived functional markers were developed which included 187 EST-SSRs, 130 potential intron polymorphisms (PIPs) and 50 expressed sequence tag polymorphisms (ESTPs). Along with these, 590 previously published markers including 385 EST-based markers and 205 genomic SSRs were utilized. Of the 957 markers tested for analysis of parental polymorphism between the two parents of the mapping population, 135 (14.64%) were found to be polymorphic. Of these, 131 polymorphic markers could be mapped to the 8 linkage groups. The linkage map had a total length of 1140.54 cM with an average marker density of 8.7 cM. The map was further used for QTL identification using composite interval mapping method (CIM). Two QTLs each for seed weight, qSW-1 and qSW-2 (explaining 11.54 and 19.24% of phenotypic variance, respectively) and plant height, qPH-1 and qPH-2 (explaining 13.98 and 12.17% of phenotypic variance, respectively) were detected. The novel set of genic markers, the intra-specific linkage map and the QTLs identified in the present study will serve as valuable genomic resources in improving the chickpea seed yield using marker-assisted selection (MAS) strategies.

Deploying QTL-seq for rapid delineation of a potential candidate gene underlying major trait-associated QTL in chickpea

A rapid high-resolution genome-wide strategy for molecular mapping of major QTL(s)/gene(s) regulating important agronomic traits is vital for in-depth dissection of complex quantitative traits and genetic enhancement in chickpea. The present study for the first time employed a NGS-based whole-genome QTL-seq strategy to identify one major genomic region harbouring a robust 100-seed weight QTL using an intra-specific 221 chickpea mapping population (desi cv. ICC 7184 × desi cv. ICC 15061). The QTL-seq-derived major SW QTL (CaqSW1.1) was further validated by single-nucleotide polymorphism (SNP) and simple sequence repeat (SSR) marker-based traditional QTL mapping (47.6% R² at higher LOD >19). This reflects the reliability and efficacy of QTL-seq as a strategy for rapid genome-wide scanning and fine mapping of major trait regulatory QTLs in chickpea. The use of QTL-seq and classical QTL mapping in combination narrowed down the 1.37 Mb (comprising 177 genes) major SW QTL (CaqSW1.1) region into a 35 kb genomic interval on desi chickpea chromosome 1 containing six genes. One coding SNP (G/A)-carrying constitutive photomorphogenic9 (COP9) signalosome complex subunit 8 (CSN8) gene of these exhibited seed-specific expression, including pronounced differential up-/down-regulation in low and high seed weight mapping parents and homozygous individuals during seed development. The coding SNP mined in this potential seed weight-governing candidate CSN8 gene was found to be present exclusively in all cultivated species/genotypes, but not in any wild species/genotypes of primary, secondary and tertiary gene pools. This indicates the effect of strong artificial and/or natural selection pressure on target SW locus during chickpea domestication. The proposed QTL-seq-driven integrated genome-wide strategy has potential to delineate major candidate gene(s) harbouring a robust trait regulatory QTL rapidly with optimal use of resources. This will further assist us to extrapolate the molecular mechanism underlying complex quantitative traits at a genome-wide scale leading to fast-paced marker-assisted genetic improvement in diverse crop plants, including 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.

Evaluation of microsatellite-based genetic diversity, protein and mineral content in chickpea accessions grown in Kyrgyzstan

The genetic diversity of 23 chickpea accessions representing Kyrgyz landraces and cultivars, ICARDA breeding lines, Spanish and Turkish cultivars was characterized using nine microsatellite (SSR) markers which generated a total of 122 alleles. The number of alleles (Na) per locus varied from 9 to 20. The observed heterozygosity (Ho) ranged between 0.05 and 0.43 (average 0.13) whereas both the expected heterozygosity (He) and polymorphic information content (PIC) ranged from 0.71 to 0.90 (average 0.83). Analysis of molecular variance (AMOVA) showed that 62% of the total genetic variation was found within accessions while the remaining 38% was found among accessions. Principal coordinate analysis (PCoA) indicated the presence of two groups. The two Kyrgyz cultivars were found apart from these groups. Cluster analysis generally confirmed the results of PCoA and also separated the Kyrgyz cultivars from the subcluster formed by Kyrgyz landraces and the subclusters formed by breeding lines from ICARDA along with landraces from Turkey and Spain. In addition, protein content and mineral concentration were determined. Protein content and mineral concentrations for Ca, S, Mg, P, K, Fe, Mn, Cu and Zn varied significantly among accessions. The results show that Kyrgyz germplasm provides a source of diversity for improvement of chickpea.

An integrated genomic approach for rapid delineation of candidate genes regulating agro-morphological traits in chickpea

The identification and fine mapping of robust quantitative trait loci (QTLs)/genes governing important agro-morphological traits in chickpea still lacks systematic efforts at a genome-wide scale involving wild Cicer accessions. In this context, an 834 simple sequence repeat and single-nucleotide polymorphism marker-based high-density genetic linkage map between cultivated and wild parental accessions (Cicer arietinum desi cv. ICC 4958 and Cicer reticulatum wild cv. ICC 17160) was constructed. This inter-specific genetic map comprising eight linkage groups spanned a map length of 949.4 cM with an average inter-marker distance of 1.14 cM. Eleven novel major genomic regions harbouring 15 robust QTLs (15.6–39.8% R² at 4.2–15.7 logarithm of odds) associated with four agro-morphological traits (100-seed weight, pod and branch number/plant and plant hairiness) were identified and mapped on chickpea chromosomes. Most of these QTLs showed positive additive gene effects with effective allelic contribution from ICC 4958, particularly for increasing seed weight (SW) and pod and branch number. One robust SW-influencing major QTL region (qSW4.2) has been narrowed down by combining QTL mapping with high-resolution QTL region-specific association analysis, differential expression profiling and gene haplotype-based association/LD mapping. This enabled to delineate a strong SW-regulating ABI3VP1 transcription factor (TF) gene at trait-specific QTL interval and consequently identified favourable natural allelic variants and superior high seed weight-specific haplotypes in the upstream regulatory region of this gene showing increased transcript expression during seed development. The genes (TFs) harbouring diverse trait-regulating QTLs, once validated and fine-mapped by our developed rapid integrated genomic approach and through gene/QTL map-based cloning, can be utilized as potential candidates for marker-assisted genetic enhancement of chickpea.

Marker Assisted Selection (MAS) for chickpea Fusarium oxysporum wilt resistant genotypes using PCR based molecular markers

The exploration of genetically superior accessions is the key source of germplasm conservation and potential breeding material for the future. To meet the demand of better yielding chickpea cultivars in Pakistan the present study was organized to select more stable and resistant lines from indigenous as well as exotic chickpea germplasm obtained from Plant Genetic Resource Institute (PGRI), National Agricultural Research Centre, Islamabad, Pakistan. For the identification and evaluation of chickpea wilt resistant lines against Fusarium oxysporum f. sp. ciceris (Schlechtends), the germplasm was tested in the field for the selection of wilt resistant lines and the PCR based molecular markers were investigated to use Marker Assisted Selection (MAS) for selection of the desirable cultivars. In field trial, 70 % accessions were resistant to wilt disease, while the remaining 30 % have shown susceptibility to the disease. A total of 5 RAPD and 15 SSR markers were screened for molecular based characterization of wilt response. The data of molecular markers were scored by the presence (1) and absence (0) of allele and subjected to statistical analysis. The analysis was based on coefficient of molecular similarity using UPGMA and sorted the germplasm into two groups based on disease response. Among the total used RAPD/SSR primers, only TA194 SSR marker showed linkage to wilt resistant locus at 85 % probability. The linkage of a marker was reconfirmed by receiver operating characteristic curve. The use of the sorted wilt resistant genotypes through SSR marker TA194 can make available ample prospect in MAS breeding for yield improvement of the crop in Pakistan.

CicArMiSatDB: the chickpea microsatellite database

Background

Chickpea (Cicer arietinum) is a widely grown legume crop in tropical, sub-tropical and temperate regions. Molecular breeding approaches seem to be essential for enhancing crop productivity in chickpea. Until recently, limited numbers of molecular markers were available in the case of chickpea for use in molecular breeding. However, the recent advances in genomics facilitated the development of large scale markers especially SSRs (simple sequence repeats), the markers of choice in any breeding program. Availability of genome sequence very recently opens new avenues for accelerating molecular breeding approaches for chickpea improvement.

Description

In order to assist genetic studies and breeding applications, we have developed a user friendly relational database named the Chickpea Microsatellite Database (CicArMiSatDB http://cicarmisatdb.icrisat.org). This database provides detailed information on SSRs along with their features in the genome. SSRs have been classified and made accessible through an easy-to-use web interface.

Conclusions

This database is expected to help chickpea community in particular and legume community in general, to select SSRs of particular type or from a specific region in the genome to advance both basic genomics research as well as applied aspects of crop improvement.

Mapping and identification of a Cicer arietinum NSP2 gene involved in nodulation pathway

For the first time the putative NSP2 gene in chickpea has been identified using pairs of NILs differing for the Rn1 / rn1 nodulation gene that was located in LG5 of chickpea genetic map. An intraspecific cross between the mutant non-nodulating genotype PM233, carrying the recessive gene rn1, and the wild-type CA2139 was used to develop two pairs of near-isogenic lines (NILs) for nodulation in chickpea. These pairs of NILs were characterized using sequence tagged microsatellite site (STMS) markers distributed across different linkage groups (LGs) of the chickpea genetic map leading to the detection of polymorphic markers located in LG5. Using this information, together with the genome annotation in Medicago truncatula, a candidate gene (NSP2) known to be involved in nodulation pathway was selected for mapping in chickpea. The full length sequence obtained in chickpea wild-type (CaNSP2) was 1,503 bp. Linkage analysis in an F3 population of 118 plants derived from the cross between the pair of NILS NIL7-2A (nod) × NIL7-2B (non-nod) revealed a co-localization between CaNSP2 and Rn1 gene. These data implicate the CaNSP2 gene as a candidate for identity to Rn1, and suggest that it could act in the nodulation signaling transduction pathway similarly to that in other legumes species.

Genetic dissection of drought tolerance in chickpea (Cicer arietinum L.)

Analysis of phenotypic data for 20 drought tolerance traits in 1-7 seasons at 1-5 locations together with genetic mapping data for two mapping populations provided 9 QTL clusters of which one present on CaLG04 has a high potential to enhance drought tolerance in chickpea improvement. Chickpea (Cicer arietinum L.) is the second most important grain legume cultivated by resource poor farmers in the arid and semi-arid regions of the world. Drought is one of the major constraints leading up to 50% production losses in chickpea. In order to dissect the complex nature of drought tolerance and to use genomics tools for enhancing yield of chickpea under drought conditions, two mapping populations-ICCRIL03 (ICC 4958 × ICC 1882) and ICCRIL04 (ICC 283 × ICC 8261) segregating for drought tolerance-related root traits were phenotyped for a total of 20 drought component traits in 1-7 seasons at 1-5 locations in India. Individual genetic maps comprising 241 loci and 168 loci for ICCRIL03 and ICCRIL04, respectively, and a consensus genetic map comprising 352 loci were constructed ( http://cmap.icrisat.ac.in/cmap/sm/cp/varshney/). Analysis of extensive genotypic and precise phenotypic data revealed 45 robust main-effect QTLs (M-QTLs) explaining up to 58.20% phenotypic variation and 973 epistatic QTLs (E-QTLs) explaining up to 92.19% phenotypic variation for several target traits. Nine QTL clusters containing QTLs for several drought tolerance traits have been identified that can be targeted for molecular breeding. Among these clusters, one cluster harboring 48% robust M-QTLs for 12 traits and explaining about 58.20% phenotypic variation present on CaLG04 has been referred as "QTL-hotspot". This genomic region contains seven SSR markers (ICCM0249, NCPGR127, TAA170, NCPGR21, TR11, GA24 and STMS11). Introgression of this region into elite cultivars is expected to enhance drought tolerance in chickpea.

Functionally relevant microsatellite markers from chickpea transcription factor genes for efficient genotyping applications and trait association mapping

We developed 1108 transcription factor gene-derived microsatellite (TFGMS) and 161 transcription factor functional domain-associated microsatellite (TFFDMS) markers from 707 TFs of chickpea. The robust amplification efficiency (96.5%) and high intra-specific polymorphic potential (34%) detected by markers suggest their immense utilities in efficient large-scale genotyping applications, including construction of both physical and functional transcript maps and understanding population structure. Candidate gene-based association analysis revealed strong genetic association of TFFDMS markers with three major seed and pod traits. Further, TFGMS markers in the 5' untranslated regions of TF genes showing differential expression during seed development had higher trait association potential. The significance of TFFDMS markers was demonstrated by correlating their allelic variation with amino acid sequence expansion/contraction in the functional domain and alteration of secondary protein structure encoded by genes. The seed weight-associated markers were validated through traditional bi-parental genetic mapping. The determination of gene-specific linkage disequilibrium (LD) patterns in desi and kabuli based on single nucleotide polymorphism-microsatellite marker haplotypes revealed extended LD decay, enhanced LD resolution and trait association potential of genes. The evolutionary history of a strong seed-size/weight-associated TF based on natural variation and haplotype sharing among desi, kabuli and wild unravelled useful information having implication for seed-size trait evolution during chickpea domestication.

Molecular mapping and characterization of genes governing time to flowering, seed weight, and plant height in an intraspecific genetic linkage map of chickpea (Cicer arietinum)

Drought is the major constraint to chickpea productivity worldwide. Utilizing early flowering genotypes and larger seed size have been suggested as strategies for breeding in drought zones. Therefore, this study aimed to identify potential markers linked to days-to-flowering, 100-seed weight, and plant height in a chickpea intraspecific F(2:3) population derived from the cross ILC3279 × ICCV2. A closely linked marker (TA117) on linkage group LG3 was identified for the days-to-flowering trait, explaining 33% of the variation. In relation to plant height, a quantitative trait loci (QTL) was located in LG3, close to the Ts5 marker, that explained 29% of phenotypic variation. A QTL for 100-seed weight located in LG4, close to TA176, explained 51% of variation. The identification of a locus linked both to high 100-seed weight and days-to-flowering may account for the correlation observed between these traits in this and other breeding attempts.

An intra-specific consensus genetic map of pigeonpea [Cajanus cajan (L.) Millspaugh] derived from six mapping populations

Pigeonpea (Cajanus cajan L.) is an important food legume crop of rainfed agriculture. Owing to exposure of the crop to a number of biotic and abiotic stresses, the crop productivity has remained stagnant for almost last five decades at ca. 750 kg/ha. The availability of a cytoplasmic male sterility (CMS) system has facilitated the development and release of hybrids which are expected to enhance the productivity of pigeonpea. Recent advances in genomics and molecular breeding such as marker-assisted selection (MAS) offer the possibility to accelerate hybrid breeding. Molecular markers and genetic maps are pre-requisites for deploying MAS in breeding. However, in the case of pigeonpea, only one inter- and two intra-specific genetic maps are available so far. Here, four new intra-specific genetic maps comprising 59-140 simple sequence repeat (SSR) loci with map lengths ranging from 586.9 to 881.6 cM have been constructed. Using these four genetic maps together with two recently published intra-specific genetic maps, a consensus map was constructed, comprising of 339 SSR loci spanning a distance of 1,059 cM. Furthermore, quantitative trait loci (QTL) analysis for fertility restoration (Rf) conducted in three mapping populations identified four major QTLs explaining phenotypic variances up to 24 %. To the best of our knowledge, this is the first report on construction of a consensus genetic map in pigeonpea and on the identification of QTLs for fertility restoration. The developed consensus genetic map should serve as a reference for developing new genetic maps as well as correlating with the physical map in pigeonpea to be developed in near future. The availability of more informative markers in the bins harbouring QTLs for sterility mosaic disease (SMD) and Rf will facilitate the selection of the most suitable markers for genetic analysis and molecular breeding applications in pigeonpea.

Development of SSR markers and construction of a linkage map in jute

Jute is an important natural fibre crop, which is only second to cotton in its importance at the global level. It is mostly grown in Indian subcontinent and has been recently used for the development of genomics resources.We recently initiated a programme to develop simple sequence repeat markers and reported a set of 2469 SSR that were developed using four SSR-enriched libraries (Mir et al. 2009). In this communication, we report an additional set of 607 novel SSR in 393 SSR containing sequences. However, primers could be designed for only 417 potentially useful SSR. Polymorphism survey was carried out for 374 primer pairs using two parental genotypes (JRO 524 and PPO4) of a mapping population developed for fibre fineness; only 66 SSR were polymorphic. Owing to a low level of polymorphism between the parental genotypes and a high degree of segregation distortion in recombinant inbred lines, genotypic data of only 53 polymorphic SSR on the mapping population consisting of 120 RIL could be used for the construction of a linkage map; 36 SSR loci were mapped on six linkage groups that covered a total genetic distance of 784.3 cM. Hopefully, this map will be enriched with more SSR loci in future and will prove useful for identification of quantitative trait loci/genes for molecular breeding involving improvement of fibre fineness and other related traits in jute.

High-throughput SNP discovery and genotyping for constructing a saturated linkage map of chickpea (Cicer arietinum L.)

The present study reports the large-scale discovery of genome-wide single-nucleotide polymorphisms (SNPs) in chickpea, identified mainly through the next generation sequencing of two genotypes, i.e. Cicer arietinum ICC4958 and its wild progenitor C. reticulatum PI489777, parents of an inter-specific reference mapping population of chickpea. Development and validation of a high-throughput SNP genotyping assay based on Illumina's GoldenGate Genotyping Technology and its application in building a high-resolution genetic linkage map of chickpea is described for the first time. In this study, 1022 SNPs were identified, of which 768 high-confidence SNPs were selected for designing the custom Oligo Pool All (CpOPA-I) for genotyping. Of these, 697 SNPs could be successfully used for genotyping, demonstrating a high success rate of 90.75%. Genotyping data of the 697 SNPs were compiled along with those of 368 co-dominant markers mapped in an earlier study, and a saturated genetic linkage map of chickpea was constructed. One thousand and sixty-three markers were mapped onto eight linkage groups spanning 1808.7 cM (centiMorgans) with an average inter-marker distance of 1.70 cM, thereby representing one of the most advanced maps of chickpea. The map was used for the synteny analysis of chickpea, which revealed a higher degree of synteny with the phylogenetically close Medicago than with soybean. The first set of validated SNPs and map resources developed in this study will not only facilitate QTL mapping, genome-wide association analysis and comparative mapping in legumes but also help anchor scaffolds arising out of the whole-genome sequencing of chickpea.

Transcriptome sequencing of wild chickpea as a rich resource for marker development

The transcriptome of cultivated chickpea (Cicer arietinum L.), an important crop legume, has recently been sequenced. Here, we report sequencing of the transcriptome of wild chickpea, C. reticulatum (PI489777), the progenitor of cultivated chickpea, by GS-FLX 454 technology. The optimized assembly of C. reticulatum transcriptome generated 37 265 transcripts in total with an average length of 946 bp. A total of 4072 simple sequence repeats (SSRs) could be identified in these transcript sequences, of which at least 561 SSRs were polymorphic between C. arietinum and C. reticulatum. In addition, a total of 36 446 single-nucleotide polymorphisms (SNPs) were identified after optimization of probability score, quality score, read depth and consensus base ratio. Several of these SSRs and SNPs could be associated with tissue-specific and transcription factor encoding transcripts. A high proportion (92-94%) of polymorphic SSRs and SNPs identified between the two chickpea species were validated successfully. Further, the estimation of synonymous substitution rates of orthologous transcript pairs suggested that the speciation event for divergence of C. arietinum and C. reticulatum may have happened approximately 0.53 million years ago. The results of our study provide a rich resource for exploiting genetic variations in chickpea for breeding programmes.

Large-scale development of cost-effective SNP marker assays for diversity assessment and genetic mapping in chickpea and comparative mapping in legumes

A set of 2486 single nucleotide polymorphisms (SNPs) were compiled in chickpea using four approaches, namely (i) Solexa/Illumina sequencing (1409), (ii) amplicon sequencing of tentative orthologous genes (TOGs) (604), (iii) mining of expressed sequence tags (ESTs) (286) and (iv) sequencing of candidate genes (187). Conversion of these SNPs to the cost-effective and flexible throughput Competitive Allele Specific PCR (KASPar) assays generated successful assays for 2005 SNPs. These marker assays have been designated as Chickpea KASPar Assay Markers (CKAMs). Screening of 70 genotypes including 58 diverse chickpea accessions and 12 BC(3) F(2) lines showed 1341 CKAMs as being polymorphic. Genetic analysis of these data clustered chickpea accessions based on geographical origin. Genotyping data generated for 671 CKAMs on the reference mapping population (Cicer arietinum ICC 4958 × Cicer reticulatum PI 489777) were compiled with 317 unpublished TOG-SNPs and 396 published markers for developing the genetic map. As a result, a second-generation genetic map comprising 1328 marker loci including novel 625 CKAMs, 314 TOG-SNPs and 389 published marker loci with an average inter-marker distance of 0.59 cM was constructed. Detailed analyses of 1064 mapped loci of this second-generation chickpea genetic map showed a higher degree of synteny with genome of Medicago truncatula, followed by Glycine max, Lotus japonicus and least with Vigna unguiculata. Development of these cost-effective CKAMs for SNP genotyping will be useful not only for genetics research and breeding applications in chickpea, but also for utilizing genome information from other sequenced or model legumes.

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.

Phenotyping chickpeas and pigeonpeas for adaptation to drought

The chickpea and pigeonpea are protein-rich grain legumes used for human consumption in many countries. Grain yield of these crops is low to moderate in the semi-arid tropics with large variation due to high GxE interaction. In the Indian subcontinent chickpea is grown in the post-rainy winter season on receding soil moisture, and in other countries during the cool and dry post winter or spring seasons. The pigeonpea is sown during rainy season which flowers and matures in post-rainy season. The rainy months are hot and humid with diurnal temperature varying between 25 and 35°C (maximum) and 20 and 25°C (minimum) with an erratic rainfall. The available soil water during post-rainy season is about 200-250 mm which is bare minimum to meet the normal evapotranspiration. Thus occurrence of drought is frequent and at varying degrees. To enhance productivity of these crops cultivars tolerant to drought need to be developed. ICRISAT conserves a large number of accessions of chickpea (>20,000) and pigeonpea (>15,000). However only a small proportion (<1%) has been used in crop improvement programs mainly due to non-availability of reliable information on traits of economic importance. To overcome this, core and mini core collections (10% of core, 1% of entire collection) have been developed. Using the mini core approach, trait-specific donor lines were identified for agronomic, quality, and stress related traits in both crops. Composite collections were developed both in chickpea (3000 accessions) and pigeonpea (1000 accessions), genotyped using SSR markers and genotype based reference sets of 300 accessions selected for each crop. Screening methods for different drought-tolerant traits such as early maturity (drought escape), large and deep root system, high water-use efficiency, smaller leaflets, reduced canopy temperature, carbon isotope discrimination, high leaf chlorophyll content (drought avoidance), and breeding strategies for improving drought tolerance have been discussed.

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.

Novel SSR markers from BAC-end sequences, DArT arrays and a comprehensive genetic map with 1,291 marker loci for chickpea (Cicer arietinum L.)

Chickpea (Cicer arietinum L.) is the third most important cool season food legume, cultivated in arid and semi-arid regions of the world. The goal of this study was to develop novel molecular markers such as microsatellite or simple sequence repeat (SSR) markers from bacterial artificial chromosome (BAC)-end sequences (BESs) and diversity arrays technology (DArT) markers, and to construct a high-density genetic map based on recombinant inbred line (RIL) population ICC 4958 (C. arietinum)×PI 489777 (C. reticulatum). A BAC-library comprising 55,680 clones was constructed and 46,270 BESs were generated. Mining of these BESs provided 6,845 SSRs, and primer pairs were designed for 1,344 SSRs. In parallel, DArT arrays with ca. 15,000 clones were developed, and 5,397 clones were found polymorphic among 94 genotypes tested. Screening of newly developed BES-SSR markers and DArT arrays on the parental genotypes of the RIL mapping population showed polymorphism with 253 BES-SSR markers and 675 DArT markers. Segregation data obtained for these polymorphic markers and 494 markers data compiled from published reports or collaborators were used for constructing the genetic map. As a result, a comprehensive genetic map comprising 1,291 markers on eight linkage groups (LGs) spanning a total of 845.56 cM distance was developed (http://cmap.icrisat.ac.in/cmap/sm/cp/thudi/). The number of markers per linkage group ranged from 68 (LG 8) to 218 (LG 3) with an average inter-marker distance of 0.65 cM. While the developed resource of molecular markers will be useful for genetic diversity, genetic mapping and molecular breeding applications, the comprehensive genetic map with integrated BES-SSR markers will facilitate its anchoring to the physical map (under construction) to accelerate map-based cloning of genes in chickpea and comparative genome evolution studies in legumes.

Large-scale transcriptome analysis in chickpea (Cicer arietinum L.), an orphan legume crop of the semi-arid tropics of Asia and Africa

Chickpea (Cicer arietinum L.) is an important legume crop in the semi-arid regions of Asia and Africa. Gains in crop productivity have been low however, particularly because of biotic and abiotic stresses. To help enhance crop productivity using molecular breeding techniques, next generation sequencing technologies such as Roche/454 and Illumina/Solexa were used to determine the sequence of most gene transcripts and to identify drought-responsive genes and gene-based molecular markers. A total of 103,215 tentative unique sequences (TUSs) have been produced from 435,018 Roche/454 reads and 21,491 Sanger expressed sequence tags (ESTs). Putative functions were determined for 49,437 (47.8%) of the TUSs, and gene ontology assignments were determined for 20,634 (41.7%) of the TUSs. Comparison of the chickpea TUSs with the Medicago truncatula genome assembly (Mt 3.5.1 build) resulted in 42,141 aligned TUSs with putative gene structures (including 39,281 predicted intron/splice junctions). Alignment of ∼37 million Illumina/Solexa tags generated from drought-challenged root tissues of two chickpea genotypes against the TUSs identified 44,639 differentially expressed TUSs. The TUSs were also used to identify a diverse set of markers, including 728 simple sequence repeats (SSRs), 495 single nucleotide polymorphisms (SNPs), 387 conserved orthologous sequence (COS) markers, and 2088 intron-spanning region (ISR) markers. This resource will be useful for basic and applied research for genome analysis and crop improvement in chickpea.

Integration of genetic and physical maps of the chickpea (Cicer arietinum L.) genome using flow-sorted chromosomes

Cultivated chickpea is the third most important legume after field bean and garden pea worldwide. Despite considerable breeding towards improved yield and resistance to biotic and abiotic stresses, the production of chickpea remained stagnant, but molecular tools are expected to increase the impact of current improvement programs. As a first step towards this goal, various genetic linkage maps have been established and markers linked to resistance genes been identified. However, until now, only one linkage group (LG) has been assigned to a specific chromosome. In the present work, mitotic chromosomes were sorted using flow cytometry and used as template for PCR with primers designed for genomic regions flanking microsatellites. These primers amplify sequence-tagged microsatellite site markers. This approach confirmed the assignment of LG8 to the smallest chromosome H. For the first time, LG5 was linked to the largest chromosome A, LG4 to a medium-sized chromosome E, while LG3 was anchored to the second largest chromosome B. Chromosomes C and D could not be flow-sorted separately and were jointly associated to LG6 and LG7. By the same token, chromosomes F and G were anchored to LG1 and LG2. To establish a set of preferably diagnostic cytogenetic markers, the genomic distribution of various probes was verified using FISH. Moreover, a partial genomic bacterial artificial chromosome (BAC) library was constructed and putative single/low-copy BAC clones were mapped cytogenetically. As a result, two clones were identified localizing specifically to chromosomes E and H, for which no cytogenetic markers were yet available.

Molecular mapping for resistance to pea rust caused by Uromyces fabae (Pers.) de-Bary

Pea rust caused by Uromyces fabae (Pers.) de-Bary is a major problem in warm humid regions causing huge economic losses. A mapping population of 136 F(6:7) recombinant inbred lines (RILs) derived from the cross between pea genotypes, HUVP 1 (susceptible) and FC 1 (resistant) was evaluated in polyhouse as well as under field conditions during two consecutive years. Infection frequency (IF) and area under disease progress curve (AUDPC) were used for evaluation of rust reaction of the RILs. A linkage map was constructed with 57 polymorphic loci selected from 148 simple sequence repeats (SSRs), 3 sequence tagged sites (STS), and 2 random amplified polymorphic (RAPD) markers covering 634 cM of genetic distance on the seven linkage groups of pea with an average interval length of 11.3 cM. Composite interval mapping (CIM) revealed one major (Qruf) and one minor (Qruf1) QTL for rust resistance on LGVII. The LOD (5.2-15.8) peak for Qruf was flanked by SSR markers, AA505 and AA446 (10.8 cM), explaining 22.2-42.4% and 23.5-58.8% of the total phenotypic variation for IF and AUDPC, respectively. The minor QTL was environment-specific, and it was detected only in the polyhouse (LOD values 4.2 and 4.8). It was flanked by SSR markers, AD146 and AA416 (7.3 cM), and explained 11.2-12.4% of the total phenotypic variation. The major QTL Qruf was consistently identified across all the four environments. Therefore, the SSR markers flanking Qruf would be useful for marker-assisted selection for pea rust (U. fabae) resistance.

Development and use of genic molecular markers (GMMs) for construction of a transcript map of chickpea (Cicer arietinum L.)

A transcript map has been constructed by the development and integration of genic molecular markers (GMMs) including single nucleotide polymorphism (SNP), genic microsatellite or simple sequence repeat (SSR) and intron spanning region (ISR)-based markers, on an inter-specific mapping population of chickpea, the third food legume crop of the world and the first food legume crop of India. For SNP discovery through allele re-sequencing, primer pairs were designed for 688 genes/expressed sequence tags (ESTs) of chickpea and 657 genes/ESTs of closely related species of chickpea. High-quality sequence data obtained for 220 candidate genic regions on 2-20 genotypes representing 9 Cicer species provided 1,893 SNPs with an average frequency of 1/35.83 bp and 0.34 PIC (polymorphism information content) value. On an average 2.9 haplotypes were present in 220 candidate genic regions with an average haplotype diversity of 0.6326. SNP2CAPS analysis of 220 sequence alignments, as mentioned above, provided a total of 192 CAPS candidates. Experimental analysis of these 192 CAPS candidates together with 87 CAPS candidates identified earlier through in silico mining of ESTs provided scorable amplification in 173 (62.01%) cases of which predicted assays were validated in 143 (82.66%) cases (CGMM). Alignments of chickpea unigenes with Medicago truncatula genome were used to develop 121 intron spanning region (CISR) markers of which 87 yielded scorable products. In addition, optimization of 77 EST-derived SSR (ICCeM) markers provided 51 scorable markers. Screening of easily assayable 281 markers including 143 CGMMs, 87 CISRs and 51 ICCeMs on 5 parental genotypes of three mapping populations identified 104 polymorphic markers including 90 markers on the inter-specific mapping population. Sixty-two of these GMMs together with 218 earlier published markers (including 64 GMM loci) and 20 other unpublished markers could be integrated into this genetic map. A genetic map developed here, therefore, has a total of 300 loci including 126 GMM loci and spans 766.56 cM, with an average inter-marker distance of 2.55 cM. In summary, this is the first report on the development of large-scale genic markers including development of easily assayable markers and a transcript map of chickpea. These resources should be useful not only for genome analysis and genetics and breeding applications of chickpea, but also for comparative legume genomics.

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.