Comparative Transcriptome Analysis Provides Insights into the Effect of Epicuticular Wax Accumulation on Salt Stress in Coconuts.

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.

Developing functional markers for vitamin E biosynthesis in oil palm.

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.

Genome-wide discovery and characterization of long noncoding RNAs in African oil palm (Elaeis guineensis Jacq.).

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.

Alternative splicing of flowering time gene FT is associated with halving of time to flowering in coconut.

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.

Identifying Vitamin E Biosynthesis Genes in Elaeis guineensis by Genome-Wide Association Study.

Elaeis guineensis is a tropical oil crop and has the highest oil yield per unit area. Palm oil has high palmitic acid content and is also rich in vitamins, including vitamin E. We conducted genome-wide association studies in a diversity panel of 161 E. guineensis accessions to identify single-nucleotide polymorphisms (SNPs) linked with vitamin E and validated candidate genes in these marker-associated intervals. Based on the SNPs reported in our previous research, 47 SNP markers were detected to be significantly associated with the variation of tocopherol and tocotrienol content at a cutoff P value of 6.3 × 10-7. A total of 656 candidate genes in the flanking regions of the 47 SNPs were identified, followed by pathway enrichment analysis. Of these candidate genes, EgHGGT (homogentisate geranylgeranyl transferase) involved in the biosynthesis of tocotrienols had a higher expression level in the mesocarp compared to other tissues. Expression of the EgHGGT gene was positively correlated with the variation in α-tocotrienol content. Induced overexpression of the gene in Arabidopsis caused a significant increase in vitamin E content and production of α-tocotrienols compared to wild Arabidopsis.

Identification and Validation of Candidate Genes Involved in Fatty Acid Content in Oil Palm by Genome-Wide Association Analysis.

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.

Development of high-resolution DNA barcodes for Dioscorea species discrimination and phylogenetic analysis.

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.

Development of High-Density SNP Markers and Their Application in Evaluating Genetic Diversity and Population Structure in Elaeis guineensis.

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.

Genome-wide identification of WRKY genes and their expression profiles under different abiotic stresses in Elaeis guineensis.

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.

[Biodiversity of Rhizobia associated with Acacia melanoxylon grown in South China].

OBJECTIVE: This study aimed to investigate the genetic diversity of 174 isolates of symbiotic bacteria associated with Acacia melanoxylon obtained from 15 sampling sites in Guangdong, Fujian and Jiangxi provinces of China. METHODS: The 16S rDNA restriction fragment length polymorphism (RFLP) and phylogenetic analyses of the 16S rDNA and housekeeping genes (recA, glnII and atpD). RESULTS: In the 16S rDNA PCR-RFLP analysis, 9 rDNA types were identified among the 174 isolates; Phylogenetic analyses based on 16S rDNA and housekeeping gene sequences indicated that 34 representative isolates belonged to the genus Bradyrhizobium, Rhizobium, Mesorhizobium in Alpha-Proteobacteria, and the most closely related strains are Bradyrhizobium liaoningense, Bradyrhizobium betae, Bradyrhizobium cytisi, Rhizobium multihospitium and Mesorhizobium plurifarium. CONCLUSION: All of the isolates could be identified to general, and Bradyrhizobium, Rhizobium or Mesorhizobium could be the dominant microsymbiont. The microsymbionts associated with Acacia melanoxylon showed relative genetic diversity.