A Multiple QTL-Seq Strategy Delineates Potential Genomic Loci Governing Flowering Time in Chickpea

Genetic dissection of domestication traits in interspecific chickpea populations

Chickpea (Cicer arietinum) is a pulse crop that provides an integral source of nutrition for human consumption. The close wild relatives Cicer reticulatum and Cicer echinospermum harbor untapped genetic diversity that can be exploited by chickpea breeders to improve domestic varieties. Knowledge of genomic loci that control important chickpea domestication traits will expedite the development of improved chickpea varieties derived from interspecific crosses. Therefore, we set out to identify genomic loci underlying key chickpea domestication traits by both association and quantitative trait locus (QTL) mapping using interspecific F2 populations. Diverse phenotypes were recorded for various agronomic traits. A total of 11 high-confidence markers were detected on chromosomes 1, 3, and 7 by both association and QTL mapping; these were associated with growth habit, flowering time, and seed traits. Furthermore, we identified candidate genes linked to these markers, which advanced our understanding of the genetic basis of domestication traits and validated known genes such as the FLOWERING LOCUS gene cluster that regulates flowering time. Collectively, this study has elucidated the genetic basis of chickpea domestication traits, which can facilitate the development of superior chickpea varieties.

Delineation of genes for a major QTL governing heat stress tolerance in chickpea

Chickpea (Cicer arietinum) is a cool season grain legume experiencing severe yield loss during heat stress due to the intensifying climate changes and its associated gradual increase of mean temperature. Hence, understanding the genetic architecture regulating heat stress tolerance has emerged as an important trait to be addressed for enhancing yield and productivity of chickpea under heat stress. The present study is intended to identify the major genomic region(s) governing heat stress tolerance in chickpea. For this, an integrated genomics-assisted breeding strategy involving NGS-based high-resolution QTL-seq assay, QTL region-specific association analysis and molecular haplotyping was deployed in a population of 206 mapping individuals and a diversity panel of 217 germplasm accessions of chickpea. This combinatorial strategy delineated a major 156.8 kb QTL genomic region, which was subsequently narrowed-down to a functional candidate gene CaHSFA5 and its natural alleles associated strongly with heat stress tolerance in chickpea. Superior natural alleles and haplotypes delineated from the CaHSFA5 gene have functional significance in regulating heat stress tolerance in chickpea. Histochemical staining, interaction studies along with differential expression profiling of CaHSFA5 and ROS scavenging genes suggest a cross talk between CaHSFA5 with ROS homeostasis pertaining to heat stress tolerance in chickpea. Heterologous gene expression followed by heat stress screening further validated the functional significance of CaHSFA5 for heat stress tolerance. The salient outcomes obtained here can have potential to accelerate multiple translational genomic analysis including marker-assisted breeding and gene editing in order to develop high-yielding heat stress tolerant chickpea varieties.

Natural alleles of Mediator subunit genes modulate plant height in chickpea

Plant height (PH) is an important plant architectural trait targeted during Green Revolution to enhance crop yields. Identification of genes and natural alleles governing plant height without compromising agronomic performance can fill the lacuna of knowledge connecting ideal plant architecture with maximum achievable yield in chickpea. Through coherent strategy involving genome‐wide association study, QTL/fine mapping, map‐based cloning, molecular haplotyping, and downstream functional genomics, the current study identified two Mediator subunit genes namely, CaMED23 and CaMED5b and their derived natural alleles/haplotypes underlying the major QTLs and trans‐acting eQTLs regulating plant height in chickpea. Differential accumulation of haplotype‐specific transcripts of these two Mediator genes in corresponding haplotype‐introgressed near‐isogenic lines (NILs) correlates negatively with the plant height trait. Quantitative as well as qualitative estimation based on histology, scanning electron microscopy, and histochemical assay unraveled the reduced lengths and cell sizes of internodes along with compromised lignin levels in dwarf/semi‐dwarf chickpea NILs introgressed with superior CaMED23 and CaMED5b gene haplotypes. This observation, supported by global transcriptome profiling‐based diminished expression of various phenylpropanoid pathway genes upstream of lignin biosynthesis in dwarf/semi‐dwarf NILs, essentially links plant height with lignin accumulation. The identified molecular signatures in the Mediator subunit genes can be efficiently utilized to develop desirable dwarf/semi‐dwarf‐type chickpea cultivars without affecting their yield per plant via modulating lignin/phenylpropanoid biosynthesis.

A superior gene allele involved in abscisic acid signaling enhances drought tolerance and yield in chickpea

Identifying potential molecular tags for drought tolerance is essential for achieving higher crop productivity under drought stress. We employed an integrated genomics-assisted breeding and functional genomics strategy involving association mapping, fine mapping, map-based cloning, molecular haplotyping and transcript profiling in the introgression lines (ILs)- and near isogenic lines (NILs)-based association panel and mapping population of chickpea (Cicer arietinum). This combinatorial approach delineated a bHLH (basic helix-loop-helix) transcription factor, CabHLH10 (Cicer arietinum bHLH10) underlying a major QTL, along with its derived natural alleles/haplotypes governing yield traits under drought stress in chickpea. CabHLH10 binds to a cis-regulatory G-box promoter element to modulate the expression of RD22 (responsive to desiccation 22), a drought/abscisic acid (ABA)-responsive gene (via a trans-expression QTL), and two strong yield-enhancement photosynthetic efficiency (PE) genes. This, in turn, upregulates other downstream drought-responsive and ABA signaling genes, as well as yield-enhancing PE genes, thus increasing plant adaptation to drought with reduced yield penalty. We showed that a superior allele of CabHLH10 introgressed into the NILs improved root and shoot biomass and PE, thereby enhancing yield and productivity during drought without compromising agronomic performance. Furthermore, overexpression of CabHLH10 in chickpea and Arabidopsis (Arabidopsis thaliana) conferred enhanced drought tolerance by improving root and shoot agro-morphological traits. These findings facilitate translational genomics for crop improvement and the development of genetically tailored, climate-resilient, high-yielding chickpea cultivars.

Whole genome resequencing identifies candidate genes and allelic diagnostic markers for resistance to Ralstonia solanacearum infection in cultivated peanut (Arachis hypogaea L.)

Bacterial wilt disease (BWD), caused by Ralstonia solanacearum is a major challenge for peanut production in China and significantly affects global peanut field productivity. It is imperative to identify genetic loci and putative genes controlling resistance to R. solanacearum (RRS). Therefore, a sequencing-based trait mapping approach termed "QTL-seq" was applied to a recombination inbred line population of 581 individuals from the cross of Yueyou 92 (resistant) and Xinhuixiaoli (susceptible). A total of 381,642 homozygous single nucleotide polymorphisms (SNPs) and 98,918 InDels were identified through whole genome resequencing of resistant and susceptible parents for RRS. Using QTL-seq analysis, a candidate genomic region comprising of 7.2 Mb (1.8-9.0 Mb) was identified on chromosome 12 which was found to be significantly associated with RRS based on combined Euclidean Distance (ED) and SNP-index methods. This candidate genomic region had 180 nonsynonymous SNPs and 14 InDels that affected 75 and 11 putative candidate genes, respectively. Finally, eight nucleotide binding site leucine rich repeat (NBS-LRR) putative resistant genes were identified as the important candidate genes with high confidence. Two diagnostic SNP markers were validated and revealed high phenotypic variation in the different resistant and susceptible RIL lines. These findings advocate the expediency of the QTL-seq approach for precise and rapid identification of candidate genomic regions, and the development of diagnostic markers that are applicable in breeding disease-resistant peanut varieties.

The Key to the Future Lies in the Past: Insights from Grain Legume Domestication and Improvement Should Inform Future Breeding Strategies

Crop domestication is a co-evolutionary process that has rendered plants and animals significantly dependent on human interventions for survival and propagation. Grain legumes have played an important role in the development of Neolithic agriculture some 12,000 years ago. Despite being early companions of cereals in the origin and evolution of agriculture, the understanding of grain legume domestication has lagged behind that of cereals. Adapting plants for human use has resulted in distinct morpho-physiological changes between the wild ancestors and domesticates, and this distinction has been the focus of several studies aimed at understanding the domestication process and the genetic diversity bottlenecks created. Growing evidence from research on archeological remains, combined with genetic analysis and the geographical distribution of wild forms, has improved the resolution of the process of domestication, diversification and crop improvement. In this review, we summarize the significance of legume wild relatives as reservoirs of novel genetic variation for crop breeding programs. We describe key legume features, which evolved in response to anthropogenic activities. Here, we highlight how whole genome sequencing and incorporation of omics-level data have expanded our capacity to monitor the genetic changes accompanying these processes. Finally, we present our perspective on alternative routes centered on de novo domestication and re-domestication to impart significant agronomic advances of novel crops over existing commodities. A finely resolved domestication history of grain legumes will uncover future breeding targets to develop modern cultivars enriched with alleles that improve yield, quality and stress tolerance.

Mechanisms of Vernalization-Induced Flowering in Legumes

Vernalization is the requirement for exposure to low temperatures to trigger flowering. The best knowledge about the mechanisms of vernalization response has been accumulated for Arabidopsis and cereals. In Arabidopsis thaliana, vernalization involves an epigenetic silencing of the MADS-box gene FLOWERING LOCUS C (FLC), which is a flowering repressor. FLC silencing releases the expression of the main flowering inductor FLOWERING LOCUS T (FT), resulting in a floral transition. Remarkably, no FLC homologues have been identified in the vernalization-responsive legumes, and the mechanisms of cold-mediated transition to flowering in these species remain elusive. Nevertheless, legume FT genes have been shown to retain the function of the main vernalization signal integrators. Unlike Arabidopsis, legumes have three subclades of FT genes, which demonstrate distinct patterns of regulation with respect to environmental cues and tissue specificity. This implies complex mechanisms of vernalization signal propagation in the flowering network, that remain largely elusive. Here, for the first time, we summarize the available information on the genetic basis of cold-induced flowering in legumes with a special focus on the role of FT genes.

Exploring Chickpea Germplasm Diversity for Broadening the Genetic Base Utilizing Genomic Resourses

Legume crops provide significant nutrition to humans as a source of protein, omega-3 fatty acids as well as specific macro and micronutrients. Additionally, legumes improve the cropping environment by replenishing the soil nitrogen content. Chickpeas are the second most significant staple legume food crop worldwide behind dry bean which contains 17%–24% protein, 41%–51% carbohydrate, and other important essential minerals, vitamins, dietary fiber, folate, β-carotene, anti-oxidants, micronutrients (phosphorus, calcium, magnesium, iron, and zinc) as well as linoleic and oleic unsaturated fatty acids. Despite these advantages, legumes are far behind cereals in terms of genetic improvement mainly due to far less effort, the bottlenecks of the narrow genetic base, and several biotic and abiotic factors in the scenario of changing climatic conditions. Measures are now called for beyond conventional breeding practices to strategically broadening of narrow genetic base utilizing chickpea wild relatives and improvement of cultivars through advanced breeding approaches with a focus on high yield productivity, biotic and abiotic stresses including climate resilience, and enhanced nutritional values. Desirable donors having such multiple traits have been identified using core and mini core collections from the cultivated gene pool and wild relatives of Chickpea. Several methods have been developed to address cross-species fertilization obstacles and to aid in inter-specific hybridization and introgression of the target gene sequences from wild Cicer species. Additionally, recent advances in “Omics” sciences along with high-throughput and precise phenotyping tools have made it easier to identify genes that regulate traits of interest. Next-generation sequencing technologies, whole-genome sequencing, transcriptomics, and differential genes expression profiling along with a plethora of novel techniques like single nucleotide polymorphism exploiting high-density genotyping by sequencing assays, simple sequence repeat markers, diversity array technology platform, and whole-genome re-sequencing technique led to the identification and development of QTLs and high-density trait mapping of the global chickpea germplasm. These altogether have helped in broadening the narrow genetic base of chickpeas.

Harnessing the hidden allelic diversity of wild Cicer to accelerate genomics-assisted chickpea crop improvement

Chickpea, commonly called Bengal gram or Garbanzo bean, faces a productivity crisis around the globe due to numerous biotic and abiotic stresses. The eroded genetic base of the cultivated Cicer gene pool is becoming a significant bottleneck in developing stress-resilient chickpea cultivars. In this scenario, the crop wild relatives (CWR) of chickpea, with the useful genomic wealth of their wild adaptation, give a ray of hope to improve the genetic background of the cultivated Cicer gene pool. To extrapolate these unearthed genomic diversities of wild, we require a thorough understanding of the pre-historic domestication episodes that are changing their shape with the expansion of the available scientific evidence. Keeping aforesaid in view, the current review article provides a glimpsed overview on several efforts done so far to reveal the mysterious origin and evolution of the Cicer gene pool, along with the constraints in their utilization for chickpea crop improvement. It encapsulates various stress-resilient CWR of chickpea and their use in several pre-breeding programs to develop numerous breeding populations for crop genetic enhancement. Further, this review will recapitulate the significant contributions of structural, functional and comparative genomics, pan-genomics and diverse genomics-assisted breeding strategy in dissecting the untapped trait-specific allelic/gene diversity and domestication pattern behind the CWR of chickpea, along with their potential and promises. We expect the newly explored genetic variations may be used in the breeding programs for re-wilding the cultigens' genomic background to open a new avenue for genetic gain and crop improvement capacity of chickpea.

Bulk segregation analysis in the NGS era: a review of its teenage years

Bulk segregation analysis (BSA) utilizes a strategy of pooling individuals with extreme phenotypes to conduct economical and rapidly linked marker screening or quantitative trait locus (QTL) mapping. With the development of next-generation sequencing (NGS) technology in the past 10 years, BSA methods and technical systems have been gradually developed and improved. At the same time, the ever-decreasing costs of sequencing accelerate NGS-based BSA application in different species, including eukaryotic yeast, grain crops, economic crops, horticultural crops, trees, aquatic animals, and insects. This paper provides a landscape of BSA methods and reviews the BSA development process in the past decade, including the sequencing method for BSA, different populations, different mapping algorithms, associated region threshold determination, and factors affecting BSA mapping. Finally, we summarize related strategies in QTL fine mapping combining BSA.

Comparative Flower Transcriptome Network Analysis Reveals DEGs Involved in Chickpea Reproductive Success during Salinity

Salinity is increasingly becoming a significant problem for the most important yet intrinsically salt-sensitive grain legume chickpea. Chickpea is extremely sensitive to salinity during the reproductive phase. Therefore, it is essential to understand the molecular mechanisms by comparing the transcriptomic dynamics between the two contrasting genotypes in response to salt stress. Chickpea exhibits considerable genetic variation amongst improved cultivars, which show better yields in saline conditions but still need to be enhanced for sustainable crop production. Based on previous extensive multi-location physiological screening, two identified genotypes, JG11 (salt-tolerant) and ICCV2 (salt-sensitive), were subjected to salt stress to evaluate their phenological and transcriptional responses. RNA-Sequencing is a revolutionary tool that allows for comprehensive transcriptome profiling to identify genes and alleles associated with stress tolerance and sensitivity. After the first flowering, the whole flower from stress-tolerant and sensitive genotypes was collected. A total of ~300 million RNA-Seq reads were sequenced, resulting in 2022 differentially expressed genes (DEGs) in response to salt stress. Genes involved in flowering time such as FLOWERING LOCUS T (FT) and pollen development such as ABORTED MICROSPORES (AMS), rho-GTPase, and pollen-receptor kinase were significantly differentially regulated, suggesting their role in salt tolerance. In addition to this, we identify a suite of essential genes such as MYB proteins, MADS-box, and chloride ion channel genes, which are crucial regulators of transcriptional responses to salinity tolerance. The gene set enrichment analysis and functional annotation of these genes in flower development suggest that they can be potential candidates for chickpea crop improvement for salt tolerance.

Genomic resources in plant breeding for sustainable agriculture

Climate change during the last 40 years has had a serious impact on agriculture and threatens global food and nutritional security. From over half a million plant species, cereals and legumes are the most important for food and nutritional security. Although systematic plant breeding has a relatively short history, conventional breeding coupled with advances in technology and crop management strategies has increased crop yields by 56 % globally between 1965-85, referred to as the Green Revolution. Nevertheless, increased demand for food, feed, fiber, and fuel necessitates the need to break existing yield barriers in many crop plants. In the first decade of the 21st century we witnessed rapid discovery, transformative technological development and declining costs of genomics technologies. In the second decade, the field turned towards making sense of the vast amount of genomic information and subsequently moved towards accurately predicting gene-to-phenotype associations and tailoring plants for climate resilience and global food security. In this review we focus on genomic resources, genome and germplasm sequencing, sequencing-based trait mapping, and genomics-assisted breeding approaches aimed at developing biotic stress resistant, abiotic stress tolerant and high nutrition varieties in six major cereals (rice, maize, wheat, barley, sorghum and pearl millet), and six major legumes (soybean, groundnut, cowpea, common bean, chickpea and pigeonpea). We further provide a perspective and way forward to use genomic breeding approaches including marker-assisted selection, marker-assisted backcrossing, haplotype based breeding and genomic prediction approaches coupled with machine learning and artificial intelligence, to speed breeding approaches. The overall goal is to accelerate genetic gains and deliver climate resilient and high nutrition crop varieties for sustainable agriculture.

QTL-seq analysis of powdery mildew resistance in a Korean cucumber inbred line

QTL mapping and RT-PCR analyses identified the CsGy5G015660 as a strong powdery mildew resistance candidate gene and natural variation of CsGy5G015660 allele was observed using 115 core germplasm. Powdery mildew (PM) is among the most serious fungal diseases encountered in the cultivation of cucurbits. The development of PM-resistant inbred lines is thus of considerable significance for cucumber breeding programs. In this study, we applied bulked segregant analysis combined with QTL-seq to identify PM resistance loci using F₂ population derived from a cross between two Korean cucumber inbred lines, PM-R (resistant) and PM-S (susceptible). Genome-wide SNP profiling using bulks of the two extreme phenotypes identified two QTLs on chromosomes 5 and 6, designated pm5.2 and pm6.1, respectively. The two PM resistance loci were validated using molecular marker-based classical QTL analysis: pm5.2 (30% R² at LOD 11) and pm6.1 (11% R² at LOD 3.2). Furthermore, reverse transcriptase-PCR analyses, using genes found to be polymorphic between PM-R and PM-S, were conducted to identify the candidate gene(s) responsible for PM resistance. We found that transcripts of the gene CsGy5G015660, encoding a putative leucine-rich repeat receptor-like serine/threonine-protein kinase (RPK2), showed specific accumulation in PM-R prior to the appearance of disease symptoms, and was accordingly considered a strong candidate gene for PM resistance. In addition, cleaved amplified polymorphic sequence markers from CsGy5G015660 were developed and used to screen 35 inbred lines. Natural variation in the CsGy5G015660 allele was also observed based on analysis of a core collection of 115 cucumber accessions. Our results provide new genetic insights for gaining a better understanding of the genetic basis of PM resistance in cucumber, and pave the way for further utilization in cucumber PM resistance breeding programs.

A high-density SNP-based linkage map using genotyping-by-sequencing and its utilization for improved genome assembly of chickpea (Cicer arietinum L.)

Genotyping-by-sequencing (GBS) allows rapid identification of markers for use in development of linkage maps, which expedite efficient breeding programs. In the present study, we have utilized GBS approach to identify and genotype single-nucleotide polymorphism (SNP) markers in an inter-specific RIL population of Cicer arietinum L. X C. reticulatum. A total of 141,639 raw SNPs were identified using the TASSEL-GBS pipeline. After stringent filtering, 8208 candidate SNPs were identified of which ~ 37% were localized in the intragenic regions followed by genic regions (~ 30%) and intergenic regions (~ 27%). We then utilized 6920 stringent selected SNPs from present study and 6714 SNPs and microsatellite markers available from previous studies for construction of linkage map. The resulting high-density linkage map comprising of eight linkage groups contained 13,590 markers which spanned 1299.14 cM of map length with an average marker density of 0.095 cM. Further, the derived linkage map was used to improve the available assembly of desi chickpea genome by anchoring 443 previously unplaced scaffolds onto eight linkage groups. The present efforts have refined anchoring of the desi chickpea genome assembly to 55.57% of the ~ 520 Mb of assembled desi genome. To the best of our knowledge, the linkage map generated in the present study represents one of the most dense linkage map developed for the crop till date. It will serve as a valuable resource for fine mapping and positional cloning of important quantitative trait loci (QTLs) associated with agronomical traits and also for anchoring and ordering of future genome sequence assemblies.

Integrated Genomic Strategies for Cereal Genetic Enhancement: Combining QTL and Association Mapping

Identification of genetic basis for important agronomic traits is essential for marker-assisted crop improvement. Linkage mapping is one of the most popular approaches utilized for identification of major quantitative trait loci (QTLs) governing important agronomic traits in cereals. However, the identified QTLs usually span large genomic intervals and very few of these are subsequently fine mapped to single major effect gene. This hinders application of these QTLs in marker-aided breeding and crop genetic enhancement. On the contrary, association mapping, another popular approach for identification of QTLs, provides very high resolution but suffers from high level of false positives. Joint linkage-association analysis provides a way to combine advantages and avoid the pitfalls associated with both these methods. In this context, we recently developed MetaQTL specific regional association analysis and demonstrated its utility to rapidly narrow down previously identified QTL intervals to few candidate genes. Here, we describe the detailed step-by-step guide for performing MetaQTL specific regional association analysis to identify important genomic regions and underlying potential major effect genes governing traits of agronomic importance in cereals.

Identification of heat-tolerance QTLs and high-temperature stress-responsive genes through conventional QTL mapping, QTL-seq and RNA-seq in tomato

BACKGROUND

High temperature is one of the major abiotic stresses in tomato and greatly reduces fruit yield and quality. Identifying high-temperature stress-responsive (HSR) genes and breeding heat-tolerant varieties is an effective way to address this issue. However, there are few reports on the fine mapping of heat-tolerance quantitative trait locus (QTL) and the identification of HSR genes in tomato. Here, we applied three heat tolerance-related physiological indexes, namely, relative electrical conductivity (REC), chlorophyll content (CC) and maximum photochemical quantum efficiency (F_v/F_m) of PSII (photosystem II), as well as the phenotypic index, the heat injury index (HII), and conventional QTL analysis combined with QTL-seq technology to comprehensively detect heat-tolerance QTLs in tomato seedlings. In addition, we integrated the QTL mapping results with RNA-seq to identify key HSR genes within the major QTLs.

RESULTS

A total of five major QTLs were detected: qHII-1-1, qHII-1-2, qHII-1-3, qHII-2-1 and qCC-1-5 (qREC-1-3). qHII-1-1, qHII-1-2 and qHII-1-3 were located, respectively, in the intervals of 1.43, 1.17 and 1.19 Mb on chromosome 1, while the interval of qHII-2-1 was located in the intervals of 1.87 Mb on chromosome 2. The locations observed with conventional QTL mapping and QTL-seq were consistent. qCC-1-5 and qREC-1-3 for CC and REC, respectively, were located at the same position by conventional QTL mapping. Although qCC-1-5 was not detected in QTL-seq analysis, its phenotypic variation (16.48%) and positive additive effect (0.22) were the highest among all heat tolerance QTLs. To investigate the genes involved in heat tolerance within the major QTLs in tomato, RNA-seq analysis was performed, and four candidate genes (SlCathB2, SlGST, SlUBC5, and SlARG1) associated with heat tolerance were finally detected within the major QTLs by DEG analysis, qRT-PCR screening and biological function analysis.

CONCLUSIONS

In conclusion, this study demonstrated that the combination of conventional QTL mapping, QTL-seq analysis and RNA-seq can rapidly identify candidate genes within major QTLs for a complex trait of interest to replace the fine-mapping process, thus greatly shortening the breeding process and improving breeding efficiency. The results have important applications for the fine mapping and identification of HSR genes and breeding for improved thermotolerance.

Transcriptional signatures modulating shoot apical meristem morphometric and plant architectural traits enhance yield and productivity in chickpea

Plant height (PH) and plant width (PW), two of the major plant architectural traits determining the yield and productivity of a crop, are defined by diverse morphometric characteristics of the shoot apical meristem (SAM). The identification of potential molecular tags from a single gene that simultaneously modulates these plant/SAM architectural traits is therefore prerequisite to achieve enhanced yield and productivity in crop plants, including chickpea. Large-scale multienvironment phenotyping of the association panel and mapping population have ascertained the efficacy of three vital SAM morphometric trait parameters, SAM width, SAM height and SAM area, as key indicators to unravel the genetic basis of the wide PW and PH trait variations observed in desi chickpea. This study integrated a genome-wide association study (GWAS); quantitative trait locus (QTL)/fine-mapping and map-based cloning with molecular haplotyping; transcript profiling; and protein-DNA interaction assays for the dissection of plant architectural traits in chickpea. These exertions delineated natural alleles and superior haplotypes from a CabHLH121 transcription factor (TF) gene within the major QTL governing PW, PH and SAM morphometric traits. A genome-wide protein-DNA interaction assay assured the direct binding of a known stem cell master regulator, CaWUS, to the WOX-homeodomain TF binding sites of a CabHLH121 gene and its constituted haplotypes. The differential expression of CaWUS and transcriptional regulation of its target CabHLH121 gene/haplotypes were apparent, suggesting their collective role in altering SAM morphometric characteristics and plant architectural traits in the contrasting near isogenic lines (NILs). The NILs introgressed with a superior haplotype of a CabHLH121 exhibited optimal PW and desirable PH as well as enhanced yield and productivity without compromising any component of agronomic performance. These molecular signatures of the CabHLH121 TF gene have the potential to regulate both PW and PH traits through the modulation of proliferation, differentiation and maintenance of the meristematic stem cell population in the SAM; therefore, these signatures will be useful in the translational genomic study of chickpea genetic enhancement. The restructured cultivars with desirable PH (semidwarf) and PW will ensure maximal planting density in a specified cultivable field area, thereby enhancing the overall yield and productivity of chickpea. This can essentially facilitate the achievement of better remunerative outputs by farmers with rational land use, therefore ensuring global food security in the present scenario of an increasing population density and shrinking per capita land area.

Identification of quantitative trait loci governing subgynoecy in cucumber

QTL-seq analysis identified three major QTLs conferring subgynoecy in cucumbers. Furthermore, sequence and expression analyses predicted candidate genes controlling subgynoecy. The cucumber (Cucumis sativus L.) is a typical monoecious having individual male and female flowers, and sex differentiation is an important developmental process that directly affects its fruit yield. Subgynoecy represents a sex form with a high degree of femaleness and would have alternative use as gynoecy under limited resource conditions. Recently, many studies have been reported that QTL-seq, which integrates the advantages of bulked segregant analysis and high-throughput whole-genome resequencing, can be a rapid and cost-effective way of mapping QTLs. Segregation analysis in the F₂ and BC₁ populations derived from a cross between subgynoecious LOSUAS and monoecious BMB suggested the quantitative nature of subgynoecy in cucumbers. Both genome-wide SNP profiling of subgynoecious and monoecious bulks constructed from F₂ and BC₁ plants consistently identified three significant genomic regions, one on chromosome 3 (sg3.1) and another two on short and long arms of chromosome 1 (sg1.1 and sg1.2). Classical QTL analysis using the F₂ confirmed sg3.1 (R² = 42%), sg1.1 (R² = 29%) and sg1.2 (R² = 18%) as major QTLs. These results revealed the unique genetic inheritance of subgynoecious line LOSUAS through two distinct major QTLs, sg3.1 and sg1.1, which mainly increase degree of femaleness, while another QTL, sg1.2, contributes to decrease it. This study demonstrated that QTL-seq allows rapid and powerful detection of QTLs using preliminary generation mapping populations such as F₂ or BC₁ population and further that the identified QTLs could be useful for molecular breeding of cucumber lines with high yield potential.

ABC Transporter-Mediated Transport of Glutathione Conjugates Enhances Seed Yield and Quality in Chickpea

The identification of functionally relevant molecular tags is vital for genomics-assisted crop improvement and enhancement of seed yield, quality, and productivity in chickpea (Cicer arietinum). The simultaneous improvement of yield/productivity as well as quality traits often requires pyramiding of multiple genes, which remains a major hurdle given various associated epistatic and pleotropic effects. Unfortunately, no single gene that can improve yield/productivity along with quality and other desirable agromorphological traits is known, hampering the genetic enhancement of chickpea. Using a combinatorial genomics-assisted breeding and functional genomics strategy, this study identified natural alleles and haplotypes of an ABCC3-type transporter gene that regulates seed weight, an important domestication trait, by transcriptional regulation and modulation of the transport of glutathione conjugates in seeds of desi and kabuli chickpea. The superior allele/haplotype of this gene introgressed in desi and kabuli near-isogenic lines enhances the seed weight, yield, productivity, and multiple desirable plant architecture and seed-quality traits without compromising agronomic performance. These salient findings can expedite crop improvement endeavors and the development of nutritionally enriched high-yielding cultivars in chickpea.

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

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

Current advances in chickpea genomics: applications and future perspectives

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