METHODS: Phylogenetic analysis was used to identify candidate Dendrobium catenatum R2R3-MYB (DcaMYB) sequences associated with pigment and cell shape development. Gene silencing of candidate DhMYBs in Dendrobium hybrid by direct application of dsRNA to developing flowers was followed by observation of gene expression level and flower phenotypes. Silencing of the structural gene chalcone synthase was used as a comparative control.
KEY RESULTS: Ten candidate flower-associated DcaMYBs were identified. Flowers treated with dsRNA of DhMYB22 and DhMYB60 sequences were less pigmented and had relatively low expression of anthocyanin biosynthetic genes (F3'H and DFR), lower total anthocyanin concentration and markedly lower levels of cyanidin-3-glucoside and cyanidin-3-rutinoside. Petals of DhMYB22-treated flowers and sepals of DhMYB60-treated flowers showed the greatest colour difference relative to the same organs in untreated flowers. DhMYB22-treated flowers had relatively narrow and constricted lips, while DhMYB60-treated flowers had narrow and constricted sepals. No significant difference in shape was observed for DhCHS-treated or untreated flowers.
CONCLUSIONS: Our results demonstrate that DhMYB22 and DhMYB60 regulate pigment intensity and floral organ shape in Dendrobium. This is a first report of MYB regulation of floral organ shape in orchids.
RESULTS: Using total RNA extracted from young roots as template, we performed a comparative transcriptome analysis on oil palm responding to 14d and 28d of Pi deprivation treatment and under adequate Pi supply. By using Illumina HiSeq4000 platform, RNA-Seq analysis was successfully conducted on 12 paired-end RNA-Seq libraries and generated more than 1.2 billion of clean reads in total. Transcript abundance estimated by fragments per kilobase per million fragments (FPKM) and differential expression analysis revealed 36 and 252 genes that are differentially regulated in Pi-starved roots at 14d and 28d, respectively. Genes possibly involved in regulating Pi homeostasis, nutrient uptake and transport, hormonal signaling and gene transcription were found among the differentially expressed genes.
CONCLUSIONS: Our results showed that the molecular response mechanism underlying Pi starvation in oil palm is complexed and involved multilevel regulation of various sensing and signaling components. This contribution would generate valuable genomic resources in the effort to develop oil palm planting materials that possess Pi-use efficient trait through molecular manipulation and breeding programs.
RESULT: SPME GC-MS analysis showed the highest terpenoid accumulation on the 6th day post-inoculation (dpi) compared to the other treatment time points (0 dpi, 3 dpi, and 9 dpi). Among the increased terpenoid compounds, α-cedrene, valencene and β-bisabolene were prominent. P. minor inoculated for 6 days was selected for miRNA library construction using next generation sequencing. Differential gene expression analysis showed that 58 miRNAs belonging to 30 families had significantly altered regulation.
Among these 58 differentially expressed genes (DEGs), 27 [corrected] miRNAs were upregulated, whereas 31 [corrected] miRNAs were downregulated. Two putative novel pre-miRNAs were identified and validated through reverse transcriptase PCR. Prediction of target transcripts potentially involved in the mevalonate pathway (MVA) was carried out by psRobot software, resulting in four miRNAs: pmi-miR530, pmi-miR6173, pmi-miR6300 and a novel miRNA, pmi-Nov_13. In addition, two miRNAs, miR396a and miR398f/g, were predicted to have their target transcripts in the non-mevalonate pathway (MEP). In addition, a novel miRNA, pmi-Nov_12, was identified to have a target gene involved in green leaf volatile (GLV) biosynthesis. RT-qPCR analysis showed that pmi-miR6173, pmi-miR6300 and pmi-nov_13 were downregulated, while miR396a and miR398f/g were upregulated. Pmi-miR530 showed upregulation at 9 dpi, and dynamic expression was observed for pmi-nov_12. Pmi-6300 and pmi-miR396a cleavage sites were detected through degradome sequence analysis. Furthermore, the relationship between miRNA metabolites and mRNA metabolites was validated using correlation analysis.
CONCLUSION: Our findings suggest that six studied miRNAs post-transcriptionally regulate terpenoid biosynthesis in P. minor. This regulatory behaviour of miRNAs has potential as a genetic tool to regulate terpenoid biosynthesis in P. minor.
RESULTS: The study generated 846,762 high quality sequence reads, with an average length of 334 bp and totalling 283 Mbp. De novo assembly generated 36,384 and 35,269 unigene sequences for M. acuminata Calcutta 4 and Cavendish Grande Naine, respectively. A total of 64.4% of the unigenes were annotated through Basic Local Alignment Search Tool (BLAST) similarity analyses against public databases.Assembled sequences were functionally mapped to Gene Ontology (GO) terms, with unigene functions covering a diverse range of molecular functions, biological processes and cellular components. Genes from a number of defense-related pathways were observed in transcripts from each cDNA library. Over 99% of contig unigenes mapped to exon regions in the reference M. acuminata DH Pahang whole genome sequence. A total of 4068 genic-SSR loci were identified in Calcutta 4 and 4095 in Cavendish Grande Naine. A subset of 95 potential defense-related gene-derived simple sequence repeat (SSR) loci were validated for specific amplification and polymorphism across M. acuminata accessions. Fourteen loci were polymorphic, with alleles per polymorphic locus ranging from 3 to 8 and polymorphism information content ranging from 0.34 to 0.82.
CONCLUSIONS: A large set of unigenes were characterized in this study for both M. acuminata Calcutta 4 and Cavendish Grande Naine, increasing the number of public domain Musa ESTs. This transcriptome is an invaluable resource for furthering our understanding of biological processes elicited during biotic stresses in Musa. Gene-based markers will facilitate molecular breeding strategies, forming the basis of genetic linkage mapping and analysis of quantitative trait loci.
RESULTS: The present report is a continuing study from our previous published transcriptomic profiling of oil palm seedlings against G. boninense. We focused on identifying differentially expressed genes (DEGs) encoding transcription factors (TFs) from the same RNA-seq data; resulting in 106 upregulated and 108 downregulated TFs being identified. The DEGs are involved in four established defense-related pathways responsible for cell wall modification, reactive oxygen species (ROS)-mediated signaling, programmed cell death (PCD) and plant innate immunity. We discovered upregulation of JUNGBRUNNEN 1 (EgJUB1) during the fungal biotrophic phase while Ethylene Responsive Factor 113 (EgERF113) demonstrated prominent upregulation when the palm switches to defense against necrotrophic phase. EgJUB1 was shown to have a binding activity to a 19 bp palindromic SNBE1 element, WNNYBTNNNNNNNAMGNHW found in the promoter region of co-expressing EgHSFC-2b. Further in silico analysis of promoter regions revealed co-expression of EgJUB1 with TFs containing SNBE1 element with single nucleotide change at either the 5th or 18th position. Meanwhile, EgERF113 binds to both GCC and DRE/CRT elements promoting plasticity in upregulating the downstream defense-related genes. Both TFs were proven to be nuclear-localized based on subcellular localization experiment using onion epidermal cells.
CONCLUSION: Our findings demonstrated unprecedented transcriptional reprogramming of specific TFs potentially to enable regulation of a specific set of genes during different infection phases of this hemibiotrophic fungal pathogen. The results propose the intricacy of oil palm defense response in orchestrating EgJUB1 during biotrophic and EgERF113 during the subsequent transition to the necrotrophic phase. Binding of EgJUB1 to SNBE motif instead of NACBS while EgERF113 to GCC-box and DRE/CRT motifs is unconventional and not normally associated with pathogen infection. Identification of these phase-specific oil palm TFs is important in designing strategies to tackle or attenuate the progress of infection.
RESULTS: A total of 12 standard cDNA libraries, representing three main developmental stages in oil palm tissue culture, were generated in this study. Random sequencing of clones from these cDNA libraries generated 17,599 expressed sequence tags (ESTs). The ESTs were analysed, annotated and assembled to generate 9,584 putative unigenes distributed in 3,268 consensi and 6,316 singletons. These unigenes were assigned putative functions based on similarity and gene ontology annotations. Cluster analysis, which surveyed the relatedness of each library based on the abundance of ESTs in each consensus, revealed that lipid transfer proteins were highly expressed in embryogenic tissues. A glutathione S-transferase was found to be highly expressed in non-embryogenic callus. Further analysis of the unigenes identified 648 non-redundant simple sequence repeats and 211 putative full-length open reading frames.
CONCLUSION: This study has provided an overview of genes expressed during oil palm tissue culture. Candidate genes with expression that are modulated during tissue culture were identified. However, in order to confirm whether these genes are suitable as early markers for embryogenesis, the genes need to be tested on earlier stages of tissue culture and a wider range of genotypes. This collection of ESTs is an important resource for genetic and genome analyses of the oil palm, particularly during tissue culture development.
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.