Transcriptome analysis of salt tolerant common bean (Phaseolus vulgaris L.) under saline conditions

Predicting the role of β-GAL genes in bean under abiotic stress and genome-wide characterization of β-GAL gene family members

Β-Gals are a subgroup of the glycoside hydrolase (GH) family of enzymes, which possess the Glyco_hydro_35 (GH35) domain. Although studies have been conducted on the β-Gal gene family in numerous plant species, no such research has been conducted on beans. The purpose of this study was to determine the gene expression levels of β-Gal genes in the leaf tissue of P. vulgaris under salt and drought stress using quantitative real-time polymerase chain reaction (qRT-PCR) and to perform a comprehensive analysis of β-Gal gene family members using bioinformatics tools. In the bean genome, 25 Pvul-βGAL proteins with amino acid numbers ranging from 291 to 1119, molecular weights from 32.94 to 126.56 kDa, and isoelectric points from 5.46 to 9.08 were identified. Both segmental and tandem duplication have occurred in β-Gal genes in the bean genome, and Pvul-BGAL genes have been subject to negative selection in the evolutionary process. For a deeper comprehension of the evolutionary proximity of Pvul-BGAL genes, a phylogenetic tree and synteny map were drawn together with Arabidopsis thaliana and Glycine max β-Gal genes. The expression profiles of β-Gal genes in different tissues of the bean were determined in silico. In addition, the expression profiles of β-Gal genes in the leaves of bean plants subjected to drought and salt stress were analyzed, and the role of β-Gal genes in salt and drought stress was estimated. In this study, the role of β-Gal gene family in abiotic stress response and the characterization of β-Gal genes in beans were determined for the first time and will provide a basis for future functional genomics studies.

Omics-Driven Strategies for Developing Saline-Smart Lentils: A Comprehensive Review

A number of consequences of climate change, notably salinity, put global food security at risk by impacting the development and production of lentils. Salinity-induced stress alters lentil genetics, resulting in severe developmental issues and eventual phenotypic damage. Lentils have evolved sophisticated signaling networks to combat salinity stress. Lentil genomics and transcriptomics have discovered key genes and pathways that play an important role in mitigating salinity stress. The development of saline-smart cultivars can be further revolutionized by implementing proteomics, metabolomics, miRNAomics, epigenomics, phenomics, ionomics, machine learning, and speed breeding approaches. All these cutting-edge approaches represent a viable path toward creating saline-tolerant lentil cultivars that can withstand climate change and meet the growing demand for high-quality food worldwide. The review emphasizes the gaps that must be filled for future food security in a changing climate while also highlighting the significant discoveries and insights made possible by omics and other state-of-the-art biotechnological techniques.

Understanding of Plant Salt Tolerance Mechanisms and Application to Molecular Breeding

Soil salinization is a widespread hindrance that endangers agricultural production and ecological security. High salt concentrations in saline soils are primarily caused by osmotic stress, ionic toxicity and oxidative stress, which have a negative impact on plant growth and development. In order to withstand salt stress, plants have developed a series of complicated physiological and molecular mechanisms, encompassing adaptive changes in the structure and function of various plant organs, as well as the intricate signal transduction networks enabling plants to survive in high-salinity environments. This review summarizes the recent advances in salt perception under different tissues, physiological responses and signaling regulations of plant tolerance to salt stress. We also examine the current knowledge of strategies for breeding salt-tolerant plants, including the applications of omics technologies and transgenic approaches, aiming to provide the basis for the cultivation of salt-tolerant crops through molecular breeding. Finally, future research on the application of wild germplasm resources and muti-omics technologies to discover new tolerant genes as well as investigation of crosstalk among plant hormone signaling pathways to uncover plant salt tolerance mechanisms are also discussed in this review.

Investigating the genetic basis of salt-tolerance in common bean: a genome-wide association study at the early vegetative stage

Salinity poses a significant challenge to global crop productivity, affecting approximately 20% of cultivated and 33% of irrigated farmland, and this issue is on the rise. Negative impact of salinity on plant development and metabolism leads to physiological and morphological alterations mainly due to high ion concentration in tissues and the reduced water and nutrients uptake. Common bean (Phaseolus vulgaris L.), a staple food crop accounting for a substantial portion of consumed grain legumes worldwide, is highly susceptible to salt stress resulting in noticeable reduction in dry matter gain in roots and shoots even at low salt concentrations. In this study we screened a common bean panel of diversity encompassing 192 homozygous genotypes for salt tolerance at seedling stage. Phenotypic data were leveraged to identify genomic regions involved in salt stress tolerance in the species through GWAS. We detected seven significant associations between shoot dry weight and SNP markers. The candidate genes, in linkage with the regions associated to salt tolerance or harbouring the detected SNP, showed strong homology with genes known to be involved in salt tolerance in Arabidopsis. Our findings provide valuable insights onto the genetic control of salt tolerance in common bean and represent a first contribution to address the challenge of salinity-induced yield losses in this species and poses the ground to eventually breed salt tolerant common bean varieties.

Characterization of Two-Component System gene (TCS) in melatonin-treated common bean under salt and drought stress

The two-component system (TCS) generally consists of three elements, namely the histidine kinase (HK), response regulator (RR), and histidine phosphotransfer (HP) gene families. This study aimed to assess the expression of TCS genes in P. vulgaris leaf tissue under salt and drought stress and perform a genome-wide analysis of TCS gene family members using bioinformatics methods. This study identified 67 PvTCS genes, including 10 PvHP, 38 PvRR, and 19 PvHK, in the bean genome. PvHK2 had the maximum number of amino acids with 1261, whilst PvHP8 had the lowest number with 87. In addition, their theoretical isoelectric points were between 4.56 (PvHP8) and 9.15 (PvPRR10). The majority of PvTCS genes are unstable. Phylogenetic analysis of TCS genes in A. thaliana, G. max, and bean found that PvTCS genes had close phylogenetic relationships with the genes of other plants. Segmental and tandem duplicate gene pairs were detected among the TCS genes and TCS genes have been subjected to purifying selection pressure in the evolutionary process. Furthermore, the TCS gene family, which has an important role in abiotic stress and hormonal responses in plants, was characterized for the first time in beans, and its expression of TCS genes in bean leaves under salt and drought stress was established using RNAseq and qRT-PCR analyses. The findings of this study will aid future functional and genomic studies by providing essential information about the members of the TCS gene family in beans.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01406-5.

Salt stress responses and alleviation strategies in legumes: a review of the current knowledge

Salinity is one of the most significant environmental factors limiting legumes development and productivity. Salt stress disturbs all developmental stages of legumes and affects their hormonal regulation, photosynthesis and biological nitrogen fixation, causing nutritional imbalance, plant growth inhibition and yield losses. At the molecular level, salt stress exposure involves large number of factors that are implicated in stress perception, transduction, and regulation of salt responsive genes' expression through the intervention of transcription factors. Along with the complex gene network, epigenetic regulation mediated by non-coding RNAs, and DNA methylation events are also involved in legumes' response to salinity. Different alleviation strategies can increase salt tolerance in legume plants. The most promising ones are Plant Growth Promoting Rhizobia, Arbuscular Mycorrhizal Fungi, seed and plant's priming. Genetic manipulation offers an effective approach for improving salt tolerance. In this review, we present a detailed overview of the adverse effect of salt stress on legumes and their molecular responses. We also provide an overview of various ameliorative strategies that have been implemented to mitigate/overcome the harmful effects of salt stress on legumes.

Genome-wide characterization of the common bean kinome: Catalog and insights into expression patterns and genetic organization

The protein kinase (PK) superfamily is one of the largest superfamilies in plants and is the core regulator of cellular signaling. Even considering this substantial importance, the kinome of common bean (Phaseolus vulgaris) has not been profiled yet. Here, we identified and characterised the complete set of kinases of common bean, performing an in-depth investigation with phylogenetic analyses and measurements of gene distribution, structural organization, protein properties, and expression patterns over a large set of RNA-Sequencing data. Being composed of 1,203 PKs distributed across all P. vulgaris chromosomes, this set represents 3.25% of all predicted proteins for the species. These PKs could be classified into 20 groups and 119 subfamilies, with a more pronounced abundance of subfamilies belonging to the receptor-like kinase (RLK)-Pelle group. In addition to provide a vast and rich reservoir of data, our study supplied insights into the compositional similarities between PK subfamilies, their evolutionary divergences, highly variable functional profile, structural diversity, and expression patterns, modeled with coexpression networks for investigating putative interactions associated with stress response.

Integrated transcriptomics and metabolomics analysis reveals key regulatory network that response to cold stress in common Bean (Phaseolus vulgaris L.)

Cold temperatures can be detrimental to crop survival and productivity. Breeding progress can be improved by understanding the molecular basis of low temperature tolerance. We investigated the key routes and critical metabolites related to low temperature resistance in cold-tolerant and -sensitive common bean cultivars 120 and 093, respectively. Many potential genes and metabolites implicated in major metabolic pathways during the chilling stress response were identified through transcriptomics and metabolomics research. Under chilling stress, the expression of many genes involved in lipid, amino acid, and flavonoid metabolism, as well as metabolite accumulation increased in the two bean types. Malondialdehyde (MDA) content was lower in 120 than in 093. Regarding amino acid metabolism, 120 had a higher concentration of acidic amino acids than 093, whereas 093 had a higher concentration of basic amino acids. Methionine accumulation was clearly higher in 120 than in 093. In addition, 120 had a higher concentration of many types of flavonoids than 093. Flavonoids, methionine and malondialdehyde could be used as biomarkers of plant chilling injury. Transcriptome analysis of hormone metabolism revealed considerably greater, expression of abscisic acid (ABA), gibberellin (GA), and jasmonic acid (JA) in 093 than in 120 during chilling stress, indicating that hormone regulation modes in 093 and 120 were different. Thus, chilling stress tolerance is different between 093 and 120 possibly due to transcriptional and metabolic regulation.

Characterization and Application of EST-SSR Markers Developed from Transcriptome Sequences in Elymus breviaristatus (Poaceae: Triticeae)

BACKGROUND

Elymus L. is the largest genus in the Triticeae tribe. Most species in this genus are highly stress resistant, with excellent forage value. Elymus breviaristatus, a rare species endemic to the Qinghai-Tibet Plateau (QTP), is declining due to habitat fragmentation. However, genetic data for E. breviaristatus are limited, with expressed sequence tag (EST) markers being particularly rare, hampering genetic studies and protection measures.

RESULTS

We obtained 9.06 Gb clean sequences from the transcriptome of E. breviaristatus, generating 171,522 unigenes, which were assembled and functionally annotated against five public databases. We identified 30,668 SSRs in the E. breviaristatus transcriptome, from which 103 EST-SSR primer pairs were randomly selected. Of these, 58 pairs of amplified products of the expected size, and 18 of the amplified products were polymorphic. Model-based Bayesian clustering, the unweighted pair group method with arithmetic average (UPGMA), and principal coordinate analysis (PCoA) of 179 wild E. breviaristatus in 12 populations using these EST-SSRs were generally consistent, grouping the 12 populations into two major clades. Analysis of molecular variance (AMOVA) found 70% of the genetic variation among the 12 populations and 30% within the populations, indicating a high level of genetic differentiation (or low gene exchange) among the 12 populations. The transferability of the 58 successful EST-SSR primers to 22 related hexaploid species was 86.2-98.3%. UPGMA analysis generally grouped species with similar genome types together.

CONCLUSIONS

Here, we developed EST-SSR markers from the transcriptome of E. breviaristatus. The transferability of these markers was evaluated, and the genetic structure and diversity of E. breviaristatus were explored. Our results provide a basis for the conservation and management of this endangered species, and the obtained molecular markers represent valuable resources for the exploration of genetic relationships among species in the Elymus genus.

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

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

Transcriptome analysis of halophyte Nitraria tangutorum reveals multiple mechanisms to enhance salt resistance

As a typical halophyte, Nitraria tangutorum Bobr. has attracted the interest of many researchers with the excellent salt tolerance. Elucidation of the mechanism of N. tangutorum salinity tolerance will facilitate the genetic improvement of productive plants faced with salinity. To reveal the molecular response to gradually accumulated salt stress in N. tangutorum, RNA-sequencing and analysis of gradually accumulated NaCl treated samples and control samples were performed, and a total of 1419 differentially expressed genes were identified, including 949 down-regulated genes and 470 up-regulated genes. Detailed analysis uncovered that the catabolism of organic compounds mainly based on oxidative phosphorylation genes was up-regulated. Additionally, various antioxidant genes, especially anthocyanin-related genes, were found to help N. tangutorum remove reactive oxygen species. Moreover, the Mitogen activated protein kinase signaling pathway and other signaling pathways co-regulated various salt tolerance activities. Additionally, intracellular ion homeostasis was maintained via regulation of osmotic regulator-related genes, cutin-related genes, and cell elongation-related genes to retain cellular water and reduce ion concentration. In particularly, simultaneous up-regulation in cytoskeleton-related genes, cell wall-related genes, and auxin-related genes, provided evidence of important role of cell expansion in plant salt tolerance. In conclusion, complex regulatory mechanisms modulated by multiple genes might contribute to the salt tolerance by N. tangutorum.

Physiological Analysis and Transcriptome Sequencing Reveal the Effects of Salt Stress on Banana (Musa acuminata cv. BD) Leaf

The salinization of soil is a widespread environmental problem. Banana (Musa acuminata L.) is a salt-sensitive plant whose growth, development, and production are constrained by salt stresses. However, the tolerance mechanism of this salt-sensitive banana to salt stress is still unclear. This study aimed to investigate the influence of NaCl treatment on phenotypic, physiological, and transcriptome changes in bananas. We found that the content of root activity, MDA, Pro, soluble sugar, soluble protein, and antioxidant enzymes activity in salt-stress treatment were significantly higher than the control in bananas. Transcriptome sequencing result identified an overall of 3,378 differentially expressed genes (DEGs) in banana leaves, and the Kyoto Encyclopedia of Genes and Genomes analysis indicated that these DEGs were involved in phenylpropanoid biosynthesis process, ribosome process, starch and sucrose metabolism, amino sugar process, and plant hormone signal transduction process that had simultaneously changed their expression under salt stress, which indicated these DEGs may play a role in promoting BD banana growth under salt treatments. The genes which were enriched in the phenylpropanoid biosynthesis process, starch and sucrose metabolism process, amino sugar process, and plant hormone signal transduction process were specifically regulated to respond to the salt stress treatments. Here, totally 48 differentially expressed transcription factors (TFs), including WRKY, MYB, NAC, and bHLH, were annotated in BD banana under salt stress. In the phenylpropane biosynthesis pathway, all transcripts encoding key enzymes were found to be significantly up-regulated, indicating that the genes in these pathways may play a significant function in the response of BD banana to salt stress. In conclusion, this study provides new insights into the mechanism of banana tolerance to salt stress, which provides a potential application for the genetic improvement of banana with salt tolerance.

Transcriptome Analysis Reveals Key Pathways and Candidate Genes Controlling Seed Development and Size in Ricebean (Vigna umbellata)

Ricebean (Vigna umbellata) is a lesser known pulse with well-recognized potential. Recently, it has emerged as a legume with endowed nutritional potential because of high concentration of quality protein and other vital nutrients in its seeds. However, the genes and pathways involved in regulating seed development and size are not understood in this crop. In our study, we analyzed the transcriptome of two genotypes with contrasting grain size (IC426787: large seeded and IC552985: small seeded) at two different time points, namely, 5 and 10 days post-anthesis (DPA). The bold seeded genotype across the time points (B5_B10) revealed 6,928 differentially expressed genes (DEGs), whereas the small seeded genotype across the time point (S5_S10) contributed to 14,544 DEGs. We have also identified several candidate genes for seed development-related traits like seed size and 100-seed weight. On the basis of similarity search and domain analysis, some candidate genes (PHO1, cytokinin dehydrogenase, A-type cytokinin, and ARR response negative regulator) related to 100-seed weight and seed size showed downregulation in the small seeded genotype. The MapMan and KEGG analysis confirmed that auxin and cytokinin pathways varied in both the contrasting genotypes and can therefore be the regulators of the seed size and other seed development-related traits in ricebeans. A total of 51 genes encoding SCF TIR1/AFB , Aux/IAA, ARFs, E3 ubiquitin transferase enzyme, and 26S proteasome showing distinct expression dynamics in bold and small genotypes were also identified. We have also validated randomly selected SSR markers in eight accessions of the Vigna species (V. umbellata: 6; Vigna radiata: 1; and Vigna mungo: 1). Cross-species transferability pattern of ricebean-derived SSR markers was higher in V. radiata (73.08%) than V. mungo (50%). To the best of our knowledge, this is the first transcriptomic study conducted in this crop to understand the molecular basis of any trait. It would provide us a comprehensive understanding of the complex transcriptome dynamics during the seed development and gene regulatory mechanism of the seed size determination in ricebeans.

The molecular events underpinning cultivar differences in melatonin counteracting salt damage in Phaseolus vulgaris

Melatonin (N -acetyl-5-methoxytryptamine) plays important roles in multiple stress responses, especially under salt stress. However, cultivar differences in melatonin mediated salt stress tolerance are unclear. Phaseolus vulgaris L. (common bean) cultivars Jiyin 1 (JY, salt-tolerant) and Xuliyabai (XL, salt-sensitive) were used in this study. Exogenous melatonin significantly improved root growth under salt stress in JY, but had little effect on XL. Physiology analysis showed significant differences in activities of antioxidant enzymes (superoxide, SOD; and catalase, CAT) and malondialdehyde content between JY and XL. Meanwhile, the change of ABA content in JY and XL root was opposite in salt plus melatonin treatment. Comparative root transcriptomes of JY and XL revealed 3505 and 668 differentially expressed genes (DEGs) regulated by salt stress and melatonin. The most enriched melatonin-responsive genes under salt stress are mainly involved in regulation of transcription, oxidation-reduction process, transcription factor activity, oxidoreductase activity. In addition, melatonin induced more obvious changes of DEGs in JY than that in XL under salt condition. Melatonin also significantly induced 41 DEGs only in JY, including signal transduction genes, transcription factors, ubiquitin protein ligases, ion homeostasis and osmotic adjustment genes etc. This study uncovered the molecular mechanism of cultivar difference of melatonin response under salt stress in common bean.

Transcriptome of Cenostigma pyramidale roots, a woody legume, under different salt stress times

Salinity stress has a significant impact on the gain of plant biomass. Our study provides the first root transcriptome of Cenostigma pyramidale, a tolerant woody legume from a tropical dry forest, under three different salt stress times (30 min, 2 h, and 11 days). The transcriptome was assembled using the RNA sequencing (RNA-Seq) de novo pipeline from GenPipes. We observed 932, 804, and 3157 upregulated differentially expressed genes (DEGs) and 164, 273, and 1332 downregulated DEGs for salt over 30 min, 2 h, and 11 days, respectively. For DEGs annotated with the Viridiplantae clade in the early stress periods, the response to salt stress was mainly achieved by stabilizing homeostasis of such ions like Na⁺ and K⁺ , signaling by Ca²⁺ , transcription factor modulation, water transport, and oxidative stress. For salt stress at 11 days, we observed a higher modulation of transcription factors including the WRKY, MYB, bHLH, NAC, HSF, and AP2-EREBP families, as well as DEGs involved in hormonal responses, water transport, sugar metabolism, proline, and reactive oxygen scavenging mechanisms. Five selected DEGs (K⁺ transporter, aquaporin, glutathione S-transferase, cyclic nucleotide-gated channel, and superoxide dismutase) were validated by qPCR. Our results indicated that C. pyramidale had an early perception of salt stress modulating ionic channels and transporters, and as the stress progressed, the focus turned to the antioxidant system, aquaporins, and complex hormone responses. The results of this first root transcriptome provide clues on how this native species modulate gene expression to achieve salt stress tolerance.

Transcriptome analysis and functional identification of GmMYB46 in soybean seedlings under salt stress

Salinity is one of the major abiotic stress that limits crop growth and productivity. We investigated the transcriptomes of salt-treated soybean seedlings versus a control using RNA-seq to better understand the molecular mechanisms of the soybean (Glycine max L.) response to salt stress. Transcriptome analysis revealed 1,235 differentially expressed genes (DEGs) under salt stress. Several important pathways and key candidate genes were identified by KEGG enrichment. A total of 116 differentially expressed transcription factors (TFs) were identified, and 17 TFs were found to belong to MYB families. Phylogenetic analysis revealed that these TFs may be involved in salt stress adaptation. Further analysis revealed that GmMYB46 was up-regulated by salt and mannitol and was localized in the nucleus. The salt tolerance of transgenic Arabidopsis overexpressing GmMYB46 was significantly enhanced compared to wild-type (WT). GmMYB46 activates the expression of salt stress response genes (P5CS1, SOD, POD, NCED3) in Arabidopsis under salt stress, indicating that the GmMYB46 protein mediates the salt stress response through complex regulatory mechanisms. This study provides information with which to better understand the molecular mechanism of salt tolerance in soybeans and to genetically improve the crop.

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

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

Alternative splicing of the dormancy-associated MADS-box transcription factor gene PpDAM1 is associated with flower bud dormancy in 'Dangshansu' pear (Pyrus pyrifolia white pear group)

Alternative splicing (AS) plays a crucial role in plant growth, development and response to various environmental changes. However, whether alternative splicing of MADS-box transcription factors contributes to the flower bud dormancy process in fruit trees still remains unknown. In this work, the AS profile of genes in the dormant flower buds of 'Dangshansu' pear tree were examined. A total number of 3661 alternatively spliced genes were identified, and three mRNA isoforms of the dormancy associated MADS box (DAM) gene, PpDAM1, derived by alternative splicing, designated as PpDAM1.1, PpDAM1.2 and PpDAM1.3, were characterized. Bimolecular fluorescence complementation (BiFC) analysis indicated that AS of PpDAM1 didn't affect the nucleus localization and homo-/heterodimerization of PpDAM1.1, PpDAM1.2 and PpDAM1.3 proteins, but disturbed the translocation of PpDAM1.1/PpDAM1.1, PpDAM1.3/PpDAM1.3, PpDAM1.1/PpDAM1.3, and PpDAM1.2/PpDAM1.3 dimers to the nucleus. Constitutive expression of PpDAM1.2, but not PpDAM1.1 and PpDAM1.3, in Arabidopsis retarded the growth and development of transgenic plants. Further comparative expression analyses of PpDAM1.1, PpDAM1.2 and PpDAM1.3 in the flower buds of 'Dangshansu' and a less dormant pear cultivar, 'Cuiguan', exhibited that the expression of all the three isoforms in 'Dangshansu' were significantly higher than in 'Cuiguan', especially PpDAM1.2, which showed a predominantly higher expression than PpDAM1.1 and PpDAM1.3 in both cultivars. Our results suggest that alternative splicing of PpDAM1 could play a crucial role in pear flower bud dormancy process.

Identification and Characterization of SRS Genes in Phaseolus vulgaris Genome and Their Responses Under Salt Stress

SHI-Related Sequence (SRS) transcription factors comprise a protein family with important roles in growth and development. However, the genome-wide study of the SRS protein family has not yet been carried out in the common bean. For this reason, the SRS family has been characterized in depth at both gene and protein levels and several bioinformatics methods have been used. As a result, 10 SRS genes have been identified and their proteins have been phylogenetically categorized into three major groups within the common bean. By investigating duplications that play a major role in the development of gene families, 19 duplication events have been identified in the SRS family (18 segmental and 1 tandem). In addition, using available RNAseq data, comparative expression analysis of Pvul-SRS genes was performed and expression changes in Pvul-SRS-1, 2, 4, 6, 7, and 10 genes were observed under both salt and drought stress. Five Pvul-SRS genes were selected based on RNAseq data (Pvul-SRS-1, 2, 4, 6, and 10) and screened with RT-qPCR in two common bean cultivars (Yakutiye 'salt-resistant' and Zulbiye 'salt-susceptible' cv.). These genes also showed different levels of expression between two common bean cultivars under salt stress conditions and this may explain the responses of Pvul-SRS genes against abiotic stress. In summary, this work is the first study in which in silico identification and characterization of Pvul-SRS genes have been examined at gene expression level. The results could therefore provide the basis for future studies of functional characterization of Pvul-SRS genes.

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.

Overexpression of A Biotic Stress-Inducible Pvgstu Gene Activates Early Protective Responses in Tobacco under Combined Heat and Drought

Drought and heat stresses are major factors limiting crop growth and productivity, and their effect is more devastating when occurring concurrently. Plant glutathione transferases (GSTs) are differentially expressed in response to different stimuli, conferring tolerance to a wide range of abiotic stresses. GSTs from drought-tolerant Phaseolus vulgaris var. "Plake Megalosperma Prespon" is expected to play an important role in the response mechanisms to combined and single heat and drought stresses. Herein, we examined wild-type N. tabacum plants (cv. Basmas Xanthi) and T1 transgenic lines overexpressing the stress-induced Pvgstu3-3 and Pvgstu2-2 genes. The overexpression of Pvgstu3-3 contributed to potential thermotolerance and greater plant performance under combined stress. Significant alterations in the primary metabolism were observed in the transgenic plants between combined stress and stress-free conditions. Stress-responsive differentially expressed genes (DEGs) and transcription factors (TFs) related to photosynthesis, signal transduction, starch and sucrose metabolism, osmotic adjustment and thermotolerance, were identified under combined stress. In contrast, induction of certain DEGs and TF families under stress-free conditions indicated that transgenic plants were in a primed state. The overexpression of the Pvgstu3-3 is playing a leading role in the production of signaling molecules, induction of specific metabolites and activation of the protective mechanisms for enhanced protection against combined abiotic stresses in tobacco.

Comparative Transcriptome, Metabolome, and Ionome Analysis of Two Contrasting Common Bean Genotypes in Saline Conditions

Soil salinity is a major abiotic stress factor that limits agricultural productivity worldwide, and this problem is expected to grow in the future. Common bean is an important protein source in developing countries however highly susceptible to salt stress. To understand the underlying mechanism of salt stress responses, transcriptomics, metabolomics, and ion content analysis were performed on both salt-tolerant and susceptible common bean genotypes in saline conditions. Transcriptomics has demonstrated increased photosynthesis in saline conditions for tolerant genotype while the susceptible genotype acted in contrast. Transcriptome also displayed active carbon and amino-acid metabolism for the tolerant genotype. Analysis of metabolites with GC-MS demonstrated the boosted carbohydrate metabolism in the tolerant genotype with increased sugar content as well as better amino-acid metabolism. Accumulation of lysine, valine, and isoleucine in the roots of the susceptible genotype suggested a halted stress response. According to ion content comparison, the tolerant genotype managed to block accumulation of Na⁺ in the leaves while accumulating significantly less Na⁺ in the roots compared to susceptible genotype. K⁺ levels increased in the leaves of both genotype and the roots of the susceptible one but dropped in the roots of the tolerant genotype. Additionally, Zn⁺² and Mn⁺² levels were dropped in the tolerant roots, while Mo⁺² levels were significantly higher in all tissues in both control and saline conditions for tolerant genotype. The results of the presented study have demonstrated the differences in contrasting genotypes and thus provide valuable information on the pivotal molecular mechanisms underlying salt tolerance.

Protein and gene integration analysis through proteome and transcriptome brings new insight into salt stress tolerance in pigeonpea (Cajanus cajan L.)

Salt stress is a major constrain to the productivity of nutritionally rich pigeonpea, an important legume of SE Asia and other parts of the world. The present study provides a comprehensive insight on integrated proteomic and transcriptomic analysis of root and shoot tissues of contrasting pigeonpea varieties (ICP1071- salt-sensitive; ICP7- salt-tolerant) to unravel salt stress induced pathways. Proteome analysis revealed 82 differentially expressed proteins (DEPs) with ≥±1.5 fold expression on 2-Dimensional (2D) gel. Of these, 25 DEPs identified through MALDI-TOF/TOF were classified using Uniprot software into functional categories. Pathways analyses using KAAS server showed the highest abundance of functional genes regulating metabolisms of carbohydrate followed by protein folding/degradation, amino acids and lipids. Expression studies on six genes (triosephosphate isomerase, oxygen evolving enhancer protein 1, phosphoribulokinase, cysteine synthase, oxygen evolving enhancer protein 2 and early nodulin like protein 2) with ≥±3 fold change were performed, and five of these showed consistency in transcript and protein expressions. Transcript analysis of root and shoot led to positive identification of 25 differentially expressed salt-responsive genes, with seven genes having ≥±5 fold change have diverse biological functions. Our combinatorial analysis suggests important role of these genes/proteins in providing salt tolerance in pigeonpea.

Comparative Transcriptome Analysis Reveals Molecular Defensive Mechanism of Arachis hypogaea in Response to Salt Stress

Abiotic stresses comprise all nonliving factors, such as soil salinity, drought, extreme temperatures, and metal toxicity, posing a serious threat to agriculture and affecting the plant production around the world. Peanut (Arachis hypogaea L.) is one of the most important crops for vegetable oil, proteins, minerals, and vitamins in the world. Therefore, it is of importance to understand the molecular mechanism of peanut against salt stress. Six transcriptome sequencing libraries including 24-hour salt treatments and control samples were constructed from the young leaves of peanut. A comprehensive analysis between two groups detected 3,425 differentially expressed genes (DEGs) including 2,013 upregulated genes and 1,412 downregulated genes. Of these DEGs, 141 transcription factors (TFs) mainly consisting of MYB, AP2/ERF, WRKY, bHLH, and HSF were identified in response to salinity stress. Further, GO categories of the DEGs highly related to regulation of cell growth, cell periphery, sustained external encapsulating structure, cell wall organization or biogenesis, antioxidant activity, and peroxidase activity were significantly enriched for upregulated DEGs. The function of downregulated DEGs was mainly enriched in regulation of metabolic processes, oxidoreductase activity, and catalytic activity. Fourteen DEGs with response to salt tolerance were validated by real-time PCR. Taken together, the identification of DEGs' response to salt tolerance of cultivated peanut will provide a solid foundation for improving salt-tolerant peanut genetic manipulation in the future.

Transcriptome Analysis of Salt-Sensitive and Tolerant Genotypes Reveals Salt-Tolerance Metabolic Pathways in Sugar Beet

Soil salinization is a common environmental problem that seriously affects the yield and quality of crops. Sugar beet (Beta vulgaris L.), one of the main sugar crops in the world, shows a strong tolerance to salt stress. To decipher the molecular mechanism of sugar beet under salt stress, we conducted transcriptomic analyses of two contrasting sugar beet genotypes. To the best of our knowledge, this is the first comparison of salt-response transcriptomes in sugar beet with contrasting genotypes. Compared to the salt-sensitive cultivar (S710), the salt-tolerant one (T710MU) showed better growth and exhibited a higher chlorophyll content, higher antioxidant enzyme activity, and increased levels of osmotic adjustment molecules. Based on a high-throughput experimental system, 1714 differentially expressed genes were identified in the leaves of the salt-sensitive genotype, and 2912 in the salt-tolerant one. Many of the differentially expressed genes were involved in stress and defense responses, metabolic processes, signal transduction, transport processes, and cell wall synthesis. Moreover, expression patterns of several genes differed between the two cultivars in response to salt stress, and several key pathways involved in determining the salt tolerance of sugar beet, were identified. Our results revealed the mechanism of salt tolerance in sugar beet and provided potential metabolic pathways and gene markers for growing salt-tolerant cultivars.

Functional annotation of mulberry (Morus spp.) transcriptome, differential expression of genes related to growth and identification of putative genic SSRs, SNPs and InDels

Growth is a complex trait associated with mulberry leaf yield and controlled by several genes. In this study, we have explored the molecular basis underlying growth using Transcriptome profiling of contrasting genotypes. A total of 66.6 Mbp of primary transcriptomes from high growth (HGG)-Jalalgarah-3 and M. laevigata (H) and, low growth genotypes (LGG)-Harmutty and Vadagaraparai-2; resulting in 24210, 27998, 28085 and 28764 final transcripts respectively. Out of the 34096 pooled transcripts, 20249 transcripts matched with at least one sequence of the non-redundant database. Functional annotation resulted in the categorization of 18970 transcripts into 3 gene ontology (GO) terms and 7440 were assigned to 23 Kyoto encyclopaedia of genes and genomes (KEGG) pathway. Based on the differentially expressed genes and gene enrichment analysis, over expression of photosynthetic related transcripts in HGG and defence related transcripts in LGG were noted. Simple sequence repeats were mined from unique transcripts and the most abundant motifs were tri- (1883) followed by di- (1710), tetra- (192), penta- (68) and hexa- (40) repeats. Further, a total of 390897 high quality SNPs and 8081 InDels were identified by mapping onto Morus notabilis reference genome. The study provides an insight into the expression of genes involved in growth and further research on utilization in gentic improvement of the crop.

Contrasting Effects of NaCl and NaHCO3 Stresses on Seed Germination, Seedling Growth, Photosynthesis, and Osmoregulators of the Common Bean (Phaseolus vulgaris L.)

The common bean (Phaseolus vulgaris L.), the most important food legume for human nutrition globally, contributes greatly to the improvement of soil fertility in semi-dry lands where most of the soil is already salinized or alkalized, such as in the Songnen Plain of China. In this study, we investigated the effects of salt stress (neutral and alkaline) on the salt-tolerant common bean. Seed germination, seedling growth, photosynthesis, and osmotic adjustment were assessed. Neutral and alkaline salt growth environments were simulated using NaCl and NaHCO3, respectively. The results indicated that at ≥60 mmol·L−1, both NaCl and NaHCO3 caused significant delays in seedling emergence and decreased seedling emergence rates. NaHCO3 stress suppressed seedling survival regardless of concentration; however, only NaCl concentrations >60 mmol·L−1 had the same effect. Alkaline salt stress remarkably suppressed photosynthesis and seedling establishment. The common bean compensated for the increase in inorganic anion concentration (influx of Na+) by synthesizing more organic acids and soluble sugars. This adaptive mechanism enabled the common bean to balance the large inflow of cations for maintaining a stable cell pH environment under alkaline salt stress.

Genome-wide identification of CAMTA gene family members in Phaseolus vulgaris L. and their expression profiling during salt stress

The calmodulin-binding transcriptional activator (CAMTA) family was first observed in tobacco (NtER1) during a screening for the CaM-binding proteins, which are known to be one of the fast response stress proteins. Due to the increased importance of plant transcription factors in recent years; genome-wide identification of CAMTA genes has been performed in several plant species, except for Phaseolus vulgaris. Therefore, our aim was to identify and characterize CAMTA genes in P. vulgaris via in silico genome-wide analysis approach. Our results showed a total of eight CAMTA genes that were identified and observed on five out of 11 chromosomes of P. vulgaris. Four gene couples were found to be segmentally-duplicated and these segmental duplication events were shown to occur from 29.97 to 92.06 MYA. The phylogenetic tree of CAMTA homologs from P. vulgaris, A. thaliana, and G. max. revealed three groups based on their homology and the intron numbers of Pvul-CAMTA genes, ranged from 11 to 12. According to the syteny analysis; CAMTA genes of P. vulgaris and G. max revealed higher similarity, because they have highly similar genomes compared to A. thaliana. All Pvul-CAMTA genes were targeted by miRNAs, which play a role in response mechanism of salt stress. To detect expression levels in different plant tissues, mRNA analysis of Pvul-CAMTA genes were performed using publicly available expression data in Phytozome v12.1. In addition, responses of Pvul-CAMTA genes to salt stress, were also examined via both RNAseq and qRT-PCR analysis. To identify and to obtain insight into biological functions of CAMTA genes in the genome of P. vulgaris, several analyses were conducted using many online and offline bioinformatic tools, genome databases and qRT-PCR analyses. Due to this study being the first in the identification of CAMTA genes in P. vulgaris, this study could be considered as an useful source for future CAMTA genes studies in either P. vulgaris or comparative different plant species.

Salinity stress response and 'omics' approaches for improving salinity stress tolerance in major grain legumes

Sustaining yield gains of grain legume crops under growing salt-stressed conditions demands a thorough understanding of plant salinity response and more efficient breeding techniques that effectively integrate modern omics knowledge. Grain legume crops are important to global food security being an affordable source of dietary protein and essential mineral nutrients to human population, especially in the developing countries. The global productivity of grain legume crops is severely challenged by the salinity stress particularly in the face of changing climates coupled with injudicious use of irrigation water and improper agricultural land management. Plants adapt to sustain under salinity-challenged conditions through evoking complex molecular mechanisms. Elucidating the underlying complex mechanisms remains pivotal to our knowledge about plant salinity response. Improving salinity tolerance of plants demand enriching cultivated gene pool of grain legume crops through capitalizing on 'adaptive traits' that contribute to salinity stress tolerance. Here, we review the current progress in understanding the genetic makeup of salinity tolerance and highlight the role of germplasm resources and omics advances in improving salt tolerance of grain legumes. In parallel, scope of next generation phenotyping platforms that efficiently bridge the phenotyping-genotyping gap and latest research advances including epigenetics is also discussed in context to salt stress tolerance. Breeding salt-tolerant cultivars of grain legumes will require an integrated "omics-assisted" approach enabling accelerated improvement of salt-tolerance traits in crop breeding programs.

Physiological, epigenetic and genetic regulation in some olive cultivars under salt stress

Cultivated olive, a typical fruit crop species of the semi-arid regions, could successfully face the new scenarios driven by the climate change through the selection of tolerant varieties to salt and drought stresses. In the present work, multidisciplinary approaches, including physiological, epigenetic and genetic studies, have been applied to clarify the salt tolerance mechanisms in olive. Four varieties (Koroneiki, Royal de Cazorla, Arbequina and Picual) and a related form (O. europaea subsp. cuspidata) were grown in a hydroponic system under different salt concentrations from zero to 200 mM. In order to verify the plant response under salt stress, photosynthesis, gas exchange and relative water content were measured at different time points, whereas chlorophyll and leaf concentration of Na+, K+ and Ca2+ ions, were quantified at 43 and 60 days after treatment, when stress symptoms became prominent. Methylation sensitive amplification polymorphism (MSAP) technique was used to assess the effects of salt stress on plant DNA methylation. Several fragments resulted differentially methylated among genotypes, treatments and time points. Real time quantitative PCR (RT-qPCR) analysis revealed significant expression changes related to plant response to salinity. Four genes (OePIP1.1, OePetD, OePI4Kg4 and OeXyla) were identified, as well as multiple retrotransposon elements usually targeted by methylation under stress conditions.

Analysis of the common bean (Phaseolus vulgaris L.) transcriptome regarding efficiency of phosphorus use

Common bean is a highly important food in tropical regions, where most production occurs on small farms with limited use of technology and, consequently, greater vulnerability to abiotic stresses such as nutritional stress. Usually phosphorus (P) is the most limiting nutrient for crop growth in these regions. The aim of this study was to characterize the gene expression profiles of the genotypes of common bean IAC Imperador (P-responsive) and DOR 364 (P-unresponsive) under different P concentrations using RNA-seq transcriptome sequencing technology. Plants were grown hydroponically, with application of two P concentrations (4.00 mg L^-1 restrictive level and 8.00 mg L^-1 control level). Differential expression analyses, annotation, and functional classification were performed comparing genotypes within each P rate administered and comparing each genotype response to the different P levels. Considering differential expression analyses within genotypes, IAC Imperador exhibited 1538 up-regulated genes under P restriction and 1679 up-regulated genes in the control, while DOR 364 exhibited 13 up-regulated genes in the control and only 2 up-regulated genes under P restriction, strongly corroborating P-unresponsiveness of this genotype. Genes related to phosphorus restriction were identified among the differentially expressed genes, including transcription factors such as WRKY, ERF, and MYB families, phosphatase related genes such as pyrophosphatase, acid phosphatase, and purple acid phosphatase, and phosphate transporters. The enrichment test for the P restriction treatment showed 123 enriched gene ontologies (GO) for IAC Imperador, while DOR 364 enriched only 24. Also, the enriched GO correlated with P metabolism, compound metabolic processes containing phosphate, nucleoside phosphate binding, phosphorylation, and also response to stresses. Thus, this study proved to be informative to phosphorus limitation in common bean showing global changes at transcript level.

Transcriptome profiling of a beach-adapted wild legume for dissecting novel mechanisms of salinity tolerance

Strophostyles helvola is a close relative to common bean (Phaseolus vulgaris) and inhabits both coastal and non-coastal regions in North America. However, the mechanism of saline adaptation in S. helvola remains unclear. A transcriptome profiling would facilitate dissecting the underlying molecular mechanisms in salinity-adapted S. helvola. In this study, we reported the RNA-seq analyses of two genotypes (a salt-tolerant beach genotype and a salt-sensitive inland genotype) of S. helvola stressed with salt. S. helvola plants were grown in pots and treated with half lethal-guided dose of NaCl solution for 3 h, 24 h, and 7d. The plants supplied with the same amount of water were used as controls. The whole roots sampled from the three time points were equally pooled as one biological replicate, and three replicates were used for library construction and transcriptome sequencing on Illumina Hiseq 2500. The comparative analyses of root transcriptomes presented here provides a valuable resource for discovery of genes and networks involved in salt tolerance in S. helvola.

Legume genetic resources and transcriptome dynamics under abiotic stress conditions

Grain legumes are an important source of nutrition and income for billions of consumers and farmers around the world. However, the low productivity of new legume varieties, due to the limited genetic diversity available for legume breeding programmes and poor policymaker support, combined with an increasingly unpredictable global climate is resulting in a large gap between current yields and the increasing demand for legumes as food. Hence, there is a need for novel approaches to develop new high-yielding legume cultivars that are able to cope with a range of environmental stressors. Next-generation technologies are providing the tools that could enable the more rapid and cost-effective genomic and transcriptomic studies for most major crops, allowing the identification of key functional and regulatory genes involved in abiotic stress resistance. In this review, we provide an overview of the recent achievements regarding abiotic stress resistance in a wide range of legume crops and highlight the transcriptomic and miRNA approaches that have been used. In addition, we critically evaluate the availability and importance of legume genetic resources with desirable abiotic stress resistance traits.

Transcriptome - Scale characterization of salt responsive bean TCP transcription factors

TEOSINTE-BRANCHED1/CYCLOIDEA/PCF (TCP) proteins are important regulators of growth and developmental processes including branching, floral organ morphogenesis and leaf growth as well as stress response. This study identified 27 TCP genes of Phaseolus vulgaris (common bean), which were divided into three clusters based on phylogenetic relationship. In addition, this study showed that some of TCP genes such as Pvul-TCP-4 and Pvul-TCP-15 located on chromosomes 3 and 7, Pvul-TCP-7 and Pvul-TCP-20 located on chromosome 7 and 9, were segmentally duplicated. On the other hand, a total of 20 Pvul-TCP genes have predicted to be targeted by microRNAs (miRNA). Most of the miRNA-target genes were Pvul-TCP-1, -11, -13 and -27, which were targeted by 13, 17, 22 and 13 plant miRNAs, respectively. miR319 was one of the highly represented regulatory miRNAs to target TCP transcripts. Promoter region analysis of TCP genes resulted that the GT-1 motif, which was related to salt stress, was found in 14 different Pvul-TCP genes. Expression profiling of 10 Pvul-TCP genes based on RNA-sequencing data further confirmed with quantitative real-time RT-PCR measurements identified that Pvul-TCP genes under salt stress are expressed in a cultivar- and tissue-specific manner.

Genome-wide analysis of Phaseolus vulgaris C2C2-YABBY transcription factors under salt stress conditions

The aim of this study was to identify and characterize the C2C2-YABBY family of genes by a genome-wide scale in common bean. Various in silico approaches were used for the study and the results were confirmed through common molecular biology techniques. Quantitative real-time PCR (qPCR) analysis was performed for identified putative PvulYABBY genes in leaf and root tissues of two common bean cultivars, namely Yakutiye and Zulbiye under salt stress condition. Eight candidate PvulYABBY proteins were discovered and the length of these proteins ranged from 173 to 256 amino acids. The isoelectric points (pIs) of YABBY proteins were between 5.18 and 9.34 and ranged from acidic to alkaline, and the molecular weight of PvulYABBYs were between 18978.4 and 28916.8 Da. Three segmentally duplicated gene couples among the identified eight PvulYABBY genes were detected. These segmentally duplicated gene couples were PvulYABBY-1/PvulYABBY-3, PvulYABBY-5/PvulYABBY-7 and PvulYABBY-6/PvulYABBY-8. The predicted number of exons among the PvulYABBY genes varied from 6 to 8 exons. Additionally, all genes found included introns within ORFs. PvulYABBY-2, -4, -5 and -7 genes were targeted by miRNAs of five plant species and a total of five miRNA families (miR5660, miR1157, miR5769, miR5286 and miR8120) were detected. According to RNA-seq analysis, all genes were up- or down-regulated except for PvulYABBY-1 and PvulYABBY-6 after salt stress treatment in leaf and root tissues of common bean. According to the qPCR analysis, six out of eight genes were expressed in the leaves but only four out of eight genes were expressed in the roots and these genes exhibited tissue- and cultivar-specific expression patterns.

Genome-wide expression profiling in leaves and roots of date palm (Phoenix dactylifera L.) exposed to salinity

BACKGROUND

Date palm, as one of the most important fruit crops in North African and West Asian countries including Oman, is facing serious growth problems due to salinity, arising from persistent use of saline water for irrigation. Although date palm is a relatively salt-tolerant plant species, its adaptive mechanisms to salt stress are largely unknown.

RESULTS

In order to get an insight into molecular mechanisms of salt tolerance, RNA was profiled in leaves and roots of date palm seedlings subjected to NaCl for 10 days. Under salt stress, photosynthetic parameters were differentially affected; all gas exchange parameters were decreased but the quantum yield of PSII was unaffected while non-photochemical quenching was increased. Analyses of gene expression profiles revealed 2630 and 4687 genes were differentially expressed in leaves and roots, respectively, under salt stress. Of these, 194 genes were identified as commonly responding in both the tissue sources. Gene ontology (GO) analysis in leaves revealed enrichment of transcripts involved in metabolic pathways including photosynthesis, sucrose and starch metabolism, and oxidative phosphorylation, while in roots genes involved in membrane transport, phenylpropanoid biosynthesis, purine, thiamine, and tryptophan metabolism, and casparian strip development were enriched. Differentially expressed genes (DEGs) common to both tissues included the auxin responsive gene, GH3, a putative potassium transporter 8 and vacuolar membrane proton pump.

CONCLUSIONS

Leaf and root tissues respond differentially to salinity stress and this study has revealed genes and pathways that are associated with responses to elevated NaCl levels and thus may play important roles in salt tolerance providing a foundation for functional characterization of salt stress-responsive genes in the date palm.

Exploring drought stress-regulated genes in senna (Cassia angustifolia Vahl.): a transcriptomic approach

De novo assembly of reads produced by next-generation sequencing (NGS) technologies offers a rapid approach to obtain expressed gene sequences for non-model organisms. Senna (Cassia angustifolia Vahl.) is a drought-tolerant annual undershrub of Caesalpiniaceae, a subfamily of Fabaceae. There are insufficient transcriptomic and genomic data in public databases for understanding the molecular mechanism underlying the drought tolerance of senna. Therefore, the main purpose of this study was to know the transcriptome profile of senna, with special reference to drought stress. RNA from two different stages of leaf development was extracted and sequenced separately using the Illumina technology. A total of 200 million reads were generated, and a de novo assembly of processed reads in the pooled transcriptome using Trinity yielded 43,413 transcripts which were further annotated using NCBI BLAST with "green plant database (txid 33090)," Swiss Prot, Kyoto Encyclopedia of Genes and Genomes (KEGG), Clusters of Orthologous Groups (COG), and Gene Ontology (GO). Out of the total transcripts, 42,280 (95.0 %) were annotated by BLASTX against the green plant database of NCBI. Senna transcriptome showed the highest similarity to Glycine max (41 %), followed by Phaseolus vulgaris (16 %), Cicer arietinum (15 %), and Medicago trancatula (5 %). The highest number of GO terms were enriched for the molecular functions category; of these "catalytic activity" (GO: 0003824) (25.10 %) and "binding activity" (GO: 0005488) (20.10 %) were most abundantly represented. We used InterProscan to see protein similarity at domain level; a total of 33,256 transcripts were annotated against the Pfam domains. The transcripts were assigned with various KEGG pathways. Coding DNA sequences (CDS) encoding various drought stress-regulated pathways such as signaling factors, protein-modifying/degrading enzymes, biosynthesis of phytohormone, phytohormone signaling, osmotically active compounds, free radical scavengers, chlorophyll metabolism, leaf cuticular wax, polyamines, and protective proteins were identified through BLASTX search. The lucine-rich repeat kinase family was the most abundantly found group of protein kinases. Orphan, bHLH, and bZIP family TFs were the most abundantly found in senna. Six genes encoding MYC2 transcription factor, 9-cis-epoxycarotenoid dioxygenase (NCED), l -ascorbate peroxidase (APX), aminocyclopropane carboxylate oxidase (ACO), abscisic acid 8'-hydroxylase (ABA), and WRKY transcription factor were confirmed through reverse transcriptase-PCR (RT-PCR) and Sanger sequencing for the first time in senna. The potential drought stress-related transcripts identified in this study provide a good start for further investigation into the drought adaptation in senna. Additionally, our transcriptome sequences are the valuable resource for accelerated genomics-assisted genetic improvement programs and facilitate manipulation of biochemical pathways for developing drought-tolerant genotypes of crop plants.

Transcriptome Analysis of the Response to NaCl in Suaeda maritima Provides an Insight into Salt Tolerance Mechanisms in Halophytes

Although salt tolerance is a feature representative of halophytes, most studies on this topic in plants have been conducted on glycophytes. Transcriptome profiles are also available for only a limited number of halophytes. Hence, the present study was conducted to understand the molecular basis of salt tolerance through the transcriptome profiling of the halophyte Suaeda maritima, which is an emerging plant model for research on salt tolerance. Illumina sequencing revealed 72,588 clustered transcripts, including 27,434 that were annotated using BLASTX. Salt application resulted in the 2-fold or greater upregulation of 647 genes and downregulation of 735 genes. Of these, 391 proteins were homologous to proteins in the COGs (cluster of orthologous groups) database, and the majorities were grouped into the poorly characterized category. Approximately 50% of the genes assigned to MapMan pathways showed homology to S. maritima. The majority of such genes represented transcription factors. Several genes also contributed to cell wall and carbohydrate metabolism, ion relation, redox responses and G protein, phosphoinositide and hormone signaling. Real-time PCR was used to validate the results of the deep sequencing for the most of the genes. This study demonstrates the expression of protein kinase C, the target of diacylglycerol in phosphoinositide signaling, for the first time in plants. This study further reveals that the biochemical and molecular responses occurring at several levels are associated with salt tolerance in S. maritima. At the structural level, adaptations to high salinity levels include the remodeling of cell walls and the modification of membrane lipids. At the cellular level, the accumulation of glycinebetaine and the sequestration and exclusion of Na+ appear to be important. Moreover, this study also shows that the processes related to salt tolerance might be highly complex, as reflected by the salt-induced enhancement of transcription factor expression, including hormone-responsive factors, and that this process might be initially triggered by G protein and phosphoinositide signaling.

Genome-wide identification of salinity responsive HSP70s in common bean

The present study is aimed to identify and characterize HSP70 (PvHSP70) genes in two different common bean cultivars under salt stress. For this purpose various in silico methods such as RNAseq data and qRT-PCR analysis were used. A total of 24 candidate PvHSP70 gene were identified. Except for chromosome 4 and 7, these candidate PvHSP70 genes were distributed on the remaining chromosomes. While the lowest number of PvHSP70 genes was determined on chromosomes 1, 3, 5, 7, 9, 10 and 11 (one HSP70 gene), the highest number of PvHSP70s was on chromosomes 6 and 8 (seven HSP70 genes each). Three genes; PvHSP70-5, -9, and -10 were found to have no-introns. In addition, four tandemly and six segmentally duplicated gene couples were detected. A total of 13 PvHSP70 genes were targeted by miRNAs of 44 plant species and the most targeted genes were PvHSP70-5 and -23. The expression profile of PvHSP70 genes based on publicly available RNA-seq data was identified and salt treated leaf tissue was found to have more gene expression levels compared to the root. qRT-PCR analysis showed that the transcript concentrations of upregulated PvHSP70 genes in leaves of Zulbiye (sensitive) were mostly higher than those of Yakutiye (resistant). The present study revealed that PvHSP70 genes might play an important role in salt stress response for common bean cultivars and variability between cultivars also suggests that these genes could be used as functional markers for salt tolerance in common bean.

Comparative Transcriptome Analysis of Resistant and Susceptible Common Bean Genotypes in Response to Soybean Cyst Nematode Infection

Soybean cyst nematode (SCN; Heterodera glycines Ichinohe) reproduces on the roots of common bean (Phaseolus vulgaris L.) and can cause reductions in plant growth and seed yield. The molecular changes in common bean roots caused by SCN infection are unknown. Identification of genetic factors associated with SCN resistance could help in development of improved bean varieties with high SCN resistance. Gene expression profiling was conducted on common bean roots infected by SCN HG type 0 using next generation RNA sequencing technology. Two pinto bean genotypes, PI533561 and GTS-900, resistant and susceptible to SCN infection, respectively, were used as RNA sources eight days post inoculation. Total reads generated ranged between ~ 3.2 and 5.7 million per library and were mapped to the common bean reference genome. Approximately 70-90% of filtered RNA-seq reads uniquely mapped to the reference genome. In the inoculated roots of resistant genotype PI533561, a total of 353 genes were differentially expressed with 154 up-regulated genes and 199 down-regulated genes when compared to the transcriptome of non- inoculated roots. On the other hand, 990 genes were differentially expressed in SCN-inoculated roots of susceptible genotype GTS-900 with 406 up-regulated and 584 down-regulated genes when compared to non-inoculated roots. Genes encoding nucleotide-binding site leucine-rich repeat resistance (NLR) proteins, WRKY transcription factors, pathogenesis-related (PR) proteins and heat shock proteins involved in diverse biological processes were differentially expressed in both resistant and susceptible genotypes. Overall, suppression of the photosystem was observed in both the responses. Furthermore, RNA-seq results were validated through quantitative real time PCR. This is the first report describing genes/transcripts involved in SCN-common bean interaction and the results will have important implications for further characterization of SCN resistance genes in common bean.

Emerging Genomic Tools for Legume Breeding: Current Status and Future Prospects

Legumes play a vital role in ensuring global nutritional food security and improving soil quality through nitrogen fixation. Accelerated higher genetic gains is required to meet the demand of ever increasing global population. In recent years, speedy developments have been witnessed in legume genomics due to advancements in next-generation sequencing (NGS) and high-throughput genotyping technologies. Reference genome sequences for many legume crops have been reported in the last 5 years. The availability of the draft genome sequences and re-sequencing of elite genotypes for several important legume crops have made it possible to identify structural variations at large scale. Availability of large-scale genomic resources and low-cost and high-throughput genotyping technologies are enhancing the efficiency and resolution of genetic mapping and marker-trait association studies. Most importantly, deployment of molecular breeding approaches has resulted in development of improved lines in some legume crops such as chickpea and groundnut. In order to support genomics-driven crop improvement at a fast pace, the deployment of breeder-friendly genomics and decision support tools seems appear to be critical in breeding programs in developing countries. This review provides an overview of emerging genomics and informatics tools/approaches that will be the key driving force for accelerating genomics-assisted breeding and ultimately ensuring nutritional and food security in developing countries.

Development of Gene-Based SSR Markers in Rice Bean (Vigna umbellata L.) Based on Transcriptome Data

Rice bean (Vigna umbellata (Thunb.) Ohwi & Ohashi) is a warm season annual legume mainly grown in East Asia. Only scarce genomic resources are currently available for this legume crop species and no simple sequence repeat (SSR) markers have been specifically developed for rice bean yet. In this study, approximately 26 million high quality cDNA sequence reads were obtained from rice bean using Illumina paired-end sequencing technology and assembled into 71,929 unigenes with an average length of 986 bp. Of these unigenes, 38,840 (33.2%) showed significant similarity to proteins in the NCBI non-redundant protein and nucleotide sequence databases. Furthermore, 30,170 (76.3%) could be classified into gene ontology categories, 25,451 (64.4%) into Swiss-Prot categories and 21,982 (55.6%) into KOG database categories (E-value < 1.0E-5). A total of 9,301 (23.5%) were mapped onto 118 pathways using the Kyoto Encyclopedia of Genes and Genome (KEGG) pathway database. A total of 3,011 genic SSRs were identified as potential molecular markers. AG/CT (30.3%), AAG/CTT (8.1%) and AGAA/TTCT (20.0%) are the three main repeat motifs. A total of 300 SSR loci were randomly selected for validation by using PCR amplification. Of these loci, 23 primer pairs were polymorphic among 32 rice bean accessions. A UPGMA dendrogram revealed three major clusters among 32 rice bean accessions. The large number of SSR-containing sequences and genic SSRs in this study will be valuable for the construction of high-resolution genetic linkage maps, association or comparative mapping and genetic analyses of various Vigna species.

Genome-wide investigation and expression analysis of AP2-ERF gene family in salt tolerant common bean

Apetala2-ethylene-responsive element binding factor (AP2-ERF) superfamily with common AP2-DNA binding domain have developmentally and physiologically important roles in plants. Since common bean genome project has been completed recently, it is possible to identify all of the AP2-ERF genes in the common bean genome. In this study, a comprehensive genome-wide in silico analysis identified 180 AP2-ERF superfamily genes in common bean (Phaseolus vulgaris). Based on the amino acid alignment and phylogenetic analyses, superfamily members were classified into four subfamilies: DREB (54), ERF (95), AP2 (27) and RAV (3), as well as one soloist. The physical and chemical characteristics of amino acids, interaction between AP2-ERF proteins, cis elements of promoter region of AP2-ERF genes and phylogenetic trees were predicted and analyzed. Additionally, expression levels of AP2-ERF genes were evaluated by in silico and qRT-PCR analyses. In silico micro-RNA target transcript analyses identified nearly all PvAP2-ERF genes as targets of by 44 different plant species' miRNAs were identified in this study. The most abundant target genes were PvAP2/ERF-20-25-62-78-113-173. miR156, miR172 and miR838 were the most important miRNAs found in targeting and BLAST analyses. Interactome analysis revealed that the transcription factor PvAP2-ERF78, an ortholog of Arabidopsis At2G28550, was potentially interacted with at least 15 proteins, indicating that it was very important in transcriptional regulation. Here we present the first study to identify and characterize the AP2-ERF transcription factors in common bean using whole-genome analysis, and the findings may serve as a references for future functional research on the transcription factors in common bean.

Wild soybean roots depend on specific transcription factors and oxidation reduction related genesin response to alkaline stress

Soil alkalinity is an important environmental problem limiting agricultural productivity. Wild soybean (Glycine soja) shows strong alkaline stress tolerance, so it is an ideal plant candidate for studying the molecular mechanisms of alkaline tolerance and identifying alkaline stress-responsive genes. However, limited information is available about G. soja responses to alkaline stress on a genomic scale. Therefore, in the present study, we used RNA sequencing to compare transcript profiles of G. soja root responses to sodium bicarbonate (NaHCO3) at six time points, and a total of 68,138,478 pairs of clean reads were obtained using the Illumina GAIIX. Expression patterns of 46,404 G. soja genes were profiled in all six samples based on RNA-seq data using Cufflinks software. Then, t12 transcription factors from MYB, WRKY, NAC, bZIP, C2H2, HB, and TIFY families and 12 oxidation reduction related genes were chosen and verified to be induced in response to alkaline stress by using quantitative real-time polymerase chain reaction (qRT-PCR). The GO functional annotation analysis showed that besides "transcriptional regulation" and "oxidation reduction," these genes were involved in a variety of processes, such as "binding" and "response to stress." This is the first comprehensive transcriptome profiling analysis of wild soybean root under alkaline stress by RNA sequencing. Our results highlight changes in the gene expression patterns and identify a set of genes induced by NaHCO3 stress. These findings provide a base for the global analyses of G. soja alkaline stress tolerance mechanisms.

De novo transcriptome sequencing of Acer palmatum and comprehensive analysis of differentially expressed genes under salt stress in two contrasting genotypes

Maple (Acer palmatum) is an important species for landscape planting worldwide. Salt stress affects the normal growth of the Maple leaf directly, leading to loss of esthetic value. However, the limited availability of Maple genomic information has hindered research on the mechanisms underlying this tolerance. In this study, we performed comprehensive analyses of the salt tolerance in two genotypes of Maple using RNA-seq. Approximately 146.4 million paired-end reads, representing 181,769 unigenes, were obtained. The N50 length of the unigenes was 738 bp, and their total length over 102.66 Mb. 14,090 simple sequence repeats and over 500,000 single nucleotide polymorphisms were identified, which represent useful resources for marker development. Importantly, 181,769 genes were detected in at least one library, and 303 differentially expressed genes (DEGs) were identified between salt-sensitive and salt-tolerant genotypes. Among these DEGs, 125 were upregulated and 178 were downregulated genes. Two MYB-related proteins and one LEA protein were detected among the first 10 most downregulated genes. Moreover, a methyltransferase-related gene was detected among the first 10 most upregulated genes. The three most significantly enriched pathways were plant hormone signal transduction, arginine and proline metabolism, and photosynthesis. The transcriptome analysis provided a rich genetic resource for gene discovery related to salt tolerance in Maple, and in closely related species. The data will serve as an important public information platform to further our understanding of the molecular mechanisms involved in salt tolerance in Maple.

Genome-wide characterization and expression analysis of common bean bHLH transcription factors in response to excess salt concentration

Members of basic helix-loop-helix (bHLH) gene family found in all eukaryotes play crucial roles in response to stress. Though, most eukaryotes carry the proteins of this family, biological functions of the most bHLH family members are not deeply evaluated in plants. In this study, we conducted a comprehensive genome-wide analysis of bHLH transcription factors in salt tolerant common bean. We identified 155 bHLH protein-encoding genes (PvbHLH) by using in silico comparative genomics tools. Based on the phylogenetic tree, PvbHLH genes were classified into 8 main groups with 21 subfamilies. Exon-intron analysis indicated that proteins belonging to same main groups exhibited a closely related gene structure. While, the PvbHLH gene family has been mainly expanded through segmental duplications, a total of 11 tandem duplication were detected. Genome-wide expression analysis of bHLH genes showed that 63 PvbHLH genes were differentially expressed in at least one tissue. Three of them displayed higher expression values in both leaf and root tissues. The in silico micro-RNA target transcript analyses revealed that totally 100 PvHLH genes targeted by 86 plant miRNAs. The most abundant transcripts, which were targeted by all 18 plant miRNA, were belonging to PvHLH-22 and PvHLH-44 genes. The expression of 16 PvbHLH genes in the root and leaf tissues of salt-stressed common bean was evaluated using qRT-PCR. Among them, two of PvbHLHs, PvbHLH-54, PvbHLH-148, were found to be up-regulated in both tissues in correlation with RNA-seq measurements. The results of this study could help improve understanding of biological functions of common bean bHLH family under salt stress. Additionally, it may provide basic resources for analyzing bHLH protein function for improving economic, agronomic and ecological benefit in common bean and other species.

The Medicago sativa gene index 1.2: a web-accessible gene expression atlas for investigating expression differences between Medicago sativa subspecies

Background

Alfalfa (Medicago sativa L.) is the primary forage legume crop species in the United States and plays essential economic and ecological roles in agricultural systems across the country. Modern alfalfa is the result of hybridization between tetraploid M. sativa ssp. sativa and M. sativa ssp. falcata. Due to its large and complex genome, there are few genomic resources available for alfalfa improvement.

Results

A de novo transcriptome assembly from two alfalfa subspecies, M. sativa ssp. sativa (B47) and M. sativa ssp. falcata (F56) was developed using Illumina RNA-seq technology. Transcripts from roots, nitrogen-fixing root nodules, leaves, flowers, elongating stem internodes, and post-elongation stem internodes were assembled into the Medicago sativa Gene Index 1.2 (MSGI 1.2) representing 112,626 unique transcript sequences. Nodule-specific and transcripts involved in cell wall biosynthesis were identified. Statistical analyses identified 20,447 transcripts differentially expressed between the two subspecies. Pair-wise comparisons of each tissue combination identified 58,932 sequences differentially expressed in B47 and 69,143 sequences differentially expressed in F56. Comparing transcript abundance in floral tissues of B47 and F56 identified expression differences in sequences involved in anthocyanin and carotenoid synthesis, which determine flower pigmentation. Single nucleotide polymorphisms (SNPs) unique to each M. sativa subspecies (110,241) were identified.

Conclusions

The Medicago sativa Gene Index 1.2 increases the expressed sequence data available for alfalfa by ninefold and can be expanded as additional experiments are performed. The MSGI 1.2 transcriptome sequences, annotations, expression profiles, and SNPs were assembled into the Alfalfa Gene Index and Expression Database (AGED) at http://plantgrn.noble.org/AGED/, a publicly available genomic resource for alfalfa improvement and legume research.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1718-7) contains supplementary material, which is available to authorized users.

Development and Validation of EST-SSR Markers from the Transcriptome of Adzuki Bean (Vigna angularis)

The adzuki bean (Vigna angularis (Ohwi) Ohwi and Ohashi) is an important grain legume of Asia. It is cultivated mainly in China, Japan and Korea. Despite its importance, few genomic resources are available for molecular genetic research of adzuki bean. In this study, we developed EST-SSR markers for the adzuki bean through next-generation sequencing. More than 112 million high-quality cDNA sequence reads were obtained from adzuki bean using Illumina paired-end sequencing technology, and the sequences were de novo assembled into 65,950 unigenes. The average length of the unigenes was 1,213 bp. Among the unigenes, 14,547 sequences contained a unique simple sequence repeat (SSR) and 3,350 sequences contained more than one SSR. A total of 7,947 EST-SSRs were identified as potential molecular markers, with mono-nucleotide A/T repeats (99.0%) as the most abundant motif class, followed by AG/CT (68.4%), AAG/CTT (30.0%), AAAG/CTTT (26.2%), AAAAG/CTTTT (16.1%), and AACGGG/CCCGTT (6.0%). A total of 500 SSR markers were randomly selected for validation, of which 296 markers produced reproducible amplicons with 38 polymorphic markers among the 32 adzuki bean genotypes selected from diverse geographical locations across China. The large number of SSR-containing sequences and EST-SSR markers will be valuable for genetic analysis of the adzuki bean and related Vigna species.