Paralogue Interference Affects the Dynamics after Gene Duplication

Homospermidine synthase evolution and the origin(s) of pyrrolizidine alkaloids in Apocynaceae

PREMISE

Enzymes that are encoded by paralogous genes and produce identical specialized metabolites in distantly related plant lineages are strong evidence of parallel phenotypic evolution. Inference of phenotypic homology for metabolites produced by orthologous genes is less straightforward, since orthologs may be recruited in parallel into novel pathways. In prior research on pyrrolizidine alkaloids (PAs), specialized metabolites of Apocynaceae, the evolution of homospermidine synthase (HSS), an enzyme of PA biosynthesis, was reconstructed and a single origin of PAs inferred because HSS enzymes of all known PA-producing Apocynaceae species are orthologous and descended from an ancestral enzyme with the motif (VXXXD) of an optimized HSS.

METHODS

We increased sampling, tested the effect of amino acid motif on HSS function, revisited motif evolution, and tested for selection to infer evolution of HSS function and its correlation with phenotype.

RESULTS

Some evidence supports a single origin of PAs: an IXXXD HSS-like gene, similar in function to VXXXD HSS, evolved in the shared ancestor of all PA-producing species; loss of HSS function occurred multiple times via pseudogenization and perhaps via evolution of an IXXXN motif. Other evidence indicates multiple origins: the VXXXD motif, highly correlated with the PA phenotype, evolved two or four times independently; the ancestral IXXXD gene was not under positive selection, while some VXXXD genes were; and substitutions at sites experiencing positive selection occurred on multiple branches in the HSS-like gene tree.

CONCLUSIONS

The complexity of the genotype-function-phenotype map confounds the inference of PA homology from HSS-like gene evolution in Apocynaceae.

An NLR paralog Pit2 generated from tandem duplication of Pit1 fine-tunes Pit1 localization and function

NLR family proteins act as intracellular receptors. Gene duplication amplifies the number of NLR genes, and subsequent mutations occasionally provide modifications to the second gene that benefits immunity. However, evolutionary processes after gene duplication and functional relationships between duplicated NLRs remain largely unclear. Here, we report that the rice NLR protein Pit1 is associated with its paralogue Pit2. The two are required for the resistance to rice blast fungus but have different functions: Pit1 induces cell death, while Pit2 competitively suppresses Pit1-mediated cell death. During evolution, the suppression of Pit1 by Pit2 was probably generated through positive selection on two fate-determining residues in the NB-ARC domain of Pit2, which account for functional differences between Pit1 and Pit2. Consequently, Pit2 lost its plasma membrane localization but acquired a new function to interfere with Pit1 in the cytosol. These findings illuminate the evolutionary trajectory of tandemly duplicated NLR genes after gene duplication.

Mutational biases favor complexity increases in protein interaction networks after gene duplication

Biological systems can gain complexity over time. While some of these transitions are likely driven by natural selection, the extent to which they occur without providing an adaptive benefit is unknown. At the molecular level, one example is heteromeric complexes replacing homomeric ones following gene duplication. Here, we build a biophysical model and simulate the evolution of homodimers and heterodimers following gene duplication using distributions of mutational effects inferred from available protein structures. We keep the specific activity of each dimer identical, so their concentrations drift neutrally without new functions. We show that for more than 60% of tested dimer structures, the relative concentration of the heteromer increases over time due to mutational biases that favor the heterodimer. However, allowing mutational effects on synthesis rates and differences in the specific activity of homo- and heterodimers can limit or reverse the observed bias toward heterodimers. Our results show that the accumulation of more complex protein quaternary structures is likely under neutral evolution, and that natural selection would be needed to reverse this tendency.

Subgenome evolutionary dynamics in allotetraploid ferns: insights from the gene expression patterns in the allotetraploid species Phegopteris decursivepinnata (Thelypteridacea, Polypodiales)

Allopolyploidization often leads to disruptive conflicts among more than two sets of subgenomes, leading to genomic modifications and changes in gene expression. Although the evolutionary trajectories of subgenomes in allopolyploids have been studied intensely in angiosperms, the dynamics of subgenome evolution remain poorly understood in ferns, despite the prevalence of allopolyploidization. In this study, we have focused on an allotetraploid fern-Phegopteris decursivepinnata-and its diploid parental species, P. koreana (K) and P. taiwaniana (T). Using RNA-seq analyses, we have compared the gene expression profiles for 9,540 genes among parental species, synthetic F1 hybrids, and natural allotetraploids. The changes in gene expression patterns were traced from the F1 hybrids to the natural allopolyploids. This study has revealed that the expression patterns observed in most genes in the F1 hybrids are largely conserved in the allopolyploids; however, there were substantial differences in certain genes between these groups. In the allopolyploids compared with the F1 hybrids, the number of genes showing a transgressive pattern in total expression levels was increased. There was a slight reduction in T-dominance and a slight increase in K-dominance, in terms of expression level dominance. Interestingly, there is no obvious bias toward the T- or K-subgenomes in the number and expression levels overall, showing the absence of subgenome dominance. These findings demonstrated the impacts of the substantial transcriptome change after hybridization and the moderate modification during allopolyploid establishment on gene expression in ferns and provided important insights into subgenome evolution in polyploid ferns.

Unusual Evolution of Cephalopod Tryptophan Indole-Lyases

Tryptophan indole-lyase (TIL), a pyridoxal-5-phosphate-dependent enzyme, catalyzes the hydrolysis of L-tryptophan (L-Trp) to indole and ammonium pyruvate. TIL is widely distributed among bacteria and bacterial TILs consist of a D2-symmetric homotetramer. On the other hand, TIL genes are also present in several metazoans. Cephalopods have two TILs, TILα and TILβ, which are believed to be derived from a gene duplication that occurred before octopus and squid diverged. However, both TILα and TILβ individually contain disruptive amino acid substitutions for TIL activity, and neither was active when expressed alone. When TILα and TILβ were coexpressed, however, they formed a heterotetramer that exhibited low TIL activity. The loss of TIL activity of the heterotetramer following site-directed mutagenesis strongly suggests that the active heterotetramer contains the TILα/TILβ heterodimer. Metazoan TILs generally have lower kcat values for L-Trp than those of bacterial TILs, but such low TIL activity may be rather suitable for metazoan physiology, where L-Trp is in high demand. Therefore, reduced activity may have been a less likely target for purifying selection in the evolution of cephalopod TILs. Meanwhile, the unusual evolution of cephalopod TILs may indicate the difficulty of post-gene duplication evolution of enzymes with catalytic sites contributed by multiple subunits, such as TIL.

How gene duplication diversifies the landscape of protein oligomeric state and function

Oligomeric proteins are central to cellular life and the duplication and divergence of their genes is a key driver of evolutionary innovations. The duplication of a gene coding for an oligomeric protein has numerous possible outcomes, which motivates questions on the relationship between structural and functional divergence. How do protein oligomeric states diversify after gene duplication? In the simple case of duplication of a homo-oligomeric protein gene, what properties can influence the fate of descendant paralogs toward forming independent homomers or maintaining their interaction as a complex? Furthermore, how are functional innovations associated with the diversification of oligomeric states? Here, we review recent literature and present specific examples in an attempt to illustrate and answer these questions.

The long road of functional recruitment-The evolution of a gene duplicate to pyrrolizidine alkaloid biosynthesis in the morning glories (Convolvulaceae)

In plants, homospermidine synthase (HSS) is a pathway-specific enzyme initiating the biosynthesis of pyrrolizidine alkaloids (PAs), which function as a chemical defense against herbivores. In PA-producing Convolvulaceae ("morning glories"), HSS originated from deoxyhypusine synthase at least >50 to 75 million years ago via a gene duplication event and subsequent functional diversification. To study the recruitment of this ancient gene duplicate to PA biosynthesis, the presence of putative hss gene copies in 11 Convolvulaceae species was analyzed. Additionally, various plant parts from seven of these species were screened for the presence of PAs. Although all of these species possess a putative hss copy, PAs could only be detected in roots of Ipomoea neei (Spreng.) O'Donell and Distimake quinquefolius (L.) A.R.Simões & Staples in this study. A precursor of PAs was detected in roots of Ipomoea alba L. Thus, despite sharing high sequence identities, the presence of an hss gene copy does not correlate with PA accumulation in particular species of Convolvulaceae. In vitro activity assays of the encoded enzymes revealed a broad spectrum of enzyme activity, further emphasizing a functional diversity of the hss gene copies. A recently identified HSS specific amino acid motif seems to be important for the loss of the ancestral protein function-the activation of the eukaryotic initiation factor 5A (eIF5A). Thus, the motif might be indicative for a change of function but allows not to predict the new function. This emphasizes the challenges in annotating functions for duplicates, even for duplicates from closely related species.

Evolution of binding preferences among whole-genome duplicated transcription factors

Throughout evolution, new transcription factors (TFs) emerge by gene duplication, promoting growth and rewiring of transcriptional networks. How TF duplicates diverge was studied in a few cases only. To provide a genome-scale view, we considered the set of budding yeast TFs classified as whole-genome duplication (WGD)-retained paralogs (~35% of all specific TFs). Using high-resolution profiling, we find that ~60% of paralogs evolved differential binding preferences. We show that this divergence results primarily from variations outside the DNA-binding domains (DBDs), while DBD preferences remain largely conserved. Analysis of non-WGD orthologs revealed uneven splitting of ancestral preferences between duplicates, and the preferential acquiring of new targets by the least conserved paralog (biased neo/sub-functionalization). Interactions between paralogs were rare, and, when present, occurred through weak competition for DNA-binding or dependency between dimer-forming paralogs. We discuss the implications of our findings for the evolutionary design of transcriptional networks.

A natural variant in ANP32B impairs influenza virus replication in human cells

Viruses require host factors to support their replication, and genetic variation in such factors can affect susceptibility to infectious disease. Influenza virus replication in human cells relies on ANP32 proteins, which are involved in assembly of replication-competent dimeric influenza virus polymerase (FluPol) complexes. Here, we investigate naturally occurring single nucleotide variants (SNV) in the human Anp32A and Anp32B genes. We note that variant rs182096718 in Anp32B is found at a higher frequency than other variants in either gene. This SNV results in a D130A substitution in ANP32B, which is less able to support FluPol activity than wild-type ANP32B and binds FluPol with lower affinity. Interestingly, ANP32B-D130A exerts a dominant negative effect over wild-type ANP32B and interferes with the functionally redundant paralogue ANP32A. FluPol activity and virus replication are attenuated in CRISPR-edited cells expressing wild-type ANP32A and mutant ANP32B-D130A. We propose a model in which the D130A mutation impairs FluPol dimer formation, thus resulting in compromised replication. We suggest that both homozygous and heterozygous carriers of rs182096718 may have some genetic protection against influenza viruses.

The integrative biology of genetic dominance

Dominance is a basic property of inheritance systems describing the link between a diploid genotype at a single locus and the resulting phenotype. Models for the evolution of dominance have long been framed as an opposition between the irreconcilable views of Fisher in 1928 supporting the role of largely elusive dominance modifiers and Wright in 1929, who viewed dominance as an emerging property of the structure of enzymatic pathways. Recent theoretical and empirical advances however suggest that these opposing views can be reconciled, notably using models investigating the regulation of gene expression and developmental processes. In this more comprehensive framework, phenotypic dominance emerges from departures from linearity between any levels of integration in the genotype‐to‐phenotype map. Here, we review how these different models illuminate the emergence and evolution of dominance. We then detail recent empirical studies shedding new light on the diversity of molecular and physiological mechanisms underlying dominance and its evolution. By reconciling population genetics and functional biology, we hope our review will facilitate cross‐talk among research fields in the integrative study of dominance evolution.

Development of an activity assay for characterizing deoxyhypusine synthase and its diverse reaction products

Deoxyhypusine synthase transfers an aminobutyl moiety from spermidine to the eukaryotic translation initiation factor 5A (eIF5A) in the first step of eIF5A activation. This exclusive post-translational modification is conserved in all eukaryotes. Activated eIF5A has been shown to be essential for cell proliferation and viability. Recent reports have linked the activation of eIF5A to several human diseases. Deoxyhypusine synthase, which is encoded by a single gene copy in most eukaryotes, was duplicated in several plant lineages during evolution, the copies being repeatedly recruited to pyrrolizidine alkaloid biosynthesis. However, the function of many of these duplicates is unknown. Notably, deoxyhypusine synthase is highly promiscuous and can catalyze various reactions, often of unknown biological relevance. To facilitate in-depth biochemical studies of this enzyme, we report here the development of a simple and robust in vitro enzyme assay. It involves precolumn derivatization of the polyamines taking part in the reaction and avoids the need for the previously used radioactively labeled tracers. The derivatized polyamines are quantified after high-performance liquid chromatography coupled to diode array and fluorescence detectors. By performing kinetic analyses of deoxyhypusine synthase and its paralog from the pyrrolizidine alkaloid-producing plant Senecio vernalis, we demonstrate that the assay unequivocally differentiates the paralogous enzymes. Furthermore, it detects and quantifies, in a single assay, the side reactions that occur in parallel to the main reaction. The presented assay thus provides a detailed biochemical characterization of deoxyhypusine synthase and its paralogs.

A hydrophobic ratchet entrenches molecular complexes

Most proteins assemble into multisubunit complexes1. The persistence of these complexes across evolutionary time is usually explained as the result of natural selection for functional properties that depend on multimerization, such as intersubunit allostery or the capacity to do mechanical work2. In many complexes, however, multimerization does not enable any known function3. An alternative explanation is that multimers could become entrenched if substitutions accumulate that are neutral in multimers but deleterious in monomers; purifying selection would then prevent reversion to the unassembled form, even if assembly per se does not enhance biological function3-7. Here we show that a hydrophobic mutational ratchet systematically entrenches molecular complexes. By applying ancestral protein reconstruction and biochemical assays to the evolution of steroid hormone receptors, we show that an ancient hydrophobic interface, conserved for hundreds of millions of years, is entrenched because exposure of this interface to solvent reduces protein stability and causes aggregation, even though the interface makes no detectable contribution to function. Using structural bioinformatics, we show that a universal mutational propensity drives sites that are buried in multimeric interfaces to accumulate hydrophobic substitutions to levels that are not tolerated in monomers. In a database of hundreds of families of multimers, most show signatures of long-term hydrophobic entrenchment. It is therefore likely that many protein complexes persist because a simple ratchet-like mechanism entrenches them across evolutionary time, even when they are functionally gratuitous.

The Evolution of Gene Duplicates in Angiosperms and the Impact of Protein-Protein Interactions and the Mechanism of Duplication

Gene duplicates, generated through either whole genome duplication (WGD) or small-scale duplication (SSD), are prominent in angiosperms and are believed to play an important role in adaptation and in generating evolutionary novelty. Previous studies reported contrasting evolutionary and functional dynamics of duplicate genes depending on the mechanism of origin, a behavior that is hypothesized to stem from constraints to maintain the relative dosage balance between the genes concerned and their interaction context. However, the mechanisms ultimately influencing loss and retention of gene duplicates over evolutionary time are not yet fully elucidated. Here, by using a robust classification of gene duplicates in Arabidopsis thaliana, Solanum lycopersicum, and Zea mays, large RNAseq expression compendia and an extensive protein-protein interaction (PPI) network from Arabidopsis, we investigated the impact of PPIs on the differential evolutionary and functional fate of WGD and SSD duplicates. In all three species, retained WGD duplicates show stronger constraints to diverge at the sequence and expression level than SSD ones, a pattern that is also observed for shared PPI partners between Arabidopsis duplicates. PPIs are preferentially distributed among WGD duplicates and specific functional categories. Furthermore, duplicates with PPIs tend to be under stronger constraints to evolve than their counterparts without PPIs regardless of their mechanism of origin. Our results support dosage balance constraint as a specific property of genes involved in biological interactions, including physical PPIs, and suggest that additional factors may be differently influencing the evolution of genes following duplication, depending on the species, time, and mechanism of origin.

A pair of E3 ubiquitin ligases compete to regulate filopodial dynamics and axon guidance

Appropriate axon guidance is necessary to form accurate neuronal connections. Axon guidance cues that stimulate cytoskeletal reorganization within the growth cone direct axon navigation. Filopodia at the growth cone periphery have long been considered sensors for axon guidance cues, yet how they respond to extracellular cues remains ill defined. Our previous work found that the filopodial actin polymerase VASP and consequently filopodial stability are negatively regulated via nondegradative TRIM9-dependent ubiquitination. Appropriate VASP ubiquitination and deubiquitination are required for axon turning in response to the guidance cue netrin-1. Here we show that the TRIM9-related protein TRIM67 outcompetes TRIM9 for interacting with VASP and antagonizes TRIM9-dependent VASP ubiquitination. The surprising antagonistic roles of two closely related E3 ubiquitin ligases are required for netrin-1-dependent filopodial responses, axon turning and branching, and fiber tract formation. We suggest a novel model in which coordinated regulation of VASP ubiquitination by a pair of interfering ligases is a critical element of VASP dynamics, filopodial stability, and axon guidance.

Paralog dependency indirectly affects the robustness of human cells

The protective redundancy of paralogous genes partly relies on the fact that they carry their functions independently. However, a significant fraction of paralogous proteins may form functionally dependent pairs, for instance, through heteromerization. As a consequence, one could expect these heteromeric paralogs to be less protective against deleterious mutations. To test this hypothesis, we examined the robustness landscape of gene loss-of-function by CRISPR-Cas9 in more than 450 human cell lines. This landscape shows regions of greater deleteriousness to gene inactivation as a function of key paralog properties. Heteromeric paralogs are more likely to occupy such regions owing to their high expression and large number of protein-protein interaction partners. Further investigation revealed that heteromers may also be under stricter dosage balance, which may also contribute to the higher deleteriousness upon gene inactivation. Finally, we suggest that physical dependency may contribute to the deleteriousness upon loss-of-function as revealed by the correlation between the strength of interactions between paralogs and their higher deleteriousness upon loss of function.

The role of structural pleiotropy and regulatory evolution in the retention of heteromers of paralogs

Gene duplication is a driver of the evolution of new functions. The duplication of genes encoding homomeric proteins leads to the formation of homomers and heteromers of paralogs, creating new complexes after a single duplication event. The loss of these heteromers may be required for the two paralogs to evolve independent functions. Using yeast as a model, we find that heteromerization is frequent among duplicated homomers and correlates with functional similarity between paralogs. Using in silico evolution, we show that for homomers and heteromers sharing binding interfaces, mutations in one paralog can have structural pleiotropic effects on both interactions, resulting in highly correlated responses of the complexes to selection. Therefore, heteromerization could be preserved indirectly due to selection for the maintenance of homomers, thus slowing down functional divergence between paralogs. We suggest that paralogs can overcome the obstacle of structural pleiotropy by regulatory evolution at the transcriptional and post-translational levels.

The effects of repeated whole genome duplication events on the evolution of cytokinin signaling pathway

BACKGROUND

It is thought that after whole-genome duplications (WGDs), a large fraction of the duplicated gene copies is lost over time while few duplicates are retained. Which factors promote survival or death of a duplicate remains unclear and the underlying mechanisms are poorly understood. According to the model of gene dosage balance, genes encoding interacting proteins are predicted to be preferentially co-retained after WGDs. Among these are genes encoding proteins involved in complexes or in signal transduction.

RESULTS

We have investigated the way that repeated WGDs during land plant evolution have affected cytokinin signaling to study patterns of gene duplicability and co-retention in this important signal transduction pathway. Through the integration of phylogenetic analyses with comparisons of genome collinearity, we have found that signal input mediated by cytokinin receptors proved to be highly conserved over long evolutionary time-scales, with receptors showing predominantly gene loss after repeated WGDs. However, the downstream elements, e,g. response regulators, were mainly retained after WGDs and thereby formed gene families in most plant lineages.

CONCLUSIONS

Gene dosage balance between the interacting components indicated by co-retention after WGDs seems to play a minor role in the evolution of cytokinin signaling pathway. Overall, core genes of cytokinin signaling show a highly heterogeneous pattern of gene retention after WGD, reflecting complex relationships between the various factors that shape the long-term fate of a duplicated gene.

Structural principles that enable oligomeric small heat-shock protein paralogs to evolve distinct functions

Oligomeric proteins assemble with exceptional selectivity, even in the presence of closely related proteins, to perform their cellular roles. We show that most proteins related by gene duplication of an oligomeric ancestor have evolved to avoid hetero-oligomerization and that this correlates with their acquisition of distinct functions. We report how coassembly is avoided by two oligomeric small heat-shock protein paralogs. A hierarchy of assembly, involving intermediates that are populated only fleetingly at equilibrium, ensures selective oligomerization. Conformational flexibility at noninterfacial regions in the monomers prevents coassembly, allowing interfaces to remain largely conserved. Homomeric oligomers must overcome the entropic benefit of coassembly and, accordingly, homomeric paralogs comprise fewer subunits than homomers that have no paralogs.

Phylogenomic analysis demonstrates a pattern of rare and long-lasting concerted evolution in prokaryotes

Concerted evolution, where paralogs in the same species show higher sequence similarity to each other than to orthologs in other species, is widely found in many species. However, cases of concerted evolution that last for hundreds of millions of years are very rare. By genome-wide analysis of a broad selection of prokaryotes, we provide strong evidence of recurrent concerted evolution in 26 genes, most of which have lasted more than ~500 million years. We find that most concertedly evolving genes are key members of important pathways, and encode proteins from the same complexes and/or pathways, suggesting coevolution of genes via concerted evolution to maintain gene balance. We also present LRCE-DB, a comprehensive online repository of long-lasting concerted evolution. Collectively, our study reveals that although most duplicated genes may diverge in sequence over a long period, on rare occasions this constraint can be breached, leading to unexpected long-lasting concerted evolution in a recurrent manner.

Sishuo Wang and Youhua Chen present an analysis of concerted evolution in prokaryotes using a new computational pipeline, iSeeCE. They find evidence in 26 genes for recurrent concerted evolution, most of which last more than ~500 million years, and provide a database, LRCE-DB, for data exploration.

Biochemical characterization of metabolism-based atrazine resistance in Amaranthus tuberculatus and identification of an expressed GST associated with resistance

Rapid detoxification of atrazine in naturally tolerant crops such as maize (Zea mays) and grain sorghum (Sorghum bicolor) results from glutathione S-transferase (GST) activity. In previous research, two atrazine-resistant waterhemp (Amaranthus tuberculatus) populations from Illinois, U.S.A. (designated ACR and MCR), displayed rapid formation of atrazine-glutathione (GSH) conjugates, implicating elevated rates of metabolism as the resistance mechanism. Our main objective was to utilize protein purification combined with qualitative proteomics to investigate the hypothesis that enhanced atrazine detoxification, catalysed by distinct GSTs, confers resistance in ACR and MCR. Additionally, candidate AtuGST expression was analysed in an F₂ population segregating for atrazine resistance. ACR and MCR showed higher specific activities towards atrazine in partially purified ammonium sulphate and GSH affinity-purified fractions compared to an atrazine-sensitive population (WCS). One-dimensional electrophoresis of these fractions displayed an approximate 26-kDa band, typical of GST subunits. Several phi- and tau-class GSTs were identified by LC-MS/MS from each population, based on peptide similarity with GSTs from Arabidopsis. Elevated constitutive expression of one phi-class GST, named AtuGSTF2, correlated strongly with atrazine resistance in ACR and MCR and segregating F₂ population. These results indicate that AtuGSTF2 may be linked to a metabolic mechanism that confers atrazine resistance in ACR and MCR.

MADS-domain transcription factors and the floral quartet model of flower development: linking plant development and evolution

The floral quartet model of floral organ specification poses that different tetramers of MIKC-type MADS-domain transcription factors control gene expression and hence the identity of floral organs during development. Here, we provide a brief history of the floral quartet model and review several lines of recent evidence that support the model. We also describe how the model has been used in contemporary developmental and evolutionary biology to shed light on enigmatic topics such as the origin of land and flowering plants. Finally, we suggest a novel hypothesis describing how floral quartet-like complexes may interact with chromatin during target gene activation and repression.

Characterization of DNA methyltransferase and demethylase genes in Fragaria vesca

DNA methylation is an epigenetic modification essential for gene regulations in plants, but understanding on how it is involved in fruit development, especially in non-climacteric fleshy fruit, is limited. The diploid woodland strawberry (Fragaria vesca) is an important model for non-climacteric fruit crops. In this study, we identified DNA methyltransferase genes and demethylase genes in Fragaria vesca and other angiosperm species. In accordance with previous studies, our phylogenetic analyses of those DNA methylation modifiers support the clustering of those genes into several classes. Our data indicate that whole-genome duplications and tandem duplications contributed to the expansion of those DNA methylation modifiers in angiosperms. We have further demonstrated that some DNA methylase and demethylase genes reach their highest expression levels in strawberry fleshy fruits when turning from white to red, suggesting that DNA methylation might undergo a dramatic change at the onset of fleshy fruit-ripening process. In addition, we have observed that expression of some DNA demethylase genes increases in response to various abiotic stresses including heat, cold, drought and salinity. Collectively, our study indicates a regulatory role of DNA methylation in the turning stage of non-climacteric fleshy fruit and responses to environment stimuli, and would facilitate functional studies of DNA methylation in the growth and development of non-climacteric fruits.