Search

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., We compared genic and intergenic regions of three Dioscorea chloroplast genomes and seven pairs of primers had 80%–100% PCR success rate. 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 majority species based on more than 100 Dioscorea individuals from 30 species.

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.

Coconut palm has two distinct types—“tall” and “dwarf”—which differ morphologically. Tall coconut varieties need 8–10 years to start flowering, while dwarf coconut varieties only require 3–5 years. We compared seedling and reproductive stage transcriptomes for both coconut types to determine potential molecular mechanisms underlying control of flowering time in coconut. Several key genes in the photoperiod pathway were differentially expressed between seedling and reproductive leaf samples in both tall and dwarf coconut. These genes included suppressor of overexpression of constans (SOC1), flowering locus T (FT), and Apetala 1 (AP1). Alternative splicing analysis of genes in the photoperiod pathway further revealed that the FT gene produces different transcripts in tall compared to dwarf coconut. The shorter alternative splice variant of FT [which included a 6 bp deletion, alternative 3′ splicing sites (A3SS)] was found to be exclusively present in dwarf coconut varieties but absent in most tall coconut varieties. Our results provide a valuable information resource as well as suggesting a probable mechanism for differentiation of flowering time onset in coconut, providing a target for future breeding work in accelerating time to flowering in this crop species.

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.