Chief EditorJ. S. Sandhu
Print ISSN 0250-5371
Online ISSN 0976-0571
NAAS Rating 6.67
Impact Factor 0.8 (2023)
Full Research Article
Prioritization of Microsatellite Markers Linked with Drought Tolerance Associated Traits in Chickpea (Cicer arietinum L.)
Background: Being a vital source of high-quality dietary protein, chickpea is an unavoidable legume. The present investigation was performed to study the applicability of the microsatellite markers linked with drought tolerance in Indian chickpea genotypes collected from different genetic background.
Methods: In Rabi 2021-22, forty chickpea genotypes including national check for drought tolerance, elite cultivars, released varieties and advanced breeding lines were screened employing forty microsatellite markers linked with drought tolerance associated traits.
Results: Among forty drought tolerance related microsatellite markers, twenty-six were found to be polymorphic and produced a total of 66 alleles, with a mean of 2.5 alleles per locus. Genotype NCPGR127 had highest numbers of alleles, gene diversity and polymorphism information content tracked by genotype TR19. These genotypes may be used for selecting drought tolerant genotypes through marker assistant selection (MAS). Model-based population structure analysis clearly distinguished the drought tolerant genotypes including ICC4958, JG74, JAKI9218, JG16, JG6, JG14 and JG11. Principle coordinate analysis (PCoA) and analysis of molecular variance (AMOVA) further confirmed these results. Findings of the present investigation have a greater potential for further utilization in breeding of drought specific chickpea cultivar(s).
Abiotic stresses, especially drought and heat reduce yield of different crops (Mishra et al., 2021; Sharma et al., 2021; Asati et al., 2022) including chickpea over 70% (Varshney et al., 2019). The application of molecular markers may shorten the crop breeding cycle and ultimately leads to crop improvement particularly in complex traits controlled by polygenes (Mishra et al., 2022; Solanki et al., 2022; Yadav et al., 2023). Several quantitative trait loci (QTLs) for drought-tolerance associated traits have been reported in chickpea. The QTL-hotspot genomic region present on CaLG04 explained >50% phenotypic variation for drought tolerance (Varshney et al., 2014). Additionally, genome-wide microsatellite markers associated with drought tolerance related traits were also reported (Li et al., 2018; Varshney et al., 2019).
Marker-assisted backcross breeding has been successfully employed for developing superior chickpea lines with improved yield and drought tolerance. Drought tolerant variety Pusa Chickpea 10216 has also been released for commercial cultivation in India (Bharadwaj et al., 2021). In the light of above facts, employing drought tolerance linked microsatellite markers, the present study was aimed to assess the pattern and level of genetic diversity and developing population structure among national check for drought tolerance, elite cultivars, released varieties and advanced breeding lines.
MATERIALS AND METHODS
Forty chickpea genotypes were obtained from the All India Coordinated Research Project (AICRP) on Chickpea, Department of Plant Breeding and Genetics, Jawaharlal Nehru Krishi Vishwa Vidyalaya, Jabalpur (Table 1). The genotypes were grown in the Net house of Biotechnology Centre, Jawaharlal Nehru Krishi Vishwa Vidyalaya, Jabalpur during Rabi 2021-22.
DNA extraction and quantification
Genomic DNA from young leaves of grown chickpea genotypes was isolated using the cetyl triethyl ammonium bromide (CTAB) method (Doyle and Doyle 1987). Two grams of leaf samples were ground in liquid nitrogen into a fine powder using sterilized pestle and mortar. About 100 mg of leaf powder was transferred to 2 ml micro-centrifuge tubes containing 1 ml of extraction buffer (0.1M Tris–HCl (pH 7.6), 0.005M EDTA and 0.2M Tris-HCl, 0.05M EDTA, 2M NaCl and 2% CTAB). Sample tubes were incubated at 65°C for one hour then centrifuged at 10000 rpm for 15 minutes with chloroform-isoamyl alcohol (24:1) mixture for protein separation. Supernatant was transferred to fresh tube and chilled isopropanol was used for the precipitation of DNA. Precipitated DNA was pelleted down through centrifugation and the DNA pellet was dried and dissolved in 50 μl of Milli-Q (MQ) water. The quantity and quality of all DNA samples were checked using agarose gel (0.8%). The working DNA sample was obtained by diluting to a final DNA concentration of 40-50 ng per microlitre (μL).
Polymerase chain reaction (PCR) and gel electrophoresis
Forty genome wide distributed drought tolerance linked microsatellite markers were selected based on polymorphic information content (PIC), allelic richness and heterozygosity from previous studies (Table 2). Primers were synthesized from Integrated DNA Technologies (Asia Pacific, Singapore). The PCR reactions were performed with a total volume of 10ìl containing 40-50 ng of DNA, 1.5 mM MgCl2, 0.2 mM dNTPs, 0.4 mM each of the forward and reverse primers and 1U Taq DNA polymerase in a Sure Cycler 8800 (Agilent Technologies) after optimizing the annealing conditions specific for each primer pair. The PCR was performed with an initial denaturation step of 5 min at 94°C, tracked by 35 cycles of 30s denaturation at 94°C, annealing at 55-60°C (depending on the primer) for 50s, initial extension at 72°C for 50s. The final extension was accomplished at 72°C for seven min. The PCR products were resolved on 2.8% agarose gel in 1 × TAE buffer with a 6 × DNA loading dye. Electrophoresis was carried out on horizontal electrophoresis with 65V for two and half hours using a standard 100 bp DNA ladder (Promega, Madison, USA). Gel pictures were snapped on a Vilber Laurmet, QUANTUM-ST5 gel electrophoresis system.
Analysis of genotypic data through a model-based approach
Scoring each polymorphic SSR marker was done with the help of Gel Analyzer V19.1 software. The allelic data scores were used for computing locus-based diversity indices with the help of GenAlEx V6.51b2 software and Power Marker V3.25 software. The Population structure for 40 genotypes was constructed using STRUCTURE version 2.3.4 software. An analysis of molecular variance (AMOVA) and Principal Coordinates Analysis (PCoA) within and among populations were performed using GenAlEx V6.51b2.
RESULTS AND DISCUSSION
Polymorphism among chickpea genotypes
The SSRs have been widely utilized for genetic diversity analysis, germplasm grouping and population structure analysis in numerous crops, including chickpea (Choudhary et al., 2012; Mishra et al., 2020; Rathore et al., 2022). The present investigation determined the suitability of selected trait linked microsatellite markers for genetic diversity analysis of chickpea genotypes. Out of 40 applied markers, only 26 showed consistent polymorphic banding patterns with genotypes, indicating their suitability for genetic diversity analysis. The remaining fourteen that produced monomorphic bands among genotypes revealing one allele at each locus were not considered for diversity analysis.
Overall, 66 alleles were identified by 26 polymorphic markers from the 40 tested genotypes with a mean of 2.5 alleles per locus (Table 3). The allelic richness (Na) per locus ranged between two (most of the markers) to four (TR 19 and NCPGR 127). The numbers of effective alleles (Ne) ranged from 1.22 (ICCeM058) to 3.57 (NCPGR 127), with an overall mean of 1.92 effective alleles per locus. The most frequent major allele frequency (0.90) was found for the marker ICCeM0058 and the lowest (0.40) for marker NCPGR127, with an average of 0.66. The highest gene diversity (0.72) and PIC value (0.67) were found for the marker NCPGR127, whilst the lowest gene diversity (0.18) and PIC value (0.16) was investigated for the marker ICCeM0058, with an average of 0.44 and 0.37, respectively. The microsatellite marker analysis result indicated moderate allelic richness per locus and relatively moderate to high PIC, Ho and He values. A high level of genetic diversity indicated the existence of higher molecular variation among the 40 tested chickpea genotypes, in agreement with the previous studies (Mir et al., 2022). During present study, the most of the loci (20) exhibited a low level of observed heterozygosity in comparison to the expected heterozygosity and high fixation index (Choudhary et al., 2012). It indicated high levels of expected inbreeding among tested chickpea genotypes because chickpea is a self-pollinated crop and only 0 to 1.58% of outcrossing is reported (Ghaffari et al., 2014). The rest of the six loci had a high level of observed heterozygosity compared to the expected heterozygosity with a low associated fixation index. It implied that higher mutation rates or inbreeding depression could be associated with these loci (Choudhary et al., 2012).
Model-based clustering method of the STRUCTURE software is valuable in showing the presence of population structure in tested genotypes, identifying different genetic populations, assigning individuals to populations and identifying pure and admixed individuals (Varshney et al., 2019). A population structure model without admixture was investigated using all 40 genotypes and 26 polymorphic markers for maximum likelihood by correlating allele frequency with a 5000-burn period with a run length of 50000 to 10 Markov Chain Monte Carlo (MCMC) with varying K from 1 to 15 with five iterations. The highest peak was observed at K = 4 with 30.78 Delta K value in Delta K analysis, indicating the presence of four populations in tested genotypes (Fig 1). At K=4, all genotypes were stratified into four populations viz., P1, P2, P3 and P4, representing 22., 22.5, 37.5 and 17.5% of genotypes used in structure analysis, respectively (Fig 2A, presents the inferred population structure). Population group 2 had clearly distingued the drought tolerant genotypes viz., ICC4958, JG74, JAKI9218, JG16, JG6, JG14 and JG11 from other studied chickpea genotypes. These distinguished drought tolerant chickpea genotypes had been already established as drought tolerant chickpea genotypes. The ICC4958 is a national check for drought tolerance in India, JG74 is an early maturing chickpea cultivar and established as most suitable chickpea cultivar for rainfed conditions of central region of India (Mannur et al., 2019), whilst JAKI9218 and JG16 are established as drought tolerant chickpea varieties. JG14 is identified as a heat tolerant chickpea variety (Dixit et al., 2019) and JG11 is established as a drought tolerant chickpea variety (Dixit, 2021).
In principal coordinate analysis (PCA) the first two Eigen vectors classified the studied genotypes into four major groups and distinguished the drought tolerant chickpea genotypes from investigated genotypes, which were comparable to observation obtained from STRUCTURE analysis (Fig 2B). The genotypes were divided into four structured populations, revealing that tested chickpea genotypes evolved from four population types demonstrating varying degrees of introgression from four types of ancestors into tested genotypes. It might be due to the resulting genotypes from independent evolutionary mechanisms including genetic drift, migration, mutation, selection and germplasm exchange that split them into discrete gene pools.
Analysis of molecular variance (AMOVA)
The four population groups generated from the structural analysis were subjected to AMOVA to estimate the percentage of genetic variation within as well as among individuals and populations (Table 4). The total genetic variance of genotypes was contributed by 7% variation within individuals, 77% among individuals and 16% among populations (Fig 3). Pairwise inbreeding coefficient within subpopulations (Fst) value (0.155) demonstrated significant variation among all the four pairs of populations, suggesting that all four groups were different from each other. The results acquired from AMOVA analysis agreed with the observations found through the basic summary statistics, structure analysis and principal coordinate analysis, confirming both the presence of a statistically moderate genetic diversity among genotypes and a high level of population structure among genotypes.
The AMOVA results demonstrated that most variation among genotypes might be credited to variation among individual genotypes compared to variation among and within the population. It suggested that individual variation of genotypes might be more valuable for chickpea improvement programs. The PCoA result depicted that the chickpea genotypes were equally spread in all four quadrants regardless of their population structure. Population one was exclusively grouped in quadrant I, showing their distinctness from the rest of the populations evenly spread in all four quadrants because of their highly diverse nature. Similar results have also been reported in various earlier studies using microsatellite markers (Mir et al., 2022).
CONFLICT OF INTEREST
- Asati, R., Tripathi, M.K., Tiwari, S., Yadav, R.K., Tripathi, N. (2022). Molecular breeding and drought tolerance in chickpea. Life. 12(11): 1846. https://doi.org/10.3390/life12111846.
- Bharadwaj, C., Tripathi, S., Soren, K.R., Thudi, M., Singh, R.K., Sheoran, S. et al. (2021). Introgression of “QTLhotspot” region enhances drought tolerance and grain yield in three elite chickpea cultivars. Plant Genome. 14(1): e20076. https://doi.org/10.1002/tpg2.20076.
- Choudhary, P., Khanna, S.M., Jain, P.K. (2012). Genetic structure and diversity analysis of the primary gene pool of chickpea using SSR markers. Genetics of Molecular Research. 11: 891-905.
- Dixit, G.P. (2021). Sustaining Chickpea growth in India: Breeder’s Perspective. Journal of Food Legumes. 34:73-75.
- Dixit, G.P., Srivastava, A.K., Singh, N.P. (2019). Marching towards self-sufficiency in chickpea. Current Science. 116: 239- 242.
- Doyle, J.J., Doyle, J.L. (1987). A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochemistry Bulletin. 19: 11-15.
- Ghaffari, P., Talebi, R., Keshavarzi, F. (2014). Genetic diversity and geographical differentiation of Iranian landrace, cultivars and exotic chickpea lines as revealed by morphological and microsatellite markers. Physiology and Molecular Biology of Plants. 20: 225-233.
- Gupta, N., Tiwari, S., Tripathi, M.K., Bhagyawant, S.S. (2021). Antinutritional and protein-based profiling of diverse desi and wild chickpea accessions. Current Journal of Applied Science and Technology. 40(6): 7-18.
- Li, Y., Ruperao, P., Batley, J., Edwards, D., Khan, T., Colmer, T.D., Sutton, T. (2018). Investigating drought tolerance in chickpea using genome-wide association mapping and genomic selection based on whole-genome resequencing data. Frontiers in Plant Science. 9: 190. doi: 10.3389/ fpls.2018.00190.
- Mannur, D.M., Babbar, A., Thudi, M., Sabbavarapu, M.M., Roorkiwal, M. et al. (2019). Super Annigeri 1 and improved JG 74: Two Fusarium wilt-resistant introgression lines developed using marker-assisted backcrossing approach in chickpea (Cicer arietinum L.). Molecular Breeding. 39: 1-13.
- Mir, A.H., Bhat, M.A., Fayaz, H., Wani, A.A., Dar, S.A., Maqbool, S., Mir, R.R. (2022). SSR markers in revealing extent of genetic diversity and phylogenetic relationships among chickpea core collection accessions for Western Himalayas. Molecular Biology Reports. 49: 11469-11479.
- Mishra, N., Tripathi, M.K., Tiwari, S., Tripathi, N., Trivedi, H.K. (2020). Morphological and molecular screening of soybean genotypes against yellow mosaic virus disease. Legume Research. 45(10): 1309-1316. DOI: 10.18805/LR-4240.
- Mishra, N., Tripathi, M.K., Tiwari, S., Tripathi, N., Gupta, N., Sharma, A. (2021). Morphological and physiological performance of Indian soybean [Glycine max (L.) Merrill] genotypes in respect to drought. Legume Research. 1-9. DOI: 10.18805/ LR-4550.
- Mishra, N., Tripathi, M.K., Tripathi, N., Tiwari, S., Gupta, N., Sharma, A. (2022). Screening of soybean genotypes against drought on the basis of gene-linked microsatellite markers. In Book: Innovations in Science and Technology. 3: 49-61.
- Ningwal, R., Tripathi, M.K., Tiwari, S., Asati, R., Yadav, R.K., Tripathi, N., Yasin, M. (2023). Identification of Polymorphic ssr marken and diversity analysis in a set of desi chickpea genotypes. Biological Forum-An International Journal. 15(3): 45-51.
- Rathore, M.S., Tiwari, S., Tripathi, M.K., Gupta, N., Yadav, S., Singh, S., Tomar, R.S. (2022). Genetic diversity analysis of groundnut germplasm lines in respect to early and late leaf spot diseases and biochemical traits. Legume Research. 1-6. DOI: 10.18805/LR-4833.
- Sahu, V.K., Tiwari, S., Gupta, N., Tripathi, M.K., Yasin, M. (2020). Evaluation of physiological and biochemical contents in desi and Kabuli chickpea. Legume Research. 45(10): 1197-1208. DOI:10.18805/LR-4265.
- Solanki, R.S., Babbar, A. and Tripathi, N. (2022). Genetic diversity analysis in kabuli chickpea (Cicer arietinum L.) genotypes based on quantitative traits and molecular markers. Bangladesh Journal of Botany. 51: 581-587.
- Varol, I.S., Kardes, Y.M., Irik, H.A., Kirnak, H., Kaplan, M. (2020). Supplementary irrigations at different physiological growth stages of chickpea (Cicer arietinum L.) change grain nutritional composition. Food Chemistry. 303: 125402. doi: 10.1016/j.foodchem.2019.125402.
- Varshney, R.K., Thudi, M., Nayak, S.N., Gaur, P.M., Kashiwagi, J., Krishnamurthy, L. et al. (2014). Genetic dissection of drought tolerance in chickpea (Cicer arietinum L.). Theoretical and Applied Genetics. 127: 445-462.
- Varshney, R.K., Thudi, M., Roorkiwal, M., He, W., Upadhyaya, H.D., Yang, W. et al. (2019). Resequencing of 429 chickpea accessions from 45 countries provides insights into genome diversity, domestication and agronomic traits. Nature Genetics. 51: 857-864.
- Yadav, R.K., Tripathi, M.K., Tiwari, S., Tripathi, N., Asati, R., Patel, V., Sikarwar, R.S., Payasi, D.K. (2023). Breeding and genomic approaches towards development of fusarium wilt resistance in chickpea. Life. 13: 988. https://doi.org/10.3390/life13040988.
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.
This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.