Oil palm genome sequence reveals divergence of interfertile species in Old and New worlds

Genetic diversity of Elaeis oleifera (HBK) Cortes populations using cross species SSRs: implication's for germplasm utilization and conservation

Background

The Elaeis oleifera genetic materials were assembled from its center of diversity in South and Central America. These materials are currently being preserved in Malaysia as ex situ living collections. Maintaining such collections is expensive and requires sizable land. Information on the genetic diversity of these collections can help achieve efficient conservation via maintenance of core collection. For this purpose, we have applied fourteen unlinked microsatellite markers to evaluate 532 E. oleifera palms representing 19 populations distributed across Honduras, Costa Rica, Panama and Colombia.

Results

In general, the genetic diversity decreased from Costa Rica towards the north (Honduras) and south-east (Colombia). Principle coordinate analysis (PCoA) showed a single cluster indicating low divergence among palms. The phylogenetic tree and STRUCTURE analysis revealed clusters based on country of origin, indicating considerable gene flow among populations within countries. Based on the values of the genetic diversity parameters, some genetically diverse populations could be identified. Further, a total of 34 individual palms that collectively captured maximum allelic diversity with reduced redundancy were also identified. High pairwise genetic differentiation (Fst > 0.250) among populations was evident, particularly between the Colombian populations and those from Honduras, Panama and Costa Rica. Crossing selected palms from highly differentiated populations could generate off-springs that retain more genetic diversity.

Conclusion

The results attained are useful for selecting palms and populations for core collection. The selected materials can also be included into crossing scheme to generate offsprings that capture greater genetic diversity for selection gain in the future.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0505-7) contains supplementary material, which is available to authorized users.

Transcriptome and functional analysis reveals hybrid vigor for oil biosynthesis in oil palm

Oil palm is the most productive oil crop in the world and composes 36% of the world production. However, the molecular mechanisms of hybrids vigor (or heterosis) between Dura, Pisifera and their hybrid progeny Tenera has not yet been well understood. Here we compared the temporal and spatial compositions of lipids and transcriptomes for two oil yielding organs mesocarp and endosperm from Dura, Pisifera and Tenera. Multiple lipid biosynthesis pathways are highly enriched in all non-additive expression pattern in endosperm, while cytokinine biosynthesis and cell cycle pathways are highly enriched both in endosperm and mesocarp. Compared with parental palms, the high oil content in Tenera was associated with much higher transcript levels of EgWRI1, homolog of Arabidopsis thaliana WRINKLED1. Among 338 identified genes in lipid synthesis, 207 (61%) has been identified to contain the WRI1 specific binding AW motif. We further functionally identified EgWRI1-1, one of three EgWRI1 orthologs, by genetic complementation of the Arabidopsis wri1 mutant. Ectopic expression of EgWRI1-1 in plant produced dramatically increased seed mass and oil content, with oil profile changed. Our findings provide an explanation for EgWRI1 as an important gene contributing hybrid vigor in lipid biosynthesis in oil palm.

Improving transcriptome de novo assembly by using a reference genome of a related species: Translational genomics from oil palm to coconut

The palms are a family of tropical origin and one of the main constituents of the ecosystems of these regions around the world. The two main species of palm represent different challenges: coconut (Cocos nucifera L.) is a source of multiple goods and services in tropical communities, while oil palm (Elaeis guineensis Jacq) is the main protagonist of the oil market. In this study, we present a workflow that exploits the comparative genomics between a target species (coconut) and a reference species (oil palm) to improve the transcriptomic data, providing a proteome useful to answer functional or evolutionary questions. This workflow reduces redundancy and fragmentation, two inherent problems of transcriptomic data, while preserving the functional representation of the target species. Our approach was validated in Arabidopsis thaliana using Arabidopsis lyrata and Capsella rubella as references species. This analysis showed the high sensitivity and specificity of our strategy, relatively independent of the reference proteome. The workflow increased the length of proteins products in A. thaliana by 13%, allowing, often, to recover 100% of the protein sequence length. In addition redundancy was reduced by a factor greater than 3. In coconut, the approach generated 29,366 proteins, 1,246 of these proteins deriving from new contigs obtained with the BRANCH software. The coconut proteome presented a functional profile similar to that observed in rice and an important number of metabolic pathways related to secondary metabolism. The new sequences found with BRANCH software were enriched in functions related to biotic stress. Our strategy can be used as a complementary step to de novo transcriptome assembly to get a representative proteome of a target species. The results of the current analysis are available on the website PalmComparomics (http://palm-comparomics.southgreen.fr/).

Plant miRNAs found in human circulating system provide evidences of cross kingdom RNAi

Background

Emerging evidence indicates that plant miRNAs can present within human circulating system through dietary intake and regulate human gene expression. Hence we deduced that comestible plants miRNAs can be identified in the public available small RNA sequencing data sets.

Results

In this study, we identified abundant plant miRNAs sequences from 410 human plasma small RNA sequencing data sets. One particular plant miRNA miR2910, conserved in fruits and vegetables, was found to present in high relative amount in the plasma samples. This miRNA, with same 6mer and 7mer-A1 target seed sequences as hsa-miR-4259 and hsa-miR-4715-5p, was predicted to target human JAK-STAT signaling pathway gene SPRY4 and transcription regulation genes.

Conclusions

Through analysis of public available plasma small RNA sequencing data, we found the supporting evidence for the plant miRNAs cross kingdom RNAi within human circulating system.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3502-3) contains supplementary material, which is available to authorized users.

Development, identification and validation of CAPS marker for SHELL trait which governs dura, pisifera and tenera fruit forms in oil palm (Elaeis guineensis Jacq.)

The oil palm fruit forms (dura, pisifera and tenera) governed by the shell thickness gene (Sh) plays a major role in identification of fruit type and also influences palm oil yield. Identification of desired fruit type is a major asset to the breeders and oil palm workers for applications in breeding, seed certification and to reduce time, space and money spent on identification of fruit form. In the present study, we developed Sh gene specific primer pairs and bulk segregant analysis was done using 300 genomic and 8 genic SSR markers. We identified one cleaved amplified polymorphic site (CAPS) marker for differentiation of oil palm fruit type which produced two alleles (280 and 250bp) in dura genotypes, three alleles in tenera genotypes (550, 280, and 250bp) and one allele in pisifera genotypes (550bp). The shell allele sequencing results showed that two SNPs were present, of which SNP2 contributed for variation of fruit forms. The nucleotide ‘A’ was present in only dura genotypes, where as ‘T’ was present only in pisifera genotypes, which in turn led to the change of amino acid lysine to aspargine. The identified CAPS marker was validated on 300 dura, 25 pisifera and 80 tenera genotypes, 80 dura/ pisifera cross progenies and 60 lines of tenera/ tenera cross progeny. Association mapping of marker data with phenotypic data of eight oil yield related traits resulted in identification of seven significant QTLs by GLM approach, four by MLM approach at a significant threshold (P) level of 0.001. Significant QTLs were identified for fruit to bunch and oil to bunch traits, which explained R² of 12.9% and 11.5% respectively. The CAPS marker used in the present study facilitate selection and timely distribution of desirable high yielding tenera sprouts to the farmers instead of waiting for 4–5 years. This saves a lot of land, time and money which will be a major breakthrough to the oil palm community.

Field Guide to Plant Model Systems

For the past several decades, advances in plant development, physiology, cell biology, and genetics have relied heavily on the model (or reference) plant Arabidopsis thaliana. Arabidopsis resembles other plants, including crop plants, in many but by no means all respects. Study of Arabidopsis alone provides little information on the evolutionary history of plants, evolutionary differences between species, plants that survive in different environments, or plants that access nutrients and photosynthesize differently. Empowered by the availability of large-scale sequencing and new technologies for investigating gene function, many new plant models are being proposed and studied.

Lineage-Specific Evolutionary Histories and Regulation of Major Starch Metabolism Genes during Banana Ripening

Starch is the most widespread and abundant storage carbohydrate in plants. It is also a major feature of cultivated bananas as it accumulates to large amounts during banana fruit development before almost complete conversion to soluble sugars during ripening. Little is known about the structure of major gene families involved in banana starch metabolism and their evolution compared to other species. To identify genes involved in banana starch metabolism and investigate their evolutionary history, we analyzed six gene families playing a crucial role in plant starch biosynthesis and degradation: the ADP-glucose pyrophosphorylases (AGPases), starch synthases (SS), starch branching enzymes (SBE), debranching enzymes (DBE), α-amylases (AMY) and β-amylases (BAM). Using comparative genomics and phylogenetic approaches, these genes were classified into families and sub-families and orthology relationships with functional genes in Eudicots and in grasses were identified. In addition to known ancestral duplications shaping starch metabolism gene families, independent evolution in banana and grasses also occurred through lineage-specific whole genome duplications for specific sub-families of AGPase, SS, SBE, and BAM genes; and through gene-scale duplications for AMY genes. In particular, banana lineage duplications yielded a set of AGPase, SBE and BAM genes that were highly or specifically expressed in banana fruits. Gene expression analysis highlighted a complex transcriptional reprogramming of starch metabolism genes during ripening of banana fruits. A differential regulation of expression between banana gene duplicates was identified for SBE and BAM genes, suggesting that part of starch metabolism regulation in the fruit evolved in the banana lineage.

Application of Genomic Technologies to the Breeding of Trees

The recent introduction of next generation sequencing (NGS) technologies represents a major revolution in providing new tools for identifying the genes and/or genomic intervals controlling important traits for selection in breeding programs. In perennial fruit trees with long generation times and large sizes of adult plants, the impact of these techniques is even more important. High-throughput DNA sequencing technologies have provided complete annotated sequences in many important tree species. Most of the high-throughput genotyping platforms described are being used for studies of genetic diversity and population structure. Dissection of complex traits became possible through the availability of genome sequences along with phenotypic variation data, which allow to elucidate the causative genetic differences that give rise to observed phenotypic variation. Association mapping facilitates the association between genetic markers and phenotype in unstructured and complex populations, identifying molecular markers for assisted selection and breeding. Also, genomic data provide in silico identification and characterization of genes and gene families related to important traits, enabling new tools for molecular marker assisted selection in tree breeding. Deep sequencing of transcriptomes is also a powerful tool for the analysis of precise expression levels of each gene in a sample. It consists in quantifying short cDNA reads, obtained by NGS technologies, in order to compare the entire transcriptomes between genotypes and environmental conditions. The miRNAs are non-coding short RNAs involved in the regulation of different physiological processes, which can be identified by high-throughput sequencing of RNA libraries obtained by reverse transcription of purified short RNAs, and by in silico comparison with known miRNAs from other species. All together, NGS techniques and their applications have increased the resources for plant breeding in tree species, closing the former gap of genetic tools between trees and annual species.

Gene coexpression network analysis of oil biosynthesis in an interspecific backcross of oil palm

Global demand for vegetable oils is increasing at a dramatic rate, while our understanding of the regulation of oil biosynthesis in plants remains limited. To gain insights into the mechanisms that govern oil synthesis and fatty acid (FA) composition in the oil palm fruit, we used a multilevel approach combining gene coexpression analysis, quantification of allele-specific expression and joint multivariate analysis of transcriptomic and lipid data, in an interspecific backcross population between the African oil palm, Elaeis guineensis, and the American oil palm, Elaeis oleifera, which display contrasting oil contents and FA compositions. The gene coexpression network produced revealed tight transcriptional coordination of fatty acid synthesis (FAS) in the plastid with sugar sensing, plastidial glycolysis, transient starch storage and carbon recapture pathways. It also revealed a concerted regulation, along with FAS, of both the transfer of nascent FA to the endoplasmic reticulum, where triacylglycerol assembly occurs, and of the production of glycerol-3-phosphate, which provides the backbone of triacylglycerols. Plastid biogenesis and auxin transport were the two other biological processes most tightly connected to FAS in the network. In addition to WRINKLED1, a transcription factor (TF) known to activate FAS genes, two novel TFs, termed NF-YB-1 and ZFP-1, were found at the core of the FAS module. The saturated FA content of palm oil appeared to vary above all in relation to the level of transcripts of the gene coding for β-ketoacyl-acyl carrier protein synthase II. Our findings should facilitate the development of breeding and engineering strategies in this and other oil crops.

Development and Validation of a High-Density SNP Genotyping Array for African Oil Palm

High-density single nucleotide polymorphism (SNP) genotyping arrays are powerful tools that can measure the level of genetic polymorphism within a population. To develop a whole-genome SNP array for oil palms, SNP discovery was performed using deep resequencing of eight libraries derived from 132 Elaeis guineensis and Elaeis oleifera palms belonging to 59 origins, resulting in the discovery of >3 million putative SNPs. After SNP filtering, the Illumina OP200K custom array was built with 170 860 successful probes. Phenetic clustering analysis revealed that the array could distinguish between palms of different origins in a way consistent with pedigree records. Genome-wide linkage disequilibrium declined more slowly for the commercial populations (ranging from 120 kb at r(2) = 0.43 to 146 kb at r(2) = 0.50) when compared with the semi-wild populations (19.5 kb at r(2) = 0.22). Genetic fixation mapping comparing the semi-wild and commercial population identified 321 selective sweeps. A genome-wide association study (GWAS) detected a significant peak on chromosome 2 associated with the polygenic component of the shell thickness trait (based on the trait shell-to-fruit; S/F %) in tenera palms. Testing of a genomic selection model on the same trait resulted in good prediction accuracy (r = 0.65) with 42% of the S/F % variation explained. The first high-density SNP genotyping array for oil palm has been developed and shown to be robust for use in genetic studies and with potential for developing early trait prediction to shorten the oil palm breeding cycle.

Draft genome sequence of an elite Dura palm and whole-genome patterns of DNA variation in oil palm

Oil palm is the world’s leading source of vegetable oil and fat. Dura, Pisifera and Tenera are three forms of oil palm. The genome sequence of Pisifera is available whereas the Dura form has not been sequenced yet. We sequenced the genome of one elite Dura palm, and re-sequenced 17 palm genomes. The assemble genome sequence of the elite Dura tree contained 10,971 scaffolds and was 1.701 Gb in length, covering 94.49% of the oil palm genome. 36,105 genes were predicted. Re-sequencing of 17 additional palm trees identified 18.1 million SNPs. We found high genetic variation among palms from different geographical regions, but lower variation among Southeast Asian Dura and Pisifera palms. We mapped 10,000 SNPs on the linkage map of oil palm. In addition, high linkage disequilibrium (LD) was detected in the oil palms used in breeding populations of Southeast Asia, suggesting that LD mapping is likely to be practical in this important oil crop. Our data provide a valuable resource for accelerating genetic improvement and studying the mechanism underlying phenotypic variations of important oil palm traits.

Identification of a Δ12 fatty acid desaturase from oil palm (Elaeis guineensis Jacq.) involved in the biosynthesis of linoleic acid by heterologous expression in Saccharomyces cerevisiae

Oil palm (Elaeis guineensis Jacq.) is one of the highest oil-yield crops in the world. A Δ12-desaturases associated with the primary steps of long-chain polyunsaturated fatty acid (LC-PUFA) biosynthesis were successfully cloned from oil palm and their functions identified. The open reading frames (ORFs) of egFAD2 (GenBank accession: KT023602) consisted of 1176bp and code for 391 amino acids. Their deduced polypeptides showed 75-93% identity to microsomal Δ12-desaturases from other higher plants, and each contained the three histidine clusters typical of the catalytic domains of such enzymes. RT-PCR experiment indicated that the egFAD2 gene exhibited the highest accumulation in the mesocarp of fruits at 120-140 DAP (i.e. the fourth period of fruit development) and, despite having different expression levels, the other four stages were at significantly lower levels compared with the fourth stage. Plasmid pYES2-egFAD2 was transformed into Saccharomyces cerevisiae strain INVSc1 using lithium acetate method for expression under the induction of galactose. Yeast cells transformed with plasmid constructs containing egFAD12 produced an appreciable amount of linoleic acids (18:2(Δ9,)(12)), not normally present in wild-type yeast cells, indicating that the genes encoded functional Δ12-desaturase enzymes.

Computational identification and comparative analysis of miRNA precursors in three palm species

In the present study, miRNA precursors in the genomes of three palm species were identified. Analyzes of sequence conservation and biological function of their putative targets contribute to understand the roles of miRNA in palm biology. MicroRNAs are small RNAs of 20-25 nucleotides in length, with important functions in the regulation of gene expression. Recent genome sequencing of the palm species Elaeis guineensis, Elaeis oleifera and Phoenix dactylifera have enabled the discovery of miRNA genes, which can be used as biotechnological tools in palm trees breeding. The goal of this study is the identification of miRNA precursors in the genomes of these species and their possible biological roles suggested by the mature miRNA-based regulation of target genes. Mature miRNA sequences from Arabidopsis thaliana, Oryza sativa, and Zea mays available at the miRBase were used to predict microRNA precursors in the palm genomes. Three hundred and thirty-eight precursors, ranging from 76 to 220 nucleotide (nt) in size and distributed in 33 families were identified. Moreover, we also identified 266 miRNA precursors of Musa acuminata, which are phylogenetically close to palms species. To understand the biological function of palm miRNAs, 374 putative miRNA targets were identified. An enrichment analysis of target-gene function was carried out using the agriGO tool. The results showed that the targets are involved in plant developmental processes, mainly regulating root development. Our findings contribute to increase the knowledge on microRNA roles in palm biology and could help breeding programs of palm trees.

Fine-mapping and cross-validation of QTLs linked to fatty acid composition in multiple independent interspecific crosses of oil palm

Background

The commercial oil palm (Elaeis guineensis Jacq.) produces a mesocarp oil (commonly called ‘palm oil’) with approximately equal proportions of saturated and unsaturated fatty acids (FAs). An increase in unsaturated FAs content or iodine value (IV) as a measure of the degree of unsaturation would help to open up new markets for the oil. One way to manipulate the fatty acid composition (FAC) in palm oil is through introgression of favourable alleles from the American oil palm, E. oleifera, which has a more unsaturated oil.

Results

In this study, a segregating E. oleifera x E. guineensis (OxG) hybrid population for FAC is used to identify quantitative trait loci (QTLs) linked to IV and various FAs. QTL analysis revealed 10 major and two putative QTLs for IV and six FAs, C14:0, C16:0, C16:1, C18:0, C18:1 and C18:2 distributed across six linkage groups (LGs), OT1, T2, T3, OT4, OT6 and T9. The major QTLs for IV and C16:0 on LGOT1 explained 60.0 – 69.0 % of the phenotypic trait variation and were validated in two independent BC₂ populations. The genomic interval contains several key structural genes in the FA and oil biosynthesis pathways such as PATE/FATB, HIBCH, BASS2, LACS4 and DGAT1 and also a relevant transcription factor (TF), WRI1. The literature suggests that some of these genes can exhibit pleiotropic effects in the regulatory networks of these traits. Using the whole genome sequence data, markers tightly linked to the candidate genes were also developed. Clustering trait values according to the allelic forms of these candidate markers revealed significant differences in the IV and FAs of the palms in the mapping and validation crosses.

Conclusions

The candidate gene approach described and exploited here is useful to identify the potential causal genes linked to FAC and can be adopted for marker-assisted selection (MAS) in oil palm.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2607-4) contains supplementary material, which is available to authorized users.

Evolution of plant genome architecture

We have witnessed an explosion in our understanding of the evolution and structure of plant genomes in recent years. Here, we highlight three important emergent realizations: (1) that the evolutionary history of all plant genomes contains multiple, cyclical episodes of whole-genome doubling that were followed by myriad fractionation processes; (2) that the vast majority of the variation in genome size reflects the dynamics of proliferation and loss of lineage-specific transposable elements; and (3) that various classes of small RNAs help shape genomic architecture and function. We illustrate ways in which understanding these organism-level and molecular genetic processes can be used for crop plant improvement.

Sequencing Crop Genomes: A Gateway to Improve Tropical Agriculture

Agricultural development in the tropics lags behind development in the temperate latitudes due to the lack of advanced technology, and various biotic and abiotic factors. To cope with the increasing demand for food and other plant-based products, improved crop varieties have to be developed. To breed improved varieties, a better understanding of crop genetics is necessary. With the advent of next-generation DNA sequencing technologies, many important crop genomes have been sequenced. Primary importance has been given to food crops, including cereals, tuber crops, vegetables, and fruits. The DNA sequence information is extremely valuable for identifying key genes controlling important agronomic traits and for identifying genetic variability among the cultivars. However, massive DNA re-sequencing and gene expression studies have to be performed to substantially improve our understanding of crop genetics. Application of the knowledge obtained from the genomes, transcriptomes, expression studies, and epigenetic studies would enable the development of improved varieties and may lead to a second green revolution. The applications of next generation DNA sequencing technologies in crop improvement, its limitations, future prospects, and the features of important crop genome projects are reviewed herein.

Genome-wide association study identifies three key loci for high mesocarp oil content in perennial crop oil palm

GWAS in out-crossing perennial crops is typically limited by insufficient marker density to account for population diversity and effects of population structure resulting in high false positive rates. The perennial crop oil palm is the most productive oil crop. We performed GWAS for oil-to-dry-mesocarp content (O/DM) on 2,045 genotyped tenera palms using 200K SNPs that were selected based on the short-range linkage disequilibrium distance, which is inherent with long breeding cycles and heterogeneous breeding populations. Eighty loci were significantly associated with O/DM (p ≤ 10⁻⁴) and three key signals were found. We then evaluated the progeny of a Deli x AVROS breeding trial and a 4% higher O/DM was observed amongst those having the beneficial genotypes at two of the three key loci (p < 0.05). We have initiated MAS and large-scale planting of elite dura and pisifera parents to generate the new commercial tenera palms with higher O/DM potential.

Genome-Wide Identification and Transferability of Microsatellite Markers between Palmae Species

The Palmae family contains 202 genera and approximately 2800 species. Except for Elaeis guineensis and Phoenix dactylifera, almost no genetic and genomic information is available for Palmae species. Therefore, this is an obstacle to the conservation and genetic assessment of Palmae species, especially those that are currently endangered. The study was performed to develop a large number of microsatellite markers which can be used for genetic analysis in different Palmae species. Based on the assembled genome of E. guineensis and P. dactylifera, a total of 814 383 and 371 629 microsatellites were identified. Among these microsatellites identified in E. guineensis, 734 509 primer pairs could be designed from the flanking sequences of these microsatellites. The majority (618 762) of these designed primer pairs had in silico products in the genome of E. guineensis. These 618 762 primer pairs were subsequently used to in silico amplify the genome of P. dactylifera. A total of 7 265 conserved microsatellites were identified between E. guineensis and P. dactylifera. One hundred and thirty-five primer pairs flanking the conserved SSRs were stochastically selected and validated to have high cross-genera transferability, varying from 16.7 to 93.3% with an average of 73.7%. These genome-wide conserved microsatellite markers will provide a useful tool for genetic assessment and conservation of different Palmae species in the future.

Non-tenera Contamination and the Economic Impact of SHELL Genetic Testing in the Malaysian Independent Oil Palm Industry

Oil palm (Elaeis guineensis) is the most productive oil bearing crop worldwide. It has three fruit forms, namely dura (thick-shelled), pisifera (shell-less) and tenera (thin-shelled), which are controlled by the SHELL gene. The fruit forms exhibit monogenic co-dominant inheritance, where tenera is a hybrid obtained by crossing maternal dura and paternal pisifera palms. Commercial palm oil production is based on planting thin-shelled tenera palms, which typically yield 30% more oil than dura palms, while pisifera palms are female-sterile and have little to no palm oil yield. It is clear that tenera hybrids produce more oil than either parent due to single gene heterosis. The unintentional planting of dura or pisifera palms reduces overall yield and impacts land utilization that would otherwise be devoted to more productive tenera palms. Here, we identify three additional novel mutant alleles of the SHELL gene, which encode a type II MADS-box transcription factor, and determine oil yield via control of shell fruit form phenotype in a manner similar to two previously identified mutant SHELL alleles. Assays encompassing all five mutations account for all dura and pisifera palms analyzed. By assaying for these variants in 10,224 mature palms or seedlings, we report the first large scale accurate genotype-based determination of the fruit forms in independent oil palm planting sites and in the nurseries that supply them throughout Malaysia. The measured non-tenera contamination rate (10.9% overall on a weighted average basis) underscores the importance of SHELL genetic testing of seedlings prior to planting in production fields. By eliminating non-tenera contamination, comprehensive SHELL genetic testing can improve sustainability by increasing yield on existing planted lands. In addition, economic modeling demonstrates that SHELL gene testing will confer substantial annual economic gains to the oil palm industry, to Malaysian gross national income and to Malaysian government tax receipts.

Natural (13) C distribution in oil palm (Elaeis guineensis Jacq.) and consequences for allocation pattern

Oil palm has now become one of the most important crops, palm oil representing nearly 25% of global plant oil consumption. Many studies have thus addressed oil palm ecophysiology and photosynthesis-based models of carbon allocation have been used. However, there is a lack of experimental data on carbon fixation and redistribution within palm trees, and important C-sinks have not been fully characterized yet. Here, we carried out extensive measurement of natural (13) C-abundance (δ(13) C) in oil palm tissues, including fruits at different maturation stages. We find a (13) C-enrichment in heterotrophic organs compared to mature leaves, with roots being the most (13) C-enriched. The δ(13) C in fruits decreased during maturation, reflecting the accumulation in (13) C-depleted lipids. We further used observed δ(13) C values to compute plausible carbon fluxes using a steady-state model of (13) C-distribution including metabolic isotope effects ((12) v/(13) v). The results suggest that fruits represent a major respiratory loss (≈39% of total tree respiration) and that sink organs such as fruits are fed by sucrose from leaves. That is, glucose appears to be a quantitatively important compound in palm tissues, but computations indicate that it is involved in dynamic starch metabolism rather that C-exchange between organs.

Effects of methylation-sensitive enzymes on the enrichment of genic SNPs and the degree of genome complexity reduction in a two-enzyme genotyping-by-sequencing (GBS) approach: a case study in oil palm (Elaeis guineensis)

Advances in next generation sequencing have facilitated a large-scale single nucleotide polymorphism (SNP) discovery in many crop species. Genotyping-by-sequencing (GBS) approach couples next generation sequencing with genome complexity reduction techniques to simultaneously identify and genotype SNPs. Choice of enzymes used in GBS library preparation depends on several factors including the number of markers required, the desired level of multiplexing, and whether the enrichment of genic SNP is preferred. We evaluated various combinations of methylation-sensitive (AatII, PstI, MspI) and methylation-insensitive (SphI, MseI) enzymes for their effectiveness in genome complexity reduction and enrichment of genic SNPs. We discovered that the use of two methylation-sensitive enzymes effectively reduced genome complexity and did not require a size selection step. On the contrary, the genome coverage of libraries constructed with methylation-insensitive enzymes was quite high, and the additional size selection step may be required to increase the overall read depth. We also demonstrated the effectiveness of methylation-sensitive enzymes in enriching for SNPs located in genic regions. When two methylation-insensitive enzymes were used, only 16% of SNPs identified were located in genes and 18% in the vicinity (± 5 kb) of the genic regions, while most SNPs resided in the intergenic regions. In contrast, a remarkable degree of enrichment was observed when two methylation-sensitive enzymes were employed. Almost two thirds of the SNPs were located either inside (32-36%) or in the vicinity (28-31%) of the genic regions. These results provide useful information to help researchers choose appropriate GBS enzymes in oil palm and other crop species.

De Novo Assembly of Coding Sequences of the Mangrove Palm (Nypa fruticans) Using RNA-Seq and Discovery of Whole-Genome Duplications in the Ancestor of Palms

Nypa fruticans (Arecaceae) is the only monocot species of true mangroves. This species represents the earliest mangrove fossil recorded. How N. fruticans adapts to the harsh and unstable intertidal zone is an interesting question. However, the 60 gene segments deposited in NCBI are insufficient for solving this question. In this study, we sequenced, assembled and annotated the transcriptome of N. fruticans using next-generation sequencing technology. A total of 19,918,800 clean paired-end reads were de novo assembled into 45,368 unigenes with a N50 length of 1,096 bp. A total of 41.35% unigenes were functionally annotated using Blast2GO. Many genes annotated to "response to stress" and 15 putative positively selected genes were identified. Simple sequence repeats were identified and compared with other palms. The divergence time between N. fruticans and other palms was estimated at 75 million years ago using the genomic data, which is consistent with the fossil record. After calculating the synonymous substitution rate between paralogs, we found that two whole-genome duplication events were shared by N. fruticans and other palms. These duplication events provided a large amount of raw material for the more than 2,000 later speciation events in Arecaceae. This study provides a high quality resource for further functional and evolutionary studies of N. fruticans and palms in general.

Authentication of medicinal plants by DNA markers

Medicinal plants have been used worldwide for centuries to maintain health and to treat diseases, more so chronic diseases. However, adulteration and use of spurious materials as substitutes have become a major concern for users and industry for reasons of safety and efficacy. Therefore, authentication of medicinal plants is of utmost importance. Morphological, anatomical, chemical and DNA markers solve the problem by differentiating the genuine material from the adulterants, substitutes and spurious drugs. DNA markers use nucleotide sequences to identify species; it takes preference over the other two markers being not age dependent, tissue specific and having a higher discriminating power. Therefore, characterization of plants with such markers is an ideal approach for identification of medicinal plant species and populations/varieties of the same species. Availability of certified taxonomic specimens in herbaria is certainly required for unambiguous confirmation through final visual comparison and analysis.

The pineapple genome and the evolution of CAM photosynthesis

Pineapple (Ananas comosus (L.) Merr.) is the most economically valuable crop possessing crassulacean acid metabolism (CAM), a photosynthetic carbon assimilation pathway with high water-use efficiency, and the second most important tropical fruit. We sequenced the genomes of pineapple varieties F153 and MD2 and a wild pineapple relative, Ananas bracteatus accession CB5. The pineapple genome has one fewer ancient whole-genome duplication event than sequenced grass genomes and a conserved karyotype with seven chromosomes from before the ρ duplication event. The pineapple lineage has transitioned from C3 photosynthesis to CAM, with CAM-related genes exhibiting a diel expression pattern in photosynthetic tissues. CAM pathway genes were enriched with cis-regulatory elements associated with the regulation of circadian clock genes, providing the first cis-regulatory link between CAM and circadian clock regulation. Pineapple CAM photosynthesis evolved by the reconfiguration of pathways in C3 plants, through the regulatory neofunctionalization of preexisting genes and not through the acquisition of neofunctionalized genes via whole-genome or tandem gene duplication.

Non-canonical structure, function and phylogeny of the Bsister MADS-box gene OsMADS30 of rice (Oryza sativa)

Bsister MADS-box genes play key roles in female reproductive organ and seed development throughout seed plants. This view is supported by their high conservation in terms of sequence, expression and function. In grasses, there are three subclades of Bsister genes: the OsMADS29-, the OsMADS30- and the OsMADS31-like genes. Here, we report on the evolution of the OsMADS30-like genes. Our analyses indicate that these genes evolved under relaxed purifying selection and are rather weakly expressed. OsMADS30, the representative of the OsMADS30-like genes from rice (Oryza sativa), shows strong sequence deviations in its 3' region when compared to orthologues from other grass species. We show that this is due to a 2.4-kbp insertion, possibly of a hitherto unknown helitron, which confers a heterologous C-terminal domain to OsMADS30. This putative helitron is not present in the OsMADS30 orthologues from closely related wild rice species, pointing to a relatively recent insertion event. Unlike other Bsister mutants O. sativa plants carrying a T-DNA insertion in the OsMADS30 gene do not show aberrant seed phenotypes, indicating that OsMADS30 likely does not have a canonical 'Bsister function'. However, imaging-based phenotyping of the T-DNA carrying plants revealed alterations in shoot size and architecture. We hypothesize that sequence deviations that accumulated during a period of relaxed selection in the gene lineage that led to OsMADS30 and the alteration of the C-terminal domain might have been a precondition for a potential neo-functionalization of OsMADS30 in O. sativa.

Genome-wide analysis of LTR-retrotransposons in oil palm

BACKGROUND

The oil palm (Elaeis guineensis Jacq.) is a major cultivated crop and the world's largest source of edible vegetable oil. The genus Elaeis comprises two species E. guineensis, the commercial African oil palm and E. oleifera, which is used in oil palm genetic breeding. The recent publication of both the African oil palm genome assembly and the first draft sequence of its Latin American relative now allows us to tackle the challenge of understanding the genome composition, structure and evolution of these palm genomes through the annotation of their repeated sequences.

METHODS

In this study, we identified, annotated and compared Transposable Elements (TE) from the African and Latin American oil palms. In a first step, Transposable Element databases were built through de novo detection in both genome sequences then the TE content of both genomes was estimated. Then putative full-length retrotransposons with Long Terminal Repeats (LTRs) were further identified in the E. guineensis genome for characterization of their structural diversity, copy number and chromosomal distribution. Finally, their relative expression in several tissues was determined through in silico analysis of publicly available transcriptome data.

RESULTS

Our results reveal a congruence in the transpositional history of LTR retrotransposons between E. oleifera and E. guineensis, especially the Sto-4 family. Also, we have identified and described 583 full-length LTR-retrotransposons in the Elaeis guineensis genome. Our work shows that these elements are most likely no longer mobile and that no recent insertion event has occurred. Moreover, the analysis of chromosomal distribution suggests a preferential insertion of Copia elements in gene-rich regions, whereas Gypsy elements appear to be evenly distributed throughout the genome.

CONCLUSIONS

Considering the high proportion of LTR retrotransposon in the oil palm genome, our work will contribute to a greater understanding of their impact on genome organization and evolution. Moreover, the knowledge gained from this study constitutes a valuable resource for both the improvement of genome annotation and the investigation of the evolutionary history of palms.

Ploidy and domestication are associated with genome size variation in Palms

PREMISE OF THE STUDY

The genome size of a species (C-value) is associated with growth, development and adaptation to environmental changes. Angiosperm C-values range 1200-fold and frequently vary within species, although little is known about the impacts of domestication on genome size. Genome size variation among related species of palms is of evolutionary significance because changes characterize clades and may be associated with polyploidy, transposon amplifications, deletions, or rearrangements. Further knowledge of genome size will provide crucial information needed for planning of whole genome sequencing and accurate annotations. We studied the genome size of Cocos nucifera and its variation among cultivars, and compared it to values for related palms from the Attaleinae subtribe.

METHODS

Flow cytometric analysis of isolated nuclei from young palm leaves was used to estimate genome sizes of 23 coconut cultivars (Talls, Dwarfs, and hybrids) worldwide and 17 Cocoseae species. Ancestral genome size was reconstructed on a maximum likelihood phylogeny of Attaleinae from seven WRKY loci.

KEY RESULTS

The coconut genome is large-averaging 5.966 pg-and shows intraspecific variation associated with domestication. Variation among Tall coconuts was significantly greater than among Dwarfs. Attaleinae genomes showed moderate size variation across genera, except polyploids Jubaeopsis caffra, Voanioala gerardii, Beccariophoenix alfredii, and Allagoptera caudescens, which had larger genomes.

CONCLUSIONS

Our results contribute to the understanding of the relationship between domestication and genome size in long-lived tree crops and provide a basis for whole-genome sequencing of the coconut and other domesticated plants. Polyploidy evolved independently in two clades within Attaleinae.

Loss of Karma transposon methylation underlies the mantled somaclonal variant of oil palm

Somaclonal variation arises in plants and animals when differentiated somatic cells are induced into a pluripotent state, but the resulting clones differ from each other and from their parents. In agriculture, somaclonal variation has hindered the micropropagation of elite hybrids and genetically modified crops, but the mechanism responsible remains unknown. The oil palm fruit 'mantled' abnormality is a somaclonal variant arising from tissue culture that drastically reduces yield, and has largely halted efforts to clone elite hybrids for oil production. Widely regarded as an epigenetic phenomenon, 'mantling' has defied explanation, but here we identify the MANTLED locus using epigenome-wide association studies of the African oil palm Elaeis guineensis. DNA hypomethylation of a LINE retrotransposon related to rice Karma, in the intron of the homeotic gene DEFICIENS, is common to all mantled clones and is associated with alternative splicing and premature termination. Dense methylation near the Karma splice site (termed the Good Karma epiallele) predicts normal fruit set, whereas hypomethylation (the Bad Karma epiallele) predicts homeotic transformation, parthenocarpy and marked loss of yield. Loss of Karma methylation and of small RNA in tissue culture contributes to the origin of mantled, while restoration in spontaneous revertants accounts for non-Mendelian inheritance. The ability to predict and cull mantling at the plantlet stage will facilitate the introduction of higher performing clones and optimize environmentally sensitive land resources.

Molecular defense response of oil palm to Ganoderma infection

Basal stem rot (BSR) of oil palm roots is due to the invasion of fungal mycelia of Ganoderma species which spreads to the bole of the stem. In addition to root contact, BSR can also spread by airborne basidiospores. These fungi are able to break down cell wall components including lignin. BSR not only decreases oil yield, it also causes the stands to collapse thus causing severe economic loss to the oil palm industry. The transmission and mode of action of Ganoderma, its interactions with oil palm as a hemibiotroph, and the molecular defence responses of oil palm to the infection of Ganoderma boninense in BSR are reviewed, based on the transcript profiles of infected oil palms. The knowledge gaps that need to be filled in oil palm-Ganoderma molecular interactions i.e. the associations of hypersensitive reaction (HR)-induced cell death and reactive oxygen species (ROS) kinetics to the susceptibility of oil palm to Ganoderma spp., the interactions of phytohormones (salicylate, jasmonate and ethylene) at early and late stages of BSR, and cell wall strengthening through increased production of guaiacyl (G)-type lignin, are also discussed.

The Genome of Dendrobium officinale Illuminates the Biology of the Important Traditional Chinese Orchid Herb

Dendrobium officinale Kimura et Migo is a traditional Chinese orchid herb that has both ornamental value and a broad range of therapeutic effects. Here, we report the first de novo assembled 1.35 Gb genome sequences for D. officinale by combining the second-generation Illumina Hiseq 2000 and third-generation PacBio sequencing technologies. We found that orchids have a complete inflorescence gene set and have some specific inflorescence genes. We observed gene expansion in gene families related to fungus symbiosis and drought resistance. We analyzed biosynthesis pathways of medicinal components of D. officinale and found extensive duplication of SPS and SuSy genes, which are related to polysaccharide generation, and that the pathway of D. officinale alkaloid synthesis could be extended to generate 16-epivellosimine. The D. officinale genome assembly demonstrates a new approach to deciphering large complex genomes and, as an important orchid species and a traditional Chinese medicine, the D. officinale genome will facilitate future research on the evolution of orchid plants, as well as the study of medicinal components and potential genetic breeding of the dendrobe.

Molecular basis of angiosperm tree architecture

The architecture of trees greatly impacts the productivity of orchards and forestry plantations. Amassing greater knowledge on the molecular genetics that underlie tree form can benefit these industries, as well as contribute to basic knowledge of plant developmental biology. This review describes the fundamental components of branch architecture, a prominent aspect of tree structure, as well as genetic and hormonal influences inferred from studies in model plant systems and from trees with non-standard architectures. The bulk of the molecular and genetic data described here is from studies of fruit trees and poplar, as these species have been the primary subjects of investigation in this field of science.

Biotechnology of oil palm: strategies towards manipulation of lipid content and composition

Oil palm is a major economic crop for Malaysia. The major challenges faced by the industry are labor shortage, availability of arable land and unstable commodity price. This has caused the industry to diversify its applications into higher value products besides increasing its yield. While conventional breeding has its limitations, biotechnology was identified as one of the tools for overcoming the above challenges. Research on biotechnology of oil palm began more than two decades ago leveraging a multidisciplinary approach involving biochemical studies, gene and promoter isolation, transformation vector construction and finally genetic transformation to produce the targeted products. The main target of oil palm biotechnology research is to increase oleic acid in the mesocarp. Other targets are stearic acid, palmitoleic acid, ricinoleic acid, lycopene (carotenoid) and biodegradable plastics. Significant achievements were reported for the biochemical studies, isolation of useful oil palm genes and characterization of important promoters. A large number of transformation constructs for various targeted products were successfully produced using the isolated oil palm genes and promoters. Finally transformation of these constructs into oil palm embryogenic calli was carried out while the regeneration of transgenic oil palm harboring the useful genes is in progress.

Genomic selection prediction accuracy in a perennial crop: case study of oil palm (Elaeis guineensis Jacq.)

Genomic selection empirically appeared valuable for reciprocal recurrent selection in oil palm as it could account for family effects and Mendelian sampling terms, despite small populations and low marker density. Genomic selection (GS) can increase the genetic gain in plants. In perennial crops, this is expected mainly through shortened breeding cycles and increased selection intensity, which requires sufficient GS accuracy in selection candidates, despite often small training populations. Our objective was to obtain the first empirical estimate of GS accuracy in oil palm (Elaeis guineensis), the major world oil crop. We used two parental populations involved in conventional reciprocal recurrent selection (Deli and Group B) with 131 individuals each, genotyped with 265 SSR. We estimated within-population GS accuracies when predicting breeding values of non-progeny-tested individuals for eight yield traits. We used three methods to sample training sets and five statistical methods to estimate genomic breeding values. The results showed that GS could account for family effects and Mendelian sampling terms in Group B but only for family effects in Deli. Presumably, this difference between populations originated from their contrasting breeding history. The GS accuracy ranged from -0.41 to 0.94 and was positively correlated with the relationship between training and test sets. Training sets optimized with the so-called CDmean criterion gave the highest accuracies, ranging from 0.49 (pulp to fruit ratio in Group B) to 0.94 (fruit weight in Group B). The statistical methods did not affect the accuracy. Finally, Group B could be preselected for progeny tests by applying GS to key yield traits, therefore increasing the selection intensity. Our results should be valuable for breeding programs with small populations, long breeding cycles, or reduced effective size.

Discovery of novel genes derived from transposable elements using integrative genomic analysis

Complex eukaryotes contain millions of transposable elements (TEs), comprising large fractions of their nuclear genomes. TEs consist of structural, regulatory, and coding sequences that are ordinarily associated with transposition, but that occasionally confer on the organism a selective advantage and may thereby become exapted. Exapted transposable element genes (ETEs) are known to play critical roles in diverse systems, from vertebrate adaptive immunity to plant development. Yet despite their evident importance, most ETEs have been identified fortuitously and few systematic searches have been conducted, suggesting that additional ETEs may await discovery. To explore this possibility, we develop a comprehensive systematic approach to searching for ETEs. We use TE-specific conserved domains to identify with high precision genes derived from TEs and screen them for signatures of exaptation based on their similarities to reference sets of known ETEs, conventional (non-TE) genes, and TE genes across diverse genetic attributes including repetitiveness, conservation of genomic location and sequence, and levels of expression and repressive small RNAs. Applying this approach in the model plant Arabidopsis thaliana, we discover a surprisingly large number of novel high confidence ETEs. Intriguingly, unlike known plant ETEs, several of the novel ETE families form tandemly arrayed gene clusters, whereas others are relatively young. Our results not only identify novel TE-derived genes that may have practical applications but also challenge the notion that TE exaptation is merely a relic of ancient life, instead suggesting that it may continue to fundamentally drive evolution.

Genome-wide SNP discovery and identification of QTL associated with agronomic traits in oil palm using genotyping-by-sequencing (GBS)

Oil palm has become one of the most important oil crops in the world. Marker-assisted selections have played a pivotal role in oil palm breeding programs. Here, we report the use of genotyping-by-sequencing (GBS) approach for a large-scale SNP discovery and genotyping of a mapping population. Reduced representation libraries of 108 F2 progeny were sequenced and a total of 524 million reads were obtained. We detected 21,471 single nucleotide substitutions, most of which (62.6%) represented transition events. Of 3417 fully informative SNP markers, we were able to place 1085 on a linkage map, which spanned 1429.6 cM and had an average of one marker every 1.26 cM. Three QTL affecting trunk height were detected on LG 10, 14 and 15, whereas a single QTL associated with fruit bunch weight was identified on LG 3. The use of GBS approach proved to be rapid, cost-effective and highly reproducible in this species.

A consensus linkage map of oil palm and a major QTL for stem height

Oil palm (Elaeis guinensis Jacquin) is the most important source of vegetable oil and fat. Several linkage maps had been constructed using dominant and co-dominant markers to facilitate mapping of QTL. However, dominant markers are not easily transferable among different laboratories. We constructed a consensus linkage map for oil palm using co-dominant markers (i.e. microsatellite and SNPs) and two F1 breeding populations generated by crossing Dura and Pisifera individuals. Four hundreds and forty-four microsatellites and 36 SNPs were mapped onto 16 linkage groups. The map length was 1565.6 cM, with an average marker space of 3.72 cM. A genome-wide scan of QTL identified a major QTL for stem height on the linkage group 5, which explained 51% of the phenotypic variation. Genes in the QTL were predicted using the palm genome sequence and bioinformatic tools. The linkage map supplies a base for mapping QTL for accelerating the genetic improvement, and will be also useful in the improvement of the assembly of the genome sequences. Markers linked to the QTL may be used in selecting dwarf trees. Genes within the QTL will be characterized to understand the mechanisms underlying dwarfing.

Oil palm natural diversity and the potential for yield improvement

African oil palm has the highest productivity amongst cultivated oleaginous crops. Species can constitute a single crop capable to fulfill the growing global demand for vegetable oils, which is estimated to reach 240 million tons by 2050. Two types of vegetable oil are extracted from the palm fruit on commercial scale. The crude palm oil and kernel palm oil have different fatty acid profiles, which increases versatility of the crop in industrial applications. Plantations of the current varieties have economic life-span around 25-30 years and produce fruits around the year. Thus, predictable annual palm oil supply enables marketing plans and adjustments in line with the economic forecasts. Oil palm cultivation is one of the most profitable land uses in the humid tropics. Oil palm fruits are the richest plant source of pro-vitamin A and vitamin E. Hence, crop both alleviates poverty, and could provide a simple practical solution to eliminate global pro-vitamin A deficiency. Oil palm is a perennial, evergreen tree adapted to cultivation in biodiversity rich equatorial land areas. The growing demand for the palm oil threatens the future of the rain forests and has a large negative impact on biodiversity. Plant science faces three major challenges to make oil palm the key element of building the future sustainable world. The global average yield of 3.5 tons of oil per hectare (t) should be raised to the full yield potential estimated at 11-18t. The tree architecture must be changed to lower labor intensity and improve mechanization of the harvest. Oil composition should be tailored to the evolving needs of the food, oleochemical and fuel industries. The release of the oil palm reference genome sequence in 2013 was the key step toward this goal. The molecular bases of agronomically important traits can be and are beginning to be understood at the single base pair resolution, enabling gene-centered breeding and engineering of this remarkable crop.

Evolutionary divergence of β-expansin structure and function in grasses parallels emergence of distinctive primary cell wall traits

Expansins are wall-loosening proteins that promote the extension of primary cell walls without the hydrolysis of major structural components. Previously, proteins from the EXPA (α-expansin) family were found to loosen eudicot cell walls but to be less effective on grass cell walls, whereas the reverse pattern was found for EXPB (β-expansin) proteins obtained from grass pollen. To understand the evolutionary and structural bases for the selectivity of EXPB action, we assessed the extension (creep) response of cell walls from diverse monocot families to EXPA and EXPB treatments. Cell walls from Cyperaceae and Juncaceae (families closely related to grasses) displayed a typical grass response ('β-response'). Walls from more distant monocots, including some species that share with grasses high levels of arabinoxylan, responded preferentially to α-expansins ('α-response'), behaving in this regard like eudicots. An expansin with selective activity for grass cell walls was detected in Cyperaceae pollen, coinciding with the expression of genes from the divergent EXPB-I branch that includes grass pollen β-expansins. The evolutionary origin of this branch was located within Poales on the basis of phylogenetic analyses and its association with the 'sigma' whole-genome duplication. Accelerated evolution in this branch has remodeled the protein surface in contact with the substrate, potentially for binding highly substituted arabinoxylan. We propose that the evolution of the divergent EXPB-I group made a fundamental change in the target and mechanism of wall loosening in the grass lineage possible, involving a new structural role for xylans and the expansins that target them.

Oil accumulation by the oleaginous diatom Fistulifera solaris as revealed by the genome and transcriptome

Oleaginous photosynthetic organisms such as microalgae are promising sources for biofuel production through the generation of carbon-neutral sustainable energy. However, the metabolic mechanisms driving high-rate lipid production in these oleaginous organisms remain unclear, thus impeding efforts to improve productivity through genetic modifications. We analyzed the genome and transcriptome of the oleaginous diatom Fistulifera solaris JPCC DA0580. Next-generation sequencing technology provided evidence of an allodiploid genome structure, suggesting unorthodox molecular evolutionary and genetic regulatory systems for reinforcing metabolic efficiencies. Although major metabolic pathways were shared with nonoleaginous diatoms, transcriptome analysis revealed unique expression patterns, such as concomitant upregulation of fatty acid/triacylglycerol biosynthesis and fatty acid degradation (β-oxidation) in concert with ATP production. This peculiar pattern of gene expression may account for the simultaneous growth and oil accumulation phenotype and may inspire novel biofuel production technology based on this oleaginous microalga.

The bimodal distribution of genic GC content is ancestral to monocot species

In grasses such as rice or maize, the distribution of genic GC content is well known to be bimodal. It is mainly driven by GC content at third codon positions (GC3 for short). This feature is thought to be specific to grasses as closely related species like banana have a unimodal GC3 distribution. GC3 is associated with numerous genomics features and uncovering the origin of this peculiar distribution will help understanding the potential roles and consequences of GC3 variations within and between genomes. Until recently, the origin of the peculiar GC3 distribution in grasses has remained unknown. Thanks to the recent publication of several complete genomes and transcriptomes of nongrass monocots, we studied more than 1,000 groups of one-to-one orthologous genes in seven grasses and three outgroup species (banana, palm tree, and yam). Using a maximum likelihood-based method, we reconstructed GC3 at several ancestral nodes. We found that the bimodal GC3 distribution observed in extant grasses is ancestral to both grasses and most monocot species, and that other species studied here have lost this peculiar structure. We also found that GC3 in grass lineages is globally evolving very slowly and that the decreasing GC3 gradient observed from 5′ to 3′ along coding sequences is also conserved and ancestral to monocots. This result strongly challenges the previous views on the specificity of grass genomes and we discuss its implications for the possible causes of the evolution of GC content in monocots.

Exploiting transcriptome data for the development and characterization of gene-based SSR markers related to cold tolerance in oil palm (Elaeis guineensis)

BACKGROUND

The oil palm (Elaeis guineensis, 2n = 32) has the highest oil yield of any crop species, as well as comprising the richest dietary source of provitamin A. For the tropical species, the best mean growth temperature is about 27°C, with a minimal growth temperature of 15°C. Hence, the plantation area is limited into the geographical ranges of 10°N to 10°S. Enhancing cold tolerance capability will increase the total cultivation area and subsequently oil productivity of this tropical species. Developing molecular markers related to cold tolerance would be helpful for molecular breeding of cold tolerant Elaeis guineensis.

RESULTS

In total, 5791 gene-based SSRs were identified in 51,452 expressed sequences from Elaeis guineensis transcriptome data: approximately one SSR was detected per 10 expressed sequences. Of these 5791 gene-based SSRs, 916 were derived from expressed sequences up- or down-regulated at least two-fold in response to cold stress. A total of 182 polymorphic markers were developed and characterized from 442 primer pairs flanking these cold-responsive SSR repeats. The polymorphic information content (PIC) of these polymorphic SSR markers across 24 lines of Elaeis guineensis varied from 0.08 to 0.65 (mean = 0.31 ± 0.12). Using in-silico mapping, 137 (75.3%) of the 182 polymorphic SSR markers were located onto the 16 Elaeis guineensis chromosomes. Total coverage of 473 Mbp was achieved, with an average physical distance of 3.4 Mbp between adjacent markers (range 96 bp - 20.8 Mbp). Meanwhile, Comparative analysis of transcriptome under cold stress revealed that one ICE1 putative ortholog, five CBF putative orthologs, 19 NAC transcription factors and four cold-induced orhologs were up-regulated at least two fold in response to cold stress. Interestingly, 5' untranslated region of both Unigene21287 (ICE1) and CL2628.Contig1 (NAC) both contained an SSR markers.

CONCLUSIONS

In the present study, a series of SSR markers were developed based on sequences differentially expressed in response to cold stress. These EST-SSR markers would be particularly useful for gene mapping and population structure analysis in Elaeis guineensis. Meanwhile, the EST-SSR loci were inducible expressed in response to low temperature, which may have potential application in identifying trait-associated markers in oil palm in the future.

Expression Comparison of Oil Biosynthesis Genes in Oil Palm Mesocarp Tissue Using Custom Array

Gene expression changes that occur during mesocarp development are a major research focus in oil palm research due to the economic importance of this tissue and the relatively rapid increase in lipid content to very high levels at fruit ripeness. Here, we report the development of a transcriptome-based 105,000-probe oil palm mesocarp microarray. The expression of genes involved in fatty acid (FA) and triacylglycerol (TAG) assembly, along with the tricarboxylic acid cycle (TCA) and glycolysis pathway at 16 Weeks After Anthesis (WAA) exhibited significantly higher signals compared to those obtained from a cross-species hybridization to the Arabidopsis (p-value < 0.01), and rice (p-value < 0.01) arrays. The oil palm microarray data also showed comparable correlation of expression (r² = 0.569, p < 0.01) throughout mesocarp development to transcriptome (RNA sequencing) data, and improved correlation over quantitative real-time PCR (qPCR) (r² = 0.721, p < 0.01) of the same RNA samples. The results confirm the advantage of the custom microarray over commercially available arrays derived from model species. We demonstrate the utility of this custom microarray to gain a better understanding of gene expression patterns in the oil palm mesocarp that may lead to increasing future oil yield.