Comparative analysis of expressed sequence tags from Sesamum indicum and Arabidopsis thaliana developing seeds

Omics technologies towards sesame improvement: a review

Genetic improvement of sesame (Sesamum indicum L.), one of the most important oilseed crops providing edible oil, proteins, minerals, and vitamins, is important to ensure a balanced diet for the growing world population. Increasing yield, seed protein, oil, minerals, and vitamins is urgently needed to meet the global demand. The production and productivity of sesame is very low due to various biotic and abiotic stresses. Therefore, various efforts have been made to combat these constraints and increase the production and productivity of sesame through conventional breeding. However, less attention has been paid to the genetic improvement of the crop through modern biotechnological methods, leaving it lagging behind other oilseed crops. Recently, however, the scenario has changed as sesame research has entered the era of "omics" and has made significant progress. Therefore, the purpose of this paper is to provide an overview of the progress made by omics research in improving sesame. This review presents a number of efforts that have been made over past decade using omics technologies to improve various traits of sesame, including seed composition, yield, and biotic and abiotic resistant varieties. It summarizes the advances in genetic improvement of sesame using omics technologies, such as germplasm development (web-based functional databases and germplasm resources), gene discovery (molecular markers and genetic linkage map construction), proteomics, transcriptomics, and metabolomics that have been carried out in the last decade. In conclusion, this review highlights future directions that may be important for omics-assisted breeding in sesame genetic improvement.

Oil plant genomes: current state of the science

Vegetable oils are an indispensable nutritional component of the human diet as well as important raw materials for a variety of industrial applications such as pharmaceuticals, cosmetics, oleochemicals, and biofuels. Oil plant genomes are highly diverse, and their genetic variation leads to a diversity in oil biosynthesis and accumulation along with agronomic traits. This review discusses plant oil biosynthetic pathways, current state of genome assembly, polyploidy and asymmetric evolution of genomes of oil plants and their wild relatives, and research progress of pan-genomics in oil plants. The availability of complete high-resolution genomes and pan-genomes has enabled the identification of structural variations in the genomes that are associated with the diversity of agronomic and environment fitness traits. These and future genomes also provide powerful tools to understand crop evolution and to harvest the rich natural variations to improve oil crops for enhanced productivity, oil quality, and adaptability to changing environments.

Characterization of the Genetic Diversity Present in a Diverse Sesame Landrace Collection Based on Phenotypic Traits and EST-SSR Markers Coupled With an HRM Analysis

A selection of sesame (Sesamum indicum L.) landraces of different eco-geographical origin and breeding history have been characterized using 28 qualitative morpho-physiological descriptors and seven expressed sequence tag-simple sequence repeat (EST-SSR) markers coupled with a high-resolution melting (HRM) analysis. The most variable qualitative traits that could efficiently discriminate landraces, as revealed by the correlation analyses, were the plant growth type and position of the branches, leaf blade width, stem pubescence, flowering initiation, capsule traits and seed coat texture. The agglomerative hierarchical clustering analysis based on a dissimilarity matrix highlighted three main groups among the sesame landraces. An EST-SSR marker analysis revealed an average polymorphism information content (PIC) value of 0.82, which indicated that the selected markers were highly polymorphic. A principal coordinate analysis and dendrogram reconstruction based on the molecular data classified the sesame genotypes into four major clades. Both the morpho-physiological and molecular analyses showed that landraces from the same geographical origin were not always grouped in the same cluster, forming heterotic groups; however, clustering patterns were observed for the Greek landraces. The selective breeding of such traits could be employed to unlock the bottleneck of local phenotypic diversity and create new cultivars with desirable traits.

High oil accumulation in tuber of yellow nutsedge compared to purple nutsedge is associated with more abundant expression of genes involved in fatty acid synthesis and triacylglycerol storage

BACKGROUND

Yellow nutsedge is a unique plant species that can accumulate up to 35% oil of tuber dry weight, perhaps the highest level observed in the tuber tissues of plant kingdom. To gain insight into the molecular mechanism that leads to high oil accumulation in yellow nutsedge, gene expression profiles of oil production pathways involved carbon metabolism, fatty acid synthesis, triacylglycerol synthesis, and triacylglycerol storage during tuber development were compared with purple nutsedge, the closest relative of yellow nutsedge that is poor in oil accumulation.

RESULTS

Compared with purple nutsedge, high oil accumulation in yellow nutsedge was associated with significant up-regulation of specific key enzymes of plastidial RubisCO bypass as well as malate and pyruvate metabolism, almost all fatty acid synthesis enzymes, and seed-like oil-body proteins. However, overall transcripts for carbon metabolism toward carbon precursor for fatty acid synthesis were comparable and for triacylglycerol synthesis were similar in both species. Two seed-like master transcription factors ABI3 and WRI1 were found to display similar transcript patterns but were expressed at 6.5- and 14.3-fold higher levels in yellow nutsedge than in purple nutsedge, respectively. A weighted gene co-expression network analysis revealed that ABI3 was in strong transcriptional coordination with WRI1 and other key oil-related genes.

CONCLUSIONS

These results implied that pyruvate availability and fatty acid synthesis in plastid, along with triacylglycerol storage in oil bodies, rather than triacylglycerol synthesis in endoplasmic reticulum, are the major factors responsible for high oil production in tuber of yellow nutsedge, and ABI3 most likely plays a critical role in regulating oil accumulation. This study is of significance with regard to understanding the molecular mechanism controlling carbon partitioning toward oil production in oil-rich tuber and provides a valuable reference for enhancing oil accumulation in non-seed tissues of crops through genetic breeding or metabolic engineering.

Identification and evaluation of anti-inflammatory properties of aqueous components extracted from sesame (Sesamum indicum) oil

We previously reported that sesame oil (SO) has anti-inflammatory, anti-atherosclerotic and lipid lowering properties in vivo. Our recent studies have shown that, an aqueous extract of sesame oil (SOAE) has also anti-inflammatory and anti-atherosclerotic properties but with no lipid lowering effects. The extent of reduction in atherosclerosis led us to identify components of SOAE and evaluate their anti-inflammatory properties in vitro. Liquid chromatography mass spectrometric method was used to detect and identify components of SOAE. Methoxyphenol derivatives, short and long chain carboxylic acids, dicarboxylic acids, hydroxy and oxo- carboxylic acids were detected. To our surprise, sesamol and its derivatives (lignans), were not present in the SOAE. Among the identified, a combination of methoxy phenol compounds were selected and tested their ability to reduce LPS induced inflammatory gene expression. Monocyte derived macrophages/RAW 264.7 macrophages were pre-treated with these compounds for 2 h, followed by LPS stimulation for 24 h and pro-inflammatory gene expressions were analyzed. These methoxyphenol derivatives showed potent anti-inflammatory properties. In conclusion, the anti-inflammatory molecules associated with SO may contribute the anti-inflammatory and anti-atherosclerotic properties. Also, our results shed light for the development of SOAE based non-pharmacological therapeutics, nutritional supplements and health products for various inflammatory diseases in the future.

The Emerging Oilseed Crop Sesamum indicum Enters the "Omics" Era

Sesame (Sesamum indicum L.) is one of the oldest oilseed crops widely grown in Africa and Asia for its high-quality nutritional seeds. It is well adapted to harsh environments and constitutes an alternative cash crop for smallholders in developing countries. Despite its economic and nutritional importance, sesame is considered as an orphan crop because it has received very little attention from science. As a consequence, it lags behind the other major oil crops as far as genetic improvement is concerned. In recent years, the scenario has considerably changed with the decoding of the sesame nuclear genome leading to the development of various genomic resources including molecular markers, comprehensive genetic maps, high-quality transcriptome assemblies, web-based functional databases and diverse daft genome sequences. The availability of these tools in association with the discovery of candidate genes and quantitative trait locis for key agronomic traits including high oil content and quality, waterlogging and drought tolerance, disease resistance, cytoplasmic male sterility, high yield, pave the way to the development of some new strategies for sesame genetic improvement. As a result, sesame has graduated from an "orphan crop" to a "genomic resource-rich crop." With the limited research teams working on sesame worldwide, more synergic efforts are needed to integrate these resources in sesame breeding for productivity upsurge, ensuring food security and improved livelihood in developing countries. This review retraces the evolution of sesame research by highlighting the recent advances in the "Omics" area and also critically discusses the future prospects for a further genetic improvement and a better expansion of this crop.

Comparative transcriptome analysis during early fruit development between three seedy citrus genotypes and their seedless mutants

Identification of genes with differential transcript abundance (GDTA) in seedless mutants may enhance understanding of seedless citrus development. Transcriptome analysis was conducted at three time points during early fruit development (Phase 1) of three seedy citrus genotypes: Fallglo (Bower citrus hybrid (Citrus reticulata×C. reticulata×C. paradisi)×Temple (C. reticulata×C. sinensis)), grapefruit (C. paradisi), Pineapple sweet orange (C. sinensis), and their seedless mutants. Seed abortion in seedless mutants was observed at 26 days post anthesis (Time point 2). Affymetrix transcriptomic analysis revealed 359 to 1077 probe sets with differential transcript abundance in the comparison of seedless versus seedy fruits for each citrus genotypes and time points. The GDTA identified by 18 microarray probe sets were validated by qPCR. Hierarchical clustering analysis revealed a range of GDTA associated with development, hormone and protein metabolism, all of which may reflect genes associated with seedless fruit development. There were 14, 9 and 12 genes found exhibiting similar abundance ratios in all three seedless versus seedy genotype comparisons at time point 1, 2 and 3, respectively. Among those genes were genes coding for an aspartic protease and a cysteine protease, which may play important roles in seedless fruit development. New insights into seedless citrus fruit development may contribute to biotech approaches to create seedless cultivars.

Tracking sesamin synthase gene expression through seed maturity in wild and cultivated sesame species--a domestication footprint

Sesamin and sesamolin are the major oil-soluble lignans present in sesame seed, having a wide range of biological functions beneficial to human health. Understanding sesame domestication history using sesamin synthase gene expression could enable delineation of the sesame putative progenitor. This report examined the functional expression of sesamin synthase (CYP81Q1) during capsule maturation (0-40 days after flowering) in three wild Sesamum species and four sesame cultivars. Among the cultivated accessions, only S. indicum (CO-1) exhibited transcript abundance of sesamin synthase along with high sesamin content similar to S. malabaricum, while the other cultivated sesame showed low expression. The sesamin synthase expression analysis, coupled with quantification of sesamin level, indicates that sesamin synthase was not positively favoured during domestication. The sesamin synthase expression pattern and lignan content, along with phylogenetic analysis suggested a close relationship of cultivated sesame and the wild species S. malabaricum. The high genetic identity between the two species S. indicum and S. malabaricum points towards the role of the putative progenitor S. malabaricum in sesame breeding programmes to broaden the genetic base of sesame cultivars. This study emphasises the need to investigate intraspecific and interspecific variation in the primary, secondary and tertiary gene pools to develop superior sesame genotypes.

Association mapping of seed oil and protein content in Sesamum indicum L. using SSR markers

Sesame is an important oil crop for the high oil content and quality. The seed oil and protein contents are two important traits in sesame. To identify the molecular markers associated with the seed oil and protein contents in sesame, we systematically performed the association mapping among 369 worldwide germplasm accessions under 5 environments using 112 polymorphic SSR markers. The general linear model (GLM) was applied with the criteria of logP≥3.0 and high stability under all 5 environments. Among the 369 sesame accessions, the oil content ranged from 27.89%–58.73% and the protein content ranged from 16.72%–27.79%. A significant negative correlation of the oil content with the protein content was found in the population. A total of 19 markers for oil content were detected with a R² value range from 4% to 29%; 24 markers for protein content were detected with a R² value range from 3% to 29%, of which 19 markers were associated with both traits. Moreover, partial markers were confirmed using mixed linear model (MLM) method, which suggested that the oil and protein contents are controlled mostly by major genes. Allele effect analysis showed that the allele associated with high oil content was always associated with low protein content, and vice versa. Of the 19 markers associated with oil content, 17 presented near the locations of the plant lipid pathway genes and 2 were located just next to a fatty acid elongation gene and a gene encoding Stearoyl-ACP Desaturase, respectively. The findings provided a valuable foundation for oil synthesis gene identification and molecular marker assistant selection (MAS) breeding in sesame.

Genome sequencing of the high oil crop sesame provides insight into oil biosynthesis

Background

Sesame, Sesamum indicum L., is considered the queen of oilseeds for its high oil content and quality, and is grown widely in tropical and subtropical areas as an important source of oil and protein. However, the molecular biology of sesame is largely unexplored.

Results

Here, we report a high-quality genome sequence of sesame assembled de novo with a contig N50 of 52.2 kb and a scaffold N50 of 2.1 Mb, containing an estimated 27,148 genes. The results reveal novel, independent whole genome duplication and the absence of the Toll/interleukin-1 receptor domain in resistance genes. Candidate genes and oil biosynthetic pathways contributing to high oil content were discovered by comparative genomic and transcriptomic analyses. These revealed the expansion of type 1 lipid transfer genes by tandem duplication, the contraction of lipid degradation genes, and the differential expression of essential genes in the triacylglycerol biosynthesis pathway, particularly in the early stage of seed development. Resequencing data in 29 sesame accessions from 12 countries suggested that the high genetic diversity of lipid-related genes might be associated with the wide variation in oil content. Additionally, the results shed light on the pivotal stage of seed development, oil accumulation and potential key genes for sesamin production, an important pharmacological constituent of sesame.

Conclusions

As an important species from the order Lamiales and a high oil crop, the sesame genome will facilitate future research on the evolution of eudicots, as well as the study of lipid biosynthesis and potential genetic improvement of sesame.

Comparative transcriptional analysis reveals differential gene expression between Sand Daffodil tissues

Sand Daffodil (Pancratium maritimum) is a world-wide endangered Amayllidaceae species and represents an important anti-cancer medicinal resource due to alkaloids production. Despite its increasing pharmaceutical importance, there are not molecular resources that can be utilized toward improving genetic traits. In our research, the suppression subtractive hybridization (SSH) method conducted to generate large-scale expressed sequence tags (EST), was designed to identify gene candidates related to the morphological and physiological differences between the two tissues, leaves and bulbs, since lycorine, the main anti-cancer compound, is there synthesized. We focused on identification of transcripts in different tissues from Sand Daffodil using PCR-based suppression SSH to identify genes involved in global pathway control. Sequencing of 2,000 differentially screened clones from the SSH libraries resulted in 136 unigenes. Functional annotation and gene ontology analysis of up-regulated EST libraries showed several known biosynthetic genes and novel transcripts that may be involved in signaling, cellular transport, or metabolism. Real time RT-PCR analysis of a set of 8 candidate genes further confirmed the differential gene expression.

Analysis of expressed sequence tags (ESTs) from avocado seed (Persea americana var. drymifolia) reveals abundant expression of the gene encoding the antimicrobial peptide snakin

Avocado is one of the most important fruits in the world. Avocado "native mexicano" (Persea americana var. drymifolia) seeds are widely used in the propagation of this plant and are the primary source of rootstocks globally for a variety of avocado cultivars, such as the Hass avocado. Here, we report the isolation of 5005 ESTs from the 5' ends of P. americana var. drymifolia seed cDNA clones representing 1584 possible unigenes. These avocado seed ESTs were compared with the avocado flower EST library, and we detected several genes that are expressed either in both tissues or only in the seed. The snakin gene, which encodes an element of the innate immune response in plants, was one of those most frequently found among the seed ESTs, and this suggests that it is abundantly expressed in the avocado seed. We expressed the snakin gene in a heterologous system, namely the bovine endothelial cell line BVE-E6E7. Conditioned media from transfected BVE-E6E7 cells showed antimicrobial activity against strains of Escherichia coli and Staphylococcus aureus. This is the first study of the function of the snakin gene in plant seed tissue, and our observations suggest that this gene might play a protective role in the avocado seed.

Development, inheritance and breeding potential of a recessive genic male sterile line D248A in Sesame (Sesamum indicum L.)

Genic male sterility (GMS) has great potential for heterosis exploitation in sesame (Sesamum indicum L.). Two spontaneous male-sterile plants were discovered in a Chinese sesame cultivar (Zhuzhi 4) in 2006. By consecutive sib mating with fertile plants from Zhuzhi 4, a new sterile line, D248A, was developed. Anatomy studies showed that D248A has thin, small and greenish anthers on which there are no or little pollen grains. The pollens are irregularly shaped and completely aborted, resulting in no germination and no formation of pollen tubes as revealed by acetocarmine stain or semi-solid suspension culture. Furthermore, D248A has a better performance in growth vigor, bloom duration and yield per plant than the other GMS lines (i.e. 95 ms-2A and 95 ms-5A). To investigate the inheritance mode of fertility, D248A was crossed and backcrossed with six varieties, and a segregating ratio of 3:1 and 1:1 for fertile and sterile plants was observed in F₂ and BC₁ populations, respectively. These results suggested that D248A is controlled by a recessive GMS gene. The average yield of four D248A-derived F₁ hybrids is as high as 1695 kg·ha^-1, which is almost twice of that of 95 ms-5A-derived F₁ hybrids. These results indicated that this newly developed recessive GMS line has great potential in sesame hybrid breeding.

Identification and testing of reference genes for Sesame gene expression analysis by quantitative real-time PCR

Sesame (Sesamum indicum L.) is an ancient and important oilseed crop. However, few sesame reference genes have been selected for quantitative real-time PCR until now. Screening and validating reference genes is a requisite for gene expression normalization in sesame functional genomics research. In this study, ten candidate reference genes, i.e., SiACT, SiUBQ6, SiTUB, Si18S rRNA, SiEF1α, SiCYP, SiHistone, SiDNAJ, SiAPT and SiGAPDH, were chosen and examined systematically in 32 sesame samples. Three qRT-PCR analysis methods, i.e., geNorm, NormFinder and BestKeeper, were evaluated systematically. Results indicated that all ten candidate reference genes could be used as reference genes in sesame. SiUBQ6 and SiAPT were the optimal reference genes for sesame plant development; SiTUB was suitable for sesame vegetative tissue development, SiDNAJ for pathogen treatment, SiHistone for abiotic stress, SiUBQ6 for bud development and SiACT for seed germination. As for hormone treatment and seed development, SiHistone, SiCYP, SiDNAJ or SiUBQ6, as well as SiACT, SiDNAJ, SiTUB or SiAPT, could be used as reference gene, respectively. To illustrate the suitability of these reference genes, we analyzed the expression variation of three functional sesame genes of SiSS, SiLEA and SiGH in different organs using the optimal qRT-PCR system for the first time. The stability levels of optimal and worst reference genes screened for seed development, anther sterility and plant development were validated in the qRT-PCR normalization. Our results provided a reference gene application guideline for sesame gene expression characterization using qRT-PCR system.

Genome sequencing of the important oilseed crop Sesamum indicum L

The Sesame Genome Working Group (SGWG) has been formed to sequence and assemble the sesame (Sesamum indicum L.) genome. The status of this project and our planned analyses are described.

Phylogenetic, Molecular, and Biochemical Characterization of Caffeic Acid o-Methyltransferase Gene Family in Brachypodium distachyon

Caffeic acid o-methyltransferase (COMT) is one of the important enzymes controlling lignin monomer production in plant cell wall synthesis. Analysis of the genome sequence of the new grass model Brachypodium distachyon identified four COMT gene homologs, designated as BdCOMT1, BdCOMT2, BdCOMT3, and BdCOMT4. Phylogenetic analysis suggested that they belong to the COMT gene family, whereas syntenic analysis through comparisons with rice and sorghum revealed that BdCOMT4 on Chromosome 3 is the orthologous copy of the COMT genes well characterized in other grass species. The other three COMT genes are unique to Brachypodium since orthologous copies are not found in the collinear regions of rice and sorghum genomes. Expression studies indicated that all four Brachypodium COMT genes are transcribed but with distinct patterns of tissue specificity. Full-length cDNAs were cloned in frame into the pQE-T7 expression vector for the purification of recombinant Brachypodium COMT proteins. Biochemical characterization of enzyme activity and substrate specificity showed that BdCOMT4 has significant effect on a broad range of substrates with the highest preference for caffeic acid. The other three COMTs had low or no effect on these substrates, suggesting that a diversified evolution occurred on these duplicate genes that not only impacted their pattern of expression, but also altered their biochemical properties.

Development and validation of genic-SSR markers in sesame by RNA-seq

Background

Sesame (Sesamum indicum L.) is one of the most important oil crops; however, a lack of useful molecular markers hinders current genetic research. We performed transcriptome sequencing of samples from different sesame growth and developmental stages, and mining of genic-SSR markers to identify valuable markers for sesame molecular genetics research.

Results

In this study, 75 bp and 100 bp paired-end RNA-seq was used to sequence 24 cDNA libraries, and 42,566 uni-transcripts were assembled from more than 260 million filtered reads. The total length of uni-transcript sequences was 47.99 Mb, and 7,324 SSRs (SSRs ≥15 bp) and 4,440 SSRs (SSRs ≥18 bp) were identified. On average, there was one genic-SSR per 6.55 kb (SSRs ≥15 bp) or 10.81 kb (SSRs ≥18 bp). Among perfect SSRs (≥18 bp), di-nucleotide motifs (48.01%) were the most abundant, followed by tri- (20.96%), hexa- (25.37%), penta- (2.97%), tetra- (2.12%), and mono-nucleotides (0.57%). The top four motif repeats were (AG/CT)n [1,268 (34.51%)], (CA/TG)n [281 (7.65%)], (AT/AT)n [215 (5.85%)], and (GAA/TTC)n [131 (3.57%)]. A total of 2,164 SSR primer pairs were identified in the 4,440 SSR-containing sequences (≥18 bp), and 300 SSR primer pairs were randomly chosen for validation. These SSR markers were amplified and validated in 25 sesame accessions (24 cultivated accessions, one wild species). 276 (92.0%) primer pairs yielded PCR amplification products in 24 cultivars. Thirty two primer pairs (11.59%) exhibited polymorphisms. Moreover, 203 primer pairs (67.67%) yielded PCR amplicons in the wild accession and 167 (60.51%) were polymorphic between species. A UPGMA dendrogram based on genetic similarity coefficients showed that the correlation between genotype and geographical source was low and that the genetic basis of sesame in China is narrow, as previously reported. The 32 polymorphic primer pairs were validated using an F₂ mapping population; 18 primer pairs exhibited polymorphisms between the parents, and 14 genic-SSRs could be integrated into 9 main linkage groups.

Conclusions

2,164 genic-SSR markers have been developed in sesame using transcriptome sequencing. 276 of 300 validated primer pairs successfully yielded PCR amplicons in 24 cultivated sesame accessions. These markers increase current SSR marker resources and will greatly benefit genetic diversity, qualitative and quantitative trait mapping and marker-assisted selection studies in sesame.

Metabolic cartography: experimental quantification of metabolic fluxes from isotopic labelling studies

For the past decade, flux maps have provided researchers with an in-depth perspective on plant metabolism. As a rapidly developing field, significant headway has been made recently in computation, experimentation, and overall understanding of metabolic flux analysis. These advances are particularly applicable to the study of plant metabolism. New dynamic computational methods such as non-stationary metabolic flux analysis are finding their place in the toolbox of metabolic engineering, allowing more organisms to be studied and decreasing the time necessary for experimentation, thereby opening new avenues by which to explore the vast diversity of plant metabolism. Also, improved methods of metabolite detection and measurement have been developed, enabling increasingly greater resolution of flux measurements and the analysis of a greater number of the multitude of plant metabolic pathways. Methods to deconvolute organelle-specific metabolism are employed with increasing effectiveness, elucidating the compartmental specificity inherent in plant metabolism. Advances in metabolite measurements have also enabled new types of experiments, such as the calculation of metabolic fluxes based on (13)CO(2) dynamic labelling data, and will continue to direct plant metabolic engineering. Newly calculated metabolic flux maps reveal surprising and useful information about plant metabolism, guiding future genetic engineering of crops to higher yields. Due to the significant level of complexity in plants, these methods in combination with other systems biology measurements are necessary to guide plant metabolic engineering in the future.

Analysis of expression sequence tags from a full-length-enriched cDNA library of developing sesame seeds (Sesamum indicum)

BACKGROUND

Sesame (Sesamum indicum) is one of the most important oilseed crops with high oil contents and rich nutrient value. However, genetic improvement efforts in sesame could not get benefit from molecular biology technology due to poor DNA and RNA sequence resources. In this study, we carried out a large scale of expressed sequence tags (ESTs) sequencing from developing sesame seeds and further conducted analysis on seed storage products-related genes.

RESULTS

A normalized and full-length enriched cDNA library from 5 ~ 30 days old immature seeds was constructed and randomly sequenced, leading to generation of 41,248 expressed sequence tags (ESTs) which then formed 4,713 contigs and 27,708 singletons with 44.9% uniESTs being putative full-length open reading frames. Approximately 26,091 of all these uniESTs have significant matches to the counterparts in Nr database of GenBank, and 21,628 of them were assigned to one or more Gene ontology (GO) terms. Homologous genes involved in oil biosynthesis were identified including some conservative transcription factors regulating oil biosynthesis such as LEAFY COTYLEDON1 (LEC1), PICKLE (PKL), WRINKLED1 (WRI1) and majority of them were found for the first time in sesame seeds. One hundred and 17 ESTs were identified possibly involved in biosynthesis of sesame lignans, sesamin and sesamolin. In total, 9,347 putative functional genes from developing seeds were identified, which accounts for one third of total genes in the sesame genome. Further analysis of the uniESTs identified 1,949 non-redundant simple sequence repeats (SSRs).

CONCLUSIONS

This study has provided an overview of genes expressed during sesame seed development. This collection of sesame full-length cDNAs covered a wide variety of genes in seeds, in particular, candidate genes involved in biosynthesis of sesame oils and lignans. These EST sequences enriched with full length will contribute to comparative genomic studies on sesame and other oilseed plants and serve as an abundant information platform for functional marker development and functional gene study.

Characterization of the sesame (Sesamum indicum L.) global transcriptome using Illumina paired-end sequencing and development of EST-SSR markers

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.

Characterization of the sesame (Sesamum indicum L.) global transcriptome using Illumina paired-end sequencing and development of EST-SSR markers

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.

Seed-expressed casein kinase I acts as a positive regulator of the SeFAD2 promoter via phosphorylation of the SebHLH transcription factor

Microsomal oleic acid desaturase (FAD2) catalyzes the first committed step of the biosynthesis of polyunsaturated fatty acids via extra-plastidial desaturation of oleic acid to linoleic acid. In the regulatory mechanism controlling seed-specific SeFAD2 expression, trans-activation of the seed-specific SeFAD2 promoter is mediated by the SebHLH transcription factor (Kim et al. in Plant Mol Biol 64:453-466, 2007). In this study, a protein interacting with SebHLH was isolated from yeast two-hybrid analysis. The protein shares approximately 80% sequence identity with other putative casein kinases and was named SeCKI (Sesame Casein Kinase I). SeCKI transcripts were predominantly expressed in developing sesame seeds and were induced approximately threefold by exogenous application of ABA. eGFP:SeCKI fusion protein was localized to the nucleus. The SeCKI protein specifically bound to SebHLH. The SeCKI protein was autophosphorylated in a calcium-independent manner and transphosphorylated the SebHLH protein. Both the SebHLH and the SeCKI genes or both the SebHLH and mutated SemCKI (K182G) genes, under the control of CaMV 35S promoter, and the GUS reporter gene driven by SeFAD2 promoter containing E- and G-Box motifs were co-expressed in developing sesame seeds. This co-expression revealed that SeCKI enhanced the SebHLH-mediated transactivation of the SeFAD2 gene promoter via phosphorylation of the SebHLH transcription factor.

Global analysis of gene expression profiles in Brassica napus developing seeds reveals a conserved lipid metabolism regulation with Arabidopsis thaliana

In order to study Brassica napus fatty acid (FA) metabolism and relevant regulatory networks, a systematic identification of fatty acid (FA) biosynthesis-related genes was conducted. Following gene identification, gene expression profiles during B. napus seed development and FA metabolism were performed by cDNA chip hybridization (>8000 EST clones from seed). The results showed that FA biosynthesis and regulation, and carbon flux, were conserved between B. napus and Arabidopsis. However, a more critical role of starch metabolism was detected for B. napus seed FA metabolism and storage-component accumulation when compared with Arabidopsis. In addition, a crucial stage for the transition of seed-to-sink tissue was 17-21 d after flowering (DAF), whereas FA biosynthesis-related genes were highly expressed primarily at 21 DAF. Hormone (auxin and jasmonate) signaling is found to be important for FA metabolism. This study helps to reveal the global regulatory network of FA metabolism in developing B. napus seeds.

Isolation and characterization of multiple abundant lipid transfer protein isoforms in developing sesame (Sesamum indicum L.) seeds

Sesame (Sesamum indicum) is an important oilseed crop; approximately 50% of the seed dry weight is storage oil. In a previous report, developing sesame seed expressed sequence tags (ESTs) revealed that ESTs encoding lipid transfer protein (LTPs) were one of the most abundant groups of sesame ESTs. LTP functions in the transfer of wax or cutin monomers and in the defense response against pathogen attack. To study the biological role of the abundant LTP isoforms in developing seeds, 122 ESTs out of 3328 sesame ESTs were analyzed against Arabidopsis and rice proteome databases. LTP fraction, which was partially purified from developing sesame seeds, actively transferred fluorescent phospholipids and bound to fatty acids. Full-length cDNAs of five out of 21 LTP isoforms were isolated and named SiLTP1-SiLTP5. The predicted amino acid sequences of the five SiLTPs harbor typical characteristics of LTPs, including conserved arrangement of cysteine residues. Northern blot analysis revealed that the five SiLTP isoforms were most abundantly expressed in developing seeds, but were also detected in flower tissues. Also, SiLTP3 and SiLTP4 transcripts were expressed in leaves and seed-pot walls, respectively. In addition, SiLTP2 and SiLTP4 transcripts were significantly induced in 6-day-old sesame seedlings by application of NaCl, mannitol, and abscisic acid (ABA). Transient expression of green fluorescent protein (GFP)-fusion constructs in Arabidopsis protoplasts revealed that SiLTP1 and SiLTP2 were secreted by different pathways. Taken together, the abundant LTPs in developing sesame seeds are involved in lipid transfer into the extracellular matrix. Possible biological roles of SiLTPs related to organ-specific expression and abiotic stresses are discussed.

The SebHLH transcription factor mediates trans-activation of the SeFAD2 gene promoter through binding to E- and G-box elements

Microsomal oleic acid desaturase (FAD2) catalyzes the first extra-plastidial desaturation in plants, converting oleic acid to linoleic acid, which is a major constituent in all cellular membranes as well as in seed storage oils. Seed-specific FAD2 (SeFAD2) produced 40% of linoleic acids in the total fatty acids of sesame (Sesamum indicum) seeds. The expression of SeFAD2 transcripts was spatially and temporally controlled during seed development. To investigate the regulatory mechanism controlling seed-specific SeFAD2 expression, we isolated a well-matched sequence homologous to the basic region/helix-loop-helix proteins by yeast one-hybrid screening and named it SebHLH. SebHLH transcripts were expressed in developing seeds and roots of sesame. SebHLH:GFP fusion protein localized in the nucleus. Recombinant SebHLH protein bound E-box (CANNTG) and G-box (CACGTG) elements in the region from -179 to -53 of the SeFAD2 gene promoter, and the external C and G nucleotides in the E- and G-box motifs were essential for SebHLH protein binding. The SebHLH gene, under the CaMV35S promoter, and the GUS reporter gene driven by E- and G-box motifs were co-expressed in developing sesame seeds and Arabidopsis transgenic leaves. This co-expression demonstrated that SebHLH protein mediates transactivation of the SeFAD2 gene promoter through binding to E- and G-box elements. E- or G-box elements frequently occur in the 5'-flanking region of genes that are involved in triacylglycerol biosynthesis and that exhibit seed-specific expression in Arabidopsis and other plants, suggesting that bHLH transcription factors play a key role in the transcriptional regulation of genes related to storage lipid biosynthesis and accumulation during seed development.

Evidence for a role of AtCAD 1 in lignification of elongating stems of Arabidopsis thaliana

The cinnamyl alcohol dehydrogenase (AtCAD) multigene family in Arabidopsis is composed of nine genes. Our previous studies focused on the two isoforms AtCAD C and AtCAD D which show a high homology to those related to lignification in other plants. This study focuses on the seven other Arabidopsis CAD for which functions are not yet elucidated. Their expression patterns were determined in different parts of Arabidopsis. Only CAD 1 protein can be detected in elongating stems, flowers, and siliques using Western-blot analysis. Tissue specific expression of CAD 1, B1, and G genes was determined using their promoters fused to the GUS reporter gene. CAD 1 expression was observed in primary xylem in accordance with a potential role in lignification. Arabidopsis T-DNA mutants knockout for the different genes CAD genes were characterized. Their stems displayed no substantial reduction of CAD activities for coniferyl and sinapyl alcohols as well as no modifications of lignin quantity and structure in mature inflorescence stems. Only a small reduction of lignin content could be observed in elongating stems of Atcad 1 mutant. These CAD genes in combination with the CAD D promoter were used to complement a CAD double mutant severely altered in lignification (cad c cad d). The expression of AtCAD A, B1, B2, F, and G had no effect on restoring a normal lignin profile of this mutant. In contrast, CAD 1 complemented partly this mutant as revealed by the partial restoration of conventional lignin units and by the decrease in the frequency of beta-O-4 linked p-OH cinnamaldehydes.

Formation of two methylenedioxy bridges by a Sesamum CYP81Q protein yielding a furofuran lignan, (+)-sesamin

(+)-Sesamin, a furofuran class lignan, is widespread in vascular plants and represented by Sesamum spp. (+)-Sesamin has been of rapidly growing interest because of its beneficial biological effects in mammals, but its biosynthesis and physiological roles in plants remain to be clarified. It is speculated to be synthesized from (+)-pinoresinol by means of (+)-piperitol by formation of two methylenedioxy bridges mediated by two distinct Sesamum indicum cytochrome P450 (SiP450) proteins. Here, we report an SiP450, CYP81Q1, that alone catalyzes (+)-sesamin biosynthesis from (+)-pinoresinol by means of (+)-piperitol by forming two methylenedioxy bridges. The CYP81Q1 gene expression profile was temporally consistent with the accumulation pattern of (+)-sesamin during seed development. The CYP81Q1-GFP chimera protein was colocalized with an endoplasmic reticulum (ER)-targeting chimera protein, indicating that (+)-sesamin biosynthesis occurs on the ER cytoplasmic surface. Moreover, we isolated two CYP81Q1 homologs from other Sesamum spp. Sesamum radiatum CYP81Q2 showed dual (+)-piperitol/(+)-sesamin synthetic activity. CYP81Q2, as well as CYP81Q1, therefore, corresponds to a (+)-piperitol/(+)-sesamin synthase in lignan biosynthesis. In contrast, Sesamum alatum CYP81Q3 showed no activity, in accord with (+)-sesamin being deficient in S. alatum. Our findings not only provide insight into lignan biosynthesis but also unravel a unique mode of cytochrome P450 action.

Primary metabolic pathways and signal transduction in sunflower (Helianthus annuus L.): comparison of transcriptional profiling in leaves and immature embryos using cDNA microarrays

The early stage of embryo development is a critical step in plant production. To identify genes with potential roles in the early sunflower seed development, a cDNA microarray approach was employed. We developed a thematic cDNA microarray containing clones representing high sequence similarities with known or predicted Arabidopsis genes implicated in different metabolic and signal transduction pathways. This 800-element cDNA array was used to compare the expression patterns in leaves and immature embryos (2 mm and 6 mm). Statistical analysis, using two-step ANOVA, revealed that 143 cDNA clones can be considered as differentially expressed. Of these, 62 clones were found to be up-regulated in leaves, 81 in embryos whereas only seven clones displayed increased level of mRNA in the 6 mm embryos when compared with 2 mm embryos. The differentially expressed clones are distributed among many metabolic and signal transduction pathways. For example, genes related to fatty acid metabolism and amino acid biosynthesis exhibited preferential expression patterns in immature embryos. Also, clones potentially encoding enzymes involved in the metabolism of ascorbate and aldarate, pyruvate, propanoate and inositol, and citrate cycle were found to be up-regulated in embryos. In contrast, cDNA clones putatively involved in energy metabolism were more abundant in leaves than embryos. Clones encoding potential signal transduction components including receptors, protein kinases, protein phosphatases, and transcription factors were also identified, with preferential expression profiles in immature embryos. The expression patterns derived from this study provide initial characterization of metabolic pathways and signalling transduction networks occurring in the early stage of sunflower seed development.

Identification of new gene expression regulators specifically expressed during plant seed maturation

A cDNA-AFLP approach on Linum usitatissimum (flax) was used to identify genes specifically expressed during the seed maturation process. Among the 20,000 cDNA-AFLP tags produced, 486 were selected for their seed-specific expression during maturation. When compared with the publicly available databases, half of them presented some significant similarity with known plant sequences. The results obtained confirmed the accuracy of the approach as numerous genes previously described as being expressed exclusively in plant seeds were identified in this screen. The focus was on sequences similar to plant regulators involved in the control of gene expression, either at the transcriptional, post-transcriptional, or post-translational levels. Using a real-time RT-PCR approach, seed-specific expression kinetics were confirmed for 13 of these regulators that were never characterized for being expressed during seed maturation. Among these, a flax gene of the non-LEC1-like HAP3 family and a flax MYB factor were shown to be expressed in specialized tissues of flax embryo using an in situ hybridization approach. By expression kinetic comparison between these flax genes and their Arabidopsis counterparts, it was found that the new HAP3 gene should be related to a ubiquitous seed maturation mechanism, while a new MYB factor appears to be related to a more seed-specific maturation mechanism. These results demonstrate the utility of the flax database in not only identifying new genes expressed during seed maturation but also in being able to highlight the distinction between conserved and non-conserved seed maturation mechanisms.

Coffee and tomato share common gene repertoires as revealed by deep sequencing of seed and cherry transcripts

An EST database has been generated for coffee based on sequences from approximately 47,000 cDNA clones derived from five different stages/tissues, with a special focus on developing seeds. When computationally assembled, these sequences correspond to 13,175 unigenes, which were analyzed with respect to functional annotation, expression profile and evolution. Compared with Arabidopsis, the coffee unigenes encode a higher proportion of proteins related to protein modification/turnover and metabolism-an observation that may explain the high diversity of metabolites found in coffee and related species. Several gene families were found to be either expanded or unique to coffee when compared with Arabidopsis. A high proportion of these families encode proteins assigned to functions related to disease resistance. Such families may have expanded and evolved rapidly under the intense pathogen pressure experienced by a tropical, perennial species like coffee. Finally, the coffee gene repertoire was compared with that of Arabidopsis and Solanaceous species (e.g. tomato). Unlike Arabidopsis, tomato has a nearly perfect gene-for-gene match with coffee. These results are consistent with the facts that coffee and tomato have a similar genome size, chromosome karyotype (tomato, n=12; coffee n=11) and chromosome architecture. Moreover, both belong to the Asterid I clade of dicot plant families. Thus, the biology of coffee (family Rubiacaeae) and tomato (family Solanaceae) may be united into one common network of shared discoveries, resources and information.

Comparative analysis of early embryonic sunflower cDNA libraries

To gain information concerning cell functions and activities during sunflower embryogenesis, an expressed sequence tag (EST) approach was used to analyse gene expression in the early stages of sunflower embryos development. Confocal microscopy observations of whole-mounted embryos allowed us to identify precisely the major steps of the zygotic embryonic development. A time-course analysis was then employed to collect the embryonic material. Three cDNA libraries were constructed from microdissected embryos, and three other cDNA libraries were created using a classical day after pollination schedule. A total of 7106 ESTs were produced and assembled. The total number of putative different genes represents about 43.1 (3064 tentative contigs and singlets) of the analysed sequences. The unigenes that showed similarity to proteins with known or predicted functions (50.3) were classified into 15 different functional categories. The functional profiles were found to be quite similar for all studied embryo stages but statistical analysis revealed that successive and coordinate sets of genes are expressed at each embryonic stage. The analysis allowed us to identify abundant and differentially expressed genes at the early stages of embryos development as well as some putatively interesting genes, showing strong similarities with genes playing key roles in plant and animal embryogenesis. The data presented in this study not only provide a first global overview of the genes expression profile during sunflower embryogenesis but also represent an original and valuable tool for developmental genomics studies on exalbuminous dicots.