Sesame is an important oil crop due to its high oil, antioxidant, and protein content. Drought stress is a major abiotic stress that affects sesame production as well as the quality of sesame seed. To reveal the adaptive mechanism of sesame in response to water deficient conditions, transcriptomic and metabolomics were applied in drought-tolerant (DT) and drought-susceptible (DS) sesame genotypes. Transcriptomic analysis reveals a set of core drought-responsive genes (684 up-regulated and 1346 down-regulated) in sesame that was robustly differently expressed in both genotypes. Most enriched drought-responsive genes are mainly involved in protein processing in endoplasmic reticulum, plant hormone signal transduction photosynthesis, lipid metabolism, and amino acid metabolism. Drought-susceptible genotype was more disturbed by drought stress at both transcriptional and metabolic levels, since more drought-responsive genes/metabolites were identified in DS. Drought-responsive genes associated with stress response, amino acid metabolism, and reactive oxygen species scavenging were more enriched or activated in DT. According to the partial least-squares discriminate analysis, the most important metabolites which were accumulated under drought stress in both genotypes includes ABA, amino acids, and organic acids. Especially, higher levels of ABA, proline, arginine, lysine, aromatic and branched chain amino acids, GABA, saccharopine, 2-aminoadipate, and allantoin were found in DT under stress condition. Combination of transcriptomic and metabolomic analysis highlights the important role of amino acid metabolism (especially saccharopine pathway) and ABA metabolism and signaling pathway for drought tolerance in sesame. The results of the present study provide valuable information for better understanding the molecular mechanism underlying drought tolerance of sesame, and also provide useful clues for the genetic improvement of drought tolerance in sesame.

Sesame (Sesamum indicum L.) is one of the world’s most important oil crops. However, it is susceptible to abiotic stresses in general, and to waterlogging and drought stresses in particular. The molecular mechanisms of abiotic stress tolerance in sesame have not yet been elucidated. The WRKY domain transcription factors play significant roles in plant growth, development, and responses to stresses. However, little is known about the number, location, structure, molecular phylogenetics, and expression of the WRKY genes in sesame. We performed a comprehensive study of the WRKY gene family in sesame and identified 71 SiWRKYs. In total, 65 of these genes were mapped to 15 linkage groups within the sesame genome. A phylogenetic analysis was performed using a related species (Arabidopsis thaliana) to investigate the evolution of the sesame WRKY genes. Tissue expression profiles of the WRKY genes demonstrated that six SiWRKY genes were highly expressed in all organs, suggesting that these genes may be important for plant growth and organ development in sesame. Analysis of the SiWRKY gene expression patterns revealed that 33 and 26 SiWRKYs respond strongly to waterlogging and drought stresses, respectively. Changes in the expression of 12 SiWRKY genes were observed at different times after the waterlogging and drought treatments had begun, demonstrating that sesame gene expression patterns vary in response to abiotic stresses. In this study, we analyzed the WRKY family of transcription factors encoded by the sesame genome. Insight was gained into the classification, evolution, and function of the SiWRKY genes, revealing their putative roles in a variety of tissues. Responses to abiotic stresses in different sesame cultivars were also investigated. The results of our study provide a better understanding of the structures and functions of sesame WRKY genes and suggest that manipulating these WRKYs could enhance resistance to waterlogging and drought.

Soil salinity is one of the major serious factors that affect agricultural productivity of almost all crops worldwide, including the important oilseed crop sesame. In order to improve salinity resistance in sesame, it is crucial to understand the molecular mechanisms underlying the adaptive response to salinity stress. In the present study, two contrasting sesame genotypes differing in salt tolerance were used to decipher the adaptive responses to salt stress based on morphological, transcriptome and metabolome characterizations. Morphological results indicated that under salt stress, the salt-tolerant (ST) genotype has enhanced capacity to withstand salinity stress, higher seed germination rate and plant survival rate, as well as better growth rate than the salt-sensitive genotype. Transcriptome analysis revealed strongly induced salt-responsive genes in sesame mainly related to amino acid metabolism, carbohydrate metabolism, biosynthesis of secondary metabolites, plant hormone signal transduction, and oxidation-reduction process. Especially, several pathways were preferably enriched with differentially expressed genes in ST genotype, including alanine, aspartate and glutamate metabolism, carotenoid biosynthesis, galactose metabolism, glycolysis/gluconeogenesis, glyoxylate and dicarboxylate metabolism, porphyrin and chlorophyll metabolism. Metabolome profiling under salt stress showed a higher accumulation degree of metabolites involved in stress tolerance in ST, and further highlighted that the amino acid metabolism, and sucrose and raffinose family oligosaccharides metabolism were enhanced in ST. These findings suggest that the candidate genes and metabolites involved in crucial biological pathways may regulate salt tolerance of sesame, and increase our understanding of the molecular mechanisms underlying the adaptation of sesame to salt stress.

Glycosyltransferases comprise a highly divergent and polyphyletic multigene family that is involved in widespread modification of plant secondary metabolites in a process called glycosylation. According to conserved domains identified in their amino acid sequences, these glycosyltransferases can be classified into a single UDP-glycosyltransferase (UGT) 1 superfamily. We performed genome-wide comparative analysis of UGT genes to trace evolutionary history in algae, bryophytes, pteridophytes, and angiosperms; then, we further investigated the expansion mechanisms and function characterization of UGT gene families in Brassica rapa and Brassica oleracea. Using Hidden Markov Model search, we identified 3, 21, 140, 200, 115, 147, and 147 UGTs in Chlamydomonas reinhardtii, Physcomitrella patens, Selaginella moellendorffii, Oryza sativa, Arabidopsis thaliana, B. rapa, and B. oleracea, respectively. Phylogenetic analysis revealed that UGT80 gene family is an ancient gene family, which is shared by all plants and UGT74 gene family is shared by ferns and angiosperms, but the remaining UGT gene families were shared by angiosperms. In dicot lineage, UGTs among three species were classified into three subgroups containing 3, 6, and 12 UGT gene families. Analysis of chromosomal distribution indicates that 98.6 and 71.4% of UGTs were located on B. rapa and B. oleracea pseudo-molecules, respectively. Expansion mechanism analyses uncovered that whole genome duplication event exerted larger influence than tandem duplication on expansion of UGT gene families in B. rapa, and B. oleracea. Analysis of selection forces of UGT orthologous gene pairs in B. rapa, and B. oleracea compared to A. thaliana suggested that orthologous genes in B. rapa, and B. oleracea have undergone negative selection, but there were no significant differences between A. thaliana –B. rapa and A. thaliana –B. oleracea lineages. Our comparisons of expression profiling illustrated that UGTs in B. rapa performed more discrete expression patterns than these in B. oleracea indicating stronger function divergence. Combing with phylogeny and expression analysis, the UGTs in B. rapa and B. oleracea experienced parallel evolution after they diverged from a common ancestor. We first traced the evolutionary history of UGT gene families in plants and revealed its evolutionary and functional characterization of UGTs in B. rapa, and B. oleracea. This study provides novel insights into the evolutionary history and functional divergence of important traits or phenotype-related gene families in plants.

The cytochrome P450 monooxygenase (P450) superfamily is involved in the biosynthesis of various primary and secondary metabolites. However, little is known about the effects of whole genome duplication (WGD) and tandem duplication (TD) events on the evolutionary history and functional divergence of P450s in Brassica after splitting from a common ancestor with Arabidopsis thaliana. Using Hidden Markov Model search and manual curation, we detected that Brassica species have nearly 1.4-fold as many P450 members as A. thaliana. Most P450s in A. thaliana and Brassica species were located on pseudo-chromosomes. The inferred phylogeny indicated that all P450s were clustered into two different subgroups. Analysis of WGD event revealed that different P450 gene families had appeared after evolutionary events of species. For the TD event analyses, the P450s from TD events in Brassica species can be divided into ancient and recent parts. Our comparison of influence of WGD and TD events on the P450 gene superfamily between A. thaliana and Brassica species indicated that the family-specific evolution in the Brassica lineage can be attributed to both WGD and TD, whereas WGD was recognized as the major mechanism for the recent evolution of the P450 super gene family. Expression analysis of P450s from A. thaliana and Brassica species indicated that WGD-type P450s showed the same expression pattern but completely different expression with TD-type P450s across different tissues in Brassica species. Selection force analysis suggested that P450 orthologous gene pairs between A. thaliana and Brassica species underwent negative selection, but no significant differences were found between P450 orthologous gene pairs in A. thaliana–B. rapa and A. thaliana–B. oleracea lineages, as well as in different subgenomes in B. rapa or B. oleracea compared with A. thaliana. This study is the first to investigate the effects of WGD and TD on the evolutionary history and functional divergence of P450 gene families in A. thaliana and Brassica species. This study provides a biology model to study the mechanism of gene family formation, particularly in the context of the evolutionary history of angiosperms, and offers novel insights for the study of angiosperm genomes.

Sesame is poised to become a major oilseed crop owing to its high oil quality and adaptation to various ecological areas. However, the seed yield of sesame is very low and the underlying genetic basis is still elusive. Here, we performed genome-wide association studies of 39 seed yield-related traits categorized into five major trait groups, in three different environments, using 705 diverse lines. Extensive variation was observed for the traits with capsule size, capsule number and seed size-related traits, found to be highly correlated with seed yield indexes. In total, 646 loci were significantly associated with the 39 traits (p < 10⁻⁷) and resolved to 547 quantitative trait loci QTLs. We identified six multi-environment QTLs and 76 pleiotropic QTLs associated with two to five different traits. By analyzing the candidate genes for the assayed traits, we retrieved 48 potential genes containing significant functional loci. Several homologs of these candidate genes in Arabidopsis are described to be involved in seed or biomass formation. However, we also identified novel candidate genes, such as SiLPT3 and SiACS8, which may control capsule length and capsule number traits. Altogether, we provided the highly-anticipated basis for research on genetics and functional genomics towards seed yield improvement in sesame.