The chromatin remodeller MdRAD5B enhances drought tolerance by coupling MdLHP1-mediated H3K27me3 in apple.

RAD5B belongs to the Rad5/16-like group of the SNF2 family, which often functions in chromatin remodelling. However, whether RAD5B is involved in chromatin remodelling, histone modification, and drought stress tolerance is largely unclear. We identified a drought-inducible chromatin remodeler, MdRAD5B, which positively regulates apple drought tolerance. Transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) analysis showed that MdRAD5B affects the expression of 466 drought-responsive genes through its chromatin remodelling function in response to drought stress. In addition, MdRAD5B interacts with and degrades MdLHP1, a crucial regulator of histone H3 trimethylation at K27 (H3K27me3), through the ubiquitin-independent 20S proteasome. Chromatin immunoprecipitation-sequencing (ChIP-seq) analysis revealed that MdRAD5B modulates the H3K27me3 deposition of 615 genes in response to drought stress. Genetic interaction analysis showed that MdRAD5B mediates the H3K27me3 deposition of drought-responsive genes through MdLHP1, which causes their expression changes under drought stress. Our results unravelled a dual function of MdRAD5B in gene expression modulation in apple in response to drought, that is, via the regulation of chromatin remodelling and H3K27me3.

Global hypermethylation of the N6-methyladenosine RNA modification associated with apple heterografting.

Grafting can facilitate better scion performance and is widely used in plants. Numerous studies have studied the involvement of mRNAs, small RNAs, and epigenetic regulations in the grafting process. However, it remains unclear whether the mRNA N6-methyladenosine (m6A) modification participates in the apple (Malus x domestica Borkh.) grafting process. Here, we decoded the landscape of m6A modification profiles in 'Golden delicious' (a cultivar, Gd) and Malus prunifolia 'Fupingqiuzi' (a unique rootstock with resistance to environmental stresses, Mp), as well as their heterografted and self-grafted plants. Interestingly, global hypermethylation of m6A occurred in both heterografted scion and rootstock compared with their self-grafting controls. Gene Ontology (GO) term enrichment analysis showed that grafting-induced differentially m6A-modified genes were mainly involved in RNA processing, epigenetic regulation, stress response, and development. Differentially m6A-modified genes harboring expression alterations were mainly involved in various stress responses and fatty acid metabolism. Furthermore, grafting-induced mobile mRNAs with m6A and gene expression alterations mainly participated in ABA synthesis and transport (e.g. carotenoid cleavage dioxygenase 1 [CCD1] and ATP-binding cassette G22 [ABCG22]) and abiotic and biotic stress responses, which might contribute to the better performance of heterografted plants. Additionally, the DNA methylome analysis also demonstrated the DNA methylation alterations during grafting. Downregulated expression of m6A methyltransferase gene MdMTA (ortholog of METTL3) in apples induced the global m6A hypomethylation and distinctly activated the expression level of DNA demethylase gene MdROS1 (REPRESSOR OF SILENCING 1) showing the possible association between m6A and 5mC methylation in apples. Our results reveal the m6A modification profiles in the apple grafting process and enhance our understanding of the m6A regulatory mechanism in plant biological processes.

HISTONE DEACETYLASE 6 interaction with ABSCISIC ACID-INSENSITIVE 5 decreases apple drought tolerance.

Understanding the molecular regulation of plant response to drought is the basis of drought-resistance improvement through molecular strategies. Here, we characterized apple (Malus × domestica) histone deacetylase 6 (MdHDA6), which negatively regulates apple drought tolerance by catalyzing deacetylation on histones associated with drought-responsive genes. Transgenic apple plants over-expressing MdHDA6 were less drought-tolerant, while those with down-regulated MdHDA6 expression were more drought-resistant than nontransgenic apple plants. Transcriptomic and histone 3 acetylation (H3ac) Chromatin immunoprecipitation-seq analyses indicated that MdHDA6 could facilitate histone deacetylation on the drought-responsive genes, repressing gene expression. Moreover, MdHDA6 interacted with the abscisic acid (ABA) signaling transcriptional factor, ABSCISIC ACID-INSENSITIVE 5 (MdABI5), forming the MdHDA6-MdABI5 complex. Interestingly, MdHDA6 facilitated histone deacetylation on the drought-responsive genes regulated by MdABI5, resulting in gene repression. Furthermore, a dual-Luc experiment showed that MdHDA6 could repress the regulation of a drought-responsive gene, RESPONSIVE TO DESICCATION 29A (MdRD29A) activated by MdABI5. On the one hand, MdHDA6 can facilitate histone deacetylation and gene repression on the positive drought-responsive genes to negatively regulate drought tolerance in apple. On the other hand, MdHDA6 directly interacts with MdABI5 and represses the expression of genes downstream of MdABI5 via histone deacetylation around these genes to reduce drought tolerance. Our study uncovers a different drought response regulatory mechanism in apple based on the MdHDA6-MdABI5 complex function and provides the molecular basis for drought-resistance improvement in apple.

Histone deacetylase MdHDA6 is an antagonist in regulation of transcription factor MdTCP15 to promote cold tolerance in apple.

Understanding the molecular regulation of plant cold response is the basis for cold resistance germplasm improvement. Here, we revealed that the apple histone deacetylase MdHDA6 can perform histone deacetylation on cold-negative regulator genes and repress their expression, leading to the positive regulation of cold tolerance in apples. Moreover, MdHDA6 directly interacts with the transcription factor MdTCP15. Phenotypic analysis of MdTCP15 transgenic apple lines and wild types reveals that MdTCP15 negatively regulates cold tolerance in apples. Furthermore, we found that MdHDA6 can facilitate histone deacetylation of MdTCP15 and repress the expression of MdTCP15, which positively contributes to cold tolerance in apples. Additionally, the transcription factor MdTCP15 can directly bind to the promoter of the cold-negative regulator gene MdABI1 and activate its expression, and it can also directly bind to the promoter of the cold-positive regulator gene MdCOR47 and repress its expression. However, the co-expression of MdHDA6 and MdTCP15 can inhibit MdTCP15-induced activation of MdABI1 and repression of MdCOR47, suggesting that MdHDA6 suppresses the transcriptional regulation of MdTCP15 on its downstream genes. Our results demonstrate that histone deacetylase MdHDA6 plays an antagonistic role in the regulation of MdTCP15-induced transcriptional activation or repression to positively regulate cold tolerance in apples, revealing a new regulatory mechanism of plant cold response.

Advances in Plant Epigenome Editing Research and Its Application in Plants.

Plant epistatic regulation is the DNA methylation, non-coding RNA regulation, and histone modification of gene sequences without altering the genome sequence, thus regulating gene expression patterns and the growth process of plants to produce heritable changes. Epistatic regulation in plants can regulate plant responses to different environmental stresses, regulate fruit growth and development, etc. Genome editing can effectively improve plant genetic efficiency by targeting the design and efficient editing of genome-specific loci with specific nucleases, such as zinc finger nucleases (ZFNs), transcription activator-like effector nucleases (TALEN), and clustered regularly interspaced short palindromic repeats/CRISPR-associated 9 (CRISPR/Cas9). As research progresses, the CRISPR/Cas9 system has been widely used in crop breeding, gene expression, and epistatic modification due to its high editing efficiency and rapid translation of results. In this review, we summarize the recent progress of CRISPR/Cas9 in epigenome editing and look forward to the future development direction of this system in plant epigenetic modification to provide a reference for the application of CRISPR/Cas9 in genome editing.

Integrating ATAC-seq and RNA-seq Reveals the Dynamics of Chromatin Accessibility and Gene Expression in Apple Response to Drought.

Drought resistance in plants is influenced by multiple signaling pathways that involve various transcription factors, many target genes, and multiple types of epigenetic modifications. Studies on epigenetic modifications of drought focus on DNA methylation and histone modifications, with fewer on chromatin remodeling. Changes in chromatin accessibility can play an important role in abiotic stress in plants by affecting RNA polymerase binding and various regulatory factors. However, the changes in chromatin accessibility during drought in apples are not well understood. In this study, the landscape of chromatin accessibility associated with the gene expression of apple (GL3) under drought conditions was analyzed by Assay for Transposase Accessible Chromatin with high-throughput sequencing (ATAC-seq) and RNA-seq. Differential analysis between drought treatment and control identified 23,466 peaks of upregulated chromatin accessibility and 2447 peaks of downregulated accessibility. The drought-induced chromatin accessibility changed genes were mainly enriched in metabolism, stimulus, and binding pathways. By combining results from differential analysis of RNA-seq and ATAC-seq, we identified 240 genes with higher chromatin accessibility and increased gene expression under drought conditions that may play important functions in the drought response process. Among them, a total of nine transcription factor genes were identified, including ATHB7, HAT5, and WRKY26. These transcription factor genes are differentially expressed with different chromatin accessibility motif binding loci that may participate in apple response to drought by regulating downstream genes. Our study provides a reference for chromatin accessibility under drought stress in apples and the results will facilitate subsequent studies on chromatin remodelers and transcription factors.

Variation burst during dedifferentiation and increased CHH-type DNA methylation after 30 years of in vitro culture of sweet orange.

Somaclonal variation arising from tissue culture may provide a valuable resource for the selection of new germplasm, but may not preserve true-to-type characteristics, which is a major concern for germplasm conservation or genome editing. The genomic changes associated with dedifferentiation and somaclonal variation during long-term in vitro culture are largely unknown. Sweet orange was one of the earliest plant species to be cultured in vitro and induced via somatic embryogenesis. We compared four sweet orange callus lines after 30 years of constant tissue culture with newly induced calli by comprehensively determining the single-nucleotide polymorphisms, copy number variations, transposable element insertions, methylomic and transcriptomic changes. We identified a burst of variation during early dedifferentiation, including a retrotransposon outbreak, followed by a variation purge during long-term in vitro culture. Notably, CHH methylation showed a dynamic pattern, initially disappearing during dedifferentiation and then more than recovering after 30 years of in vitro culture. We also analyzed the effects of somaclonal variation on transcriptional reprogramming, and indicated subgenome dominance was evident in the tetraploid callus. We identified a retrotransposon insertion and DNA modification alternations in the potential regeneration-related gene CLAVATA3/EMBRYO SURROUNDING REGION-RELATED 16. This study provides the foundation to harness in vitro variation and offers a deeper understanding of the variation introduced by tissue culture during germplasm conservation, somatic embryogenesis, gene editing, and breeding programs.

Genome-wide analysis of SET-domain group histone methyltransferases in apple reveals their role in development and stress responses.

BACKGROUND: Histone lysine methylation plays an important role in plant development and stress responses by activating or repressing gene expression. Histone lysine methylation is catalyzed by a class of SET-domain group proteins (SDGs). Although an increasing number of studies have shown that SDGs play important regulatory roles in development and stress responses, the functions of SDGs in apple remain unclear. RESULTS: A total of 67 SDG members were identified in the Malus×domestica genome. Syntenic analysis revealed that most of the MdSDG duplicated gene pairs were associated with a recent genome-wide duplication event of the apple genome. These 67 MdSDG members were grouped into six classes based on sequence similarity and the findings of previous studies. The domain organization of each MdSDG class was characterized by specific patterns, which was consistent with the classification results. The tissue-specific expression patterns of MdSDGs among the 72 apple tissues in the different apple developmental stages were characterized to provide insight into their potential functions in development. The expression profiles of MdSDGs were also investigated in fruit development, the breaking of bud dormancy, and responses to abiotic and biotic stress; the results indicated that MdSDGs might play a regulatory role in development and stress responses. The subcellular localization and putative interaction network of MdSDG proteins were also analyzed. CONCLUSIONS: This work presents a fundamental comprehensive analysis of SDG histone methyltransferases in apple and provides a basis for future studies of MdSDGs involved in apple development and stress responses.

Integrative Analyses of Widely Targeted Metabolic Profiling and Transcriptome Data Reveals Molecular Insight into Metabolomic Variations during Apple (Malus domestica) Fruit Development and Ripening.

The apple is a favorite fruit for human diet and is one of the most important commercial fruit crops around the world. Investigating metabolic variations during fruit development can provide a better understanding on the formation of fruit quality. The present study applied a widely targeted LC-MS-based metabolomics approach with large-scale detection, identification and quantification to investigate the widespread metabolic changes during "Pinova" apple development and ripening. A total of 462 primary and secondary metabolites were simultaneously detected, and their changes along with the four fruit-development stages were further investigated. The results indicated that most of the sugars presented increasing accumulation levels while organic acid, including Tricarboxylic acid cycle (TCA) intermediates, showed a distinct decreasing trend across the four fruit-development stages. A total of 207 secondary metabolites consisted of 104 flavonoids and 103 other secondary metabolites. Many flavonoids maintained relatively high levels in the early fruit stage and then rapidly decreased their levels at the following developmental stages. Further correlation analyses of each metabolite-metabolite pair highlighted the cross talk between the primary and secondary metabolisms across fruit development and ripening, indicating the significant negative correlations between sugars and secondary metabolites. Moreover, transcriptome analysis provided the molecular basis for metabolic variations during fruit development. The results showed that most differentially expressed genes (DEGs) involved in the TCA cycle were upregulated from the early fruit stage to the preripening stage. The extensive downregulation of controlling genes involved in the flavonoid pathway is probably responsible for the rapid decrease of flavonoid content at the early fruit stage. These data provide a global view of the apple metabolome and a comprehensive analysis on metabolomic variations during fruit development, providing a broader and better understanding on the molecular and metabolic basis of important fruit quality traits in commercial apples.

Single-base methylome analysis reveals dynamic epigenomic differences associated with water deficit in apple.

Cytosine methylation is an essential feature of epigenetic regulation and is involved in various biological processes. Although cytosine methylation has been analysed at the genomic scale for several plant species, there is a general lack of understanding of the dynamics of global and genic DNA methylation in plants growing in environments challenged with biotic and abiotic stresses. In this study, we mapped cytosine methylation at single-base resolution in the genome of commercial apple (Malus x domestica), and analysed changes in methylation patterns associated with water deficit in representative drought-sensitive and drought-tolerant cultivars. We found that the apple genome exhibits ~54%, ~38% and ~8.5% methylation at CG, CHG and CHH sequence contexts, respectively. We additionally documented changes in gene expression associated with water deficit in an attempt to link methylation and gene expression changes. Global methylation and transcription analysis revealed that promoter-unmethylated genes showed higher expression levels than promoter-methylated genes. Gene body methylation appears to be positively correlated with gene expression. Water deficit stress was associated with changes in methylation at a multitude of genes, including those encoding transcription factors (TFs) and transposable elements (TEs). These results present a methylome map of the apple genome and reveal widespread DNA methylation alterations in response to water deficit stress. These data will be helpful for understanding potential linkages between DNA methylation and gene expression in plants growing in natural environments and challenged with abiotic and biotic stresses.

Dynamic changes in methylome and transcriptome patterns in response to methyltransferase inhibitor 5-azacytidine treatment in citrus.

DNA methylation is known to play an important role in various developmental processes in plants. However, there is a general lack of understanding about the possible functions of DNA methylation in fruit trees. Using callus as a model, methylome, transcriptome and metabolite changes were assessed after treatment with the DNA methyltransferase inhibitor 5-azacytidine (5azaC). Genome-wide methylome analysis revealed the demethylation of a diverse of genes, including many genes encoding transcription factors (TFs), genes involved in biological processes, and the up-regulation of a wide range of transposable elements (TEs). Combined with the RNA-seq data, we observed no obvious genome-wide correlation between the changes in methylation status and expression levels. Furthermore, 5azaC treatment induced carotenoid degradation along with strong activation of carotenoid cleavage dioxygenases 1 (CpCCD1). Functional complementation analysis in bacterial system showed that CpCCD1 exhibited strong catalytic activities toward zeaxanthin, ß-carotene and lycopene. In summary, 5azaC treatments induced carotenoid degradation by CpCCD1 activation and led to a genome-wide demethylation effect.

Genomic analyses of primitive, wild and cultivated citrus provide insights into asexual reproduction.

The emergence of apomixis-the transition from sexual to asexual reproduction-is a prominent feature of modern citrus. Here we de novo sequenced and comprehensively studied the genomes of four representative citrus species. Additionally, we sequenced 100 accessions of primitive, wild and cultivated citrus. Comparative population analysis suggested that genomic regions harboring energy- and reproduction-associated genes are probably under selection in cultivated citrus. We also narrowed the genetic locus responsible for citrus polyembryony, a form of apomixis, to an 80-kb region containing 11 candidate genes. One of these, CitRWP, is expressed at higher levels in ovules of polyembryonic cultivars. We found a miniature inverted-repeat transposable element insertion in the promoter region of CitRWP that cosegregated with polyembryony. This study provides new insights into citrus apomixis and constitutes a promising resource for the mining of agriculturally important genes.

Improved hybrid de novo genome assembly of domesticated apple (Malus x domestica).

BACKGROUND: Domesticated apple (Malus × domestica Borkh) is a popular temperate fruit with high nutrient levels and diverse flavors. In 2012, global apple production accounted for at least one tenth of all harvested fruits. A high-quality apple genome assembly is crucial for the selection and breeding of new cultivars. Currently, a single reference genome is available for apple, assembled from 16.9 × genome coverage short reads via Sanger and 454 sequencing technologies. Although a useful resource, this assembly covers only ~89 % of the non-repetitive portion of the genome, and has a relatively short (16.7 kb) contig N50 length. These downsides make it difficult to apply this reference in transcriptive or whole-genome re-sequencing analyses. FINDINGS: Here we present an improved hybrid de novo genomic assembly of apple (Golden Delicious), which was obtained from 76 Gb (~102 × genome coverage) Illumina HiSeq data and 21.7 Gb (~29 × genome coverage) PacBio data. The final draft genome is approximately 632.4 Mb, representing ~ 90 % of the estimated genome. The contig N50 size is 111,619 bp, representing a 7 fold improvement. Further annotation analyses predicted 53,922 protein-coding genes and 2,765 non-coding RNA genes. CONCLUSIONS: The new apple genome assembly will serve as a valuable resource for investigating complex apple traits at the genomic level. It is not only suitable for genome editing and gene cloning, but also for RNA-seq and whole-genome re-sequencing studies.

Genome-wide identification of sweet orange (Citrus sinensis) histone modification gene families and their expression analysis during the fruit development and fruit-blue mold infection process.

In eukaryotes, histone acetylation and methylation have been known to be involved in regulating diverse developmental processes and plant defense. These histone modification events are controlled by a series of histone modification gene families. To date, there is no study regarding genome-wide characterization of histone modification related genes in citrus species. Based on the two recent sequenced sweet orange genome databases, a total of 136 CsHMs (Citrus sinensis histone modification genes), including 47 CsHMTs (histone methyltransferase genes), 23 CsHDMs (histone demethylase genes), 50 CsHATs (histone acetyltransferase genes), and 16 CsHDACs (histone deacetylase genes) were identified. These genes were categorized to 11 gene families. A comprehensive analysis of these 11 gene families was performed with chromosome locations, phylogenetic comparison, gene structures, and conserved domain compositions of proteins. In order to gain an insight into the potential roles of these genes in citrus fruit development, 42 CsHMs with high mRNA abundance in fruit tissues were selected to further analyze their expression profiles at six stages of fruit development. Interestingly, a numbers of genes were expressed highly in flesh of ripening fruit and some of them showed the increasing expression levels along with the fruit development. Furthermore, we analyzed the expression patterns of all 136 CsHMs response to the infection of blue mold (Penicillium digitatum), which is the most devastating pathogen in citrus post-harvest process. The results indicated that 20 of them showed the strong alterations of their expression levels during the fruit-pathogen infection. In conclusion, this study presents a comprehensive analysis of the histone modification gene families in sweet orange and further elucidates their behaviors during the fruit development and the blue mold infection responses.

Comprehending crystalline ß-carotene accumulation by comparing engineered cell models and the natural carotenoid-rich system of citrus.

Genetic manipulation of carotenoid biosynthesis has become a recent focus for the alleviation of vitamin A deficiency. However, the genetically modified phenotypes often challenge the expectation, suggesting the incomplete comprehension of carotenogenesis. Here, embryogenic calli were engineered from four citrus genotypes as engineered cell models (ECMs) by over-expressing a bacterial phytoene synthase gene (CrtB). Ripe flavedos (the coloured outer layer of citrus fruits), which exhibit diverse natural carotenoid patterns, were offered as a comparative system to the ECMs. In the ECMs, carotenoid patterns showed diversity depending on the genotypes and produced additional carotenoids, such as lycopene, that were absent from the wild-type lines. Especially in the ECMs from dark-grown culture, there emerged a favoured ß,ß-pathway characterized by a striking accumulation of ß-carotene, which was dramatically different from those in the wild-type calli and ripe flavedos. Unlike flavedos that contained a typical chromoplast development, the ECMs sequestered most carotenoids in the amyloplasts in crystal form, which led the amyloplast morphology to show a chromoplast-like profile. Transcriptional analysis revealed a markedly flavedo-specific expression of the ß-carotene hydroxylase gene (HYD), which was suppressed in the calli. Co-expression of CrtB and HYD in the ECMs confirmed that HYD predominantly mediated the preferred carotenoid patterns between the ECMs and flavedos, and also revealed that the carotenoid crystals in the ECMs were mainly composed of ß-carotene. In addition, a model is proposed to interpret the common appearance of a favoured ß,ß-pathway and the likelihood of carotenoid degradation potentially mediated by photo-oxidation and vacuolar phagocytosis in the ECMs is discussed.

Phylogenetic and evolutionary analysis of NBS-encoding genes in Rutaceae fruit crops.

The nucleotide-binding site leucine-rich repeat (NBS-LRR) genes are the largest class of disease resistance genes in plants. However, our understanding of the evolution of NBS-LRR genes in Rutaceae fruit crops is rather limited. We report an evolutionary study of 103 NBS-encoding genes isolated from Poncirus trifoliata (trifoliate orange), Citrus reticulata (tangerine) and their F(1) progeny. In all, 58 of the sequences contained a continuous open reading frame. Phylogenetic analysis classified the 58 NBS genes into nine clades, eight of which were genus specific. This was taken to imply that most of the ancestors of these NBS genes evolved after the genus split. The motif pattern of the 58 NBS-encoding genes was consistent with their phylogenetic profile. An extended phylogenetic analysis, incorporating citrus NBS genes from the public database, classified 95 citrus NBS genes into six clades, half of which were genus specific. RFLP analysis showed that citrus NBS-encoding genes have been evolving rapidly, and that they are unstable when passed through an intergeneric cross. Of 32 NBS-encoding genes tracked by gene-specific PCR, 24 showed segregation distortion among a set of 94 F(1) individuals. This study provides new insight into the evolution of Rutaceae NBS genes and their behaviour following an intergeneric cross.