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

Insight into a region of chickpea (Cicer arietinum L.) Chromosome 2 revealed potential candidate genes linked to Foc4 Fusarium wilt resistance

Two markers on Chromosome 2 of chickpea (Cicer arietinum ) are reportedly associated with resistance to race 4 Fusarium wilt, and are frequently used in breeding. However, the genes in this region that actually confer wilt resistance are unknown. We aimed to characterise them using both in silico approaches and marker trait association (MTA) analysis. Of the 225 protein-encoding genes in this region, 51 showed significant differential expression in two contrasting chickpea genotypes under wilt, with potential involvement in stress response. From a diverse set of 244 chickpea genotypes, two sets of 40 resistant and 40 susceptible genotypes were selected based on disease incidence and amplification pattern of the TA59 marker. All cultivars were further genotyped with 1238 single nucleotide polymorphisms (SNPs) specific to the 51 genes; only seven SNPs were significantly correlated with disease. SNP Ca2_24099002, specific to the LOC101498008 (Transmembrane protein 87A) gene, accounted for the highest phenotypic variance for disease incidence at 16.30%, whereas SNPs Ca2_25166118 and Ca2_27029215, specific to the LOC101494644 (β-glucosidase BoGH3B-like) and LOC101505289 (Putative tRNA pseudouridine synthase) genes, explained 10.51% and 10.50% of the variation, respectively, in the sets with contrasting disease susceptibility. Together with the TA59 and TR19 markers, these SNPs can be used in a chickpea breeding scheme to develop wilt resistance.

Meta QTL analysis for dissecting abiotic stress tolerance in chickpea

BACKGROUND

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

RESULTS

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

De novo transcriptomic analysis of Doum Palm (Hyphaene compressa) revealed an insight into its potential drought tolerance

BACKGROUND

Doum palms (Hyphaene compressa) perform a crucial starring role in the lives of Kenya's arid and semi-arid people for empowerment and sustenance. Despite the crop's potential for economic gain, there is a lack of genetic resources and detailed information about its domestication at the molecular level. Given the doum palm's vast potential as a widely distributed plant in semi-arid and arid climates and a source of many applications, coupled with the current changing climate scenario, it is essential to understand the molecular processes that provide drought resistance to this plant.

RESULTS

Assembly of the first transcriptome of doum palms subjected to water stress generated about 39.97 Gb of RNA-Seq data. The assembled transcriptome revealed 193,167 unigenes with an average length of 1655 bp, with 128,708 (66.63%) successfully annotated in seven public databases. Unigenes exhibited significant differentially expressed genes (DEGs) in well-watered and stressed-treated plants, with 45071 and 42457 accounting for up-regulated and down-regulated DEGs, respectively. GO term, KEGG, and KOG analysis showed that DEGs were functionally enriched cellular processes, metabolic processes, cellular and catalytic activity, metabolism, genetic information processing, signal transduction mechanisms, and posttranslational modification pathways. Transcription factors (TF), such as the MYB, WRKY, NAC family, FAR1, B3, bHLH, and bZIP, were the prominent TF families identified as doum palm DEGs encoding drought stress tolerance.

CONCLUSIONS

This study provides a complete understanding of DEGs involved in drought stress at the transcriptome level in doum palms. This research is, therefore, the foundation for the characterization of potential genes, leading to a clear understanding of its drought stress responses and providing resources for improved genetic modification.

Transcriptomic analysis reveals candidate genes associated with salinity stress tolerance during the early vegetative stage in fababean genotype, Hassawi-2

Abiotic stresses are a significant constraint to plant production globally. Identifying stress-related genes can aid in the development of stress-tolerant elite genotypes and facilitate trait and crop manipulation. The primary aim of this study was to conduct whole transcriptome analyses of the salt-tolerant faba bean genotype, Hassawi-2, under different durations of salt stress (6 h, 12 h, 24 h, 48 h, and 72 h) at the early vegetative stage, to better understand the molecular basis of salt tolerance. After de novo assembly, a total of 140,308 unigenes were obtained. The up-regulated differentially expressed genes (DEGs) were 2380, 2863, 3057, 3484, and 4820 at 6 h, 12 h, 24 h, 48 h, and 72 h of salt stress, respectively. Meanwhile, 1974, 3436, 2371, 3502, and 5958 genes were downregulated at 6 h, 12 h, 24 h, 48 h, and 72 h of salt stress, respectively. These DEGs encoded various regulatory and functional proteins, including kinases, plant hormone proteins, transcriptional factors (TFs) basic helix-loop-helix (bHLH), Myeloblastosis (MYB), and (WRKY), heat shock proteins (HSPs), late embryogenesis abundant (LEA) proteins, dehydrin, antioxidant enzymes, and aquaporin proteins. This suggests that the faba bean genome possesses an abundance of salinity resistance genes, which trigger different adaptive mechanisms under salt stress. Some selected DEGs validated the RNA sequencing results, thus confirming similar gene expression levels. This study represents the first transcriptome analysis of faba bean leaves subjected to salinity stress offering valuable insights into the mechanisms governing salt tolerance in faba bean during the vegetative stage. This comprehensive investigation enhances our understanding of precise gene regulatory mechanisms and holds promise for the development of novel salt-tolerant faba bean salt-tolerant cultivars.

Genome-wide discovery of di-nucleotide SSR markers based on whole genome re-sequencing data of Cicer arietinum L. and Cicer reticulatum Ladiz

Simple sequence repeats (SSRs) are valuable genetic markers due to their co-dominant inheritance, multi-allelic and reproducible nature. They have been largely used for exploiting genetic architecture of plant germplasms, phylogenetic analysis, and mapping studies. Among the SSRs, di-nucleotide repeats are the most frequent of the simple repeats distributed throughout the plant genomes. In present study, we aimed to discover and develop di-nucleotide SSR markers by using the whole genome re-sequencing (WGRS) data from Cicer arietinum L. and C. reticulatum Ladiz. A total of 35,329 InDels were obtained in C. arietinum, whereas 44,331 InDels in C. reticulatum. 3387 InDels with 2 bp length were detected in C. arietinum, there were 4704 in C. reticulatum. Among 8091 InDels, 58 di-nucleotide regions that were polymorphic between two species were selected and used for validation. We tested primers for evaluation of genetic diversity in 30 chickpea genotypes including C. arietinum, C. reticulatum, C. echinospermum P.H. Davis, C. anatolicum Alef., C. canariense A. Santos & G.P. Lewis, C. microphyllum Benth., C. multijugum Maesen, C. oxyodon Boiss. & Hohen. and C. songaricum Steph ex DC. A total of 244 alleles were obtained for 58 SSR markers giving an average of 2.36 alleles per locus. The observed heterozygosity was 0.08 while the expected heterozygosity was 0.345. Polymorphism information content was found to be 0.73 across all loci. Phylogenetic tree and principal coordinate analysis clearly divided the accessions into four groups. The SSR markers were also evaluated in 30 genotypes of a RIL population obtained from an interspecific cross between C. arietinum and C. reticulatum. Chi-square (χ²) test revealed an expected 1:1 segregation ratio in the population. These results demonstrated the success of SSR identification and marker development for chickpea with the use of WGRS data. The newly developed 58 SSR markers are expected to be useful for chickpea breeders.

Transcriptome analysis of ovules offers early developmental clues after fertilization in Cicer arietinum L

Chickpea (Cicer arietinum L.) seeds are valued for their nutritional scores and limited information on the molecular mechanisms of chickpea fertilization and seed development is available. In the current work, comparative transcriptome analysis was performed on two different stages of chickpea ovules (pre- and post-fertilization) to identify key regulatory transcripts. Two-staged transcriptome sequencing was generated and over 208 million reads were mapped to quantify transcript abundance during fertilization events. Mapping to the reference genome showed that the majority (92.88%) of high-quality Illumina reads were aligned to the chickpea genome. Reference-guided genome and transcriptome assembly yielded a total of 28,783 genes. Of these, 3399 genes were differentially expressed after the fertilization event. These involve upregulated genes including a protease-like secreted in CO(2) response (LOC101500970), amino acid permease 4-like (LOC101506539), and downregulated genes MYB-related protein 305-like (LOC101493897), receptor like protein 29 (LOC101491695). WGCNA analysis and pairwise comparison of datasets, successfully constructed four co-expression modules. Transcription factor families including bHLH, MYB, MYB-related, C2H2 zinc finger, ERF, WRKY and NAC transcription factor were also found to be activated after fertilization. Activation of these genes and transcription factors results in the accumulation of carbohydrates and proteins by enhancing their trafficking and biosynthesis. Total 17 differentially expressed genes, were randomly selected for qRT-PCR for validation of transcriptome analysis and showed statistically significant correlations with the transcriptome data. Our findings provide insights into the regulatory mechanisms underlying changes in fertilized chickpea ovules. This work may come closer to a comprehensive understanding of the mechanisms that initiate developmental events in chickpea seeds after fertilization.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03599-8.

Weed community changes in saffron+chickpea intercropping under different irrigation management

Saffron (Crocus sativus L.) is among the world's most expensive crops; nevertheless, it struggles to compete with weeds. Non-chemical farming practices, such as intercropping and reduced irrigation, can help to decrease weed problems. Therefore, this study aimed to evaluate the changes in the weed density, biomass and weed diversity under saffron-chickpea intercropping system with two irrigation regimes. The study's treatments included two irrigation regimes, namely one-time irrigation and conventional irrigation (carried out four times from October through May), and six planting ratios of saffron and chickpea, namely saffron sole-crop (C1), chickpea sole-crop (C2) in eight rows, 1:1 (C3), 2:2 (C4), 2:1 (C5), and 3:1 (C6)] as main and sub-plots, respectively. The result showed that the conventional irrigation regimes increased weed diversity, however, it didn't affect the Pielou index. Intercropping ratios decreased weed diversity compared to saffron and chickpea mono-cropping systems. The interaction effect of treatments was significant for weed density and weed biomass. In most intercropping ratios, weed density and weed biomass decreased under one-time irrigation regimes. The lowest values for weed density and biomass were observed with an average of 15.5 plants/m2 and 37.51 g/m2, respectively, under the one-time irrigation regime with C4 intercropping systems. This intercropping system did not show a significant difference with C3. Overall, the results indicate that a one-time irrigation regime and intercropping with chickpea, specifically with a 1:1 saffron-chickpea ratio (C3) and a 2:2 saffron-chickpea ratio (C4), could be effective strategies for weed management in saffron in semiarid cropping systems.

Underutilized legumes: nutrient status and advanced breeding approaches for qualitative and quantitative enhancement

Underutilized/orphan legumes provide food and nutritional security to resource-poor rural populations during periods of drought and extreme hunger, thus, saving millions of lives. The Leguminaceae, which is the third largest flowering plant family, has approximately 650 genera and 20,000 species and are distributed globally. There are various protein-rich accessible and edible legumes, such as soybean, cowpea, and others; nevertheless, their consumption rate is far higher than production, owing to ever-increasing demand. The growing global urge to switch from an animal-based protein diet to a vegetarian-based protein diet has also accelerated their demand. In this context, underutilized legumes offer significant potential for food security, nutritional requirements, and agricultural development. Many of the known legumes like Mucuna spp., Canavalia spp., Sesbania spp., Phaseolus spp., and others are reported to contain comparable amounts of protein, essential amino acids, polyunsaturated fatty acids (PUFAs), dietary fiber, essential minerals and vitamins along with other bioactive compounds. Keeping this in mind, the current review focuses on the potential of discovering underutilized legumes as a source of food, feed and pharmaceutically valuable chemicals, in order to provide baseline data for addressing malnutrition-related problems and sustaining pulse needs across the globe. There is a scarcity of information about underutilized legumes and is restricted to specific geographical zones with local or traditional significance. Around 700 genera and 20,000 species remain for domestication, improvement, and mainstreaming. Significant efforts in research, breeding, and development are required to transform existing local landraces of carefully selected, promising crops into types with broad adaptability and economic viability. Different breeding efforts and the use of biotechnological methods such as micro-propagation, molecular markers research and genetic transformation for the development of underutilized crops are offered to popularize lesser-known legume crops and help farmers diversify their agricultural systems and boost their profitability.

Breeding and Genomic Approaches towards Development of Fusarium Wilt Resistance in Chickpea

Chickpea is an important leguminous crop with potential to provide dietary proteins to both humans and animals. It also ameliorates soil nitrogen through biological nitrogen fixation. The crop is affected by an array of biotic and abiotic factors. Among different biotic stresses, a major fungal disease called Fusarium wilt, caused by Fusarium oxysporum f. sp. ciceris (FOC), is responsible for low productivity in chickpea. To date, eight pathogenic races of FOC (race 0, 1A, and 1B/C, 2-6) have been reported worldwide. The development of resistant cultivars using different conventional breeding methods is very time consuming and depends upon the environment. Modern technologies can improve conventional methods to solve these major constraints. Understanding the molecular response of chickpea to Fusarium wilt can help to provide effective management strategies. The identification of molecular markers closely linked to genes/QTLs has provided great potential for chickpea improvement programs. Moreover, omics approaches, including transcriptomics, metabolomics, and proteomics give scientists a vast viewpoint of functional genomics. In this review, we will discuss the integration of all available strategies and provide comprehensive knowledge about chickpea plant defense against Fusarium wilt.

Chickpea (Cicer arietinum L.) Biology and Biotechnology: From Domestication to Biofortification and Biopharming

Chickpea (Cicer arietinum L.), the world's second most consumed legume crop, is cultivated in more than 50 countries around the world. It is a boon for diabetics and is an excellent source of important nutrients such as vitamins A, C, E, K, B1-B3, B5, B6, B9 and minerals (Fe, Zn, Mg and Ca) which all have beneficial effects on human health. By 2050, the world population can cross 9 billion, and in order to feed the teaming millions, chickpea production should also be increased, as it is a healthy alternative to wheat flour and a boon for diabetics. Moreover, it is an important legume that is crucial for food, nutrition, and health security and the livelihood of the small-scale farmers with poor resources, in developing countries. Although marvelous improvement has been made in the development of biotic and abiotic stress-resistant varieties, still there are many lacunae, and to fulfill that, the incorporation of genomic technologies in chickpea breeding (genomics-assisted breeding, high-throughput and precise-phenotyping and implementation of novel breeding strategies) will facilitate the researchers in developing high yielding, climate resilient, water use efficient, salt-tolerant, insect/pathogen resistant varieties, acceptable to farmers, consumers, and industries. This review focuses on the origin and distribution, nutritional profile, genomic studies, and recent updates on crop improvement strategies for combating abiotic and biotic stresses in chickpea.

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.

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.

De novo Transcriptome Analysis of Drought-Adapted Cluster Bean (Cultivar RGC-1025) Reveals the Wax Regulatory Genes Involved in Drought Resistance

Cluster bean (Cyamopsis tetragonoloba L.) is one of the multipurpose underexplored crops grown as green vegetable and for gum production in dryland areas. Cluster bean is known as relatively tolerant to drought and salinity stress. To elucidate the molecular mechanisms involved in the drought tolerance of cluster bean cultivar RGC-1025, RNA sequencing (RNA-seq) of the drought-stressed and control samples was performed. De novo assembly of the reads resulted in 66,838 transcripts involving 203 pathways. Among these transcripts, differentially expressed gene (DEG) analysis resulted in some of the drought-responsive genes expressing alpha dioxygenase 2, low temperature-induced 65 kDa protein (LDI65), putative vacuolar amino acid transporter, and late embryogenesis abundant protein (LEA 3). The analysis also reported drought-responsive transcription factors (TFs), such as NAC, WRKY, GRAS, and MYB families. The relative expression of genes by qRT-PCR revealed consistency with the DEG analysis. Key genes involved in the wax biosynthesis pathway were mapped using the DEG data analysis. These results were positively correlated with epicuticular wax content and the wax depositions on the leaf surfaces, as evidenced by scanning electron microscope (SEM) image analysis. Further, these findings support the fact that enhanced wax deposits on the leaf surface had played a crucial role in combating the drought stress in cluster beans under drought stress conditions. In addition, this study provided a set of unknown genes and TFs that could be a source of engineering tolerance against drought stress in cluster beans.

A Comprehensive Review on Chickpea (Cicer arietinum L.) Breeding for Abiotic Stress Tolerance and Climate Change Resilience

Chickpea is one of the most important pulse crops worldwide, being an excellent source of protein. It is grown under rain-fed conditions averaging yields of 1 t/ha, far from its potential of 6 t/ha under optimum conditions. The combined effects of heat, cold, drought, and salinity affect species productivity. In this regard, several physiological, biochemical, and molecular mechanisms are reviewed to confer tolerance to abiotic stress. A large collection of nearly 100,000 chickpea accessions is the basis of breeding programs, and important advances have been achieved through conventional breeding, such as germplasm introduction, gene/allele introgression, and mutagenesis. In parallel, advances in molecular biology and high-throughput sequencing have allowed the development of specific molecular markers for the genus Cicer, facilitating marker-assisted selection for yield components and abiotic tolerance. Further, transcriptomics, proteomics, and metabolomics have permitted the identification of specific genes, proteins, and metabolites associated with tolerance to abiotic stress of chickpea. Furthermore, some promising results have been obtained in studies with transgenic plants and with the use of gene editing to obtain drought-tolerant chickpea. Finally, we propose some future lines of research that may be useful to obtain chickpea genotypes tolerant to abiotic stress in a scenario of climate change.

Transcriptome analysis of Auricularia fibrillifera fruit-body responses to drought stress and rehydration

BACKGROUND

Drought stress severely restricts edible fungus production. The genus Auricularia has a rare drought tolerance, a rehydration capability, and is nutrient rich.

RESULTS

The key genes and metabolic pathways involved in drought-stress and rehydration were investigated using a transcriptome analysis to clarify the relevant molecular mechanisms. In total, 173.93 Mb clean reads, 26.09 Gb of data bulk, and 52,954 unigenes were obtained. Under drought-stress and rehydration conditions, 14,235 and 8539 differentially expressed genes, respectively, were detected. 'Tyrosine metabolic', 'caffeine metabolism', 'ribosome', 'phagosome', and 'proline and arginine metabolism', as well as 'peroxisome' and 'mitogen-activated protein kinase signaling' pathways, had major roles in A. fibrillifera responses to drought stress. 'Tyrosine' and 'caffeine metabolism' might reveal unknown mechanisms for the antioxidation of A. fibrillifera under drought-stress conditions. During the rehydration process, 'diterpenoid biosynthesis', 'butanoate metabolism', 'C5-branched dibasic acid', and 'aflatoxin biosynthesis' pathways were significantly enriched. Gibberellins and γ-aminobutyric acid were important in the recovery of A. fibrillifera growth after rehydration. Many genes related to antibiotics, vitamins, and other health-related ingredients were found in A. fibrillifera.

CONCLUSION

These findings suggested that the candidate genes and metabolites involved in crucial biological pathways might regulate the drought tolerance or rehydration of Auricularia, shedding light on the corresponding mechanisms and providing new potential targets for the breeding and cultivation of drought-tolerant fungi.

Genomics Associated Interventions for Heat Stress Tolerance in Cool Season Adapted Grain Legumes

Cool season grain legumes occupy an important place among the agricultural crops and essentially provide multiple benefits including food supply, nutrition security, soil fertility improvement and revenue for farmers all over the world. However, owing to climate change, the average temperature is steadily rising, which negatively affects crop performance and limits their yield. Terminal heat stress that mainly occurred during grain development phases severely harms grain quality and weight in legumes adapted to the cool season, such as lentils, faba beans, chickpeas, field peas, etc. Although, traditional breeding approaches with advanced screening procedures have been employed to identify heat tolerant legume cultivars. Unfortunately, traditional breeding pipelines alone are no longer enough to meet global demands. Genomics-assisted interventions including new-generation sequencing technologies and genotyping platforms have facilitated the development of high-resolution molecular maps, QTL/gene discovery and marker-assisted introgression, thereby improving the efficiency in legumes breeding to develop stress-resilient varieties. Based on the current scenario, we attempted to review the intervention of genomics to decipher different components of tolerance to heat stress and future possibilities of using newly developed genomics-based interventions in cool season adapted grain legumes.

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.

Identification of stable heat tolerance QTLs using inter-specific recombinant inbred line population derived from GPF 2 and ILWC 292

Heat stress during reproductive stages has been leading to significant yield losses in chickpea (Cicer arietinum L.). With an aim of identifying the genomic regions or QTLs responsible for heat tolerance, 187 F8 recombinant inbred lines (RILs) derived from the cross GPF 2 (heat tolerant) × ILWC 292 (heat sensitive) were evaluated under late-sown irrigated (January-May) and timely-sown irrigated environments (November-April) at Ludhiana and Faridkot in Punjab, India for 13 heat tolerance related traits. The pooled ANOVA for both locations for the traits namely days to germination (DG), days to flowering initiation (DFI), days to 50% flowering (DFF), days to 100% flowering (DHF), plant height (PH), pods per plant (NPP), biomass (BIO), grain yield (YLD), 100-seed weight (HSW), harvest index (HI), membrane permeability index (MPI), relative leaf water content (RLWC) and pollen viability (PV)) showed a highly significant difference in RILs. The phenotyping data coupled with the genetic map comprising of 1365 ddRAD-Seq based SNP markers were used for identifying the QTLs for heat tolerance. Composite interval mapping provided a total of 28 and 23 QTLs, respectively at Ludhiana and Faridkot locations. Of these, 13 consensus QTLs for DG, DFI, DFF, DHF, PH, YLD, and MPI have been identified at both locations. Four QTL clusters containing QTLs for multiple traits were identified on the same genomic region at both locations. Stable QTLs for days to flowering can be one of the major factors for providing heat tolerance as early flowering has an advantage of more seed setting due to a comparatively longer reproductive period. Identified QTLs can be used in genomics-assisted breeding to develop heat stress-tolerant high yielding chickpea cultivars.

Drought and Heat Stress in Cool-Season Food Legumes in Sub-Tropical Regions: Consequences, Adaptation, and Mitigation Strategies

Drought and heat stress are two major abiotic stresses that challenge the sustainability of agriculture to a larger extend. The changing and unpredictable climate further aggravates the efforts made by researchers as well as farmers. The stresses during the terminal stage of cool-season food legumes may affect numerous physiological and biochemical reactions that may result in poor yield. The plants possess a good number of adaptative and avoiding mechanisms to sustain the adverse situation. The various agronomic and breeding approaches may help in stress-induced alteration. The physiological and biochemical response of crops to any adverse situation is very important to understand to develop mechanisms and approaches for tolerance in plants. Agronomic approaches like altering the planting time, seed priming, foliar application of various macro and micro nutrients, and the application of rhizobacteria may help in mitigating the adverse effect of heat and drought stress to some extent. Breeding approaches like trait-based selection, inheritance studies of marker-based selection, genetic approaches using the transcriptome and metabolome may further pave the way to select and develop crops with better heat and drought stress adaptation and mitigation.

A comprehensive RNA-Seq-based gene expression atlas of the summer squash (Cucurbita pepo) provides insights into fruit morphology and ripening mechanisms

BACKGROUND

Summer squash (Cucurbita pepo: Cucurbitaceae) are a popular horticultural crop for which there is insufficient genomic and transcriptomic information. Gene expression atlases are crucial for the identification of genes expressed in different tissues at various plant developmental stages. Here, we present the first comprehensive gene expression atlas for a summer squash cultivar, including transcripts obtained from seeds, shoots, leaf stem, young and developed leaves, male and female flowers, fruits of seven developmental stages, as well as primary and lateral roots.

RESULTS

In total, 27,868 genes and 2352 novel transcripts were annotated from these 16 tissues, with over 18,000 genes common to all tissue groups. Of these, 3812 were identified as housekeeping genes, half of which assigned to known gene ontologies. Flowers, seeds, and young fruits had the largest number of specific genes, whilst intermediate-age fruits the fewest. There also were genes that were differentially expressed in the various tissues, the male flower being the tissue with the most differentially expressed genes in pair-wise comparisons with the remaining tissues, and the leaf stem the least. The largest expression change during fruit development was early on, from female flower to fruit two days after pollination. A weighted correlation network analysis performed on the global gene expression dataset assigned 25,413 genes to 24 coexpression groups, and some of these groups exhibited strong tissue specificity.

CONCLUSIONS

These findings enrich our understanding about the transcriptomic events associated with summer squash development and ripening. This comprehensive gene expression atlas is expected not only to provide a global view of gene expression patterns in all major tissues in C. pepo but to also serve as a valuable resource for functional genomics and gene discovery in Cucurbitaceae.

Genome-wide simple sequence repeats (SSR) markers discovered from whole-genome sequence comparisons of multiple spinach accessions

The availability of well-assembled genome sequences and reduced sequencing costs have enabled the resequencing of many additional accessions in several crops, thus facilitating the rapid discovery and development of simple sequence repeat (SSR) markers. Although the genome sequence of inbred spinach line Sp75 is available, previous efforts have resulted in a limited number of useful SSR markers. Identification of additional polymorphic SSR markers will support genetics and breeding research in spinach. This study aimed to use the available genomic resources to mine and catalog a large number of polymorphic SSR markers. A search for SSR loci on six chromosome sequences of spinach line Sp75 using GMATA identified a total of 42,155 loci with repeat motifs of two to six nucleotides in the Sp75 reference genome. Whole-genome sequences (30x) of additional 21 accessions were aligned against the chromosome sequences of the reference genome and in silico genotyped using the HipSTR program by comparing and counting repeat numbers variation across the SSR loci among the accessions. The HipSTR program generated SSR genotype data were filtered for monomorphic and high missing loci, and a final set of the 5986 polymorphic SSR loci were identified. The polymorphic SSR loci were present at a density of 12.9 SSRs/Mb and were physically mapped. Out of 36 randomly selected SSR loci for validation, two failed to amplify, while the remaining were all polymorphic in a set of 48 spinach accessions from 34 countries. Genetic diversity analysis performed using the SSRs allele score data on the 48 spinach accessions showed three main population groups. This strategy to mine and develop polymorphic SSR markers by a comparative analysis of the genome sequences of multiple accessions and computational genotyping of the candidate SSR loci eliminates the need for laborious experimental screening. Our approach increased the efficiency of discovering a large set of novel polymorphic SSR markers, as demonstrated in this report.

Genome-wide simple sequence repeats (SSR) markers discovered from whole-genome sequence comparisons of multiple spinach accessions

The availability of well-assembled genome sequences and reduced sequencing costs have enabled the resequencing of many additional accessions in several crops, thus facilitating the rapid discovery and development of simple sequence repeat (SSR) markers. Although the genome sequence of inbred spinach line Sp75 is available, previous efforts have resulted in a limited number of useful SSR markers. Identification of additional polymorphic SSR markers will support genetics and breeding research in spinach. This study aimed to use the available genomic resources to mine and catalog a large number of polymorphic SSR markers. A search for SSR loci on six chromosome sequences of spinach line Sp75 using GMATA identified a total of 42,155 loci with repeat motifs of two to six nucleotides in the Sp75 reference genome. Whole-genome sequences (30x) of additional 21 accessions were aligned against the chromosome sequences of the reference genome and in silico genotyped using the HipSTR program by comparing and counting repeat numbers variation across the SSR loci among the accessions. The HipSTR program generated SSR genotype data were filtered for monomorphic and high missing loci, and a final set of the 5986 polymorphic SSR loci were identified. The polymorphic SSR loci were present at a density of 12.9 SSRs/Mb and were physically mapped. Out of 36 randomly selected SSR loci for validation, two failed to amplify, while the remaining were all polymorphic in a set of 48 spinach accessions from 34 countries. Genetic diversity analysis performed using the SSRs allele score data on the 48 spinach accessions showed three main population groups. This strategy to mine and develop polymorphic SSR markers by a comparative analysis of the genome sequences of multiple accessions and computational genotyping of the candidate SSR loci eliminates the need for laborious experimental screening. Our approach increased the efficiency of discovering a large set of novel polymorphic SSR markers, as demonstrated in this report.

Ion homeostasis for salinity tolerance in plants: a molecular approach

Soil salinity is one of the major environmental stresses faced by the plants. Sodium chloride is the most important salt responsible for inducing salt stress by disrupting the osmotic potential. Due to various innate mechanisms, plants adapt to the sodic niche around them. Genes and transcription factors regulating ion transport and exclusion such as salt overly sensitive (SOS), Na⁺ /H⁺ exchangers (NHXs), high sodium affinity transporter (HKT) and plasma membrane protein (PMP) are activated during salinity stress and help in alleviating cells of ion toxicity. For salt tolerance in plants signal transduction and gene expression is regulated via transcription factors such as NAM (no apical meristem), ATAF (Arabidopsis transcription activation factor), CUC (cup-shaped cotyledon), Apetala 2/ethylene responsive factor (AP2/ERF), W-box binding factor (WRKY) and basic leucine zipper domain (bZIP). Cross-talk between all these transcription factors and genes aid in developing the tolerance mechanisms adopted by plants against salt stress. These genes and transcription factors regulate the movement of ions out of the cells by opening various membrane ion channels. Mutants or knockouts of all these genes are known to be less salt-tolerant compared to wild-types. Using novel molecular techniques such as analysis of genome, transcriptome, ionome and metabolome of a plant, can help in expanding the understanding of salt tolerance mechanism in plants. In this review, we discuss the genes responsible for imparting salt tolerance under salinity stress through transport dynamics of ion balance and need to integrate high-throughput molecular biology techniques to delineate the issue.

Molecular mapping of quantitative trait loci for ascochyta blight and botrytis grey mould resistance in an inter-specific cross in chickpea (Cicer arietinum L.) using genotyping by sequencing

Ascochyta blight (AB) and botrytis grey mould (BGM) are the most devastating fungal diseases of chickpea worldwide. The wild relative of chickpea, C. reticulatum acc. ILWC 292 was found resistant to BGM whereas, GPF2 (Cicer arietinum L.) is resistant to AB. A total of 187 F₈ Recombinant Inbred Lines (RILs) developed from an inter-specific cross of GPF2 × C. reticulatum acc. ILWC 292 were used to identify quantitative trait loci (QTLs) responsible for resistance to AB and BGM. RILs along with parents were evaluated under artificial epiphytotic field/laboratory conditions for two years. Highly significant differences (P < 0.001) were observed for reaction to both pathogens in both years. Parents and RILs were genotyped-by-sequencing to identify genome wide single nucleotide polymorphism (SNPs). A total of 1365 filtered and parental polymorphic SNPs were used for linkage map construction, of which, 673 SNPs were arranged on eight linkage groups. Composite interval mapping revealed three QTLs for AB and four QTLs for BGM resistance. Out of which, two QTLs for AB and three QTLs for BGM were consistent in both years. These QTLs can be targeted for further fine mapping for deployment of resistance to AB and BGM in elite chickpea cultivars using marker-assisted-selection.

Identification of Differentially Expressed Drought-Responsive Genes in Guar [Cyamopsis tetragonoloba (L.) Taub]

Drought remains one of the most serious environmental stresses because of the continuous reduction in soil moisture, which requires the improvement of crops with features such as drought tolerance. Guar [Cyamopsis tetragonoloba (L.) Taub], a forage and industrial crop, is a nonthirsty plant. However, the information on the transcriptome changes that occur under drought stress in guar is very limited; therefore, a gene expression analysis is necessary in this context. Here, we studied the differentially expressed genes (DEGs) in response to drought stress and their metabolic pathways. RNA-Seq via an expectation-maximization algorithm was used to estimate gene abundance. Subsequently, an Empirical Analysis of Digital Gene Expression Data in the R Bioconductor package was used to identify DEGs. Blast2GO, InterProScan, and the Kyoto Encyclopedia of Genes and Genomes were used to explore functional annotation, protein analysis, enzymes, and metabolic pathways. Transcription factors were identified using the PlantTFDB database. Our study identified 499 upregulated and 191 downregulated genes in response to drought stress. Of those, 32 upregulated and six downregulated genes were deemed as novel genes exclusive to guar. An aggregate of 137 protein families, 306 domains, 12 repeats, and two sites were upregulated. The proton-dependent oligopeptide transporter family and transferase, aquaporin transporter, calcium/calmodulin-dependent/calcium-dependent protein kinase, aspartic peptidase A1 family, UDP-glucuronosyl/UDP-glucosyltransferase, and major intrinsic protein were the most upregulated protein families. The upregulated unigenes were associated with 88 enzymes and 77 KEGG pathways. Finally, the MYB-related, MYB, and ERF transcription factor families were upregulated. These data may be useful for understanding the plant molecular response to drought stress.

Exploitation of Heterosis in Pearl Millet: A Review

The phenomenon of heterosis has fascinated plant breeders ever since it was first described by Charles Darwin in 1876 in the vegetable kingdom and later elaborated by George H Shull and Edward M East in maize during 1908. Heterosis is the phenotypic and functional superiority manifested in the F1 crosses over the parents. Various classical complementation mechanisms gave way to the study of the underlying potential cellular and molecular mechanisms responsible for heterosis. In cereals, such as maize, heterosis has been exploited very well, with the development of many single-cross hybrids that revolutionized the yield and productivity enhancements. Pearl millet (Pennisetum glaucum (L.) R. Br.) is one of the important cereal crops with nutritious grains and lower water and energy footprints in addition to the capability of growing in some of the harshest and most marginal environments of the world. In this highly cross-pollinating crop, heterosis was exploited by the development of a commercially viable cytoplasmic male-sterility (CMS) system involving a three-lines breeding system (A-, B- and R-lines). The first set of male-sterile lines, i.e., Tift 23A and Tift18A, were developed in the early 1960s in Tifton, Georgia, USA. These provided a breakthrough in the development of hybrids worldwide, e.g., Tift 23A was extensively used by Punjab Agricultural University (PAU), Ludhiana, India, for the development of the first single-cross pearl millet hybrid, named Hybrid Bajra 1 (HB 1), in 1965. Over the past five decades, the pearl millet community has shown tremendous improvement in terms of cytoplasmic and nuclear diversification of the hybrid parental lines, which led to a progressive increase in the yield and adaptability of the hybrids that were developed, resulting in significant genetic gains. Lately, the whole genome sequencing of Tift 23D2B1 and re-sequencing of circa 1000 genomes by a consortium led by the International Crops Research Institute for the Semi-Arid Tropics (ICRISAT) has been a significant milestone in the development of cutting-edge genetic and genomic resources in pearl millet. Recently, the application of genomics and molecular technologies has provided better insights into genetic architecture and patterns of heterotic gene pools. Development of whole-genome prediction models incorporating heterotic gene pool models, mapped traits and markers have the potential to take heterosis breeding to a new level in pearl millet. This review discusses advances and prospects in various fronts of heterosis for pearl millet.

Characterization of Genetic and Allelic Diversity Amongst Cultivated and Wild Lentil Accessions for Germplasm Enhancement

Intensive breeding of cultivated lentil has resulted in a relatively narrow genetic base, which limits the options to increase crop productivity through selection. Assessment of genetic diversity in the wild gene pool of lentil, as well as characterization of useful and novel alleles/genes that can be introgressed into elite germplasm, presents new opportunities and pathways for germplasm enhancement, followed by successful crop improvement. In the current study, a lentil collection consisting of 467 wild and cultivated accessions that originated from 10 diverse geographical regions was assessed, to understand genetic relationships among different lentil species/subspecies. A total of 422,101 high-confidence SNP markers were identified against the reference lentil genome (cv. CDC Redberry). Phylogenetic analysis clustered the germplasm collection into four groups, namely, Lens culinaris/Lens orientalis, Lens lamottei/Lens odemensis, Lens ervoides, and Lens nigricans. A weak correlation was observed between geographical origin and genetic relationship, except for some accessions of L. culinaris and L. ervoides. Genetic distance matrices revealed a comparable level of variation within the gene pools of L. culinaris (Nei’s coefficient 0.01468–0.71163), L. ervoides (Nei’s coefficient 0.01807–0.71877), and L. nigricans (Nei’s coefficient 0.02188–1.2219). In order to understand any genic differences at species/subspecies level, allele frequencies were calculated from a subset of 263 lentil accessions. Among all cultivated and wild lentil species, L. nigricans exhibited the greatest allelic differentiation across the genome compared to all other species/subspecies. Major differences were observed on six genomic regions with the largest being on Chromosome 1 (c. 1 Mbp). These results indicate that L. nigricans is the most distantly related to L. culinaris and additional structural variations are likely to be identified from genome sequencing studies. This would provide further insights into evolutionary relationships between cultivated and wild lentil germplasm, for germplasm improvement and introgression.

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.

Updates on Genomic Resources in Chickpea for Crop Improvement

In recent years, rapid advancement has been done in generation of genomic resources for the important legume crop chickpea. Here, we provide an update on important advancements made on availability of genomic resources for this crop. The availability of reference genome and transcriptome sequences, and resequencing of several accessions have enabled the discovery of gene space and molecular markers in chickpea. These resources have helped in elucidating evolutionary relationship and identification of quantitative trait loci for important agronomic traits. Gene expression in different tissues/organs during development and under abiotic/biotic stresses has been interrogated. In addition, single-base resolution DNA methylation patterns in different organs have been analyzed to understand gene regulation. Overall, we provide a consolidated overview of available genomic resources of chickpea that may help in fulfilling the promises for improvement of this important crop.

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.

Transcriptome sequencing, de novo assembly, characterisation of wild accession of blackgram (Vigna mungo var. silvestris) as a rich resource for development of molecular markers and validation of SNPs by high resolution melting (HRM) analysis

BACKGROUND

Blackgram [Vigna mungo (L.) Hepper], is an important legume crop of Asia with limited genomic resources. We report a comprehensive set of genic simple sequence repeat (SSR) and single nucleotide polymorphism (SNPs) markers using Illumina MiSeq sequencing of transcriptome and its application in genetic variation analysis and mapping.

RESULTS

Transcriptome sequencing of immature seeds of wild blackgram, V. mungo var. silvestris by Illumina MiSeq technology generated 1.9 × 10⁷ reads, which were assembled into 40,178 transcripts (TCS) with an average length of 446 bp covering 2.97 GB of the genome. A total of 38,753 CDS (Coding sequences) were predicted from 40,178 TCS and 28,984 CDS were annotated through BLASTX and mapped to GO and KEGG database resulting in 140 unique pathways. The tri-nucleotides were most abundant (39.9%) followed by di-nucleotide (30.2%). About 60.3 and 37.6% of SSR motifs were present in the coding sequences (CDS) and untranslated regions (UTRs) respectively. Among SNPs, the most abundant substitution type were transitions (Ts) (61%) followed by transversions (Tv) type (39%), with a Ts/Tv ratio of 1.58. A total of 2306 DEGs were identified by RNA Seq between wild and cultivar and validation was done by quantitative reverse transcription polymerase chain reaction. In this study, we genotyped SNPs with a validation rate of 78.87% by High Resolution Melting (HRM) Assay.

CONCLUSION

In the present study, 1621genic-SSR and 1844 SNP markers were developed from immature seed transcriptome sequence of blackgram and 31 genic-SSR markers were used to study genetic variations among different blackgram accessions. Above developed markers contribute towards enriching available genomic resources for blackgram and aid in breeding programmes.

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.

Transcript expression profiling in two contrasting cultivars and molecular cloning of a SKP-1 like gene, a component of SCF-ubiquitin proteasome system from mungbean Vigna radiate L

Protein degradation and turnover under various environmental stresses is basically regulated by ubiquitin-proteasome system (UPS), of which SKP1 is a very essential component. Isolation and cloning of an identified potential stress responsive candidate gene SKP1, was successfully done for the first time to fathom the role of SKP1 in drought tolerance at genetic level in drought tolerant mungbean cultivar Pratap, which was screened after a detailed physio-biochemical screening amongst seven popular mungbean cultivars. The cloned gene SKP1 (accession number KX881912) is 550 bp in length, encodes 114 amino acids. It shows high sequence homology with SKP1 from Zea mays (NP_001148633). The protein expression of isolated SKP1 was confirmed by GUS fused expression using a Histochemical assay under control as well as under drought stress. Further, up-regulation in relative expression level of SKP1 in different plant parts under drought stress confirmed its utility as a potential drought responsive candidate gene certainly demanding extensive genetic research for further incorporation in breeding programs. Moreover, the structure of VrSKP1 (Vigna radiata SKP1) has been modelled, validated and an Essential Dynamics (ED) was done on the Molecular Dynamics (MD) simulation trajectories for filtering large-scale concerted motions. Free-energy calculations on the ED revealed a complex free-energy landscape (FEL) implying the conformational diversity of the modelled VrSPK1 protein.

Research Progress and Perspective on Drought Stress in Legumes: A Review

Climate change, food shortage, water scarcity, and population growth are some of the threatening challenges being faced in today's world. Drought stress (DS) poses a constant challenge for agricultural crops and has been considered a severe constraint for global agricultural productivity; its intensity and severity are predicted to increase in the near future. Legumes demonstrate high sensitivity to DS, especially at vegetative and reproductive stages. They are mostly grown in the dry areas and are moderately drought tolerant, but severe DS leads to remarkable production losses. The most prominent effects of DS are reduced germination, stunted growth, serious damage to the photosynthetic apparatus, decrease in net photosynthesis, and a reduction in nutrient uptake. To curb the catastrophic effect of DS in legumes, it is imperative to understand its effects, mechanisms, and the agronomic and genetic basis of drought for sustainable management. This review highlights the impact of DS on legumes, mechanisms, and proposes appropriate management approaches to alleviate the severity of water stress. In our discussion, we outline the influence of water stress on physiological aspects (such as germination, photosynthesis, water and nutrient uptake), growth parameters and yield. Additionally, mechanisms, various management strategies, for instance, agronomic practices (planting time and geometry, nutrient management), plant growth-promoting Rhizobacteria and arbuscular mycorrhizal fungal inoculation, quantitative trait loci (QTLs), functional genomics and advanced strategies (CRISPR-Cas9) are also critically discussed. We propose that the integration of several approaches such as agronomic and biotechnological strategies as well as advanced genome editing tools is needed to develop drought-tolerant legume cultivars.

Overexpression of Chickpea Defensin Gene Confers Tolerance to Water-Deficit Stress in Arabidopsis thaliana

Plant defensins are mainly known for their antifungal activity. However, limited information is available regarding their function in abiotic stresses. In this study, a defensin gene, Ca-AFP, from Cicer arietinum, commonly known as chickpea, was cloned and transformed in Arabidopsis thaliana for its functional characterization under simulated water-deficit conditions. Under simulated water-deficit conditions (mannitol and polyethylene glycol-6000 induced), the transgenic A. thaliana plants had higher accumulation of the Ca-AFP transcript compared to that under non-stress condition and showed higher germination rate, root length, and biomass than the wild-type (WT) plants. To get further insights into the role of Ca-AFP in conferring tolerance to water-deficit stress, we determined various physiological parameters and found significant reduction in the transpiration rate and stomatal conductance whereas the net photosynthesis and water use efficiency was increased in the transgenic plants compared to that in the WT plants under water deficit conditions. The transgenic plants showed enhanced superoxide dismutase, ascorbate peroxidase, and catalase activities, had higher proline, chlorophyll, and relative water content, and exhibited reduced ion leakage and malondialdehyde content under water-deficit conditions. Overall, our results indicate that overexpression of Ca-AFP could be an efficient approach for conferring tolerance to water-deficit stress in plants.

The RNA-Seq-based high resolution gene expression atlas of chickpea (Cicer arietinum L.) reveals dynamic spatio-temporal changes associated with growth and development

Chickpea is one of the world's largest cultivated food legumes and is an excellent source of high-quality protein to the human diet. Plant growth and development are controlled by programmed expression of a suite of genes at the given time, stage, and tissue. Understanding how the underlying genome sequence translates into specific plant phenotypes at key developmental stages, information on gene expression patterns is crucial. Here, we present a comprehensive Cicer arietinum Gene Expression Atlas (CaGEA) across different plant developmental stages and organs covering the entire life cycle of chickpea. One of the widely used drought tolerant cultivars, ICC 4958 has been used to generate RNA-Seq data from 27 samples at 5 major developmental stages of the plant. A total of 816 million raw reads were generated and of these, 794 million filtered reads after quality control (QC) were subjected to downstream analysis. A total of 15,947 unique number of differentially expressed genes across different pairwise tissue combinations were identified. Significant differences in gene expression patterns contributing in the process of flowering, nodulation, and seed and root development were inferred in this study. Furthermore, differentially expressed candidate genes from "QTL-hotspot" region associated with drought stress response in chickpea were validated.

The early transcriptome response of cassava (Manihot esculenta Crantz) to mealybug (Phenacoccus manihoti) feeding

The mealybug, Phenacoccus manihoti, is a leading pest of cassava (Manihot esculenta Crantz), damaging this crop globally. Although the biological control of this mealybug using natural predators has been established, resistance breeding remains an important means of control. Understanding plant responses to insect herbivory, by determining and identifying differentially expressed genes (DEGs), is a vital step towards the understanding of molecular mechanisms of defence responses in plants and the development of resistant cultivars by gene editing. Morphological and molecular analysis confirmed the mealybug identity as Phenacoccus manihoti (Matile-Ferrero). The transcriptome response of the green mite resistant cassava genotype AR23.1 was compared to P40/1 with no known resistance at 24 and 72 hours of mealybug infestation compared to non-infested mock. A total of 301 and 206 genes were differentially expressed at 24 and 72 of mealybug infestation for AR23.1 and P40/1 genotypes respectively, using a log2 fold change and P-value ≤ 0.05. Gene ontology functional classification revealed an enrichment of genes in the secondary metabolic process category in AR23.1 in comparison with P40/1, while genes in the regulation of molecular function, cellular component biogenesis and electron carrier categories were more significantly enriched in P40/1 than in AR23.1. Biological pathway analysis, based on KEGG, revealed a significant enrichment of plant-pathogen interaction and plant hormonal signal transduction pathways for a cohort of up-regulated and down-regulated DEGs in both genotypes. Defence-related genes such as 2-oxogluterate, gibberellin oxidase and terpene synthase proteins were only induced in genotype AR23.1 and not in P40/1, and subsequently validated by RT-qPCR. The study revealed a difference in response to mealybug infestation in the two genotypes studied, with AR23.1 showing a higher number of differentially expressed transcripts post mealybug infestation at 24 and 72 hours. Candidate defence-related genes that were overexpressed in the AR23.1 genotype post mealybug infestation will be useful in future functional studies towards the control of mealybugs.

Improvement of Drought Tolerance in Rice (Oryza sativa L.): Genetics, Genomic Tools, and the WRKY Gene Family

Drought tolerance is an important quantitative trait with multipart phenotypes that are often further complicated by plant phenology. Different types of environmental stresses, such as high irradiance, high temperatures, nutrient deficiencies, and toxicities, may challenge crops simultaneously; therefore, breeding for drought tolerance is very complicated. Interdisciplinary researchers have been attempting to dissect and comprehend the mechanisms of plant tolerance to drought stress using various methods; however, the limited success of molecular breeding and physiological approaches suggests that we rethink our strategies. Recent genetic techniques and genomics tools coupled with advances in breeding methodologies and precise phenotyping will likely reveal candidate genes and metabolic pathways underlying drought tolerance in crops. The WRKY transcription factors are involved in different biological processes in plant development. This zinc (Zn) finger protein family, particularly members that respond to and mediate stress responses, is exclusively found in plants. A total of 89 WRKY genes in japonica and 97 WRKY genes in O. nivara (OnWRKY) have been identified and mapped onto individual chromosomes. To increase the drought tolerance of rice (Oryza sativa L.), research programs should address the problem using a multidisciplinary strategy, including the interaction of plant phenology and multiple stresses, and the combination of drought tolerance traits with different genetic and genomics approaches, such as microarrays, quantitative trait loci (QTLs), WRKY gene family members with roles in drought tolerance, and transgenic crops. This review discusses the newest advances in plant physiology for the exact phenotyping of plant responses to drought to update methods of analysing drought tolerance in rice. Finally, based on the physiological/morphological and molecular mechanisms found in resistant parent lines, a strategy is suggested to select a particular environment and adapt suitable germplasm to that environment.

The leaf transcriptome of fennel (Foeniculum vulgare Mill.) enables characterization of the t-anethole pathway and the discovery of microsatellites and single-nucleotide variants

Fennel is a plant species of both agronomic and pharmaceutical interest that is characterized by a shortage of genetic and molecular data. Taking advantage of NGS technology, we sequenced and annotated the first fennel leaf transcriptome using material from four different lines and two different bioinformatic approaches: de novo and genome-guided transcriptome assembly. A reference transcriptome for assembly was produced by combining these two approaches. Among the 79,263 transcripts obtained, 47,775 were annotated using BLASTX analysis performed against the NR protein database subset with 11,853 transcripts representing putative full-length CDS. Bioinformatic analyses revealed 1,011 transcripts encoding transcription factors, mainly from the BHLH, MYB-related, C2H2, MYB, and ERF families, and 6,411 EST-SSR regions. Single-nucleotide variants of SNPs and indels were identified among the 8 samples at a frequency of 0.5 and 0.04 variants per Kb, respectively. Finally, the assembled transcripts were screened to identify genes related to the biosynthesis of t-anethole, a compound well-known for its nutraceutical and medical properties. For each of the 11 genes encoding structural enzymes in the t-anethole biosynthetic pathway, we identified at least one transcript showing a significant match. Overall, our work represents a treasure trove of information exploitable both for marker-assisted breeding and for in-depth studies on thousands of genes, including those involved in t-anethole biosynthesis.

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.

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

Drought is the most important constraint that effects chickpea production globally. RNA-Seq has great potential to dissect the molecular mechanisms of tolerance to environmental stresses. Transcriptome profiles in roots and shoots of two contrasting Iranian kabuli chickpea genotypes (Bivanij and Hashem) were investigated under water-limited conditions at early flowering stage using RNA-Seq approach. A total of 4,572 differentially expressed genes (DEGs) were identified. Of these, 261 and 169 drought stress responsive genes were identified in the shoots and the roots, respectively, and 17 genes were common in the shoots and the roots. Gene Ontology (GO) analysis revealed several sub-categories related to the stress, including response to stress, defense response and response to stimulus in the tolerant genotype Bivanij as compared to the sensitive genotype Hashem under drought stress. In addition, several Transcription factors (TFs) were identified in major metabolic pathways such as, ABA, proline and flavonoid biosynthesis. Furthermore, a number of the DEGs were observed in "QTL-hotspot" regions which were reported earlier in chickpea. Drought tolerance dissection in the genotypes revealed that the genes and the pathways involved in shoots of Bivanij were the most important factor to make a difference between the genotypes for drought tolerance. The identified TFs in the experiment, particularly those which were up-regulated in shoots of Bivanij during drought stress, were potential candidates for enhancing tolerance to drought.

De novo transcriptome assembly and identification of salt-responsive genes in sugar beet M14

Sugar beet (Beta vulgaris) is an important crop of sugar production in the world. Previous studies reported that sugar beet monosomic addition line M14 obtained from the intercross between Beta vulgaris L. (cultivated species) and B. corolliflora Zoss (wild species) exhibited tolerance to salt (up to 0.5 M NaCl) stress. To estimate a broad spectrum of genes involved in the M14 salt tolerance will help elucidate the molecular mechanisms underlying salt stress. Comparative transcriptomics was performed to monitor genes differentially expressed in the leaf and root samples of the sugar beet M14 seedlings treated with 0, 200 and 400 mM NaCl, respectively. Digital gene expression revealed that 3856 unigenes in leaves and 7157 unigenes in roots were differentially expressed under salt stress. Enrichment analysis of the differentially expressed genes based on GO and KEGG databases showed that in both leaves and roots genes related to regulation of redox balance, signal transduction, and protein phosphorylation were differentially expressed. Comparison of gene expression in the leaf and root samples treated with 200 and 400 mM NaCl revealed different mechanisms for coping with salt stress. In addition, the expression levels of nine unigenes in the reactive oxygen species (ROS) scavenging system exhibited significant differences in the leaves and roots. Our transcriptomics results have provided new insights into the salt-stress responses in the leaves and roots of sugar beet.

Comparative transcriptome reveal the potential adaptive evolutionary genes in Andrias davidianus

To search the evidence of molecular evolution mechanism for aquatic and cave habitat in Andrias davidianus, the evolution analysis was carried out among several species transcriptome data. The transcriptome data of Notophthalmus viridescens, Xenopus tropicalis, Cynops pyrrhogaster, Hynobius chinensis and A. davidianus were obtained from the Genbank and reassembled except Xenopus tropicalis. The BLAST search of transcriptome data obtained 1244 single-copy orthologous genes among five species. A phylogenetic tree showed A. davidianus to have the closest relationship to H. chinensis. Fourteen positively selected genes were detected in A. davidianus and N. vridescens group and fifteen in A. davidianus and H. chinensis group. Five genes were shared in the both groups which involved in the immune system, suggesting that A. davidianus adaptation to an aquatic and cave environment required rapid evolution of the immune system compared to N. viridescens and H. chinensis.

Mining and Development of Novel SSR Markers Using Next Generation Sequencing (NGS) Data in Plants

Microsatellites, or simple sequence repeats (SSRs), are one of the most informative and multi-purpose genetic markers exploited in plant functional genomics. However, the discovery of SSRs and development using traditional methods are laborious, time-consuming, and costly. Recently, the availability of high-throughput sequencing technologies has enabled researchers to identify a substantial number of microsatellites at less cost and effort than traditional approaches. Illumina is a noteworthy transcriptome sequencing technology that is currently used in SSR marker development. Although 454 pyrosequencing datasets can be used for SSR development, this type of sequencing is no longer supported. This review aims to present an overview of the next generation sequencing, with a focus on the efficient use of de novo transcriptome sequencing (RNA-Seq) and related tools for mining and development of microsatellites in plants.

Improving Salt Tolerance of Chickpea Using Modern Genomics Tools and Molecular Breeding

INTRODUCTION

The high protein value, essential minerals, dietary fibre and notable ability to fix atmospheric nitrogen make chickpea a highly remunerative crop, particularly in low-input food production systems. Of the variety of constraints challenging chickpea productivity worldwide, salinity remains of prime concern owing to the intrinsic sensitivity of the crop. In view of the projected expansion of chickpea into arable and salt-stressed land by 2050, increasing attention is being placed on improving the salt tolerance of this crop. Considerable effort is currently underway to address salinity stress and substantial breeding progress is being made despite the seemingly highly-complex and environment-dependent nature of the tolerance trait.

CONCLUSION

This review aims to provide a holistic view of recent advances in breeding chickpea for salt tolerance. Initially, we focus on the identification of novel genetic resources for salt tolerance via extensive germplasm screening. We then expand on the use of genome-wide and cost-effective techniques to gain new insights into the genetic control of salt tolerance, including the responsive genes/QTL(s), gene(s) networks/cross talk and intricate signalling cascades.

Food Legumes and Rising Temperatures: Effects, Adaptive Functional Mechanisms Specific to Reproductive Growth Stage and Strategies to Improve Heat Tolerance

Ambient temperatures are predicted to rise in the future owing to several reasons associated with global climate changes. These temperature increases can result in heat stress- a severe threat to crop production in most countries. Legumes are well-known for their impact on agricultural sustainability as well as their nutritional and health benefits. Heat stress imposes challenges for legume crops and has deleterious effects on the morphology, physiology, and reproductive growth of plants. High-temperature stress at the time of the reproductive stage is becoming a severe limitation for production of grain legumes as their cultivation expands to warmer environments and temperature variability increases due to climate change. The reproductive period is vital in the life cycle of all plants and is susceptible to high-temperature stress as various metabolic processes are adversely impacted during this phase, which reduces crop yield. Food legumes exposed to high-temperature stress during reproduction show flower abortion, pollen and ovule infertility, impaired fertilization, and reduced seed filling, leading to smaller seeds and poor yields. Through various breeding techniques, heat tolerance in major legumes can be enhanced to improve performance in the field. Omics approaches unravel different mechanisms underlying thermotolerance, which is imperative to understand the processes of molecular responses toward high-temperature stress.