Search

The coconut is an important tropical economical crop and exhibits high tolerance to various types of salinity stress. However, little is known about the molecular mechanism underlying its salt tolerance. In this study, RNA-Seq was applied to examine the different genes expressed in four coconut varieties when exposed to a salt environment, resulting in the generation of data for 48 transcriptomes. Comparative transcriptome analysis showed that some genes involved in cutin and wax biosynthesis were significantly upregulated in salt treatment compared to the control, including CYP86A4, HTH, CER1, CER2, CER3, DCR, GPAT4, LTP3, LTP4, and LTP5. In particular, the expression of CER2 was induced more than sixfold, with an RPKM value of up to 205 ten days after salt treatment in Hainan Tall coconut, demonstrating superior capacity in salt tolerance compared to dwarf coconut varieties. However, for yellow dwarf and red dwarf coconut varieties, the expression level of the CER2 gene was low at four different time points after exposure to salt treatment, suggesting that this gene may contribute to the divergence in salt tolerance between tall and dwarf coconut varieties. Cytological evidence showed a higher abundance of cuticle accumulation in tall coconut and severe damage to cuticular wax in dwarf coconut.

Abstract The genus Dioscorea is widely distributed in tropical and subtropical regions, and is economically important in terms of food supply and pharmaceutical applications. However, DNA barcodes are relatively unsuccessful in discriminating between Dioscorea species, with the highest discrimination rate (23.26%) derived from matK sequences. In this study, we compared genic and intergenic regions of three Dioscorea chloroplast genomes and found that the density of SNPs and indels in intergenic sites was about twice and seven times higher than that of SNPs and indels in the genic regions, respectively. A total of 52 primer pairs covering highly variable regions were designed and seven pairs of primers had 80%–100% PCR success rate. PCR amplicons of 73 Dioscorea individuals and assembled sequences of 47 Dioscorea SRAs were used for estimating intraspecific and interspecific divergence for the seven loci: The rpoB‐trnC locus had the highest interspecific divergence. Automatic barcoding gap discovery (ABGD), Poisson tree processes (PTP), and generalized mixed Yule coalescence (GMYC) analysis were applied for species delimitation based on the seven loci and successfully identified the majority of species, except for species in the Enantiophyllum section. Phylogenetic analysis of 51 Dioscorea individuals (28 species) showed that most individuals belonging to the same species tended to cluster in the same group. Our results suggest that the variable loci derived from comparative analysis of plastid genome sequences could be good DNA barcode candidates for taxonomic analysis and species delimitation.

Vitamin E is essential for human health and plays positive roles in anti-oxidation. Previously, we detected large variation in vitamin E content among 161 oil palm accessions. In this study, twenty oil palm accessions with distinct variation in vitamin E contents (171.30 to 1 258.50 ppm) were selected for genetic variation analysis and developing functional markers associated with vitamin E contents. Thirty-seven homologous genes in oil palm belonging to vitamin E biosynthesis pathway were identified via BLASTP analysis, the lengths of which ranged from 426 to 25 717 bp (average 7 089 bp). Multiplex PCR sequencing for the 37 genes found 1 703 SNPs and 85 indels among the 20 oil palm accessions, with 226 SNPs locating in the coding regions. Clustering analysis for these polymorphic loci showed that the 20 oil palm accessions could be divided into five groups. Among these groups, group I included eight oil palm accessions whose vitamin E content (mean value: 893.50 ppm) was far higher than other groups (mean value 256.29 to 532.94 ppm). Correlation analysis between the markers and vitamin E traits showed that 134 SNP and 7 indel markers were significantly (p < 0.05) related with total vitamin E content. Among these functional markers, the indel EgTMT-1-24 was highly correlated with variation in vitamin E content, especially tocotrienol content. Our study identified a number of candidate function associated markers and provided clues for further research into molecular breeding for high vitamin E content oil palm.

Long noncoding RNAs (lncRNAs) are an important class of genes and play important roles in a range of biological processes. However, few reports have described the identification of lncRNAs in oil palm. In this study, we applied strand specific RNA-seq with rRNA removal to identify 1,363 lncRNAs from the equally mixed tissues of oil palm spear leaf and six different developmental stages of mesocarp (8–24 weeks). Based on strand specific RNA-seq data and 18 released oil palm transcriptomes, we systematically characterized the expression patterns of lncRNA loci and their target genes. A total of 875 uniq target genes for natural antisense lncRNAs (NAT-lncRNA, 712), long intergenic noncoding RNAs (lincRNAs, 92), intronic-lncRNAs (33), and sense-lncRNAs (52) were predicted. A majority of lncRNA loci (77.8%–89.6%) had low expression in 18 transcriptomes, while only 89 lncRNA loci had medium to high expression in at least one transcriptome. Coexpression analysis between lncRNAs and their target genes indicated that 6% of lncRNAs had expression patterns positively correlated with those of target genes. Based on single nucleotide polymorphism (SNP) markers derived from our previous research, 6,882 SNPs were detected for lncRNAs and 28 SNPs belonging to 21 lncRNAs were associated with the variation of fatty acid contents. Moreover, seven lncRNAs showed expression patterns positively correlated expression pattern with those of genes in de novo fatty acid synthesis pathways. Our study identified a collection of lncRNAs for oil palm and provided clues for further research into lncRNAs that may regulate mesocarp development and lipid metabolism.

Oil palm (Elaeis guineensis) is the highest yielding oil crop per unit area worldwide, but its oil is considered unhealthy for human consumption due to its high palmitic acid content (C16:0). In order to facilitate breeding for fatty acid content in oil palm, genome-wide association analysis (GWAS) was used to identify and validate single-nucleotide polymorphism (SNP) markers and underlying candidate genes associated with fatty acid content in a diversity panel of 200 oil palm individuals. A total of 1,261,501 SNP markers previously developed using SLAF-seq (specific locus amplified fragment sequencing) were used for GWAS. Based on this analysis, 62 SNP markers were significantly associated with fatty acid composition, and 223 candidate genes were identified in the flanking regions of these SNPs. We found one gene (acyl-ACP thioesterase B genes) that was involved in fatty acid biosynthesis and that was associated with high palmitic acid content in the mesocarp. Over-expression of this gene caused a significant increase in palmitic acid content. Our study provides key loci that can be used for breeding oil palm cultivars with low palmitic acid content.

High-density single nucleotide polymorphisms (SNPs) are used as highly favored makers to analyze genetic diversity and population structure, to construct high-density genetic maps and provide genotypes for genome-wide association analysis. In order to develop genome-wide SNP markers in oil palm (Elaeis guineensis), single locus amplified fragment sequencing (SLAF-seq) technology was performed in a diversity panel of 200 oil palm individuals and 1,261,501 SNPs were identified with minor allele frequency > 0.05 and integrity > 1. Among them, only 17.81% can be mapped within the genic region and the remaining was located into the intergenic region. A positive correlation was detected between the distribution of SNP markers and retrotransposons [transposable elements (TEs)]. Population structure analysis showed that the 200 individuals of oil palm can be divided into five subgroups based on cross-validation errors. However, the subpopulations divided for the 200 oil palm individuals based on the SNP markers were not accurately related to their geographical origins and 80 oil palm individuals from Malaysia showed highest genetic diversity. In addition, the physical distance of linkage disequilibrium (LD) decay in the analyzed oil palm population was 14.516 kb when r2 = 0.1. The LD decay distances for different chromosomes varied from 3.324 (chromosome 15) to 19.983 kb (chromosome 7). Our research provides genome-wide SNPs for future targeted breeding in palm oil.

African oil palm (Elaeis guineensis) is an important oil crop grown in tropical region and sensitive to low temperature along with high tolerance to salt and drought stresses. Since the WRKY transcription factor family plays central roles in the regulation of plant stress tolerance, 95 genes belonging to the WRKY family were identified and characterized in oil palm genome. Gene structure analysis showed that EgWRKY genes have considerable variation in intron number (0 to 12) and gene length (477bp to 89,167 bp). Duplicated genes identification indicated 32 EgWRKY genes originated from segmental duplication and two from tandem duplication. Based on transcriptome data, most EgWRKY genes showed tissue-specific expression patterns and their expression could be induced under cold stress. Furthermore, six EgWRKY genes with more than two-folded increased expression level under cold stress were validated by RT-qPCR, which has higher expression level in cold, drought and high salinity treatment. The identification and characterization of WRKY gene family showed that EgWRKY were associated with a wide range of abiotic stress responses in Elaeis guineensis and some EgWRKY members with high expression levels could be selected for further research in analyzing their functions in the stress response in African oil palm.