The current study employed a high-throughput genome-wide next-generation sequencing-led multiple QTL-seq (mQTL-seq) strategy in two inter- and intra-specific recombinant inbred line (RIL) mapping populations to identify the major genomic regions underlying robust quantitative trait loci (QTLs) regulating plant height in chickpea. The whole genome resequencing discovered 446,475 and 150,434 high-quality homozygous single nucleotide polymorphisms (SNPs) exhibiting polymorphism between tall and dwarf/semi-dwarf mapping parents and bulk/homozygous individuals selected from each of two chickpea RIL populations. These SNP-led mQTL-seq assays in RIL mapping populations scaled-down two longer major genomic regions (1.26–1.34 Mb) underlying robust plant height QTLs into the shorter high-resolution QTL intervals (653.2–756.3 kb) on chickpea chromosomes 3 and 8. This essentially delineated regulatory novel natural SNP allelic variants from brassinosteroid insensitive 1-receptor kinase 1 (BAK1) and gibberellin (GA) 20-oxidase genes governing plant height in chickpea. A strong impact of evolutionary bottlenecks including strong artificial/natural selection on two plant height gene loci during chickpea domestication was observed. The shoot apical meristem-specific expression aside from down-regulation of two plant height genes especially in dwarf/semi-dwarf as compared to tall parents and homozygous mapping individuals of two aforementioned RIL populations was apparent. The integrated genomics-assisted breeding strategy combining mQTL-seq with differential gene expression profiling and functional allelic diversity-based trait domestication study collectively identified potential natural allelic variants of candidate genes underlying major plant height QTLs in chickpea. These functionally relevant molecular signatures can be of immense use for marker-aided genetic enhancement to develop high seed- and pod-yielding non-lodging cultivars restructured with desirable plant height in chickpea.
Read full abstract