Chromosome-scale reference genome provides insights into the genetic origin and grafting-mediated stress tolerance of Malus prunifolia

Genome-wide identification of starch phosphorylase gene family in Rosa chinensis and expression in response to abiotic stress

Chinese rose (Rosa chinensis) is an important ornamental plant, with economic, cultural, and symbolic significance. During the application of outdoor greening, adverse environments such as high temperature and drought are often encountered, which affect its application scope and ornamental quality. The starch phosphorylase (Pho) gene family participate in the synthesis and decomposition of starch, not only related to plant energy metabolism, but also plays an important role in plant stress resistance. The role of Pho in combating salinity and high temperature stress in R. chinensis remains unknown. In this work, 4 Phos from R. chinensis were detected with Pfam number of Pho (PF00343.23) and predicted by homolog-based prediction (HBP). The Phos are characterized by sequence lengths of 821 to 997 bp, and the proteins are predicted to subcellularly located in the plastid and cytoplasm. The regulatory regions of the Phos contain abundant stress and phytohormone-responsive cis-acting elements. Based on transcriptome analysis, the Phos were found to respond to abiotic stress factors such as drought, salinity, high temperature, and plant phytohormone of jasmonic acid and salicylic acid. The response of Phos to abiotic stress factors such as salinity and high temperature was confirmed by qRT-PCR analysis. To evaluate the genetic characteristics of Phos, a total of 69 Phos from 17 species were analyzed and then classified into 3 groups in phylogenetic tree. The collinearity analysis of Phos in R. chinensis and other species was conducted for the first time. This work provides a view of evolution for the Pho gene family and indicates that Phos play an important role in abiotic stress response of R. chinensis.

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.

Chromosome-level phased genome assembly of "Antonovka" identified candidate apple scab-resistance genes highly homologous to HcrVf2 and HcrVf1 on linkage group 1

Apple scab, a fungal disease caused by Venturia inaequalis, leads to losses in both yield and fruit quality of apples (Malus domestica Borkh.). Most commercial apple cultivars, including those containing the well-characterized Rvi6-scab-resistance locus on linkage group (LG) 1, are susceptible to scab. HcrVf2 and HcrVf1 are considered the main paralogs of the Rvi6 locus. The major apple scab-resistance loci Vhc1 in "Honeycrisp" and Rvi17 in "Antonovka," were identified in close proximity to HcrVf2. In this study, we used long-read sequencing and in silico gene sequence characterization to identify candidate resistance genes homologous to HcrVf2 and HcrVf1 in Honeycrisp and Antonovka. Previously published chromosome-scale phased assembly of Honeycrisp and a newly assembled phased genome of Antonovka 172670-B were used to identify HcrVf2 and HcrVf1 homologs spanning Vhc1 and Rvi17 loci. In combination with 8 available Malus assemblies, 43 and 46 DNA sequences highly homologous to HcrVf2 and HcrVf1, respectively, were identified on LG 1 and 6, with identity and coverage ranging between 87-95 and 81-95%, respectively. Among these homologs, 2 candidate genes in Antonovka and Honeycrisp haplome A are located in close physical proximity to the scab-resistance marker Ch-Vf1 on LG 1. They showed the highest identity and coverage (95%) of HcrVf2 and only minor changes in the protein motifs. They were identical by state between each other, but not with HcrVf2. This study offers novel genomic resources and insights into the Vhc1 and Rvi17 loci on LG 1 and identifies candidate genes for further resistance characterization.

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.

A haplotype resolved chromosome-scale assembly of North American wild apple Malus fusca and comparative genomics of the fire blight Mfu10 locus

The Pacific crabapple (Malus fusca) is a wild relative of the commercial apple (Malus × domestica). With a range extending from Alaska to Northern California, M. fusca is extremely hardy and disease resistant. The species represents an untapped genetic resource for the development of new apple cultivars with enhanced stress resistance. However, gene discovery and utilization of M. fusca have been hampered by the lack of genomic resources. Here, we present a high‐quality, haplotype‐resolved, chromosome‐scale genome assembly and annotation for M. fusca. The genome was assembled using high‐fidelity long‐reads and scaffolded using genetic maps and high‐throughput chromatin conformation capture sequencing, resulting in one of the most contiguous apple genomes to date. We annotated the genome using public transcriptomic data from the same species taken from diverse plant structures and developmental stages. Using this assembly, we explored haplotypic structural variation within the genome of M. fusca, identifying thousands of large variants. We further showed high sequence co‐linearity with other domesticated and wild Malus species. Finally, we resolve a known quantitative trait locus associated with resistance to fire blight (Erwinia amylovora). Insights gained from the assembly of a reference‐quality genome of this hardy wild apple relative will be invaluable as a tool to facilitate DNA‐informed introgression breeding.

Insights into the molecular mechanisms underlying responses of apple trees to abiotic stresses

Apple (Malus[Formula: see text]domestica) is a popular temperate fruit crop worldwide. However, its growth, productivity, and quality are often adversely affected by abiotic stresses such as drought, extreme temperature, and high salinity. Due to the long juvenile phase and highly heterozygous genome, the conventional breeding approaches for stress-tolerant cultivars are time-consuming and resource-intensive. These issues may be resolved by feasible molecular breeding techniques for apples, such as gene editing and marker-assisted selection. Therefore, it is necessary to acquire a more comprehensive comprehension of the molecular mechanisms underpinning apples' response to abiotic stress. In this review, we summarize the latest research progress in the molecular response of apples to abiotic stressors, including the gene expression regulation, protein modifications, and epigenetic modifications. We also provide updates on new approaches for improving apple abiotic stress tolerance, while discussing current challenges and future perspectives for apple molecular breeding.

Molecular Structure and Variation Characteristics of the Plastomes from Six Malus baccata (L.) Borkh. Individuals and Comparative Genomic Analysis with Other Malus Species

Malus baccata (L.) Borkh. is an important wild species of Malus. Its rich variation types and population history are not well understood. Chloroplast genome mining plays an active role in germplasm identification and genetic evolution. In this study, by assembly and annotation, six complete cp genome sequences, ranging in size from 160,083 to 160,295 bp, were obtained. The GC content of stable IR regions (42.7%) was significantly higher than that of full length (36.5%) and SC regions (LSC-34.2%, SSC-30.4%). Compared with other Malus species, it was found that there were more sites of polymorphisms and hotspots of variation in LSC and SSC regions, with high variation sites including trnR/UCU-atpA, trnT/UGU-trnL/UAA, ndhF-rpl32 and ccsA-ndhD. The intraspecific and interspecific collinearity was good, and no structural rearrangement was observed. A large number of repeating elements and different boundary expansions may be involved in shaping the cp genome size. Up to 77 or 78 coding genes were annotated in the cp genomes of M. baccata, and high frequency codons such as UUA (Leu), GCU (Ala) and AGA (Arg) were identified by relative synonymous codon usage analysis. Phylogeographic analysis showed that 12 individuals of M. baccata clustered into three different groups with complex structure, whereas variant xiaojinensis (M.H. Cheng & N.G. Jiang) was not closely related to M. baccata evolutionarily. The phylogenetic analysis suggested that two main clades of different M. baccata in the genus Malus were formed and that I and II diverged about 9.7 MYA. In conclusion, through cp genome assembly and comparison, the interspecific relationships and molecular variations of M. baccata were further elucidated, and the results of this study provide valuable information for the phylogenetic evolution and germplasm conservation of M. baccata and Malus.

The role of DNA methylation in the maintenance of phenotypic variation induced by grafting chimerism in Brassica

Grafting facilitates the interaction between heterologous cells with different genomes, resulting in abundant phenotypic variation, which provides opportunities for crop improvement. However, how grafting-induced variation occurs and is transmitted to progeny remains elusive. A graft chimera, especially a periclinal chimera, which has genetically distinct cell layers throughout the plant, is an excellent model to probe the molecular mechanisms of grafting-induced variation maintenance. Here we regenerated a plant from the T-cell layer of a periclinal chimera, TCC (where the apical meristem was artificially divided into three cell layers - from outside to inside, L1, L2, and L3; T = Tuber mustard, C = red Cabbage), named rTTT0 (r = regenerated). Compared with the control (rsTTT, s = self-grafted), rTTT0 had multiple phenotypic variations, especially leaf shape variation, which could be maintained in sexual progeny. Transcriptomes were analyzed and 58 phenotypic variation-associated genes were identified. Whole-genome bisulfite sequencing analyses revealed that the methylome of rTTT0 was changed, and the CG methylation level was significantly increased by 8.74%. In rTTT0, the coding gene bodies are hypermethylated in the CG context, while their promoter regions are hypomethylated in the non-CG context. DNA methylation changes in the leaf shape variation-associated coding genes, ARF10, IAA20, ROF1, and TPR2, were maintained for five generations of rTTT0. Interestingly, grafting chimerism also affected transcription of the microRNA gene (MIR), among which the DNA methylation levels of the promoters of three MIRs associated with leaf shape variation were changed in rTTT0, and the DNA methylation modification of MIR319 was maintained to the fifth generation of selfed progeny of rTTT0 (rTTT5). These findings demonstrate that DNA methylation of coding and non-coding genes plays an important role in heterologous cell interaction-induced variation formation and its transgenerational inheritance.

De novo transcriptome assembly and functional analysis reveal a dihydrochalcone 3-hydroxylase(DHC3H) of wild Malus species that produces sieboldin in vivo

Sieboldin is a specialised secondary metabolite of the group of dihydrochalcones (DHC), found in high concentrations only in some wild Malus species, closely related to the domesticated apple (Malus × domestica L.). To date, the first committed step towards the biosynthesis of sieboldin remains unknown. In this study, we combined transcriptomic analysis and a de novo transcriptome assembly to identify two putative 3-hydroxylases in two wild Malus species (Malus toringo (K. Koch) Carriere syn. sieboldii Rehder, Malus micromalus Makino) whose DHC profile is dominated by sieboldin. We assessed the in vivo activity of putative candidates to produce 3-hydroxyphloretin and sieboldin by de novo production in Saccharomyces cerevisiae. We found that CYP98A proteins of wild Malus accessions (CYP98A195, M. toringo and CYP98A196, M. micromalus) were able to produce 3-hydroxyphloretin, ultimately leading to sieboldin accumulation by co-expression with PGT2. CYP98A197-198 genes of M. × domestica, however, were unable to hydroxylate phloretin in vivo. CYP98A195-196 proteins exerting 3-hydroxylase activity co-localised with an endoplasmic reticulum marker. CYP98A protein model from wild accessions showed mutations in key residues close to the ligand pocket predicted using phloretin for protein docking modelling. These mutations are located within known substrate recognition sites of cytochrome P450s, which could explain the acceptance of phloretin in CYP98A protein of wild accessions. Screening a Malus germplasm collection by HRM marker analysis for CYP98A genes identified three clusters that correspond to the alleles of domesticated and wild species. Moreover, CYP98A isoforms identified in M. toringo and M. micromalus correlate with the accumulation of sieboldin in other wild and hybrid Malus genotypes. Taken together, we provide the first evidence of an enzyme producing sieboldin in vivo that could be involved in the key hydroxylation step towards the synthesis of sieboldin in Malus species.

Genome-wide identification and stress response analysis of cyclophilin gene family in apple (Malus × domestica)

BACKGROUND

Cyclophilin (CYP) belongs to the immunophilin family and has peptidyl-prolyl cis-trans isomerase (PPIase) activity, which catalyzes the cis-trans isomerization process of proline residues. CYPs widely exist in eukaryotes and prokaryotes, and contain a conserved cyclophilin-like domain (CLD). Plant cyclophilins are widely involved in a range of biological processes including stress response, metabolic regulation, and growth and development.

RESULT

In this study, 30 cyclophilin genes on 15 chromosomes were identified from the 'Golden Delicious' apple (M. domestica) genome. Phylogenetic analysis showed that the cyclophilin family genes can be divided into three clades in Malus. Collinear analysis showed that ten gene pairs were the result of segmental duplication. Analysis of gene and protein structure further supported the phylogenetic tree and collinearity analysis. The expression of MdCYPs in different organs was higher in leaves, flowers, and fruits. Ten and eight CYPs responded to drought and salt stress, respectively. MdCYP16, a nuclear-localized MD CYP, was screened from the intersection of the two expression profiling datasets and was highly sensitive to drought and salt stress. GUS staining of transgenic Arabidopsis indicated that MdCYP16 may be involved in the regulation of abiotic stress.

CONCLUSION

This study systematically analyzed members of the apple cyclophilin family and confirmed the involvement of MdCYP16 as a nuclear-localized MD cyclophilin that acts in response to salt and drought stress in apple. Our work identifies members of the apple cyclophilin gene family, and provides an important theoretical basis for in-depth study of cyclophilin function. Additionally, the analysis provides candidate genes that may be involved in stress response in apple.