RNA-Seq analysis revealed genes associated with drought stress response in kabuli chickpea (Cicer arietinum L.)

Global transcriptome analysis identifies critical functional modules associated with multiple abiotic stress responses in microalgae Chromochloris zofingiensis

In the current study, systems biology approach was applied to get a deep insight regarding the regulatory mechanisms of Chromochloris zofingiensis under overall stress conditions. Meta-analysis was performed using p-values combination of differentially expressed genes. To identify the informative models related to stress conditions, two distinct weighted gene co-expression networks were constructed and preservation analyses were performed using medianRankand Zsummary algorithms. Moreover, functional enrichment analysis of non-preserved modules was performed to shed light on the biological performance of underlying genes in the non-preserved modules. In the next step, the gene regulatory networks between top hub genes of non-preserved modules and transcription factors were inferred using ensemble of trees algorithm. Results showed that the power of beta = 7 was the best soft-thresholding value to ensure a scale-free network, leading to the determination of 12 co-expression modules with an average size of 128 genes. Preservation analysis showed that the connectivity pattern of the six modules including the blue, black, yellow, pink, greenyellow, and turquoise changed during stress condition which defined as non-preserved modules. Examples of enriched pathways in non-preserved modules were Oxidative phosphorylation", "Vitamin B6 metabolism", and "Arachidonic acid metabolism". Constructed regulatory network between identified TFs and top hub genes of non-preserved module such as Cz06g10250, Cz03g12130 showed that some specific TFs such as C3H and SQUAMOSA promoter binding protein (SBP) specifically regulates the specific hubs. The current findings add substantially to our understanding of the stress responsive underlying mechanism of C. zofingiensis for future studies and metabolite production programs.

Integrated multi-omics analysis reveals drought stress response mechanism in chickpea (Cicer arietinum L.)

Drought is one of the major constraints limiting chickpea productivity. To unravel complex mechanisms regulating drought response in chickpea, we generated transcriptomics, proteomics, and metabolomics datasets from root tissues of four contrasting drought-responsive chickpea genotypes: ICC 4958, JG 11, and JG 11+ (drought-tolerant), and ICC 1882 (drought-sensitive) under control and drought stress conditions. Integration of transcriptomics and proteomics data identified enriched hub proteins encoding isoflavone 4'-O-methyltransferase, UDP-d-glucose/UDP-d-galactose 4-epimerase, and delta-1-pyrroline-5-carboxylate synthetase. These proteins highlighted the involvement of pathways such as antibiotic biosynthesis, galactose metabolism, and isoflavonoid biosynthesis in activating drought stress response mechanisms. Subsequently, the integration of metabolomics data identified six metabolites (fructose, galactose, glucose, myoinositol, galactinol, and raffinose) that showed a significant correlation with galactose metabolism. Integration of root-omics data also revealed some key candidate genes underlying the drought-responsive "QTL-hotspot" region. These results provided key insights into complex molecular mechanisms underlying drought stress response in chickpea.

Chickpea: Its Origin, Distribution, Nutrition, Benefits, Breeding, and Symbiotic Relationship with Mesorhizobium Species

Chickpea (Cicer arietinum L.), encompassing the desi and kabuli varieties, is a beloved pulse crop globally. Its cultivation spans over fifty countries, from the Indian subcontinent and southern Europe to the Middle East, North Africa, the Americas, Australia, and China. With a rich composition of carbohydrates and protein, constituting 80% of its dry seed mass, chickpea is also touted for its numerous health benefits, earning it the title of a 'functional food'. In the past two decades, research has extensively explored the rhizobial diversity associated with chickpea and its breeding in various countries across Europe, Asia, and Oceania, aiming to understand its impact on the sustainable yield and quality of chickpea crops. To date, four notable species of Mesorhizobium-M. ciceri, M. mediterraneum, M. muleiense, and M. wenxiniae-have been reported, originally isolated from chickpea root nodules. Other species, such as M. amorphae, M. loti, M. tianshanense, M. oportunistum, M. abyssinicae, and M. shonense, have been identified as potential symbionts of chickpea, possibly acquiring symbiotic genes through lateral gene transfer. While M. ciceri and M. mediterraneum are widely distributed and studied across chickpea-growing regions, they remain absent in China, where M. muleiense and M. wenxiniae are the sole rhizobial species associated with chickpea. The geographic distribution of chickpea rhizobia is believed to be influenced by factors such as genetic characteristics, competitiveness, evolutionary adaptation to local soil conditions, and compatibility with native soil microbes. Inoculating chickpea with suitable rhizobial strains is crucial when introducing the crop to new regions lacking indigenous chickpea rhizobia. The introduction of a novel chickpea variety, coupled with the effective use of rhizobia for inoculation, offers the potential not only to boost the yield and seed quality of chickpeas, but also to enhance crop productivity within rotation and intercropped systems involving chickpea and other crops. Consequently, this advancement holds the promise to drive forward the cause of sustainable agriculture on a global scale.

Comparative transcriptome analysis of two contrasting genotypes provides new insights into the drought response mechanism in pigeon pea (Cajanus cajan L. Millsp.)

BACKGROUND

Despite plant's ability to adapt and withstand challenging environments, drought poses a severe threat to their growth and development. Although pigeon pea is already quite resistant to drought, the prolonged dehydration induced by the aberrant climate poses a serious threat to their survival and productivity.

OBJECTIVE

Comparative physiological and transcriptome analyses of drought-tolerant (CO5) and drought-sensitive (CO1) pigeon pea genotypes subjected to drought stress were carried out in order to understand the molecular basis of drought tolerance in pigeon pea.

METHODS

The transcriptomic analysis allowed us to examine how drought affects the gene expression of C. cajan. Using bioinformatics tools, the unigenes were de novo assembled, annotated, and functionally evaluated. Additionally, a homology-based sequence search against the droughtDB database was performed to identify the orthologs of the DEGs.

RESULTS

1102 potential drought-responsive genes were found to be differentially expressed genes (DEGs) between drought-tolerant and drought-sensitive genotypes. These included Abscisic acid insensitive 5 (ABI5), Nuclear transcription factor Y subunit A-7 (NF-YA7), WD40 repeat-containing protein 55 (WDR55), Anthocyanidin reductase (ANR) and Zinc-finger homeodomain protein 6 (ZF-HD6) and were highly expressed in the tolerant genotype. Further, GO analysis revealed that the most enriched classes belonged to biosynthetic and metabolic processes in the biological process category, binding and catalytic activity in the molecular function category and nucleus and protein-containing complex in the cellular component category. Results of KEGG pathway analysis revealed that the DEGs were significantly abundant in signalling pathways such as plant hormone signal transduction and MAPK signalling pathways. Consequently, in our investigation, we have identified and validated by qPCR a group of genes involved in signal reception and propagation, stress-specific TFs, and basal regulatory genes associated with drought response.

CONCLUSION

In conclusion, our comprehensive transcriptome dataset enabled the discovery of candidate genes connected to pathways involved in pigeon pea drought response. Our research uncovered a number of unidentified genes and transcription factors that could be used to understand and improve susceptibility to drought.

Identification and expression profile of the SMAX/SMXL family genes in chickpea and lentil provide important players of biotechnological interest involved in plant branching

MAIN CONCLUSION

SMAX/SMXL family genes were successfully identified and characterized in the chickpea and lentil and gene expression data revealed several genes associated with the modulation of plant branching and powerful targets for use in transgenesis and genome editing. Strigolactones (SL) play essential roles in plant growth, rooting, development, and branching, and are associated with plant resilience to abiotic and biotic stress conditions. Likewise, karrikins (KAR) are "plant smoke-derived molecules" that act in a hormonal signaling pathway similar to SL playing an important role in seed germination and hairy root elongation. The SMAX/SMXL family genes are part of these two signaling pathways, in addition to some of these members acting in a still little known SL- and KAR-independent signaling pathway. To date, the identification and functional characterization of the SMAX/SMXL family genes has not been performed in the chickpea and lentil. In this study, nine SMAX/SMXL genes were systematically identified and characterized in the chickpea and lentil, and their expression profiles were explored under different unstressless or different stress conditions. After a comprehensive in silico characterization of the genes, promoters, proteins, and protein-protein interaction network, the expression profile for each gene was determined using a meta-analysis from the RNAseq datasets and complemented with real-time PCR analysis. The expression profiles of the SMAX/SMXL family genes were very dynamic in different chickpea and lentil organs, with some genes assuming a tissue-specific expression pattern. In addition, these genes were significantly modulated by different stress conditions, indicating that SMAX/SMXL genes, although working in three distinct signaling pathways, can act to modulate plant resilience. Most CaSMAX/SMXL and partner genes such as CaTiE1 and CaLAP1, have a positive correlation with the plant branching level, while most LcSMAX/SMXL genes were less correlated with the plant branching level. The SMXL6, SMXL7, SMXL8, TiE1, LAP1, BES1, and BRC1 genes were highlighted as powerful targets for use in transgenesis and genome editing aiming to develop chickpea and lentil cultivars with improved architecture. Therefore, this study presented a detailed characterization of the SMAX/SMXL genes in the chickpea and lentil, and provided new insights for further studies focused on each SMAX/SMXL gene.

Expression profile of the NCED/CCD genes in chickpea and lentil during abiotic stress reveals a positive correlation with increased plant tolerance

Carotenoid cleavage dioxygenase (CCD) gene family is organized in two subfamilies: (i) 9-cis epoxycarotenoid dioxygenase (NCED) genes and (ii) CCD genes. NCED genes are essential for catalyzing the first step of the abscisic-acid (ABA) biosynthesis, while CCD genes produce precursors of the strigolactones hormone. The functional characterization of these gene subfamilies has not been yet performed in chickpea and lentil. Herein, were identified and systematically characterized two NCED and five CCD genes in the chickpea and two NCED and six CCD genes in lentil. After in silico sequence analysis and phylogeny, the expression profile of the NCED/CCD genes was determined by meta-analysis and real-time PCR in plants under different stress conditions. Sequence data revealed that NCED/CCD genes are highly conserved between chickpea and lentil. This conservation was observed both at gene and protein sequence levels and phylogenetic relationships. Analysis of the promoter sequences revealed that all NCED/CCD genes have a considerable number of cis-regulatory elements responsive to biotic and abiotic stress. Protein sequence analysis evidenced that NCED/CCD genes share several conserved motifs and that they have a highly interconnected interaction network. Furthermore, the three-dimensional structure of these proteins was determined and indicated that some proteins have structures with considerable similarity. The meta-analysis revealed that NCED/CCD genes are dynamically modulated in different organs and under different stress conditions, but they have a positive correlation with plant tolerance. In accordance, real-time PCR data showed that both NCED and CCD genes are differentially modulated in plants under drought stress. In particular, CaNCED2, CaCCD5, LcNCED2, LcCCD1, and LcCCD2 genes have a positive correlation with improved plant tolerance to drought stress. Therefore, this study presented a detailed characterization of the chickpea and lentil NCED/CCD genes and provided new insights to improve abiotic stress tolerance in these two important crops.

Genetic mapping of quantitative trait loci associated with drought tolerance in chickpea (Cicer arietinum L.)

Elucidation of the genetic basis of drought tolerance is vital for genomics-assisted breeding of drought tolerant crop varieties. Here, we used genotyping-by-sequencing (GBS) to identify single nucleotide polymorphisms (SNPs) in recombinant inbred lines (RILs) derived from a cross between a drought tolerant chickpea variety, Pusa 362 and a drought sensitive variety, SBD 377. The GBS identified a total of 35,502 SNPs and subsequent filtering of these resulted in 3237 high-quality SNPs included in the eight linkage groups. Fifty-one percent of these SNPs were located in the genic regions distributed throughout the genome. The high density linkage map has total map length of 1069 cm with an average marker interval of 0.33 cm. The linkage map was used to identify 9 robust and consistent QTLs for four drought related traits viz. membrane stability index, relative water content, seed weight and yield under drought, with percent variance explained within the range of 6.29%-90.68% and LOD scores of 2.64 to 6.38, which were located on five of the eight linkage groups. A genomic region on LG 7 harbors quantitative trait loci (QTLs) explaining > 90% phenotypic variance for membrane stability index, and > 10% PVE for yield. This study also provides the first report of major QTLs for physiological traits such as membrane stability index and relative water content for drought stress in chickpea. A total of 369 putative candidate genes were identified in the 6.6 Mb genomic region spanning these QTLs. In-silico expression profiling based on the available transcriptome data revealed that 326 of these genes were differentially expressed under drought stress. KEGG analysis resulted in reduction of candidate genes from 369 to 99, revealing enrichment in various signaling pathways. Haplotype analysis confirmed 5 QTLs among the initially identified 9 QTLs. Two QTLs, qRWC1.1 and qYLD7.1, were chosen based on high SNP density. Candidate gene-based analysis revealed distinct haplotypes in qYLD7.1 associated with significant phenotypic differences, potentially linked to pathways for secondary metabolite biosynthesis. These identified candidate genes bolster defenses through flavonoids and phenylalanine-derived compounds, aiding UV protection, pathogen resistance, and plant structure.The study provides novel genomic regions and candidate genes which can be utilized in genomics-assisted breeding of superior drought tolerant chickpea cultivars.

Transcriptomic analysis of bell pepper (Capsicum annuum L.) revealing key mechanisms in response to low temperature stress

BACKGROUND

Bell pepper (Capsicum annuum L.) is one of the most economically and nutritionally important vegetables worldwide. However, its production can be affected by various abiotic stresses, such as low temperature. This causes various biochemical, morphological and molecular changes affecting membrane lipid composition, photosynthetic pigments, accumulation of free sugars and proline, secondary metabolism, as well as a change in gene expression. However, the mechanism of molecular response to this type of stress has not yet been elucidated.

METHODS AND RESULTS

To further investigate the response mechanism to this abiotic stress, we performed an RNA-Seq transcriptomic analysis to obtain the transcriptomic profile of Capsicum annuum exposed to low temperature stress, where libraries were constructed from reads of control and low temperature stress samples, varying on average per treatment from 22,952,190.5-27,305,327 paired reads ranging in size from 30 to 150 bp. The number of differentially expressed genes (DEGs) for each treatment was 388, 417 and 664 at T-17 h, T-22 h and T-41 h, respectively, identifying 58 up-regulated genes and 169 down-regulated genes shared among the three exposure times. Likewise, 23 DEGs encoding TFs were identified at T-17 h, 30 DEGs at T-22 h and 47 DEGs at T-42 h, respectively. GO analysis revealed that DEGs were involved in catalytic activity, response to temperature stimulus, oxidoreductase activity, stress response, phosphate ion transport and response to abscisic acid. KEGG pathway analysis identified that DEGs were related to flavonoid biosynthesis, alkaloid biosynthesis and plant circadian rhythm pathways in the case of up-regulated genes, while in the case of down-regulated genes, they pertained to MAPK signaling and plant hormone signal transduction pathways, present at all the three time points of low temperature exposure. Validation of the transcriptomic method was performed by evaluation of five DEGs by quantitative polymerase chain reaction (q-PCR).

CONCLUSIONS

The data obtained in the present study provide new insights into the transcriptome profiles of Capsicum annuum stem in response to low temperature stress. The data generated may be useful for the identification of key candidate genes and molecular mechanisms involved in response to this type of stress.

Transcriptome profiling reveals the expression and regulation of genes associated with Fusarium wilt resistance in chickpea (Cicer arietinum L.)

Fusarium wilt (FW) is one of the most significant biotic stresses limiting chickpea production worldwide. To dissect the molecular mechanism of FW resistance in chickpea, comparative transcriptome analyses of contrasting resistance sources of chickpea genotypes under control and Fusarium oxysporum f. sp. ciceris (Foc) inoculated conditions were performed. The high-throughput transcriptome sequencing generated about 1137 million sequencing reads from 24 samples representing two resistant genotypes, two susceptible genotypes, and two near-isogenic lines under control and stress conditions at two-time points (7th- and 12th-day post-inoculation). The analysis identified 5182 differentially expressed genes (DEGs) between different combinations of chickpea genotypes. Functional annotation of these genes indicated their involvement in various biological processes such as defense response, cell wall biogenesis, secondary metabolism, and disease resistance. A significant number (382) of transcription factor encoding genes exhibited differential expression patterns under stress. Further, a considerable number of the identified DEGs (287) co-localized with previously reported quantitative trait locus for FW resistance. Several resistance/susceptibility-related genes, such as SERINE/THREONINE PROTEIN KINASE, DIRIGENT, and MLO exhibiting contrasting expression patterns in resistant and susceptible genotypes upon Foc inoculation, were identified. The results presented in the study provide valuable insights into the transcriptional dynamics associated with FW stress response in chickpea and provide candidate genes for the development of disease-resistant chickpea cultivars.

Unravelling the molecular mechanism underlying drought stress response in chickpea via integrated multi-omics analysis

Drought stress affects growth and productivity significantly in chickpea. An integrated multi-omics analysis can provide a better molecular-level understanding of drought stress tolerance. In the present study, comparative transcriptome, proteome and metabolome analyses of two chickpea genotypes with contrasting responses to drought stress, ICC 4958 (drought-tolerant, DT) and ICC 1882 (drought-sensitive, DS), was performed to gain insights into the molecular mechanisms underlying drought stress response/tolerance. Pathway enrichment analysis of differentially abundant transcripts and proteins suggested the involvement of glycolysis/gluconeogenesis, galactose metabolism, and starch and sucrose metabolism in the DT genotype. An integrated multi-omics analysis of transcriptome, proteome and metabolome data revealed co-expressed genes, proteins and metabolites involved in phosphatidylinositol signaling, glutathione metabolism and glycolysis/gluconeogenesis pathways, specifically in the DT genotype under drought. These stress-responsive pathways were coordinately regulated by the differentially abundant transcripts, proteins and metabolites to circumvent the drought stress response/tolerance in the DT genotype. The QTL-hotspot associated genes, proteins and transcription factors may further contribute to improved drought tolerance in the DT genotype. Altogether, the multi-omics approach provided an in-depth understanding of stress-responsive pathways and candidate genes involved in drought tolerance in chickpea.

Genome-wide transcriptional profiling provides clues to molecular mechanisms underlying cold tolerance in chickpea

Chickpea is an important food legume cultivated in several countries. A sudden drop in autumn temperature, freezing winter temperature, and late spring cold events result in significant losses in chickpea production. The current study used RNA sequencing of two cold tolerant (Saral) and sensitive (ILC533) Kabuli chickpea genotypes to identify cold tolerance-associated genes/pathways. A total of 200.85 million raw reads were acquired from the leaf samples by Illumina sequencing, and around 86% of the clean reads (199 million) were mapped to the chickpea reference genome. The results indicated that 3710 (1980 up- and 1730 down-regulated) and 3473 (1972 up- and 1501 down-regulated) genes were expressed differentially under cold stress in the tolerant and sensitive genotypes, respectively. According to the GO enrichment analysis of uniquely down-regulated genes under cold stress in ILC533, photosynthetic membrane, photosystem II, chloroplast part, and photosystem processes were enriched, revealing that the photosynthesis is severely sensitive to cold stress in this sensitive genotype. Many remarkable transcription factors (CaDREB1E, CaMYB4, CaNAC47, CaTCP4, and CaWRKY33), signaling/regulatory genes (CaCDPK4, CaPP2C6, CaMKK2, and CaHSFA3), and protective genes (CaCOR47, CaLEA3, and CaGST) were identified among the cold-responsive genes of the tolerant genotype. These findings would help improve cold tolerance across chickpea genotypes by molecular breeding or genetic engineering.

Molecular Breeding and Drought Tolerance in Chickpea

Cicer arietinum L. is the third greatest widely planted imperative pulse crop worldwide, and it belongs to the Leguminosae family. Drought is the utmost common abiotic factor on plants, distressing their water status and limiting their growth and development. Chickpea genotypes have the natural ability to fight drought stress using certain strategies viz., escape, avoidance and tolerance. Assorted breeding methods, including hybridization, mutation, and marker-aided breeding, genome sequencing along with omics approaches, could be used to improve the chickpea germplasm lines(s) against drought stress. Root features, for instance depth and root biomass, have been recognized as the greatest beneficial morphological factors for managing terminal drought tolerance in the chickpea. Marker-aided selection, for example, is a genomics-assisted breeding (GAB) strategy that can considerably increase crop breeding accuracy and competence. These breeding technologies, notably marker-assisted breeding, omics, and plant physiology knowledge, underlined the importance of chickpea breeding and can be used in future crop improvement programmes to generate drought-tolerant cultivars(s).

Recent advancement in OMICS approaches to enhance abiotic stress tolerance in legumes

The world is facing rapid climate change and a fast-growing global population. It is believed that the world population will be 9.7 billion in 2050. However, recent agriculture production is not enough to feed the current population of 7.9 billion people, which is causing a huge hunger problem. Therefore, feeding the 9.7 billion population in 2050 will be a huge target. Climate change is becoming a huge threat to global agricultural production, and it is expected to become the worst threat to it in the upcoming years. Keeping this in view, it is very important to breed climate-resilient plants. Legumes are considered an important pillar of the agriculture production system and a great source of high-quality protein, minerals, and vitamins. During the last two decades, advancements in OMICs technology revolutionized plant breeding and emerged as a crop-saving tool in wake of the climate change. Various OMICs approaches like Next-Generation sequencing (NGS), Transcriptomics, Proteomics, and Metabolomics have been used in legumes under abiotic stresses. The scientific community successfully utilized these platforms and investigated the Quantitative Trait Loci (QTL), linked markers through genome-wide association studies, and developed KASP markers that can be helpful for the marker-assisted breeding of legumes. Gene-editing techniques have been successfully proven for soybean, cowpea, chickpea, and model legumes such as Medicago truncatula and Lotus japonicus. A number of efforts have been made to perform gene editing in legumes. Moreover, the scientific community did a great job of identifying various genes involved in the metabolic pathways and utilizing the resulted information in the development of climate-resilient legume cultivars at a rapid pace. Keeping in view, this review highlights the contribution of OMICs approaches to abiotic stresses in legumes. We envisage that the presented information will be helpful for the scientific community to develop climate-resilient legume cultivars.

Genetic mapping of QTLs for drought tolerance in chickpea (Cicer arietinum L.)

Chickpea yield is severely affected by drought stress, which is a complex quantitative trait regulated by multiple small-effect genes. Identifying genomic regions associated with drought tolerance component traits may increase our understanding of drought tolerance mechanisms and assist in the development of drought-tolerant varieties. Here, a total of 187 F8 recombinant inbred lines (RILs) developed from an interspecific cross between drought-tolerant genotype GPF 2 (Cicer arietinum) and drought-sensitive accession ILWC 292 (C. reticulatum) were evaluated to identify quantitative trait loci (QTLs) associated with drought tolerance component traits. A total of 21 traits, including 12 morpho-physiological traits and nine root-related traits, were studied under rainfed and irrigated conditions. Composite interval mapping identified 31 QTLs at Ludhiana and 23 QTLs at Faridkot locations for morphological and physiological traits, and seven QTLs were identified for root-related traits. QTL analysis identified eight consensus QTLs for six traits and five QTL clusters containing QTLs for multiple traits on linkage groups CaLG04 and CaLG06. The identified major QTLs and genomic regions associated with drought tolerance component traits can be introgressed into elite cultivars using genomics-assisted breeding to enhance drought tolerance in chickpea.

Epigenomics as Potential Tools for Enhancing Magnitude of Breeding Approaches for Developing Climate Resilient Chickpea

Epigenomics has become a significant research interest at a time when rapid environmental changes are occurring. Epigenetic mechanisms mainly result from systems like DNA methylation, histone modification, and RNA interference. Epigenetic mechanisms are gaining importance in classical genetics, developmental biology, molecular biology, cancer biology, epidemiology, and evolution. Epigenetic mechanisms play important role in the action and interaction of plant genes during development, and also have an impact on classical plant breeding programs, inclusive of novel variation, single plant heritability, hybrid vigor, plant-environment interactions, stress tolerance, and performance stability. The epigenetics and epigenomics may be significant for crop adaptability and pliability to ambient alterations, directing to the creation of stout climate-resilient elegant crop cultivars. In this review, we have summarized recent progress made in understanding the epigenetic mechanisms in plant responses to biotic and abiotic stresses and have also tried to provide the ways for the efficient utilization of epigenomic mechanisms in developing climate-resilient crop cultivars, especially in chickpea, and other legume crops.

RNA-Seq Analysis Demonstrates Different Strategies Employed by Tiger Nuts (Cyperus esculentus L.) in Response to Drought Stress

Drought stress, an important abiotic stress, has affected global agricultural production by limiting the yield and the quality of crops. Tiger nuts (Cyperus esculentus L.) are C4 crops in the Cyperaceae family, which have high-quality wholesome ingredients. However, data on mechanisms underlying the response of tiger nuts to drought stress are few. Here, the variety of Jisha 1 and 15% polyethylene glycol (PEG; a drought stress simulator) were used to study the mechanisms of stress response in tiger nuts. Our evaluation of the changes in physiological indicators such as electrolyte leakage (El), malondialdehyde (MDA), hydrogen peroxide (H₂O₂), superoxide anion (O₂⁻) and activities of reactive oxygen species (ROS) showed that 12 h was the most suitable time point to harvest and analyze the response to drought stress. Thereafter, we performed transcriptome (RNA-Seq) analysis in the control (CK) and stress treatment groups and showed that there was a total of 1760 differentially expressed genes (DEGs). Gene Ontology (GO) analysis showed that the DEGs were enriched in abscisic acid (ABA) terms, and pathways such as starch and sucrose metabolism (ko00500), phenylpropanoid biosynthesis (ko00940) and plant hormone signal transduction (ko04075) were significantly enriched in the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. In addition, quantitative real-time PCR (qRT-PCR) analysis of the DEGs demonstrated an upregulation of ABA and lignin content, as well as enzyme activities in enriched pathways, which validated the RNA-Seq data. These results revealed the pathways and mechanisms adopted by the tiger nuts in response to drought stress.

Identification of Candidate Genes Regulating Drought Tolerance in Pearl Millet

Pearl millet is an important crop of the arid and semi-arid ecologies to sustain food and fodder production. The greater tolerance to drought stress attracts us to examine its cellular and molecular mechanisms via functional genomics approaches to augment the grain yield. Here, we studied the drought response of 48 inbreds representing four different maturity groups at the flowering stage. A set of 74 drought-responsive genes were separated into five major phylogenic groups belonging to eight functional groups, namely ABA signaling, hormone signaling, ion and osmotic homeostasis, TF-mediated regulation, molecular adaptation, signal transduction, physiological adaptation, detoxification, which were comprehensively studied. Among the conserved motifs of the drought-responsive genes, the protein kinases and MYB domain proteins were the most conserved ones. Comparative in-silico analysis of the drought genes across millet crops showed foxtail millet had most orthologs with pearl millet. Of 698 haplotypes identified across millet crops, MyC2 and Myb4 had maximum haplotypes. The protein–protein interaction network identified ABI2, P5CS, CDPK, DREB, MYB, and CYP707A3 as major hub genes. The expression assay showed the presence of common as well as unique drought-responsive genes across maturity groups. Drought tolerant genotypes in respective maturity groups were identified from the expression pattern of genes. Among several gene families, ABA signaling, TFs, and signaling proteins were the prospective contributors to drought tolerance across maturity groups. The functionally validated genes could be used as promising candidates in backcross breeding, genomic selection, and gene-editing schemes in pearl millet and other millet crops to increase the yield in drought-prone arid and semi-arid ecologies.

Distinctive Physio-Biochemical Properties and Transcriptional Changes Unfold the Mungbean Cultivars Differing by Their Response to Drought Stress at Flowering Stage

Mungbean is a nutritionally and economically important pulse crop cultivated around Asia, mainly in India. The crop is sensitive to drought at various developmental stages of its growing period. However, there is limited or almost no research on a comparative evaluation of mung-bean plants at the flowering stage under drought conditions. Hence, the aim of this research was to impose the drought stress on two mungbean cultivars VRM (Gg) 1 and CO6 at the flowering stage and assess the physio-biochemical and transcriptional changes. After imposing the drought stress, we found that VRM (Gg) 1 exhibited a low reduction in physiological traits (Chlorophyll, relative water content, and plant dry mass) and high proline content than CO6. Additionally, VRM (Gg) 1 has a low level of H2O2 and MDA contents and higher antioxidant enzymes (SOD, POD, and CAT) activity than CO6 during drought stress. The transcriptional analysis of photosynthesis (PS II-PsbP, PS II-LHC, PS I-PsaG/PsaK, and PEPC 3), antioxidant (SOD 2, POD, CAT 2), and drought-responsive genes (HSP-90, DREB2C, NAC 3 and AREB 2) show that VRM (Gg) 1 had increased transcripts more than CO6 under drought stress. Taken together, VRM (Gg) 1 had a better photosynthetic performance which resulted in fewer reductions in chlorophyll, relative water content, and plant dry mass during drought stress. In addition, higher antioxidative enzyme activities led to lower H2O2 and MDA levels, limiting oxidative damage in VRM (Gg) 1. This was positively correlated with increased transcripts of photosynthesis and antioxidant-related genes in VRM (Gg) 1. Further, the increased transcripts of drought-responsive genes indicate that VRM (Gg) 1 has a better genetic basis against drought stress than CO6. These findings help to understand the mungbean response to drought stress and will aid in the development of genotypes with greater drought tolerance by utilizing natural genetic variants.

Insights into the gene and protein structures of the CaSWEET family members in chickpea (Cicer arietinum), and their gene expression patterns in different organs under various stress and abscisic acid treatments

'Sugars Will Eventually be Exported Transporters' (SWEETs) are a group of sugar transporters that play crucial roles in various biological processes, particularly plant stress responses. However, no information is available yet for the CaSWEET family in chickpea. Here, we identified all putative CaSWEET members in chickpea, and obtained their major characteristics, including physicochemical patterns, chromosomal distribution, subcellular localization, gene organization, conserved motifs and three-dimensional protein structures. Subsequently, we explored available transcriptome data to compare spatiotemporal transcript abundance of CaSWEET genes in various major organs. Finally, we studied the changes in their transcript levels in leaves and/or roots following dehydration and exogenous abscisic acid treatments using RT-qPCR to obtain valuable information underlying their potential roles in chickpea responses to water-stress conditions. Our results provide the first insights into the characteristics of the CaSWEET family members and a foundation for further functional characterizations of selected candidate genes for genetic engineering of chickpea.

Genome-wide identification and expression analysis of the GRAS gene family in response to drought stress in chickpea (Cicer arietinum L.)

The GRAS (gibberellic acid insensitive, repressor of GAI and scarecrow) transcription factors (TFs) regulate diverse biological processes involved in plant growth and development. These TFs are also known to regulate gene expression in response to various abiotic stress factors like cold, drought, etc. In chickpea one of the most devastating abiotic stress factors is terminal drought. The GRAS TF family has not been characterized in chickpea (Cicer arietinum L.) until now. In this study, we report 46 GRAS TF genes (CaGRAS genes) in the chickpea genome. The CaGRAS proteins were categorized into nine subfamilies based on their phylogenetic relationship with known GRAS members of Arabidopsis and soybean. The PAT subfamily was the largest consisting of ten CaGRAS members whereas the LAS subfamily was the smallest with only one member. Gene duplication analysis revealed that segmental duplication was the primary reason for the expansion of this gene family within the chickpea genome. The gene expression levels of CaGRAS genes were analysed using two different chickpea varieties contrasting for drought tolerance trait, i.e., ICC 4958 (drought tolerant) and ICC 1882 (drought sensitive). On exposure to drought stress, the two chickpea genotypes, exhibited differential drought response, which was quantified and estimated in terms of differences in leaf relative water content (RWC). The well-watered or control plants of the drought tolerant variety were able to maintain a higher leaf RWC by the end of the drought stress period, whereas the control plants of the drought sensitive variety continued to show a decline in leaf RWC. The two genotypes also differed in their root morphologies, under well-watered and drought stress conditions. The gene expression analysis revealed a potential role of PAT, SCR, SCL3 and SHR GRAS members in the regulation of differential response to drought, in the root tissues, for both the genotypes. CaGRAS 12 (SCR) was identified as a drought-responsive GRAS TF gene, which could serve as a potential candidate gene for utilization in developing chickpea varieties with improved drought tolerance. This study demonstrates the drought-responsive expression of CaGRAS genes in chickpea and also describes the morpho-physiological response of chickpea plants to drought stress conditions.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-03104-z.

Using mathematical models to evaluate germination rate and seedlings length of chickpea seed (Cicer arietinum L.) to osmotic stress at cardinal temperatures

Cicer arietinum is the 3^rd most important cool season legume crop growing in vast arid and semi-arid regions of the world. A lab experiment was designed using hydrothermal time model (HTT) to investigate the chickpea seed germination (SG) behavior, cardinal temperatures and germination responses across fluctuating temperatures (T_s) and water potentials (Ψ_s). Seeds of chickpea var. NIFA 1995 were germinated at six constant T_s (7, 14, 21, 28, 35 and 42°C) each having the following five water potentials: 0, -0.2, -0.4–0.6 and -0.8 MPa. Germination percentage (G%) decreased significantly at (*P ≤ 0.05) from 86.7% at 28°C in -0.2 MPa to 10% in -0.2 MPa at 7°C. The germination rate (GR = 1/t₅₀) against different T percentiles exhibited that linear increase was observed in the GR pattern above and below the T_o. Based on the confidence intervals of the model coefficients and (R²: 0.96), the average cardinal temperatures were 4.7, 23 and 44.2°C for the base (T_b), optimal (T_o) and ceiling (T_c) temperatures respectively. θT1 value was observed maximum at 28°C in -0.2 MPa and decreases with decreasing Ψ (-0.8 MPa). In comparison with control, the θT2 value was also highest in -0.2 MPa at 28°C. The thermal time (TT) concept is well fitted to germination fraction data in distilled water with an R² value increasing 0.972. The hydro time constant (θH) increased with an increase in T to T_o and then decreased when T>T_o. The ѱ_b(50) irregularly varied with increasing T, σΨ_b was also recorded lowest (0.166 MPa) at 28°C and highest (0.457 MPa) at 7°C. Based on the statistical analysis, cardinal temperatures, hydrothermal time constant (θHTT) and germination findings the HTT gives an insight into the interactive effect of T and Ψ on seed germination time courses under varying environmental conditions.

Comparative transcriptomic and metabolic profiling provides insight into the mechanism by which the autophagy inhibitor 3-MA enhances salt stress sensitivity in wheat seedlings

BACKGROUND

Salt stress hinders plant growth and production around the world. Autophagy induced by salt stress helps plants improve their adaptability to salt stress. However, the underlying mechanism behind this adaptability remains unclear. To obtain deeper insight into this phenomenon, combined metabolomics and transcriptomics analyses were used to explore the coexpression of differentially expressed-metabolite (DEM) and gene (DEG) between control and salt-stressed wheat roots and leaves in the presence or absence of the added autophagy inhibitor 3-methyladenine (3-MA).

RESULTS

The results indicated that 3-MA addition inhibited autophagy, increased ROS accumulation, damaged photosynthesis apparatus and impaired the tolerance of wheat seedlings to NaCl stress. A total of 14,759 DEGs and 554 DEMs in roots and leaves of wheat seedlings were induced by salt stress. DEGs were predominantly enriched in cellular amino acid catabolic process, response to external biotic stimulus, regulation of the response to salt stress, reactive oxygen species (ROS) biosynthetic process, regulation of response to osmotic stress, ect. The DEMs were mostly associated with amino acid metabolism, carbohydrate metabolism, phenylalanine metabolism, carbapenem biosynthesis, and pantothenate and CoA biosynthesis. Further analysis identified some critical genes (gene involved in the oxidative stress response, gene encoding transcription factor (TF) and gene involved in the synthesis of metabolite such as alanine, asparagine, aspartate, glutamate, glutamine, 4-aminobutyric acid, abscisic acid, jasmonic acid, ect.) that potentially participated in a complex regulatory network in the wheat response to NaCl stress. The expression of the upregulated DEGs and DEMs were higher, and the expression of the down-regulated DEGs and DEMs was lower in 3-MA-treated plants under NaCl treatment.

CONCLUSION

3-MA enhanced the salt stress sensitivity of wheat seedlings by inhibiting the activity of the roots and leaves, inhibiting autophagy in the roots and leaves, increasing the content of both H₂O₂ and O₂•-, damaged photosynthesis apparatus and changing the transcriptome and metabolome of salt-stressed wheat seedlings.

RNA-Seq reveals different responses to drought in Neotropical trees from savannas and seasonally dry forests

BACKGROUND

Water is one of the main limiting factors for plant growth and crop productivity. Plants constantly monitor water availability and can rapidly adjust their metabolism by altering gene expression. This leads to phenotypic plasticity, which aids rapid adaptation to climate changes. Here, we address phenotypic plasticity under drought stress by analyzing differentially expressed genes (DEG) in four phylogenetically related neotropical Bignoniaceae tree species: two from savanna, Handroanthus ochraceus and Tabebuia aurea, and two from seasonally dry tropical forests (SDTF), Handroanthus impetiginosus and Handroanthus serratifolius. To the best of our knowledge, this is the first report of an RNA-Seq study comparing tree species from seasonally dry tropical forest and savanna ecosystems.

RESULTS

Using a completely randomized block design with 4 species × 2 treatments (drought and wet) × 3 blocks (24 plants) and an RNA-seq approach, we detected a higher number of DEGs between treatments for the SDTF species H. serratifolius (3153 up-regulated and 2821 down-regulated under drought) and H. impetiginosus (332 and 207), than for the savanna species. H. ochraceus showed the lowest number of DEGs, with only five up and nine down-regulated genes, while T. aurea exhibited 242 up- and 96 down-regulated genes. The number of shared DEGs among species was not related to habitat of origin or phylogenetic relationship, since both T. aurea and H impetiginosus shared a similar number of DEGs with H. serratifolius. All four species shared a low number of enriched gene ontology (GO) terms and, in general, exhibited different mechanisms of response to water deficit. We also found 175 down-regulated and 255 up-regulated transcription factors from several families, indicating the importance of these master regulators in drought response.

CONCLUSION

Our findings show that phylogenetically related species may respond differently at gene expression level to drought stress. Savanna species seem to be less responsive to drought at the transcriptional level, likely due to morphological and anatomical adaptations to seasonal drought. The species with the largest geographic range and widest edaphic-climatic niche, H. serratifolius, was the most responsive, exhibiting the highest number of DEG and up- and down-regulated transcription factors (TF).

Fab Advances in Fabaceae for Abiotic Stress Resilience: From 'Omics' to Artificial Intelligence

Legumes are a better source of proteins and are richer in diverse micronutrients over the nutritional profile of widely consumed cereals. However, when exposed to a diverse range of abiotic stresses, their overall productivity and quality are hugely impacted. Our limited understanding of genetic determinants and novel variants associated with the abiotic stress response in food legume crops restricts its amelioration. Therefore, it is imperative to understand different molecular approaches in food legume crops that can be utilized in crop improvement programs to minimize the economic loss. 'Omics'-based molecular breeding provides better opportunities over conventional breeding for diversifying the natural germplasm together with improving yield and quality parameters. Due to molecular advancements, the technique is now equipped with novel 'omics' approaches such as ionomics, epigenomics, fluxomics, RNomics, glycomics, glycoproteomics, phosphoproteomics, lipidomics, regulomics, and secretomics. Pan-omics-which utilizes the molecular bases of the stress response to identify genes (genomics), mRNAs (transcriptomics), proteins (proteomics), and biomolecules (metabolomics) associated with stress regulation-has been widely used for abiotic stress amelioration in food legume crops. Integration of pan-omics with novel omics approaches will fast-track legume breeding programs. Moreover, artificial intelligence (AI)-based algorithms can be utilized for simulating crop yield under changing environments, which can help in predicting the genetic gain beforehand. Application of machine learning (ML) in quantitative trait loci (QTL) mining will further help in determining the genetic determinants of abiotic stress tolerance in pulses.

Transcriptome analysis reveals key genes associated with root-lesion nematode Pratylenchus thornei resistance in chickpea

The root-lesion nematode, Pratylenchus thornei, is one of the major plant-parasitic nematode species causing significant yield losses in chickpea (Cicer arietinum). In order to identify the underlying mechanisms of resistance to P. thornei, the transcriptomes of control and inoculated roots of three chickpea genotypes viz. D05253 > F3TMWR2AB001 (resistant advanced breeding line), PBA HatTrick (moderately resistant cultivar), and Kyabra (susceptible cultivar) were studied at 20 and 50 days post inoculation using the RNA-seq approach. On analyzing the 633.3 million reads generated, 962 differentially expressed genes (DEGs) were identified. Comparative analysis revealed that the majority of DEGs upregulated in the resistant genotype were downregulated in the moderately resistant and susceptible genotypes. Transcription factor families WRKY and bZIP were uniquely expressed in the resistant genotype. The genes Cysteine-rich receptor-like protein kinase 10, Protein lifeguard-like, Protein detoxification, Bidirectional sugar transporter Sugars Will Eventually be Exported Transporters1 (SWEET1), and Subtilisin-like protease were found to play cross-functional roles in the resistant chickpea genotype against P. thornei. The identified candidate genes for resistance to P. thornei in chickpea can be explored further to develop markers and accelerate the introgression of P. thornei resistance into elite chickpea cultivars.

Root-omics for drought tolerance in cool-season grain legumes

Root traits can be exploited to increase the physiological efficiency of crop water use under drought. Root length, root hairs, root branching, root diameter, and root proliferation rate are genetically defined traits that can help to improve the water productivity potential of crops. Recently, high-throughput phenotyping techniques/platforms have been used to screen the germplasm of major cool-season grain legumes for root traits and their impact on different physiological processes, including nutrient uptake and yield potential. Advances in omics approaches have led to the dissection of genomic, proteomic, and metabolomic structures of these traits. This knowledge facilitates breeders to improve the water productivity and nutrient uptake of cultivars under limited soil moisture conditions in major cool-season grain legumes that usually face terminal drought. This review discusses the advances in root traits and their potential for developing drought-tolerant cultivars in cool-season grain legumes.

Time-course transcriptomic analysis of Petunia ×hybrida leaves under water deficit stress using RNA sequencing

Water deficit limits plant growth and development, resulting in quality loss of horticultural crops. However, there is limited information on gene regulation and signaling pathways related to water deficit stress response at multiple time points. The objective of this research was to investigate global gene expression patterns under water deficit stress to provide an insight into how petunia (Petunia ×hybrida 'Mitchell Diploid') responded in the process of stress. Nine-week-old petunias were irrigated daily or placed under water stress by withholding water. Stressed plants reduced stomatal conductance after five days of water deficit, indicating they perceived stress and initiated stress response mechanisms. To analyze transcriptomic changes at the early stage of water deficit, leaf tissue samples were collected 1, 3, and 5 days after water was withheld for RNA sequencing. Under water deficit stress, 154, 3611, and 980 genes were upregulated and 41, 2806, and 253 genes were downregulated on day 1, 3, and 5, respectively. Gene Ontology analysis revealed that redox homeostasis processes through sulfur and glutathione metabolism pathways, and hormone signal transduction, especially abscisic acid and ethylene, were enriched under water deficit stress. Thirty-four transcription factor families were identified, including members of AP2/ERF, NAC, MYB-related, C2H2, and bZIP families, and TFs in AP2/ERF family was the most abundant in petunia. Interestingly, only one member of GRFs was upregulated on day 1, while most of TFs were differentially expressed on day 3 and/or 5. The transcriptome data from this research will provide valuable molecular resources for understanding the early stages of water stress-responsive networks as well as engineering petunia with enhanced water stress tolerance.

Physiological changes and transcriptome profiling in Saccharum spontaneum L. leaf under water stress and re-watering conditions

As the polyploidy progenitor of modern sugarcane, Saccharum spontaneum is considered to be a valuable resistance source to various biotic and abiotic stresses. However, little has been reported on the mechanism of drought tolerance in S. spontaneum. Herein, the physiological changes of S. spontaneum GXS87-16 at three water-deficit levels (mild, moderate, and severe) and after re-watering during the elongation stage were investigated. RNA sequencing was utilized for global transcriptome profiling of GXS87-16 under severe drought and re-watered conditions. There were significant alterations in the physiological parameters of GXS87-16 in response to drought stress and then recovered differently after re-watering. A total of 1569 differentially expressed genes (DEGs) associated with water stress and re-watering were identified. Notably, the majority of the DEGs were induced by stress. GO functional annotations and KEGG pathway analysis assigned the DEGs to 47 GO categories and 93 pathway categories. The pathway categories were involved in various processes, such as RNA transport, mRNA surveillance, plant hormone signal transduction, and plant-pathogen interaction. The reliability of the RNA-seq results was confirmed by qRT-PCR. This study shed light on the regulatory processes of drought tolerance in S. spontaneum and identifies useful genes for genetic improvement of drought tolerance in sugarcane.

Introgression of "QTL-hotspot" region enhances drought tolerance and grain yield in three elite chickpea cultivars

With an aim of enhancing drought tolerance using a marker-assisted backcrossing (MABC) approach, we introgressed the "QTL-hotspot" region from ICC 4958 accession that harbors quantitative trait loci (QTLs) for several drought-tolerance related traits into three elite Indian chickpea (Cicer arietinum L.) cultivars: Pusa 372, Pusa 362, and DCP 92-3. Of eight simple sequence repeat (SSR) markers in the QTL-hotspot region, two to three polymorphic markers were used for foreground selection with respective cross-combinations. A total of 47, 53, and 46 SSRs were used for background selection in case of introgression lines (ILs) developed in genetic backgrounds of Pusa 372, Pusa 362, and DCP 92-3, respectively. In total, 61 ILs (20 BC₃ F₃ in Pusa 372; 20 BC₂ F₃ in Pusa 362, and 21 BC₃ F₃ in DCP 92-3), with >90% recurrent parent genome recovery were developed. Six improved lines in different genetic backgrounds (e.g. BGM 10216 in Pusa 372; BG 3097 and BG 4005 in Pusa 362; IPC(L4-14), IPC(L4-16), and IPC(L19-1) in DCP 92-3) showed better performance than their respective recurrent parents. BGM 10216, with 16% yield gain over Pusa 372, has been released as Pusa Chickpea 10216 by the Central Sub-Committees on Crop Standards, Notification and Release of Varieties of Agricultural Crops, Ministry of Agriculture and Farmers Welfare, Government of India, for commercial cultivation in India. In summary, this study reports introgression of the QTL-hotspot for enhancing yield under rainfed conditions, development of several introgression lines, and release of Pusa Chickpea 10216 developed through molecular breeding in India.

Differential Regulation of Drought Responses in Two Phaseolus vulgaris Genotypes

Drought is probably the most harmful stress affecting common bean crops. Domestication, worldwide spread and local farming practices has entailed the development of a wide variety of common bean genotypes with different degrees of resistance to water stress. In this work, physiological and molecular responses to water stress have been compared in two common bean accessions, PHA-0683 and PMB-0220, previously identified as highly and moderately resistant to water stress, respectively. Our hypothesis was that only quantitative differences in the expression patterns of key genes should be found if molecular mechanisms regulating drought resistance are similar in the two accessions. However, results presented here indicate that the resistance to drought in PMB-0220 and PHA-0683 common bean accessions is regulated by different molecular mechanisms. Differential regulation of ABA synthesis and ABA signaling related genes among the two genotypes, and the control of the drought-induced senescence have a relevant contribution to the higher resistance level of PHA-0683 accession. Our results also suggest that expression patterns of key senescence-related transcription factors could be considered in the screening for drought resistance in common bean germplasm collections.

De-novo transcriptome analysis unveils differentially expressed genes regulating drought and salt stress response in Panicum sumatrense

Screening the transcriptome of drought tolerant variety of little millet (Panicum sumatrense), a marginally cultivated, nutritionally rich, susbsistent crop, can identify genes responsible for its hardiness and enable identification of new sources of genetic variation which can be used for crop improvement. RNA-Seq generated ~ 230 million reads from control and treated tissues, which were assembled into 86,614 unigenes. In silico differential gene expression analysis created an overview of patterns of gene expression during exposure to drought and salt stress. Separate gene expression profiles for leaf and root tissue revealed the differences in regulatory mechanisms operating in these tissues during exposure to abiotic stress. Several transcription factors were identified and studied for differential expression. 61 differentially expressed genes were found to be common to both tissues under drought and salinity stress and were further validated using qRT-PCR. Transcriptome of P. sumatrense was also used to mine for genic SSR markers relevant to abiotic stress tolerance. This study is first report on a detailed analysis of molecular mechanisms of drought and salinity stress tolerance in a little millet variety. Resources generated in this study can be used as potential candidates for further characterization and to improve abiotic stress tolerance in food crops.

Integrating genomics for chickpea improvement: achievements and opportunities

Integration of genomic technologies with breeding efforts have been used in recent years for chickpea improvement. Modern breeding along with low cost genotyping platforms have potential to further accelerate chickpea improvement efforts. The implementation of novel breeding technologies is expected to contribute substantial improvements in crop productivity. While conventional breeding methods have led to development of more than 200 improved chickpea varieties in the past, still there is ample scope to increase productivity. It is predicted that integration of modern genomic resources with conventional breeding efforts will help in the delivery of climate-resilient chickpea varieties in comparatively less time. Recent advances in genomics tools and technologies have facilitated the generation of large-scale sequencing and genotyping data sets in chickpea. Combined analysis of high-resolution phenotypic and genetic data is paving the way for identifying genes and biological pathways associated with breeding-related traits. Genomics technologies have been used to develop diagnostic markers for use in marker-assisted backcrossing programmes, which have yielded several molecular breeding products in chickpea. We anticipate that a sequence-based holistic breeding approach, including the integration of functional omics, parental selection, forward breeding and genome-wide selection, will bring a paradigm shift in development of superior chickpea varieties. There is a need to integrate the knowledge generated by modern genomics technologies with molecular breeding efforts to bridge the genome-to-phenome gap. Here, we review recent advances that have led to new possibilities for developing and screening breeding populations, and provide strategies for enhancing the selection efficiency and accelerating the rate of genetic gain in chickpea.

Comparative Root Transcriptomics Provide Insights into Drought Adaptation Strategies in Chickpea (Cicer arietinum L.)

Drought adversely affects crop production across the globe. The root system immensely contributes to water management and the adaptability of plants to drought stress. In this study, drought-induced phenotypic and transcriptomic responses of two contrasting chickpea (Cicer arietinum L.) genotypes were compared at the vegetative, reproductive transition, and reproductive stages. At the vegetative stage, drought-tolerant genotype maintained higher root biomass, length, and surface area under drought stress as compared to sensitive genotype. However, at the reproductive stage, root length and surface area of tolerant genotype was lower but displayed higher root diameter than sensitive genotype. The shoot biomass of tolerant genotype was overall higher than the sensitive genotype under drought stress. RNA-seq analysis identified genotype- and developmental-stage specific differentially expressed genes (DEGs) in response to drought stress. At the vegetative stage, a total of 2161 and 1873 DEGs, and at reproductive stage 4109 and 3772 DEGs, were identified in the tolerant and sensitive genotypes, respectively. Gene ontology (GO) analysis revealed enrichment of biological categories related to cellular process, metabolic process, response to stimulus, response to abiotic stress, and response to hormones. Interestingly, the expression of stress-responsive transcription factors, kinases, ROS signaling and scavenging, transporters, root nodulation, and oxylipin biosynthesis genes were robustly upregulated in the tolerant genotype, possibly contributing to drought adaptation. Furthermore, activation/repression of hormone signaling and biosynthesis genes was observed. Overall, this study sheds new insights on drought tolerance mechanisms operating in roots with broader implications for chickpea improvement.

Genome-wide transcriptome and physiological analyses provide new insights into peanut drought response mechanisms

Drought is one of the main constraints in peanut production in West Texas and eastern New Mexico regions due to the depletion of groundwater. A multi-seasonal phenotypic analysis of 10 peanut genotypes revealed C76-16 (C-76) and Valencia-C (Val-C) as the best and poor performers under deficit irrigation (DI) in West Texas, respectively. In order to decipher transcriptome changes under DI, RNA-seq was performed in C-76 and Val-C. Approximately 369 million raw reads were generated from 12 different libraries of two genotypes subjected to fully irrigated (FI) and DI conditions, of which ~329 million (90.2%) filtered reads were mapped to the diploid ancestors of peanut. The transcriptome analysis detected 4,508 differentially expressed genes (DEGs), 1554 genes encoding transcription factors (TFs) and a total of 514 single nucleotide polymorphisms (SNPs) among the identified DEGs. The comparative analysis between the two genotypes revealed higher and integral tolerance in C-76 through activation of key genes involved in ABA and sucrose metabolic pathways. Interestingly, one SNP from the gene coding F-box protein (Araip.3WN1Q) and another SNP from gene coding for the lipid transfer protein (Aradu.03ENG) showed polymorphism in selected contrasting genotypes. These SNPs after further validation may be useful for performing early generation selection for selecting drought-responsive genotypes.

Developing Climate-Resilient Chickpea Involving Physiological and Molecular Approaches With a Focus on Temperature and Drought Stresses

Chickpea is one of the most economically important food legumes, and a significant source of proteins. It is cultivated in more than 50 countries across Asia, Africa, Europe, Australia, North America, and South America. Chickpea production is limited by various abiotic stresses (cold, heat, drought, salt, etc.). Being a winter-season crop in northern south Asia and some parts of the Australia, chickpea faces low-temperature stress (0-15°C) during the reproductive stage that causes substantial loss of flowers, and thus pods, to inhibit its yield potential by 30-40%. The winter-sown chickpea in the Mediterranean, however, faces cold stress at vegetative stage. In late-sown environments, chickpea faces high-temperature stress during reproductive and pod filling stages, causing considerable yield losses. Both the low and the high temperatures reduce pollen viability, pollen germination on the stigma, and pollen tube growth resulting in poor pod set. Chickpea also experiences drought stress at various growth stages; terminal drought, along with heat stress at flowering and seed filling can reduce yields by 40-45%. In southern Australia and northern regions of south Asia, lack of chilling tolerance in cultivars delays flowering and pod set, and the crop is usually exposed to terminal drought. The incidences of temperature extremes (cold and heat) as well as inconsistent rainfall patterns are expected to increase in near future owing to climate change thereby necessitating the development of stress-tolerant and climate-resilient chickpea cultivars having region specific traits, which perform well under drought, heat, and/or low-temperature stress. Different approaches, such as genetic variability, genomic selection, molecular markers involving quantitative trait loci (QTLs), whole genome sequencing, and transcriptomics analysis have been exploited to improve chickpea production in extreme environments. Biotechnological tools have broadened our understanding of genetic basis as well as plants' responses to abiotic stresses in chickpea, and have opened opportunities to develop stress tolerant chickpea.

Legume genomics and transcriptomics: From classic breeding to modern technologies

Legumes are essential and play a significant role in maintaining food standards and augmenting physiochemical soil properties through the biological nitrogen fixation process. Biotic and abiotic factors are the main factors limiting legume production. Classical breeding methodologies have been explored extensively about the problem of truncated yield in legumes but have not succeeded at the desired rate. Conventional breeding improved legume genotypes but with more resources and time. Recently, the invention of next-generation sequencing (NGS) and high-throughput methods for genotyping have opened new avenues for research and developments in legume studies. During the last decade, genome sequencing for many legume crops documented. Sequencing and re-sequencing of important legume species have made structural variation and functional genomics conceivable. NGS and other molecular techniques such as the development of markers; genotyping; high density genetic linkage maps; quantitative trait loci (QTLs) identification, expressed sequence tags (ESTs), single nucleotide polymorphisms (SNPs); and transcription factors incorporated into existing breeding technologies have made possible the accurate and accelerated delivery of information for researchers. The application of genome sequencing, RNA sequencing (transcriptome sequencing), and DNA sequencing (re-sequencing) provide considerable insights for legume development and improvement programs.

Moreover, RNA-Seq helps to characterize genes, including differentially expressed genes, and can be applied for functional genomics studies, especially when there is limited information available for the studied genomes. Genome-based crop development studies and the availability of genomics data as well as decision-making gears look be specific for breeding programs. This review mainly presents an overview of the path from classical breeding to new emerging genomics tools, which will trigger and accelerate genomics-assisted breeding for recognition of novel genes for yield and quality characters for sustainable legume crop production.

Resistance to Plant-Parasitic Nematodes in Chickpea: Current Status and Future Perspectives

Plant-parasitic nematodes constrain chickpea (Cicer arietinum) production, with annual yield losses estimated to be 14% of total global production. Nematode species causing significant economic damage in chickpea include root-knot nematodes (Meloidogyne artiella, M. incognita, and M. javanica), cyst nematode (Heterodera ciceri), and root-lesion nematode (Pratylenchus thornei). Reduced functionality of roots from nematode infestation leads to water stress and nutrient deficiency, which in turn lead to poor plant growth and reduced yield. Integration of resistant crops with appropriate agronomic practices is recognized as the safest and most practical, economic and effective control strategy for plant-parasitic nematodes. However, breeding for resistance to plant-parasitic nematodes has numerous challenges that originate from the narrow genetic diversity of the C. arietinum cultigen. While levels of resistance to M. artiella, H. ciceri, and P. thornei have been identified in wild Cicer species that are superior to resistance levels in the C. arietinum cultigen, barriers to interspecific hybridization restrict the use of these crop wild relatives, as sources of nematode resistance. Wild Cicer species of the primary genepool, C. reticulatum and C. echinospermum, are the only species that have been used to introgress resistance genes into the C. arietinum cultigen. The availability of genomic resources, including genome sequence and re-sequence information, the chickpea reference set and mini-core collections, and new wild Cicer collections, provide unprecedented opportunities for chickpea improvement. This review surveys progress in the identification of novel genetic sources of nematode resistance in international germplasm collections and recommends genome-assisted breeding strategies to accelerate introgression of nematode resistance into elite chickpea cultivars.

Global Transcriptome Profiling of Multiple Porcine Organs Reveals Toxoplasma gondii-Induced Transcriptional Landscapes

We characterized the porcine tissue transcriptional landscapes that follow Toxoplasma gondii infection. RNAs were isolated from liver, spleen, cerebral cortex, lung, and mesenteric lymph nodes (MLNs) of T. gondii-infected and uninfected (control) pigs at days 6 and 18 postinfection, and were analyzed using next-generation sequencing (RNA-seq). T. gondii altered the expression of 178, 476, 199, 201, and 362 transcripts at 6 dpi and 217, 223, 347, 119, and 161 at 18 dpi in the infected brain, liver, lung, MLNs and spleen, respectively. The differentially expressed transcripts (DETs) were grouped into five expression patterns and 10 sub-clusters. Gene Ontology enrichment and pathway analysis revealed that immune-related genes dominated the overall transcriptomic signature and that metabolic processes, such as steroid biosynthesis, and metabolism of lipid and carboxylic acid, were downregulated in infected tissues. Co-expression network analysis identified transcriptional modules associated with host immune response to infection. These findings not only show how T. gondii infection alters porcine transcriptome in a tissue-specific manner, but also offer a gateway for testing new hypotheses regarding human response to T. gondii infection.

Towards Exploitation of Adaptive Traits for Climate-Resilient Smart Pulses

Pulses are the main source of protein and minerals in the vegetarian diet. These are primarily cultivated on marginal lands with few inputs in several resource-poor countries of the world, including several in South Asia. Their cultivation in resource-scarce conditions exposes them to various abiotic and biotic stresses, leading to significant yield losses. Furthermore, climate change due to global warming has increased their vulnerability to emerging new insect pests and abiotic stresses that can become even more serious in the coming years. The changing climate scenario has made it more challenging to breed and develop climate-resilient smart pulses. Although pulses are climate smart, as they simultaneously adapt to and mitigate the effects of climate change, their narrow genetic diversity has always been a major constraint to their improvement for adaptability. However, existing genetic diversity still provides opportunities to exploit novel attributes for developing climate-resilient cultivars. The mining and exploitation of adaptive traits imparting tolerance/resistance to climate-smart pulses can be accelerated further by using cutting-edge approaches of biotechnology such as transgenics, genome editing, and epigenetics. This review discusses various classical and molecular approaches and strategies to exploit adaptive traits for breeding climate-smart pulses.

Integrated transcriptome, small RNA and degradome sequencing approaches provide insights into Ascochyta blight resistance in chickpea

Ascochyta blight (AB) is one of the major biotic stresses known to limit the chickpea production worldwide. To dissect the complex mechanisms of AB resistance in chickpea, three approaches, namely, transcriptome, small RNA and degradome sequencing were used. The transcriptome sequencing of 20 samples including two resistant genotypes, two susceptible genotypes and one introgression line under control and stress conditions at two time points (3rd and 7th day post inoculation) identified a total of 6767 differentially expressed genes (DEGs). These DEGs were mainly related to pathogenesis-related proteins, disease resistance genes like NBS-LRR, cell wall biosynthesis and various secondary metabolite synthesis genes. The small RNA sequencing of the samples resulted in the identification of 651 miRNAs which included 478 known and 173 novel miRNAs. A total of 297 miRNAs were differentially expressed between different genotypes, conditions and time points. Using degradome sequencing and in silico approaches, 2131 targets were predicted for 629 miRNAs. The combined analysis of both small RNA and transcriptome datasets identified 12 miRNA-mRNA interaction pairs that exhibited contrasting expression in resistant and susceptible genotypes and also, a subset of genes that might be post-transcriptionally silenced during AB infection. The comprehensive integrated analysis in the study provides better insights into the transcriptome dynamics and regulatory network components associated with AB stress in chickpea and, also offers candidate genes for chickpea improvement.

Identification, Structural Characterization and Gene Expression Analysis of Members of the Nuclear Factor-Y Family in Chickpea (Cicer arietinum L.) under Dehydration and Abscisic Acid Treatments

In plants, the Nuclear Factor-Y (NF-Y) transcription factors (TFs), which include three distinct types of NF-YA, NF-YB, and NF-YC TFs, have been identified to play key roles in the regulation of various plant growth and developmental processes under both normal and environmental stress conditions. In this work, a total of 40 CaNF-Y-encoding genes, including eight CaNF-YAs, 21 CaNF-YBs, and 11 CaNF-YCs, were identified in chickpea, and their major gene and protein characteristics were subsequently obtained using various web-based tools. Of our interest, a phylogenetically-based analysis predicted 18 CaNF-Ys (eight CaNF-YAs, seven CaNF-YBs, and three CaNF-YCs) that potentially play roles in chickpea responses to dehydration according to their close relationship with the well-characterized GmNF-Ys in soybean. These results were in good agreement with the enrichment of drought-responsive cis-regulatory motifs and expression patterns obtained from in silico analyses using publically available transcriptome data. Most of the phylogenetically predicted drought-responsive CaNF-Y genes (15 of 18) were quantitatively validated to significantly respond to dehydration treatment in leaves and/or roots, further supporting the results of in silico analyses. Among these CaNF-Y genes, the transcript levels of CaNF-YA01 and CaNF-YC10 were the most highly accumulated in leaves (by approximately eight-fold) and roots (by approximately 18-fold), respectively, by dehydration. Furthermore, 12 of the 18 CaNF-Y genes were found to be responsive to the most well-known stress hormone, namely abscisic acid (ABA), in leaves and/or roots, suggesting that these genes may act in chickpea response to dehydration in ABA-dependent manner. Taken together, our study has provided a comprehensive and fundamental information for further functional analyses of selected CaNF-Y candidate genes, ultimately leading to the improvement of chickpea growth under water-limited conditions.