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 3412 (2001 annotated) gene candidates were found to be significantly differentially expressed between high- and low-yielding palms at at least one of the different stages of mesocarp development evaluated. Gene Ontologies (GO) enrichment analysis identified 28 significantly enriched GO terms, including regulation of transcription, fatty acid biosynthesis and metabolic processes. These differentially expressed genes comprise several transcription factors, such as, bHLH, Dof zinc finger proteins and MADS box proteins. Several genes involved in glycolysis, TCA, and fatty acid biosynthesis pathways were also found up-regulated in high-yielding oil palm, among them; pyruvate dehydrogenase E1 component Subunit Beta (PDH), ATP-citrate lyase, β- ketoacyl-ACP synthases I (KAS I), β- ketoacyl-ACP synthases III (KAS III) and ketoacyl-ACP reductase (KAR). Sucrose metabolism-related genes such as Invertase, Sucrose Synthase 2 and Sucrose Phosphatase 2 were found to be down-regulated in high-yielding oil palms, compared to the lower yield palms.
CONCLUSIONS: Our findings indicate that a higher carbon flux (channeled through down-regulation of the Sucrose Synthase 2 pathway) was being utilized by up-regulated genes involved in glycolysis, TCA and fatty acid biosynthesis leading to enhanced oil production in the high-yielding oil palm. These findings are an important stepping stone to understand the processes that lead to production of high-yielding oil palms and have implications for breeding to maximize oil production.
RESULTS: Using two independent gene-prediction pipelines, Fgenesh++ and Seqping, 26,059 oil palm genes with transcriptome and RefSeq support were identified from the oil palm genome. These coding regions of the genome have a characteristic broad distribution of GC3 (fraction of cytosine and guanine in the third position of a codon) with over half the GC3-rich genes (GC3 ≥ 0.75286) being intronless. In comparison, only one-seventh of the oil palm genes identified are intronless. Using comparative genomics analysis, characterization of conserved domains and active sites, and expression analysis, 42 key genes involved in FA biosynthesis in oil palm were identified. For three of them, namely EgFABF, EgFABH and EgFAD3, segmental duplication events were detected. Our analysis also identified 210 candidate resistance genes in six classes, grouped by their protein domain structures.
CONCLUSIONS: We present an accurate and comprehensive annotation of the oil palm genome, focusing on analysis of important categories of genes (GC3-rich and intronless), as well as those associated with important functions, such as FA biosynthesis and disease resistance. The study demonstrated the advantages of having an integrated approach to gene prediction and developed a computational framework for combining multiple genome annotations. These results, available in the oil palm annotation database ( http://palmxplore.mpob.gov.my ), will provide important resources for studies on the genomes of oil palm and related crops.
REVIEWERS: This article was reviewed by Alexander Kel, Igor Rogozin, and Vladimir A. Kuznetsov.