De novo assembly and transcriptome analysis of five major tissues of Jatropha curcas L. using GS FLX titanium platform of 454 pyrosequencing

Genome-Wide Identification and Comparative Analysis of WOX Genes in Four Euphorbiaceae Species and Their Expression Patterns in Jatropha curcas

The WUSCHEL-related homeobox (WOX) proteins are widely distributed in plants and play important regulatory roles in growth and development processes such as embryonic development and organ development. Here, series of bioinformatics methods were utilized to unravel the structural basis and genetic hierarchy of WOX genes, followed by regulation of the WOX genes in four Euphorbiaceae species. A genome-wide survey identified 59 WOX genes in Hevea brasiliensis (H. brasiliensis: 20 genes), Jatropha curcas (J. curcas: 10 genes), Manihot esculenta (M. esculenta: 18 genes), and Ricinus communis (R. communis: 11 genes). The phylogenetic analysis revealed that these WOX members could be clustered into three close proximal clades, such as namely ancient, intermediate and modern/WUS clades. In addition, gene structures and conserved motif analyses further validated that the WOX genes were conserved within each phylogenetic clade. These results suggested the relationships among WOX members in the four Euphorbiaceae species. We found that WOX genes in H. brasiliensis and M. esculenta exhibit close genetic relationship with J. curcas and R. communis. Additionally, the presence of various cis-acting regulatory elements in the promoter of J. curcas WOX genes (JcWOXs) reflected distinct functions. These speculations were further validated with the differential expression profiles of various JcWOXs in seeds, reflecting the importance of two JcWOX genes (JcWOX6 and JcWOX13) during plant growth and development. Our quantitative real-time PCR (qRT-PCR) analysis demonstrated that the JcWOX11 gene plays an indispensable role in regulating plant callus. Taken together, the present study reports the comprehensive characteristics and relationships of WOX genes in four Euphorbiaceae species, providing new insights into their characterization.

Genome-Wide Identification, Expression Patterns and Sugar Transport of the Physic Nut SWEET Gene Family and a Functional Analysis of JcSWEET16 in Arabidopsis

The Sugars Will Eventually be Exported Transporters (SWEET) family is a class of sugar transporters that play key roles in phloem loading, seed filling, pollen development and the stress response in plants. Here, a total of 18 JcSWEET genes were identified in physic nut (Jatropha curcas L.) and classified into four clades by phylogenetic analysis. These JcSWEET genes share similar gene structures, and alternative splicing of messenger RNAs was observed for five of the JcSWEET genes. Three (JcSWEET1/4/5) of the JcSWEETs were found to possess transport activity for hexose molecules in yeast. Real-time quantitative PCR analysis of JcSWEETs in different tissues under normal growth conditions and abiotic stresses revealed that most are tissue-specifically expressed, and 12 JcSWEETs responded to either drought or salinity. The JcSWEET16 gene responded to drought and salinity stress in leaves, and the protein it encodes is localized in both the plasma membrane and the vacuolar membrane. The overexpression of JcSWEET16 in Arabidopsis thaliana modified the flowering time and saline tolerance levels but not the drought tolerance of the transgenic plants. Together, these results provide insights into the characteristics of SWEET genes in physic nut and could serve as a basis for cloning and further functional analysis of these genes.

Extended mining of the oil biosynthesis pathway in biofuel plant Jatropha curcas by combined analysis of transcriptome and gene interactome data

BACKGROUND

Jatropha curcas L. is an important non-edible oilseed crop with a promising future in biodiesel production. However, little is known about the molecular biology of oil biosynthesis in this plant when compared with other established oilseed crops, resulting in the absence of agronomically improved varieties of Jatropha. To extensively discover the potentially novel genes and pathways associated with the oil biosynthesis in J. curcas, new strategy other than homology alignment is on the demand.

RESULTS

In this study, we proposed a multi-step computational framework that integrates transcriptome and gene interactome data to predict functional pathways in non-model organisms in an extended process, and applied it to study oil biosynthesis pathway in J. curcas. Using homologous mapping against Arabidopsis and transcriptome profile analysis, we first constructed protein-protein interaction (PPI) and co-expression networks in J. curcas. Then, using the homologs of Arabidopsis oil-biosynthesis-related genes as seeds, we respectively applied two algorithm models, random walk with restart (RWR) in PPI network and negative binomial distribution (NBD) in co-expression network, to further extend oil-biosynthesis-related pathways and genes in J. curcas. At last, using k-nearest neighbors (KNN) algorithm, the predicted genes were further classified into different sub-pathways according to their possible functional roles.

CONCLUSIONS

Our method exhibited a highly efficient way of mining the extended oil biosynthesis pathway of J. curcas. Overall, 27 novel oil-biosynthesis-related gene candidates were predicted and further assigned to 5 sub-pathways. These findings can help better understanding of the oil biosynthesis pathway of J. curcas, as well as paving the way for the following J. curcas breeding application.

Genome-wide characterization of bZIP transcription factors and their expression patterns in response to drought and salinity stress in Jatropha curcas

The basic leucine zipper (bZIP) family is one of the largest families of transcription factors (TFs) in plants and is responsible for various functions, including regulating development and responses to abiotic/biotic stresses. However, the roles of bZIPs in the regulation of responses to drought stress and salinity stress remain poorly understood in Jatropha curcas L., a biodiesel crop. In the present study, 50 JcbZIP genes were identified and classified into ten groups. Cis-element analysis indicated that JcbZIP genes are associated with abiotic stress. Gene expression patterns and quantitative real-time PCR (qRT-PCR) showed that four JcbZIP genes (JcbZIPs 34, 36, 49 and 50) are key resistance-related genes under both drought and salinity stress conditions. On the basis of the results of cis-element and phylogenetic analyses, JcbZIP49 and JcbZIP50 are likely involved in responses to drought and salinity stress; moreover, JcbZIP34 and JcbZIP36 might also play important roles in seed development and response to abiotic stress. These findings advance our understanding of the comprehensive characteristics of JcbZIP genes and provide new insights for functional validation in the further.

Comparative gene expression profile analysis of ovules provides insights into Jatropha curcas L. ovule development

Jatropha curcas, an economically important biofuel feedstock with oil-rich seeds, has attracted considerable attention among researchers in recent years. Nevertheless, valuable information on the yield component of this plant, particularly regarding ovule development, remains scarce. In this study, transcriptome profiles of anther and ovule development were established to investigate the ovule development mechanism of J. curcas. In total, 64,325 unigenes with annotation were obtained, and 1723 differentially expressed genes (DEGs) were identified between different stages. The DEG analysis showed the participation of five transcription factor families (bHLH, WRKY, MYB, NAC and ERF), five hormone signaling pathways (auxin, gibberellic acid (GA), cytokinin, brassinosteroids (BR) and jasmonic acid (JA)), five MADS-box genes (AGAMOUS-2, AGAMOUS-1, AGL1, AGL11, and AGL14), SUP and SLK3 in ovule development. The role of GA and JA in ovule development was evident with increases in flower buds during ovule development: GA was increased approximately twofold, and JA was increased approximately sevenfold. In addition, the expression pattern analysis using qRT-PCR revealed that CRABS CLAW and AGAMOUS-2 were also involved in ovule development. The upregulation of BR signaling genes during ovule development might have been regulated by other phytohormone signaling pathways through crosstalk. This study provides a valuable framework for investigating the regulatory networks of ovule development in J. curcas.

Genomics analysis of genes encoding respiratory burst oxidase homologs (RBOHs) in jatropha and the comparison with castor bean

Respiratory burst oxidase homologs (RBOHs), which catalyze the production of superoxide from oxygen and NADPH, play key roles in plant growth and development, hormone signaling, and stress responses. Compared with extensive studies in model plants arabidopsis and rice, little is known about RBOHs in other species. This study presents a genome-wide analysis of Rboh family genes in jatropha (Jatropha curcas) as well as the comparison with castor bean (Ricinus communis), another economically important non-food oilseed crop of the Euphorbiaceae family. The family number of seven members identified from the jatropha genome is equal to that present in castor bean, and further phylogenetic analysis assigned these genes into seven groups named RBOHD, -C, -B, -E, -F, -N, and -H. In contrast to a high number of paralogs present in arabidopsis and rice that experienced several rounds of recent whole-genome duplications, no duplicate was identified in both jatropha and castor bean. Conserved synteny and one-to-one orthologous relationship were observed between jatropha and castor bean Rboh genes. Although exon-intron structures are usually highly conserved between orthologs, loss of certain introns was observed for JcRbohB, JcRbohD, and RcRbohN, supporting their divergence. Global gene expression profiling revealed diverse patterns of JcRbohs over various tissues. Moreover, expression patterns of JcRbohs during flower development as well as various stresses were also investigated. These findings will not only improve our knowledge on species-specific evolution of the Rboh gene family, but also provide valuable information for further functional analysis of Rboh genes in jatropha.

Transcriptome Analysis Reveals Dynamic Fat Accumulation in the Walnut Kernel

Walnut (Juglans regia L.) is an important woody oilseed species cultivated throughout the world. In this study, comparative transcript profiling was performed using high-throughput RNA sequencing technology at the following three stages of walnut fat synthesis in the “Lvling” walnut cultivar: the initial developmental stage (L1), the fast developing stage (L2), and the last developing stage (L3). A total of 68.18 GB of data were obtained on the three developmental stages, and 92% to 94% of clean data were able to be located to the reference genome. Further comparisons of the transcripts in the three libraries revealed that 724, 2027, and 4817 genes were differentially expressed between the L2 and L1 (L2vsL1), L3 and L2 (L3vsL2), and L3 and L1 (L3vsL1) samples, respectively. Through the GO gene enrichment analysis, differentially expressed genes (DEGs) in L2vsL1, L3vsL2, and L3vsL1 were enriched into 3, 0, and 2 functional categories, respectively. According to the KEGG enrichment analysis, DEGs in L2vsL1, L3vsL2, and L3vsL1 were annotated into 77, 110, and 3717 taxonomic metabolic pathways in the KEGG database, respectively. Next, we analyzed expression levels of genes related to fat synthesis. Our results indicated that ACCase, LACS, and FAD7 were the key genes related to fat synthesis. The high-throughput transcriptome sequencing of walnut in different developmental stages has greatly enriched the current genomic available resources. The comparison of DEGs under different developmental stages identified a wealth of candidate genes involved in fat synthesis, which will facilitate further genetic improvement and molecular studies of the walnut.

Genome-wide comparative analysis of papain-like cysteine protease family genes in castor bean and physic nut

Papain-like cysteine proteases (PLCPs) are a class of proteolytic enzymes involved in many plant processes. Compared with the extensive research in Arabidopsis thaliana, little is known in castor bean (Ricinus communis) and physic nut (Jatropha curcas), two Euphorbiaceous plants without any recent whole-genome duplication. In this study, a total of 26 or 23 PLCP genes were identified from the genomes of castor bean and physic nut respectively, which can be divided into nine subfamilies based on the phylogenetic analysis: RD21, CEP, XCP, XBCP3, THI, SAG12, RD19, ALP and CTB. Although most of them harbor orthologs in Arabidopsis, several members in subfamilies RD21, CEP, XBCP3 and SAG12 form new groups or subgroups as observed in other species, suggesting specific gene loss occurred in Arabidopsis. Recent gene duplicates were also identified in these two species, but they are limited to the SAG12 subfamily and were all derived from local duplication. Expression profiling revealed diverse patterns of different family members over various tissues. Furthermore, the evolution characteristics of PLCP genes were also compared and discussed. Our findings provide a useful reference to characterize PLCP genes and investigate the family evolution in Euphorbiaceae and species beyond.

Cloning and gene expression analysis of two cDNA of cysteine proteinase genes involved in programmed cell death in the inner integument from developing seeds of Jatropha curcas L

In this paper, two cysteine proteinases were cloned from Jatropha curcas seeds. The full length cDNAs obtained from cloning of Jc-CysEP1 and Jc-CysEP2 genes were 1.516bp and 1500 pb, respectively. The Jc-CysEP1 contained a 1083bp open reading frame (ORF) coding for 360 amino acids. The JcCysEP1 protein sequence had an estimated native molecular weight of 36.89 kDa, with a predicted isoelectric point of 4.55. The average lengths of JcCysEP1 5' UTR and 3' UTR were 269 bp and 167bp, respectively. The Jc-CysEP2 contained a 1077 pb open reading frame (ORF) that encoded 358 amino acids. We also identified UTRs with lengths of 229 pb (5'UTR) and 194 pb (3'UTR). The Jc-CysEP2 sequence had a native molecular weight of 39.94 kDa, with a predicted isoelectric point of 6.19. Real-time PCR analyses of developing seeds (stages I-VII) showed that most cysteine proteinase genes were expressed at stage IV (middle stage) revealing peculiar spatio-temporal differences. JcCysEP2 was the cysteine proteinase gene with the highest expression in inner integument tissue, while JcCysEP1 was expressed in lower levels. Our results suggest that JcCysEP2 could be the major cysteine proteinase gene involved in PCD events in inner integument tissue, playing a critical role in PCD events during seed development, while Jc-CyEP1 and JcCysEP2 genes act cooperatively in stages IV-VII. JcCysEP2 is important to complete their participation in PCD until development of seeds.

Seed development of Jatropha curcas L. (Euphorbiaceae): integrating anatomical, ultrastructural and molecular studies

This work provides a detailed histological analysis of the development of Jatropha curcas seeds, together with an assessment of the role of programmed cell death in this process. Seeds of Jatropha curcas are a potential source of raw material for the production of biodiesel, but very little is known about how the architecture of the seeds is shaped by the coordinated development of the embryo, endosperm and maternal tissues, namely integuments and nucellus. This study used standard anatomical and ultrastructural techniques to evaluate seed development and programmed cell death (PCD) in the inner integument was monitored by qPCR. In these studies, we found that the embryo sac formation is of the Polygonum type. We also found that embryogenesis is a slow process and the embryo is nourished by the suspensor at earlier stages and by nutrients remobilized from the lysis of the inner integument at later stages. Two types of programmed cell death contribute to the differentiation of the inner integument that begins at early stages of seed development. In addition, the mature embryo presents features of adaptation to dry environments such as the presence of four seminal roots, water absorbing stomata in the root zone and already differentiated protoxylem elements. The findings in this study fill in gaps related to the ontogeny of J. curcas seed development and provide novel insights regarding the types of PCD occurring in the inner integument.

Enriching Genomic Resources and Marker Development from Transcript Sequences of Jatropha curcas for Microgravity Studies

Jatropha (Jatropha curcas L.) is an economically important species with a great potential for biodiesel production. To enrich the jatropha genomic databases and resources for microgravity studies, we sequenced and annotated the transcriptome of jatropha and developed SSR and SNP markers from the transcriptome sequences. In total 1,714,433 raw reads with an average length of 441.2 nucleotides were generated. De novo assembling and clustering resulted in 115,611 uniquely assembled sequences (UASs) including 21,418 full-length cDNAs and 23,264 new jatropha transcript sequences. The whole set of UASs were fully annotated, out of which 59,903 (51.81%) were assigned with gene ontology (GO) term, 12,584 (10.88%) had orthologs in Eukaryotic Orthologous Groups (KOG), and 8,822 (7.63%) were mapped to 317 pathways in six different categories in Kyoto Encyclopedia of Genes and Genome (KEGG) database, and it contained 3,588 putative transcription factors. From the UASs, 9,798 SSRs were discovered with AG/CT as the most frequent (45.8%) SSR motif type. Further 38,693 SNPs were detected and 7,584 remained after filtering. This UAS set has enriched the current jatropha genomic databases and provided a large number of genetic markers, which can facilitate jatropha genetic improvement and many other genetic and biological studies.

Integrated analysis of 454 and Illumina transcriptomic sequencing characterizes carbon flux and energy source for fatty acid synthesis in developing Lindera glauca fruits for woody biodiesel

BACKGROUND

Lindera glauca fruit with high quality and quantity of oil has emerged as a novel potential source of biodiesel in China, but the molecular regulatory mechanism of carbon flux and energy source for oil biosynthesis in developing fruits is still unknown. To better develop fruit oils of L. glauca as woody biodiesel, a combination of two different sequencing platforms (454 and Illumina) and qRT-PCR analysis was used to define a minimal reference transcriptome of developing L. glauca fruits, and to construct carbon and energy metabolic model for regulation of carbon partitioning and energy supply for FA biosynthesis and oil accumulation.

RESULTS

We first analyzed the dynamic patterns of growth tendency, oil content, FA compositions, biodiesel properties, and the contents of ATP and pyridine nucleotide of L. glauca fruits from seven different developing stages. Comprehensive characterization of transcriptome of the developing L. glauca fruit was performed using a combination of two different next-generation sequencing platforms, of which three representative fruit samples (50, 125, and 150 DAF) and one mixed sample from seven developing stages were selected for Illumina and 454 sequencing, respectively. The unigenes separately obtained from long and short reads (201, and 259, respectively, in total) were reconciled using TGICL software, resulting in a total of 60,031 unigenes (mean length = 1061.95 bp) to describe a transcriptome for developing L. glauca fruits. Notably, 198 genes were annotated for photosynthesis, sucrose cleavage, carbon allocation, metabolite transport, acetyl-CoA formation, oil synthesis, and energy metabolism, among which some specific transporters, transcription factors, and enzymes were identified to be implicated in carbon partitioning and energy source for oil synthesis by an integrated analysis of transcriptomic sequencing and qRT-PCR. Importantly, the carbon and energy metabolic model was well established for oil biosynthesis of developing L. glauca fruits, which could help to reveal the molecular regulatory mechanism of the increased oil production in developing fruits.

CONCLUSIONS

This study presents for the first time the application of an integrated two different sequencing analyses (Illumina and 454) and qRT-PCR detection to define a minimal reference transcriptome for developing L. glauca fruits, and to elucidate the molecular regulatory mechanism of carbon flux control and energy provision for oil synthesis. Our results will provide a valuable resource for future fundamental and applied research on the woody biodiesel plants.

Transcriptome analysis revealed the dynamic oil accumulation in Symplocos paniculata fruit

BACKGROUND

Symplocos paniculata, asiatic sweetleaf or sapphire berry, is a widespread shrub or small tree from Symplocaceae with high oil content and excellent fatty acid composition in fruit. It has been used as feedstocks for biodiesel and cooking oil production in China. Little transcriptome information is available on the regulatory molecular mechanism of oil accumulation at different fruit development stages.

RESULTS

The transcriptome at four different stages of fruit development (10, 80,140, and 170 days after flowering) of S. paniculata were analyzed. Approximately 28 million high quality clean reads were generated. These reads were trimmed and assembled into 182,904 non-redundant putative transcripts with a mean length of 592.91 bp and N50 length of 785 bp, respectively. Based on the functional annotation through Basic Local Alignment Search Tool (BLAST) with public protein database, the key enzymes involved in lipid metabolism were identified, and a schematic diagram of the pathway and temporal expression patterns of lipid metabolism was established. About 13,939 differentially expressed unigenes (DEGs) were screened out using differentially expressed sequencing (DESeq) method. The transcriptional regulatory patterns of the identified enzymes were highly related to the dynamic oil accumulation along with the fruit development of S. paniculata. In addition, quantitative real-time PCR (qRT-PCR) of six vital genes was significantly correlated with DESeq data.

CONCLUSIONS

The transcriptome sequences obtained and deposited in NCBI would enrich the public database and provide an unprecedented resource for the discovery of the genes associated with lipid metabolism pathway in S. paniculata. Results in this study will lay the foundation for exploring transcriptional regulatory profiles, elucidating molecular regulatory mechanisms, and accelerating genetic engineering process to improve the yield and quality of seed oil of S. paniculata.

Unravelling molecular mechanisms from floral initiation to lipid biosynthesis in a promising biofuel tree species, Pongamia pinnata using transcriptome analysis

Pongamia pinnata (L.) (Fabaceae) is a promising biofuel tree species which is underexploited in the areas of both fundamental and applied research, due to the lack of information either on transcriptome or genomic data. To investigate the possible metabolic pathways, we performed whole transcriptome analysis of Pongamia through Illumina NextSeq platform and generated 2.8 GB of paired end sequence reads. The de novo assembly of raw reads generated 40,000 contigs and 35,000 transcripts, representing leaf, flower and seed unigenes. Spatial and temporal expression profiles of photoperiod and floral homeotic genes in Pongamia, identified GIGANTEA (GI) - CONSTANS (CO) - FLOWERING LOCUS T (FT) as active signal cascade for floral initiation. Four prominent stages of seed development were selected in a high yielding Pongamia accession (TOIL 1) to follow the temporal expression patterns of key fatty acid biosynthetic genes involved in lipid biosynthesis and accumulation. Our results provide insights into an array of molecular events from flowering to seed maturity in Pongamia which will provide substantial basis for modulation of fatty acid composition and enhancing oil yields which should serve as a potential feedstock for biofuel production.

De novo transcriptome sequencing facilitates genomic resource generation in Tinospora cordifolia

Tinospora cordifolia is known for its medicinal properties owing to the presence of useful constituents such as terpenes, glycosides, steroids, alkaloids, and flavonoids belonging to secondary metabolism origin. However, there is little information available pertaining to critical genomic elements (ESTs, molecular markers) necessary for judicious exploitation of its germplasm. We employed 454 GS-FLX pyrosequencing of entire transcripts and altogether ∼25 K assembled transcripts or Expressed sequence tags (ESTs) were identified. As the interest in T. cordifolia is primarily due to its secondary metabolite constituents, the ESTs pertaining to terpenoids biosynthetic pathway were identified in the present study. Additionally, several ESTs were assigned to different transcription factor families. To validate our transcripts dataset, the novel EST-SSR markers were generated to assess the genetic diversity among germplasm of T. cordifolia. These EST-SSR markers were found to be polymorphic and the dendrogram based on dice similarity index revealed three distinct clustering of accessions. The present study demonstrates effectiveness in using both NEWBLER and MIRA sequence read assembler software for enriching transcript-dataset and thus enables better exploitation of EST resources for mining candidate genes and designing molecular markers.

Hypobaric Treatment Effects on Chilling Injury, Mitochondrial Dysfunction, and the Ascorbate-Glutathione (AsA-GSH) Cycle in Postharvest Peach Fruit

In this study, hypobaric treatment effects were investigated on chilling injury, mitochondrial dysfunction, and the ascorbate-glutathione (AsA-GSH) cycle in peach fruit stored at 0 °C. Internal browning of peaches was dramatically reduced by applying 10-20 kPa pressure. Hypobaric treatment markedly inhibited membrane fluidity increase, whereas it kept mitochondrial permeability transition pore (MPTP) concentration and cytochrome C oxidase (CCO) and succinic dehydrogenase (SDH) activity relatively high in mitochondria. Similarly, 10-20 kPa pressure treatment reduced the level of decrease observed in AsA and GSH concentrations, while it enhanced ascorbate peroxidase (APX), glutathione reductase (GR), and monodehydroascorbate reductase (MDHAR) activities related to the AsA-GSH cycle. Furthermore, comparative transcriptomic analysis showed that differentially expressed genes (DEGs) associated with the metabolism of glutathione, ascorbate, and aldarate were up-regulated in peaches treated with 10-20 kPa for 30 days at 0 °C. Genes encoding GR, MDHAR, and APX were identified and exhibited higher expression in fruits treated with low pressure than in fruits treated with normal atmospheric pressure. Our findings indicate that the alleviation of chilling injury by hypobaric treatment was associated with preventing mitochondrial dysfunction and triggering the AsA-GSH cycle by the transcriptional up-regulation of related enzymes.

Genome-Wide Identification of Jatropha curcas Aquaporin Genes and the Comparative Analysis Provides Insights into the Gene Family Expansion and Evolution in Hevea brasiliensis

Aquaporins (AQPs) are channel-forming integral membrane proteins that transport water and other small solutes across biological membranes. Despite the vital role of AQPs, to date, little is known in physic nut (Jatropha curcas L., Euphorbiaceae), an important non-edible oilseed crop with great potential for the production of biodiesel. In this study, 32 AQP genes were identified from the physic nut genome and the family number is relatively small in comparison to 51 in another Euphorbiaceae plant, rubber tree (Hevea brasiliensis Muell. Arg.). Based on the phylogenetic analysis, the JcAQPs were assigned to five subfamilies, i.e., nine plasma membrane intrinsic proteins (PIPs), nine tonoplast intrinsic proteins (TIPs), eight NOD26-like intrinsic proteins (NIPs), two X intrinsic proteins (XIPs), and four small basic intrinsic proteins (SIPs). Like rubber tree and other plant species, functional prediction based on the aromatic/arginine selectivity filter, Froger's positions, and specificity-determining positions showed a remarkable difference in substrate specificity among subfamilies of JcAQPs. Genome-wide comparative analysis revealed the specific expansion of PIP and TIP subfamilies in rubber tree and the specific gene loss of the XIP subfamily in physic nut. Furthermore, by analyzing deep transcriptome sequencing data, the expression evolution especially the expression divergence of duplicated HbAQP genes was also investigated and discussed. Results obtained from this study not only provide valuable information for future functional analysis and utilization of Jc/HbAQP genes, but also provide a useful reference to survey the gene family expansion and evolution in Euphorbiaceae plants and other plant species.

De novo assembly and functional annotation of Myrciaria dubia fruit transcriptome reveals multiple metabolic pathways for L-ascorbic acid biosynthesis

Background

Myrciaria dubia is an Amazonian fruit shrub that produces numerous bioactive phytochemicals, but is best known by its high L-ascorbic acid (AsA) content in fruits. Pronounced variation in AsA content has been observed both within and among individuals, but the genetic factors responsible for this variation are largely unknown. The goals of this research, therefore, were to assemble, characterize, and annotate the fruit transcriptome of M. dubia in order to reconstruct metabolic pathways and determine if multiple pathways contribute to AsA biosynthesis.

Results

In total 24,551,882 high-quality sequence reads were de novo assembled into 70,048 unigenes (mean length = 1150 bp, N50 = 1775 bp). Assembled sequences were annotated using BLASTX against public databases such as TAIR, GR-protein, FB, MGI, RGD, ZFIN, SGN, WB, TIGR_CMR, and JCVI-CMR with 75.2 % of unigenes having annotations. Of the three core GO annotation categories, biological processes comprised 53.6 % of the total assigned annotations, whereas cellular components and molecular functions comprised 23.3 and 23.1 %, respectively. Based on the KEGG pathway assignment of the functionally annotated transcripts, five metabolic pathways for AsA biosynthesis were identified: animal-like pathway, myo-inositol pathway, L-gulose pathway, D-mannose/L-galactose pathway, and uronic acid pathway. All transcripts coding enzymes involved in the ascorbate-glutathione cycle were also identified. Finally, we used the assembly to identified 6314 genic microsatellites and 23,481 high quality SNPs.

Conclusions

This study describes the first next-generation sequencing effort and transcriptome annotation of a non-model Amazonian plant that is relevant for AsA production and other bioactive phytochemicals. Genes encoding key enzymes were successfully identified and metabolic pathways involved in biosynthesis of AsA, anthocyanins, and other metabolic pathways have been reconstructed. The identification of these genes and pathways is in agreement with the empirically observed capability of M. dubia to synthesize and accumulate AsA and other important molecules, and adds to our current knowledge of the molecular biology and biochemistry of their production in plants. By providing insights into the mechanisms underpinning these metabolic processes, these results can be used to direct efforts to genetically manipulate this organism in order to enhance the production of these bioactive phytochemicals.

The accumulation of AsA precursor and discovery of genes associated with their biosynthesis and metabolism in M. dubia is intriguing and worthy of further investigation. The sequences and pathways produced here present the genetic framework required for further studies. Quantitative transcriptomics in concert with studies of the genome, proteome, and metabolome under conditions that stimulate production and accumulation of AsA and their precursors are needed to provide a more comprehensive view of how these pathways for AsA metabolism are regulated and linked in this species.

Electronic supplementary material

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

An efficient protocol for Agrobacterium-mediated transformation of the biofuel plant Jatropha curcas by optimizing kanamycin concentration and duration of delayed selection

Jatropha curcas is considered a potential biodiesel feedstock crop. Currently, the value of J. curcas is limited because its seed yield is generally low. Transgenic modification is a promising approach to improve the seed yield of J. curcas. Although Agrobacterium-mediated genetic transformation of J. curcas has been pursued for several years, the transformation efficiency remains unsatisfying. Therefore, a highly efficient and simple Agrobacterium-mediated genetic transformation method for J. curcas should be developed. We examined and optimized several key factors that affect genetic transformation of J. curcas in this study. The results showed that the EHA105 strain was superior to the other three Agrobacterium tumefaciens strains for infecting J. curcas cotyledons, and the supplementation of 100 mM acetosyringone slightly increased the transient transformation frequency. Use of the appropriate inoculation method, optimal kanamycin concentration and appropriate duration of delayed selection also improved the efficiency of stable genetic transformation of J. curcas. The percentage of β-glucuronidase positive J. curcas shoots reached as high as 56.0 %, and 1.70 transformants per explant were obtained with this protocol. Furthermore, we optimized the root-inducing medium to achieve a rooting rate of 84.9 %. Stable integration of the T-DNA into the genomes of putative transgenic lines was confirmed by PCR and Southern blot analysis. Using this improved protocol, a large number of transgenic J. curcas plantlets can be routinely obtained within approximately 4 months. The detailed information provided here for each step of J. curcas transformation should enable successful implementation of this transgenic technology in other laboratories.

Molecular marker development from transcript sequences and germplasm evaluation for cultivated peanut (Arachis hypogaea L.)

Molecular markers are important tools for genotyping in genetic studies and molecular breeding. The SSR and SNP are two commonly used marker systems developed from genomic or transcript sequences. The objectives of this study were to: (1) assemble and annotate the publicly available ESTs in Arachis and the in-house short reads, (2) develop and validate SSR and SNP markers, and (3) investigate the genetic diversity and population structure of the peanut breeding lines and the U.S. peanut mini core collection using developed SSR markers. An NCBI EST dataset with 252,951 sequences and an in-house 454 RNAseq dataset with 288,701 sequences were assembled separately after trimming. Transcript sequence comparison and phylogenetic analysis suggested that peanut is closer to cowpea and scarlet bean than to soybean, common bean and Medicago. From these two datasets, 6455 novel SSRs and 11,902 SNPs were identified. Of the discovered SSRs, 380 representing various SSR types were selected for PCR validation. The amplification rate was 89.2 %. Twenty-two (6.5 %) SSRs were polymorphic between at least one pair of four genotypes. Sanger sequencing of PCR products targeting 110 SNPs revealed 13 true SNPs between tetraploid genotypes and 193 homoeologous SNPs within genotypes. Eight out of the 22 polymorphic SSR markers were selected to evaluate the genetic diversity of Florida peanut breeding lines and the U.S. peanut mini core collection. This marker set demonstrated high discrimination power by displaying an average polymorphism information content value of 0.783, a combined probability of identity of 10(-11), and a combined power of exclusion of 0.99991. The structure analysis revealed four sub-populations among the peanut accessions and lines evaluated. The results of this study enriched the peanut genomic resources, provided over 6000 novel SSR markers and the credentials for true peanut SNP marker development, and demonstrated the power of newly developed SSR markers in genotyping peanut germplasm and breeding materials.

Regulation of FA and TAG biosynthesis pathway genes in endosperms and embryos of high and low oil content genotypes of Jatropha curcas L

The rising demand for biofuels has raised concerns about selecting alternate and promising renewable energy crops which do not compete with food supply. Jatropha (Jatropha curcas L.), a non-edible energy crop of the family euphorbiaceae, has the potential of providing biodiesel feedstock due to the presence of high proportion of unsaturated fatty acids (75%) in seed oil which is mainly accumulated in endosperm and embryo. The molecular basis of seed oil biosynthesis machinery has been studied in J. curcas, however, what genetic differences contribute to differential oil biosynthesis and accumulation in genotypes varying for oil content is poorly understood. We investigated expression profile of 18 FA and TAG biosynthetic pathway genes in different developmental stages of embryo and endosperm from high (42%) and low (30%) oil content genotypes grown at two geographical locations. Most of the genes showed relatively higher expression in endosperms of high oil content genotype, whereas no significant difference was observed in endosperms versus embryos of low oil content genotype. The promoter regions of key genes from FA and TAG biosynthetic pathways as well as other genes implicated in oil accumulation were analyzed for regulatory elements and transcription factors specific to oil or lipid accumulation in plants such as Dof, CBF (LEC1), SORLIP, GATA and Skn-1_motif etc. Identification of key genes from oil biosynthesis and regulatory elements specific to oil deposition will be useful not only in dissecting the molecular basis of high oil content but also improving seed oil content through transgenic or molecular breeding approaches.

Proteomic Analysis of the Endosperm Ontogeny of Jatropha curcas L. Seeds

Seeds of Jatropha curcas L. represent a potential source of raw material for the production of biodiesel. However, this use is hampered by the lack of basic information on the biosynthetic pathways associated with synthesis of toxic diterpenes, fatty acids, and triacylglycerols, as well as the pattern of deposition of storage proteins during seed development. In this study, we performed an in-depth proteome analysis of the endosperm isolated from five developmental stages which resulted in the identification of 1517, 1256, 1033, 752, and 307 proteins, respectively, summing up 1760 different proteins. Proteins with similar label free quantitation expression pattern were grouped into five clusters. The biological significance of these identifications is discussed with special focus on the analysis of seed storage proteins, proteins involved in the metabolism of fatty acids, carbohydrates, toxic components and proteolytic processing. Although several enzymes belonging to the biosynthesis of diterpenoid precursors were identified, we were unable to find any terpene synthase/cyclase, indicating that the synthesis of phorbol esters, the main toxic diterpenes, does not occur in seeds. The strategy used enabled us to provide a first in depth proteome analysis of the developing endosperm of this biodiesel plant, providing an important glimpse into the enzymatic machinery devoted to the production of C and N sources to sustain seed development.

Transcriptome Analysis Comparison of Lipid Biosynthesis in the Leaves and Developing Seeds of Brassica napus

Brassica napus seed is a lipid storage organ containing approximately 40% oil, while its leaves contain many kinds of lipids for many biological roles, but the overall amounts are less than in seeds. Thus, lipid biosynthesis in the developing seeds and the leaves is strictly regulated which results the final difference of lipids. However, there are few reports about the molecular mechanism controlling the difference in lipid biosynthesis between developing seeds and leaves. In this study, we tried to uncover this mechanism by analyzing the transcriptome data for lipid biosynthesis. The transcriptome data were de novo assembled and a total of 47,216 unigenes were obtained, which had an N50 length and median of 1271 and 755 bp, respectively. Among these unigenes, 36,368 (about 77.02%) were annotated and there were 109 up-regulated unigenes and 72 down-regulated unigenes in the developing seeds lipid synthetic pathway after comparing with leaves. In the oleic acid pathway, 23 unigenes were up-regulated and four unigenes were down-regulated. During triacylglycerol (TAG) synthesis, the key unigenes were all up-regulated, such as phosphatidate phosphatase and diacylglycerol O-acyltransferase. During palmitic acid, palmitoleic acid, stearic acid, linoleic acid and linolenic acid synthesis in leaves, the unigenes were nearly all up-regulated, which indicated that the biosynthesis of these particular fatty acids were more important in leaves. In the developing seeds, almost all the unigenes in the ABI3VP1, RKD, CPP, E2F-DP, GRF, JUMONJI, MYB-related, PHD and REM transcript factor families were up-regulated, which helped us to discern the regulation mechanism underlying lipid biosynthesis. The differential up/down-regulation of the genes and TFs involved in lipid biosynthesis in developing seeds and leaves provided direct evidence that allowed us to map the network that regulates lipid biosynthesis, and the identification of new TFs that are up-regulated in developing seeds will help us to further elucidate the lipids biosynthesis pathway in developing seeds and leaves.

Transcriptome of the inflorescence meristems of the biofuel plant Jatropha curcas treated with cytokinin

Background

Jatropha curcas, whose seed content is approximately 30–40% oil, is an ideal feedstock for producing biodiesel and bio-jet fuels. However, Jatropha plants have a low number of female flowers, which results in low seed yield that cannot meet the needs of the biofuel industry. Thus, increasing the number of female flowers is critical for the improvement of Jatropha seed yield. Our previous findings showed that cytokinin treatment can increase the flower number and female to male ratio and also induce bisexual flowers in Jatropha. The mechanisms underlying the influence of cytokinin on Jatropha flower development and sex determination, however, have not been clarified.

Results

This study examined the transcriptional levels of genes involved in the response to cytokinin in Jatropha inflorescence meristems at different time points after cytokinin treatment by 454 sequencing, which gave rise to a total of 294.6 Mb of transcript sequences. Up-regulated and down-regulated annotated and novel genes were identified, and the expression levels of the genes of interest were confirmed by qRT-PCR. The identified transcripts include those encoding genes involved in the biosynthesis, metabolism, and signaling of cytokinin and other plant hormones, flower development and cell division, which may be related to phenotypic changes of Jatropha in response to cytokinin treatment. Our analysis indicated that Jatropha orthologs of the floral organ identity genes known as ABCE model genes, JcAP1,2, JcPI, JcAG, and JcSEP1,2,3, were all significantly repressed, with an exception of one B-function gene JcAP3 that was shown to be up-regulated by BA treatment, indicating different mechanisms to be involved in the floral organ development of unisexual flowers of Jatropha and bisexual flowers of Arabidopsis. Several cell division-related genes, including JcCycA3;2, JcCycD3;1, JcCycD3;2 and JcTSO1, were up-regulated, which may contribute to the increased flower number after cytokinin treatment.

Conclusions

This study presents the first report of global expression patterns of cytokinin-regulated transcripts in Jatropha inflorescence meristems. This report laid the foundation for further mechanistic studies on Jatropha and other non-model plants responding to cytokinin. Moreover, the identification of functional candidate genes will be useful for generating superior varieties of high-yielding transgenic Jatropha.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-974) contains supplementary material, which is available to authorized users.

Multi-Spectroscopic Analysis of Seed Quality and 13C-Stable-Iotopologue Monitoring in Initial Growth Metabolism of Jatropha curcas L

In the present study, we applied nuclear magnetic resonance (NMR), as well as near-infrared (NIR) spectroscopy, to Jatropha curcas to fulfill two objectives: (1) to qualitatively examine the seeds stored at different conditions, and (2) to monitor the metabolism of J. curcas during its initial growth stage under stable-isotope-labeling condition (until 15 days after seeding). NIR spectra could non-invasively distinguish differences in storage conditions. NMR metabolic analysis of water-soluble metabolites identified sucrose and raffinose family oligosaccharides as positive markers and gluconic acid as a negative marker of seed germination. Isotopic labeling patteren of metabolites in germinated seedlings cultured in agar-plate containg 13C-glucose and 15N-nitrate was analyzed by zero-quantum-filtered-total correlation spectroscopy (ZQF-TOCSY) and 13C-detected 1H-13C heteronuclear correlation spectroscopy (HETCOR). 13C-detected HETOCR with 13C-optimized cryogenic probe provided high-resolution 13C-NMR spectra of each metabolite in molecular crowd. The 13C-13C/12C bondmer estimated from 1H-13C HETCOR spectra indicated that glutamine and arginine were the major organic compounds for nitrogen and carbon transfer from roots to leaves.

Proteomic analysis of oil bodies in mature Jatropha curcas seeds with different lipid content

To reveal the difference among three mature Jatropha curcas seeds (JcVH, variant with high lipid content; JcW, wild type and JcVL, variant with low lipid content) with different lipid content, comparative proteomics was employed to profile the changes of oil body (OB) associated protein species by using gels-based proteomic technique. Eighty-three protein species were successfully identified through LTQ-ES-MS/MS from mature JcW seeds purified OBs. Two-dimensional electrophoresis analysis of J. curcas OB associated protein species revealed they had essential interactions with other organelles and demonstrated that oleosin and caleosin were the most abundant OB structural protein species. Twenty-eight OB associated protein species showed significant difference among JcVH, JcW and JcVL according to statistical analysis. Complementary transient expression analysis revealed that calcium ion binding protein (CalBP) and glycine-rich RNA binding protein (GRP) were well targeted in OBs apart from the oleosins. This study demonstrated that ratio of lipid content to caleosins abundance was involved in the regulation of OB size, and the mutant induced by ethylmethylsulfone treatment might be related to the caleosin like protein species. These findings are important for biotechnological improvement with the aim to alter the lipid content in J. curcas seeds.

BIOLOGICAL SIGNIFICANCE

The economic value of Jatropha curcas largely depends on the lipid content in seeds which are mainly stored in the special organelle called oil bodies (OBs). In consideration of the biological importance and applications of J. curcas OB in seeds, it is necessary to further explore the components and functions of J. curcas OBs. Although a previous study concerning the J. curcas OB proteome revealed oleosins were the major OB protein component and additional protein species were similar to those in other oil seed plants, these identified OB associated protein species were corresponding to the protein bands instead of protein spots in the electrophoresis gels. Furthermore, the interaction of OB associated protein species and their contribution to OB formation and stabilization are still blank. In this study, with the overall object of profiling OB protein species from mature J. curcas seeds with different lipid content, we provided a setting of comparative OB proteomics with biochemical data and transient expression to explore the core of OB associated protein species involved in the regulation of OB size and lipid accumulation. The results were important for biotechnological improvement with the aim to a global modification of lipid storage in J. curcas seeds. Meanwhile, this study gave insight into possible associations between OBs and other organelles in mature J. curcas seeds. It may represent new aspects of the biological functions of the OBs during the oil mobilization. Combined the technique of transient transformation, a newly reported protein species, glycine-rich RNA binding protein (GRP) was successfully targeted in OBs. Therefore, further molecular analysis of these protein species is warranted to verify this association and what role they have in OBs.

Qualitative de novo analysis of full length cDNA and quantitative analysis of gene expression for common marmoset (Callithrix jacchus) transcriptomes using parallel long-read technology and short-read sequencing

The common marmoset (Callithrix jacchus) is a non-human primate that could prove useful as human pharmacokinetic and biomedical research models. The cytochromes P450 (P450s) are a superfamily of enzymes that have critical roles in drug metabolism and disposition via monooxygenation of a broad range of xenobiotics; however, information on some marmoset P450s is currently limited. Therefore, identification and quantitative analysis of tissue-specific mRNA transcripts, including those of P450s and flavin-containing monooxygenases (FMO, another monooxygenase family), need to be carried out in detail before the marmoset can be used as an animal model in drug development. De novo assembly and expression analysis of marmoset transcripts were conducted with pooled liver, intestine, kidney, and brain samples from three male and three female marmosets. After unique sequences were automatically aligned by assembling software, the mean contig length was 718 bp (with a standard deviation of 457 bp) among a total of 47,883 transcripts. Approximately 30% of the total transcripts were matched to known marmoset sequences. Gene expression in 18 marmoset P450- and 4 FMO-like genes displayed some tissue-specific patterns. Of these, the three most highly expressed in marmoset liver were P450 2D-, 2E-, and 3A-like genes. In extrahepatic tissues, including brain, gene expressions of these monooxygenases were lower than those in liver, although P450 3A4 (previously P450 3A21) in intestine and P450 4A11- and FMO1-like genes in kidney were relatively highly expressed. By means of massive parallel long-read sequencing and short-read technology applied to marmoset liver, intestine, kidney, and brain, the combined next-generation sequencing analyses reported here were able to identify novel marmoset drug-metabolizing P450 transcripts that have until now been little reported. These results provide a foundation for mechanistic studies and pave the way for the use of marmosets as model animals for drug development in the future.

Transcriptome analysis by Illumina high-throughout paired-end sequencing reveals the complexity of differential gene expression during in vitro plantlet growth and flowering in Amaranthus tricolor L

Amaranthus tricolor L. is a C4 plant, which is consumed as a major leafy vegetable in some tropical countries. Under conditions of high temperature and short daylight, Am. tricolor readily bolts and blooms, degrading leaf quality. A preliminary in vitro flowering study demonstrated that the flowering control pathway in Am. tricolor may differ from that of Arabidopsis. Nevertheless, no transcriptome analysis of the flowering process in Amaranthus has been conducted. To study Am. tricolor floral regulatory mechanisms, we conducted a large-scale transcriptome analysis--based on Illumina HiSeq sequencing of cDNA libraries generated from Am. tricolor at young seedling (YSS), adult seedling (ASS), flower bud (FBS), and flowering (FS) stages. A total of 99,312 unigenes were obtained. Using BLASTX, 43,088 unigenes (43.39%) were found to have significant similarity with accessions in Nr, Nt, and Swiss-Prot databases. Of these unigenes, 11,291 were mapped to 266 KEGG pathways. Further analysis of the four digital transcriptomes revealed that 735, 17,184, 274, and 206 unigenes were specifically expressed during YSS, ASS, FBS, and FS, respectively, with 59,517 unigenes expressed throughout the four stages. These unigenes were involved in many metabolic pathways related to in vitro flowering. Among these pathways, 259 unigenes were associated with ubiquitin-mediated proteolysis, indicating its importance for in vitro flowering in Am. tricolor. Other pathways, such as circadian rhythm and cell cycle, also had important roles. Finally, 26 unigenes were validated by qRT-PCR in samples from Am. tricolor at YSS, ASS, FBS, and FS; their differential expressions at the various stages indicate their possible roles in Am. tricolor growth and development, but the results were somewhat similar to Arabidopsis. Because unigenes involved in many metabolic pathways or of unknown function were revealed to regulate in vitro plantlet growth and flowering in Am. tricolor, the process appears to be highly complex in this species.

Rapid development of microsatellite markers for Callosobruchus chinensis using Illumina paired-end sequencing

Background

The adzuki bean weevil, Callosobruchus chinensis L., is one of the most destructive pests of stored legume seeds such as mungbean, cowpea, and adzuki bean, which usually cause considerable loss in the quantity and quality of stored seeds during transportation and storage. However, a lack of genetic information of this pest results in a series of genetic questions remain largely unknown, including population genetic structure, kinship, biotype abundance, and so on. Co-dominant microsatellite markers offer a great resolving power to determine these events. Here, we report rapid microsatellite isolation from C. chinensis via high-throughput sequencing.

Principal Findings

In this study, 94,560,852 quality-filtered and trimmed reads were obtained for the assembly of genome using Illumina paired-end sequencing technology. In total, the genome with total length of 497,124,785 bp, comprising 403,113 high quality contigs was generated with de novo assembly. More than 6800 SSR loci were detected and a suit of 6303 primer pair sequences were designed and 500 of them were randomly selected for validation. Of these, 196 pair of primers, i.e. 39.2%, produced reproducible amplicons that were polymorphic among 8 C. chinensis genotypes collected from different geographical regions. Twenty out of 196 polymorphic SSR markers were used to analyze the genetic diversity of 18 C. chinensis populations. The results showed the twenty SSR loci were highly polymorphic among these populations.

Conclusions

This study presents a first report of genome sequencing and de novo assembly for C. chinensis and demonstrates the feasibility of generating a large scale of sequence information and SSR loci isolation by Illumina paired-end sequencing. Our results provide a valuable resource for C. chinensis research. These novel markers are valuable for future genetic mapping, trait association, genetic structure and kinship among C. chinensis.

Global analysis of gene expression profiles in physic nut (Jatropha curcas L.) seedlings exposed to salt stress

BACKGROUND

Salt stress interferes with plant growth and production. Plants have evolved a series of molecular and morphological adaptations to cope with this abiotic stress, and overexpression of salt response genes reportedly enhances the productivity of various crops. However, little is known about the salt responsive genes in the energy plant physic nut (Jatropha curcas L.). Thus, excavate salt responsive genes in this plant are informative in uncovering the molecular mechanisms for the salt response in physic nut.

METHODOLOGY/PRINCIPAL FINDINGS

We applied next-generation Illumina sequencing technology to analyze global gene expression profiles of physic nut plants (roots and leaves) 2 hours, 2 days and 7 days after the onset of salt stress. A total of 1,504 and 1,115 genes were significantly up and down-regulated in roots and leaves, respectively, under salt stress condition. Gene ontology (GO) analysis of physiological process revealed that, in the physic nut, many "biological processes" were affected by salt stress, particular those categories belong to "metabolic process", such as "primary metabolism process", "cellular metabolism process" and "macromolecule metabolism process". The gene expression profiles indicated that the associated genes were responsible for ABA and ethylene signaling, osmotic regulation, the reactive oxygen species scavenging system and the cell structure in physic nut.

CONCLUSIONS/SIGNIFICANCE

The major regulated genes detected in this transcriptomic data were related to trehalose synthesis and cell wall structure modification in roots, while related to raffinose synthesis and reactive oxygen scavenger in leaves. The current study shows a comprehensive gene expression profile of physic nut under salt stress. The differential expression genes detected in this study allows the underling the salt responsive mechanism in physic nut with the aim of improving its salt resistance in the future.

Jatropha curcas, a biofuel crop: functional genomics for understanding metabolic pathways and genetic improvement

Jatropha curcas is currently attracting much attention as an oilseed crop for biofuel, as Jatropha can grow under climate and soil conditions that are unsuitable for food production. However, little is known about Jatropha, and there are a number of challenges to be overcome. In fact, Jatropha has not really been domesticated; most of the Jatropha accessions are toxic, which renders the seedcake unsuitable for use as animal feed. The seeds of Jatropha contain high levels of polyunsaturated fatty acids, which negatively impact the biofuel quality. Fruiting of Jatropha is fairly continuous, thus increasing costs of harvesting. Therefore, before starting any improvement program using conventional or molecular breeding techniques, understanding gene function and the genome scale of Jatropha are prerequisites. This review presents currently available and relevant information on the latest technologies (genomics, transcriptomics, proteomics and metabolomics) to decipher important metabolic pathways within Jatropha, such as oil and toxin synthesis. Further, it discusses future directions for biotechnological approaches in Jatropha breeding and improvement.

De novo transcriptome analysis of the Siberian apricot (Prunus sibirica L.) and search for potential SSR markers by 454 pyrosequencing

The Siberian apricot, an economically and ecologically important plant in China, contains seeds high in oil and can grow on marginal land. Although this species has multiple purposes and may be a feedstock of biofuel in China, transcriptome information and molecular research on this species remain limited. RNA-Seq technology has been widely applied to transcriptomics, genomics and the development of molecular markers, and functional gene studies. In this study, we obtained 1,243,067 high-quality reads with a mean size of 425 bp in a single run, totaling 528.4 Mb of sequence data using 454 GS FLX Titanium sequencing. All reads were assembled de novo into 46,940 unigenes with a mean size of 651 bp (range: 45-5566 bp). Assembled unigenes were annotated in multiple public databases based on similarity alignments to genes and proteins. 191 unigenes involving in lipid biosynthesis and metabolism were found, among them, expression patterns of two desaturase enzymes were analyzed by quantitative real-time polymerase chain reaction (qRT-PCR), based on six tissues from Siberian apricot, the seeds had the highest expression. 7304 simple sequence repeats (SSR) were identified from 6509 unigenes, a total of 9930 primer pairs were designed, 50 primer pairs were randomly selected to validate of the usefulness, and 24 (48%) primer pairs produced bands of the expected size. These data provide a base of sequence information to improve agronomic characters and molecular marker-assisted breeding to alter the composition of fatty acids in seeds from this plant, and hence, facilitate its utilization as a future biodiesel feedstock.

Metatranscriptome analysis of fungal strains Penicillium camemberti and Geotrichum candidum reveal cheese matrix breakdown and potential development of sensory properties of ripened Camembert-type cheese

BACKGROUND

Camembert-type cheese ripening is driven mainly by fungal microflora including Geotrichum candidum and Penicillium camemberti. These species are major contributors to the texture and flavour of typical bloomy rind cheeses. Biochemical studies showed that G. candidum reduces bitterness, enhances sulphur flavors through amino acid catabolism and has an impact on rind texture, firmness and thickness, while P. camemberti is responsible for the white and bloomy aspect of the rind, and produces enzymes involved in proteolysis and lipolysis activities. However, very little is known about the genetic determinants that code for these activities and their expression profile over time during the ripening process.

RESULTS

The metatranscriptome of an industrial Canadian Camembert-type cheese was studied at seven different sampling days over 77 days of ripening. A database called CamemBank01 was generated, containing a total of 1,060,019 sequence tags (reads) assembled in 7916 contigs. Sequence analysis revealed that 57% of the contigs could be affiliated to molds, 16% originated from yeasts, and 27% could not be identified. According to the functional annotation performed, the predominant processes during Camembert ripening include gene expression, energy-, carbohydrate-, organic acid-, lipid- and protein- metabolic processes, cell growth, and response to different stresses. Relative expression data showed that these functions occurred mostly in the first two weeks of the ripening period.

CONCLUSIONS

These data provide further advances in our knowledge about the biological activities of the dominant ripening microflora of Camembert cheese and will help select biological markers to improve cheese quality assessment.

Peroxidase gene discovery from the horseradish transcriptome

BACKGROUND

Horseradish peroxidases (HRPs) from Armoracia rusticana have long been utilized as reporters in various diagnostic assays and histochemical stainings. Regardless of their increasing importance in the field of life sciences and suggested uses in medical applications, chemical synthesis and other industrial applications, the HRP isoenzymes, their substrate specificities and enzymatic properties are poorly characterized. Due to lacking sequence information of natural isoenzymes and the low levels of HRP expression in heterologous hosts, commercially available HRP is still extracted as a mixture of isoenzymes from the roots of A. rusticana.

RESULTS

In this study, a normalized, size-selected A. rusticana transcriptome library was sequenced using 454 Titanium technology. The resulting reads were assembled into 14871 isotigs with an average length of 1133 bp. Sequence databases, ORF finding and ORF characterization were utilized to identify peroxidase genes from the 14871 isotigs generated by de novo assembly. The sequences were manually reviewed and verified with Sanger sequencing of PCR amplified genomic fragments, resulting in the discovery of 28 secretory peroxidases, 23 of them previously unknown. A total of 22 isoenzymes including allelic variants were successfully expressed in Pichia pastoris and showed peroxidase activity with at least one of the substrates tested, thus enabling their development into commercial pure isoenzymes.

CONCLUSIONS

This study demonstrates that transcriptome sequencing combined with sequence motif search is a powerful concept for the discovery and quick supply of new enzymes and isoenzymes from any plant or other eukaryotic organisms. Identification and manual verification of the sequences of 28 HRP isoenzymes do not only contribute a set of peroxidases for industrial, biological and biomedical applications, but also provide valuable information on the reliability of the approach in identifying and characterizing a large group of isoenzymes.

Single-nucleotide polymorphism markers from de-novo assembly of the pomegranate transcriptome reveal germplasm genetic diversity

Pomegranate is a valuable crop that is grown commercially in many parts of the world. Wild species have been reported from India, Turkmenistan and Socotra. Pomegranate fruit has a variety of health-beneficial qualities. However, despite this crop's importance, only moderate effort has been invested in studying its biochemical or physiological properties or in establishing genomic and genetic infrastructures. In this study, we reconstructed a transcriptome from two phenotypically different accessions using 454-GS-FLX Titanium technology. These data were used to explore the functional annotation of 45,187 fully annotated contigs. We further compiled a genetic-variation resource of 7,155 simple-sequence repeats (SSRs) and 6,500 single-nucleotide polymorphisms (SNPs). A subset of 480 SNPs was sampled to investigate the genetic structure of the broad pomegranate germplasm collection at the Agricultural Research Organization (ARO), which includes accessions from different geographical areas worldwide. This subset of SNPs was found to be polymorphic, with 10.7% loci with minor allele frequencies of (MAF<0.05). These SNPs were successfully used to classify the ARO pomegranate collection into two major groups of accessions: one from India, China and Iran, composed of mainly unknown country origin and which was more of an admixture than the other major group, composed of accessions mainly from the Mediterranean basin, Central Asia and California. This study establishes a high-throughput transcriptome and genetic-marker infrastructure. Moreover, it sheds new light on the genetic interrelations between pomegranate species worldwide and more accurately defines their genetic nature.

RNA sequencing read depth requirement for optimal transcriptome coverage in Hevea brasiliensis

BACKGROUND

One of the concerns of assembling de novo transcriptomes is determining the amount of read sequences required to ensure a comprehensive coverage of genes expressed in a particular sample. In this report, we describe the use of Illumina paired-end RNA-Seq (PE RNA-Seq) reads from Hevea brasiliensis (rubber tree) bark to devise a transcript mapping approach for the estimation of the read amount needed for deep transcriptome coverage.

FINDINGS

We optimized the assembly of a Hevea bark transcriptome based on 16 Gb Illumina PE RNA-Seq reads using the Oases assembler across a range of k-mer sizes. We then assessed assembly quality based on transcript N50 length and transcript mapping statistics in relation to (a) known Hevea cDNAs with complete open reading frames, (b) a set of core eukaryotic genes and (c) Hevea genome scaffolds. This was followed by a systematic transcript mapping process where sub-assemblies from a series of incremental amounts of bark transcripts were aligned to transcripts from the entire bark transcriptome assembly. The exercise served to relate read amounts to the degree of transcript mapping level, the latter being an indicator of the coverage of gene transcripts expressed in the sample. As read amounts or datasize increased toward 16 Gb, the number of transcripts mapped to the entire bark assembly approached saturation. A colour matrix was subsequently generated to illustrate sequencing depth requirement in relation to the degree of coverage of total sample transcripts.

CONCLUSIONS

We devised a procedure, the "transcript mapping saturation test", to estimate the amount of RNA-Seq reads needed for deep coverage of transcriptomes. For Hevea de novo assembly, we propose generating between 5-8 Gb reads, whereby around 90% transcript coverage could be achieved with optimized k-mers and transcript N50 length. The principle behind this methodology may also be applied to other non-model plants, or with reads from other second generation sequencing platforms.

Comparative de novo transcriptome analysis and metabolic pathway studies of Citrus paradisi flavedo from naive stage to ripened stage

Grapefruit (Citrus pardisi) is a popular citrus fruit that is a cross between a sweet orange and pummelo. This research article focuses on an in silico approach for comparative analysis of C. paradisi green flavedo (GF) and ethylene treated flavedo (ETF) transcriptome data. Our pathway analysis provides comprehensive information of genes playing significant role in different stages of ripening in fruit. De novo assembly was carried out using six different assemblers namely GS assembler, SeqMan NGEN, Velvet/Oases, CLC, iAssembler and Cortex followed by subsequent meta-assembly, annotation and pathway analysis. We conclude that de novo transcriptome assembly using meta-assembly approach is used to increase assembly quality in comparison to single assembler.

De novo transcriptome analysis of an imminent biofuel crop, Camelina sativa L. using Illumina GAIIX sequencing platform and identification of SSR markers

Camelina sativa L. is an emerging biofuel crop with potential applications in industry, medicine, cosmetics and human nutrition. The crop is unexploited owing to very limited availability of transcriptome and genomic data. In order to analyse the various metabolic pathways, we performed de novo assembly of the transcriptome on Illumina GAIIX platform with paired end sequencing for obtaining short reads. The sequencing output generated a FastQ file size of 2.97 GB with 10.83 million reads having a maximum read length of 101 nucleotides. The number of contigs generated was 53,854 with maximum and minimum lengths of 10,086 and 200 nucleotides respectively. These trancripts were annotated using BLAST search against the Aracyc, Swiss-Prot, TrEMBL, gene ontology and clusters of orthologous groups (KOG) databases. The genes involved in lipid metabolism were studied and the transcription factors were identified. Sequence similarity studies of Camelina with the other related organisms indicated the close relatedness of Camelina with Arabidopsis. In addition, bioinformatics analysis revealed the presence of a total of 19,379 simple sequence repeats. This is the first report on Camelina sativa L., where the transcriptome of the entire plant, including seedlings, seed, root, leaves and stem was done. Our data established an excellent resource for gene discovery and provide useful information for functional and comparative genomic studies in this promising biofuel crop.

Comparative analyses of two Geraniaceae transcriptomes using next-generation sequencing

BACKGROUND

Organelle genomes of Geraniaceae exhibit several unusual evolutionary phenomena compared to other angiosperm families including accelerated nucleotide substitution rates, widespread gene loss, reduced RNA editing, and extensive genomic rearrangements. Since most organelle-encoded proteins function in multi-subunit complexes that also contain nuclear-encoded proteins, it is likely that the atypical organellar phenomena affect the evolution of nuclear genes encoding organellar proteins. To begin to unravel the complex co-evolutionary interplay between organellar and nuclear genomes in this family, we sequenced nuclear transcriptomes of two species, Geranium maderense and Pelargonium x hortorum.

RESULTS

Normalized cDNA libraries of G. maderense and P. x hortorum were used for transcriptome sequencing. Five assemblers (MIRA, Newbler, SOAPdenovo, SOAPdenovo-trans [SOAPtrans], Trinity) and two next-generation technologies (454 and Illumina) were compared to determine the optimal transcriptome sequencing approach. Trinity provided the highest quality assembly of Illumina data with the deepest transcriptome coverage. An analysis to determine the amount of sequencing needed for de novo assembly revealed diminishing returns of coverage and quality with data sets larger than sixty million Illumina paired end reads for both species. The G. maderense and P. x hortorum transcriptomes contained fewer transcripts encoding the PLS subclass of PPR proteins relative to other angiosperms, consistent with reduced mitochondrial RNA editing activity in Geraniaceae. In addition, transcripts for all six plastid targeted sigma factors were identified in both transcriptomes, suggesting that one of the highly divergent rpoA-like ORFs in the P. x hortorum plastid genome is functional.

CONCLUSIONS

The findings support the use of the Illumina platform and assemblers optimized for transcriptome assembly, such as Trinity or SOAPtrans, to generate high-quality de novo transcriptomes with broad coverage. In addition, results indicated no major improvements in breadth of coverage with data sets larger than six billion nucleotides or when sampling RNA from four tissue types rather than from a single tissue. Finally, this work demonstrates the power of cross-compartmental genomic analyses to deepen our understanding of the correlated evolution of the nuclear, plastid, and mitochondrial genomes in plants.

Global analysis of transcriptome responses and gene expression profiles to cold stress of Jatropha curcas L

BACKGROUND

Jatropha curcas L., also called the Physic nut, is an oil-rich shrub with multiple uses, including biodiesel production, and is currently exploited as a renewable energy resource in many countries. Nevertheless, because of its origin from the tropical MidAmerican zone, J. curcas confers an inherent but undesirable characteristic (low cold resistance) that may seriously restrict its large-scale popularization. This adaptive flaw can be genetically improved by elucidating the mechanisms underlying plant tolerance to cold temperatures. The newly developed Illumina Hiseq™ 2000 RNA-seq and Digital Gene Expression (DGE) are deep high-throughput approaches for gene expression analysis at the transcriptome level, using which we carefully investigated the gene expression profiles in response to cold stress to gain insight into the molecular mechanisms of cold response in J. curcas.

RESULTS

In total, 45,251 unigenes were obtained by assembly of clean data generated by RNA-seq analysis of the J. curcas transcriptome. A total of 33,363 and 912 complete or partial coding sequences (CDSs) were determined by protein database alignments and ESTScan prediction, respectively. Among these unigenes, more than 41.52% were involved in approximately 128 known metabolic or signaling pathways, and 4,185 were possibly associated with cold resistance. DGE analysis was used to assess the changes in gene expression when exposed to cold condition (12°C) for 12, 24, and 48 h. The results showed that 3,178 genes were significantly upregulated and 1,244 were downregulated under cold stress. These genes were then functionally annotated based on the transcriptome data from RNA-seq analysis.

CONCLUSIONS

This study provides a global view of transcriptome response and gene expression profiling of J. curcas in response to cold stress. The results can help improve our current understanding of the mechanisms underlying plant cold resistance and favor the screening of crucial genes for genetically enhancing cold resistance in J. curcas.

Functional characterization of two microsomal fatty acid desaturases from Jatropha curcas L

Linoleic acid (LA, C18:2) and α-linolenic acid (ALA, C18:3) are polyunsaturated fatty acids (PUFAs) and major storage compounds in plant seed oils. Microsomal ω-6 and ω-3 fatty acid (FA) desaturases catalyze the synthesis of seed oil LA and ALA, respectively. Jatropha curcas L. seed oils contain large proportions of LA, but very little ALA. In this study, two microsomal desaturase genes, named JcFAD2 and JcFAD3, were isolated from J. curcas. Both deduced amino acid sequences possessed eight histidines shown to be essential for desaturases activity, and contained motif in the C-terminal for endoplasmic reticulum localization. Heterologous expression in Saccharomyces cerevisiae and Arabidopsis thaliana confirmed that the isolated JcFAD2 and JcFAD3 proteins could catalyze LA and ALA synthesis, respectively. The results indicate that JcFAD2 and JcFAD3 are functional in controlling PUFA contents of seed oils and could be exploited in the genetic engineering of J. curcas, and potentially other plants.

Linkage mapping in the oilseed crop Jatropha curcas L. reveals a locus controlling the biosynthesis of phorbol esters which cause seed toxicity

Current efforts to grow the tropical oilseed crop Jatropha curcas L. economically are hampered by the lack of cultivars and the presence of toxic phorbol esters (PE) within the seeds of most provenances. These PE restrict the conversion of seed cake into animal feed, although naturally occurring 'nontoxic' provenances exist which produce seed lacking PE. As an important step towards the development of genetically improved varieties of J. curcas, we constructed a linkage map from four F₂ mapping populations. The consensus linkage map contains 502 codominant markers, distributed over 11 linkage groups, with a mean marker density of 1.8 cM per unique locus. Analysis of the inheritance of PE biosynthesis indicated that this is a maternally controlled dominant monogenic trait. This maternal control is due to biosynthesis of the PE occurring only within maternal tissues. The trait segregated 3 : 1 within seeds collected from F₂ plants, and QTL analysis revealed that a locus on linkage group 8 was responsible for phorbol ester biosynthesis. By taking advantage of the draft genome assemblies of J. curcas and Ricinus communis (castor), a comparative mapping approach was used to develop additional markers to fine map this mutation within 2.3 cM. The linkage map provides a framework for the dissection of agronomic traits in J. curcas, and the development of improved varieties by marker-assisted breeding. The identification of the locus responsible for PE biosynthesis means that it is now possible to rapidly breed new nontoxic varieties.

Transcriptome analysis of yellow horn (Xanthoceras sorbifolia Bunge): a potential oil-rich seed tree for biodiesel in China

Background

Yellow horn (Xanthoceras sorbifolia Bunge) is an oil-rich seed shrub that grows well in cold, barren environments and has great potential for biodiesel production in China. However, the limited genetic data means that little information about the key genes involved in oil biosynthesis is available, which limits further improvement of this species. In this study, we describe sequencing and de novo transcriptome assembly to produce the first comprehensive and integrated genomic resource for yellow horn and identify the pathways and key genes related to oil accumulation. In addition, potential molecular markers were identified and compiled.

Methodology/Principal Findings

Total RNA was isolated from 30 plants from two regions, including buds, leaves, flowers and seeds. Equal quantities of RNA from these tissues were pooled to construct a cDNA library for 454 pyrosequencing. A total of 1,147,624 high-quality reads with total and average lengths of 530.6 Mb and 462 bp, respectively, were generated. These reads were assembled into 51,867 unigenes, corresponding to a total of 36.1 Mb with a mean length, N50 and median of 696, 928 and 570 bp, respectively. Of the unigenes, 17,541 (33.82%) were unmatched in any public protein databases. We identified 281 unigenes that may be involved in de novo fatty acid (FA) and triacylglycerol (TAG) biosynthesis and metabolism. Furthermore, 6,707 SSRs, 16,925 SNPs and 6,201 InDels with high-confidence were also identified in this study.

Conclusions

This transcriptome represents a new functional genomics resource and a foundation for further studies on the metabolic engineering of yellow horn to increase oil content and modify oil composition. The potential molecular markers identified in this study provide a basis for polymorphism analysis of Xanthoceras, and even Sapindaceae; they will also accelerate the process of breeding new varieties with better agronomic characteristics.

Analysis of the global transcriptome of longan (Dimocarpus longan Lour.) embryogenic callus using Illumina paired-end sequencing

Background

Longan is a tropical/subtropical fruit tree of great economic importance in Southeast Asia. Progress in understanding molecular mechanisms of longan embryogenesis, which is the primary influence on fruit quality and yield, is slowed by lack of transcriptomic and genomic information. Illumina second generation sequencing, which is suitable for generating enormous numbers of transcript sequences that can be used for functional genomic analysis of longan.

Results

In this study, a longan embryogenic callus (EC) cDNA library was sequenced using an Illumina HiSeq 2000 system. A total of 64,876,258 clean reads comprising 5.84 Gb of nucleotides were assembled into 68,925 unigenes of 448-bp mean length, with unigenes ≥1000 bp accounting for 8.26% of the total. Using BLASTx, 40,634 unigenes were found to have significant similarity with accessions in Nr and Swiss- Prot databases. Of these, 38,845 unigenes were assigned to 43 GO sub-categories and 17,118 unigenes were classified into 25 COG sub-groups. In addition, 17,306 unigenes mapped to 199 KEGG pathways, with the categories of Metabolic pathways, Plant-pathogen interaction, Biosynthesis of secondary metabolites, and Genetic information processing being well represented. Analyses of unigenes ≥1000 bp revealed 328 embryogenesis-related unigenes as well as numerous unigenes expressed in EC associated with functions of reproductive growth, such as flowering, gametophytogenesis, and fertility, and vegetative growth, such as root and shoot growth. Furthermore, 23 unigenes related to embryogenesis and reproductive and vegetative growth were validated by quantitative real time PCR (qPCR) in samples from different stages of longan somatic embryogenesis (SE); their differentially expressions in the various embryogenic cultures indicated their possible roles in longan SE.

Conclusions

The quantity and variety of expressed EC genes identified in this study is sufficient to serve as a global transcriptome dataset for longan EC and to provide more molecular resources for longan functional genomics.

Proteomic analysis of the seed development in Jatropha curcas: from carbon flux to the lipid accumulation

To characterize the metabolic signatures of lipid accumulation in Jatropha curcas seeds, comparative proteomic technique was employed to profile protein changes during the seed development. Temporal changes in comparative proteome were examined using gels-based proteomic technique at six developmental stages for lipid accumulation. And 104 differentially expressed proteins were identified by MALDI-TOF/TOF tandem mass spectrometry. These protein species were classified into 10 functional categories, and the results demonstrated that protein species related to energy and metabolism were notably accumulated and involved in the carbon flux to lipid accumulation that occurs primarily from early to late stage in seed development. Glycolysis and oxidative pentose phosphate pathways were the major pathways of producing carbon flux, and the glucose-6-phosphate and triose-phosphate are the major carbon source for fatty acid synthesis. Lipid analysis revealed that fatty acid accumulation initiated 25days after flowering at the late stage of seed development of J. curcas. Furthermore, C16:0 was initially synthesized as the precursor for the elongation to C18:1 and C18:2 in the developing seeds of J. curcas. Together, the metabolic signatures on protein changes in seed development provide profound knowledge and perspective insights into understanding lipid network in J. curcas.

BIOLOGICAL SIGNIFICANCE

Due to the abundant oil content in seeds, Jatropha curcas seeds are being considered as the ideal materials for biodiesel. Although several studies had carried out the transcriptomic project to study the genes expression profiles in seed development of J. curcas, these ESTs hadn't been confirmed by qRT-PCR. Yet, the seed development of J. curcas had been described for a pool of developing seeds instead of being characterized systematically. Moreover, cellular metabolic events are also controlled by protein-protein interactions, posttranslational protein modifications, and enzymatic activities which cannot be described by transcriptional profiling approaches alone. In this study, within the overall objective of profiling differential protein abundance in developing J. curcas seeds, we provide a setting of physiological data with dynamic proteomic and qRT-PCR analysis to characterize the metabolic pathways and the relationship between mRNA and protein patterns from early stage to seed filling during the seed development of J. curcas. The construction of J. curcas seed development proteome profiles will significantly increase our understanding of the process of seed development and provide a foundation to examine the dynamic changes of the metabolic network during seed development process and certainly suggest some clues to improve the lipid content of J. curcas seeds.

Biofuels as a sustainable energy source: an update of the applications of proteomics in bioenergy crops and algae

Sustainable energy is the need of the 21st century, not because of the numerous environmental and political reasons but because it is necessary to human civilization's energy future. Sustainable energy is loosely grouped into renewable energy, energy conservation, and sustainable transport disciplines. In this review, we deal with the renewable energy aspect focusing on the biomass from bioenergy crops to microalgae to produce biofuels to the utilization of high-throughput omics technologies, in particular proteomics in advancing our understanding and increasing biofuel production. We look at biofuel production by plant- and algal-based sources, and the role proteomics has played therein. This article is part of a Special Issue entitled: Translational Plant Proteomics.

De novo sequence assembly and characterization of Lycoris aurea transcriptome using GS FLX titanium platform of 454 pyrosequencing

Background

Lycoris aurea, also called Golden Magic Lily, is an ornamentally and medicinally important species of the Amaryllidaceae family. To date, the sequencing of its whole genome is unavailable as a non-model organism. Transcriptomic information is also scarce for this species. In this study, we performed de novo transcriptome sequencing to produce the first comprehensive expressed sequence tag (EST) dataset for L. aurea using high-throughput sequencing technology.

Methodology and Principal Findings

Total RNA was isolated from leaves with sodium nitroprusside (SNP), salicylic acid (SA), or methyl jasmonate (MeJA) treatment, stems, and flowers at the bud, blooming, and wilting stages. Equal quantities of RNA from each tissue and stage were pooled to construct a cDNA library. Using 454 pyrosequencing technology, a total of 937,990 high quality reads (308.63 Mb) with an average read length of 329 bp were generated. Clustering and assembly of these reads produced a non-redundant set of 141,111 unique sequences, comprising 24,604 contigs and 116,507 singletons. All of the unique sequences were involved in the biological process, cellular component and molecular function categories by GO analysis. Potential genes and their functions were predicted by KEGG pathway mapping and COG analysis. Based on our sequence analysis and published literatures, many putative genes involved in Amaryllidaceae alkaloids synthesis, including PAL, TYDC OMT, NMT, P450, and other potentially important candidate genes, were identified for the first time in this Lycoris. Furthermore, 6,386 SSRs and 18,107 high-confidence SNPs were identified in this EST dataset.

Conclusions

The transcriptome provides an invaluable new data for a functional genomics resource and future biological research in L. aurea. The molecular markers identified in this study will provide a material basis for future genetic linkage and quantitative trait loci analyses, and will provide useful information for functional genomic research in future.

Identification of the Hevea brasiliensis AP2/ERF superfamily by RNA sequencing

Background

Rubber tree (Hevea brasiliensis) laticifers are the source of natural rubber. Rubber production depends on endogenous and exogenous ethylene (ethephon). AP2/ERF transcription factors, and especially Ethylene-Response Factors, play a crucial role in plant development and response to biotic and abiotic stresses. This study set out to sequence transcript expressed in various tissues using next-generation sequencing and to identify AP2/ERF superfamily in the rubber tree.

Results

The 454 sequencing technique was used to produce five tissue-type transcript libraries (leaf, bark, latex, embryogenic tissues and root). Reads from all libraries were pooled and reassembled to improve mRNA lengths and produce a global library. One hundred and seventy-three AP2/ERF contigs were identified by in silico analysis based on the amino acid sequence of the conserved AP2 domain from the global library. The 142 contigs with the full AP2 domain were classified into three main families (20 AP2 members, 115 ERF members divided into 11 groups, and 4 RAV members) and 3 soloist members. Fifty-nine AP2/ERF transcripts were found in latex. Alongside the microRNA172 already described in plants, eleven additional microRNAs were predicted to inhibit Hevea AP2/ERF transcripts.

Conclusions

Hevea has a similar number of AP2/ERF genes to that of other dicot species. We adapted the alignment and classification methods to data from next-generation sequencing techniques to provide reliable information. We observed several specific features for the ERF family. Three HbSoloist members form a group in Hevea. Several AP2/ERF genes highly expressed in latex suggest they have a specific function in Hevea. The analysis of AP2/ERF transcripts in Hevea presented here provides the basis for studying the molecular regulation of latex production in response to abiotic stresses and latex cell differentiation.