Genome-wide dissection of AP2/ERF and HSP90 gene families in five legumes and expression profiles in chickpea and pigeonpea

Genome-wide identification, classification and expression analysis of the heat shock transcription factor family in Garlic (Allium sativum L.)

BACKGROUND

The heat shock transcription factor (HSF) plays a crucial role in the regulatory network by coordinating responses to heat stress as well as other stress signaling pathways. Despite extensive studies on HSF functions in various plant species, our understanding of this gene family in garlic, an important crop with nutritional and medicinal value, remains limited. In this study, we conducted a comprehensive investigation of the entire garlic genome to elucidate the characteristics of the AsHSF gene family.

RESULTS

In this study, we identified a total of 17 AsHSF transcription factors. Phylogenetic analysis classified these transcription factors into three subfamilies: Class A (9 members), Class B (6 members), and Class C (2 members). Each subfamily was characterized by shared gene structures and conserved motifs. The evolutionary features of the AsHSF genes were investigated through a comprehensive analysis of chromosome location, conserved protein motifs, and gene duplication events. These findings suggested that the evolution of AsHSF genes is likely driven by both tandem and segmental duplication events. Moreover, the nucleotide diversity of the AsHSF genes decreased by only 0.0002% from wild garlic to local garlic, indicating a slight genetic bottleneck experienced by this gene family during domestication. Furthermore, the analysis of cis-acting elements in the promoters of AsHSF genes indicated their crucial roles in plant growth, development, and stress responses. qRT-PCR analysis, co-expression analysis, and protein interaction prediction collectively highlighted the significance of Asa6G04911. Subsequent experimental investigations using yeast two-hybridization and yeast induction experiments confirmed its interaction with HSP70/90, reinforcing its significance in heat stress.

CONCLUSIONS

This study is the first to unravel and analyze the AsHSF genes in garlic, thereby opening up new avenues for understanding their functions. The insights gained from this research provide a valuable resource for future investigations, particularly in the functional analysis of AsHSF genes.

The Brachypodium distachyon DREB transcription factor BdDREB-39 confers oxidative stress tolerance in transgenic tobacco

Key message BdDREB-39 is a DREB/CBF transcription factor, localized in the nucleus with transactivation activity, and BdDREB-39-overexpressing transgenic yeasts and tobacco enhanced the tolerance to oxidative stress.Abstract The DREB/CBF transcription factors are generally recognized to play an important factor in plant growth, development and response to various abiotic stresses. However, the mechanism of DREB/CBFs in oxidative stress response is largely unknown. This study isolated a DREB/CBF gene BdDREB-39 from Brachypodium distachyon (B. distachyon). Multiple sequence alignment and phylogenetic analysis showed that BdDREB-39 was closely related to the DREB proteins of oats, barley, wheat and rye and therefore its study can provide a reference for the excavation and genetic improvement of BdDREB-39 or its homologs in its closely related species. The transcript levels of BdDREB-39 were significantly up-regulated under H2O2 stress. BdDREB-39 was localised in the nucleus and functioned as a transcriptional activator. Overexpression of BdDREB-39 enhanced H2O2 tolerance in yeast. Transgenic tobaccos with BdDREB-39 had higher germination rates, longer root, better growth status, lesser reactive oxygen species (ROS) and malondialdehyde (MDA), and higher superoxide dismutase (SOD) and peroxidase (POD) activities than wild type (WT). The expression levels of ROS-related and stress-related genes were improved by BdDREB-39. In summary, these results revealed that BdDREB-39 can improve the viability of tobacco by regulating the expression of ROS and stress-related genes, allowing transgenic tobacco to accumulate lower levels of ROS and reducing the damage caused by ROS to cells. The BdDREB-39 gene has the potential for developing plant varieties tolerant to stress.

The cold-responsive C-repeat binding factors in Betula platyphylla Suk. positively regulate cold tolerance

Cold stress is one of the most destructive abiotic stresses limiting plant growth and development. CBF (C-repeat binding factor) transcription factors and their roles in cold response have been identified in Arabidopsis as well as several other plant species. However, the biological functions and related molecular mechanisms of CBFs in birch (Betula platyphylla Suk.) remain undetermined. In this study, five cold-responsive BpCBF genes, BpCBF1, BpCBF2, BpCBF7, BpCBF10 and BpCBF12 were cloned. Via protoplast transformation, BpCBF7 was found to be localized in nucleus. The result of yeast one hybrid assay validated the binding of BpCBF7 to the CRT/DRE (C-repeat/dehydration responsive element) elements in the promoter of BpERF1.1 gene. By overexpressing and repressing BpCBFs in birch plants, it was proven that BpCBFs play positive roles in the cold tolerance. At the metabolic level, BpCBFs OE lines had lower ROS accumulation, as well as higher activities of antioxidant enzymes (SOD, POD and CAT) and higher accumulation of protective substances (soluble sugar, soluble protein and proline). Via yeast one hybrid and co-transformation of effector and reporter vectors assay, it was proven that BpCBF7 can regulate the expression of BpERF5 and BpZAT10 genes by directly binding to their promoters. An interacting protein of BpCBF7, BpWRKY17, was identified by yeast two hybrid library sequencing and the interaction was validated with in vivo methods. These results indicates that BpCBFs can increase the cold tolerance of birch plants, partly by gene regulation and protein interaction. This study provides a reference for the research on CBF transcription factors and genetic improvement of forest trees upon abiotic stresses.

Deciphering the possible role of RNA-helicase genes mechanism in response to abiotic stresses in rapeseed (Brassica napus L.)

BACKGROUND

Plants mediate several defense mechanisms to withstand abiotic stresses. Several gene families respond to stress as well as multiple transcription factors to minimize abiotic stresses without minimizing their effects on performance potential. RNA helicase (RH) is one of the foremost critical gene families that can play an influential role in tolerating abiotic stresses in plants. However, little knowledge is present about this protein family in rapeseed (canola). Here, we performed a comprehensive survey analysis of the RH protein family in rapeseed (Brassica napus L.).

RESULTS

A total of 133 BnRHs genes have been discovered in this study. By phylogenetic analysis, RHs genes were divided into one main group and a subgroup. Examination of the chromosomal position of the identified genes showed that most of the genes (27%) were located on chromosome 3. All 133 identified sequences contained the main DEXDC domain, the HELICC domain, and a number of sub-domains. The results of biological process studies showed that about 17% of the proteins acted as RHs, 22% as ATP binding, and 14% as mRNA binding. Each part of the conserved motifs, communication network, and three-dimensional structure of the proteins were examined separately. The results showed that the RWC in leaf tissue decreased with higher levels of drought stress and in both root and leaf tissues sodium concentration was increased upon increased levels of salt stress treatments. The proline content were found to be increased in leaf and root with the increased level of stress treatment. Finally, the expression patterns of eight selected RHs genes that have been exposed to drought, salinity, cold, heat and cadmium stresses were investigated by qPCR. The results showed the effect of genes under stress. Examination of gene expression in the Hayola #4815 cultivar showed that all primers except primer #79 had less expression in both leaves and roots than the control level.

CONCLUSIONS

New finding from the study have been presented new insights for better understanding the function and possible mechanism of RH in response to abiotic stress in rapeseed.

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.

Genome-wide identification of 2-oxoglutarate and Fe (II)-dependent dioxygenase family genes and their expression profiling under drought and salt stress in potato

The 2-Oxoglutatrate-dependent dioxygenases (2OGDs) comprise the 2-Oxoglutatrate and Fe(II)-dependent dioxygenases (2ODD) enzyme families that facilitate the biosynthesis of various compounds like gibberellin, ethylene, etc. The 2OGDs are also involved in various catabolism pathways, such as auxin and salicylic acid catabolism. Despite their important roles, 2ODDs have not been studied in potato, which is the third most important crop globally. In this study, a comprehensive genome wide analysis was done to identify all 2ODDs in potatoes, and the putative genes were analysed for the presence of the signature 2OG-FeII_Oxy (PF03171) domain and the conserved DIOX_N (PF14226) domain. A total of 205 St2ODDs were identified and classified into eight groups based on their function. The physiochemical properties, gene structures, and motifs were analysed, and gene duplication events were also searched for St2ODDs. The active amino acid residues responsible for binding with 2-oxoglutarate and Fe (II) were conserved throughout the St2ODDs. The three-dimensional (3D) structures of the representative members of flavanol synthase (FNS), 1-aminocyclopropane-1-carboxylic acid oxidases (ACOs), and gibberellin oxidases (GAOXs) were made and docked with their respective substrates, and the potential interactions were visualised. The expression patterns of the St2ODDs under abiotic stressors such as heat, salt, and drought were also analysed. We found altered expression levels of St2ODDs under abiotic stress conditions, which was further confirmed for drought and salt stress using qRT-PCR. The expression levels of St2ODD115, St2ODD34, and St2ODD99 were found to be upregulated in drought stress with 2.2, 1.8, and 2.6 fold changes, respectively. After rewatering, the expression levels were normal. In salt stress, the expression levels of St2ODD151, St2ODD76, St2ODD91, and St2ODD34 were found to be upregulated after 24 hours (h), 48 hours (h), 72 hours (h), and 96 hours (h). Altogether, the elevated expression levels suggest the importance of St2ODDs under abiotic stresses, i.e., drought and salt. Overall, our study provided a knowledge base for the 2ODD gene family in potato, which can be used further to study the important roles of 2ODDs in potato plants.

Functional Characterization of AP2/ERF Transcription Factors during Flower Development and Anthocyanin Biosynthesis Related Candidate Genes in Lycoris

The APETALA2/ethylene-responsive transcription factor (AP2/ERF) family has been extensively investigated because of its significant involvement in plant development, growth, fruit ripening, metabolism, and plant stress responses. To date, there has been little investigation into how the AP2/ERF genes influence flower formation and anthocyanin biosynthesis in Lycoris. Herein, 80 putative LrAP2/ERF transcription factors (TFs) with complete open reading frames (ORFs) were retrieved from the Lycoris transcriptome sequence data, which could be divided into five subfamilies dependent on their complete protein sequences. Furthermore, our findings demonstrated that genes belonging to the same subfamily had structural similarities and conserved motifs. LrAP2/ERF genes were analyzed for playing an important role in plant growth, water deprivation, and flower formation by means of gene ontology (GO) enrichment analysis. The expression pattern of the LrAP2/ERF genes differed across tissues and might be important for Lycoris growth and flower development. In response to methyl jasmonate (MeJA) exposure and drought stress, the expression of each LrAP2/ERF gene varied across tissues and time. Moreover, a total of 20 anthocyanin components were characterized using ultra-performance liquid chromatography-electrospray ionization tandem mass spectrometry (UPLC-ESI-MS/MS) analysis, and pelargonidin-3-O-glucoside-5-O-arabinoside was identified as the major anthocyanin aglycone responsible for the coloration of the red petals in Lycoris. In addition, we mapped the relationships between genes and metabolites and found that LrAP2/ERF16 is strongly linked to pelargonidin accumulation in Lycoris petals. These findings provide the basic conceptual groundwork for future research into the molecular underpinnings and regulation mechanisms of AP2/ERF TFs in anthocyanin accumulation and Lycoris floral development.

Comprehensive investigation of cucumber heat shock proteins under abiotic stress conditions: A multi-omics survey

Heat-shock proteins (Hsps) are a family of proteins essential in preserving the vitality and functionality of proteins under stress conditions. Cucumber (Cucumis sativus) is a widely grown plant with high nutritional value and is used as a model organism in many studies. This study employed a genomics, transcriptomics, and metabolomics approach to investigate cucumbers' Hsps against abiotic stress conditions. Bioinformatics methods were used to identify six Hsp families in the cucumber genome and to characterize family members. Transcriptomics data from the Sequence Read Archive (SRA) database was also conducted to select CsHsp genes for further study. Real-time PCR was used to evaluate gene expression levels under different stress conditions, revealing that CssHsp-08 was a vital gene for resistance to stress conditions; including drought, salinity, cold, heat stresses, and ABA application. Gas Chromatography-Mass Spectrometry (GC-MS) analysis of plant extracts revealed that amino acids accumulate in leaves under high temperatures and roots under drought, while sucrose accumulates in both tissues under applied most stress factors. The study provides valuable insights into the structure, organization, evolution, and expression profiles of the Hsp family and contributes to a better understanding of plant stress mechanisms. These findings have important implications for developing crops that can withstand environmental stress conditions better.

Transcript-wide identification and expression pattern analysis to comprehend the roles of AP2/ERF genes under development and abiotic stress in Trichosanthes kirilowii

BACKGROUND

The APETALA 2/ ethylene-responsive element binding factors (AP2/ERF), are thought to be associated with plant abiotic stress response, and involved in some plant hormone signaling pathways. Trichosanthes kirilowii is an important edible and medicinal crop, so far no research has been conducted on the TkAP2/ERF genes.

RESULT

In this study, a total of 135 TkERFs were identified, these genes were divided into 4 subfamilies and clustered into 13 groups. Moreover, 37 paralogous pairs were identified, with only two having Ka/Ks values greater than 1, proving that most TkERF genes underwent purifying selection during evolution. Co-expression networks constructed using transcriptome data at various flowering stages revealed that 50, 64, and 67 AP2/ERF genes correlated with members of the ethylene, gibberellin, and abscisic acid signaling pathways, respectively. When tissue cultured seedlings were treated with ETH, GA3 and ABA, 11, 12 and 17 genes were found to be up-regulated, respectively, suggesting that some members of the TkERF gene family may be involved in plant hormone signaling pathways. And under 4 ℃, PEG and NaCl treatment, 15, 20 and 19 genes were up-regulated, respectively, this suggested that these selected genes might be involved in plant abiotic stresses.

CONCLUSIONS

Overall, we identified 135 AP2/ERF family members, a comprehensive analysis of AP2/ERF gene expression patterns by RNA-seq and qRT-PCR showed that they played important roles in flower development and abiotic stress. This study provided a theoretical basis for the functional study of TkAP2/ERF genes and the genetic improvement of T. kirilowii.

RNA-seq analysis and reconstruction of gene networks involved in response to salinity stress in quinoa (cv. Titicaca)

To better understand the mechanisms involved in salinity stress, the adaptability of quinoa cv. Titicaca-a halophytic plant-was investigated at the transcriptome level under saline and non-saline conditions. RNA-sequencing analysis of leaf tissue at the four-leaf stage by Illumina paired-end method was used to compare salt stress treatment (four days after stress at 13.8 dsm^-1) and control. Among the obtained 30,846,354 transcripts sequenced, 30,303 differentially expressed genes from the control and stress treatment samples were identified, with 3363 genes expressed ≥ 2 and false discovery rate (FDR) of < 0.001. Six differential expression genes were then selected and qRT-PCR was used to confirm the RNA-seq results. Some of the genes (Include; CML39, CBSX5, TRX1, GRXC9, SnRKγ1 and BAG6) and signaling pathways discussed in this paper not been previously studied in quinoa. Genes with ≥ 2 were used to design the gene interaction network using Cytoscape software, and AgriGO software and STRING database were used for gene ontology. The results led to the identification of 14 key genes involved in salt stress. The most effective hub genes involved in salt tolerance were the heat shock protein gene family. The transcription factors that showed a significant increase in expression under stress conditions mainly belonged to the WRKY, bZIP and MYB families. Ontology analysis of salt stress-responsive genes and hub genes revealed that metabolic pathways, binding, cellular processes and cellular anatomical entity are among the most effective processes involved in salt stress.

Genome-wide characterization of DcHsp90 gene family in carnation (Dianthus caryophyllus L.) and functional analysis of DcHsp90-6 in heat tolerance

Plant heat shock protein 90 (Hsp90) participates in various physiological processes including protein folding, degradation, and signal transduction. However, the DcHsp90 gene family in carnation (Dianthus caryophyllus L.) has not been systematically analyzed. We thoroughly examined and comprehensively analyzed the carnation DcHsp90 gene family in this study and discovered 9 DcHsp90 genes. Based on the phylogenetic examination, DcHsp90 proteins may be divided into two groups. DcHsp90 structural features were similar but varied between groups. Promoter analysis revealed the presence of many cis-acting elements, most of which were connected to growth and development, hormones, and stress. DcHsp90 genes may play distinct functions in heat stress response, according to gene expression analyses. The DcHsp90-6 was isolated, and its role in the reaction to heat stress was studied. Thermotolerance and superoxide dismutase activity in transgenic seedlings were enhanced by Arabidopsis overexpression of DcHsp90-6. After heat stress, transgenic plants' electrolyte leakage and malondialdehyde levels were much lower than wild-type plants. Furthermore, overexpression of DcHsp90-6 altered the expressions of stress-responsive genes such as AtHsp101, AtHsp90, AtGolS1, AtRS4/5, and AtHsfB1. This study provides comprehensive information on the DcHsp90 gene family and suggests that overexpressed DcHsp90-6 positively regulates thermotolerance highlighting the adaptation mechanism of carnation under heat stress.

Genome-Wide Analysis of DREB Family Genes and Characterization of Cold Stress Responses in the Woody Plant Prunus nana

Dehydration response element binding factor (DREB) is a family of plant-specific transcription factors, whose members participate in the regulation of plant responses to various abiotic stresses. Prunus nana, also known as the wild almond, is a member of the Rosaceae family that is rare and found to grow in the wild in China. These wild almond trees are found in hilly regions in northern Xinjiang, and exhibit greater drought and cold stress resistance than cultivated almond varieties. However, the response of P. nana DREBs (PnaDREBs) under low temperature stress is still unclear. In this study, 46 DREB genes were identified in the wild almond genome, with this number being slightly lower than that in the sweet almond (Prunus dulcis cultivar ‘Nonpareil’). These DREB genes in wild almond were separated into two classes. All PnaDREB genes were located on six chromosomes. PnaDREB proteins that were classified in the same groups contained specific shared motifs, and promoter analyses revealed that PnaDREB genes harbored a range of stress-responsive elements associated with drought, low-temperature stress, light responsivity, and hormone-responsive cis-regulatory elements within their promoter regions. MicroRNA target site prediction analyses also suggested that 79 miRNAs may regulate the expression of 40 of these PnaDREB genes, with PnaDREB2. To examine if these identified PnaDREB genes responded to low temperature stress, 15 of these genes were selected including seven homologous to Arabidopsis C-repeat binding factor (CBFs), and their expression was assessed following incubation for 2 h at 25 °C, 5 °C, 0 °C, −5 °C, or −10 °C. In summary, this analysis provides an overview of the P. nana PnaDREB gene family and provides a foundation for further studies of the ability of different PnaDREB genes to regulate cold stress responses in almond plants.

Genome-Wide Analysis of AP2/ERF Gene Superfamily in Ramie (Boehmeria nivea L.) Revealed Their Synergistic Roles in Regulating Abiotic Stress Resistance and Ramet Development

AP2/ERF transcription factors (TFs) are one of the largest superfamilies in plants, and play vital roles in growth and response to biotic/abiotic stresses. Although the AP2/ERF family has been extensively characterized in many species, very little is known about this family in ramie (Boehmeria nivea L.). In this study, 138 AP2/ERF TFs were identified from the ramie genome and were grouped into five subfamilies, including the AP2 (19), RAV (5), Soloist (1), ERF (77), and DREB (36). Unique motifs were found in the DREB/ERF subfamily members, implying significance to the AP2/ERF TF functions in these evolutionary branches. Segmental duplication events were found to play predominant roles in the BnAP2/ERF TF family expansion. Light-, stress-, and phytohormone-responsive elements were identified in the promoter region of BnAP2/ERF genes, with abscisic acid response elements (ABRE), methyl jasmonate response elements, and the dehydration response element (DRE) being dominant. The integrated transcriptome and quantitative real-time PCR (qPCR) revealed 12 key BnAP2/ERF genes positively responding to waterlogging. Five of the genes are also involved in ramet development, with two (BnERF-30 and BnERF-32) further showing multifunctional roles. The protein interaction prediction analysis further verified their crosstalk mechanism in coordinating waterlogging resistance and ramet development. Our study provides new insights into the presence of AP2/ERF TFs in ramie, and provides candidate AP2/ERF TFs for further studies on breeding varieties with coupling between water stress tolerance and high yield.

Expression Analysis of Molecular Chaperones Hsp70 and Hsp90 on Development and Metabolism of Different Organs and Testis in Cattle (Cattle-yak and Yak)

Hsp70 and Hsp90 play an important role in testis development and spermatogenesis regulation, but the exact connection between Hsp70 and Hsp90 and metabolic stress in cattle is unclear. Here, we focused on the male cattle−yak and yak, investigated the expression and localization of Hsp70 and Hsp90 in their tissues, and explored the influence of these factors on development and metabolism. In our study, a total of 54 cattle (24 cattle−yaks and 30 yaks; aged 1 day to 10 years) were examined. The Hsp90 mRNA of the cattle−yak was first cloned and compared with that of the yak, and variation in the amino acid sequence was found, which led to differences in protein spatial structure. Using real-time quantitative PCR (RT-qPCR) and Western blot (WB) techniques, we investigated whether the expression of Hsp70 and Hsp90 mRNA and protein are different in the cattle−yak and yak. We found a disparity in Hsp70 and Hsp90 mRNA and protein expression in different non-reproductive organs and in testicular tissues at different stages of development, while high expression was observed in the testes of both juveniles and adults. Moreover, it was intriguing to observe that Hsp70 expression was significantly high in the yak, whereas Hsp90 was high in the cattle−yak (p < 0.01). We also examined the location of Hsp70 and Hsp90 in the testis by immunohistochemical (IHC) and immunofluorescence (IF) techniques, and the results showed that Hsp70 and Hsp90 were positive in the epithelial cells, spermatogenic cells, and mesenchymal cells. In summary, our study proved that Hsp70 and Hsp90 expressions were different in different tissues (kidney, heart, cerebellum, liver, lung, spleen, and testis), and Hsp90 expression was high in the testis of the cattle−yak, suggesting that dysplasia of the cattle−yak may correlate with an over-metabolism of Hsp90.

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.

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.

Mining of Cloned Disease Resistance Gene Homologs (CDRHs) in Brassica Species and Arabidopsis thaliana

Simple Summary

Developing cultivars with resistance genes (R genes) is an effective strategy to support high yield and quality in Brassica crops. The availability of clone R gene and genomic sequences in Brassica species and Arabidopsis thaliana provide the opportunity to compare genomic regions and survey R genes across genomic databases. In this paper, we aim to identify genes related to cloned genes through sequence identity, providing a repertoire of species-wide related R genes in Brassica crops. The comprehensive list of candidate R genes can be used as a reference for functional analysis.

Abstract

Various diseases severely affect Brassica crops, leading to significant global yield losses and a reduction in crop quality. In this study, we used the complete protein sequences of 49 cloned resistance genes (R genes) that confer resistance to fungal and bacterial diseases known to impact species in the Brassicaceae family. Homology searches were carried out across Brassica napus, B. rapa, B. oleracea, B. nigra, B. juncea, B. carinata and Arabidopsis thaliana genomes. In total, 660 cloned disease R gene homologs (CDRHs) were identified across the seven species, including 431 resistance gene analogs (RGAs) (248 nucleotide binding site-leucine rich repeats (NLRs), 150 receptor-like protein kinases (RLKs) and 33 receptor-like proteins (RLPs)) and 229 non-RGAs. Based on the position and distribution of specific homologs in each of the species, we observed a total of 87 CDRH clusters composed of 36 NLR, 16 RLK and 3 RLP homogeneous clusters and 32 heterogeneous clusters. The CDRHs detected consistently across the seven species are candidates that can be investigated for broad-spectrum resistance, potentially providing resistance to multiple pathogens. The R genes identified in this study provide a novel resource for the future functional analysis and gene cloning of Brassicaceae R genes towards crop improvement.

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.

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

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

Genome-wide identification and expression profiles of AP2/ERF transcription factor family in mung bean (Vigna radiata L.)

Mung bean (Vigna radiata L. Wilczek) is an economically important grain legume crop in Asia, with high nutritional quality and potential in other parts of the world particularly arid and semiarid regions. Considering the potential adverse effects of drought, high salt, and other abiotic stresses on crop yield, significant efforts have been made to understand the underlying molecular mechanisms of tolerance to these abiotic stresses in legumes. In this study, a total of 186 putative AP2/ERF genes were identified, which were named VrERF1-186. These VrERF genes were classified into four main subfamilies according to the number of AP2 domains and sequence similarity, including 24 AP2 gene members, 81 ERF gene members, 79 DREB gene members, and 2 RAV members. VrERF genes are scattered across all 11 chromosomes and form small gene clusters on chromosomes due to segmental or tandem duplication. Promoter analysis revealed various cis-acting elements related to light, hormones, and stress responsiveness processes. The expression profiles of the VrERF genes in tissues during development and in response to abiotic stresses were assessed by transcriptome sequencing, and the selected reference genes were validated by qRT-PCR. A total of 174 VrERF genes were expressed in at least one of five tissues, while others showed distinct expression patterns in different tissues or under specific abiotic stress treatments, which indicates that VrERF genes are involved in developmental and environmental stress responses in V. radiata. In conclusion, the genome localization, genome-wide characterization, gene duplication, phylogenetic relationships, and expression pattern of VrERF genes in V. radiata were analyzed, and these results will lay the foundation for further functional analysis of these genes and improve stress tolerance to adverse conditions in plants.

Genome-wide investigation of AP2/ERF gene family in the desert legume Eremosparton songoricum: Identification, classification, evolution, and expression profiling under drought stress

Eremosparton songoricum (Litv.) Vass. is a rare leafless legume shrub endemic to central Asia which grows on bare sand. It shows extreme drought tolerance and is being developed as a model organism for investigating morphological, physiological, and molecular adaptations to harsh desert environments. APETALA2/Ethylene Responsive Factor (AP2/ERF) is a large plant transcription factor family that plays important roles in plant responses to various biotic and abiotic stresses and has been extensively studied in several plants. However, our knowledge on the AP2/ERF family in legume species is limited, and no respective study was conducted so far on the desert shrubby legume E. songoricum. Here, 153 AP2/ERF genes were identified based on the E. songoricum genome data. EsAP2/ERFs covered AP2 (24 genes), DREB (59 genes), ERF (68 genes), and Soloist (2 genes) subfamilies, and lacked canonical RAV subfamily genes based on the widely used classification method. The DREB and ERF subfamilies were further divided into A1-A6 and B1-B6 groups, respectively. Protein motifs and exon-intron structures of EsAP2/ERFs were also examined, which matched the subfamily/group classification. Cis-acting element analysis suggested that EsAP2/ERF genes shared many stress- and hormone-related cis-regulatory elements. Moreover, the gene numbers and the ratio of each subfamily and the intron-exon structures were systematically compared with other model plants ranging from algae to angiosperms, including ten legumes. Our results supported the view that AP2 and ERF evolved early and already existed in algae, whereas RAV and DREB began to appear in moss species. Almost all plant AP2 and Soloist genes contained introns, whereas most DREB and ERF genes did not. The majority of EsAP2/ERFs were induced by drought stress based on RNA-seq data, EsDREBs were highly induced and had the largest number of differentially expressed genes in response to drought. Eight out of twelve representative EsAP2/ERFs were significantly up-regulated as assessed by RT-qPCR. This study provides detailed insights into the classification, gene structure, motifs, chromosome distribution, and gene expression of AP2/ERF genes in E. songoricum and lays a foundation for better understanding of drought stress tolerance mechanisms in legume plants. Moreover, candidate genes for drought-resistant plant breeding are proposed.

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.

The key regulator LcERF056 enhances salt tolerance by modulating reactive oxygen species-related genes in Lotus corniculatus

BACKGROUND

The APETALA2/ethylene response factor (AP2/ERF) family are important regulatory factors involved in plants' response to environmental stimuli. However, their roles in salt tolerance in Lotus corniculatus remain unclear.

RESULTS

Here, the key salt-responsive transcription factor LcERF056 was cloned and characterised. LcERF056 belonging to the B3-1 (IX) subfamily of ERFs was considerably upregulated by salt treatment. LcERF056-fused GFP was exclusively localised to nuclei. Furthermore, LcERF056- overexpression (OE) transgenic Arabidopsis and L. corniculatus lines exhibited significantly high tolerance to salt treatment compared with wild-type (WT) or RNA interference expression (RNAi) transgenic lines at the phenotypic and physiological levels. Transcriptome analysis of OE, RNAi, and WT lines showed that LcERF056 regulated the downstream genes involved in several metabolic pathways. Chromatin immunoprecipitation-quantitative polymerase chain reaction (ChIP-qPCR) and yeast one-hybrid (Y1H) assay demonstrated that LcERF056 could bind to cis-element GCC box or DRE of reactive oxygen species (ROS)-related genes such as lipid-transfer protein, peroxidase and ribosomal protein.

CONCLUSION

Our results suggested that the key regulator LcERF056 plays important roles in salt tolerance in L. corniculatus by modulating ROS-related genes. Therefore, it may be a useful target for engineering salt-tolerant L. corniculatus or other crops.

Comprehensive RNAseq analysis for identification of genes expressed under heat stress in lentil

Lentils are highly sensitive to abrupt increases in temperature during the mid to late reproductive stages, leading to severe biomass and seed yield reduction. Therefore, we carried out an RNAseq analysis between IG4258 (heat tolerant) and IG3973 (heat sensitive) lentil genotypes at the reproductive stage under both normal and heat stress conditions in the field. It resulted in 209,549 assembled transcripts and among these 161,809 transcripts had coding regions, of which 94,437 transcripts were annotated. The differential gene expression analysis showed upregulation of 678 transcripts and downregulation of 680 transcripts between the tolerant and sensitive genotypes at the early reproductive stage. While 76 transcripts were upregulated and 47 transcripts were downregulated at the late reproductive stage under heat stress conditions. The validation of 12 up-or downregulated transcripts through RT-PCR corresponded well with the expression analysis data of RNAseq, with a correlation of R²  = 0.89. Among these transcripts, the DN364_c1_g1_i9 and DN2218_c0_g1_i5 transcripts encoded enzymes involved in the tryptophan pathway, indicating that tryptophan biosynthesis plays a role under heat stress in lentil. Moreover, KEGG pathways enrichment analysis identified transcripts associated with genes encoding proteins/regulating factors related to different metabolic pathways including signal transduction, fatty acid biosynthesis, rRNA processing, ribosome biogenesis, gibberellin (GA) biosynthesis, and riboflavin biosynthesis. This analysis also identified 6852 genic-SSRs leading to the development of 4968 SSR primers that are potential genomic resources for molecular mapping of heat-tolerant genes in lentil.

Chickpea tolerance to temperature stress: Status and opportunity for improvement

Chickpea is a globally important commercial crop and a key source of protein for vegetarian populations. Though chickpea was domesticated at least 3000 years ago, research into abiotic stress tolerance has been limited compared to cereal crops such as wheat. This review investigates the impacts of heat stress on chickpea, focusing on reproductive development. The fertilisation process is particularly sensitive to environmental stress, such as drought and heat that can reduce yields by up to 70%. Current research has largely focused on breeding cultivars that reach maturity faster to avoid stress rather than true thermotolerance and little is known of the impact of heat on cellular processes. This review suggests that there is ample variation within the chickpea gene pool for selective breeding to achieve improved abiotic stress tolerance. Rates of genetic progress will improve once key QTL are identified and the link between thermotolerance and pollen viability confirmed. Other benefits may arise from better understanding of heat shock proteins and molecular chaperones and their role in the protection of reproductive processes.

Genome-wide survey and identification of AP2/ERF genes involved in shoot and leaf development in Liriodendron chinense

BACKGROUND

Liriodendron chinense is a distinctive ornamental tree species due to its unique leaves and tulip-like flowers. The discovery of genes involved in leaf development and morphogenesis is critical for uncovering the underlying genetic basis of these traits. Genes in the AP2/ERF family are recognized as plant-specific transcription factors that contribute to plant growth, hormone-induced development, ethylene response factors, and stress responses.

RESULTS

In this study, we identified 104 putative AP2/ERF genes in the recently released L. chinense genome and transcriptome database. In addition, all 104 genes were grouped into four subfamilies, the AP2, ERF, RAV, and Soloist subfamilies. This classification was further supported by the results of gene structure and conserved motif analyses. Intriguingly, after application of a series test of cluster analysis, three AP2 genes, LcERF 94, LcERF 96, and LcERF 98, were identified as tissue-specific in buds based on the expression profiles of various tissues. These results were further validated via RT-qPCR assays and were highly consistent with the STC analysis. We further investigated the dynamic changes of immature leaves by dissecting fresh shoots into seven discontinuous periods, which were empirically identified as shoot apical meristem (SAM), leaf primordia and tender leaf developmental stages according to the anatomic structure. Subsequently, these three candidates were highly expressed in SAM and leaf primordia but rarely in tender leaves, indicating that they were mainly involved in early leaf development and morphogenesis. Moreover, these three genes displayed nuclear subcellular localizations through the transient transformation of tobacco epidermal cells.

CONCLUSIONS

Overall, we identified 104 AP2/ERF family members at the genome-wide level and discerned three candidate genes that might participate in the development and morphogenesis of the leaf primordium in L. chinense.

Transcriptome-based identification and expression profiling of AP2/ERF members in Caragana intermedia and functional analysis of CiDREB3

BACKGROUND

The AP2/ERF transcription factor family plays important roles in regulation of plant growth and development as well as the response of plants to stress. However, there are currently few studies focusing on the function of the AP2/ERF-type transcription factors in Caragana intermedia Kuang et H. C. Fu. Here, the expression pattern of AP2/ERF transcription factors family in different tissues and under four stress treatments were evaluated, and the function of CiDREB3 was examined.

METHODS AND RESULTS

In this study, the genes encoding the AP2/ERF family of transcription factors were screened from the C. intermedia drought transcriptome database and subjected to bioinformatic analysis using the online tool and software. The expression pattern of the members of AP2/ERF transcription factors in C. intermedia were detected via quantitative real-time PCR (qRT-PCR). The function of CiDREB3 on growth, development and drought tolerance was evaluated by transgenic Arabidopsis. As a result, 22 sequences with complete ORFs were obtained and all sequences were divided into 13 sub-groups. Most of the AP2/ERF transcription factors exhibited tissue-specific expression and were induced by cold, heat, NaCl and mannitol treatments. Furthermore, heterologous expression of CiDREB3 altered the morphology of the transgenic Arabidopsis thaliana L. Heynh and improved its drought tolerance during seedlings development.

CONCLUSIONS

Taken together, the results of the present study helped to better understand the function of the AP2/ERF family transcription factors in response to multiple abiotic stresses and uncovered the role of CiDREB3 in affecting the morphology and abiotic stress tolerance of Arabidopsis.

Investigation of Distinctive Morpho-Physio and Biochemical Alterations in Desi Chickpea at Seedling Stage Under Irrigation, Heat, and Combined Stress

Global climatic instabilities have become the main reason for drastic yield losses in chickpea. This shift in climate could be a great threat in the future for food security in developing countries. Chickpea production is badly hampered by heat stress coupled with drought stress, and these factors can reduce yields by 40-45%. To mitigate yield losses due these abiotic factors, irrigation supplementation could be the best strategy. The present study aimed to (i) investigate the tolerance response of 9 desi chickpea genotypes against heat stress (H), irrigation (I), and a combination of both (I+H) through morphophysiological and biochemical indices at early growth stage, and (ii) assess yield performance across multiple locations of the country. Results revealed that under irrigation treatment, all genotypes perform well, but the genotypes D-09027 and D-09013 showed best performance because, as compared to control, they retained root length, seedling fresh weight, root fresh weight, root dry weight, esterase activity, Malondialdehyde (MDA) content, total chlorophyll, and total carotenoids. Shoot length and total phenolic contents (TPC) increased in both genotypes. Superoxide dismutase (SOD) and peroxidase (POD) increased in D-09027 and retained in D-09013. Catalase activity increased in D-09013 and retained in D-09027. Protease activity, total water potential and osmotic potential decreased in both genotypes and depicted high yield potential with 27 and 30% increase in yield over Bhakhar-2011 (check), respectively. In case of heat stress, maximum tolerance was found in genotypes CH104/06 and D-09013 with no change in shoot and root length, seedling dry weight, shoot fresh and dry weight, root dry weight, relative water content, turgor water potential, catalase (CAT) activity, esterase activity, increased root fresh weight, peroxidase activity (POD), ascorbate peroxidase activity (APX), and lycopene with low accumulation of protease and Malondialdehyde content (MDA). Both genotypes depicted high yield potential with 30 and 43% increase in yield over check across multiple locations of the country. Under the combined treatment, most genotypes showed good performance, while CH104/06 was selected as best performer genotype because significant of its increased root fresh weight, lycopene content, chlorophyll b, total carotenoids, total chlorophyll, retained shoot length, root length, seedling fresh and dry weight, total water potential, osmotic potential, relative water content, peroxidase activity (POD), catalase, esterase, and its ascorbate peroxidase (APX) activity and total soluble proteins (TSP) showed highest yield potential with 43% increase over check. Identified best performing and tolerant genotypes can further be employed for breeding climate-smart chickpea genotypes for sustainable production under changing climate.

Genome-Wide Identification and Characterization of HSP90-RAR1-SGT1-Complex Members From Arachis Genomes and Their Responses to Biotic and Abiotic Stresses

The molecular chaperone complex HSP90-RAR1-SGT1 (HRS) plays important roles in both biotic and abiotic stress responses in plants. A previous study showed that wild peanut Arachis diogoi SGT1 (AdSGT1) could enhance disease resistance in transgenic tobacco and peanut. However, no systematic analysis of the HRS complex in Arachis has been conducted to date. In this study, a comprehensive analysis of the HRS complex were performed in Arachis. Nineteen HSP90, two RAR1 and six SGT1 genes were identified from the allotetraploid peanut Arachis hypogaea, a number close to the sum of those from the two wild diploid peanut species Arachis duranensis and Arachis ipaensis. According to phylogenetic and chromosomal location analyses, thirteen orthologous gene pairs from Arachis were identified, all of which except AhHSP90-A8, AhHSP90-B9, AdHSP90-9, and AiHSP90-9 were localized on the syntenic locus, and they shared similar exon-intron structures, conserved motifs and expression patterns. Phylogenetic analysis showed that HSP90 and RAR1 from dicot and monocot plants diverged into different clusters throughout their evolution. Chromosomal location analysis indicated that AdSGT1 (the orthologous gene of AhSGT1-B3 in this study) might provide resistance to leaf late spot disease dependent on the orthologous genes of AhHSP90-B10 and AhRAR1-B in the wild peanut A. diogoi. Several HRS genes exhibited tissue-specific expression patterns, which may reflect the sites where they perform functions. By exploring published RNA-seq data, we found that several HSP90 genes play major roles in both biotic and abiotic stress responses, especially salt and drought responses. Autoactivation assays showed that AhSGT1-B1 could not be used as bait for yeast two-hybrid (Y2H) library screening. AhRAR1 and AhSGT1 could strongly interact with each other and interact with AhHSP90-B8. The present study represents the first systematic analysis of HRS complex genes in Arachis and provides valuable information for functional analyses of HRS complex genes. This study also offers potential stress-resistant genes for peanut improvement.

Genome-wide analysis of microRNA156 and its targets, the genes encoding SQUAMOSA promoter-binding protein-like (SPL) transcription factors, in the grass family Poaceae

MicroRNAs (miRNAs) are endogenous small non-coding RNAs that play an important role in post-transcriptional gene regulation in plants and animals by targeting messenger RNAs (mRNAs) for cleavage or repressing translation of specific mRNAs. The first miRNA identified in plants, miRNA156 (miR156), targets the SQUAMOSA promoter-binding protein-like (SPL) transcription factors, which play critical roles in plant phase transition, flower and plant architecture, and fruit development. We identified multiple copies of MIR156 and SPL in the rice, Brachypodium, sorghum, maize, and foxtail millet genomes. Sequence and chromosomal synteny analysis showed that both MIR156s and SPLs are conserved across species in the grass family. Analysis of expression data of the SPLs in eleven juvenile and adult rice tissues revealed that four non-miR156-targeted genes were highly expressed and three miR156-targeted genes were only slightly expressed in all tissues/developmental stages. The remaining SPLs were highly expressed in the juvenile stage, but their expression was lower in the adult stage. It has been proposed that under strong selective pressure, non-miR156-targeted mRNA may be able to re-structure to form a miRNA-responsive element. In our analysis, some non-miR156-targeted SPLs (SPL5/8/10) had gene structure and gene expression patterns similar to those of miR156-targeted genes, suggesting that they could diversify into miR156-targeted genes. DNA methylation profiles of SPLs and MIR156s in different rice tissues showed diverse methylation patterns, and hypomethylation of non-CG sites was observed in rice endosperm. Our findings suggested that MIR156s and SPLs had different origination and evolutionary mechanisms: the SPLs appear to have resulted from vertical evolution, whereas MIR156s appear to have resulted from strong evolutionary selection on mature sequences.

Genome wide identification, classification and functional characterization of heat shock transcription factors in cultivated and ancestral cottons (Gossypium spp.)

Heat shock transcription factors (HSF) have been demonstrated to play a significant transcriptional regulatory role in plants and considered as an integral part of signal transduction pathways against environmental stresses especially heat stress. Despite of their importance, HSFs have not yet been identified and characterized in all cotton species. In this study, we report the identification of 42, 39, 67, and 79 non-redundant HSF genes from diploid cottons G. arboreum (A2) and G. raimondii (D5), and tetraploid cottons G. barbadense (AD2) and G. hirsutum (AD1) respectively. The chromosome localization of identified HSFs revealed their random distribution on all the 13 chromosomes of A and D genomes of cotton with few regions containing HSFs in clusters. The genes structure and conserved domain analysis revealed the family-specific conservation of intron/exon organization and conserved domains in HSFs. Various abiotic stress-related cis-regulatory elements were identified from the putative promoter regions of cotton HSFs suggesting their possible role in mediating abiotic stress tolerance. The combined phylogenetic analysis of all the cotton HSFs grouped them into three subfamilies; with 145 HSFs belong to class A, 85 to class B, and 17 to class C subfamily. Moreover, a detailed analysis of HSF gene family in four species of cotton elucidated the role of allopolyploid and hybridization during evolutionary cascade of allotetraploid cotton. Comparatively, existence of more orthologous genes in cotton species than Arabidopsis, advocated that polyploidization produced new cotton specific orthologous gene clusters. Phylogenetic, collinearity and multiple synteny analyses exhibited dispersed, segmental, proximal, and tandem gene duplication events in HSF gene family. Duplication of gene events suggests that HSF gene family of cotton evolution was under strong purifying selection. Expression analysis revealed that GarHSF04 were found to be actively involved in PEG and salinity tolerance in G. arboreum. GhiHSF14 upregulated in heat and downregulated in salinity whilst almost illustrated similar behavior under cold and PEG treatments and GhiHSF21 exhibited down regulation almost across all the stresses in G. hirsutum. Overwhelmingly, present study paves the way to better understand the evolution of cotton HSF TFs and lays a foundation for future investigation of HSFs in improving abiotic stress tolerance in cotton.

Major QTLs and Potential Candidate Genes for Heat Stress Tolerance Identified in Chickpea (Cicer arietinum L.)

In the context of climate change, heat stress during the reproductive stages of chickpea (Cicer arietinum L.) leads to significant yield losses. In order to identify the genomic regions responsible for heat stress tolerance, a recombinant inbred line population derived from DCP 92-3 (heat sensitive) and ICCV 92944 (heat tolerant) was genotyped using the genotyping-by-sequencing approach and evaluated for two consecutive years (2017 and 2018) under normal and late sown or heat stress environments. A high-density genetic map comprising 788 single-nucleotide polymorphism markers spanning 1,125 cM was constructed. Using composite interval mapping, a total of 77 QTLs (37 major and 40 minor) were identified for 12 of 13 traits. A genomic region on CaLG07 harbors quantitative trait loci (QTLs) explaining >30% phenotypic variation for days to pod initiation, 100 seed weight, and for nitrogen balance index explaining >10% PVE. In addition, we also reported for the first time major QTLs for proxy traits (physiological traits such as chlorophyll content, nitrogen balance index, normalized difference vegetative index, and cell membrane stability). Furthermore, 32 candidate genes in the QTL regions that encode the heat shock protein genes, heat shock transcription factors, are involved in flowering time regulation as well as pollen-specific genes. The major QTLs reported in this study, after validation, may be useful in molecular breeding for developing heat-tolerant superior lines or varieties.

Genome-Wide Identification and Evolution of Receptor-Like Kinases (RLKs) and Receptor like Proteins (RLPs) in Brassica juncea

Brassica juncea, an allotetraploid species, is an important germplasm resource for canola improvement, due to its many beneficial agronomic traits, such as heat and drought tolerance and blackleg resistance. Receptor-like kinase (RLK) and receptor-like protein (RLP) genes are two types of resistance gene analogues (RGA) that play important roles in plant innate immunity, stress response and various development processes. In this study, genome wide analysis of RLKs and RLPs is performed in B. juncea. In total, 493 RLKs (LysM-RLKs and LRR-RLKs) and 228 RLPs (LysM-RLPs and LRR-RLPs) are identified in the genome of B. juncea, using RGAugury. Only 13.54% RLKs and 11.79% RLPs are observed to be grouped within gene clusters. The majority of RLKs (90.17%) and RLPs (52.83%) are identified as duplicates, indicating that gene duplications significantly contribute to the expansion of RLK and RLP families. Comparative analysis between B. juncea and its progenitor species, B. rapa and B. nigra, indicate that 83.62% RLKs and 41.98% RLPs are conserved in B. juncea, and RLPs are likely to have a faster evolution than RLKs. This study provides a valuable resource for the identification and characterisation of candidate RLK and RLP genes.

Alternative Splicing Diversified the Heat Response and Evolutionary Strategy of Conserved Heat Shock Protein 90s in Hexaploid Wheat (Triticum aestivum L.)

Crops are challenged by the increasing high temperature. Heat shock protein 90 (HSP90), a molecular chaperone, plays a critical role in the heat response in plants. However, the evolutionary conservation and divergence of HSP90s homeologs in polyploidy crops are largely unknown. Using the newly released hexaploid wheat reference sequence, we identified 18 TaHSP90s that are evenly distributed as homeologous genes among three wheat subgenomes, and were highly conserved in terms of sequence identity and gene structure among homeologs. Intensive time-course transcriptomes showed uniform expression and transcriptional response profiles among the three TaHSP90 homeologs. Based on the comprehensive isoforms generated by combining full-length single-molecule sequencing and Illumina short read sequencing, 126 isoforms, including 90 newly identified isoforms of TaHSP90s, were identified, and each TaHSP90 generated one to three major isoforms. Intriguingly, the numbers and the splicing modes of the major isoforms generated by three TaHSP90 homeologs were obviously different. Furthermore, the quantified expression profiles of the major isoforms generated by three TaHSP90 homeologs are also distinctly varied, exhibiting differential alternative splicing (AS) responses of homeologs. Our results showed that the AS diversified the heat response of the conserved TaHSP90s and provided a new perspective for understanding about functional conservation and divergence of homologous genes.

Genome-wide analysis of HSP90 gene family in the Mediterranean olive (Olea europaea subsp. europaea) provides insight into structural patterns, evolution and functional diversity

Plants regularly experience multiple abiotic and biotic pressures affecting their normal development. The 90-kDa heat shock protein (HSP90) plays a dynamic role in countering abiotic and biotic stresses via a plethora of functional mechanisms. The HSP90 has been investigated in many plant species. However, there is little information available about this gene family in the cultivated Mediterranean olive tree, Olea europaea subsp. europaea var. europaea. In the current study, we systematically performed genome-wide identification and characterization of the HSP90 gene family in O. europaea var. europaea (OeHSP90s). Twelve regular OeHSP90s were identified, which were phylogenetically grouped into two major clusters and four sub-clusters, showing five paralogous gene pairs evolving under purifying selection. All of the 12 proteins contained a Histidine kinase-like ATPase (HATPase_c) domain, justifying the role played by HSP90 proteins in ATP binding and hydrolysis. The predicted 3D structure of OeHSP90 proteins provided information to understand their functions at the biochemical level. Consistent with their phylogenetic relationships, OeHSP90 members were predicted to be localized in different cellular compartments, suggesting their involvement in various subcellular processes. In consonance with their spatial organization, olive HSP90 family members were found to share similar motif arrangements and similar number of exons. We found that OeHSP90 promoters contained various cis-acting elements associated with light responsiveness, hormone signaling pathways and reaction to various stress conditions. In addition, expression sequence tags (ESTs) analysis offered a view of OeHSP90 tissue- and developmental stage specific pattern of expression. Proteins interacting with OeHSP90s were predicted and their potential roles were discussed. Overall, our results offer premises for further investigation of the implication of HSP90 genes in the physiological processes of the olive and its adaptation to stresses.

Identification and Characterization of Contrasting Genotypes/Cultivars for Developing Heat Tolerance in Agricultural Crops: Current Status and Prospects

Rising global temperatures due to climate change are affecting crop performance in several regions of the world. High temperatures affect plants at various organizational levels, primarily accelerating phenology to limit biomass production and shortening reproductive phase to curtail flower and fruit numbers, thus resulting in severe yield losses. Besides, heat stress also disrupts normal growth, development, cellular metabolism, and gene expression, which alters shoot and root structures, branching patterns, leaf surface and orientation, and anatomical, structural, and functional aspects of leaves and flowers. The reproductive growth stage is crucial in plants' life cycle, and susceptible to high temperatures, as reproductive processes are negatively impacted thus reducing crop yield. Genetic variation exists among genotypes of various crops to resist impacts of heat stress. Several screening studies have successfully phenotyped large populations of various crops to distinguish heat-tolerant and heat-sensitive genotypes using various traits, related to shoots (including leaves), flowers, fruits (pods, spikes, spikelets), and seeds (or grains), which have led to direct release of heat-tolerant cultivars in some cases (such as chickpea). In the present review, we discuss examples of contrasting genotypes for heat tolerance in different crops, involving many traits related to thermotolerance in leaves (membrane thermostability, photosynthetic efficiency, chlorophyll content, chlorophyll fluorescence, stomatal activity), flowers (pollen viability, pollen germination, fertilization, ovule viability), roots (architecture), biomolecules (antioxidants, osmolytes, phytohormones, heat-shock proteins, other stress proteins), and "omics" (phenomics, transcriptomics, genomics) approaches. The traits linked to heat tolerance can be introgressed into high yielding but heat-sensitive genotypes of crops to enhance their thermotolerance. Involving these traits will be useful for screening contrasting genotypes and would pave the way for characterizing the underlying molecular mechanisms, which could be valuable for engineering plants with enhanced thermotolerance. Wherever possible, we discussed breeding and biotechnological approaches for using these traits to develop heat-tolerant genotypes of various food crops.

Genome-wide analysis of tandem duplicated genes and their contribution to stress resistance in pigeonpea (Cajanus cajan)

Pigeonpea is the main protein source for more than one billion people, and it shows a strong adaptation to biotic stress and abiotic stress. Gene duplication is a fundamental process in genome evolution. Although the draft sequence of the pigeonpea genome has been available since 2011, further analysis of tandem duplicated genes (TDGs) and their contribution to the evolution of pigeonpea has not been reported. In this study, we identify 3211 TDGs in the pigeonpea genome and KEGG enrichment analysis of these genes shows that the TDGs are significantly enriched in resistance-related pathways. In addition, we find that TDGs are more abundant in retrotransposon-related genes in pigeonpea than in the other species included in our study. These results indicate that stress resistance in pigeonpea may be ascribed to resistance-related pathways and retrotransposons originating from tandem duplications. Our study will provide an important basis for further research in pigeonpea breeding.

Third-generation sequencing found LncRNA associated with heat shock protein response to heat stress in Populus qiongdaoensis seedlings

Background

As air temperatures increase globally, more and more plants are exposed to heat-stress conditions. Although many studies have explored regulation networks in plants with the aim of improving their heat-stress tolerance, only few have revealed them in trees. Here, individuals of Populus qiongdaoensis seedlings, which grows naturally in tropical areas, exposed to heat at 40 °C and the non-coding regulation networks were explored using the PacBio RSII and the Illumina sequencing platform.

Results

In total, we obtained 88,161 full-length transcripts representing 39,343 genes using 5,498,988 long reads and 350,026,252 clean reads, and also 216 microRNAs (miRNAs) via 95,794,107 reads. We then identified 928 putative long non-coding RNAs (lncRNAs), consisting of 828 sense lncRNAs (89.22%), 34 long intergenic non-coding RNAs (3.66%), 16 antisense (1.72%), and 50 sense intronic lncRNAs (5.39%). Under the dual criteria of |log₂fold-change| ≥ 1 and P-value < 0.05, 1690 genes, 25 lncRNAs, and 15 miRNAs were found differentially expressed under the heat stress treatment. Furthermore, 563 and 595 mRNAs were detected as target genes of 14 differently expressed miRNAs and 26 differentially expressed lncRNAs. Functional annotation analysis of these target genes demonstrated they were related to cell membrane stability, plant hormone signal transduction, antioxidation, and aldarate metabolism. Lastly, we uncovered a key interaction network of lncRNAs, miRNAs and mRNAs that consisted of miR1444d, miR482a.1, miR530a, lncHSP18.2, HSP18.1, and HSP18.2. Expression level analysis showed that miRNAs in the network were up-regulated, while mRNAs and lncRNA were down-regulated, and also found that lncHSP18.2 may cis-regulate HSP18.2.

Conclusions

Functional enrichment analysis of target genes of miRNAs and lncRNAs indicated that miRNAs and lncRNAs play an important role in the response to heat stress P. qiongdaoensis. Lastly, by investigating the miRNA–lncRNA–mRNA network of this species, we revealed that miRNAs may negatively regulate both lncRNAs and mRNAs in tree responses to heat stress, and found that lncHSP18.2 may cis-regulate HSP18.2.

Genome-Wide Identification and Functional Characterization of the Heat Shock Factor Family in Eggplant (Solanum melongena L.) under Abiotic Stress Conditions

Plant heat shock factors (Hsfs) play crucial roles in various environmental stress responses. Eggplant (Solanum melongena L.) is an agronomically important and thermophilic vegetable grown worldwide. Although the functions of Hsfs under environmental stress conditions have been characterized in the model plant Arabidopsis thaliana and tomato, their roles in responding to various stresses remain unclear in eggplant. Therefore, we characterized the eggplant SmeHsf family and surveyed expression profiles mediated by the SmeHsfs under various stress conditions. Here, using reported Hsfs from other species as queries to search SmeHsfs in the eggplant genome and confirming the typical conserved domains, we identified 20 SmeHsf genes. The SmeHsfs were further classified into 14 subgroups on the basis of their structure. Additionally, quantitative real-time PCR revealed that SmeHsfs responded to four stresses—cold, heat, salinity and drought—which indicated that SmeHsfs play crucial roles in improving tolerance to various abiotic stresses. The expression pattern of SmeHsfA6b exhibited the most immediate response to the various environmental stresses, except drought. The genome-wide identification and abiotic stress-responsive expression pattern analysis provide clues for further analysis of the roles and regulatory mechanism of SmeHsfs under environmental stresses.

Drought tolerance in Triticum aestivum L. genotypes associated with enhanced antioxidative protection and declined lipid peroxidation

Drought is one of the major constraints in wheat production and causes a huge loss at grain-filling stage. In this study we highlighted the response of different wheat genotypes under drought stress at the grain-filling stage. Field experiments were conducted to evaluate 72 wheat (Triticum aestivum L.) genotypes under two water regimes: irrigated and drought condition. Four wheat genotypes, two each of drought tolerant (IC36761A, IC128335) and drought-susceptible category (IC335732 and IC138852) were selected on the basis of agronomic traits and drought susceptibility index (DSI), to understand their morphological, biochemical and molecular basis of drought stress tolerance. Among agronomic traits, productive tillers followed by biomass had high percent reduction under drought stress, thus drought stress had a great impact. Antioxidant activity (AO), total phenolic and proline content were found to be significantly higher in IC128335 genotype. Differential expression pattern of transcription factors of ten genes revealed that transcription factor qTaWRKY2 followed by qTaDREB, qTaEXPB23 and qTaAPEX might be utilized for developing wheat varieties resistant to drought stress. Understanding the role of TFs would be helpful to decipher the molecular mechanism involved in drought stress. Identified genotypes (IC128335 and IC36761A) may be useful as parental material for future breeding program to generate new drought-tolerant varieties.

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.

Characterization of the AP2/ERF Transcription Factor Family and Expression Profiling of DREB Subfamily under Cold and Osmotic Stresses in Ammopiptanthus nanus

APETALA2/ethylene-responsive factor (AP2/ERF) is one of the largest transcription factor (TF) families in plants, which play important roles in regulating plant growth, development, and response to environmental stresses. Ammopiptanthus nanus, an unusual evergreen broad-leaved shrub in the arid region in the northern temperate zone, demonstrates a strong tolerance to low temperature and drought stresses, and AP2/ERF transcription factors may contribute to the stress tolerance of A. nanus. In the current study, 174 AP2/ERF family members were identified from the A. nanus genome, and they were divided into five subfamilies, including 92 ERF members, 55 dehydration-responsive element binding (DREB) members, 24 AP2 members, 2 RAV members, and 1 Soloist member. Compared with the other leguminous plants, A. nanus has more members of the DREB subfamily and the B1 group of the ERF subfamily, and gene expansion in the AP2/ERF family is primarily driven by tandem and segmental duplications. Promoter analysis showed that many stress-related cis-acting elements existed in promoter regions of the DREB genes, implying that MYB, ICE1, and WRKY transcription factors regulate the expression of DREB genes in A. nanus. Expression profiling revealed that the majority of DREB members were responsive to osmotic and cold stresses, and several DREB genes such as EVM0023336.1 and EVM0013392.1 were highly induced by cold stress, which may play important roles in cold response in A. nanus. This study provided important data for understanding the evolution and functions of AP2/ERF and DREB transcription factors in A. nanus.

The chaperone MeHSP90 recruits MeWRKY20 and MeCatalase1 to regulate drought stress resistance in cassava

The 90 kDa heat shock protein (HSP90) is widely involved in various developmental processes and stress responses in plants. However, the molecular chaperone HSP90-constructed protein complex and its function in cassava remain elusive. In this study, we report that HSP90 is essential for drought stress resistance in cassava by regulating abscisic acid (ABA) and hydrogen peroxide (H₂ O₂ ) using two specific protein inhibitors of HSP90 (geldanamycin (GDA) and radicicol (RAD)). Among 10 MeHSP90s, the transcript of MeHSP90.9 is largely induced during drought stress. Further investigation identifies MeWRKY20 and MeCatalase1 as MeHSP90.9-interacting proteins. MeHSP90.9-, MeWRKY20-, or MeCatalase1-silenced plants through virus-induced gene silencing display drought sensitivity in cassava, indicating that they are important to drought stress response. MeHSP90.9 can promote the direct transcriptional activation of MeWRKY20 on the W-box element of MeNCED5 promoter, encoding a key enzyme in ABA biosynthesis. Moreover, MeHSP90.9 positively regulates the activity of MeCatalase1, and MeHSP90.9-silenced cassava leaves accumulate more H₂ O₂ under drought stress. Taken together, we demonstrate that the MeHSP90.9 chaperone complex is a regulator of drought stress resistance in cassava.

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.

The expanding roles of APETALA2/Ethylene Responsive Factors and their potential applications in crop improvement

Understanding the molecular basis of the gene-regulatory networks underlying agronomic traits or plant responses to abiotic/biotic stresses is very important for crop improvement. In this context, transcription factors, which either singularly or in conjugation directly control the expression of many target genes, are suitable candidates for improving agronomic traits via genetic engineering. In this regard, members of one of the largest class of plant-specific APETALA2/Ethylene Response Factor (AP2/ERF) superfamily, which is implicated in various aspects of development and plant stress adaptation responses, are considered high-value targets for crop improvement. Besides their long-known regulatory roles in mediating plant responses to abiotic stresses such as drought and submergence, the novel roles of AP2/ERFs during fruit ripening or secondary metabolites production have also recently emerged. The astounding functional plasticity of AP2/ERF members is considered to be achieved by their interplay with other regulatory networks and signalling pathways. In this review, we have integrated the recently accumulated evidence from functional genomics studies and described their newly emerged functions in plants. The key structural features of AP2/ERF proteins and the modes of their action are briefly summarized. The importance of AP2/ERFs in plant development and stress responses and a summary of the event of their successful applications in crop improvement programs are also provided. Altogether, we envisage that the synthesized information presented in this review will be useful to design effective strategies for improving agronomic traits in crop plants.

Impact of Combined Heat and Drought Stress on the Potential Growth Responses of the Desert Grass Artemisia sieberi alba: Relation to Biochemical and Molecular Adaptation

Artemisia sieberi alba is one of the important plants frequently encountered by the combined effect of drought and heat stress. In the present study, we investigated the individual and combined effect of drought and heat stress on growth, photosynthesis, oxidative damage, and gene expression in A. sieberi alba. Drought and heat stress triggered oxidative damage by increasing the accumulation of hydrogen peroxide, and therefore electrolyte leakage. The accumulation of secondary metabolites, such as phenol and flavonoids, and proline, mannitol, inositol, and sorbitol, was increased due to drought and heat stress exposure. Photosynthetic attributes including chlorophyll synthesis, stomatal conductance, transpiration rate, photosynthetic efficiency, and chlorophyll fluorescence parameters were drastically reduced due to drought and heat stress exposure. Relative water content declined significantly in stressed plants, which was evident by the reduced leaf water potential and the water use efficiency, therefore, affecting the overall growth performance. Relative expression of aquaporin (AQP), dehydrin (DHN1), late embryogenesis abundant (LEA), osmotin (OSM-34), and heat shock proteins (HSP70) were significantly higher in stressed plants. Drought triggered the expression of AQP, DHN1, LEA, and OSM-34 more than heat, which improved the HSP70 transcript levels. A. sieberi alba responded to drought and heat stress by initiating key physio-biochemical and molecular responses, which were distinct in plants exposed to a combination of drought and heat stress.

A comprehensive analysis of the B3 superfamily identifies tissue-specific and stress-responsive genes in chickpea (Cicer arietinum L.)

The aim of this study was to provide a comprehensive analysis of the plant-specific B3 domain-containing transcription factors (TFs) in chickpea. Scanning of the chickpea genome resulted in the identification of 51 B3 domain-containing TFs that were located on seven out of eight chickpea chromosomes. Based on the presence of additional domains other than the B3 domain, the candidates were classified into four subfamilies, i.e., ARF (24), REM (19), LAV (6) and RAV (2). Phylogenetic analysis classified them into four groups in which members of the same group had similar intron-exon organization and motif composition. Genome duplication analysis of the candidate B3 genes revealed an event of segmental duplication that was instrumental in the expansion of the B3 gene family. Ka/Ks analysis showed that the B3 gene family was under purifying selection. Further, chickpea B3 genes showed maximum orthology with Medicago followed by soybean and Arabidopsis. Promoter analyses of the B3 genes led to the identification of several tissue-specific and stress-responsive cis-regulatory elements. Expression profiling of the candidate B3 genes using publicly available RNA-seq data of several chickpea tissues indicated their putative role in plant development and abiotic stress response. These findings were further validated by real-time expression analysis. Overall, this study provides a comprehensive analysis of the B3 domain-containing proteins in chickpea that would aid in devising strategies for crop manipulation in chickpea.