The plant shoot system consists of reproductive organs such as inflorescences, buds and fruits, and the vegetative leaves and stems. In this study, the reproductive part of the Jatropha curcas shoot system, which includes the aerial shoots, shoots bearing the inflorescence and inflorescence were investigated in regard to gene-to-gene interactions underpinning yield-related biological processes. An RNA-seq based sequencing of shoot tissues performed on an Illumina HiSeq. 2500 platform generated 18 transcriptomes. Using the reference genome-based mapping approach, a total of 64 361 genes was identified in all samples and the data was annotated against the non-redundant database by the BLAST2GO Pro. Suite. After removing the outlier genes and samples, a total of 12 734 genes across 17 samples were subjected to gene co-expression network construction using petal, an R library. A gene co-expression network model built with scale-free and small-world properties extracted four vicinity networks (VNs) with putative involvement in yield-related biological processes as follow; heat stress tolerance, floral and shoot meristem differentiation, biosynthesis of chlorophyll molecules and laticifers, cell wall metabolism and epigenetic regulations. Our VNs revealed putative key players that could be adapted in breeding strategies for J. curcas shoot system improvements.
Jatropha curcas is an oil-rich seed crop with huge potentials for bioenergy production. The inflorescence carries a number of processes that are likely to affect the overall yield potentials; floral development, male-to-female flower ratio, floral abscission and fruit set. In this study, a weighted gene co-expression network analysis which integrates the transcriptome, physical and simple sugar data of J. curcas inflorescence was performed and nine modules were identified by means of hierarchical clustering. Among them, four modules (green4, antiquewhite2, brown2 and lightskyblue4) showed significant correlation to yield factors at p≤0.01. The four modules are categorized into two clusters; cluster 1 of green4 and antiquewhite2 modules correspond to number of flowers/inflorescence, total seed weight/plant, number of seeds/plant, and number of fruits/plant, whereas cluster 2 of brown2 and lightskyblue4 modules correspond to glucose and fructose. Descriptive characterizations of cluster 1 show putative involvement in gibberellin signaling and responses, whereas cluster 2 may have been involved in sugar signaling, signal transductions and regulation of flowerings. Our findings present a list of hub genes for J. curcas yield improvement and reproductive biology enhancement strategies.
Rafflesia is a biologically enigmatic species that is very rare in occurrence and possesses an extraordinary morphology. This parasitic plant produces a gigantic flower up to one metre in diameter with no leaves, stem or roots. However, little is known about the floral biology of this species especially at the molecular level. In an effort to address this issue, we have generated and characterised the transcriptome of the Rafflesia cantleyi flower, and performed a comparison with the transcriptome of its floral bud to predict genes that are expressed and regulated during flower development. Approximately 40 million sequencing reads were generated and assembled de novo into 18,053 transcripts with an average length of 641 bp. Of these, more than 79% of the transcripts had significant matches to annotated sequences in the public protein database. A total of 11,756 and 7,891 transcripts were assigned to Gene Ontology categories and clusters of orthologous groups respectively. In addition, 6,019 transcripts could be mapped to 129 pathways in Kyoto Encyclopaedia of Genes and Genomes Pathway database. Digital abundance analysis identified 52 transcripts with very high expression in the flower transcriptome of R. cantleyi. Subsequently, analysis of differential expression between developing flower and the floral bud revealed a set of 105 transcripts with potential role in flower development. Our work presents a deep transcriptome resource analysis for the developing flower of R. cantleyi. Genes potentially involved in the growth and development of the R. cantleyi flower were identified and provide insights into biological processes that occur during flower development.
Rafflesia possesses unique biological features and known primarily for producing the world's largest and existing as a single flower. However, to date, little is known about key regulators participating in Rafflesia flower development. In order to further understand the molecular mechanism that regulates Rafflesia cantleyi flower development, RNA-seq data from three developmental stages of floral bud, representing the floral organ primordia initiation, floral organ differentiation, and floral bud outgrowth, were analysed. A total of 89,890 transcripts were assembled of which up to 35% could be annotated based on homology search. Advanced transcriptome analysis using K-mean clustering on the differentially expressed genes (DEGs) was able to identify 12 expression clusters that reflect major trends and key transitional states, which correlate to specific developmental stages. Through this, comparative gene expression analysis of different floral bud stages identified various transcription factors related to flower development. The members of WRKY, NAC, bHLH, and MYB families are the most represented among the DEGs, suggesting their important function in flower development. Furthermore, pathway enrichment analysis also revealed DEGs that are involved in various phytohormone signal transduction events such as auxin and auxin transport, cytokinin and gibberellin biosynthesis. Results of this study imply that transcription factors and phytohormone signalling pathways play major role in Rafflesia floral bud development. This study provides an invaluable resource for molecular studies of the flower development process in Rafflesia and other plant species.
Rafflesia is a unique plant species existing as a single flower and produces the largest flower in the world. While Rafflesia buds take up to 21 months to develop, its flowers bloom and wither within about a week. In this study, transcriptome analysis was carried out to shed light on the molecular mechanism of senescence in Rafflesia. A total of 53.3 million high quality reads were obtained from two Rafflesia cantleyi flower developmental stages and assembled to generate 64,152 unigenes. Analysis of this dataset showed that 5,166 unigenes were differentially expressed, in which 1,073 unigenes were identified as genes involved in flower senescence. Results revealed that as the flowers progress to senescence, more genes related to flower senescence were significantly over-represented compared to those related to plant growth and development. Senescence of the R. cantleyi flower activates senescence-associated genes in the transcription activity (members of the transcription factor families MYB, bHLH, NAC, and WRKY), nutrient remobilization (autophagy-related protein and transporter genes), and redox regulation (CATALASE). Most of the senescence-related genes were found to be differentially regulated, perhaps for the fine-tuning of various responses in the senescing R. cantleyi flower. Additionally, pathway analysis showed the activation of genes such as ETHYLENE RECEPTOR, ETHYLENE-INSENSITIVE 2, ETHYLENE-INSENSITIVE 3, and ETHYLENE-RESPONSIVE TRANSCRIPTION FACTOR, indicating the possible involvement of the ethylene hormone response pathway in the regulation of R. cantleyi senescence. Our results provide a model of the molecular mechanism underlying R. cantleyi flower senescence, and contribute essential information towards further understanding the biology of the Rafflesiaceae family.
Dehydration-responsive element binding (DREB) transcription factor plays an important role in controlling the expression of abiotic stress responsive genes. An intronless oil palm EgDREB1 was isolated and confirmed to be a nuclear localized protein. Electrophoretic mobility shift and yeast one-hybrid assays validated its ability to interact with DRE/CRT motif. Its close evolutionary relation to the dicot NtDREB2 suggests a universal regulatory role. In order to determine its involvement in abiotic stress response, functional characterization was performed in oil palm seedlings subjected to different levels of drought severity and in EgDREB1 transgenic tomato seedlings treated by abiotic stresses. Its expression in roots and leaves was compared with several antioxidant genes using quantitative real-time PCR. Early accumulation of EgDREB1 in oil palm roots under mild drought suggests possible involvement in the initiation of signaling communication from root to shoot. Ectopic expression of EgDREB1 in T1 transgenic tomato seedlings enhanced expression of DRE/CRT and non-DRE/CRT containing genes, including tomato peroxidase (LePOD), ascorbate peroxidase (LeAPX), catalase (LeCAT), superoxide dismutase (LeSOD), glutathione reductase (LeGR), glutathione peroxidase (LeGP), heat shock protein 70 (LeHSP70), late embryogenesis abundant (LeLEA), metallothionine type 2 (LeMET2), delta 1-pyrroline-5- carboxylate synthetase (LePCS), ABA-aldehyde oxidase (LeAAO) and 9-cis- Epoxycarotenoid dioxygenase (LeECD) under PEG treatment and cold stress (4 °C). Altogether, these findings suggest that EgDREB1 is a functional regulator in enhancing tolerance to drought and cold stress.
Hevea brasiliensis remains the primary crop commercially exploited to obtain latex, which is produced from the articulated secondary laticifer. Here, we described the transcriptional events related to jasmonic acid (JA)- and linolenic acid (LA)-induced secondary laticifer differentiation (SLD) in H. brasiliensis clone RRIM 600 based on RNA-seq approach. Histochemical approach proved that JA- and LA-treated samples resulted in SLD in H. brasiliensis when compared to ethephon and untreated control. RNA-seq data resulted in 86,614 unigenes, of which 2,664 genes were differentially expressed in JA and LA-induced secondary laticifer harvested from H. brasiliensis bark samples. Among these, 450 genes were unique to JA and LA as they were not differentially expressed in ethephon-treated samples compared with the untreated samples. Most transcription factors from the JA- and LA-specific dataset were classified under MYB, APETALA2/ethylene response factor (AP2/ERF), and basic-helix-loop-helix (bHLH) gene families that were involved in tissue developmental pathways, and we proposed that Bel5-GA2 oxidase 1-KNOTTED-like homeobox complex are likely involved in JA- and LA-induced SLD in H. brasiliensis. We also discovered alternative spliced transcripts, putative novel transcripts, and cis-natural antisense transcript pairs related to SLD event. This study has advanced understanding on the transcriptional regulatory network of SLD in H. brasiliensis.
Atmospheric CO2 level is one of the most important factors which affect plant growth and crop production. Although many crucial genes and pathways have been identified in response to atmospheric CO2 changes, the integrated and precise mechanisms of plant CO2 response are not well understood. Alternative splicing (AS) is an important gene regulation process that affects many biological processes in plants. However, the AS pattern changes in plants in response to elevated CO2 levels have not yet been investigated. Here, we used RNA-Seq data of Arabidopsis thaliana grown under different CO2 concentration to analyze the global changes in AS. We found that AS increased with the rise in CO2 concentration. Additionally, we identified 345 differentially expressed (DE) genes and 251 differentially alternative splicing (DAS) genes under the elevated CO2 condition. Moreover, the results showed that the expression of most of the DAS genes did not change significantly, indicating that AS can serve as an independent mechanism for gene regulation in response to elevated CO2. Furthermore, our analysis of function categories revealed that the DAS genes were associated mainly with the stimulus response. Overall, this the first study to explore the changes of AS in plants in response to elevated CO2.
Morus alba is an important medicinal plant that is used to treat human diseases. The leaf, branch, and root of Morus can be applied as antidiabetic, antioxidant, and anti-inflammatory medicines, respectively. To explore the molecular mechanisms underlying the various pharmacological functions within different parts of Morus, organ-specific proteomics were performed. Protein profiles of the Morus leaf, branch, and root were determined using a gel-free/label-free proteomic technique. In the Morus leaf, branch, and root, a total of 492, 414, and 355 proteins were identified, respectively, including 84 common proteins. In leaf, the main function was related to protein degradation, photosynthesis, and redox ascorbate/glutathione metabolism. In branch, the main function was related to protein synthesis/degradation, stress, and redox ascorbate/glutathione metabolism. In root, the main function was related to protein synthesis/degradation, stress, and cell wall. Additionally, organ-specific metabolites and antioxidant activities were analyzed. These results revealed that flavonoids were highly accumulated in Morus root compared with the branch and leaf. Accordingly, two root-specific proteins named chalcone flavanone isomerase and flavonoid 3,5-hydroxylase were accumulated in the flavonoid pathway. Consistent with this finding, the content of the total flavonoids was higher in root compared to those detected in branch and leaf. These results suggest that the flavonoids in Morus root might be responsible for its biological activity and the root is the main part for flavonoid biosynthesis in Morus.
Ultraviolet-B (UVB) irradiation induces oxidative stress in plant cells due to the generation of excessive reactive oxygen species. Morus alba L. (M. abla) is an important medicinal plant used for the treatment of human diseases. Also, its leaves are widely used as food for silkworms. In our previous research, we found that a high level of UVB irradiation with dark incubation led to the accumulation of secondary metabolites in M. abla leaf. The aim of the present study was to describe and compare M. alba leaf transcriptomics with different treatments (control, UVB, UVB+dark). Leaf transcripts from M. alba were sequenced using an Illumina Hiseq 2000 system, which produced 14.27Gb of data including 153,204,462 paired-end reads among the three libraries. We de novo assembled 133,002 transcripts with an average length of 1270bp and filtered 69,728 non-redundant unigenes. A similarity search was performed against the non-redundant National Center of Biotechnology Information (NCBI) protein database, which returned 41.08% hits. Among the 20,040 unigenes annotated in UniProtKB/SwissProt database, 16,683 unigenes were assigned 102,232 gene ontology terms and 6667 unigenes were identified in 287 known metabolic pathways. Results of differential gene expression analysis together with real-time quantitative PCR tests indicated that UVB irradiation with dark incubation enhanced the flavonoid biosynthesis in M. alba leaf. Our findings provided a valuable proof for a better understanding of the metabolic mechanism under abiotic stresses in M. alba leaf.
We have characterized an oil palm (Elaeis guineensis Jacq.) constitutive promoter that is derived from a translationally control tumor protein (TCTP) gene. The TCTP promoter was fused transcriptionally with the gusA reporter gene and transferred to monocot and dicot systems in order to study its regulatory role in a transient expression study. It was found that the 5' region of TCTP was capable of driving the gusA expression in all the oil palm tissues tested, including immature embryo, embryogenic callus, embryoid, young leaflet from mature palm, green leaf, mesocarp and stem. It could also be used in dicot systems as it was also capable of driving gusA expression in tobacco leaves. The results indicate that the TCTP promoter could be used for the production of recombinant proteins that require constitutive expression in the plant system.
Drought is a major abiotic stress that influences rice production. Although the transcriptomic data of rice against drought is widely available, the regulation of small open reading frames (sORFs) in response to drought stress in rice is yet to be investigated. Different levels of drought stress have different regulatory mechanisms in plants. In this study, drought stress was imposed on four-leaf stage rice, divided into two treatments, 40% and 30% soil moisture content (SMC). The RNAs of the samples were extracted, followed by the RNA sequencing analysis on their sORF expression changes under 40%_SMC and 30%_SMC, and lastly, the expression was validated through NanoString. A total of 122 and 143 sORFs were differentially expressed (DE) in 40%_SMC and 30%_SMC, respectively. In 40%_SMC, 69 sORFs out of 696 (9%) DEGs were found to be upregulated. On the other hand, 69 sORFs out of 449 DEGs (11%) were significantly downregulated. The trend seemed to be higher in 30%_SMC, where 112 (12%) sORFs were found to be upregulated from 928 significantly upregulated DEGs. However, only 8% (31 sORFs out of 385 DEGs) sORFs were downregulated in 30%_SMC. Among the identified sORFs, 110 sORFs with high similarity to rice proteome in the PsORF database were detected in 40%_SMC, while 126 were detected in 30%_SMC. The Gene Ontology (GO) enrichment analysis of DE sORFs revealed their involvement in defense-related biological processes, such as defense response, response to biotic stimulus, and cellular homeostasis, whereas enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways indicated that DE sORFs were associated with tryptophan and phenylalanine metabolisms. Several DE sORFs were identified, including the top five sORFs (OsisORF_3394, OsisORF_0050, OsisORF_3007, OsisORF_6407, and OsisORF_7805), which have yet to be characterised. Since these sORFs were responsive to drought stress, they might hold significant potential as targets for future climate-resilient rice development.
Elucidating the physiological mechanisms of the irregular yet concerted flowering rhythm of mass flowering tree species in the tropics requires long-term monitoring of flowering phenology, exogenous and endogenous environmental factors, as well as identifying interactions and dependencies among these factors. To investigate the proximate factors for floral initiation of mast seeding trees in the tropics, we monitored the expression dynamics of two key flowering genes, meteorological conditions and endogenous resources over two flowering events of Shorea curtisii and Shorea leprosula in the Malay Peninsula. Comparisons of expression dynamics of genes studied indicated functional conservation of FLOWERING LOCUS T (FT) and LEAFY (LFY) in Shorea. The genes were highly expressed at least 1 month before anthesis for both species. A mathematical model considering the synergistic effect of cool temperature and drought on activation of the flowering gene was successful in predicting the observed gene expression patterns. Requirement of both cool temperature and drought for floral transition suggested by the model implies that flowering phenologies of these species are sensitive to climate change. Our molecular phenology approach in the tropics sheds light on the conserved role of flowering genes in plants inhabiting different climate zones and can be widely applied to dissect the flowering processes in other plant species.
Impatiens balsamina L. is a tropical ornamental and traditional medicinal herb rich in natural compounds, especially 2-methoxy-1,4-naphthoquinone (MNQ) which is a bioactive compound with tested anticancer activities. Characterization of key genes involved in the shikimate and 1,4-dihydroxy-2-naphthoate (DHNA) pathways responsible for MNQ biosynthesis and their expression profiles in I. balsamina will facilitate adoption of genetic/metabolic engineering or synthetic biology approaches to further increase production for pre-commercialization. In this study, HPLC analysis showed that MNQ was present in significantly higher quantities in the capsule pericarps throughout three developmental stages (early-, mature- and postbreaker stages) whilst its immediate precursor, 2-hydroxy-1,4-naphthoquinone (lawsone) was mainly detected in mature leaves. Transcriptomes of I. balsamina derived from leaf, flower, and three capsule developmental stages were generated, totalling 59.643 Gb of raw reads that were assembled into 94,659 unigenes (595,828 transcripts). A total of 73.96% of unigenes were functionally annotated against seven public databases and 50,786 differentially expressed genes (DEGs) were identified. Expression profiles of 20 selected genes from four major secondary metabolism pathways were studied and validated using qRT-PCR method. Majority of the DHNA pathway genes were found to be significantly upregulated in early stage capsule compared to flower and leaf, suggesting tissue-specific synthesis of MNQ. Correlation analysis identified 11 candidate unigenes related to three enzymes (NADH-quinone oxidoreductase, UDP-glycosyltransferases and S-adenosylmethionine-dependent O-methyltransferase) important in the final steps of MNQ biosynthesis based on genes expression profiles consistent with MNQ content. This study provides the first molecular insight into the dynamics of MNQ biosynthesis and accumulation across different tissues of I. balsamina and serves as a valuable resource to facilitate further manipulation to increase production of MNQ.
A natural shift is taking place in the approaches being adopted by plant scientists in response to the accessibility of systems-based technology platforms. Metabolomics is one such field, which involves a comprehensive non-biased analysis of metabolites in a given cell at a specific time. This review briefly introduces the emerging field and a range of analytical techniques that are most useful in metabolomics when combined with computational approaches in data analyses. Using cases from Arabidopsis and other selected plant systems, this review highlights how information can be integrated from metabolomics and other functional genomics platforms to obtain a global picture of plant cellular responses. We discuss how metabolomics is enabling large-scale and parallel interrogation of cell states under different stages of development and defined environmental conditions to uncover novel interactions among various pathways. Finally, we discuss selected applications of metabolomics.
Lignification of the plant cell wall could serve as the first line of defense against pathogen attack, but the molecular mechanisms of virulence and disease between oil palm and Ganoderma boninense are poorly understood. This study presents the biochemical, histochemical, enzymology and gene expression evidences of enhanced lignin biosynthesis in young oil palm as a response to G. boninense (GBLS strain). Comparative studies with control (T1), wounded (T2) and infected (T3) oil palm plantlets showed significant accumulation of total lignin content and monolignol derivatives (syringaldehyde and vanillin). These derivatives were deposited on the epidermal cell wall of infected plants. Moreover, substantial differences were detected in the activities of enzyme and relative expressions of genes encoding phenylalanine ammonia lyase (EC 4.3.1.24), cinnamate 4-hydroxylase (EC 1.14.13.11), caffeic acid O-methyltransferase (EC 2.1.1.68) and cinnamyl alcohol dehydrogenase (CAD, EC 1.1.1.195). These enzymes are key intermediates dedicated to the biosynthesis of lignin monomers, the guaicyl (G), syringyl (S) and ρ-hydroxyphenyl (H) subunits. Results confirmed an early, biphasic and transient positive induction of all gene intermediates, except for CAD enzyme activities. These differences were visualized by anatomical and metabolic changes in the profile of lignin in the oil palm plantlets such as low G lignin, indicating a potential mechanism for enhanced susceptibility toward G. boninense infection.
The R2R3-MYB genes comprise one of the largest transcription factor gene families in plants, playing regulatory roles in plant-specific developmental processes, defense responses and metabolite accumulation. To date MYB family genes have not yet been comprehensively identified in the major staple fruit crop banana. In this study, we present a comprehensive, genome-wide analysis of the MYB genes from Musa acuminata DH-Pahang (A genome). A total of 285 R2R3-MYB genes as well as genes encoding three other classes of MYB proteins containing multiple MYB repeats were identified and characterised with respect to structure and chromosomal organisation. Organ- and development-specific expression patterns were determined from RNA-Seq data. For 280 M. acuminata MYB genes for which expression was found in at least one of the analysed samples, a variety of expression patterns were detected. The M. acuminata R2R3-MYB genes were functionally categorised, leading to the identification of seven clades containing only M. acuminata R2R3-MYBs. The encoded proteins may have specialised functions that were acquired or expanded in Musa during genome evolution. This functional classification and expression analysis of the MYB gene family in banana establishes a solid foundation for future comprehensive functional analysis of MaMYBs and can be utilized in banana improvement programmes.
One of the targets in oil palm genetic engineering programme is the production of polyhydroxybutyrate (PHB) and polyhydroxybutyrate-co-valerate (PHBV) in the oil palm leaf tissues. Production of PHB requires the use of phbA (beta-ketothiolase type A), phbB (acetoacetyl-CoA reductase) and phbC (PHB synthase) genes of Ralstonia eutropha, whereas bktB (beta-ketothiolase type B), phbB, phbC genes of R. eutropha and tdcB (threonine dehydratase) gene of Escherichia coli were used for PHBV production. Each of these genes was fused with a transit peptide (Tp) of oil palm acyl-carrier-protein (ACP) gene, driven by an oil palm leaf-specific promoter (LSP1) to genetically engineer the PHB/PHBV pathway to the plastids of the leaf tissues. In total, four transformation vectors, designated pLSP15 (PHB) and pLSP20 (PHBV), and pLSP13 (PHB) and pLSP23 (PHBV), were constructed for transformation in Arabidopsis thaliana and oil palm, respectively. The phosphinothricin acetyltransferase gene (bar) driven by CaMV35S promoter in pLSP15 and pLSP20, and ubiquitin promoter in pLSP13 and pLSP23 were used as the plant selectable markers. Matrix attachment region of tobacco (RB7MAR) was also included in the vectors to stabilize the transgene expression and to minimize silencing due to positional effect. Restriction digestion, PCR amplification and/or sequencing were carried out to ensure sequence integrity and orientation.
The somatic embryogenesis tissue culture process has been utilized to propagate high yielding oil palm. Due to the low callogenesis and embryogenesis rates, molecular studies were initiated to identify genes regulating the process, and their expression levels are usually quantified using reverse transcription quantitative real-time PCR (RT-qPCR). With the recent release of oil palm genome sequences, it is crucial to establish a proper strategy for gene analysis using RT-qPCR. Selection of the most suitable reference genes should be performed for accurate quantification of gene expression levels.
KEY MESSAGE: Transcript profiling during the early induction phase of oil palm tissue culture and RNAi studies in a model somatic embryogenesis system showed that EgENOD93 expression is essential for somatic embryogenesis. Micropropagation of oil palm through tissue culture is vital for the generation of superior and uniform elite planting materials. Studies were carried out to identify genes to distinguish between leaf explants with the potential to develop into embryogenic or non-embryogenic callus. Oil palm cDNA microarrays were co-hybridized with cDNA probes of reference tissue, separately with embryo forming (media T527) and non-embryo (media T694) forming leaf explants sampled at Day 7, Day 14 and Day 21. Analysis of the normalized datasets has identified 77, 115 and 127 significantly differentially expressed genes at Day 7, Day 14, and Day 21, respectively. An early nodulin 93 protein gene (ENOD93), was highly expressed at Day 7, Day 14, and Day 21 and in callus (media T527), as assessed by RT-qPCR. Validation of EgENOD93 across tissue culture lines of different genetic background and media composition showed the potential of this gene as an embryogenic marker. In situ RNA hybridization and functional characterization in Medicago truncatula provided additional evidence that ENOD93 is essential for somatic embryogenesis. This study supports the suitability of EgENOD93 as a marker to predict the potential of leaf explants to produce embryogenic callus. Crosstalk among stresses, auxin, and Nod-factor like signalling molecules likely induces the expression of EgENOD93 for embryogenic callus formation.