# Transcriptome Analysis of Cell Wall and NAC Domain Transcription Factor Genes during Elaeis guineensis Fruit Ripening: Evidence for Widespread Conservation within Monocot and Eudicot Lineages

Timothy J Tranbarger, Kim Fooyontphanich, Peerapat Roongsattham, Maxime Pizot, Myriam Collin, Chatchawan Jantasuriyarat, Potjamarn Suraninpong, Somvong Tragoonrung, Stéphane Dussert, Jean-Luc Verdeil, et al.

<figure><img src="image_1.png" /></figure>

To cite this version:

Timothy J Tranbarger, Kim Fooyontphanich, Peerapat Roongsattham, Maxime Pizot, Myriam Collin, et al. Transcriptome Analysis of Cell Wall and NAC Domain Transcription Factor Genes during Elaeis guineensis Fruit Ripening: Evidence for Widespread Conservation within Monocot and Eudicot Lineages. Frontiers in Plant Science, 2017, 8, pp.613. 10.3389/fpls.2017.00603. hal-03204464

HAL Id: hal-03204464

https://hal.science/hal-03204464

Submitted on 21 Apr 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

# Transcriptome Analysis of Cell Wall and NAC Domain Transcription Factor Genes during *Elaeis guineensis* Fruit Ripening: Evidence for Widespread Conservation within Monocot and Eudicot Lineages

OPEN ACCESS

**Edited by:**
Antonio Ferrante,
Università degli Studi di Milano, Italy

**Reviewed by:**
Xiaohong Kou,
Tianjin University, China
Anna Mensuali,
Sant'Anna School of Advanced
Studies, Italy

**\*Correspondence:**
Timothy J. Tranbarger
timothy.tranbarger@ird.fr

**†Present Address:**
Peerapat Roongsattham,
Department of Biology, Faculty of
Science, Prince of Songkla University,
Songkhla, Thailand
Kim Fooyontphanich,
Kasetsart University and Winrock
International, Bangkok, Thailand

**Specialty section:**
This article was submitted to
Plant Physiology,
a section of the journal
Frontiers in Plant Science

Received: 09 December 2016

Accepted: 03 April 2017

Published: 25 April 2017

**Citation:**
Tranbarger TJ, Fooyontphanich K,
Roongsattham P, Pizot M, Collin M,
Jantasuriyarat C, Suraninpong P,
Tragoonrung S, Dussert S, Verdeil J-L
and Morcillo F (2017) Transcriptome
Analysis of Cell Wall and NAC Domain
Transcription Factor Genes during
Elaeis guineensis Fruit Ripening:
Evidence for Widespread
Conservation within Monocot and
Eudicot Lineages.
Front. Plant Sci. 8:603.
doi: 10.3389/fpls.2017.00603

Timothy J. Tranbarger<fnref n="1" />*, Kim Fooyontphanich<fnref n="1" />†, Peerapat Roongsattham<fnref n="1" />†,
Maxime Pizot<fnref n="1" />, Myriam Collin<fnref n="1" />, Chatchawan Jantasuriyarat², Potjamarn Suraninpong³,
Somvong Tragoonrung⁴, Stéphane Dussert<fnref n="1" />, Jean-Luc Verdeil⁵ and Fabienne Morcillo⁶

The oil palm (*Elaeis guineensis*), a monocotyledonous species in the family Arecaceae, has an extraordinarily oil rich fleshy mesocarp, and presents an original model to examine the ripening processes and regulation in this particular monocot fruit. Histochemical analysis and cell parameter measurements revealed cell wall and middle lamella expansion and degradation during ripening and in response to ethylene. Cell wall related transcript profiles suggest a transition from synthesis to degradation is under transcriptional control during ripening, in particular a switch from cellulose, hemicellulose, and pectin synthesis to hydrolysis and degradation. The data provide evidence for the transcriptional activation of expansin, polygalacturonase, mannosidase, beta-galactosidase, and xyloglucan endotransglucosylase/hydrolase proteins in the ripening oil palm mesocarp, suggesting widespread conservation of these activities during ripening for monocotyledonous and eudicotyledonous fruit types. Profiling of the most abundant oil palm polygalacturonase (*EgPG4*) and 1-aminocyclopropane-1-carboxylic acid oxidase (ACO) transcripts during development and in response to ethylene demonstrated both are sensitive markers of ethylene production and inducible gene expression during mesocarp ripening, and provide evidence for a conserved regulatory module between ethylene and cell wall pectin degradation. A comprehensive analysis of NAC transcription factors confirmed at least 10 transcripts from diverse NAC domain clades are expressed in the mesocarp during ripening, four of which are induced by ethylene treatment, with the two most inducible (*EgNAC6* and *EgNAC7*) phylogenetically similar to the tomato NAC-NOR master-ripening regulator. Overall, the

<footnote n="1">¹Institut de Recherche pour le Développement, IRD, UMR DIADE, Montpellier, France, ²Department of Genetics, Kasetsart University, Bangkok, Thailand, ³Department of Plant Science, Institute of Agricultural Technology, Walailak University, Nakhon Si Thammarat, Thailand, ⁴Genome Institute, National Center for Genetic Engineering and Biotechnology, Pathumthani, Thailand, ⁵Centre de Coopération Internationale en Recherche Agronomique pour le Développement, UMR AGAP, Montpellier, France, ⁶Centre de Coopération Internationale en Recherche Agronomique pour le Développement, UMR DIADE, Montpellier, France</footnote>
results provide evidence that despite the phylogenetic distance of the oil palm within the family Arecaceae from the most extensively studied monocot banana fruit, it appears ripening of divergent monocot and eudicot fruit lineages are regulated by evolutionarily conserved molecular physiological processes.

Keywords: ripening, cell wall, NAC domain, oil palm, monocotyledon, mesocarp, ethylene

INTRODUCTION

Fruit ripening is a biological character unique to flowering plants, not only of central importance to seed dispersal and the reproductive success of plants, but also essential for both human and animal diets. The color, metabolic, and textural transitions that occur in the ripening fleshy fruit tissues have economic and nutritional consequences for humans. Ethylene is known to play a central role in the ripening of climacteric fruit in which there is a shift from a basal non-catalytic (system 1) to an autocatalytic (system 2) increase in ethylene production (Lelievre et al., 1997; Klee and Giovannoni, 2011; Liu et al., 2015). Studies with the climacteric fruit model tomato indicate a central role for ethylene in the transcriptional coordination of metabolic processes that occur during ripening (Alba et al., 2005; Osorio et al., 2011). The role of ethylene is conserved across many eudicot species and is likely to play similar roles in other climacteric fleshy fruits such apple, peach, pear, mango, and banana (Prasanna et al., 2007; Bapat et al., 2010; Seymour et al., 2013a,b). Furthermore, recent studies indicate that ethylene signal transduction is also important for the transcriptional regulation or ripening process in non-climacteric fruit that lack the large climacteric increase in ethylene, such as pepper, grape, and strawberry (Trainotti et al., 2005; Chervin et al., 2008; Bapat et al., 2010; Osorio et al., 2012). While ethylene may be viewed as a common regulator for fruit ripening in general, how the ripening related regulatory networks have evolved and adapted during plant evolution is unknown. Indeed, most of the research on fruit ripening focuses on the important cultivated eudicotyledonous species, in particular tomato. However, research with the monocotyledonous banana and the basal angiosperm *Persea americana* (avocado) support the notion that the role of ethylene has been conserved through plant evolution, but given the large phylogenetic distances of these species, diversification can be expected to be discovered (Chanderbali et al., 2008, 2009; Elitzur et al., 2010, 2016; Jourda et al., 2014).

One of the most important processes regulated by ethylene during fruit ripening are the textural changes that occur due to cell wall disassembly, indeed, transgenic tomato with modified ethylene biosynthesis or perception have provided evidence that ethylene regulated cell wall metabolism is central to climacteric fruit ripening (Bennett and Labavitch, 2008). However, while extensive work has been done with the ripening tomato as a eudicot model, our understanding of how the cell wall is modified for softening to occur is still incomplete. In general, it appears that textural changes during ripening are controlled simultaneously by many genes, and may be species and/or fruit type

dependent (Bennett and Labavitch, 2008; Seymour et al., 2013a).

NAC-NOR is an example of a master regulator that controls ethylene production and downstream changes in the cell wall of the ripening tomato (Giovannoni, 2004; Seymour et al., 2013a). The *non-ripening* (*nor*) tomato mutation affects a NAC domain transcription factor (TF) and results in fruit that do not ripening normally, similar to both *ripening-inhibitor* (*rin*) and *never-ripe* (*nr*) (Giovannoni, 2004; Osorio et al., 2011). All these mutants share a strong down regulation of the tomato fruit polygalacturonase (TFPG) transcript (*PG2*), indicating that the transcriptional activation of the tomato *PG2* gene is an important control point regulating the expression of *PG2* during ripening (Dellapenna et al., 1986, 1987, 1989). While PG can hydrolyze the backbone of pectin homogalacturonan polymers, the down-regulation or knockout of the TFPG transcript results in a decrease in pectin depolymerization, but no change in fruit softening, and it is generally thought that a large number of downstream cell wall related target genes are activated transcriptionally during ripening and are necessary for softening (Sheehy et al., 1988; Smith et al., 1988, 1990; Cooley and Yoder, 1998; Brummell and Harpster, 2001; Brummell, 2006; Seymour et al., 2013a). NAC domain proteins have emerged as important transcriptional regulators of tomato fruit ripening. For example, the silencing of two NAC domain TFs *SNAC4* and *SNAC9* results in decreased expression of the tomato ethylene biosynthesis genes *LeACS2*, *LeACS4* (encoding 1-amino-cyclopropane-1-carboxylic acid or ACC synthases), and *LeACO1* (tomato ACC oxidase), which consequently inhibition the ripening process (Kou et al., 2016). Indeed, the SNAC TFs bind to the promoter of these ethylene biosynthesis genes *in vitro*, while the silencing of *LeACS4*, *LeACO1*, and *LeERF2* reduced the expression of *SNAC4* and *SNAC9*, suggesting a possible feedback system between these SNAC TFs and ethylene biosynthesis gene transcription. In contrast, less is known about monocot fruit ripening. Recent studies found NAC domain transcripts were differentially expressed in the ripening oil palm mesocarp, suggesting ethylene-regulated expression, while a study with banana identified at least 6 NAC domain TFs expressed during banana ripening, two of which are induced by ethylene (Tranbarger et al., 2011; Shan et al., 2012). Together these results suggest that NAC domain TFs may have conserved central roles during ethylene induced fruit ripening of both eudicots and monocots.

The monocot-eudicot divergence occurred 130-135 million years ago (D'hont et al., 2012; Magallón et al., 2015). The best-studied monocot fruit is banana in which ethylene plays a conserved role during ripening (Elitzur et al., 2010, 2016; Jourda et al., 2014). Arecaceae (the monocot palm family) displays a very large amount of fruit diversity and is separated from its closest
monocot lineage by approximately 98 million years (Magallón et al., 2015). Ethylene appears to be involved in fruit ripening of both oil palm (Elaeis guineensis) and date palm (Phoenix dactylifera L), while oil palm has twice the number of ripening and abscission genes and more ACO genes than date palm (Abbas and Ibrahim, 1996; Tranbarger et al., 2011; Singh et al., 2013; Nualwijit and Lerslerwong, 2014). Recent studies by our group identified a polygalacturonase (EgPG4) highly induced by ethylene in oil palm fruit abscission zone cells and associated with the cell separation and fruit abscission (Roongsattham et al., 2012, 2016). Whereas, EgPG4 is not expressed or induced in the adjacent pedicel tissue, it is also induced by ethylene and highly expressed in the mesocarp cells (Roongsattham et al., 2012). In addition, quantitative changes in the mesocarp cell walls occur during development and ripening, including thickness increases in the primary cell wall and middle lamella, while ethylene treatment of ripe fruit results in a decrease in middle lamella thickness (Roongsattham et al., 2016).

The majority of our knowledge about ethylene function during fruit ripening in higher plants comes from studies with eudicot model plants, in particular with the tomato fruit, and to a lesser extent with the monocot model fruit banana. To what extent this knowledge extends to non-model and phylogenetically distant species is not fully known. The current study focuses on the ethylene and cell wall related processes that take place during the ripening of the oil palm fruit mesocarp in comparison with what is known from the main eudicot and monocot fruit models, banana, and tomato.

# MATERIALS AND METHODS

## Plant Material, Treatments, and RNA Extraction

Oil palm fruits for the mesocarp transcriptome were harvested at Pobè CRA-PP Station (INRAB) Benin, from a *dura* parent of Deli Dabou origin, within the same self-progeny of a single palm as described previously (Tranbarger et al., 2011). Oil palm spikelets with ripe fruit 180 days after pollination (DAP), at which time fruit are typically harvested and transported for oil extraction, were sampled from a tenera clone (clone C) at the Krabi Golden Tenera plantation, Thailand as previously described (Roongsattham et al., 2012). Spikelets were treated for 3, 6, and 9 h with 10 µl l⁻¹ ethylene, a quantity of ethylene that was previously found to induce changes to the cell walls of the fruit abscission zone cells that results in 25% fruit abscission after 9 h (Roongsattham et al., 2012, 2016). For RNA extractions, mesocarp fruit tissue samples were collected and frozen immediately in liquid nitrogen. Total RNA from mesocarp was extracted as previously described (Morcillo et al., 2006). Total RNA (1 µg) was used to synthesize cDNA using the first-strand cDNA synthesis kit (ImProm-II™ Reverse Transcription System, Promega).

## Quantitative Real-Time RT-PCR

qPCR was conducted as previously described (Roongsattham et al., 2012). The analysis was conducted on a LightCycler 480 (Roche) in 96-well plates in a volume of 10 µl containing

2 µl of cDNA diluted 1/100, 1.5 µl of primer forward (2 µM), 1.5 µl of reverse primer (2 µM), and 5 µl SYBR® Green Mastermix (Roche). Primers used for EgPG4 and ACO amplification were previously published (Fooyontphanich et al., 2016). PCR was initiated by denaturation at 95°C for 10 min, followed by 45 cycles of 95°C for 15 s, 60°C for 15 s, and a final extension at 70°C for 1 min. All expression was normalized to the *EgEfα1* (accession number: AY550990) mRNA from *E. guineensis*, and mRNA abundance for each experiment was calculated relative to the sample with the lowest amount of transcript present determined with the formula as described previously (Pfaffl, 2001). No change of *EgEfα1* transcript accumulation was found in the fruit tissues treated or not treated with ethylene. Control using RNA matrices were also conducted to validate the absence of DNA in each sample. Each time point was replicated three times from two independent biological samples, and all amplified cDNA fragments were sequenced by Beckman-Cogenics (https://www. genewiz.com) to check the specificity of the amplified products. Gene abundance is expressed as mean and standard error bars are calculated from the technical replicates of one of the biological repetitions.

## Mesocarp Transcriptome Data Mining

Transcriptome data of the developing mesocarp previously clustered (clusters A, B, C, and D) was searched for transcripts with expression profiles that either increase or decrease during the burst of ethylene production observed during oil palm fruit ripening between 100 and 160 DAP (Supplementary Table 1; Tranbarger et al., 2011). The BLAST2GO and InterProScan web services with BLASTX using an E-value cutoff of 1e-5 were used to annotate the gene sets (Altschul et al., 1990; Zdobnov and Apweiler, 2001; Götz et al., 2008). Cell wall sequences were identified by searching the GO annotated sequences for InterPro accessions and key words related to cell wall processes, and by searching (TBLASTX) the 454 sequence database with known candidates related to cell wall biosynthesis and degradation.

## RNA-Seq Data Processing

Nine RNAseq datasets generated by Orion Genomics (http://www.oriongenomics.com/index.php/en/) for selected *E. guineensis* tissues were downloaded from NCBI databases: 17 cm long leaf (SRX278048), spear leaf (SRX278049), 2.5 cm (SRX278052), and 20 cm long female inflorescence (SRX278053), pollen (SRX278051), root (SRX278062), shoot apex (SRX278055), and kernel at 10 and 15 weeks after anthesis (SRX278021 and SRX278018, respectively). Reads from these datasets in addition to RNAseq data (100-nt Illumina reads) for the *E. guineensis* ripening mesocarp (Guerin et al., 2016) and embryonic cell suspensions (kindly provided by Thierry Beulé) were then mapped on the *E. guineensis* CDS reference (NCBI website, GCF_000442705.1_EG5_rna.fna; January 2015) and the RPKM (number of reads per kilobase and million reads) for each locus were compared.
## Transcription Factor Prediction

The NAC domain gene family of *E. guineensis* was retrieved from the Plant Transcription Factor Database (PlantTFDB, http://planttfdb.cbi.pku.edu.cn/prediction.php) 4.0 using the *E. guineensis* CDS reference (NCBI website, GCF_000442705.1_EG5_rna.fna; January 2015; Jin et al., 2015, 2016). From a total of 44,360 oil palm putative coding sequences (CDS), 175 NAC domain-containing sequences were identified in the oil palm reference CDS, which corresponded to a total of 124 non-redundant loci. The same result was obtained using the iTAK database (Zheng et al., 2016; http://bioinfo.bti.cornell.edu/cgi-bin/itak/online_itak.cgi).

## Phylogenetic Analysis

Amino acid alignments were performed with ClustalW (Larkin et al., 2007) using default settings (http://www.ebi.ac.uk/). The amino acid sequences containing the conserved NAC subdomains (A-E) were used for the phylogenetic analysis (Ooka et al., 2003) using the default settings without the G-blocks step (Dereeper et al., 2008). Branch support values are based on the aLRT statistical test (Anisimova and Gascuel, 2006).

## Histological Analysis and Cell Measurements

Mesocarp fruit samples were collected from ethylene treated fruit and fixed in 0.2 M phosphate buffer containing 2% (w/v) paraformaldehyde, 1% (w/v) caffeine, and 0.5% (v/v) glutaraldehyde or a minimum of 2 days at 4°C as previously described (Buffard-Morel et al., 1992). Serial dehydration with ethanol from 30 to 100%, then 100% butanol/100% ethanol (v/v), and finally 100% butanol was performed for each sample and followed by impregnation and embedding in Technovit 7100 resin (Heraeus Kulzer). Semi-thin sections of 3 µm were cut using a microtome. Each section was stained with toluidine blue or ruthenium red. Microphotographs were taken with a Leica camera (DFC 300 FX) on a Leica (LEITZ DMRB) light microscope (x20/0.5; x40/0.7; and x100/1.3). Cell parameter measurements were performed for cell width, cell wall width, and middle lamella as previously reported (Roongsattham et al., 2016).

# RESULTS

## EgACO and Polygalacturonase EgPG4 Transcript Profiling during Mesocarp Development and in Response to Ethylene Treatment of Ripe Fruit

To provide molecular indicators of ethylene production in the mesocarp, we examined the quantitative expression profiles of the most highly abundant ACO previously observed in the oil palm mesocarp (CL1Contig999, *EgACO*), in addition to the most highly expressed polygalacturonase (CL1Contig5616, *EgPG4*) in the mesocarp induced by ethylene (Tranbarger et al., 2011; Roongsattham et al., 2012; Fooyontphanich et al., 2016) (Supplementary Tables 1, 2). The *EgACO* nucleotide sequence aligns most significantly (E-value 0) to LOC105039092

(Chromosome 2) annotated as an ACO3 in the oil palm genome (Supplementary Table 2). During mesocarp development, the *EgACO* transcript is expressed in the pollinated flower at the onset of fruit development (10 DAP) and is undetectable or barely detectable until 120 DAP, and then increases 79 and 199 fold at 140 and 160 DAP respectively (Figure 1A). Illumina sequencing confirmed that the *EgACO3* transcript increased between 120 and 160 DAP (Figure 1A, Supplementary Table 2). The pattern of *EgACO* expression correlates to the increase in ethylene observed between 120 and 160 DAP, which corresponds to the transition from system 1 to system 2 ethylene production (Tranbarger et al., 2011). To determine whether the *EgACO* transcript could be regulated by ethylene, the amount of *EgACO* transcript was quantified in the mesocarp of ripe fruit (150 and 180 DAP) treated with ethylene for 3-9 h (Figure 1B). After 3 h of ethylene treatment, in 150 DAP fruit the *EgACO* transcript increased 25 fold, while in the 180 DAP fruit, the transcript increased 9 fold. A further increase was observed at 9 h while at 6 h the *EgACO* transcript amounts remained stable (150 DAP) or decreased (180 DAP). In contrast to *EgACO*, the *EgPG4* transcript is undetectable or present at very low amounts from 10 to 140 DAP, then increases dramatically over 3,000 fold at 160 DAP, which corresponds to the autocatalytic system 2 ethylene production (Figure 1B). *EgPG4* is annotated at LOC105034919 (Chromosome 2, Supplementary Table 2). Illumina sequencing confirmed an increase of *EgPG4* transcript increased between 140 and 160 DAP (Figure 1A, Supplementary Table 2). The *EgPG4* transcript increased in both 150 and 180 DAP fruit after 3 h of ethylene treatment, but the magnitude was much higher in 150 DAP (500 fold) than in the older 180 DAP fruit (10 fold). The *EgPG4* transcript continued to increase after 6 and 9 h of ethylene treatment and eventually reached over 5,000 fold by 9 h in 150 DAP fruit. While both the *EgACO* and *EgPG4* transcripts increase rapidly (3 h) in response to ethylene treatment, the induction was less in 180 DAP fruit for both genes, while the magnitude of induction was consistently higher for *EgPG4*. Furthermore, *EgPG4* expression is limited to the ripest fruit in contrast to *EgACO*, which is also expressed 10 days after the flower is pollinated.

## Histological Analysis of the Mesocarp in Response to Ethylene in Ripening Fruit

Histological analysis of abscission zone and mesocarp cell walls revealed changes in integrity during development and in response to ethylene treatments (Roongsattham et al., 2016). While the previous article focused on cellular changes in abscission zone cells, here we focus on changes that occur in the mesocarp cells during ripening and in response to ethylene (Supplementary Figures 1, 2). To examine the cellular processes in the mesocarp, longitudinal tissue sections were analyzed from samples of mesocarp at key stages of development and ripening, and in response to ethylene. At 30 DAP, anticlinal cell divisions have taken place, cells have not fully elongated, and the formation of intercellular spaces is observed (Roongsattham et al., 2016; Supplementary Figure 1). By 120 DAP, cells have
<figure><img src="image_4.png" /><figcaption>FIGURE 1 | *EgACO* and *EgPG4* expression in the oil palm mesocarp increase during fruit ripening and are induced by external ethylene treatments at two different ripening stages. qPCR analysis of *EgACO* and *EgPG4* transcript abundance in the mesocarp during development (A) and in response to exogenous ethylene (10 µL L⁻¹) treatments (B). Inserts are RNA-Seq based RPKM profiles of the mesocarp at 100, 120, 140, and 160 DAP.</figcaption></figure>

elongated and a significant increase in the cell wall width is observed, with the average cell wall width more than 1 µm while the middle lamella is approximately 0.6 µm wide and the presence of intercellular spaces is also observed at 120 DAP. At 180 DAP, the cell wall width increased, and cell separation between adjacent cells is apparent. By 180 DAP a 2-fold increase in cell wall width is observed, and the middle lamella width also increases significantly to ~0.9 µm. Based on these measurements, the cell wall accounts for approximately 4, 7, and 15% of the total cell width at 30, 120, and 180 DAP respectively. After 9 h of ethylene treatment of 180 DAP fruit,

the toluidine blue staining increases between adjacent cells, while ruthenium red stained sections reveal non-continuous dark strands that corresponds to an apparent breakdown of the middle lamella that is measurable after 9 h of ethylene treatment (Figure 2). Finally, in ethylene treated cells, intercellular material accumulates based on both toluidine blue and ruthenium red staining (Figure 2). Based on these data, a significant amount of cell wall and the middle lamella expansion occurs between 30 to 180 DAP, while the integrity of the middle lamella in the mesocarp cells is most affected by ethylene treatment of 180 DAP fruit (Figure 2).
<figure><img src="image_5.png" /><figcaption>**FIGURE 2** | Mesocarp cell characteristics change during development and after 9 h ethylene treatment. A longitudinal sections of mesocarp cells toward the base of the fruit stained with toluidine blue (upper panels) and ruthenium red (lower panels). cs, cell separation; dc, divided cell; IM, intercellular material; IS, intercellular space; ml, middle lamella; Scale bar, 10 µm.</figcaption></figure>

## Cell Wall Associated Transcriptional Activity during Mesocarp Development and Ripening

Based on the histological analysis, an increase in cell wall width occurs during mesocarp development and ripening, while exogenous ethylene treatment results in a loss of middle lamella integrity (Figure 2; Roongsattham et al., 2016). To provide insight into the transcriptional basis for these processes, proteins associated with cell wall assembly, metabolism, and modification were searched for within the differentially abundant transcripts previously identified in the mesocarp transcriptome at 100, 120, 140, and 160 DAP (Supplementary Table 1; Tranbarger et al., 2011). A total of 75 transcripts for cell wall related activities were found to be differentially expressed in the mesocarp during ripening, 63% of which have expression peaks at 140 and 160 DAP (including *EgPG4*) concomitant with the ethylene burst as measured previously (Tranbarger et al., 2011). Illumina sequencing at 100, 120, 140, and 160 DAP mesocarp validated the expression 54 of these genes annotated in the oil palm genome with roles in cell wall component biosynthesis, assembly, modification and hydrolysis during fruit ripening (Figure 3; Supplementary Table 2; Guerin et al., 2016). A global comparison of each transcript type expressed as a percentage of all cell wall transcripts gives an estimation of the transcriptional contribution of each cell wall activity during these stages of

<figure><img src="image_6.png" /></figure>

fruit ripening (Figure 4; Supplementary Table 3). Globally, the most abundant transcripts encode sequences similar to glucan endo-1,3-beta-glucosidase, with 25% of the total read counts from the 4 stages of development, and a peak of 43% at 160
<figure><img src="image_7.png" /><figcaption>**FIGURE 4** | Global cell wall related transcript abundance during mesocarp development and ripening. The transcripts and their encoded products can be grouped into two categories including biosynthesis (upper panel) and hydrolysis (lower panel) related functions. The upper black arrow indicates the stage at which the transition from ethylene system 1 to system 2 occurs while the lower arrow indicates the stage at which the ethylene burst occurs (Tranbarger et al., 2011). See **Supplementary Tables 1, 2** for details.</figcaption></figure>

DAP. Chitinase and pathogenesis-related proteins accounted for 13% of the total, with 16% and 15% at 140 DAP and 160 DAP respectively, while xyloglucan biosynthesis or modification proteins accounted for 12% of the total with 32% at 120 DAP (Figure 4, Supplementary Table 3). Xyloglucan biosynthesis or modification proteins were also highly abundant (27%) along with cellulose synthase and cellulose synthase-like transcripts (26%) at 100 DAP. Transcripts for cellulose synthase and cellulose-synthase like proteins decreased progressively during the later stages to 12% (120 DAP), 5% (140 DAP), and 2% (160 DAP). Transcripts for uridine diphosphate (UDP) glucose 6-dehydrogenase (9%), involved in both pectin and hemicellulose biogenesis, were more abundant at 120 DAP, but also present throughout the developmental stages examined (Klinghammer and Tenhaken, 2007). Pectin related transcripts included pectin methyltransferases (PMTs, 13% 120 DAP), a pectin methylesterase (PME, 3% 120 DAP) and the transcript

EgPG4, which represents 20% of the total cell wall related transcripts at 160 DAP (Figure 3). Overall, there is a progressive change from cell wall biosynthesis to hydrolysis during the transition from ethylene system 1 to 2 (between 120 and 140 DAP) that includes transcripts encoding glucan endo-1,3-beta-glucosidase, endoglucanase, beta- and alpha-galactosidase, chitinase and pathogenesis-related, mannan endo-1,4-beta-mannosidase and expansin proteins. Between 140 and 160 DAP, the transcripts that increased the most were EgPG4 (210 fold) and glucan endo-1,3-beta-glucosidase-like (LOC105033964, ~350 fold; Supplementary Table 3). In addition, one transcript similar to a xyloglucan endotransglucosylase/hydrolase protein 23 (LOC105039293) transcript was expressed at relatively high amounts exclusively at 160 DAP. The profile similarity of EgPG4 with other transcripts suggest certain cell wall activities may be coordinated by ethylene during mesocarp ripening (Figure 4, Supplementary Tables 2, 3).
# Analysis of NAC Transcription Factor Transcripts Expressed in the Mesocarp and in Response to Ethylene

NAC domain containing TFs were found previously to be one of the most highly represented TFs expressed in the mesocarp (Tranbarger et al., 2011). To respond to the question whether any of these putative NAC TFs could be regulated by ethylene, first a genome wide search for NAC domain TFs was performed, then expression profiles of selected candidates were examined by qPCR in fruit treated with ethylene at two stages of ripening (Figure 5, Supplementary Table 4). A total of 124 putative NAC domain encoding genes were identified in the oil palm genome, 33 of which had total read counts higher than 50 in the ripening mesocarp (Supplementary Table 4). A selection was made that consisted of the top five mesocarp expressed genes in addition to five additional genes with either down regulated or up regulated expression profiles in the ripening mesocarp. Gene specific primer sets amplified 10 non-redundant NAC sequences and confirmed their expression in the mesocarp (Supplementary Table 5). As with the studies on *EgACO* and *EgPG4*, ripe fruit at two stages (150 and 180 DAP) were treated with ethylene and the NAC transcripts were quantified. All the transcripts examined either increased or decreased at one time point in 150 or 180 DAP fruit. Of the 10 NAC transcripts examined, 7 NAC transcripts (*EgNAC2*, 3, 5, 6, 7, 8, and 13) were significantly induced at one time point in either 150 or 180 DAP fruit. Notable, six of these (*EgNAC2*, 3, 5, 6, 7, and 8) were induced after 3 h exposure to ethylene in the 180 DAP fruit, while two of those (*EgNAC5* and 6) were also induced at 3 h in 150 DAP fruit. In addition, the *EgNAC5*, 6, 7, and 8 increased in abundance after 9 h of ethylene exposure in both 150 and 180 DAP fruit. *EgNAC6* and 7 were the most highly expressed transcripts expressed in response to ethylene in the mesocarp at both stages of ripening (150 and 180 DAP).

To provide insight into the diversification of NAC domain TFs, and to determine whether there is a relation between NAC domain structure and ethylene response, a phylogenetic analysis of the oil palm NAC domain sequences along with those of *Arabidopsis* and other selected NAC TFs with was performed (Supplementary Figure 2). Strong branch support values provided evidence for at least two major subclades of NAC domain sequences that resolve into a number of subgroups previously identified (Jensen et al., 2010). The oil palm mesocarp NAC domain sequences are distributed throughout the two major subclades, while six are found within the subgroups III-2 and III-3 (Supplementary Figure 2). NAC1, 2, and 3 group closely in subgroup III-3 and are also less abundant in the 150 DAP mesocarp after 3-6 h treatments with ethylene (Figure 5). In contrast, two of the most highly ethylene inducible NAC transcripts (*EgNAC6* and 7) encode proteins found within subgroup III-2, along with the tomato NAC-NOR domain (Figures 5; Supplementary Figure 2). Among the other highly inducible transcripts; EgNAC5 separates in a less resolved portion of the cladogram, while EgNAC8 separates within subgroup II-3. It appears that an expansion of oil palm NAC domains occurred within the subgroups III-2 and III-3, where six sequences of

the current study are found. In particular, two highly ethylene inducible genes, *EgNAC6* and *EgNAC7*, appear to be paralogous and encode NAC domain proteins the most homologous to the tomato NAC-NOR.

## DISCUSSION

### The Transition from Ethylene Production System 1 to 2 Occurs Concomitantly with a Transcriptional Activation of *EgACO*

The transition from basal auto inhibited ethylene production (system 1) to the burst of autocatalytic ethylene production (system 2) represents a central regulatory process that coordinates the expression of ripening related genes involved in pigmentation, aroma, carbohydrate, and cell wall metabolism (Seymour et al., 2013a; Liu et al., 2015). In tomato, the transition from system 1 to 2 is controlled both by developmental ethylene-independent and ethylene dependent factors and involves the transcriptional activation of specific ACS and ACO genes by the ripening related TFs Ripening-Inhibitor (RIN) and HB1 respectively (Nakatsuka et al., 1998; Barry and Giovannoni, 2007; Cara and Giovannoni, 2008; Lin et al., 2008; Yokotani et al., 2009; Klee and Giovannoni, 2011; Martel et al., 2011; Seymour et al., 2013a,b; Liu et al., 2015). The *EgACO* transcript is present in low amounts in the pollinated flowers (10 DAP), barely detectable or undetectable in at early stages of mesocarp development (30-100 DAP), and most highly expressed during the stages of fruit ripening (120-160 DAP), consistent with a role in both flower senescence and ripening similar to tomato (Blume and Grierson, 1997). However, the transcript profile differs from that observed during tomato ripening, in which the main ripening related *LeACO1* is expressed during both ethylene system 1 and system 2 stages of fruit development, with the highest expression during ripening and the autocatalytic system 2 (Holdsworth et al., 1987; Barry et al., 1996; Blume and Grierson, 1997; Nakatsuka et al., 1998; Jafari et al., 2013). Indeed, *LeACO1* is thought to participate in both system 1 and 2, and is not specific to ripening (Barry and Giovannoni, 2007; Cara and Giovannoni, 2008). Similarly in banana, the principle *MA-ACO1* transcript is also present in system 1 stage fruit, increases during the transition from system 1 to 2, and remains high at later stages of ripening (Liu et al., 1999). Our previous analysis revealed a burst of ethylene occurred in the oil palm fruit between 120 and 160 DAP, along with changes in transcript profiles related to ethylene biosynthesis and response between 100 and 160 DAP that reflect a transition from system 1 to system 2 ethylene production (Tranbarger et al., 2011). In contrast to tomato and banana, *EgACO* expression in the mesocarp is consistent with functions during the transition from system 1 to 2 (between 100 DAP and 140 DAP), and during ripening system 2 (between 140 and 160 DAP), but not during system 1 (100 DAP).

As in tomato, the increase in *EgACO* transcript during the transition from system 1 to 2 could be controlled either by ethylene independent or dependent regulation (Blume and Grierson, 1997). A comparison of *EgACO* and *EgPG4* transcript profiles during development and ripening and in response to
<figure><img src="image_8.png" /><figcaption>FIGURE 5 | qPCR analysis of NAC transcripts expressed in the mesocarp and induction in response to ethylene treatments. Asterisks indicate significant difference and increase or decrease in relative abundance by at least 2 fold (unpaired student's t-test *P* > 0.05) compared to the control at 0 h. Profile curves (right panels) are the corresponding *E. guineensis* gene locus expression during ripening at 100, 120, 140, and 160 DAP from RNA-Seq Illumina data.</figcaption></figure>
ethylene suggests different modes of regulation in the mesocarp. *EgPG4* is barely detectable at the end of the transition from system 1 to system 2 ethylene production (approximated at 140 DAP), then increases dramatically at the autocatalytic ethylene stage (160 DAP), suggesting *EgPG4* expression may require a higher ethylene concentration threshold and/or does not respond to the same signals as *EgACO* during the transition from system 1 to 2. In contrast, expression of the principle tomato fruit TFPG increases more synchronistically with the increase in ethylene production (Dellapenna et al., 1986). In addition, there are at least two PGs expressed during banana ripening, one of which is the most highly expressed during ripening, and appears to be induced by ethylene (Asif and Nath, 2005; Asif et al., 2014). However, the ethylene treatments were for days and not hours as in the current study, so it is difficult to compare the results. In oil palm, the fact that *EgPG4* is consistently more highly induced by ethylene than *EgACO* at all treatment time points and stages examined suggests *EgPG4* expression is controlled differently by ethylene than *EgACO*, and argues against a higher ethylene threshold requirement for *EgPG4* regulation. The expression of *EgACO* during the transition from system 1 to 2 suggests developmental ethylene independent signals could induce *EgACO* during the transition, and not ethylene alone. One example of this in tomato is the LeHB-1 homeobox protein tomato that binds to the *LeACO1* promoter, while the inhibition of LeHB-1 mRNA results in reduced *LeACO1* mRNA (Lin et al., 2008). A search of the oil palm mesocarp transcriptome for sequences similar to LeHB-1 resulted in sequences with very low e-values (1e-26 to 1e-29) and low read counts, and therefore no evidence for a similar module was observed. Overall, while gene expression of key genes during the transition from system 1 to 2 in the oil palm mesocarp shares some similarities to that of banana and tomato, the examples of *EgACO* and *EgPG4* provided here suggest differences.

## The Mesocarp Undergoes a Transcriptional Transition from Cell Wall Synthesis to Degradation during Ripening

The cell wall is an extremely complex structure made up of interlaced cellulose and hemicellulose polysaccharide polymers with structural glycoproteins embedded in a pectin matrix including HG, rhamnogalacturonan I and rhamnogalacturonan II. Unsurprisingly, many cell wall modifying proteins appear to function during ripening related cell wall disassembly and fruit softening. Attempts to silence single genes encoding only one type of modification have had little success to stop softening in the tomato fruit (Brummell, 2006; Seymour et al., 2013a). The identity of cell wall related transcripts identified in the current study that are differentially expressed during mesocarp ripening suggest a transcriptional based transition from synthesis to hydrolysis of the major cell wall polymers occurs in the mesocarp during ripening. Importantly, a single PG transcript (*EgPG4*) accounts for a large percentage of the total transcripts identified, while *EgPG4* expressed at a later stage of ripening, and is also induced by ethylene, similar to TFPG *PG2* transcript of tomato (Dellapenna et al.,

1987). In addition, a number of other transcripts encoding various cell wall modifying activities were observed to have similar profiles to *EgPG4*, including glucan endo-1,3-beta-glucosidase-like, lichenase-2-like, chitinase 2-like, xyloglucan endotransglucosylase/hydrolase protein 23, UDP-glucose 6-dehydrogenase 5-like proteins were identified (Figure 3). In contrast, there are also a number of cell wall related transcripts with expression profiles similar to the ethylene biosynthesis transcript *EgACO*, also induced by ethylene, including, beta-galactosidase-like, alpha-glucosidase, and expansin-A10-like. Indeed, while *EgACO* and *EgPG4* are both induced by ethylene and are expressed in a ripening dependent manner, their ethylene induction and ripening expression profiles are different. That *EgACO* and *EgPG4* share similar expression profiles with various cell wall modifying protein transcripts suggests ethylene may transcriptionally regulate cell wall activities differentially as proposed previously (Ireland et al., 2014). Previous work on cell wall related gene expression during oil palm mesocarp ripening monitored only very few genes (Teh et al., 2014). In contrast, the current work includes a more comprehensive analysis of 56 differentially expressed genes from at least 20 different cell wall related gene families. The results suggest a close relation between the transcriptional activation of genes encoding expansin, polygalacturonase, mannosidases, beta-galactosidase, and xyloglucan endotransglucosylase/hydrolase proteins in the ripening oil palm mesocarp, comparable to that observed during tomato fruit ripening (Seymour et al., 2013a). While there is less available data with monocot fruit species, expansin and xyloglucan endotransglucosylase/hydrolase protein genes were also highly expressed during banana ripening (Asif et al., 2014). Overall, these observations provide evidence for the conservation of a large number of cell wall related activities that may be regulated by ethylene and other fruit ripening specific regulatory factors in diverse lineages of both eudicots and monocots.

## Paralogous Members of the NAC Domain Family of Transcription Factors Are Regulated in the Mesocarp during Ripening and Induced by Ethylene

The NAC gene family is a very large group that encodes plant specific TFs involved in many processes including secondary cell wall development, stress response, and leaf senescence (Puranik et al., 2012; Kim et al., 2016; Monniaux and Hay, 2016). In addition, NAC TFs are upstream of ethylene related transcriptional regulation during both eudicot and monocot fruit ripening, and the NAC-NOR is a master regulator of tomato fruit ripening (Giovannoni, 2004; Shan et al., 2012; Seymour et al., 2013a; Kou et al., 2016). From a previous study we observed a number of NAC domain family TF genes that were differentially expressed during oil palm fruit ripening suggesting these TFs may function in response to and/or during the ethylene burst (Tranbarger et al., 2011). In the current genome wide study, we identified at least 31 NAC domain transcripts expressed in the oil palm mesocarp during the stages when ethylene production increases. Experiments with ethylene treated fruit revealed at least four have enhanced expression in the presence
Asif, M. H., and Nath, P. (2005). Expression of multiple forms of polygalacturonase gene during ripening in banana fruit. *Plant Physiol. Biochem.* 43, 177-184. doi: 10.1016/j.plaphy.2005.01.011

Bapat, V. A., Trivedi, P. K., Ghosh, A., Sane, V. A., Ganapathi, T. R., and Nath, P. (2010). Ripening of fleshy fruit: molecular insight and the role of ethylene. *Biotechnol. Adv.* 28, 94-107. doi: 10.1016/j.biotechnadv.2009.10.002

Barry, C. S., Blume, B., Bouzayen, M., Cooper, W., Hamilton, A. J., and Grierson, D. (1996). Differential expression of the 1-aminocyclopropane-1-carboxylate oxidase gene family of tomato. *Plant J.* 9, 525-535. doi: 10.1046/j.1365-313X.1996.09040525.x

Barry, C. S., and Giovannoni, J. J. (2007). Ethylene and fruit ripening. *J. Plant Growth Regul.* 26, 143-159. doi: 10.1007/s00344-007-9002-y

Bennett, A. B., and Labavitch, J. M. (2008). Ethylene and ripening-regulated expression and function of fruit cell wall modifying proteins. *Plant Sci.* 175, 130-136. doi: 10.1016/j.plantsci.2008.03.004

Blume, B., and Grierson, D. (1997). Expression of ACC oxidase promoters in tomato and Nicotiana plumbaginifolia regulated by developmental and environmental stimuli. *Plant J.* 12, 731-746. doi: 10.1046/j.1365-313X.1997.12040731.x

Brummell, D. A. (2006). Cell wall disassembly in ripening fruit. *Funct. Plant Biol.* 33, 103-119. doi: 10.1071/FP05234

Brummell, D. A., and Harpster, M. H. (2001). Cell wall metabolism in fruit softening and quality and its manipulation in transgenic plants. *Plant Mol. Biol.* 47, 311-340. doi: 10.1023/A:1010656104304

Buffard-Morel, J., Verdeil, J. L., and Pannetier, C. (1992). Embryogenèse somatique du cocotier (Cocos nucifera-L.) à partir d'explants foliarires: étude histologique. *Can. J. Bot.* 70, 735-741. doi: 10.1139/b92-094

Cara, B., and Giovannoni, J. J. (2008). Molecular biology of ethylene during tomato fruit development and maturation. *Plant Sci.* 175, 106-113. doi: 10.1016/j.plantsci.2008.03.021

Chanderbali, A. S., Albert, V. A., Ashworth, V. E., Clegg, M. T., Litz, R. E., Soltis, D. E., et al. (2008). *Persea americana* (avocado): bringing ancient flowers to fruit in the genomics era. *Bioessays* 30, 386-396. doi: 10.1002/bies.20721

Chanderbali, A. S., Albert, V. A., Leebens-Mack, J., Altman, N. S., Soltis, D. E., and Soltis, P. S. (2009). Transcriptional signatures of ancient floral developmental genetics in avocado (*Persea americana*; Lauraceae). *Proc. Natl. Acad. Sci. U.S.A.* 106, 8929-8934. doi: 10.1073/pnas.0811476106

Chervin, C., Tira-Umphon, A., Terrier, N., Zouine, M., Severac, D., and Roustan, J. P. (2008). Stimulation of the grape berry expansion by ethylene and effects on related gene transcripts, over the ripening phase. *Physiol. Plant.* 134, 534-546. doi: 10.1111/j.1399-3054.2008.01158.x

Cooley, M. B., and Yoder, J. I. (1998). Insertional inactivation of the tomato polygalacturonase gene. *Plant Mol. Biol.* 38, 521-530. doi: 10.1023/A:1006086004262

Dellapenna, D., Alexander, D. C., and Bennett, A. B. (1986). Molecular cloning of tomato fruit polygalacturonase: analysis of polygalacturonase mRNA levels during ripening. *Proc. Natl. Acad. Sci. U.S.A.* 83, 6420-6424. doi: 10.1073/pnas.83.17.6420

Dellapenna, D., Kates, D. S., and Bennett, A. B. (1987). Polygalacturonase gene expression in rutgers, rin, nor, and Nr tomato fruits. *Plant Physiol.* 85, 502-507. doi: 10.1104/pp.85.2.502

Dellapenna, D., Lincoln, J. E., Fischer, R. L., and Bennett, A. B. (1989). Transcriptional analysis of polygalacturonase and other ripening associated genes in rutgers, rin, nor, and Nr tomato fruit. *Plant Physiol.* 90, 1372-1377. doi: 10.1104/pp.90.4.1372

Dereeper, A., Guignon, V., Blanc, G., Audic, S., Buffet, S., Chevenet, F., et al. (2008). Phylogeny.fr: robust phylogenetic analysis for the non-specialist. *Nucleic Acids Res.* 36, W465-W469. doi: 10.1093/nar/gkn180

D'hont, A., Denoeud, F., Aury, J. M., Baurens, F. C., Carreel, F., Garsmeur, O., et al. (2012). The banana (*Musa acuminata*) genome and the evolution of monocotyledonous plants. *Nature* 488, 213-217. doi: 10.1038/nature11241

Elitzur, T., Vrebalov, J., Giovannoni, J. J., Goldschmidt, E. E., and Friedman, H. (2010). The regulation of MADS-box gene expression during ripening of banana and their regulatory interaction with ethylene. *J. Exp. Bot.* 61, 1523-1535. doi: 10.1093/jxb/erq017

Elitzur, T., Yakir, E., Quansah, L., Zhangjun, F., Vrebalov, J., Khayat, E., et al. (2016). Banana MaMADS transcription factors are necessary for fruit ripening

and molecular tools to promote shelf-life and food security. *Plant Physiol.* 171, 380-391. doi: 10.1104/pp.15.01866

Fooyontphanich, K., Morcillo, F., Amblard, P., Collin, M., Jantasuriyarat, C., Tangphatsornruang, S., et al. (2016). A phenotypic test for delay of abscission and non-abscission oil palm fruit and validation by abscission marker gene expression analysis. *Acta Hortic.* 1119, 97-104. doi: 10.17660/ActaHortic.2016.1119.13

Giovannoni, J. J. (2004). Genetic regulation of fruit development and ripening. *Plant Cell* 16, S170-S180. doi: 10.1105/tpc.019158

Götz, S., García-Gómez, J. M., Terol, J., Williams, T. D., Nagaraj, S. H., Nueda, M. J., et al. (2008). High-throughput functional annotation and data mining with the Blast2GO suite. *Nucleic Acids Res.* 36, 3420-3435. doi: 10.1093/nar/gkn176

Guerin, C., Joët, T., Serret, J., Lasermes, P., Vaissayre, V., Agbessi, M. D., et al. (2016). Gene coexpression network analysis of oil biosynthesis in an interspecific backcross of oil palm. *Plant J.* 87, 423-441. doi: 10.1111/tpj.13208

Holdsworth, M. J., Bird, C. R., Ray, J., Schuch, W., and Grierson, D. (1987). Structure and expression of an ethylene-related mRNA from tomato. *Nucleic Acids Res.* 15, 731-739. doi: 10.1093/nar/15.2.731

Ireland, H. S., Gunaseelan, K., Muddumage, R., Tacken, E. J., Putterill, J., Johnston, J. W., et al. (2014). Ethylene regulates Apple (Malus x domestica) fruit softening through a dose x time-dependent mechanism and through differential sensitivities and dependencies of cell wall-modifying genes. *Plant Cell Physiol.* 55, 1005-1016. doi: 10.1093/pcp/pcu034

Jafari, Z., Haddad, R., Hosseini, R., and Garoosi, G. (2013). Cloning, identification and expression analysis of ACC oxidase gene involved in ethylene production pathway. *Mol. Biol. Rep.* 40, 1341-1350. doi: 10.1007/s11033-012-2178-7

Jensen, M. K., Kjaersgaard, T., Nielsen, M. M., Galberg, P., Petersen, K., O'shea, C., et al. (2010). The *Arabidopsis thaliana* NAC transcription factor family: structure-function relationships and determinants of ANAC019 stress signalling. *Biochem. J.* 426, 183-196. doi: 10.1042/BJ20091234

Jin, J., He, K., Tang, X., Li, Z., Lv, L., Zhao, Y., et al. (2015). An Arabidopsis transcriptional regulatory map reveals distinct functional and evolutionary features of novel transcription factors. *Mol. Biol. Evol.* 32, 1767-1773. doi: 10.1093/molbev/msv058

Jin, J., Tian, F., Yang, D. C., Meng, Y. Q., Kong, L., Luo, J. C., et al. (2016). PlantTFDB 4.0: toward a central hub for transcription factors and regulatory interactions in plants. *Mol. Biol. Evol.* 32, 1767-1773. doi: 10.1093/molbev/msv058

Jourda, C., Cardi, C., Mbéguie, A. M. D., Bocs, S., Garsmeur, O., D'hont, A., et al. (2014). Expansion of banana (*Musa acuminata*) gene families involved in ethylene biosynthesis and signalling after lineage-specific whole-genome duplications. *New Phytol.* 202, 986-1000. doi: 10.1111/nph.12710

Kim, H. J., Nam, H. G., and Lim, P. O. (2016). Regulatory network of NAC transcription factors in leaf senescence. *Curr. Opin. Plant Biol.* 33, 48-56. doi: 10.1016/j.pbi.2016.06.002

Klee, H. J., and Giovannoni, J. J. (2011). Genetics and control of tomato fruit ripening and quality attributes. *Annu. Rev. Genet.* 45, 41-59. doi: 10.1146/annurev-genet-110410-132507

Klinghammer, M., and Tenhaken, R. (2007). Genome-wide analysis of the UDP-glucose dehydrogenase gene family in Arabidopsis, a key enzyme for matrix polysaccharides in cell walls. *J. Exp. Bot.* 58, 3609-3621. doi: 10.1093/jxb/erm209

Kou, X., Liu, C., Han, L., Wang, S., and Xue, Z. (2016). NAC transcription factors play an important role in ethylene biosynthesis, reception and signaling of tomato fruit ripening. *Mol. Genet. Genomics* 291, 1205-1217. doi: 10.1007/s00438-016-1177-0

Larkin, M. A., Blackshields, G., Brown, N. P., Chenna, R., Mcgettigan, P. A., McWilliam, H., et al. (2007). Clustal W and clustal X version 2.0. *Bioinformatics* 23, 2947-2948. doi: 10.1093/bioinformatics/btm404

Lelievre, J. M., Latche, A., Jones, B., Bouzayen, M., and Pech, J. C. (1997). Ethylene and fruit ripening. *Physiol. Plant.* 101, 727-739. doi: 10.1111/j.1399-3054.1997.tb01057.x

Lin, Z. F., Hong, Y. G., Yin, M. G., Li, C. Y., Zhang, K., and Grierson, D. (2008). A tomato HD-Zip homeobox protein, LeHB-1, plays an important role in floral organogenesis and ripening. *Plant J.* 55, 301-310. doi: 10.1111/j.1365-313X.2008.03505.x

Liu, M. C., Pirrello, J., Chervin, C., Roustan, J. P., and Bouzayen, M. (2015). Ethylene control of fruit ripening: revisiting the complex
network of transcriptional regulation. *Plant Physiol.* 169, 2380-2390.
doi: 10.1104/pp.15.01361

Liu, X., Shiomi, S., Nakatsuka, A., Kubo, Y., Nakamura, R., and Inaba, A. (1999).
Characterization of ethylene biosynthesis associated with ripening in banana
fruit. *Plant Physiol.* 121, 1257-1266. doi: 10.1104/pp.121.4.1257

Magallón, S., Gómez-Acevedo, S., Sánchez-Reyes, L. L., and Hernández-Hernández, T. (2015). A metacalibrated time-tree documents the early rise of flowering plant phylogenetic diversity. *New Phytol.* 207, 437-453.
doi: 10.1111/nph.13264

Martel, C., Vrebalov, J., Tafelmeyer, P., and Giovannoni, J. J. (2011). The tomato MADS-box transcription factor RIPENING INHIBITOR interacts with promoters involved in numerous ripening processes in a COLORLESS NONRIPENING-dependent manner. *Plant Physiol.* 157, 1568-1579.
doi: 10.1104/pp.111.181107

Monniaux, M., and Hay, A. (2016). Cells, walls, and endless forms. *Curr. Opin. Plant Biol.* 34, 114-121. doi: 10.1016/j.pbi.2016.10.010

Morcillo, F., Gagneur, C., Adam, H., Richaud, F., Singh, R., Cheah, S. C., et al. (2006). Somaclonal variation in micropropagated oil palm. Characterization of two novel genes with enhanced expression in epigenetically abnormal cell lines and in response to auxin. *Tree Physiol.* 26, 585-594.
doi: 10.1093/treephys/26.5.585

Nakatsuka, A., Murachi, S., Okunishi, H., Shiomi, S., Nakano, R., Kubo, Y., et al. (1998). Differential expression and internal feedback regulation of 1-aminocyclopropane-1-carboxylate synthase, 1-aminocyclopropane-1-carboxylate oxidase, and ethylene receptor genes in tomato fruit during development and ripening. *Plant Physiol.* 118, 1295-1305.
doi: 10.1104/pp.118.4.1295

Nualwijit, N., and Lerslerwong, L. (2014). Post harvest ripening of oil palm fruit is accelerated by application of exogenous ethylene. *Songklanakarin J. Sci. Technol.* 36, 255-259.

Ooka, H., Satoh, K., Doi, K., Nagata, T., Otomo, Y., Murakami, K., et al. (2003).
Comprehensive analysis of NAC family genes in *Oryza sativa* and *Arabidopsis
thaliana*. *DNA Res.* 10, 239-247. doi: 10.1093/dnares/10.6.239

Osorio, S., Alba, R., Damasceno, C. M., Lopez-Casado, G., Lohse, M., Zanor, M. I., et al. (2011). Systems biology of tomato fruit development: combined transcript, protein, and metabolite analysis of tomato transcription factor (nor, rin) and ethylene receptor (Nr) mutants reveals novel regulatory interactions. *Plant Physiol.* 157, 405-425. doi: 10.1104/pp.111.175463

Osorio, S., Alba, R., Nikoloski, Z., Kochevenko, A., Fernie, A. R., and Giovannoni, J. J. (2012). Integrative comparative analyses of transcript and metabolite profiles from pepper and tomato ripening and development stages uncovers species-specific patterns of network regulatory behaviour. *Plant Physiol.* 159, 1713-1729. doi: 10.1104/pp.112.199711

Pfaffl, M. W. (2001). A new mathematical model for relative quantification in real-time RT-PCR. *Nucleic Acids Res.* 29:e45. doi: 10.1093/nar/29.9.e45

Prasanna, V., Prabha, T. N., and Tharanathan, R. N. (2007). Fruit ripening phenomena-an overview. *Crit. Rev. Food Sci.* 47, 1-19.
doi: 10.1080/10408390600976841

Puranik, S., Sahu, P. P., Srivastava, P. S., and Prasad, M. (2012). NAC proteins: regulation and role in stress tolerance. *Trends Plant Sci.* 17, 369-381.
doi: 10.1016/j.tplants.2012.02.004

Roongsattham, P., Morcillo, F., Fooyontphanich, K., Jantasuriyat, C., Tragoonrung, S., Amblard, P., et al. (2016). Cellular and pectin dynamics during abscission zone development and ripe fruit abscission of the monocot oil palm. *Front. Plant Sci.* 7:540. doi: 10.3389/fpls.2016.00540

Roongsattham, P., Morcillo, F., Jantasuriyat, C., Pizot, M., Moussu, S., Jayaweera, D., et al. (2012). Temporal and spatial expression of polygalacturonase gene family members reveals divergent regulation during fleshy fruit ripening and abscission in the monocot species oil palm. *BMC Plant Biol.* 12:150.
doi: 10.1186/1471-2229-12-150

Seymour, G. B., Chapman, N. H., Chew, B. L., and Rose, J. K. (2013a). Regulation of ripening and opportunities for control in tomato and other fruits. *Plant Biotechnol. J.* 11, 269-278. doi: 10.1111/j.1467-7652.2012.00738.x

Seymour, G. B., Ostergaard, L., Chapman, N. H., Knapp, S., and Martin, C. (2013b). Fruit development and ripening. *Annu. Rev. Plant Biol.* 64, 219-241.
doi: 10.1146/annurev-arplant-050312-120057

Shan, W., Kuang, J. F., Chen, L., Xie, H., Peng, H. H., Xiao, Y. Y., et al. (2012).
Molecular characterization of banana NAC transcription factors and their
interactions with ethylene signalling component EIL during fruit ripening. *J.
Exp. Bot.* 63, 5171-5187. doi: 10.1093/jxb/ers178

Sheehy, R. E., Kramer, M., and Hiatt, W. R. (1988). Reduction of polygalacturonase activity in tomato fruit by antisense RNA. *Proc. Natl. Acad. Sci. U.S.A.* 85, 8805-8809. doi: 10.1073/pnas.85.23.8805

Singh, R., Ong-Abdullah, M., Low, E. T., Manaf, M. A., Rosli, R., Sambanthamurthi, R., et al. (2013). Oil palm genome sequence reveals divergence of interfertile species in old and new worlds. *Nature* 500, 335-339.
doi: 10.1038/nature12309

Smith, C. J. S., Watson, C. F., Ray, J., Bird, C. R., Morris, P. C., Schuch, W., et al. (1988). Antisense RNA inhibition of polygalacturonase gene-expression in transgenic tomatoes. *Nature* 334, 724-726. doi: 10.1038/334 724a0

Smith, C. J., Watson, C. F., Morris, P. C., Bird, C. R., Seymour, G. B., Gray, J. E., et al. (1990). Inheritance and effect on ripening of antisense polygalacturonase genes in transgenic tomatoes. *Plant Mol. Biol.* 14, 369-379.
doi: 10.1007/BF00028773

Teh, H. F., Neoh, B. K., Wong, Y. C., Kwong, Q. B., Ooi, T. E. K., Ng, T. L., et al. (2014). Hormones, polyamines, and cell wall metabolism during oil palm fruit mesocarp development and ripening. *J. Agric. Food Chem.* 62, 8143-8152.
doi: 10.1021/jf500975h

Trainotti, L., Pavanello, A., and Casadoro, G. (2005). Different ethylene receptors show an increased expression during the ripening of strawberries: does such an increment imply a role for ethylene in the ripening of these non-climacteric fruits? *J. Exp. Bot.* 56, 2037-2046. doi: 10.1093/jxb/e ri202

Tranbarger, T. J., Dussert, S., Joet, T., Argout, X., Summo, M., Champion, A., et al. (2011). Regulatory mechanisms underlying oil palm fruit mesocarp maturation, ripening, and functional specialization in lipid and carotenoid metabolism. *Plant Physiol.* 156, 564-584. doi: 10.1104/pp.111.175141

Yokotani, N., Nakano, R., Imanishi, S., Nagata, M., Inaba, A., and Kubo, Y. (2009). Ripening-associated ethylene biosynthesis in tomato fruit is autocatalytically and developmentally regulated. *J. Exp. Bot.* 60, 3433-3442.
doi: 10.1093/jxb/erp185

Zdobnov, E. M., and Apweiler, R. (2001). InterProScan-an integration platform for the signature-recognition methods in interpro. *Bioinformatics* 17, 847-848.
doi: 10.1093/bioinformatics/17.9.847

Zheng, Y., Jiao, C., Sun, H., Rosli, H. G., Pombo, M. A., Zhang, P., et al. (2016). iTAK: a program for genome-wide prediction and classification of plant transcription factors, transcriptional regulators, and protein kinases. *Mol. Plant* 9, 1667-1670. doi: 10.1016/j.molp.2016.09.014

Zhu, M., Chen, G., Zhou, S., Tu, Y., Wang, Y., Dong, T., et al. (2014). A new tomato NAC (NAM/ATAF1/2/CUC2) transcription factor, SINAC4, functions as a positive regulator of fruit ripening and carotenoid accumulation. *Plant Cell Physiol.* 55, 119-135. doi: 10.1093/pcp/pct162

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2017 Tranbarger, Fooyontphanich, Roongsattham, Pizot, Collin, Jantasuriyat, Suraninpong, Tragoonrung, Dussert, Verdeil and Morcillo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.