Background: Sesame is an important oil crop, but limited transcriptomic and genomic data are currently available. This information is essential to clarify the fatty acid and lignan biosynthesis molecular mechanism. In addition, a shortage of sesame molecular markers limits the efficiency and accuracy of genetic breeding. High-throughput transcriptomic sequencing is essential to generate a large transcriptome sequence dataset for gene discovery and molecular marker development. Results: Sesame transcriptomes from five tissues were sequenced using Illumina paired-end sequencing technology. The cleaned raw reads were assembled into a total of 86,222 unigenes with an average length of 629 bp. Of the unigenes, 46,584 (54.03%) had significant similarity with proteins in the NCBI nonredundant protein database and Swiss-Prot database (E-value < 10^-5). Of these annotated unigenes, 10,805 and 27,588 unigenes were assigned to gene ontology categories and clusters of orthologous groups, respectively. In total, 22,003 (25.52%) unigenes were mapped onto 119 pathways using the Kyoto Encyclopedia of Genes and Genomes Pathway database (KEGG). Furthermore, 44,750 unigenes showed homology to 15,460 Arabidopsis genes based on BLASTx analysis against The Arabidopsis Information Resource (TAIR, Version 10) and revealed relatively high gene coverage. In total, 7,702 unigenes were converted into SSR markers (EST-SSR). Dinucleotide SSRs were the dominant repeat motif (67.07%, 5,166), followed by trinucleotide (24.89%, 1,917), tetranucleotide (4.31%, 332), hexanucleotide (2.62%, 202), and pentanucleotide (1.10%, 85) SSRs. AG/CT (46.29%) was the dominant repeat motif, followed by AC/GT (16.07%), AT/AT (10.53%), AAG/CTT (6.23%), and AGG/CCT (3.39%). Fifty EST-SSRs were randomly selected to validate amplification and to determine the degree of polymorphism in the genomic DNA pools. Forty primer pairs successfully amplified DNA fragments and detected significant amounts of polymorphism among 24 sesame accessions. Conclusions: This study demonstrates that Illumina paired-end sequencing is a fast and cost-effective approach to gene discovery and molecular marker development in non-model organisms. Our results provide a comprehensive sequence resource for sesame research.

Abiotic stresses, such as drought and salt, are major environmental stresses, affecting plant growth and crop productivity. Plant bZIP transcription factors (bZIPs) confer stress resistances in harsh environments and play important roles in each phase of plant growth processes. In this research, 15 soybean bZIP family members were identified from drought-induced de novo transcriptomic sequences of soybean, which were unevenly distributed across 12 soybean chromosomes. Promoter analysis showed that these 15 genes were rich in ABRE, MYB and MYC cis-acting elements which were reported to be involved in abiotic stress responses. Quantitative real-time polymerase chain reaction (qRT-PCR) analysis indicated that 15 GmbZIP genes could be induced by drought and salt stress. GmbZIP2 was significantly upregulated under stress conditions and thus was selected for further study. Subcellular localization analysis revealed that the GmbZIP2 protein was located in the cell nucleus. qRT-PCR results show that GmbZIP2 can be induced by multiple stresses. The overexpression of GmbZIP2 in Arabidopsis and soybean hairy roots could improve plant resistance to drought and salt stresses. The result of differential expression gene analysis shows that the overexpression of GmbZIP2 in soybean hairy roots could enhance the expression of the stress responsive genes GmMYB48, GmWD40, GmDHN15, GmGST1 and GmLEA. These results indicate that soybean bZIPs played pivotal roles in plant resistance to abiotic stresses.

GRAS genes, which form a plant-specific transcription factor family, play an important role in plant growth and development and stress responses. However, the functions of GRAS genes in soybean (Glycine max) remain largely unknown. Here, 117 GRAS genes distributed on 20 chromosomes were identified in the soybean genome and were classified into 11 subfamilies. Of the soybean GRAS genes, 80.34% did not have intron insertions, and 54 pairs of genes accounted for 88.52% of duplication events (61 pairs). RNA-seq analysis demonstrated that most GmGRASs were expressed in 14 different soybean tissues examined and responded to multiple abiotic stresses. Results from quantitative real-time PCR analysis of six selected GmGRASs suggested that GmGRAS37 was significantly upregulated under drought and salt stress conditions and abscisic acid and brassinosteroid treatment; therefore, this gene was selected for further study. Subcellular localization analysis revealed that the GmGRAS37 protein was located in the plasma membrane, nucleus, and cytosol. Soybean hairy roots overexpressing GmGRAS37 had improved resistance to drought and salt stresses. In addition, these roots showed increased transcript levels of several drought‐ and salt-related genes. The results of this study provide the basis for comprehensive analysis of GRAS genes and insight into the abiotic stress response mechanism in soybean.

Sesame (Sesamum indicum L.) is one of the four major oil crops in China. A sesame core collection (CC) was established in China in 2000, but no complete study on its genetic diversity has been carried out at either the phenotypic or molecular level. To provide technical guidance, a theoretical basis for further collection, effective protection, reasonable application, and a complete analysis of sesame genetic resources, a genetic diversity assessment of the sesame CC in China was conducted using phenotypic and molecular data and by extracting a sesame mini-core collection (MC). Results from a genetic diversity assessment of sesame CC in China were significantly inconsistent at the phenotypic and molecular levels. A Mantel test revealed the insignificant correlation between phenotype and molecular marker information (r = 0.0043, t = 0.1320, P = 0.5525). The Shannon-Weaver diversity index (I) and Nei genetic diversity index (h) were higher (I = 0.9537, h = 0.5490) when calculated using phenotypic data from the CC than when using molecular data (I = 0.3467, h = 0.2218). A mini-core collection (MC) containing 184 accessions was extracted based on both phenotypic and molecular data, with a low mean difference percentage (MD, 1.64%), low variance difference percentage (VD, 22.58%), large variable rate of coefficient of variance (VR, 114.86%), and large coincidence rate of range (CR, 95.76%). For molecular data, the diversity indices and the polymorphism information content (PIC) for the MC were significantly higher than for the CC. Compared to an alternative random sampling strategy, the advantages of capturing genetic diversity and validation by extracting a MC using an advanced maximization strategy were proven. This study provides a comprehensive characterization of the phenotypic and molecular genetic diversities of the sesame CC in China. A MC was extracted using both phenotypic and molecular data. Low MD% and VD%, and large VR% and CR% suggested that the MC provides a good representation of the genetic diversity of the original CC. The MC was more genetically diverse with higher diversity indices and a higher PIC value than the CC. A MC may aid in reasonably and efficiently selecting materials for sesame breeding and for genotypic biological studies, and may also be used as a population for association mapping in sesame.

Abiotic stresses, such as drought and salinity, severely affects the growth, development and productivity of the plants. The Catharanthus roseus RLK1-like (CrRLK1L) protein kinase family is involved in several processes in the plant life cycle. However, there have been few studies addressing the functions of CrRLK1L proteins in soybean. In this study, 38 CrRLK1L genes were identified in the soybean genome (Glycine max Wm82.a2.v1). Phylogenetic analysis demonstrated that soybean CrRLK1L genes were grouped into clusters, cluster I, II, III. The chromosomal mapping demonstrated that 38 CrRLK1L genes were located in 14 of 20 soybean chromosomes. None were discovered on chromosomes 1, 4, 6, 7, 11, and 14. Gene structure analysis indicated that 73.6% soybean CrRLK1L genes were characterized by a lack of introns.15.7% soybean CrRLK1L genes only had one intron and 10.5% soybean CrRLK1L genes had more than one intron. Five genes were obtained from soybean drought- and salt-induced transcriptome databases and were found to be highly up-regulated. GmCrRLK1L20 was notably up-regulated under drought and salinity stresses, and was therefore studied further. Subcellular localization analysis revealed that the GmCrRLK1L20 protein was located in the cell membrane. The overexpression of the GmCrRLK1L20 gene in soybean hairy roots improved both drought tolerance and salt stresses and enhanced the expression of the stress-responsive genes GmMYB84, GmWRKY40, GmDREB-like, GmGST15, GmNAC29, and GmbZIP78. These results indicated that GmCrRLK1L20 could play a vital role in defending against drought and salinity stresses in soybean.

Vascular plant one-zinc-finger (VOZ) transcription factor, a plant specific one-zinc-finger-type transcriptional activator, is involved in regulating numerous biological processes such as floral induction and development, defense against pathogens, and response to multiple types of abiotic stress. Six VOZ transcription factor-encoding genes (GmVOZs) have been reported to exist in the soybean (Glycine max) genome. In spite of this, little information is currently available regarding GmVOZs. In this study, GmVOZs were cloned and characterized. GmVOZ genes encode proteins possessing transcriptional activation activity in yeast cells. GmVOZ1E, GmVOZ2B, and GmVOZ2D gene products were widely dispersed in the cytosol, while GmVOZ1G was primarily located in the nucleus. GmVOZs displayed a differential expression profile under dehydration, salt, and salicylic acid (SA) stress conditions. Among them, GmVOZ1G showed a significantly induced expression in response to all stress treatments. Overexpression of GmVOZ1G in soybean hairy roots resulted in a greater tolerance to drought and salt stress. In contrast, RNA interference (RNAi) soybean hairy roots suppressing GmVOZ1G were more sensitive to both of these stresses. Under drought treatment, soybean composite plants with an overexpression of hairy roots had higher relative water content (RWC). In response to drought and salt stress, lower malondialdehyde (MDA) accumulation and higher peroxidase (POD) and superoxide dismutase (SOD) activities were observed in soybean composite seedlings with an overexpression of hairy roots. The opposite results for each physiological parameter were obtained in RNAi lines. In conclusion, GmVOZ1G positively regulates drought and salt stress tolerance in soybean hairy roots. Our results will be valuable for the functional characterization of soybean VOZ transcription factors under abiotic stress.