Transcriptomic analysis of salt tolerance-associated genes and diversity analysis using indel markers in yardlong bean (Vigna unguiculata ssp. sesquipedialis)
BMC Genomic Data volume 22, Article number: 34 (2021)
High salinity is a devastating abiotic stresses for crops. To understand the molecular basis of salinity stress in yardlong bean (Vigna unguiculata ssp. sesquipedalis), and to develop robust markers for improving this trait in germplasm, whole transcriptome RNA sequencing (RNA-seq) was conducted to compare the salt-tolerant variety Suzi 41 and salt-sensitive variety Sujiang 1419 under normal and salt stress conditions.
Compared with controls, 417 differentially expressed genes (DEGs) were identified under exposure to high salinity, including 42 up- and 11 down-regulated DEGs in salt-tolerant Suzi 41 and 186 up- and 197 down-regulated genes in salt-sensitive Sujiang 1419, validated by qRT-PCR. DEGs were enriched in “Glycolysis/Gluconeogenesis” (ko00010), “Cutin, suberine and wax biosynthesis” (ko00073), and “phenylpropanoid biosynthesis” (ko00940) in Sujiang 1419, although “cysteine/methionine metabolism” (ko00270) was the only pathway significantly enriched in salt-tolerant Suzi 41. Notably, AP2/ERF, LR48, WRKY, and bHLH family transcription factors (TFs) were up-regulated under high salt conditions. Genetic diversity analysis of 84 yardlong bean accessions using 26 InDel markers developed here could distinguish salt-tolerant and salt-sensitive varieties.
These findings show a limited set of DEGs, primarily TFs, respond to salinity stress in V. unguiculata, and that these InDels associated with salt-inducible loci are reliable for diversity analysis.
The legume cowpea (Vigna unguiculata L. Walp.) is the fifth most widely consumed plant-based source of protein and soluble fiber , and the sesquipedalis subspecies, i.e., asparagus bean or ‘yardlong’ bean, is cultivated as a prized vegetable among eastern and southern Asian countries [2, 3]. Abiotic stress induced by high salinity can lead to major reductions in growth, yield, and quality, so improvement to salt tolerance represents an urgent priority for yardlong bean breeding programs. Uncovering the molecular mechanisms underlying plant response to salt stress can enable development of salt-tolerant yardlong bean cultivars. To date, several mechanisms have been identified across a range of model plants and crops for their role in tolerance to high salinity, including modulation of ion and osmotic homeostasis, stress-induced cellular repair pathways, and alternative growth regulatory pathways that circumvent stress response signaling .
Salt-tolerant plants characteristically exhibit adaptive maintenance of intracellular ion homeostasis, and in particular, the salt overly sensitive (SOS) pathway has been implicated in maintaining a K+/ Na+ ratio essential for growth under high salinity conditions. The SOS pathway involves regulation of ion transport by the SOS1 Na+/H+ plasma membrane antiporter, which is activated via the SOS3 calcium sensor and SOS2 Ser/Thr protein kinase [5, 6]. Other known regulators of ion transport and exclusion include Arabidopsis K+ transporter1 (AKT1), Na+/H+ exchangers (NHXs), high sodium affinity transporter (HKT), and other plasma membrane proteins (PMP), all of which may be activated under exposure to high salinity to ensure an ion balance that allows continued cellular function [4, 7]. In addition to transporters, transcription factors from several families participate in ion homeostasis and salt tolerance through regulation of signal transduction pathways and downstream transporters, such as apetala2/ethylene responsive factor (AP2/ERF), dehydration responsive element binding protein (DREB), basic leucine zipper domain (bZIP), WRKY, and MYB, among others [8,9,10,11,12].
Osmotic homeostasis is also reportedly regulated also by different MAP kinase (MAPK) signal cascade-mediated programmed responses that control osmotic homeostasis, for example through vacuolar Na+ sequestration or via synthesis and accumulation of biocompatible osmolytes [13,14,15]. In addition, salt stress is typically accompanied by reactive oxygen species (ROS) burst that can disrupt metabolic activity or damage lipid membranes and DNA . Plants have thus evolved multiple enzymes to scavenge and detoxify ROSs, minimize their damage, and enhance repair of cellular damage including superoxide dismutase, ascorbate peroxidase, catalase, guaiacolperoxidase, and others. Furthermore, plants synthesize metabolites and small molecules that also function as antioxidants, such as ascorbic acid, alkaloids, carotenoids, flavonoids, phenolic compounds, and tocopherol, etc. [17, 18].
Since its introduction from Africa, yardlong bean has been increasingly selected for stress-resistant phenotypes suitable for cultivation in Asia. Chen et al. (2007) and Murilloamador et al. (2006) both identified salt tolerant sequipedalis genotypes [19, 20], while more recently, Xu and colleagues used genome-wide association study (GWAS) to reveal thirty-nine SNP loci associated with drought resistance . Tan et al. (2016) identified several genes that were differentially expressed genes (DEGs) between cold-tolerant and -sensitive yardlong bean cultivars, while Pan et al. (2019) found 216 and 127 salt stress-associated DEGs in roots and leaves, respectively, six of which were linked to 17 salt tolerance-associated SNP markers [22, 23]. More recently, other QTLs associated with salt tolerance in yardlong bean were mapped using a population generated by crossing Suzi 41 (salt tolerant) × Sujiang 1419 (salt sensitive) . Completion of the yardlong bean genome and the relatively low cost of re-sequencing have enabled further development of SNP/InDel markers in sesquipedalis for use in breeding and genetic analysis, as has been widely reported in common bean and mungbean [25,26,27].
In this study, RNA-seq analysis was used to compare transcriptional responses of two varieties of yardlong bean, Suzi 41 (salt-tolerant) and Sujiang 1419 (salt-sensitive), to identify the regulatory and metabolic pathways mediating salt stress response in this high value crop. A set of DEGs encoding transcription factors were identified which could regulate downstream pathways necessary for salt tolerance. In addition, KEGG and GO analysis were performed to predict the putative functions of DEGs, and then the differences in transcriptional regulation between the salt-tolerant and -sensitive varieties were compared. Importantly, this study developed a set of informative and reliable salt stress-specific InDel markers based on high throughput sequencing and revealed considerable genetic diversity in V. unguiculata. This work provides insight into the basic mechanisms underlying salt tolerance, as well as tools for applied research necessary for improvement of yardlong bean varieties for cultivation in high saline soils.
Transcriptome sequencing and discovery of novel transcripts
The Illumina HiSeq™ 2000 platform was used to sequence Suzi 41 and Sujiang 1419 transcriptomes in yardlong bean that were treated under high salt stress conditions (41S and 1419S) to compare with those of unstressed plants (41C and 1419C) in order to identify differences in their transcriptional responses to high salinity by these two phenotypically different varieties. After removing low-quality sequences and trimming adapter sequences, ~ 24 million paired-end reads were generated from each of the cDNA libraries with an average GC content of 45.65%. All clean reads were matched to the Vigna unguiculata reference genome by TopHat software. As a result, about 43 million mapped reads were obtained for each line of Suzi 41 and Sujiang 1419, with an average matching rate of 89.83% (Supplementary Table S1). Most (99.56–99.72%) of the reads with matches were unique reads (matching only one yardlong bean locus), while the remainder (~ 0.28–0.44%) were non-unique (matching more than one yardlong bean locus) or unaligned. For more detailed investigation of gene expression in the different treatments, only unique reads were used in the analysis. In both control and salt stress treatments, the numbers of mapped genes in Suzi 41 (19,606 and 19,737 genes) were found to be similar to those in Sujiang 1419 (19,433 and 19,594 genes, respectively). The mapped genes among the four treatments (41C, 41S, 1419C, and 1419S) were further compared, and ~ 95% of them were present in at least two treatments (Fig. 1).
Identification of novel transcript isoforms has emerged as one of the major advantages of RNA-seq analysis. This study revealed a total of 563 novel transcript isoforms in Suzi 41 and Sujiang 1419 yardlong bean varieties. Comparison of transcriptomic reads with the Vigna unguiculata reference genome revealed that most new genes (562; 99.82%) were annotated by nr, followed by GO (361; 64.12%) and Swissprot (319; 56.66%). Only 64 (11.37%) DEGs were annotated with COG (Supplementary Table S2). Respectively, 299 (53.11%), 243 (43.16%) and 90 (15.99%) DEGs were annotated with Pfam, KOG and KEGG. Although the novel transcript isoforms will be validated in future experiments, they were included in further analyses for preliminary functional characterization and investigation of their putative role in abiotic stress responses.
Differential gene expression in response to salt-stress treatments
Differential gene expression analysis of Suzi 41 and Sujiang 1419 genotypes revealed 390 differentially expressed genes (DEGs) in the salt stress vs control comparison (Fig. 2, Supplementary Table S3). There were 42 and 11 genes identified as being up- and down-regulated in the salt-tolerant genotype Suzi 41, respectively, and173 and 183 genes identified in the salt-sensitive genotype Sujiang 1419, respectively. There were more DEGs in Sujiang 1419 than in Suzi 41.
Under salt tolerance, a number of genes were expressed only in the salt-tolerant genotype. In Suzi 41, there were 32 and 2 DEGs were identified as being up- and down-regulated that were not also differentially expressed in Sujiang 1419 (Supplementary Table S4). Interestingly, the most highly up-regulated of these are LR48 transcription factors including (with Log2 fold change) Vigun02g152900 (2.00), Vigun10g012000 (1.91), Vigun10g011500 (1.82), and Vigun10g011900 (1.60), as well as a PR-4-like pathogenesis-related protein Vigun06g113200 (1.64). The two differentially down-regulated genes include a LR48 protein, Vigun05g219700 (− 1.17), and a WAT1-related protein Vigun06g228300 (− 1.02). The prevalence of transcription factors among DEGs strongly suggests that these genes are responsible for promoting salt tolerance in Suzi 41, and may serve as potentially strong candidates for further elucidation of the mechanisms underlying salt tolerance in yardlong bean.
Functional annotation and classification of DEGs
Ten DEGs were up-regulated and 9 DEGs were down-regulated in both varieties (Supplementary Table S5), suggesting that these genes are differentially expressed specifically under salt stress in both varieties. Among the most highly up-regulated overlapping DEGs are a hypothetical bZIP LR48 gene (annotated in Phytozome as senescence-associated protein PF06911) (Vigun11g188200), as well as several other transcription factors, and a predicted peroxidase 21 (Vigun07g080600). Among the down-regulated genes found in both varieties, metalloendoproteinase 1-like protein (Vigun02g070900) and an alcohol dehydrogenase (Vigun09g123700), both with a Log2 fold change of ~ − 1.0 in susceptible and tolerant yardlong bean varieties, and several hypothetical LR48 transcription factors were identified. These genes, which were differentially expressed under salt stress in both varieties may serve as a basis for identifying target genes for molecular breeding to improve salt tolerance in yardlong bean.
Next, GO analysis was conducted to predict the potential functions or biological roles of these salt-induced DEGs. GO terms that were enriched among the 52 most significantly up- or down-regulated DEGs (14 in Suzi 41; 40 in Sujiang 1419) under salt stress indicated that these genes were likely related to biological processes and molecular functions. More specifically, the DEGs in Suzi 41 were enriched in biological processes such as “hydrogen peroxide catabolic process” (GO:0042744) and “response to oxidative stress” (GO:0006979), while molecular function-associated terms included “peroxidase activity” (GO:0004601). By contrast, the DEGs in Sujiang 1419 were enriched in “suberin biosynthetic process” (GO: 0010345), “transition metal ion transport” (0000041), and “peptidyl-proline hydroxylation” (GO:0019511) biological processes, while molecular function-associated terms included “alcohol dehydrogenase (NAD) activity” (GO:0004022) and “potassium ion binding” (GO:0030955). It is noteworthy that “alcohol dehydrogenase (NAD) activity” (GO:0004022) was the only term enriched in both varieties (Fig. 3, Supplementary Table S6).
To identify the metabolic pathways in which the DEGs were involved and enriched, KEGG analysis was also performed . The pathways enriched for the most highly up- or down-regulated significant DEGs are listed in Fig. 4. Among these pathways, “Glycolysis/Gluconeogenesis” (ko00010, p-value = 2.32E-05), “Cutin, suberine and wax biosynthesis” (ko00073, p-value = 0.0004), and “phenylpropanoid biosynthesis” (ko00940, p-value = 0.0126) etc., were enriched in Sujiang 1419. The only significantly enriched biological pathway found in Suzi 41 was “cysteine and methionine metabolism” (ko00270, p-value = 0.0032), in which an ethylene biosynthetic enzyme-encoding gene (Vigun02g178400), 1-aminocyclopropane-1-carboxylate synthase (ACS, EC:18.104.22.168), was significantly up-regulated (Supplementary Table S7). This finding thus suggested that ethylene signaling may contribute a major role in tolerance to salt stress for V. unguiculata subsp. sesquipedalis.
Differential expression of transcription factors between the two varieties under salt stress
Transcription factors play crucial roles in regulating the expression of stress response genes during exposure to high salinity. A total of 224 differentially expressed TFs (47 families) were identified under salt stress in Suzi 41 and Sujiang 1419 (Supplementary Table S8). These TFs include MYB, B3, NAC, AP2/ERF, MADS, GNAT, plant basic helix–loop–helix (bHLH), C2H2, and WRKY. MYB composed the largest percentage (19 TFs, 16.38%), followed by B3 (15 TFs, 12.93%), NAC (13 TFs, 11.21%), and AP2/ERF (13 TFs, 11.21%), indicating that these TFs may be major determinants controlling the mechanisms of salt stress tolerance in yardlong bean. Several of the transcription factors, such as NAC and MYB, which are known to be induced by exposure to high salt conditions in Arabidopsis thaliana, halophytic Suaeda liaotungensis, wheat, and rice [29,30,31,32], were highly expressed under salinity stress in both Suzi 41 and Sujiang 1419.
Specifically, 17 transcription factors were found to be significantly up- or down-regulated only in salt tolerant Suzi 41 at a Log2 fold change of 1.0 or higher (Table 1). Three out of six of the most up-regulated TFs were hypothetical LR48 proteins including Vigun06g162600, Vigun08g102200, and Vigun02g140100, which were up-regulated 1.23, 1.27, and 1.47, respectively. The two most up-regulated transcription factors found only in Suzi 41, Vigun06g141200 and Vigun11g052100 (Log2 fold change = 1.47 and 1.36, respectively), are both MADS-M-Type TFs, the latter of which is annotated as a probable TAT2 aminotransferase. Two AP2/ERF genes were a WAT1-related AP2/ERF family protein (Vigun06g228300) and Ethylene-responsive transcription factor (Vigun07g178200), which had Log2 fold lower expression = − 1.02 and 1.04 in salt-stressed plants compared to non-stressed plants. These genes may play an important role in plant response to salt stress. In addition, three other families, ABC transporter (Vigun02g026100), WRKY (Vigun01g071800), and bHLH (Vigun03g388100), showed 1.14, 1.07, and 1.00, Log2 fold up-regulation, respectively, implying that these TFs may also participate in salt-regulated pathways. The differential expression of these genes exclusively in Suzi 41 during salt exposure suggests these genes may be apt targets for molecular breeding for increased salt tolerance in yardlong bean.
Twelve transcription factors were significantly up- or down-regulated in both Suzi 41 and Sujiang 1419 (Supplementary Table S9). Notably, Vigun08g132400, a hypothetical MADS-M-Type TF showed Log2 fold increases of 2.67 and 1.26 in salt tolerant and salt sensitive varieties, respectively. Similarly, 305-like MYB-related protein Vigun07g057300 exhibited Log2 fold up-regulation of 1.50 and 2.11 in the tolerant and sensitive lines, respectively. Among differentially down-regulated TFs, MADS-M-Type Vigun06g084400 was Log2 fold lower in both Suzi 41 and Sujiang 1419 at − 1.12 and − 1.93, respectively. While these genes are not indicative of a salt stress response unique to Suzi 41, their differential expression in both varieties during salt exposure suggests that they belong to potentially universal stress response pathways across yardlong bean varieties.
qRT-PCR validation of differentially expressed genes
To validate the results of DEG identification in the RNAseq data, qRT-PCR was used to confirm the differential up- or down-regulation for 12 randomly selected DEGs (Fig. 5A). Six of the 12 DEGs (Vigun03g411300, Vigun04g067600, Vigun07g057300, Vigun07g182700, Vigun08g132400, Vigun11g212700) were up-regulated in Suzi 41 and Sujiang 1419 during NaCl treatment, while the other six DEGs (Vigun02g070900, Vigun06g084400, Vigun07g122000, Vigun08g068600, Vigun08g068700, Vigun09g123700) were down-regulated in the two genotypes. Although the gene expression values are different using both techniques, all of these genes displayed similar expression trend. Comparison of the qRT-PCR and RNAseq data revealed a high correlation for the selected unigenes (Pearson r = 0.7835, Fig. 5B), ultimately reflecting consistency between the qRT-PCR and transcriptomic results.
Several DEGs associated with salt tolerance are located on chromosome 11 of the yardlong bean genome
Based on our previous studies, six major QTLs associated with salt tolerance were detected in the region between VUIn584 ~ VUIn724, between VUIn282 ~ VUIn815, and between VUIn675 ~ VUIn578 on chromosomes 9, 11, and 8 . In the present study, nine DEGs in the region of VUIn282 ~ VUIn815 on Chr.11 were found, eight of which were differentially expressed only in Suzi 41 under salt treatment (Table 2). Three of these DEGs encoded hypothetical LR48 proteins, Vigun11g159900, Vigun11g188200 (the only one differentially expressed in both Suzi 41 and Sujiang 1419), and Vigun11g170200, which were up- or down-regulated 1.01, 1.10, and − 1.20, respectively. Furthermore, the two most up-regulated DEGs found only in Suzi 41, Vigun11g186800 and Vigun11g174500 (Log2 fold change = 1.69 and 1.54, respectively), are respectively annotated as a probable BOI-related E3 ubiquitin-protein ligase and an organ-specific protein. In addition, four other DEGs, aldehyde dehydrogenase family (Vigun11g160800), hypothetical protein (Vigun11g170300), myosin heavy chain kinase (Vigun11g170300), and putative lipid-transfer protein DIR1 (Vigun11g172800), showed Log2 fold lower expression = − 1.31, − 1.28, 1.08, and − 1.22 in salt-stressed plants compared to non-stressed plants (Table 2). The positions of these DEGs have not been previously reported, and their mapping to on Chr.11 thus provides a valuable resource for studying salt stress-induced genes in yardlong bean.
26 polymorphic InDels in DEGs between Suzi 41 and Sujiang 1419 show allelic differences and high genetic diversity among yardlong bean germplasm
Genome-wide identification of insertion/deletion polymorphisms was conducted via TopHat 2.0 software. In total, 175 InDels located in DEGs were identified in the RNAseq data (Fig. 6, Supplementary Table S10) in Suzi41 and Sujiang 1419, including 134 InDels in Suzi41 which were distributed across all the eleven chromosomes, varying from 17 on Chr.03 to six each on Chr.04 and Chr.05. At the same time, 147 InDels were identified in Sujiang 1419 that ranged from 20 on Chr.07 to 7 on Chr.05 (Supplementary Table S11). Among these, the largest InDel was 42 bp, although InDels smaller than 3 bp were prevalent and accounted for about 80% of the total. The proportions of InDels less than 10 bp were 94.78 and 94.56% in Suzi 41 and Sujiang 1419, respectively, and InDels smaller than 6 bp accounted for 87.31 and 89.80%, respectively (Supplementary Table S12). The genome sequences of two different varieties were used to validate the 175 InDel polymorphisms identified through RNA-seq. Among the tested markers, 26 primer pairs (14.86%) produced clear amplicons with the expected sizes in both yardlong bean varieties, while 127 primer pairs (72.57%) produced amplicons in only one genotype and therefore were not suitable for genetic analysis, 7 (4.00%) were monomorphic, and 15 (8.57%) failed to amplify altogether (Fig. 7). To further confirm whether these 26 InDel polymorphisms could be used as markers in other yardlong bean accessions, these primer pairs were validated in 84 additional germplasm accessions of yardlong bean (Supplementary Table S13). A total of 58 alleles were detected and scored (Supplementary Table S14). The genetic diversity index ranged from 0.0887 to 0.5907, and the average was 0.4366. The polymorphic information content (PIC) values ranged from 0.0848 to 0.5103, and the average was 0.3494. The markers used in this study represent a broad level of diversity in V. unguiculata and can effectively discriminate between many publicly available germplasm accessions. Overall, these InDel markers developed here reveal a high number of alleles for these DEGs and substantial genetic diversity among accessions.
In this study, the whole transcriptomes of salt tolerant variety Suzi 41 with that of the salt sensitive variety Sujiang 1419 were compared under high salinity conditions to determine which genes confer tolerant phenotype in V. unguiculata subsp. sesquipedalis. Using RNA-seq, transcriptional profiles were developed for salt stress responsive gene networks in salt-tolerant and -susceptible varieties. Under high salinity, stress tolerant plants apparently had fewer DEGs than sensitive plants compared with control plants. These results aligned with other transcriptomic studies of response to drought and salt stress in sorghum and Corchorus spp. [33, 34]. Consistent with this finding, fewer GO categories were enriched with DEGs in Suzi 41 than in Sujiang 1419 under salt stress. In Suzi 41, five and four GO terms were enriched with DEGs in the biological process and molecular function categories, respectively. Among the enriched biological processes in Suzi 41 were “hydrogen peroxide catabolic process” (GO: 0042744) and “response to oxidative stress” (GO: 0006979). Hydrogen peroxide plays an important role in stress response by coordinating intracellular and systemic signaling systems that facilitate plant acclimatization and tolerance to stress . Oxidative stress results from plant inability to mitigate free radical damage with sufficient antioxidant activity under salt stress . These two GO terms are only associated with the salt-tolerant yardlong bean variety and were primarily enriched with upregulated DEGs. These results suggest that salt-tolerant varieties can acclimate to salinity stress by relieving oxidative stress through the stimulation of salt resistance-related gene expression via ROS signals and regulation of redox pathways. However, there are 10, 1, and 17 GO categories enriched in biological process, cellular component and molecular function respectively in Sujiang 1419. “Alcohol dehydrogenase activity” (GO: 0004022) is the only GO category enriched (entirely with downregulated DEGs) in both Suzi 41 and Sujiang 1419. Alcohol dehydrogenase (ADH) activity is considered a necessary condition for plant survival under anaerobic conditions. The role of plant ADH gene in anaerobic stress response has been a long-standing focus of research [37, 38]. Our study shows that alcohol dehydrogenase may participate in the salinity stress response by yardlong bean.
By mapping these DEGs, three major QTL regions and a set of 26 InDel markers were successfully identified that could reliably genotype a panel of 84 V. unguiculata germplasm accessions. In addition, LR48 and ERF family transcription factors were found to potentially play a major role in tolerance to salt stress given their prevalence among the highest up-regulated DEGs in Suzi 41. Finally, genetic diversity was assessed through PIC analysis and validated InDel markers for salt response using a wide panel of yardlong bean germplasm accessions.
Although few significant DEGs were correlated with high salt treatment in the Suzi 41 variety, comparisons between high salinity treatment and control plants revealed strong induction of Vigun07g178200 (ethylene-responsive transcription factor 1-ERF1), suggesting its participation in a tolerant response to salinity. ERF1 up-regulation has been widely described as a step in the ethylene-mediated salt stress response signal cascade [39,40,41]. Importantly, AP2/ERFs are widely known to regulate hormone signaling and abiotic stress response pathways, as reviewed by Xie et al. 2019 . For example, transgenic expression of GmERF3 from soybean in Nicotiana tabacum enhanced its tolerance to salinity, as well as other abiotic stresses . In addition, Vigun02g026100 (ABC transporter G family member 21) was differentially up-regulated in Suzi 41 (Table 1). In plants, ABC transporters are typically associated with heavy metal detoxification  or with auxin transport [45, 46]. Another recent investigation of salt stress response in V. unguiculata reported similar numbers of DEGs for tolerant and sensitive varieties (i.e., 13 DEGs from six different TF families) and identified 17 SNP markers associated with six salt-induced DEGs . Similarly, our study also found six major QTLs across chromosomes 8, 9, and 11, one of which contained three differentially expressed LR48 family TFs.
Interestingly, the association of these three salt-inducible LR48 transcription factors (Vigun11g159900, Vigun11g170200, and Vigun11g188200) suggested that these genes could contribute a potentially important role in tolerance to salt stress in yardlong bean. Although surprisingly little has been reported on the structure, domains, or mechanistic function of the LR48 gene, it is commonly used in marker-assisted breeding to confer hypersensitive response-mediated resistance to Leaf Rust (Puccinia triticina) in wheat [47,48,49]. In addition to potential disease resistance, which requires further study, the preponderance of LR48-like genes in our dataset strongly suggests that they could function in abiotic stress response in V. unguiculata. In addition to the LR48 and ERF transcription factors, a bHLH35 transcription factor (Vigun03g388100), a WRKY75 transcription factor (Vigun01g071800), and an AP2-ERF ethylene-responsive transcription factor (Vigun07g178200) were identified among the highly significant, Suzi 41-specific salt-inducible DEGs.
In agreement with our findings that AP2-ERF, WRKY and bHLH transcription factors for their critical role in osmotic stress signaling mediated by salt [8, 50,51,52]. Previous studies in Arabidopsis have revealed the detailed regulatory role WRKY8 in modulating tolerance to salt. Specifically, WRKY8 was found to bind downstream stress response genes under salt exposure, and its knockout resulted in an increased Na+/K+ ratio, hypersensitivity to salt, and other developmental abnormalities . In addition, bHLH transcription factors, such as AtbHLH112 in Arabidopsis, mediate tolerance to high salinity. In this example, AtbHLH was shown to localize to the nucleus and bind GCG- and E-box motifs in target gene promoters during salt and drought treatment. Its functionality was correlated with enhanced salt tolerance, higher proline levels, and elevated expression of POD and SOD genes to mitigate ROS damage . Notably, bHLH genes have been shown to affect plant physiological response to stress, such as SlbHLH22, which was found to be elevated under high salinity- or D-mannitol-induced stress in tomato . The prevalence of LR48, bHLH, WRKY, and AP2-ERF TFs among the stress-induced DEGs suggests that these genes may serve as promising targets for improvement of salt tolerance in Vigna unguiculata.
In addition to marker development and transcriptomic profiling, genetic diversity analysis was conducted as a necessary step in accessing the full potential of these genomic and transcriptomic resources. Molecular markers are effective tools to evaluate the genetic diversity, and disclose the evolution history of cowpea resources. For example, Asare et al. (2010) analyzed the genetic diversity of 141 cowpea in Ghana using 20 pairs of SSR primers; Xiong et al. (2016) investigated the genetic polymorphism of 784 cowpea genotypes worldwide using SNP markers, and predicted the migration and domestication history of cowpea [56, 57]. In our study, the InDel markers based on high-throughput transcriptome sequencing are different from the molecular makers based on DNA polymorphism, since they locate in coding sequences that directly related to the phenotype characters. A total of 175 InDel markers were developed in our research, 26 of which were used to evaluate the genetic diversity of 84 yardlong bean accessions. A high level of diversity in V. unguiculata was found by using the InDel markers developed here. Collectively, these results indicated that the InDels not only can serve as effective markers, comparable to SSRs, for estimating genetic diversity in yardlong bean, but also provide a basis for future research on gene function.
This transcriptomic analysis of salt-inducible genes in yardlong bean revealed a suite of candidate DEGs largely comprised of ERF, LR48, WRKY, and bHLH transcription factors. A robust set of 26 salt stress-related InDel markers were developed that can be used for improvement of salt tolerance in yardlong bean germplasm. Finally, genetic diversity analysis in a wide panel of accessions was performed which demonstrated the effectiveness and reliability of these InDel markers. This work therefore provides a genetic resource for improving yields in low quality soils, and offers foundational insights into the basic mechanisms underlying abiotic stress response in V. unguiculata.
The seeds for this study contained salt-tolerant Suzi 41 and salt-sensitive Sujiang 1419 , and 84 other yardlong bean accessions, among which 77 were collected from the National Infrastructure for Vegetable Crop Germplasm Resources (NIVCGR), and 7 from Jiangsu Academy of Agricultural Sciences (JAAS) (Supplementary Table S13).
Plant growth and salt stress treatments
Two yardlong bean genotypes: salt-tolerant Suzi 41 and salt-sensitive Sujiang 1419  were employed to examine differences in the expression of genes involved in a tolerant response to high salinity. Ten seeds of Suzi 41 and Sujiang 1419 were sown in cups (9.5 × 16 cm) filled with vermiculite, with three replicates per genotype. Following germination, four plants of each variety were placed in a plastic tank (50 × 40 × 20 cm) filled with aerated half-strength Hoagland nutrient solution  and allowed to grow until the true leaves fully expanded. Four seedlings were then placed in half-strength Hoagland nutrient solution containing 50 mM NaCl for 1 day. The concentration was raised by adding 50 mM NaCl each day until a final concentration of 150 mM was reached. At the same time, untreated seedlings were transferred to a tank filled with half-strength Hoagland nutrient solution without added salt to serve as the control. The roots were harvested separately after 3 days of treatment. Each sample was derived from at least four individual plants with three biological replicates per genotype for each treatment. The plant roots were frozen in liquid nitrogen and kept at − 80 °C.
Total RNA was isolated using the RNAprep Pure Plant Kit (Tiangen Biotech Co., Ltd., Beijing, China) according to the manufacturer’s protocols. RNA samples were then visualized on a 1% agarose gel and quantified with a NanoDrop ND-1000 spectrophotometer (Thermo Fisher Scientific, Inc. Waltham, MA, USA). Twelve root samples (Suzi 41 salt-stressed, 41S; Sujiang 1419 salt-stressed, 1419S; Suzi 41 control, 41C; Sujiang 1419 control, 1419C) were used for transcriptome sequencing by the Biomarker Biotechnology Corporation (Beijing, China). For each sample, three independent biological replicates were performed.
RNA-seq libraries were prepared using the paired-end strategy. In detail, (1) Poly (A) mRNA was enriched using the NEBNext® Poly (A) mRNA Magnetic Isolation Module (New England Biolabs, Ipswich, MA, USA), and then it was fragmented into short pieces chemically. (2) The first- and second-strand cDNA were synthetized using the short mRNA as template and then subjected to end-repair and phosphorylation using T4-DNA polymerase and Klenow DNA polymerase. The repaired cDNA fragments were inserted ‘A’ bases as overhangs at the 3′ ends and connected with sequencing adapters. (3) The suitable fragments were selected for the PCR amplification as templates after agarose gel electrophoresis. Finally, the twelve libraries were sequenced using Illumina HiSeq™ 2000 sequencing system.
After RNA-seq, the raw data were purified by trimming adapters and removing low-quality sequencing to get clean reads. At the same time, Q20, Q30 and GC content of the clean data were calculated. Reference genome and gene annotation files were downloaded from the website (https://phytozome.jgi.doe.gov). All clean reads were matched to the Vigna unguiculata reference genome by TopHat v2.0 .
Analysis of differentially expressed genes
The expression levels of genes were calculated using the FPKMs (fragments per kilo-base of exon per million fragments) method, which were computed by summing the FPKMs of transcripts in each gene group . Here, only fold change with an absolute value of P ≤ 0.01 and |log2(ratio)| ≥ 1 were used as the threshold for significant differential expression and for subsequent analysis.
Differential gene functional annotation
The KEGG pathways were analyzed for differentially expressed genes and the corresponding ko numbers were predicted using KOBAS software . A statistical analysis of the GO (Gene Ontology) term for genes in the biological process, cellular component, and molecular function classifications was implemented by the GOseq R package (1.10.1) , in which gene length bias was corrected. GO terms with p < 0.05 were considered significantly enriched by differentially expressed genes.
Validation by real-time PCR (qRT-PCR)
In order to validate the reliability of RNA-seq experiments, a total of 12 DEGs were randomly selected for qRT-PCR analysis of relative expression. Sequences of the specific primers used for qRT-PCR are given in Supplementary Table S15. A total of 0.5 μg of DNaseI-treated total RNA was converted into single-stranded cDNA using a Prime-Script 1st Strand cDNA Synthesis Kit (TaKaRa, Dalian, China). The cDNA templates were then diluted 20-fold before use. The quantitative reaction was performed on a CFX96 Real-Time PCR Detection System (Bio-Rad, Singapore) using SYBR Premix Ex Taq™ (TaKaRa, Dalian, China). PCR amplification was performed under the following conditions: 30 s at 95 °C, followed by 40 cycles of 95 °C for 5 s, 60 °C for 20 s and then 72 °C for 20 s. The relative expression levels of the selected DEGs normalized to an internal reference gene were calculated using the 2-ΔΔCt method . The UBC9 housekeeping gene (Vigun05g084500) was used as the internal reference, and all analyses were performed with three technical and three biological replicates.
Identification of short InDels and primer design
The transcriptomic sequencing results of salt-tolerant Suzi 41 and salt-sensitive Sujiang 1419 were mapped to the cowpea reference genome (https://phytozome.jgi.doe.gov) using TopHat 2.0 . With reads from each genotype and treatment mapped to the reference genome, GATK v3.5  was used to call SNPs and InDels for each sample. After filtering out unreliable sites, the final set of SNPs and InDels in VCF format were obtained. To develop single-product, amplifiable InDel markers, primers were designed to amplify each InDel with 150 bp 5′ and 3′ flanking sequence using Primer 3.0  with a total amplicon lengths between 100 and 250 bp. The primers were constrained by a required melting temperature between 58 and 65 °C, primer optimal length of 20 nt (18–25 nt), minimum GC of 40%, maximum GC of 65%, and optimized GC at 50%. To avoid selection of microsatellites around InDels, only those InDels in which none of the flanks contained microsatellites identified by MISA  with default parameters, were used for further primer design.
DNA extraction and PCR
Genomic DNA of each yardlong bean accession was extracted from young leaves using the Hexadecyltrimethylammonium bromide (CTAB) method . A Nano-Drop 2000 1spectrophotometer (Nano Drop Technologies, USA) was used to evaluate the quality and concentration of all DNA. DNA samples were diluted to 25 ng/μL. PCR was performed in a total volume of 10 μl containing 50 ng of genomic DNA, 0.4 U of Taq DNA polymerase (Dingguo Biological Technology Development Co., Ltd., Beijing, China), 10× Taq Buffer II, 25 mM MgCl2, 2 mM of dNTPs, and 1 mM each forward and reverse primer. PCR conditions were as follows: 94 °C for 4 min, 30 cycles of 30 s at 94 °C, 52 ~ 58 °C for 30 s, 72 °C for 30 s, and 1 cycle at 72 °C for 2 min. The PCR products were separated on an 8.0% non-denaturing polyacrylamide gel electrophoresis (PAGE) gel and then visualized by silver staining. DL2000 Marker DNA ladder (TsingKe, Biological Technology Co., Ltd., Nanjing, China) was used as the standard size marker.
The chromosomal location of InDel markers was acquired from the cowpea genome database (https://phytozome.jgi.doe.gov) as a reference genome, and the InDel markers were mapped onto chromosomes using MapDraw .
PCR products were scored manually, and a 0/1 binary matrix was set according to the presence or absence of corresponding amplified bands. Genotypic genetic diversity analysis used PowerMarker V3.25 (http://www.powermarker.net)  to obtain allele frequency, allele number, gene diversity index, and genotype polymorphism information content (PIC).
Availability of data and materials
We have deposited our data in Sequence Read Archive (SRA) (http://www.ncbi.nlm.nih.gov/sra/), the accession number for our submissions are: PRJNA388018.
Arabidopsis K+ transporter1
High sodium affinity transporter
Plasma membrane proteins
Apetala2/ethylene responsive factor
Dehydration responsive element binding protein
Basic leucine zipper domain
Genome-wide association study
Suzi 41 salt-stressed
Sujiang 1419 salt-stressed
Suzi 41 control
Sujiang 1419 control
Polymorphic information content
Fragments per kilo-base of exon per million fragments
Quantitative real-time PCR
National Infrastructure for Vegetable Crop Germplasm Resources
Jiangsu Academy of Agricultural Sciences
Polyacrylamide gel electrophoresis
Polymorphism information content
Gillaspie AJ, Hopkins M, Dean R. Determining genetic diversity between lines of Vigna unguiculata subspecies by AFLP and SSR markers. Genet Resour Crop Evol. 2005;52(3):245–7. https://doi.org/10.1007/s10722-004-6693-9.
Ehlers JD, Hall AE. Cowpea (Vigna unguiculata L. Walp). Field Crops Res. 1997;53(1-3):187–204. https://doi.org/10.1016/S0378-4290(97)00031-2.
Xu P, Wu XH, Wang BG, Liu YH, Qin DH, et al. Development and polymorphism of Vigna unguiculata ssp. unguiculata microsatellite markers used for phylogenetic analysis in asparagus bean (Vigna unguiculata ssp. sesquipedialis (L.) Verdc.). Mol Breeding. 2010;25:675–84.
Zhu JK. Salt and drought stress signal transduction in plants. Annu Rev Plant Biol. 2002;53(1):247–73. https://doi.org/10.1146/annurev.arplant.53.091401.143329.
Mahajan S, Pandey GK, Tuteja N. Calcium- and salt-stress signaling in plants: shedding light on SOS pathway. Arch Biochem Biophys. 2008;471(2):146–58. https://doi.org/10.1016/j.abb.2008.01.010.
Qiu QS, Guo Y, Dietrich MA, Schumaker KS, Zhu JK. Regulation of SOS1, a plasma membrane Na+/H+ exchanger in Arabidopsis thaliana, by SOS2 and SOS3. Proc Natl Acad Sci. 2002;99(12):8436–41. https://doi.org/10.1073/pnas.122224699.
Yang YQ, Guo Y. Elucidating the molecular mechanisms mediating plant salt-stress responses. New Phytol. 2018;217(2):523–39. https://doi.org/10.1111/nph.14920.
Zhu Q, Zhang JT, Gao XS, Tong JH, Xiao LT, Li WB, et al. The Arabidopsis AP2/ERF transcription factor RAP2.6 participates in ABA, salt and osmotic stress responses. Gene. 2010;457(1-2):1–12. https://doi.org/10.1016/j.gene.2010.02.011.
Niu X, Luo TL, Zhao HY, Su YL, Ji WQ, Li HF. Identification of wheat DREB genes and functional characterization of TaDREB3 in response to abiotic stresses. Gene. 2020;740:144514.
Wang LF, Zhu JF, Li XM, Wang SM, Wu J. Salt and drought stress and ABA responses related to bZIP genes from V radiata and V angularis. Gene. 2018;651:152–60. https://doi.org/10.1016/j.gene.2018.02.005.
Shi WY, Du YT MJ, Min DH, Jin LG, Chen J, et al. The WRKY transcription factor GmWRKY12 confers drought and salt tolerance in soybean. Int J Mol Sci. 2018;19(12):4087.
Wu JD, Jiang YL, Liang YN, Chen L, Chen WJ, Cheng BJ. Expression of the maize MYB transcription factor ZmMYB3R enhances drought and salt stress tolerance in transgenic plants. Plant Physiol and Bioch. 2019;137:179–88. https://doi.org/10.1016/j.plaphy.2019.02.010.
Chinnusamy V, Jagendorf A, Zhu JK. Understanding and improving salt tolerance in plants. Crop Sci. 2005;45(2):437–48. https://doi.org/10.2135/cropsci2005.0437.
Munnik T, Meijer HJ. Osmotic stress activates distinct lipid and MAPK signaling pathways in plants. FEBS Lett. 2001;498(2-3):172–8. https://doi.org/10.1016/S0014-5793(01)02492-9.
Kultz D, Burg M. Evolution of osmotic stress signaling via MAP kinase cascades. J Exp Biol. 1998;201(22):3015–21. https://doi.org/10.1242/jeb.201.22.3015.
Miller G, Shulaev V, Mittler R. Reactive oxygen signaling and abiotic stress. Physiol Plant. 2008;133(3):481–9. https://doi.org/10.1111/j.1399-3054.2008.01090.x.
Huang CH, He WL, Guo JK, Chang XX, Su PX, Zhang LX. Increased sensitivity to salt stress in an ascorbate-deficient Arabidopsis mutant. J Exp Bot. 2005;56(422):3041–9. https://doi.org/10.1093/jxb/eri301.
Bandehagh A, Salekdeh GH, Toorchi M, Mohammadi A, Komatsu S. Comparative proteomic analysis of canola leavesunder salinity stress. Proteomics. 2011;11(10):1965–75. https://doi.org/10.1002/pmic.201000564.
Chen CY, Tao CX, Peng H, Ding Y. Genetic analysis of salt stress responses in asparagus bean (Vigna unguiculata (L.) ssp. sesquipedalis Verdc.). J Hered. 2007;98(7):655–65. https://doi.org/10.1093/jhered/esm084.
Murillo-Amador B, Troyo-Die’guez E, Garcia-Hernandez JL, Lopez-Aguilar RA, Vila-serrano NY, Zamora-Salgado S, et al. Effect of salinity in the genotype variation of cowpea (Vigna unguiculata) during early vegetative growth. Sci Hortic. 2006;108(4):423–31. https://doi.org/10.1016/j.scienta.2006.02.010.
Xu P, Moshelion M, Wu XH, Halperin O, Wang BG, Luo J, et al. Natural variation and gene regulatory basis for the responses of asparagus beans to soil drought. Front Plant Sci. 2015;6:1–14. https://doi.org/10.3389/fpls.2015.00891.
Tan HQ, Huang HT, Tie MM, Tang Y, Lai YS, Li HX, et al. Transcriptome profiling of two asparagus bean (Vigna unguiculata subsp. sesquipedalis) cultivars differing in chilling tolerance under cold stress. PLoS One. 2016;11(3):e0151105.
Pan L, Yu XL, Shao JJ, Liu ZC, Gao T, Zheng Y, et al. Transcriptomic profiling and analysis of differentially expressed genes in asparagus bean (Vigna unguiculata ssp. sesquipedalis) under salt stress. PLoS One. 2019;14:1–23.
Zhang H, Xu W, Chen H, Chen J, Chen X, Yang S. Evaluation and QTL mapping of salt tolerance in yardlong bean [Vigna unguiculata (L.) Walp. Subsp. unguiculata Sesquipedalis group] seedlings. Plant Mol Biol Rep. 2020;38(2):294–304. https://doi.org/10.1007/s11105-020-01194-2.
Chen JB, Somta P, Chen X, Cui XY, Yuan XX, Srinives P. Gene mapping of a mutant mungbean (Vigna radiata L.) using new molecular markers suggests a gene encoding a YUC4-like protein regulates the chasmogamous flower trait. Front. Plant Sci. 2016;7:830. https://doi.org/10.3389/fpls.2016.00830.
Moghaddam SM, Song Q, Mamidi S, Schmutz J, Lee R, Cregan P, et al. Developing market class specific InDel markers from next generation sequence data in Phaseolus vulgaris L. Front Plant Sci. 2014;5:1–13. https://doi.org/10.3389/fpls.2014.00185.
Yundaeng C, Somta P, Chen J, Yuan X, Chankaew S, Chen X. Fine mapping of QTL conferring Cercospora leaf spot disease resistance in mungbean revealed TAF5 as candidate gene for the resistance. Theor Appl Genet. 2020;134(2):1–14. https://doi.org/10.1007/s00122-020-03724-8.
Zhai RR, Feng Y, Wang HM, Zhan XD, Shen XH, Wu WM, et al. Transcriptome analysis of rice root heterosis by RNA-Seq. BMC Genomics. 2013;14(1):19. https://doi.org/10.1186/1471-2164-14-19.
Wang B, Zhong ZH, Zhang HH, Wang X, Liu BL, Yang LJ, et al. Targeted mutagenesis of NAC transcription factor gene, OsNAC041, leading to salt sensitivity in rice. Rice Sci. 2019;26(2):98–108. https://doi.org/10.1016/j.rsci.2018.12.005.
Wu DD, Sun YH, Wang HF, Shi H, Su MX, Shan HY, et al. The SlNAC8 gene of the halophyte Suaeda liaotungensis enhances drought and salt stress tolerance in transgenic Arabidopsis thaliana. Gene. 2018;662:10–20. https://doi.org/10.1016/j.gene.2018.04.012.
Yu YH, Ni ZY, Chen QJ, Qu YY. The wheat salinity-induced R2R3-MYB transcription factor TaSIM confers salt stress tolerance in Arabidopsis thaliana. Biochem Bioph Res Co. 2017;491(3):642–8. https://doi.org/10.1016/j.bbrc.2017.07.150.
Zhu N, Cheng SF, Liu XY, Du H, Dai MQ, Zhou DX, et al. The R2R3-type MYB gene OsMYB91 has a function in coordinating plant growth and salt stress tolerance in rice. Plant Sci. 2015;236:146–56. https://doi.org/10.1016/j.plantsci.2015.03.023.
Fracasso A, Trindade LM, Amaducci S. Drought stress tolerance strategies revealed by RNA-Seq in two sorghum genotypes with contrasting WUE. BMC Plant Biol. 2016;16(1):115. https://doi.org/10.1186/s12870-016-0800-x.
Yang ZM, Lu RK, Dai ZG, Yan A, Tang Q, Cheng CH, et al. Salt-stress response mechanisms using de Novo transcriptome sequencing of salt-tolerant and sensitive Corchorus spp genotypes. Genes. 2017;8:226.
Gupta K, Sengupta A, Chakraborty M, Gupta B. Hydrogen peroxide and polyamines act as double edged swords in plant abiotic stress responses. Front Plant Sci. 2016;7:1343. https://doi.org/10.3389/fpls.2016.01343.
Farooq MA, Niazi AK, Akhtar J, Farooq M, Souri Z, Karimi N, et al. Acquiring control: the evolution of ROS-induced oxidative stress and redox signaling pathways in plant stress responses. Plant Physiol Bioch. 2019;141:353–69. https://doi.org/10.1016/j.plaphy.2019.04.039.
Bailey-Serres J, Voesenek LACJ. Flooding stress: acclimations and genetic diversity. Annu Rev Plant Biol. 2008;59(1):313–39. https://doi.org/10.1146/annurev.arplant.59.032607.092752.
Xu X, Wang H, Qi X, Xu Q, Chen X. Waterlogging-induced increase in fermentation and related gene expression in the root of cucumber (Cucumis sativus L.). Sci Hortic. 2014;179:388–95.
Cao YR, Chen SY, Zhang JS. Ethylene signaling regulates salt stress response: an overview. Plant Signal Behav. 2008;3(10):761–3. https://doi.org/10.4161/psb.3.10.5934.
Quan RD, Wang J, Yang DX, Zhang HW, Zhang ZJ, Huang RF. EIN3 and SOS2 synergistically modulate plant salt tolerance. Sci Rep. 2017;7(1):44637. https://doi.org/10.1038/srep44637.
Xu ZS, Xia LQ, Chen M, Cheng XG, Zhang RY, Li LC, et al. Isolation and molecular characterization of the Triticum aestivum L. ethylene-responsive factor 1 (TaERF1) that increases multiple stress tolerance. Plant Mol Biol. 2007;65(6):719–32. https://doi.org/10.1007/s11103-007-9237-9.
Xie ZL, Nolan TM, Jiang H, Yin YH. AP2/ERF transcription factor regulatory networks in hormone and abiotic stress responses in Arabidopsis. Front Plant Sci. 2019;10:228. https://doi.org/10.3389/fpls.2019.00228.
Zhang GY, Chen M, Li LC, Xu ZS, Chen XP, Guo JM, et al. Overexpression of the soybean GmERF3 gene, an AP2/ERF type transcription factor for increased tolerances to salt, drought, and diseases in transgenic tobacco. J Exp Bot. 2009;60(13):3781–96. https://doi.org/10.1093/jxb/erp214.
Martinoia E, Klein M, Geisler M, Bovet L, Forestier C, Kolukisaoglu U, et al. Multifunctionality of plant ABC transporters - more than just detoxifiers. Planta. 2002;214(3):345–55. https://doi.org/10.1007/s004250100661.
Geisler M, Murphy AS. The ABC of auxin transport: the role of p-glycoproteins in plant development. FEBS Lett. 2006;580(4):1094–102. https://doi.org/10.1016/j.febslet.2005.11.054.
Lane TS, Rempe CS, Davitt J, Staton ME, Peng YH, Soltis DE, et al. Diversity of ABC transporter genes across the plant kingdom and their potential utility in biotechnology. BMC Biotechnol. 2016;16(1):47. https://doi.org/10.1186/s12896-016-0277-6.
Bansal UK, Hayden MJ, Venkata BP, Khanna R, Saini RG, Bariana HS. Genetic mapping of adult plant leaf rust resistance genes Lr48 and Lr49 in common wheat. Theor Appl Genet. 2008;117(3):307–12. https://doi.org/10.1007/s00122-008-0775-6.
Dhariwal R, Gahlaut V, Govindraj BR, Singh D, Mathur S, Vyas S, et al. Stage-specific reprogramming of gene expression characterizes Lr48-mediated adult plant leaf rust resistance in wheat. Funct Integr Genomics. 2015;15(2):233–45. https://doi.org/10.1007/s10142-014-0416-x.
Saini RG, Kaur M, Singh B, Sharma S, Nanda GS, Nayar SK, et al. Genes Lr48 and Lr49 for hypersensitive adult plant leaf rust resistance in wheat (Triticum aestivum L.). Euphytica. 2002;124(3):365–70. https://doi.org/10.1023/A:1015762812907.
Liang QY, Wu YH, Wang K, Bai ZY, Liu QL, Pan YZ, et al. Chrysanthemum WRKY gene DgWRKY5 enhances tolerance to salt stress in transgenic chrysanthemum. Sci Rep. 2017;7(1):4799. https://doi.org/10.1038/s41598-017-05170-x.
Yousfi FE, Makhloufi E, Marande W, Ghorbel AW, Bouzayen M, Berges H. Comparative analysis of WRKY genes potentially involved in salt stress responses in Triticum turgidum L ssp durum. Front Plant Sci. 2016;7:2034. https://doi.org/10.3389/fpls.2016.02034.
Zheng KJ, Wang YT, Wang SC. The non-DNA binding bHLH transcription factor Paclobutrazol resistances are involved in the regulation of ABA and salt responses in Arabidopsis. Plant Physiol Biochem. 2019;139:239–45. https://doi.org/10.1016/j.plaphy.2019.03.026.
Hu Y, Chen L, Wang H, Zhang L, Wang F, Yu D. Arabidopsis transcription factor WRKY8 functions antagonistically with its interacting partner VQ9 to modulate salinity stress tolerance. Plant J. 2013;74(5):730–45. https://doi.org/10.1111/tpj.12159.
Liu YJ, Ji XY, Nie XG, Qu M, Zheng L, Tan ZL, et al. Arabidopsis AtbHLH112 regulates the expression of genes involved in abiotic stress tolerance by binding to their E-box and GCG-box motifs. New Phytol. 2015;207(3):692–709. https://doi.org/10.1111/nph.13387.
Waseem M, Rong X, Li Z. Dissecting the role of a basic Helix-loop-Helix transcription factor, SlbHLH22, under salt and drought stresses in transgenic Solanum lycopersicum L. Front Plant Sci. 2019;10:734. https://doi.org/10.3389/fpls.2019.00734.
Asare AT, Gowda BS, Galyuon IKA, Aboagye LL, Takrama JF, Timko MP. Assessment of the genetic diversity in cowpea (Vigna unguiculata L. Walp.) germplasm from Ghana using simplesequence repeat markers. Plant Genet Resour-C. 2010;8(2):142–50. https://doi.org/10.1017/S1479262110000092.
Xiong HZ, Shi AN, Mou BQ, Qin J, Motes D, Lu WG, et al. Genetic diversity and population structure of cowpea (Vigna unguiculata L. Walp). PLoS One. 2016;11(8):e0160941. https://doi.org/10.1371/journal.pone.0160941.
Hoagland DR, Arnon DI. The water-culture method for growing plants without soil. Calif Agric Exp Stn Circ. 1950;347:1–32.
Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and cufflinks. Nat Protoc. 2012;7(3):562–78. https://doi.org/10.1038/nprot.2012.016.
Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28(5):511–5. https://doi.org/10.1038/nbt.1621.
Mao XZ, Cai T, Olyarchuk JG, Wei LP. Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics. 2005;21(19):3787–93. https://doi.org/10.1093/bioinformatics/bti430.
Wang LK, Feng ZX, Wang X, Wang XW, Zhang XG. DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010;26(1):136–8. https://doi.org/10.1093/bioinformatics/btp612.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2-△△CT method. Methods. 2001;25(4):402–8. https://doi.org/10.1006/meth.2001.1262.
Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25(9):1105–11. https://doi.org/10.1093/bioinformatics/btp120.
McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303. https://doi.org/10.1101/gr.107524.110.
Rozen S, Skaletsky H. Primer3 on the www for general users and for biologist programme. Methods Mol Biol. 2000;132:365–86. https://doi.org/10.1385/1-59259-192-2:365.
Beier S, Thiel T, Munch T, Scholz U, Mascher M. MISA-web: a web server for microsatellite prediction. Bioinformatics. 2017;33(16):2583–5. https://doi.org/10.1093/bioinformatics/btx198.
Englen MD, Kelley LC. A rapid DNA isolation procedure for the identification of campylobacter jejuni by the polymerase chain reaction. Lett Appl Microbiol. 2000;31(6):421–6. https://doi.org/10.1046/j.1365-2672.2000.00841.x.
Liu RH, Meng JL. MapDraw: a microsoft excel macro for drawing genetic linkage maps based on given genetic linkage data. Hereditas. 2003;25(3):317–21.
Liu KJ, Muse SV. PowerMarker: an integrated analysis environment for genetic marker analysis. Bioinformatics. 2005;21(9):2128–9. https://doi.org/10.1093/bioinformatics/bti282.
We are grateful to Shan Meng and Wei Zhang for their suggestions on experimental design.
This work was supported by Agriculture Research System of China (CARS-09) and Jiangsu Agricultural Science and Technology Innovation Fund (CX(20)3161). The funding bodies had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Number of reads obtained by RNA sequencing and their matches in the Vigna unguiculata genome.
Summary of novel transcript isoforms found in V. unguiculata genomic data.
Differentially expressed genes (DEGs) in Suzi 41 and Sujiang 1419 under salt stress.
Significant DEGs exclusive to Suzi 41 and Sujiang 1419 during salt stress.
Differentially expressed genes that are up- or down-regulated in both Suzi 41 and Sujiang 1419 during salt stress.
GO analysis of DEGs induced by salt stress in Suzi 41 and Sujiang 1419.
KEGG pathways enriched with DEGs under high salinity in Suzi 41 and Sujiang 1419.
TF genes identified as DEGs in Suzi 41 and Sujiang 1419 under salt stress.
TFs that are up- or down-regulated in both Suzi 41 and Sujiang 1419.
InDels identified on individual chromosomes of yardlong bean.
The number and distribution ratios of InDels identified in yardlong bean.
Characteristics of 175 InDels derived from RNA-seq data of salt-stressed yardlong bean.
List of 84 yardlong bean germplasm accessions used in this study.
Genetic diversity of yardlong bean based on 26 InDel markers.
Primer used for qRT-PCR analysis.
About this article
Cite this article
Zhang, H., Xu, W., Chen, H. et al. Transcriptomic analysis of salt tolerance-associated genes and diversity analysis using indel markers in yardlong bean (Vigna unguiculata ssp. sesquipedialis). BMC Genom Data 22, 34 (2021). https://doi.org/10.1186/s12863-021-00989-w
- LR48 transcription factor
- Salt stress
- Yardlong bean