High-quality de novo assembly of the apple genome and methylome dynamics of early fruit development

Upregulation of PECTATE LYASE5 by a NAC transcription factor promotes fruit softening in apple

Flesh firmness is a critical breeding trait that determines consumer selection, shelf life, and transportation. The genetic basis controlling firmness in apple (Malus × domestica Borkh.) remains to be fully elucidated. We aimed to decipher genetic variance for firmness at harvest and develop potential molecular markers for marker-assisted breeding. Maturity firmness for 439 F1 hybrids from a cross of "Cripps Pink" and "Fuji" was determined in 2016 and 2017. The phenotype segregated extensively, with a Gaussian distribution. In a combined bulked segregant analysis (BSA) and RNA-sequencing analysis, 84 differentially expressed genes were screened from the 10 quantitative trait loci regions. Interestingly, next-generation re-sequencing analysis revealed a Harbinger-like transposon element insertion upstream of the candidate gene PECTATE LYASE5 (MdPL5); the genotype was associated with flesh firmness at harvest. The presence of this transposon repressed MdPL5 expression and was closely linked to the extra-hard phenotype. MdPL5 was demonstrated to promote softening in apples and tomatoes. Subsequently, using the MdPL5 promoter as bait, MdNAC1-L was identified as a transcription activator that positively regulates ripening and softening in the developing fruit. We also demonstrated that MdNAC1-L could induce the up-regulation of MdPL5, MdPG1, and the ethylene-related genes MdACS1 and MdACO1. Our findings provide insight into TE-related genetic variation and the PL-mediated regulatory network for the firmness of apple fruit.

Genome-wide analyses of Ariadne family genes reveal their involvement in abiotic stress responses in apple

E3 ligases are essential for ubiquitination and play a role in regulating various aspects of eukaryotic life. Ariadne (ARI) proteins, a subfamily of RBR (RING-between-RING) proteins, have been recognized as a new class of RING-finger E3 ligases. Recent research has shed light on their potential involvement in plants' responses to abiotic stress. However, comprehensive studies on ARI genes in apple (Malus domestica) are still lacking. This study identified ten MdARI genes in the apple genome, and examined intraspecific and interspecific collinearity to explore the evolutionary relationships of ARI family members. Phylogenetic analyses classified MdARIs into two subfamilies (A and B), and by integrating gene structure, conserved motifs, and sequence comparison results, subfamily B was further divided into two subgroups (I and II). Tissue expression analyses revealed varied expression patterns of MdARI genes in different tissues, and subcellular localization showed that MdARI1-1, MdARI1-2, and MdARI9-1 were located in the nucleus, while the other seven MdARIs were distributed throughout the cell. Analyses of promoter cis-elements and expression patterns under cold, salt, and drought treatments indicated the involvement of MdARIs in abiotic stress responses. Several proteins crucial to the plant stress response were predicted to be potential MdARIs-interacting proteins based on the protein interaction network. Additionally, the interaction between UBC11 (E2) and MdARI7-2 was confirmed by a yeast two-hybrid (Y2H) experiment, suggesting that MdARI7-2 may function as an E3. These findings will greatly benefit future research on the role and mechanisms of ARI proteins in apple stress response.

Haplotype-resolved genome assembly and resequencing provide insights into the origin and breeding of modern rose

Modern rose (Rosa hybrida) is a recently formed interspecific hybrid and has become one of the most important and widely cultivated ornamentals. Here we report the haplotype-resolved chromosome-scale genome assembly of the tetraploid R. hybrida 'Samantha' ('JACmantha') and a genome variation map of 233 Rosa accessions involving various wild species, and old and modern cultivars. Homologous chromosomes of 'Samantha' exhibit frequent homoeologous exchanges. Population genomic and genomic composition analyses reveal the contributions of wild Rosa species to modern roses and highlight that R. odorata and its derived cultivars are important contributors to modern roses, much like R. chinensis 'Old Blush'. Furthermore, selective sweeps during modern rose breeding associated with major agronomic traits, including continuous and recurrent flowering, double flower, flower senescence and disease resistance, are identified. This study provides insights into the genetic basis of modern rose origin and breeding history, and offers unprecedented genomic resources for rose improvement.

Telomere-to-telomere Genome Assembly of two representative Asian and European pear cultivars

As the third most important temperate fruit, Pear (Pyrus spp.) exhibits a remarkable genetic diversity and is classified into two mainly categories known as Asian pear and European pear. Although several pear genomes are available, most of the released versions are fragmented and not chromosome-level high-quality. In this study, we report two high-quality genomes for Pyrus bretschneideri Rhed. cv. 'Danshansuli' (DS) and Pyrus communis L. cv. 'Conference' (KFL), which represent the predominant Asian and European cultivars, respectively, with nearly telomere-to-telomere (T2T) gap-free level. The finally assembled genome sizes for DS and KFL were 510.98 Mb and 510.71 Mb, respectively, with Contig N50 of 29.47 Mb and 30.47 Mb, where each chromosome was represented by a single contig. The DS and KFL genomes yielded a total of 46,394 and 44,702 protein-coding genes, respectively. Among these genes, the functional annotation accounted for 96.47% and 96.46% in the DS and KFL genomes. The two novels nearly T2T genomic information offers an invaluable resource for comparative genomics, genetic diversity analysis, molecular breeding strategies, and functional exploration.

Integration of mRNA-miRNA Reveals the Possible Role of PyCYCD3 in Increasing Branches Through Bud-Notching in Pear (Pyrus bretschneideri Rehd.)

Bud-notching in pear varieties with weak-branches enhances branch development, hormone distribution, and germination, promoting healthier growth and improving early yield. To examine the regulatory mechanisms of endogenous hormones on lateral bud germination in Pyrus spp. (cv. 'Huangguan') (Pyrus bretschneideri Rehd.), juvenile buds were collected from 2-year-old pear trees. Then, a comprehensive study, including assessments of endogenous hormones, germination and branching rates, RNA-seq analysis, and gene function analysis in these lateral buds was conducted. The results showed that there was no significant difference in germination rate between the control and bud-notching pear trees, but the long branch rate was significantly increased in bud-notching pear trees compared to the control (p < 0.05). After bud-notching, there was a remarkable increase in IAA and BR levels in the pruned section of shoots, specifically by 141% and 93%, respectively. However, the content of ABA in the lateral buds after bud-notching was not significantly different from the control. Based on RNA-seq analysis, a notable proportion of the differentially expressed genes (DEGs) were linked to the plant hormone signal transduction pathway. Notably, the brassinosteroid signaling pathway seemed to have the closest connection with the branching ability of pear with the related genes encoding BRI1 and CYCD3, which showed significant differences between lateral buds. Finally, the heterologous expression of PyCYCD3 has a positive regulatory effect on the increased Arabidopsis growth and branching numbers. Therefore, the PyCYCD3 was identified as an up-regulated gene that is induced via brassinosteroid (BR) and could act as a conduit, transforming bud-notching cues into proliferative signals, thereby governing lateral branching mechanisms in pear trees.

Comprehensive insights on association mapping in perennial fruit crops breeding - Its implications, current status and future perspectives

In order to provide food and nutritional security for the world's rapidly expanding population, fruit crop researchers have identified two critical priorities: increasing production and preserving fruit quality during the pre- and post-harvest periods. The genetic basis of these complex, commercially important fruit traits which are uniquely regulated by polygenes or multi-allelic genes that interact with one another and the environment can be analyzed with the aid of trait mapping tools. The most interesting trait mapping approach that offers the genetic level investigation for marker-trait associations (MTAs) for these complex fruit traits, without the development of mapping population, is association mapping. This approach was used during the genetic improvement program, emphasizing the obstacles (breeding strategies adopted, generation interval, and their genomic status) pertaining to perennial fruit crops. This method of studying population diversity and linkage disequilibrium in perennial fruit crops has been made possible by recent developments in genotyping, phenotyping, and statistical analysis. Thus, the purpose of this review is to provide an overview of different trait mapping techniques, with a focus on association mapping (method, essential components, viability, constraints, and future perspective) and its advantages, disadvantages, and possibilities for breeding perennial fruit crops.

MdNAC5: a key regulator of fructose accumulation in apple fruit

The sweetness of apple fruit is a key factor in the improvement of apple varieties, with fructose being the sweetest of the soluble sugars, playing a crucial role in determining the overall sweetness of the apple. Therefore, uncovering the key genes controlling fructose accumulation and deciphering the regulatory mechanisms of fructose are vitally important for the improvement of apple varieties. In this study, through BSA-seq and transcriptome analysis of the 'Changfu 2' × 'Golden Delicious' F1 hybrid population, MdNAC5 was identified as a key regulatory gene for fructose content. MdNAC5 was shown to significantly influence fructose accumulation in both apples and tomatoes. Furthermore, we conducted a detailed identification of sugar transporters and metabolic enzymes in apples, discovering that MdNAC5 can enhance fructose accumulation in vacuoles and the conversion of sucrose to fructose by binding to and activating the promoters of the vacuolar sugar transporter MdTST2 and the neutral invertase MdNINV6. Additionally, MdNAC5 regulated the MdEIN3.4-MdSWEET15a module, strengthening the unloading of sucrose in the phloem of the fruit. Our results reveal a new mechanism by which MdNAC5 regulates fructose accumulation in apples and provide theoretical foundations for improving apple sweetness through genetic modification.

Advancing apple genetics research: Malus coronaria and Malus ioensis genomes and a gene family-based pangenome of native North American apples

Wild Malus species flourished in North America long before Europeans introduced domesticated apples. Malus coronaria and M. ioensis are native to the mid-western and eastern United States, while M. angustifolia and M. fusca grow in the southeast and west, respectively. They offer disease resistance, climate and soil adaptability, and horticultural traits for apple breeding. However, their utilization remains limited due to insufficient genomic resources and specific genetics. We report high-quality phased chromosome-scale assemblies of M. coronaria and M. ioensis, generated using long-read and conformation capture sequencing. Phylogenetic and synteny analysis indicated high relatedness between these 2 genomes and previously published genome of M. angustifolia, and lower relatedness with M. fusca. Gene family-based pangenome of North American Malus identified 60,211 orthogroups containing 340,087 genes. Genes involved in basic cellular and metabolic processes, growth, and development were core to the existence of these species, whereas genes involved in secondary metabolism, stress response, and interactions with other organisms were accessory and are likely associated with adaptation to specific environments. Structural variation hotspots were mostly overlapping with high gene density. This study offers novel native North American Malus genome resources that can be used to identify genes for apple breeding and understand their evolution and adaptation.

The identification and analysis of meristematic mutations within the apple tree that developed the RubyMac sport mutation

BACKGROUND

Understanding the molecular basis of sport mutations in fruit trees has the potential to accelerate generation of improved cultivars.

RESULTS

For this, we analyzed the genome of the apple tree that developed the RubyMac phenotype through a sport mutation that led to the characteristic fruit coloring of this variety. Overall, we found 46 somatic mutations that distinguished the mutant and wild-type branches of the tree. In addition, we found 54 somatic gene conversions (i.e., loss-of-heterozygosity mutations) that also distinguished the two parts of the tree. Approximately 20% of the mutations were specific to individual cell lineages, suggesting that they originated from the corresponding meristematic layers. Interestingly, the de novo mutations were enriched for GC =  > AT transitions while the gene conversions showed the opposite bias for AT =  > GC transitions, suggesting that GC-biased gene conversions have the potential to counteract the AT-bias of de novo mutations. By comparing the gene expression patterns in fruit skins from mutant and wild-type branches, we found 56 differentially expressed genes including 18 involved in anthocyanin biosynthesis. While none of the differently expressed genes harbored a somatic mutation, we found that some of them in regions of the genome that were recently associated with natural variation in fruit coloration.

CONCLUSION

Our analysis revealed insights in the characteristics of somatic change, which not only included de novo mutations but also gene conversions. Some of these somatic changes displayed strong candidate mutations for the change in fruit coloration in RubyMac.

Chromosome-scale genome assembly of Docynia delavayi provides new insights into the α-farnesene biosynthesis

Docynia delavayi is an economically significant fruit species with a high market potential due to the special aroma of its fruit. Here, a 653.34 Mb high-quality genome of D. delavayi was first reported, of which 93.8 % of the sequences (612.98 Mb) could be anchored to 17 chromosomes, containing 48,325 protein-coding genes. Ks analysis proved that two whole genome duplication (WGD) events occurred in D. delavayi, resulting in the expansion of genes associated with terpene biosynthesis, which promoted its fruit-specific aroma production. Combined multi-omics analysis, α-farnesene was detected as the most abundant aroma substance emitted by D. delavayi fruit during storage, meanwhile one α-farnesene synthase gene (AFS) and 15 transcription factors (TFs) were identified as the candidate genes potentially involved in α-farnesene biosynthesis. Further studies for the regulation network of α-farnesene biosynthesis revealed that DdebHLH, DdeERF1 and DdeMYB could activate the transcription of DdeAFS. To our knowledge, it is the first report that MYB TF plays a regulatory role in α-farnesene biosynthesis, which will greatly facilitate future breeding programs for D. delavayi.

Dormancy regulator Prunus mume DAM6 promotes ethylene-mediated leaf senescence and abscission

Leaf senescence and abscission in autumn are critical phenological events in deciduous woody perennials. After leaf fall, dormant buds remain on deciduous woody perennials, which then enter a winter dormancy phase. Thus, leaf fall is widely believed to be linked to the onset of dormancy. In Rosaceae fruit trees, DORMANCY-ASSOCIATED MADS-box (DAM) transcription factors control bud dormancy. However, apart from their regulatory effects on bud dormancy, the biological functions of DAMs have not been thoroughly characterized. In this study, we revealed a novel DAM function influencing leaf senescence and abscission in autumn. In Prunus mume, PmDAM6 expression was gradually up-regulated in leaves during autumn toward leaf fall. Our comparative transcriptome analysis using two RNA-seq datasets for the leaves of transgenic plants overexpressing PmDAM6 and peach (Prunus persica) DAM6 (PpeDAM6) indicated Prunus DAM6 may up-regulate the expression of genes involved in ethylene biosynthesis and signaling as well as leaf abscission. Significant increases in 1-aminocyclopropane-1-carboxylate accumulation and ethylene emission in DEX-treated 35S:PmDAM6-GR leaves reflect the inductive effect of PmDAM6 on ethylene biosynthesis. Additionally, ethephon treatments promoted autumn leaf senescence and abscission in apple and P. mume, mirroring the changes due to PmDAM6 overexpression. Collectively, these findings suggest that PmDAM6 may induce ethylene emission from leaves, thereby promoting leaf senescence and abscission. This study clarified the effects of Prunus DAM6 on autumn leaf fall, which is associated with bud dormancy onset. Accordingly, in Rosaceae, DAMs may play multiple important roles affecting whole plant growth during the tree dormancy induction phase.

The conserved AvrE family of bacterial effectors: functions and targets during pathogenesis

The AvrE family of type III secreted effectors are highly conserved among many agriculturally important phytopathogenic bacteria. Despite their critical roles in the pathogenesis of phytopathogenic bacteria, the molecular functions and virulence mechanisms of these effectors have been largely unknown. However, recent studies have identified host-interacting proteins and demonstrated that AvrE family effectors can form water-permeable channels in the plant plasma membrane (PM) to create a hydrated and nutrient-rich extracellular space (apoplast) required for disease establishment. Here, we summarize these recent discoveries and highlight open questions related to AvrE-targeted host proteins.

The SET-Domain-Containing Protein MdSDG26 Negatively Regulates Alternaria alternata Resistance in Apple

Apple leaf spot is one of the most devastating diseases in the apple industry, caused by Alternaria alternata f. sp mali (A. alternata). SET-domain group (SDG) proteins function as the histone methyltransferases and participate in plant development and stress responses. However, whether SDG proteins are associated with A. alternata resistance is largely unclear. Here, we describe the pathogen-inducible MdSDG26 gene in apple (Malus × domestica). MdSDG26 has two transcript variants that function similarly in catalyzing histone methylation and A. alternata resistance. Transient overexpression of MdSDG26 increased the global levels of H3K4me3 and H3K36me3, whereas knockdown of MdSDG26 only reduced the H3K36me3 level. Transcriptome analysis revealed that MdSDG26 affected the genome-wide transcriptome changes in response to A. alternata infection. ChIP-qPCR analysis demonstrated that MdSDG26 modulates the levels of H3K36me3 and H3K4me3 at both the promoter and exon regions of MdNTL9. As a negative regulator of A. alternata resistance in apples, MdNTL9 plays a pivotal role in MdSDG26-mediated resistance to A. alternata. Therefore, our findings provide compelling evidence for the regulatory function of MdSDG26 in histone methylation and its molecular role in conferring resistance to A. alternata.

Whole-genome resequencing identifies candidate genes and allelic variation in the MdNADP-ME promoter that regulate fruit malate and fructose contents in apple

Soluble sugar and organic acids are key determinants of fruit organoleptic quality and directly affect the commodity value and economic returns of fruit crops. We performed whole-genome sequencing of the apple varieties Gala and Xiahongrou, along with their F1 hybrids, to construct a high-density bin map. Our quantitative genetic analysis pinpointed 53 quantitative trait loci (QTLs) related to 11 sugar and acid traits. We identified a candidate gene, MdNADP-ME, responsible for malate degradation, in a stable QTL on linkage group 15. Sequence analysis revealed an A/C SNP in the promoter region (MEp-799) that influences binding of the MdMYB2 transcription factor, thereby affecting MdNADP-ME expression. In our study of various apple genotypes, this SNP has been demonstrated to be linked to malate and fructose levels. We also developed a dCAPS marker associated with fruit fructose content. These results substantiate the role of MdNADP-ME in maintaining the equilibrium between sugar and acid contents in apple fruits.

Gene Duplication and Functional Diversification of MADS-Box Genes in Malus × domestica following WGD: Implications for Fruit Type and Floral Organ Evolution

The evolution of the MADS-box gene family is essential for the rapid differentiation of floral organs and fruit types in angiosperms. Two key processes drive the evolution of gene families: gene duplication and functional differentiation. Duplicated copies provide the material for variation, while advantageous mutations can confer new functions on gene copies. In this study, we selected the Rosaceae family, which includes a variety of fruit types and flower organs, as well as species that existed before and after whole-genome duplication (WGD). The results indicate that different fruit types are associated with different copies of MADS-box gene family duplications and WGD events. While most gene copies derived from WGD have been lost, MADS-box genes not only retain copies derived from WGD but also undergo further gene duplication. The sequences, protein structures, and expression patterns of these gene copies have undergone significant differentiation. This work provides a clear example of MADS-box genes in the context of gene duplication and functional differentiation, offering new insights into the evolution of fruit types and floral organs.

Movement of ACC oxidase 3 mRNA from seeds to flesh promotes fruit ripening in apple

Xenia, the phenomenon in which the pollen genotype directly affects the phenotypic characteristics of maternal tissues (i.e., fruit ripening), has applications in crop production and breeding. However, the underlying molecular mechanism has yet to be elucidated. Here, we investigated whether mobile mRNAs from the pollen affect the ripening and quality-related characteristics of the fruit using cross-pollination between distinct Malus domestica (apple) cultivars. We demonstrated that hundreds of mobile mRNAs originating from the seeds are delivered to the fruit. We found that the movement of one of these mRNAs, ACC oxidase 3 (MdACO3), is coordinated with fruit ripening. Salicylic acid treatment, which can cause plasmodesmal closure, blocks MdACO3 movement, indicating that MdACO3 transcripts may move through the plasmodesmata. To assess the role of mobile MdACO3 transcripts in apple fruit, we created MdACO3-GFP-expressing apple seeds using MdACO3-GFP-overexpressing pollen for pollination and showed that MdACO3 transcripts in the transgenic seeds move to the flesh, where they promote fruit ripening. Furthermore, we demonstrated that MdACO3 can be transported from the seeds to fruit in the fleshy-fruited species tomato and strawberry. These results underscore the potential of mobile mRNAs from seeds to influence fruit characteristics, providing an explanation for the xenia phenomenon. Notably, our findings highlight the feasibility of leveraging diverse pollen genomic resources, without resorting to genome editing, to improve fruit quality.

Four glycosyltransferase genes are responsible for synthesis and accumulation of different flavonol glycosides in apple tissues

Flavonols are widely synthesized throughout the plant kingdom, playing essential roles in plant physiology and providing unique health benefits for humans. Their glycosylation plays significant role in improving their stability and solubility, thus their accumulation and function. However, the genes encoding the enzymes catalyze this glycosylation remain largely unknown in apple. This study utilized a combination of methods to identify genes encoding such enzymes. Initially, candidate genes were selected based on their potential to encode UDP-dependent glycosyltransferases (UGTs) and their expression patterns in response to light induction. Subsequently, through testing the in vitro enzyme activity of the proteins produced in Escherichia coli cells, four candidates were confirmed to encode a flavonol 3-O-galactosyltransferase (UGT78T6), flavonol 3-O-glucosyltransferase (UGT78S1), flavonol 3-O-xylosyltransferase/arabinosyltransferase (UGT78T5), and flavonol 3-O-rhamnosyltransferase (UGT76AE22), respectively. Further validation of these genes' functions was conducted by modulating their expression levels in stably transformed apple plants. As anticipated, a positive correlation was observed between the expression levels of these genes and the content of specific flavonol glycosides corresponding to each gene. Moreover, overexpression of a flavonol synthase gene, MdFLS, resulted in increased flavonol glycoside content in apple roots and leaves. These findings provide valuable insights for breeding programs aimed at enriching apple flesh with flavonols and for identifying flavonol 3-O-glycosyltransferases of other plant species.

Transcriptomic dataset of Malus domestica young leaves in response to acibenzolar-S-methyl (ASM) and/or nitrogen nutrition

Plant resistance inducers (PRIs) and nitrogen (N) nutrition are both known to affect plant defence but their interaction has not been well described. We addressed this question in apple (Malus domestica) by generating a transcriptomic data set of young leaves from seedlings grown in subirrigation systems allowing variations in nitrate supply as the sole nitrogen source. Plants under three contrasting N status (high; limited for 10 days; or just resupplied after a 12 days limitation) received foliar applications of the chemical elicitor acibenzolar-S-methyl (ASM), a functional analog of salicylic acid, or water. Two days later, the youngest developed leaves were sampled for total RNA extraction and sequencing analysis (RNAseq). The current dataset includes 1) a detailed protocol of plant sample production and 2) transcriptomic profile description of young leaves as normalized counts obtained from sequence mapping against the Malus domestica GDDH13v1.1 reference transcriptome. The raw data files and processed data are available at the Gene Expression Omnibus (GEO) repository under the accession number GSE264541. This dataset is a valuable resource to investigate further the molecular mechanisms underlying the role of nitrogen and/or ASM treatment in Malus domestica.

A Functional InDel in the WRKY10 Promoter Controls the Degree of Flesh Red Pigmentation in Apple

MYB transcription factors have been linked to anthocyanin synthesis and various color phenotypes in plants. In apple, MYB10 confers a red-flesh phenotype due to a minisatellite insertion in its R6 promoter, but R6:MYB10 genotypes exhibit various degrees of red pigmentation in the flesh, suggesting the involvement of other genetic factors. Here, it is shown that MdWRKY10, a transcription factor identified via DNA pull-down trapping, binds to the promoter of MdMYB10 and activates its transcription. MdWRKY10 specifically interacts with the WDR protein MdTTG1 to join the apple MYB-bHLH-WDR (MBW) complex, which significantly enhances its transcriptional activation activity. A 163-bp InDel detected in the promoter region of the alleles of MdWRKY10 in a hybrid population of identical heterozygous genotypes regarding R6 by structural variation analysis, contains a typical W-box element that MdWRKY10 binds to for transactivation. This leads to increased transcript levels of MdWRKY10 and MdMYB10 and enhanced anthocyanin synthesis in the flesh, largely accounting for the various degrees of flesh red pigmentation in the R6 background. These findings reveal a novel regulatory role of the WRKY-containing protein complex in the formation of red flesh apple phenotypes and provide broader insights into the molecular mechanism governing anthocyanin synthesis in plants.

Why we thrive beneath a northern sky - genomic signals of selection in apple for adaptation to northern Sweden

Good understanding of the genomic regions underlying adaptation of apple to boreal climates is needed to facilitate efficient breeding of locally adapted apple cultivars. Proper infrastructure for phenotyping and evaluation is essential for identification of traits responsible for adaptation, and dissection of their genetic composition. However, such infrastructure is costly and currently not available for the boreal zone of northern Sweden. Therefore, we used historical pomological data on climate adaptation of 59 apple cultivars and whole genome sequencing to identify genomic regions that have undergone historical selection among apple cultivars recommended for cultivation in northern Sweden. We found the apple collection to be composed of two ancestral groups that are largely concordant with the grouping into 'hardy' and 'not hardy' cultivars based on the pomological literature. Using a number of genome-wide scans for signals of selection, we obtained strong evidence of positive selection at a genomic region around 29 MbHFTH1 of chromosome 1 among apple cultivars in the 'hardy' group. Using phased genotypic data from the 20 K apple Infinium® SNP array, we identified haplotypes associated with the two cultivar groups and traced transmission of these haplotypes through the pedigrees of some apple cultivars. This demonstrates that historical data from pomological literature can be analyzed by population genomic approaches as a step towards revealing the genomic control of a key property for a horticultural niche market. Such knowledge is needed to facilitate efficient breeding strategies for development of locally adapted apple cultivars in the future. The current study illustrates the response to a very strong selective pressure imposed on tree crops by climatic factors, and the importance of genetic research on this topic and feasibility of breeding efforts in the light of the ongoing climate change.

Bacillus B2 promotes root growth and enhances phosphorus absorption in apple rootstocks by affecting MhMYB15

Due to the chelation of phosphorus in the soil, it becomes unavailable for plant growth and development. The mechanisms by which phosphorus-solubilizing bacteria activate immobilized phosphorus to promote the growth and development of woody plants, as well as the intrinsic molecular mechanisms, are not clear. Through the analysis of microbial communities in the rhizosphere 16S V3-V4 and a homologous gene encoding microbial alkaline phosphomonoesterase (phoD) in phosphate-efficient (PE) and phosphate-inefficient apple rootstocks, it was found that PE significantly enriched beneficial rhizobacteria. The best phosphorus-solubilizing bacteria, Bacillus sp. strain 7DB1 (B2), was isolated, purified, and identified from the rhizosphere soil of PE rootstocks. Incubating with Bacillus B2 into the rhizosphere of apple rootstocks significantly increased the soluble phosphorus and flavonoid content in the rhizosphere soil. Simultaneously, this process stimulates the root development of the rootstocks and enhances plant phosphorus uptake. After root transcriptome sequencing, candidate transcription factor MhMYB15, responsive to Bacillus B2, was identified through heatmap and co-expression network analysis. Yeast one-hybrid, electrophoretic mobility shift assay, and LUC assay confirmed that MhMYB15 can directly bind to the promoter regions of downstream functional genes, including chalcone synthase MhCHS2 and phosphate transporter MhPHT1;15. Transgenic experiments with MhMYB15 revealed that RNAi-MhMYB15 silenced lines failed to induce an increase in flavonoid content and phosphorus levels in the roots under the treatment of Bacillus B2, and plant growth was slower than the control. In conclusion, MhMYB15 actively responds to Bacillus B2, regulating the accumulation of flavonoids and the uptake of phosphorus, thereby influencing plant growth and development.

Genomic prediction and genome-wide association study using combined genotypic data from different genotyping systems: application to apple fruit quality traits

With advances in next-generation sequencing technologies, various marker genotyping systems have been developed for genomics-based approaches such as genomic selection (GS) and genome-wide association study (GWAS). As new genotyping platforms are developed, data from different genotyping platforms must be combined. However, the potential use of combined data for GS and GWAS has not yet been clarified. In this study, the accuracy of genomic prediction (GP) and the detection power of GWAS increased for most fruit quality traits of apples when using combined data from different genotyping systems, Illumina Infinium single-nucleotide polymorphism array and genotyping by random amplicon sequencing-direct (GRAS-Di) systems. In addition, the GP model, which considered the inbreeding effect, further improved the accuracy of the seven fruit traits. Runs of homozygosity (ROH) islands overlapped with the significantly associated regions detected by the GWAS for several fruit traits. Breeders may have exploited these regions to select promising apples by breeders, increasing homozygosity. These results suggest that combining genotypic data from different genotyping platforms benefits the GS and GWAS of fruit quality traits in apples. Information on inbreeding could be beneficial for improving the accuracy of GS for fruit traits of apples; however, further analysis is required to elucidate the relationship between the fruit traits and inbreeding depression (e.g. decreased vigor).

Abscisic acid activates transcription factor module MdABI5-MdMYBS1 during carotenoid-derived apple fruit coloration

Carotenoids are major pigments contributing to fruit coloration. We previously reported that the apple (Malus domestica Borkh.) mutant fruits of "Beni Shogun" and "Yanfu 3" show a marked difference in fruit coloration. However, the regulatory mechanism underlying this phenomenon remains unclear. In this study, we determined that carotenoid is the main factor influencing fruit flesh color. We identified an R1-type MYB transcription factor (TF), MdMYBS1, which was found to be highly associated with carotenoids and abscisic acid (ABA) contents of apple fruits. Overexpression of MdMYBS1 promoted, and silencing of MdMYBS1 repressed, β-branch carotenoids synthesis and ABA accumulation. MdMYBS1 regulates carotenoid biosynthesis by directly activating the major carotenoid biosynthetic genes encoding phytoene synthase (MdPSY2-1) and lycopene β-cyclase (MdLCYb). 9-cis-epoxycarotenoid dioxygenase 1 (MdNCED1) contributes to ABA biosynthesis, and MdMYBS1 enhances endogenous ABA accumulation by activating the MdNCED1 promoter. In addition, the basic leucine zipper domain TF ABSCISIC ACID-INSENSITIVE5 (MdABI5) was identified as an upstream activator of MdMYBS1, which promotes carotenoid and ABA accumulation. Furthermore, ABA promotes carotenoid biosynthesis and enhances MdMYBS1 and MdABI5 promoter activities. Our findings demonstrate that the MdABI5-MdMYBS1 cascade activated by ABA regulates carotenoid-derived fruit coloration and ABA accumulation in apple, providing avenues in breeding and planting for improvement of fruit coloration and quality.

A haplotype-resolved genome assembly of Malus domestica 'Red Fuji'

The 'Red Fuji' apple (Malus domestica), is one of the most important and popular economic crops worldwide in the fruit industry. Using PacBio HiFi long reads and Hi-C reads, we assembled a high-quality haplotype-resolved genome of 'Red Fuji', with sizes of 668.7 and 668.8 Mb, and N50 sizes of 34.1 and 31.4 Mb. About 97.2% of sequences were anchored in 34 chromosomes. We annotated both haploid genomes, identifying a total of 95,439 protein-coding genes in the two haplotype genomes, with 98% functional annotation. The haplotype-resolved genome of 'Red Fuji' apple stands as a precise benchmark for an array of analyses, such as comparative genomics, transcriptomics, and allelic expression studies. This comprehensive resource is paramount in unraveling variations in allelic expression, advancing quality improvements, and refining breeding efforts.

Regulatory role of Prunus mume DAM6 on lipid body accumulation and phytohormone metabolism in the dormant vegetative meristem

Bud dormancy is a crucial process in the annual growth cycle of woody perennials. In Rosaceae fruit tree species, DORMANCY-ASSOCIATED MADS-box (DAM) transcription factor genes regulating bud dormancy have been identified, but their molecular roles in meristematic tissues have not been thoroughly characterized. In this study, molecular and physiological analyses of transgenic apple plants overexpressing the Japanese apricot DAM6 gene (PmDAM6) and Japanese apricot cultivars and F1 individuals with contrasting dormancy characteristics revealed the metabolic pathways controlled by PmDAM6. Our transcriptome analysis and transmission electron microscopy examination demonstrated that PmDAM6 promotes the accumulation of lipid bodies and inhibits cell division in the dormant vegetative meristem by down-regulating the expression of lipid catabolism genes (GDSL ESTERASE/LIPASE and OIL BODY LIPASE) and CYCLIN genes, respectively. Our findings also indicate PmDAM6 promotes abscisic acid (ABA) accumulation and decreases cytokinin (CTK) accumulation in vegetative buds by up-regulating the expression of the ABA biosynthesis gene ARABIDOPSIS ALDEHYDE OXIDASE and the CTK catabolism gene CYTOKININ DEHYDROGENASE, while also down-regulating the expression of the CTK biosynthesis genes ISOPENTENYL TRANSFERASE (IPT) and CYP735A. Additionally, PmDAM6 modulates gibberellin (GA) metabolism by up-regulating GA2-OXIDASE expression and down-regulating GA3-OXIDASE expression. Furthermore, PmDAM6 may also indirectly promote lipid accumulation and restrict cell division by limiting the accumulation of CTK and GA in buds. In conclusion, using our valuable genetic platform, we clarified how PmDAM6 modifies diverse cellular processes, including lipid catabolism, phytohormone (ABA, CTK, and GA) biosynthesis and catabolism, and cell division, in the dormant vegetative meristem.

Transposable elements in Rosaceae: insights into genome evolution, expression dynamics, and syntenic gene regulation

Transposable elements (TEs) exert significant influence on plant genomic structure and gene expression. Here, we explored TE-related aspects across 14 Rosaceae genomes, investigating genomic distribution, transposition activity, expression patterns, and nearby differentially expressed genes (DEGs). Analyses unveiled distinct long terminal repeat retrotransposon (LTR-RT) evolutionary patterns, reflecting varied genome size changes among nine species over the past million years. In the past 2.5 million years, Rubus idaeus showed a transposition rate twice as fast as Fragaria vesca, while Pyrus bretschneideri displayed significantly faster transposition compared with Crataegus pinnatifida. Genes adjacent to recent TE insertions were linked to adversity resistance, while those near previous insertions were functionally enriched in morphogenesis, enzyme activity, and metabolic processes. Expression analysis revealed diverse responses of LTR-RTs to internal or external conditions. Furthermore, we identified 3695 pairs of syntenic DEGs proximal to TEs in Malus domestica cv. 'Gala' and M. domestica (GDDH13), suggesting TE insertions may contribute to varietal trait differences in these apple varieties. Our study across representative Rosaceae species underscores the pivotal role of TEs in plant genome evolution within this diverse family. It elucidates how these elements regulate syntenic DEGs on a genome-wide scale, offering insights into Rosaceae-specific genomic evolution.

Recent Advances in Studies of Genomic DNA Methylation and Its Involvement in Regulating Drought Stress Response in Crops

As global arid conditions worsen and groundwater resources diminish, drought stress has emerged as a critical impediment to plant growth and development globally, notably causing declines in crop yields and even the extinction of certain cultivated species. Numerous studies on drought resistance have demonstrated that DNA methylation dynamically interacts with plant responses to drought stress by modulating gene expression and developmental processes. However, the precise mechanisms underlying these interactions remain elusive. This article consolidates the latest research on the role of DNA methylation in plant responses to drought stress across various species, focusing on methods of methylation detection, mechanisms of methylation pattern alteration (including DNA de novo methylation, DNA maintenance methylation, and DNA demethylation), and overall responses to drought conditions. While many studies have observed significant shifts in genome-wide or gene promoter methylation levels in drought-stressed plants, the identification of specific genes and pathways involved remains limited. This review aims to furnish a reference for detailed research into plant responses to drought stress through epigenetic approaches, striving to identify drought resistance genes regulated by DNA methylation, specific signaling pathways, and their molecular mechanisms of action.

A nearly gapless, highly contiguous reference genome for a doubled haploid line of Populus ussuriensis, enabling advanced genomic studies

Populus species, particularly P. trichocarpa, have long served as model trees for genomics research, owing to fully sequenced genomes. However, the high heterozygosity, and the presence of repetitive regions, including centromeres and ribosomal RNA gene clusters, have left 59 unresolved gaps, accounting for approximately 3.32% of the P. trichocarpa genome. In this study, the callus induction method was improved to derive a doubled haploid (DH) callus line from P. ussuriensis anthers. Leveraging long-read sequencing, we successfully assembled a nearly gap-free, telomere-to-telomere (T2T) P. ussuriensis genome spanning 412.13 Mb. This genome assembly contains only seven gaps and has a contig N50 length of 19.50 Mb. Annotation revealed 34,953 protein-coding genes in this genome, which is 465 more than that of P. trichocarpa. Notably, centromeric regions are characterized by higher-order repeats, we identified and annotated centromere regions in all DH genome chromosomes, a first for poplars. The derived DH genome exhibits high collinearity with P. trichocarpa and significantly fills gaps present in the latter's genome. This T2T P. ussuriensis reference genome will not only enhance our understanding of genome structure, and functions within the poplar genus but also provides valuable resources for poplar genomic and evolutionary studies.

The genetic basis of apple shape and size unraveled by digital phenotyping

Great diversity of shape, size, and skin color is observed among the fruits of different apple genotypes. These traits are critical for consumers and therefore interesting targets for breeding new apple varieties. However, they are difficult to phenotype and their genetic basis, especially for fruit shape and ground color, is largely unknown. We used the FruitPhenoBox to digitally phenotype 525 genotypes of the apple reference population (apple REFPOP) genotyped for 303,148 single nucleotide polymorphism (SNP) markers. From the apple images, 573 highly heritable features describing fruit shape and size as well as 17 highly heritable features for fruit skin color were extracted to explore genotype-phenotype relationships. Out of these features, seven principal components (PCs) and 16 features with the Pearson's correlation r < 0.75 (selected features) were chosen to carry out genome-wide association studies (GWAS) for fruit shape and size. Four PCs and eight selected features were used in GWAS for fruit skin color. In total, 69 SNPs scattered over all 17 apple chromosomes were significantly associated with round, conical, cylindrical, or symmetric fruit shapes and fruit size. Novel associations with major effect on round or conical fruit shapes and fruit size were identified on chromosomes 1 and 2. Additionally, 16 SNPs associated with PCs and selected features related to red overcolor as well as green and yellow ground color were found on eight chromosomes. The identified associations can be used to advance marker-assisted selection in apple fruit breeding to systematically select for desired fruit appearance.

The MdCBF1/2-MdTST1/2 module regulates sugar accumulation in response to low temperature in apple

Soluble sugar content is a key component in controlling fruit flavor, and its accumulation in fruit is largely determined by sugar metabolism and transportation. When the diurnal temperature range is greater, the fleshy fruits accumulated more soluble sugars and become more sweeter. However, the molecular mechanism underlying this response remains largely unknown. In this study, we verified that low-temperature treatment promoted soluble sugar accumulation in apple fruit and found that this was due to the upregulation of the Tonoplast Sugar Transporter genes MdTST1/2. A combined strategy using assay for transposase-accessible chromatin (ATAC) sequencing and gene expression and cis-acting elements analyses, we identified two C-repeat Binding Factors, MdCBF1 and MdCBF2, that were induced by low temperature and that might be upstream transcription factors of MdTST1/2. Further studies established that MdCBF1/2 could bind to the promoters of MdTST1/2 and activate their expression. Overexpression of MdCBF1 or MdCBF2 in apple calli and fruit significantly upregulated MdTST1/2 expression and increased the concentrations of glucose, fructose, and sucrose. Suppression of MdTST1 and/or MdTST2 in an MdCBF1/2-overexpression background abolished the positive effect of MdCBF1/2 on sugar accumulation. In addition, simultaneous silencing of MdCBF1/2 downregulated MdTST1/2 expression and apple fruits failed to accumulate more sugars under low-temperature conditions, indicating that MdCBF1/2-mediated sugar accumulation was dependent on MdTST1/2 expression. Hence, we concluded that the MdCBF1/2-MdTST1/2 module is crucial for sugar accumulation in apples in response to low temperatures. Our findings provide mechanistic components coordinating the relationship between low temperature and sugar accumulation as well as new avenues to improve fruit quality.

Adaptive regulation of miRNAs/milRNAs in tissue-specific interaction between apple and Valsa mali

In plant-pathogen interactions, pathogens display tissue specificity, infecting and causing disease in particular tissues. However, the involvement of microRNAs/microRNA-like RNAs (miRNAs/milRNAs) in tissue-specific regulation during plant-pathogen interactions remains largely unexplored. This study investigates the differential expression of miRNAs/milRNAs, as well as their corresponding target genes, in interactions between Valsa mali (Vm) and different apple tissues. The results demonstrated that both apple miRNAs and Vm milRNAs exhibited distinct expression profiles when Vm infected bark and leaves, with functionally diverse corresponding target genes. Furthermore, one apple miRNA (Mdo-miR482a) and one Vm milRNA (Vm-milR57) were identified as exhibiting tissue-specific expression in interactions between Vm and apple bark or leaves. Mdo-miR482a was exclusively up-regulated in response to Vm infection in bark and target a nucleotide-binding leucine-rich repeat (NLR) gene of apple. When Mdo-miR482a was transiently over-expressed or silenced, the resistance was significantly reduced or improved. Similarly, transient expression of the NLR gene also showed an increase in resistance. Vm-milR57 could target two essential pathogenicity-related genes of Vm. During Vm infection in bark, the expression of Vm-milR57 was down-regulated to enhance the expression of the corresponding target gene to improve the pathogenicity. The study is the first to reveal tissue-specific characteristics of apple miRNAs and Vm milRNAs in interactions between Vm and different apple tissues, providing new insights into adaptive regulation in tissue-specific interactions between plants and fungi.

Identification of two key genes involved in flavonoid catabolism and their different roles in apple resistance to biotic stresses

Biosynthesis of flavonoid aglycones and glycosides is well established. However, key genes involved in their catabolism are poorly understood, even though the products of hydrolysis and oxidation play important roles in plant resistance to biotic stress. Here, we report on catabolism of dihydrochalcones (DHCs), the most abundant flavonoids in domesticated apple and wild Malus. Two key genes, BGLU13.1 and PPO05, were identified by activity-directed protein purification. BGLU13.1-A hydrolyzed phlorizin, (the most abundant DHC in domesticated apple) to produce phloretin which was then oxidized by PPO05. The process differed in some wild Malus, where trilobatin (a positional isomer of phlorizin) was mainly oxidized by PPO05. The effects of DHC catabolism on apple resistance to biotic stresses was investigated using transgenic plants. Either directly or indirectly, phlorizin hydrolysis affected resistance to the phytophagous pest two-spotted spider mite, while oxidation of trilobatin was involved in resistance to the biotrophic fungus Podosphaera leucotricha. DHC catabolism did not affect apple resistance to necrotrophic pathogens Valsa mali and Erwinia amylovara. These results suggest that different DHC catabolism pathways play different roles in apple resistance to biotic stresses. The role of DHC catabolism on apple resistance appeared closely related to the mode of invasion/damage used by pathogen/pest.

MYB transcription factors encoded by diversified tandem gene clusters cause varied Morella rubra fruit color

Chinese bayberry (Morella rubra) is a fruit tree with a remarkable variation in fruit color, ranging from white to dark red as determined by anthocyanin content. In dark red "Biqi" (BQ), red "Dongkui" (DK), pink "Fenhong" (FH), and white "Shuijing" (SJ), we identified an anthocyanin-related MYB transcription factor-encoding gene cluster of four members, i.e. MrMYB1.1, MrMYB1.2, MrMYB1.3, and MrMYB2. Collinear analysis revealed that the MYB tandem cluster may have occurred in a highly conserved region of many eudicot genomes. Two alleles of MrMYB1.1 were observed; MrMYB1.1-1 (MrMYB1.1n) was a full-length allele and homozygous in "BQ", MrMYB1.1-2 (MrMYB1.1d) was a nonfunctional allele with a single base deletion and homozygous in "SJ", and MrMYB1.1n/MrMYB1.1d were heterozygous in "DK" and "FH". In these four cultivars, expression of MrMYB1.1, MrMYB1.2, and MrMYB2 was enhanced during ripening. Both alleles were equally expressed in MrMYB1.1n/MrMYB1.1d heterozygous cultivars as revealed by a cleaved amplified polymorphic sequence marker. Expression of MrMYB1.3 was restricted to some dark red cultivars only. Functional characterization revealed that MrMYB1.1n and MrMYB1.3 can induce anthocyanin accumulation while MrMYB1.1d, MrMYB1.2, and MrMYB2 cannot. DNA-protein interaction assays indicated that MrMYB1.1n and MrMYB1.3 can directly bind to and activate the promoters of anthocyanin-related genes via interaction with a MYC-like basic helix-loop-helix protein MrbHLH1. We concluded that the specific genotype of MrMYB1.1 alleles, as well as the exclusive expression of MrMYB1.3 in some dark red cultivars, contributes to fruit color variation. The study provides insights into the mechanisms for regulation of plant anthocyanin accumulation by MYB tandem clusters.

Phylogenomic profiles of whole-genome duplications in Poaceae and landscape of differential duplicate retention and losses among major Poaceae lineages

Poaceae members shared a whole-genome duplication called rho. However, little is known about the evolutionary pattern of the rho-derived duplicates among Poaceae lineages and implications in adaptive evolution. Here we present phylogenomic/phylotranscriptomic analyses of 363 grasses covering all 12 subfamilies and report nine previously unknown whole-genome duplications. Furthermore, duplications from a single whole-genome duplication were mapped to multiple nodes on the species phylogeny; a whole-genome duplication was likely shared by woody bamboos with possible gene flow from herbaceous bamboos; and recent paralogues of a tetraploid Oryza are implicated in tolerance of seawater submergence. Moreover, rho duplicates showing differential retention among subfamilies include those with functions in environmental adaptations or morphogenesis, including ACOT for aquatic environments (Oryzoideae), CK2β for cold responses (Pooideae), SPIRAL1 for rapid cell elongation (Bambusoideae), and PAI1 for drought/cold responses (Panicoideae). This study presents a Poaceae whole-genome duplication profile with evidence for multiple evolutionary mechanisms that contribute to gene retention and losses.

Phenomics and transcriptomic profiling of fruit development in distinct apple varieties

Apple is one of the most economically important and popular temperate fruit trees. The domestication of apple has resulted in substantial phenotypic differences, particularly between wild and cultivated varieties. However, the relationship between gene expression and phenotypic variations in apple remains poorly understood. Here, we present a comprehensive dataset featuring five distinct apple varieties, including two wild varieties and three representative cultivated varieties. The dataset comprises of both phenomics data, encompassing twelve fruit quality-related traits continuously measured over two years, and transcriptomic data obtained at different developmental stages with three biological replicates. We performed basic quality control process, gene expression normalization and differential gene expression analysis to demonstrate the utility and reliability of the dataset. Our findings indicate that gene expression strongly related with phenotypic variations in apple. This dataset serves as a valuable resource, encompassing phenomics and transcriptomic data in multiple formats, thereby facilitating further exploration of the relationships between gene expression and phenotypic traits in apple.

Phaseolus vulgaris STP13.1 is an H+-coupled monosaccharide transporter, present in source leaves and seed coats, with higher substrate affinity at depolarized potentials

Sugar transport proteins (STPs) are high-affinity H+-coupled hexose symporters. Recently, the contribution of STP13 to bacterial and fungal pathogen resistance across multiple plant species has garnered significant interest. Quantitative PCR analysis of source leaves, developing embryos, and seed coats of Phaseolus vulgaris L. (common bean) revealed that PvSTP13.1 was expressed in source leaves and seed coats throughout seed development. In contrast, PvSTP13.1 transcripts were detected at exceedingly low levels in developing embryos. To characterize the transport mechanism, PvSTP13.1 was expressed in Xenopus laevis oocytes, and inward-directed currents were analyzed using two-electrode voltage clamping. PvSTP13.1 was shown to function as an H+-coupled monosaccharide symporter exhibiting a unique high affinity for hexoses and aldopentoses at depolarized membrane potentials. Specifically, of the 31 assessed substrates, which included aldohexoses, deoxyhexoses, fructose, 3-O-methyl-D-glucose, aldopentoses, polyols, glycosides, disaccharides, trisaccharides, and glucuronic acid, PvSTP13.1 displayed the highest affinity (K 0.5) for glucose (43 μM), mannose (92 μM), galactose (145 μM), fructose (224 μM), xylose (1.0 mM), and fucose (3.7 mM) at pH 5.6 at a depolarized membrane potential of -40 mV. The results presented here suggest PvSTP13.1 contributes to retrieval of hexoses from the apoplasmic space in source leaves and coats of developing seeds.

Transcriptomic investigation of the interaction between a biocontrol yeast, Papiliotrema terrestris strain PT22AV, and the postharvest fungal pathogen Penicillium expansum on apple

Biocontrol strategies offer a promising alternative to control plant pathogens achieving food safety and security. In this study we apply a RNAseq analysis during interaction between the biocontrol agent (BCA) Papiliotrema terrestris, the pathogen Penicillium expansum, and the host Malus domestica. Analysis of the BCA finds overall 802 upregulated DEGs (differentially expressed genes) when grown in apple tissue, with the majority being involved in nutrients uptake and oxidative stress response. This suggests that these processes are crucial for the BCA to colonize the fruit wounds and outcompete the pathogen. As to P. expansum analysis, 1017 DEGs are upregulated when grown in apple tissue, with the most represented GO categories being transcription, oxidation reduction process, and transmembrane transport. Analysis of the host M. domestica finds a higher number of DEGs in response to the pathogen compared to the BCA, with overexpression of genes involved in host defense signaling pathways in the presence of both of them, and a prevalence of pattern-triggered immunity (PTI) and effector-triggered immunity (ETI) only during interaction with P. expansum. This analysis contributes to advance the knowledge on the molecular mechanisms that underlie biocontrol activity and the tritrophic interaction of the BCA with the pathogen and the host.

The structure, function, and evolution of plant centromeres

Centromeres are essential regions of eukaryotic chromosomes responsible for the formation of kinetochore complexes, which connect to spindle microtubules during cell division. Notably, although centromeres maintain a conserved function in chromosome segregation, the underlying DNA sequences are diverse both within and between species and are predominantly repetitive in nature. The repeat content of centromeres includes high-copy tandem repeats (satellites), and/or specific families of transposons. The functional region of the centromere is defined by loading of a specific histone 3 variant (CENH3), which nucleates the kinetochore and shows dynamic regulation. In many plants, the centromeres are composed of satellite repeat arrays that are densely DNA methylated and invaded by centrophilic retrotransposons. In some cases, the retrotransposons become the sites of CENH3 loading. We review the structure of plant centromeres, including monocentric, holocentric, and metapolycentric architectures, which vary in the number and distribution of kinetochore attachment sites along chromosomes. We discuss how variation in CENH3 loading can drive genome elimination during early cell divisions of plant embryogenesis. We review how epigenetic state may influence centromere identity and discuss evolutionary models that seek to explain the paradoxically rapid change of centromere sequences observed across species, including the potential roles of recombination. We outline putative modes of selection that could act within the centromeres, as well as the role of repeats in driving cycles of centromere evolution. Although our primary focus is on plant genomes, we draw comparisons with animal and fungal centromeres to derive a eukaryote-wide perspective of centromere structure and function.

Multi-omics provide insights into the regulation of DNA methylation in pear fruit metabolism

BACKGROUND

Extensive research has been conducted on fruit development in crops, but the metabolic regulatory networks underlying perennial fruit trees remain poorly understood. To address this knowledge gap, we conduct a comprehensive analysis of the metabolome, proteome, transcriptome, DNA methylome, and small RNAome profiles of pear fruit flesh at 11 developing stages, spanning from fruitlet to ripening. Here, we systematically investigate the metabolic landscape and regulatory network involved.

RESULTS

We generate an association database consisting of 439 metabolites and 14,399 genes to elucidate the gene regulatory network of pear flesh metabolism. Interestingly, we detect increased DNA methylation in the promoters of most genes within the database during pear flesh development. Application of a DNA methylation inhibitor to the developing fruit represses chlorophyll degradation in the pericarp and promotes xanthophyll, β-carotene, and abscisic acid (ABA) accumulation in the flesh. We find the gradual increase in ABA production during pear flesh development is correlated with the expression of several carotenoid pathway genes and multiple transcription factors. Of these transcription factors, the zinc finger protein PbZFP1 is identified as a positive mediator of ABA biosynthesis in pear flesh. Most ABA pathway genes and transcription factors are modified by DNA methylation in the promoters, although some are induced by the DNA methylation inhibitor. These results suggest that DNA methylation inhibits ABA accumulation, which may delay fruit ripening.

CONCLUSION

Our findings provide insights into epigenetic regulation of metabolic regulatory networks during pear flesh development, particularly with regard to DNA methylation.

Methylomes as key features for predicting recombination in some plant species

Knowing how chromosome recombination works is essential for plant breeding. It enables the design of crosses between different varieties to combine desirable traits and create new ones. This is because the meiotic crossovers between homologous chromatids are not purely random, and various strategies have been developed to describe and predict such exchange events. Recent studies have used methylation data to predict chromosomal recombination in rice using machine learning models. This approach proved successful due to the presence of a positive correlation between the CHH context cytosine methylation and recombination rates in rice chromosomes. This paper assesses the question if methylation can be used to predict recombination in four plant species: Arabidopsis, maize, sorghum, and tomato. The results indicate a positive association between CHH context methylation and recombination rates in certain plant species, with varying degrees of strength in their relationships. The CG and CHG methylation contexts show negative correlation with recombination. Methylation data was key effectively in predicting recombination in sorghum and tomato, with a mean determination coefficient of 0.65 ± 0.11 and 0.76 ± 0.05, respectively. In addition, the mean correlation values between predicted and experimental recombination rates were 0.83 ± 0.06 for sorghum and 0.90 ± 0.05 for tomato, confirming the significance of methylomes in both monocotyledonous and dicotyledonous species. The predictions for Arabidopsis and maize were not as accurate, likely due to the comparatively weaker relationships between methylation contexts and recombination, in contrast to sorghum and tomato, where stronger associations were observed. To enhance the accuracy of predictions, further evaluations using data sets closely related to each other might prove beneficial. In general, this methylome-based method holds great potential as a reliable strategy for predicting recombination rates in various plant species, offering valuable insights to breeders in their quest to develop novel and improved varieties.

Towards identification of postharvest fruit quality transcriptomic markers in Malus domestica

Gene expression is highly impacted by the environment and can be reflective of past events that affected developmental processes. It is therefore expected that gene expression can serve as a signal of a current or future phenotypic traits. In this paper we identify sets of genes, which we call Prognostic Transcriptomic Biomarkers (PTBs), that can predict firmness in Malus domestica (apple) fruits. In apples, all individuals of a cultivar are clones, and differences in fruit quality are due to the environment. The apples transcriptome responds to these differences in environment, which makes PTBs an attractive predictor of future fruit quality. PTBs have the potential to enhance supply chain efficiency, reduce crop loss, and provide higher and more consistent quality for consumers. However, several questions must be addressed. In this paper we answer the question of which of two common modeling approaches, Random Forest or ElasticNet, outperforms the other. We answer if PTBs with few genes are efficient at predicting traits. This is important because we need few genes to perform qPCR, and we answer the question if qPCR is a cost-effective assay as input for PTBs modeled using high-throughput RNA-seq. To do this, we conducted a pilot study using fruit texture in the 'Gala' variety of apples across several postharvest storage regiments. Fruit texture in 'Gala' apples is highly controllable by post-harvest treatments and is therefore a good candidate to explore the use of PTBs. We find that the RandomForest model is more consistent than an ElasticNet model and is predictive of firmness (r2 = 0.78) with as few as 15 genes. We also show that qPCR is reasonably consistent with RNA-seq in a follow up experiment. Results are promising for PTBs, yet more work is needed to ensure that PTBs are robust across various environmental conditions and storage treatments.

Introduction of a diverse genetic background of Pyrus into Malus through intergeneric hybridization

Wide hybridizations across species and genera have been employed to enhance agriculturally important traits in crops. Within the tribe Maleae of the Rosaceae family, different genera and species exhibit several traits useful for increasing diversity and gene pool through hybridization. This study aimed to develop and characterize intergeneric hybrid individuals between Malus and Pyrus. Through seed germination, shoot multiplication, and rooting in vitro, acclimatized seedlings showing vegetative growth on their own roots were obtained from crosses of Malus × domestica pollinated by Pyrus communis, P. bretschneideri, and the Pyrus interspecific hybrid (P. communis × P. pyrifolia). Comparative analysis of leaf morphology, flow cytometry, and molecular genotyping confirmed the hybrid status of the individuals. Genome-wide genotyping revealed that all the hybrid individuals inherited genomic fragments symmetrically from the Malus and Pyrus parents. To the best of our knowledge, this is the first report on the development of intergeneric hybrid seedlings between Malus × domestica and P. bretschneideri. Furthermore, the Pyrus interspecific hybrid individual served as a bridge plant for introducing the genetic background of P. pyrifolia into Malus × domestica. The results of this study provided a crucial foundation for breeding through intergeneric hybridization between Malus and Pyrus, facilitating the incorporation of valuable traits from diverse gene pools.

Near-gapless and haplotype-resolved apple genomes provide insights into the genetic basis of rootstock-induced dwarfing

Dwarfing rootstocks have transformed the production of cultivated apples; however, the genetic basis of rootstock-induced dwarfing remains largely unclear. We have assembled chromosome-level, near-gapless and haplotype-resolved genomes for the popular dwarfing rootstock 'M9', the semi-vigorous rootstock 'MM106' and 'Fuji', one of the most commonly grown apple cultivars. The apple orthologue of auxin response factor 3 (MdARF3) is in the Dw1 region of 'M9', the major locus for rootstock-induced dwarfing. Comparing 'M9' and 'MM106' genomes revealed a 9,723-bp allele-specific long terminal repeat retrotransposon/gypsy insertion, DwTE, located upstream of MdARF3. DwTE is cosegregated with the dwarfing trait in two segregating populations, suggesting its prospective utility in future dwarfing rootstock breeding. In addition, our pipeline discovered mobile mRNAs that may contribute to the development of dwarfed scion architecture. Our research provides valuable genomic resources and applicable methodology, which have the potential to accelerate breeding dwarfing rootstocks for apple and other perennial woody fruit trees.

Identification of Hsp20 gene family in Malus domestica and functional characterization of Hsp20 class I gene MdHsp18.2b

Heat shock protein 20 (Hsp20) is a small molecule heat shock protein that plays an important role in plant growth, development, and stress resistance. Little is known about the function of Hsp20 family genes in apple (Malus domestica). Here, we performed a genome-wide analysis of the apple Hsp20 gene family, and a total of 49 Hsp20s genes were identified from the apple genome. Phylogenetic analysis revealed that the 49 genes were divided into 11 subfamilies, and MdHsp18.2b, a member located in the CI branch, was selected as a representative member for functional characterization. Treatment with NaCl and Botryosphaeria dothidea (B. dothidea), the causal agent of apple ring rot disease, significantly induced MdHsp18.2b transcription level. Further analysis revealed that overexpressing MdHsp18.2b reduced the resistance to salt stress but enhanced the resistance to B. dothidea infection in apple calli. Moreover, MdHsp18.2b positively regulated anthocyanin accumulation in apple calli. Physiology assays revealed that MdHsp18.2b promoted H2O2 production, even in the absence of stress factors, which might contribute to its functions in response to NaCl and B. dothidea infection. Hsps usually function as homo- or heterooligomers, and we found that MdHsp18.2b could form a heterodimer with MdHsp17.9a and MdHsp17.5, two members from the same branch with MdHsp18.2b in the phylogenetic tree. Therefore, we identified 49 Hsp20s genes from the apple genome and found that MdHsp18.2b was involved in regulating plant resistance to salt stress and B. dothidea infection, as well as in regulating anthocyanin accumulation in apple calli.

Epigenomic mechanism regulating the quality and ripeness of apple fruit with differing harvest maturity

Harvest maturity significantly affects the quality of apple fruit in post-harvest storage process. Although the regulatory mechanisms underlying fruit ripening have been studied, the associated epigenetic modifications remain unclear. Thus, we compared the DNA methylation changes and the transcriptional responses of mature fruit (MF) and immature fruit (NF). There were significant correlations between DNA methylation and gene expression. Moreover, the sugar contents (sucrose, glucose, and fructose) were higher in MF than in NF, whereas the opposite pattern was detected for the starch content. The expression-level differences were due to DNA methylations and ultimately resulted in diverse fruit textures and ripeness. Furthermore, the higher ethylene, auxin, and abscisic acid levels in MF than in NF, which influenced the fruit texture and ripening, were associated with multiple differentially expressed genes in hormone synthesis, signaling, and response pathways (ACS, ACO, ZEP, NCED, and ABA2) that were regulated by DNA methylations. Multiple transcription factor genes involved in regulating fruit ripening and quality via changes in DNA methylation were identified, including MIKCC-type MADS-box genes and fruit ripening-related genes (NAP, SPL, WRKY, and NAC genes). These findings reflect the diversity in the epigenetic regulation of gene expression and may be relevant for elucidating the epigenetic regulatory mechanism underlying the ripening and quality of apple fruit with differing harvest maturity.

Constitutive expression of apple endo-POLYGALACTURONASE1 in fruit induces early maturation, alters skin structure and accelerates softening

During fruit ripening, polygalacturonases (PGs) are key contributors to the softening process in many species. Apple is a crisp fruit that normally exhibits only minor changes to cell walls and limited fruit softening. Here, we explore the effects of PG overexpression during fruit development using transgenic apple lines overexpressing the ripening-related endo-POLYGALACTURONASE1 gene. MdPG1-overexpressing (PGox) fruit displayed early maturation/ripening with black seeds, conversion of starch to sugars and ethylene production occurring by 80 days after pollination (DAP). PGox fruit exhibited a striking, white-skinned phenotype that was evident from 60 DAP and most likely resulted from increased air spaces and separation of cells in the hypodermis due to degradation of the middle lamellae. Irregularities in the integrity of the epidermis and cuticle were also observed. By 120 DAP, PGox fruit cracked and showed lenticel-associated russeting. Increased cuticular permeability was associated with microcracks in the cuticle around lenticels and was correlated with reduced cortical firmness at all time points and extensive post-harvest water loss from the fruit, resulting in premature shrivelling. Transcriptomic analysis suggested that early maturation was associated with upregulation of genes involved in stress responses, and overexpression of MdPG1 also altered the expression of genes involved in cell wall metabolism (e.g. β-galactosidase, MD15G1221000) and ethylene biosynthesis (e.g. ACC synthase, MD14G1111500). The results show that upregulation of PG not only has dramatic effects on the structure of the fruit outer cell layers, indirectly affecting water status and turgor, but also has unexpected consequences for fruit development.

Genome assembly provides insights into the genome evolution of Baccaurea ramiflora Lour

Baccaurea ramiflora Lour., an evergreen tree of the Baccaurea genus of the Phyllanthaceae family, is primarily distributed in South Asia, Southeast Asia, and southern China, including southern Yunnan Province. It is a wild or semi-cultivated tree species with ornamental, edible, and medicinal value, exhibiting significant development potential. In this study, we present the whole-genome sequencing of B. ramiflora, employing a combination of PacBio SMRT and Illumina HiSeq 2500 sequencing techniques. The assembled genome size was 975.8 Mb, with a contig N50 of 509.33 kb and the longest contig measuring 7.74 Mb. The genome comprises approximately 73.47% highly repetitive sequences, of which 52.1% are long terminal repeat-retrotransposon sequences. A total of 29,172 protein-coding genes were predicted, of which 25,980 (89.06%) have been annotated, Additionally, 3452 non-coding RNAs were identified. Comparative genomic analysis revealed a close relationship between B. ramiflora and the Euphorbiaceae family, with both being sister groups that diverged approximately 59.9 million years ago. During the evolutionary process, B. ramiflora exhibited positive selection in 278 candidate genes. Synonymous substitution rate and collinearity analysis demonstrated that B. ramiflora underwent a single ancient genome-wide triploidization event, without recent genome-wide duplication events. This high-quality B. ramiflora genome provides a valuable resource for basic research and tree improvement programs focusing on the Phyllanthaceae family.

The transcription factor MdBPC2 alters apple growth and promotes dwarfing by regulating auxin biosynthesis

Auxin plays important roles throughout plant growth and development. However, the mechanisms of auxin regulation of plant structure are poorly understood. In this study, we identified a transcription factor (TF) of the BARLEY B RECOMBINANT/BASIC PENTACYSTEINE (BBR/BPC) family in apple (Malus × domestica), MdBPC2. It was highly expressed in dwarfing rootstocks, and it negatively regulated auxin biosynthesis. Overexpression of MdBPC2 in apple decreased plant height, altered leaf morphology, and inhibited root system development. These phenotypes were due to reduced auxin levels and were restored reversed after exogenous indole acetic acid (IAA) treatment. Silencing of MdBPC2 alone had no obvious phenotypic effect, while silencing both Class I and Class II BPCs in apple significantly increased auxin content in plants. Biochemical analysis demonstrated that MdBPC2 directly bound to the GAGA-rich element in the promoters of the auxin synthesis genes MdYUC2a and MdYUC6b, inhibiting their transcription and reducing auxin accumulation in MdBPC2 overexpression lines. Further studies established that MdBPC2 interacted with the polycomb group (PcG) protein LIKE HETEROCHROMATIN PROTEIN 1 (LHP1) to inhibit MdYUC2a and MdYUC6b expression via methylation of histone 3 lysine 27 (H3K27me3). Silencing MdLHP1 reversed the negative effect of MdBPC2 on auxin accumulation. Our results reveal a dwarfing mechanism in perennial woody plants involving control of auxin biosynthesis by a BPC transcription factor, suggesting its use for genetic improvement of apple rootstock.

Target enrichment sequencing coupled with GWAS identifies MdPRX10 as a candidate gene in the control of budbreak in apple

The timing of floral budbreak in apple has a significant effect on fruit production and quality. Budbreak occurs as a result of a complex molecular mechanism that relies on accurate integration of external environmental cues, principally temperature. In the pursuit of understanding this mechanism, especially with respect to aiding adaptation to climate change, a QTL at the top of linkage group (LG) 9 has been identified by many studies on budbreak, but the genes underlying it remain elusive. Here, together with a dessert apple core collection of 239 cultivars, we used a targeted capture sequencing approach to increase SNP resolution in apple orthologues of known or suspected A. thaliana flowering time-related genes, as well as approximately 200 genes within the LG9 QTL interval. This increased the 275 223 SNP Axiom® Apple 480 K array dataset by an additional 40 857 markers. Robust GWAS analyses identified MdPRX10, a peroxidase superfamily gene, as a strong candidate that demonstrated a dormancy-related expression pattern and down-regulation in response to chilling. In-silico analyses also predicted the residue change resulting from the SNP allele associated with late budbreak could alter protein conformation and likely function. Late budbreak cultivars homozygous for this SNP allele also showed significantly up-regulated expression of C-REPEAT BINDING FACTOR (CBF) genes, which are involved in cold tolerance and perception, compared to reference cultivars, such as Gala. Taken together, these results indicate a role for MdPRX10 in budbreak, potentially via redox-mediated signaling and CBF gene regulation. Moving forward, this provides a focus for developing our understanding of the effects of temperature on flowering time and how redox processes may influence integration of external cues in dormancy pathways.