Apple fruit acidity is genetically diversified by natural variations in three hierarchical epistatic genes: MdSAUR37, MdPP2CH and MdALMTII

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).

The MdMYB44-MdTPR1 repressive complex inhibits MdCCD4 and MdCYP97A3 expression through histone deacetylation to regulate carotenoid biosynthesis in apple

Carotenoids are photosynthetic pigments and antioxidants that contribute to different plant colors. However, the involvement of TOPLESS (TPL/TPR)-mediated histone deacetylation in the modulation of carotenoid biosynthesis through ethylene-responsive element-binding factor-associated amphiphilic repression (EAR)-containing transcription factors (TFs) in apple (Malus domestica Borkh.) is poorly understood. MdMYB44 is a transcriptional repressor that contains an EAR repression motif. In the present study, we used functional analyses and molecular assays to elucidate the molecular mechanisms through which MdMYB44-MdTPR1-mediated histone deacetylation influences carotenoid biosynthesis in apples. We identified two carotenoid biosynthetic genes, MdCCD4 and MdCYP97A3, that were confirmed to be involved in MdMYB44-mediated carotenoid biosynthesis. MdMYB44 enhanced β-branch carotenoid biosynthesis by repressing MdCCD4 expression, whereas MdMYB44 suppressed lutein level by repressing MdCYP97A3 expression. Moreover, MdMYB44 partially influences carotenoid biosynthesis by interacting with the co-repressor TPR1 through the EAR motif to inhibit MdCCD4 and MdCYP97A3 expression via histone deacetylation. Our findings indicate that the MdTPR1-MdMYB44 repressive cascade regulates carotenoid biosynthesis, providing profound insights into the molecular basis of histone deacetylation-mediated carotenoid biosynthesis in plants. These results also provide evidence that the EAR-harboring TF/TPL repressive complex plays a universal role in histone deacetylation-mediated inhibition of gene expression in various plants.

MdMYB66 Is Associated with Anthocyanin Biosynthesis via the Activation of the MdF3H Promoter in the Fruit Skin of an Apple Bud Mutant

Skin color is an important trait that is mainly determined by the content and composition of anthocyanins in apples. In this study, a new bud mutant (RM) from 'Oregon Spur II' (OS) of Red Delicious apple was obtained to reveal the mechanism underlying red color formation. Results showed that the total anthocyanin content in RM was significantly higher than that in OS with the development of fruit. Through widely-targeted metabolomics, we found that cyanidin-3-O-galactoside was significantly accumulated in the fruit skin of RM. Transcriptome analysis revealed that the structural gene MdF3H and MdMYB66 transcription factor were significantly up-regulated in the mutant. Overexpression of MdMYB66 in apple fruit and apple callus significantly promoted anthocyanin accumulation and significantly increased the expression level of MdMYB66 and structural genes related to anthocyanin synthesis. Y1H and LUC analysis verified that MdMYB66 could specifically bind to the promoter of MdF3H. The results of the double luciferase activity test showed that MdMYB66 activated MdF3H 3.8 times, which led to increased anthocyanin contents. This might explain the phenotype of red color in RM at the early stage. Taken together, these results suggested that MdMYB66 was involved in regulating the anthocyanin metabolic pathways through precise regulation of gene expression. The functional characterization of MdMYB66 provides insight into the biosynthesis and regulation of anthocyanins.

MdAIL5 overexpression promotes apple adventitious shoot regeneration by regulating hormone signaling and activating the expression of shoot development-related genes

Adventitious shoot (AS) regeneration is a significant factor in the genetic transformation of horticultural plants. It is also a noteworthy approach to their vegetative propagation. AS regeneration remains highly dependent on the genotype or maturity of explants. We here found that the AS regeneration abilities of apple leaves were positively correlated with MdAIL5 expression. MdAIL5 overexpression dramatically increased AS regeneration efficiency. Notably, MdAIL5 overexpression could restore the AS formation ability of explants to a certain extent, which was lost with an increase in maturity. Endogenous hormone detection revealed that MdAIL5 overexpression changed the contents of auxin, cytokinin (CK), and other hormones in apple leaves. Transcriptome analysis revealed that many genes related to auxin, CK, and brassinolide signaling pathways were significantly and differentially expressed between MdAIL5-overexpressing transgenic apple and wild-type apple plants. Yeast one-hybrid assays, the electrophoretic mobility shift assay, and the dual-luciferase reporter assay revealed that MdAIL5 directly binds to MdARF9 and MdHB14 promoters and positively affects their expression. We here established a model of MdAIL5 regulating AS formation, which acts as a theoretical basis for facilitating genotype- or explant maturity-independent AS regeneration in the future.

Main flavor compounds and molecular regulation mechanisms in fruits and vegetables

Fruits and vegetables (F&V) are an indispensable part of a healthy diet. The volatile and nonvolatile compounds present in F&V constitute unique flavor substances. This paper reviews the main flavor substances present in F&V, as well as the biosynthetic pathways and molecular regulation mechanisms of these compounds. A series of compounds introduced include aromatic substances, soluble sugars and organic acids, which constitute the key flavor substances of F&V. Esters, phenols, alcohols, amino acids and terpenes are the main volatile aromatic substances, and nonvolatile substances are represented by amino acids, fatty acids and carbohydrates; The combination of these ingredients is the cause of the sour, sweet, bitter, astringent and spicy taste of these foods. This provides a theoretical basis for the study of the interaction between volatile and nonvolatile substances in F&V, and also provides a research direction for the healthy development of food in the future.

Abscisic acid and regulation of the sugar transporter gene MdSWEET9b promote apple sugar accumulation

Enhancing fruit sugar contents, especially for high-flavonoid apples with a sour taste, is one of the main goals of horticultural crop breeders. This study analyzed sugar accumulation and the underlying mechanisms in the F2 progenies of a hybridization between the high-sugar apple (Malus × domestica) variety "Gala" and high-flavonoid apple germplasm "CSR6R6". We revealed that MdSWEET9b (sugars will eventually be exported transporter) helps mediate sugar accumulation in fruits. Functional characterization of MdSWEET9b in yeast mutants lacking sugar transport as well as in overexpressing and CRISPR/Cas9 knockdown apple calli revealed MdSWEET9b could transport sucrose specifically, ultimately promoting normal yeast growth and accumulation of total sugar contents. Moreover, MdWRKY9 bound to the MdSWEET9b promoter and regulated its activity, which responded to abscisic acid (ABA) signaling. Furthermore, MdWRKY9 interacted with MdbZIP23 (basic leucine zipper) and MdbZIP46, key ABA signal transducers, at the protein and DNA levels to enhance its regulatory effect on MdSWEET9b expression, thereby influencing sugar accumulation. Based on the contents of ABA in lines with differing sugar contents and the effects of ABA treatments on fruits and calli, we revealed ABA as one of the main factors responsible for the diversity in apple fruit sugar content. The results of this study have clarified how MdSWEET9b influences fruit sugar accumulation, while also further elucidating the regulatory effects of the ABA-signaling network on fruit sugar accumulation. This work provides a basis for future explorations of the crosstalk between hormone and sugar metabolism pathways.

ABRE-BINDING FACTOR3-WRKY DNA-BINDING PROTEIN44 module promotes salinity-induced malate accumulation in pear

Malate impacts fruit acidity and plays a vital role in stress tolerance. Malate accumulation is induced by salinity in various plants as a metabolite in coping with this stress. However, the exact molecular mechanism responsible for salinity-induced malate accumulation remains unclear. Here, we determined that salinity treatment induces malate accumulation in pear (Pyrus spp.) fruit, calli, and plantlets compared to the control. Genetic and biochemical analyses established the key roles of PpWRKY44 and ABRE-BINDING FACTOR3 (PpABF3) transcription factors in promoting malate accumulation in response to salinity. We found that PpWRKY44 is involved in salinity-induced malate accumulation by directly binding to a W-box on the promoter of the malate-associated gene aluminum-activated malate transporter 9 (PpALMT9) to activate its expression. A series of in-vivo and in-vitro assays revealed that the G-box cis-element in the promoter of PpWRKY44 was targeted by PpABF3, which further enhanced salinity-induced malate accumulation. Taken together, these findings suggest that PpWRKY44 and PpABF3 play positive roles in salinity-induced malate accumulation in pears. This research provides insights into the molecular mechanism by which salinity affects malate accumulation and fruit quality.

Physiological and transcriptome analyses of the effects of excessive water deficit on malic acid accumulation in apple

Acidity is a determinant of the organoleptic quality of apple, whereas its regulatory mechanism under water stress remains obscure. Fruit from apple 'Yanfu 3' of Fuji trees grown under normal water irrigation (CK), excessive water deficit treatment (DRT) and excessive water irrigation treatment (WAT) were sampled at 85, 100, 115, 130, 145, 160 and 175 days after full bloom designated stages S1, S2, S3, S4, S5, S6 and S7, respectively. DRT treatment reduced the individual fruit weight and fruit moisture content, and increased fruit firmness. The malate content of DRT treatment was higher than that of CK and WAT from stages S1 to S7. RNA sequencing (RNA-seq) analysis of the transcriptome at stages S4, S6 and S7 indicated that malate anabolism was associated with cysteine and methionine, auxin signaling, glyoxylate and dicarboxylate and pyruvate metabolism. Overexpression of MdPEPC4 increased the malate content in apple calli induced by 4% PEG. Our study provides novel insights into the effects of water stress on the molecular mechanism underlying apple fruit acidity.

Organic Acid Accumulation and Associated Dynamic Changes in Enzyme Activity and Gene Expression during Fruit Development and Ripening of Common Loquat and Its Interspecific Hybrid

Loquats have gained increasing attention from consumers and growers for their essential nutrients and unusual phenology, which could help plug a gap period at market in early spring. Fruit acid is a critical contributor to fruit quality. The dynamic changes in organic acid (OA) during fruit development and ripening of common loquat (Dawuxing, DWX) and its interspecific hybrid (Chunhua, CH) were compared, as well as the corresponding enzyme activity and gene expression. At harvest, titratable acid was significantly lower (p ≤ 0.01) in CH (0.11%) than in DWX loquats (0.35%). As the predominant OA compound, malic acid accounted for 77.55% and 48.59% of the total acid of DWX and CH loquats at harvest, followed by succinic acid and tartaric acid, respectively. PEPC and NAD-MDH are key enzymes that participate in malic acid metabolism in loquat. The OA differences in DWX loquat and its interspecific hybrid could be attributed to the coordinated regulation of multiple genes and enzymes associated with OA biosynthesis, degradation, and transport. The data obtained in this work will serve as a fundamental and important basis for future loquat breeding programs and even for improvements in loquat cultural practices.

Transcriptomic and Metabolomic Analysis Reveals a Protein Module Involved in Pre-harvest Apple Peel Browning

Peel browning is a natural phenomenon that adversely affects the appearance of fruits. Research on the regulation of browning in apples (Malus × domestica Borkh.) has mainly focused on post-harvest storage, while studies at the pre-harvest stage are relatively rare. Apple is an economically important horticultural crop prone to peel browning during growth, especially when the fruits are bagged (dark conditions). The present study's integrated transcriptomics and metabolomics analysis revealed that pre-harvest apple peel browning was primarily due to changes in phenolics and flavonoids. The detailed analysis identified MdLAC7's (laccase 7) role in the pre-harvest apple peel browning process. Transient injection, overexpression, and CRISPR/Cas9 knockout of the MdLAC7 gene in apple fruit and calli identified vallinic acid, anthocyanidin, tannic acid, sinapic acid, and catechinic acid as its catalytic substrates. In addition, yeast one-hybrid (Y1H) assay, electrophoretic mobility shift assay (EMSA), luciferase (LUC) reporter assay, and ChIP-PCR analysis revealed that MdWRKY31 binds to the promoter of MdLAC7 and positively regulates its activity to promote peel browning of bagged fruits (dark conditions). Interestingly, upon light exposure, the light-responsive transcription factor MdHY5 (ELONGATED HYPOCOTYL 5) bound to the promoter of MdWRKY31 and inhibited the gene's expression, thereby indirectly inhibiting the function of MdLAC7. Subsequent analysis showed that MdHY5 binds to the MdLAC7 promoter at the G-box1/2 site and directly inhibits its expression in vivo. Thus, the study revealed the MdLAC7-mediated mechanism regulating pre-harvest apple peel browning and demonstrated the role of light in inhibiting MdLAC7 activity and subsequently reducing peel browning. These results provide theoretical guidance for producing high-quality apple fruits.

A single QTL harboring multiple genetic variations leads to complicated phenotypic segregation in apple flesh firmness and crispness

Within a QTL, the genetic recombination and interactions among five and two functional variations at MdbHLH25 and MdWDR5A caused much complicated phenotype segregation in apple FFR and FCR. The storability of climacteric fruit like apple is a quantitative trait. We previously identified 62 quantitative trait loci (QTLs) associating flesh firmness retainability (FFR) and flesh crispness retainability (FCR), but only a few functional genetic variations were identified and validated. The genetic variation network controlling fruit storability is far to be understood and diagnostic markers are needed for molecular breeding. We previously identified overlapped QTLs F16.1/H16.2 for FFR and FCR using an F1 population derived from 'Zisai Pearl' × 'Red Fuji'. In this study, five and two single-nucleotide polymorphisms (SNPs) were identified on the candidate genes MdbHLH25 and MdWDR5A within the QTL region. The SNP1 A allele at MdbHLH25 promoter reduced the expression and SNP2 T allele and/or SNP4/5 GT alleles at the exons attenuated the function of MdbHLH25 by downregulating the expression of the target genes MdACS1, which in turn led to a reduction in ethylene production and maintenance of higher flesh crispness. The SNPs did not alter the protein-protein interaction between MdbHLH25 and MdWDR5A. The joint effect of SNP genotype combinations by the SNPs on MdbHLH25 (SNP1, SNP2, and SNP4) and MdWDR5A (SNPi and SNPii) led to a much broad spectrum of phenotypic segregation in FFR and FCR. Together, the dissection of these genetic variations contributes to understanding the complicated effects of a QTL and provides good potential for marker development in molecular breeding.

Characterization of the genetic and regulatory networks associated with sugar and acid metabolism in apples via an integrated strategy

Although sugars and acids have a substantial influence on the taste of apple fruits, the genetic and regulatory networks underlying their metabolism in fruit remain insufficiently determined. To fully decipher the genetic basis of the accumulation of sugars and acids in apple fruits, we adopted an integrated strategy that included time-course RNA-seq, QTL mapping, and whole-genome sequencing to examine two typical cultivars ('HanFu' and 'Huahong') characterized by distinctive flavors. Whole-genome sequencing revealed substantial genetic variation between the two cultivars, thereby providing an indication of the genetic basis of the distinct phenotypes. Constructed co-expression networks yielded information regarding the intra-relationships among the accumulation of different types of metabolites, and also revealed key regulatory nodes associated with the accumulation of sugars and acids, including the genes MdEF2, MdPILS5, and MdGUN8. Additionally, on the basis of QTL mapping using a high-density genetic map, we identified a series of QTLs and functional genes underlying vital traits, including sugar and acid contents. Collectively, our methodology and observations will provide an important reference for further studies focusing on the flavor of apples.

MbHY5-MbYSL7 mediates chlorophyll synthesis and iron transport under iron deficiency in Malus baccata

Iron (Fe) plays an important role in cellular respiration and catalytic reactions of metalloproteins in plants and animals. Plants maintain iron homeostasis through absorption, translocation, storage, and compartmentalization of iron via a cooperative regulative network. Here, we showed different physiological characteristics in the leaves and roots of Malus baccata under Fe sufficiency and Fe deficiency conditions and propose that MbHY5 (elongated hypocotyl 5), an important transcription factor for its function in photomorphogenesis, participated in Fe deficiency response in both the leaves and roots of M. baccata. The gene co-expression network showed that MbHY5 was involved in the regulation of chlorophyll synthesis and Fe transport pathway under Fe-limiting conditions. Specifically, we found that Fe deficiency induced the expression of MbYSL7 in root, which was positively regulated by MbHY5. Overexpressing or silencing MbYSL7 influenced the expression of MbHY5 in M. baccata.

MdWOX4-2 modulated MdLBD41 functioning in adventitious shoot of apple (Malus domestica)

Apple (Malus domestica Borkh.) is not only an important fruit crop distributed worldwide, but also a common model plant. However, the lack of efficient genetic transformation procedures for apples limits the in-depth studies of their gene functions. Although leaf-regenerated adventitious shoots (LRAS) are a prerequisite for successful genetic transformation of apple, little is known about the underlying molecular mechanism of LRAS. Here, we identified the WUSCHEL-related homeobox (WOX) transcription factor in apple, MdWOX4-2, which was a transcriptional activator. Gene expression as well as morphological and histological observations revealed that MdWOX4-2 is involved in the development of LRAS. Overexpression of MdWOX4-2 conferred higher regenerative capacity in transgenic tobacco (Nicotiana tabacum) as compared to the wild type (WT). The combined results of the yeast one-hybrid (Y1H), electrophoretic mobility shift assay (EMSA), dual luciferase assays, and transient transactivation assay, revealed that MdWOX4-2 directly bound to and activated the MdLBD41 promoter. Moreover, transgenic experiments further demonstrated that MdLBD41 could significantly enhance the formation of adventitious shoot in transgenic tobacco. Collectively, our findings demonstrate that MdWOX4-2 is important for regulating the LRAS development by activating MdLBD41.

Novel Role of AaMYBC1 in Regulating Actinidia arguta Vine Architecture by Elongating Internode Based on Multi-Omics Analysis of Transgenic Tobacco

The internode length affects the status of fruiting branches and shapes the vine architecture. MYB TFs (transcription factors) have been widely studied and reported to control many biological processes including secondary metabolism, abiotic stresses, growth and development, etc. However, the roles of MYB TFs in regulating internode length remain poorly understood. Here, we demonstrated that a secondary metabolism-related R2R3-MYB TF AaMYBC1 from Actinidia arguta was involved in the regulation of internode length by combined analysis of transcriptome and metabolome of transgenic tobacco plants. The metabolome analysis of OE (over-expressed tobacco) and WT (wild-typed tobacco) showed that there were a total of 1000 metabolites, 176 of which had significant differences. A key metabolite pme1651 annotated as indole 3-acetic acid belonged to phytohormone that was involved in internode length regulation. The RNA-seq analysis presented 446 differentially expressed genes (DEGs) between OE and WT, 14 of which were common DEGs in KEGG and GO enrichment. Through the combined analysis of metabolome and transcriptome in transgenic and wild-type tobacco, three key genes including two SAUR and a GH3 gene were possibly involved in internode elongation. Finally, a regulatory module was deduced to show the role of AaMYBC1 in internode elongation. Our results proposed a molecular mechanism of AaMYBC1 regulating internode length by mediated auxin signaling, implying the potential role in regulating the vine architecture.

MdWRKY126 modulates malate accumulation in apple fruit by regulating cytosolic malate dehydrogenase (MdMDH5)

The content of organic acids greatly influences the taste and storage life of fleshy fruit. Our current understanding of the molecular mechanism of organic acid accumulation in apple (Malus domestica) fruit focuses on the aluminum-activated malate transporter 9/Ma1 gene. In this study, we identified a candidate gene, MdWRKY126, for controlling fruit acidity independent of Ma1 using homozygous recessive mutants of Ma1, namely Belle de Boskoop "BSKP" and Aifeng "AF." Analyses of transgenic apple calli and flesh and tomato (Solanum lycopersicum) fruit demonstrated that MdWRKY126 was substantially associated with malate content. MdWRKY126 was directly bound to the promoter of the cytoplasmic NAD-dependent malate dehydrogenase MdMDH5 and promoted its expression, thereby enhancing the malate content of apple fruit. In MdWRKY126 overexpressing calli, the mRNA levels of malate-associated transporters and proton pump genes also significantly increased, which contributed to the transport of malate accumulated in the cytoplasm to the vacuole. These findings demonstrated that MdWRKY126 regulates malate anabolism in the cytoplasm and coordinates the transport between cytoplasm and vacuole to regulate malate accumulation. Our study provides useful information to improve our understanding of the complex mechanism regulating apple fruit acidity.

A bulked segregant analysis tool for out-crossing species (BSATOS) and QTL-based genomics-assisted prediction of complex traits in apple

INTRODUCTION

Genomic heterozygosity, self-incompatibility, and rich-in somatic mutations hinder the molecular breeding efficiency of outcrossing plants.

OBJECTIVES

We attempted to develop an efficient integrated strategy to identify quantitative trait loci (QTLs) and trait-associated genes, to develop gene markers, and to construct genomics-assisted prediction (GAP) modes.

METHODS

A novel protocol, bulked segregant analysis tool for out-crossing species (BSATOS), is presented here, which is characterized by taking full advantage of all segregation patterns (including AB × AB markers) and haplotype information. To verify the effectiveness of the protocol in dealing with the complex traits of outbreeding species, three apple cross populations with 9,654 individuals were adopted.

RESULTS

By using BSATOS, 90, 60, and 77 significant QTLs were identified successfully and candidate genes were predicted for apple fruit weight (FW), fruit ripening date (FRD), and fruit soluble solid content (SSC), respectively. The gene-based markers were developed and genotyped for 1,396 individuals in a training population, including 145 Malus accessions and 1,251 F1 plants of the three full-sib families. GAP models were trained using marker genotype effect estimates of the training population. The prediction accuracy was 0.7658, 0.6455, and 0.3758 for FW, FRD, and SSC, respectively.

CONCLUSION

The BSATOS and GAP models provided a convenient and efficient methodology for candidate gene mining and molecular breeding in out-crossing plant species. The BSATOS pipeline can be freely downloaded from: https://github.com/maypoleflyn/BSATOS.

Variation in the promoter of the sorbitol dehydrogenase gene MdSDH2 affects binding of the transcription factor MdABI3 and alters fructose content in apple fruit

Fructose (Fru) content is a key determinant of fruit sweetness and quality. An F1 hybrid population of the apple cultivars 'Honeycrisp' × 'Qinguan' was used to investigate the quantitative trait locus (QTL) regions and genes controlling Fru content in fruit. A stable QTL on linkage group (LG) 01 in 'Honeycrisp' was detected on the single nucleotide polymorphism (SNP) genetic linkage maps. In this region, a sorbitol dehydrogenase (SDH) gene, MdSDH2, was detected and showed promoter variations and differential expression patterns between 'Honeycrisp' and 'Qinguan' fruits as well as their hybrids. A SNP variant (A/G) in the MdSDH2 promoter region (SDH2p-491) affected the binding ability of the transcription factor MdABI3, which can affect the expression of MdSDH2. Promoter sequences with an A nucleotide at SDH2p-491 had stronger binding affinity for MdABI3 than those with a G. Among 27 domesticated apple cultivars and wild relatives, this SNP (A/G) was associated with Fru content. Our results indicate that MdSDH2 can alter Fru content as the major regulatory gene and that ABA signaling might be involved in Fru content accumulation in apple fruit.

Stable QTL for malate levels in ripe fruit and their transferability across Vitis species

Malate is a major contributor to the sourness of grape berries (Vitis spp.) and their products, such as wine. Excessive malate at maturity, commonly observed in wild Vitis grapes, is detrimental to grape and wine quality and complicates the introgression of valuable disease resistance and cold hardy genes through breeding. This study investigated an interspecific Vitis family that exhibited strong and stable variation in malate at ripeness for five years and tested the separate contribution of accumulation, degradation, and dilution to malate concentration in ripe fruit in the last year of study. Genotyping was performed using transferable rhAmpSeq haplotype markers, based on the Vitis collinear core genome. Three significant QTL for ripe fruit malate on chromosomes 1, 7, and 17, accounted for over two-fold and 6.9 g/L differences, and explained 40.6% of the phenotypic variation. QTL on chromosomes 7 and 17 were stable in all and in three out of five years, respectively. Variation in pre-veraison malate was the major contributor to variation in ripe fruit malate (39%), and based on two and five years of data, respectively, their associated QTL overlapped on chromosome 7, indicating a common genetic basis. However, use of transferable markers on a closely related Vitis family did not yield a common QTL across families. This suggests that diverse physiological mechanisms regulate the levels of this key metabolite in the Vitis genus, a conclusion supported by a review of over a dozen publications from the past decade, showing malate-associated genetic loci on all 19 chromosomes.

Multi-omics approaches identify a key gene, PpTST1, for organic acid accumulation in peach

Organic acid content in fruit is an important determinant of peach organoleptic quality, which undergoes considerable variations during development and maturation. However, its molecular mechanism remains largely unclear. In this study, an integrative approach of genome-wide association studies and comparative transcriptome analysis were applied to identify candidate genes involved in organic acid accumulation in peach. A key gene PpTST1, encoding tonoplast sugar transporter, was identified and the genotype of PpTST1 with a single-base transversion (G1584T) in the third exon which leads to a single amino acid substitution (Q528H) was associated with low level of organic acid content in peach. Overexpression of PpTST1His resulted in reduced organic acid content along with increased sugar content both in peach and tomato fruits, suggesting its dual function in sugar accumulation and organic acid content reduction. Two V-type proton ATPases interact with PpTST1 in yeast two-hybridization assay. In addition, the G1584T transversion appeared and gradually accumulated during domestication and improvement, which indicated that PpTST1 was under selection. The identification and characterization of PpTST1 would facilitate the improvement of peach fruit quality.

MrMYB44-Like Negatively Regulates Anthocyanin Biosynthesis and Causes Spring Leaf Color of Malus 'Radiant' to Fade From Red to Green

The "Spring-red-leaf" crabapple cultivar has young red leaves and mature green leaves. However, the mechanism of anthocyanin biosynthesis in crabapple leaves in spring remains unknown. In this study, Illumina RNA sequencing (RNA-Seq) was performed on Malus 'Radiant' leaf tissues in different stages of development. Twenty-two genes in the anthocyanin biosynthesis pathway and 44 MYB transcription factors (TFs) were significantly enriched among differentially expressed genes (DEGs). Three R2R3-MYB TFs in subgroup 22 of the MYB TF family, MrMYB44-like1, MrMYB44-like2, and MrMYB44-like3, were highly expressed in green leaves according to RNA-Seq and quantitative real-time quantitative PCR results. Their expression levels were negatively correlated with anthocyanin content. In transient assays, overexpression of MrMYB44-like1, MrMYB44-like2, or MrMYB44-like3 inhibited anthocyanin accumulation and reduced pigment in leaf disks of M. 'Radiant' and fruit peels of M. domestica 'Fuji.' When the conserved region of the three MrMYB44-likes was silenced, the anthocyanin biosynthesis pathway was activated and pigments increased in both tissues. Moreover, bimolecular fluorescence complementation assays showed MrMYB44-likes interacted with MrWRKY6 to form protein complexes that regulated anthocyanin biosynthesis.

PpMYB36 Encodes a MYB-Type Transcription Factor That Is Involved in Russet Skin Coloration in Pear (Pyrus pyrifolia)

Fruit color is one of the most important external qualities of pear (Pyrus pyrifolia) fruits. However, the mechanisms that control russet skin coloration in pear have not been well characterized. Here, we explored the molecular mechanisms that determine the russet skin trait in pear using the F₁ population derived from a cross between russet skin ('Niitaka') and non-russet skin ('Dangshansu') cultivars. Pigment measurements indicated that the lignin content in the skin of the russet pear fruits was greater than that in the non-russet pear skin. Genetic analysis revealed that the phenotype of the russet skin pear is associated with an allele of the PpRus gene. Using bulked segregant analysis combined with the genome sequencing (BSA-seq), we identified two simple sequence repeat (SSR) marker loci linked with the russet-colored skin trait in pear. Linkage analysis showed that the PpRus locus maps to the scaffold NW_008988489.1: 53297-211921 on chromosome 8 in the pear genome. In the mapped region, the expression level of LOC103929640 was significantly increased in the russet skin pear and showed a correlation with the increase of lignin content during the ripening period. Genotyping results demonstrated that LOC103929640 encoding the transcription factor MYB36 is the causal gene for the russet skin trait in pear. Particularly, a W-box insertion at the PpMYB36 promoter of russet skin pears is essential for PpMYB36-mediated regulation of lignin accumulation and russet coloration in pear. Overall, these results show that PpMYB36 is involved in the regulation of russet skin trait in pear.

Mechanisms and regulation of organic acid accumulation in plant vacuoles

In fleshy fruits, organic acids are the main source of fruit acidity and play an important role in regulating osmotic pressure, pH homeostasis, stress resistance, and fruit quality. The transport of organic acids from the cytosol to the vacuole and their storage are complex processes. A large number of transporters carry organic acids from the cytosol to the vacuole with the assistance of various proton pumps and enzymes. However, much remains to be explored regarding the vacuolar transport mechanism of organic acids as well as the substances involved and their association. In this review, recent advances in the vacuolar transport mechanism of organic acids in plants are summarized from the perspectives of transporters, channels, proton pumps, and upstream regulators to better understand the complex regulatory networks involved in fruit acid formation.

An apple (Malus domestica) AP2/ERF transcription factor modulates carotenoid accumulation

Color is an important trait for horticultural crops. Carotenoids are one of the main pigments for coloration and have important implications for photosynthesis in plants and benefits for human health. Here, we identified an APETALA2 (AP2)/ETHYLENE RESPONSE FACTOR (ERF) transcription factor named MdAP2-34 in apple (Malus domestica Borkh.). MdAP2-34 expression exhibited a close correlation with carotenoid content in 'Benin Shogun' and 'Yanfu 3' fruit flesh. MdAP2-34 promotes carotenoid accumulation in MdAP2-34-OVX transgenic apple calli and fruits by participating in the carotenoid biosynthesis pathway. The major carotenoid contents of phytoene and β-carotene were much higher in overexpressing MdAP2-34 transgenic calli and fruit skin, yet the predominant compound of lutein showed no obvious difference, indicating that MdAP2-34 regulates phytoene and β-carotene accumulation but not lutein. MdPSY2-1 (phytoene synthase 2) is a major gene in the carotenoid biosynthesis pathway in apple fruit, and the MdPSY2-1 gene is directly bound and transcriptionally activated by MdAP2-34. In addition, overexpressing MdPSY2-1 in apple calli mainly increases phytoene and total carotenoid contents. Our findings will advance and extend our understanding of the complex molecular mechanisms of carotenoid biosynthesis in apple, and this research is valuable for accelerating the apple breeding process.

Combined Profiling of Transcriptome and DNA Methylome Reveal Genes Involved in Accumulation of Soluble Sugars and Organic Acid in Apple Fruits

Organic acids and soluble sugars are the major determinants of fruit organoleptic quality. Additionally, DNA methylation has crucial regulatory effects on various processes. However, the epigenetic modifications in the regulation of organic acid and soluble sugar accumulation in apple fruits remain uncharacterized. In this study, DNA methylation and the transcriptome were compared between ‘Honeycrisp’ and ‘Qinguan’ mature fruits, which differ significantly regarding soluble sugar and organic acid contents. In both ‘Honeycrisp’ and ‘Qinguan’ mature fruits, the CG context had the highest level of DNA methylation, and then CHG and CHH contexts. The number and distribution of differentially methylated regions (DMRs) varied among genic regions and transposable elements. The DNA methylation levels in all three contexts in the DMRs were significantly higher in ‘Honeycrisp’ mature fruits than in ‘Qinguan’ mature fruits. A combined methylation and transcriptome analysis revealed a negative correlation between methylation levels and gene expression in DMRs in promoters and gene bodies in the CG and CHG contexts and in gene bodies in the CHH context. Two candidate genes (MdTSTa and MdMa11), which encode tonoplast-localized proteins, potentially associated with fruit soluble sugar contents and acidity were identified based on expression and DNA methylation levels. Overexpression of MdTSTa in tomato increased the fruit soluble sugar content. Moreover, transient expression of MdMa11 in tobacco leaves significantly decreased the pH value. Our results reflect the diversity in epigenetic modifications influencing gene expression and will facilitate further elucidating the complex mechanism underlying fruit soluble sugar and organic acid accumulation.

Unraveling a genetic roadmap for improved taste in the domesticated apple

Although taste is an important aspect of fruit quality, an understanding of its genetic control remains elusive in apple and other fruit crops. In this study, we conducted genomic sequence analysis of 497 Malus accessions and revealed erosion of genetic diversity caused by apple breeding and possible independent domestication events of dessert and cider apples. Signatures of selection for fruit acidity and size, but not for fruit sugar content, were detected during the processes of both domestication and improvement. Furthermore, we found that single mutations in major genes affecting fruit taste, including Ma1, MdTDT, and MdSOT2, dramatically decrease malate, citrate, and sorbitol accumulation, respectively, and correspond to important domestication events. Interestingly, Ma1 was identified to have pleiotropic effects on both organic acid content and sugar:acid ratio, suggesting that it plays a vital role in determining fruit taste. Fruit taste is unlikely to have been negatively affected by linkage drag associated with selection for larger fruit that resulted from the pyramiding of multiple genes with minor effects on fruit size. Collectively, our study provides new insights into the genetic basis of fruit quality and its evolutionary roadmap during apple domestication, pinpointing several candidate genes for genetic manipulation of fruit taste in apple.

Transcriptome and Resequencing Analyses Provide Insight into Differences in Organic Acid Accumulation in Two Pear Varieties

Fruit acidity is one of the main determinants of fruit flavor and a target trait in fruit breeding. However, the genomic mechanisms governing acidity variation among different pear varieties remain poorly understood. In this study, two pear varieties with contrasting organic acid levels, ‘Dangshansuli’ (low-acidity) and ‘Amute’ (high-acidity), were selected, and a combination of transcriptome and population genomics analyses were applied to characterize their patterns of gene expression and genetic variation. Based on RNA-seq data analysis, differentially expressed genes (DEGs) involved in organic acid metabolism and accumulation were identified. Weighted correlation network analysis (WGCNA) revealed that nine candidate TCA (tricarboxylic acid)-related DEGs and three acid transporter-related DEGs were located in three key modules. The regulatory networks of the above candidate genes were also predicted. By integrating pear resequencing data, two domestication-related genes were found to be upregulated in ‘Amute’, and this trend was further validated for other pear varieties with high levels of organic acid, suggesting distinct selective sweeps during pear dissemination and domestication. Collectively, this study provides insight into organic acid differences related to expression divergence and domestication in two pear varieties, pinpointing several candidate genes for the genetic manipulation of acidity in pears.

An HD-ZIP transcription factor, MxHB13, integrates auxin-regulated and juvenility-determined control of adventitious rooting in Malus xiaojinensis

Adventitious root (AR) formation is a critical factor in the vegetative propagation of forestry and horticultural plants. Competence for AR formation declines in many species during the miR156/SPL-mediated vegetative phase change. Auxin also plays a regulatory role in AR formation. In apple rootstock, both high miR156 expression and exogenous auxin application are prerequisites for AR formation. However, the mechanism by which the miR156/SPL module interacts with auxin in controlling AR formation is unclear. In this paper, leafy cuttings of juvenile (Mx-J) and adult (Mx-A) phase Malus xiaojinensis were used in an RNA-sequencing experiment. The results revealed that numerous genes involved in phytohormone signaling, carbohydrate metabolism, cell dedifferentiation, and reactivation were downregulated in Mx-A cuttings in response to indole butyric acid treatment. Among the differentially expressed genes, an HD-ZIP transcription factor gene, MxHB13, was found to be under negative regulation of MdSPL26 by directly binding to MxHB13 promoter. MxTIFY9 interacts with MxSPL26 and may play a role in co-repressing the expression of MxHB13. The expression of MxTIFY9 was induced by exogenous indole butyric acid. MxHB13 binds to the promoter of MxABCB19-2 and positively affects the expression. A model is proposed in which MxHB13 links juvenility-limited and auxin-limited AR recalcitrance mechanisms in Mx-A.

Genetic variation in the promoter of an R2R3-MYB transcription factor determines fruit malate content in apple (Malus domestica Borkh.)

Deciphering the mechanism of malate accumulation in apple (Malus domestica Borkh.) fruits can help to improve their flavor quality and enhance their benefits for human health. Here, we analyzed malate content as a quantitative trait that is determined mainly by genetic effects. In a previous study, we identified an R2R3-MYB transcription factor named MdMYB44 that was a candidate gene in qtl08.1 (quantitative trait locus mapped to chromosome 8) of fruit malate content. In the present study, we established that MdMYB44 negatively regulates fruit malate accumulation by repressing the promoter activity of the malate-associated genes Ma1 (Al-Activated Malate Transporter 9), Ma10 (P-type ATPase 10), MdVHA-A3 (V-type ATPase A3), and MdVHA-D2 (V-type ATPase D2). Two single-nucleotide polymorphisms (SNPs) in the MdMYB44 promoter, SNP A/G and SNP T/-, were experimentally shown to associate with fruit malate content. The TATA-box in the MdMYB44 promoter in the presence of SNP A enhances the basal activity of the MdMYB44 promoter. The binding of a basic-helix-loop-helix transcription factor MdbHLH49 to the MdMYB44 promoter was enhanced by the presence of SNP T, leading to increased MdMYB44 transcript levels and reduced malate accumulation. Furthermore, MdbHLH49 interacts with MdMYB44 and enhances MdMYB44 activity. The two SNPs could be used in combination to select for sour or non-sour apples, providing a valuable tool for the selection of fruit acidity by the apple breeding industry. This research is important for understanding the complex molecular mechanisms of fruit malate accumulation and accelerating the development of germplasm innovation in apple species and cultivars.

Role of MdERF3 and MdERF118 natural variations in apple flesh firmness/crispness retainability and development of QTL-based genomics-assisted prediction

Retention of flesh texture attributes during cold storage is critical for the long-term maintenance of fruit quality. The genetic variations determining flesh firmness and crispness retainability are not well understood. The objectives of this study are to identify gene markers based on quantitative trait loci (QTLs) and to develop genomics-assisted prediction (GAP) models for apple flesh firmness and crispness retainability. Phenotype data of 2664 hybrids derived from three Malus domestica cultivars and a M. asiatica cultivar were collected in 2016 and 2017. The phenotype segregated considerably with high broad-sense heritability of 83.85% and 83.64% for flesh firmness and crispness retainability, respectively. Fifty-six candidate genes were predicted from the 62 QTLs identified using bulked segregant analysis and RNA-seq. The genotype effects of the markers designed on each candidate gene were estimated. The genomics-predicted values were obtained using pyramiding marker genotype effects and overall mean phenotype values. Fivefold cross-validation revealed that the prediction accuracy was 0.5541 and 0.6018 for retainability of flesh firmness and crispness, respectively. An 8-bp deletion in the MdERF3 promoter disrupted MdDOF5.3 binding, reduced MdERF3 expression, relieved the inhibition on MdPGLR3, MdPME2, and MdACO4 expression, and ultimately decreased flesh firmness and crispness retainability. A 3-bp deletion in the MdERF118 promoter decreased its expression by disrupting the binding of MdRAVL1, which increased MdPGLR3 and MdACO4 expression and reduced flesh firmness and crispness retainability. These results provide insights regarding the genetic variation network regulating flesh firmness and crispness retainability, and the GAP models can assist in apple breeding.

Natural variations in a pectin acetylesterase gene, MdPAE10, contribute to prolonged apple fruit shelf life

Room-temperature shelf life is a key factor in fresh market apple (Malus domestica Borkh.) quality and commercial value. To investigate the genetic and molecular mechanism underlying apple shelf life, quantitative trait loci (QTL) were identified using bulked segregant analysis via sequencing (BSA-seq). Ethylene emission, flesh firmness, or crispness of apple fruit from 1,273 F₁ plants of M. asiatica Nakai 'Zisai Pearl' × M. domestica 'Golden Delicious' were phenotyped prior to and during 6 wk of room-temperature storage. Segregation of ethylene emission and the flesh firmness or crispness traits was detected in the population. Thirteen QTL, including three major ones, were identified on chromosome 03, 08, and 16. A candidate gene encoding pectin acetylesterase, MdPAE10, from the QTL Z16.1 negatively affected fruit shelf life. A 379-bp deletion in the coding sequence of MdPAE10 disrupted its function. A single nucleotide polymorphism (SNP) in the MdPAE10 promoter region reduced its transcription activity. These findings provided insight into the genetic control of fruit shelf life and can be potentially used in apple marker-assisted selection.

Multiomics-based dissection of citrus flavonoid metabolism using a Citrus reticulata × Poncirus trifoliata population

Deciphering the genetic basis of plant secondary metabolism will provide useful insights for genetic improvement and enhance our fundamental understanding of plant biological processes. Although citrus plants are among the most important fruit crops worldwide, the genetic basis of secondary metabolism in these plants is largely unknown. Here, we use a high-density linkage map to dissect large-scale flavonoid metabolic traits measured in different tissues (young leaf, old leaf, mature pericarp, and mature pulp) of an F₁ pseudo-testcross citrus population. We detected 80 flavonoids in this population and identified 138 quantitative trait loci (QTLs) for 57 flavonoids in these four tissues. Based on transcriptional profiling and functional annotation, twenty-one candidate genes were identified, and one gene encoding flavanone 3-hydroxylase (F3H) was functionally verified to result in naturally occurring variation in dihydrokaempferol content through genetic variations in its promoter and coding regions. The abundant data resources collected for diverse citrus germplasms here lay the foundation for complete characterization of the citrus flavonoid biosynthetic pathway and will thereby promote efficient utilization of metabolites in citrus quality improvement.

A candidate PpRPH gene of the D locus controlling fruit acidity in peach

A candidate gene, designate PpRPH, in the D locus was identified to control fruit acidity in peach. Fruit acidity has a strong impact on organoleptic quality of fruit. Peach fruit acidity is controlled by a large-effect D locus on chromosome 5. In this study, the D locus was mapped to a 509-kb interval, with two markers, 5dC720 and 5C1019, co-segregating with the non-acid/acid trait of peach fruit. Within this interval, a candidate gene encoding a putative small protein, designated PpRPH, showed a consistency between gene expression and fruit acidity, with up- and down-regulation in non-acidic and acidic fruits, respectively. Transient ectopic expression of PpRPH in tobacco leaves caused an increase of pH by approximately 40% compared to the control transformed with empty vector. Whereas, the concentrations of citrate and malate decreased significantly by 22% and 37%, respectively, with respect to the empty vector control. All these results suggest that PpRPH is a strong candidate gene of the D locus. These findings contribute to our overall understanding of the complex mechanism underlying fruit acidity in peach as well as that in other fruit crops.

Rapid Mining of Candidate Genes for Verticillium Wilt Resistance in Cotton Based on BSA-Seq Analysis

Cotton is a globally important cash crop. Verticillium wilt (VW) is commonly known as "cancer" of cotton and causes serious loss of yield and fiber quality in cotton production around the world. Here, we performed a BSA-seq analysis using an F_2:3 segregation population to identify the candidate loci involved in VW resistance. Two QTLs (qvw-D05-1 and qvw-D05-2) related to VW resistance in cotton were identified using two resistant/susceptible bulks from the F₂ segregation population constructed by crossing the resistant cultivar ZZM2 with the susceptible cultivar J11. A total of 30stop-lost SNPs and 42 stop-gained SNPs, which included 17 genes, were screened in the qvw-D05-2 region by SnpEff analysis. Further analysis of the transcriptome data and qRT-PCR revealed that the expression level of Ghir_D05G037630 (designated as GhDRP) varied significantly at certain time points after infection with V. dahliae. The virus-induced gene silencing of GhDRP resulted in higher susceptibility of the plants to V. dahliae than the control, suggesting that GhDRP is involved in the resistance to V. dahlia infection. This study provides a method for rapid mining of quantitative trait loci and screening of candidate genes, as well as enriches the genomic information and gene resources for the molecular breeding of disease resistance in cotton.

Quantitative trait loci-based genomics-assisted prediction for the degree of apple fruit cover color

Apple fruit cover color is an important appearance trait determining fruit quality, high degree of fruit cover color or completely red fruit skin is also the ultimate breeding goal. MdMYB1 has repeatedly been reported as a major gene controlling apple fruit cover color. There are also multiple minor-effect genes affecting degree of fruit cover color (DFC). This study was to identify genome-wide quantitative trait loci (QTLs) and to develop genomics-assisted prediction for apple DFC. The DFC phenotype data of 9,422 hybrids from five full-sib families of Malus asiatica 'Zisai Pearl', M. domestica 'Red Fuji', 'Golden Delicious', and 'Jonathan' were collected in 2014-2017. The phenotype varied considerably among hybrids with the same MdMYB1 genotype. Ten QTLs for DFC were identified using MapQTL and bulked segregant analysis via sequencing. From these QTLs, ten candidate genes were predicted, including MdMYB1 from a year-stable QTL on chromosome 9 of 'Zisai Pearl' and 'Red Fuji'. Then, kompetitive allele-specific polymerase chain reaction (KASP) markers were designed on these candidate genes and 821 randomly selected hybrids were genotyped. The genotype effects of the markers were estimated. MdMYB1-1 (represented by marker H162) exhibited a partial dominant allelic effect on MdMYB1-2 and showed non-allelic epistasis on markers H1245 and G6. Finally, a non-additive QTL-based genomics assisted prediction model was established for DFC. The Pearson's correlation coefficient between the genomic predicted value and the observed phenotype value was 0.5690. These results can be beneficial for apple genomics-assisted breeding and may provide insights for understanding the mechanism of fruit coloration.

Spatial variation of the soil bacterial community in major apple producing regions of China

AIMS

In China, apple production areas are largely from the coastal to inland areas and across varied climate zones. However, the relationship among soil micro-organisms, environmental factors and fruit quality has not been clearly confirmed in orchards. Here we attempted to identify the variation of soil bacteria in the main apple producing regions and reveal the relationship among climatic factor, soil properties, soil bacterial community and fruit quality.

METHODS AND RESULTS

Sixty soil samples were collected from six main apple producing areas in China. We examined the soil bacteria using bacterial 16S rRNA gene amplicon profiling. The results show that the soil bacterial diversity of apple orchards varied from the Bohai Bay Region to the Loess Plateau Region. Proteobacteria, Acidobacteria and Actinobacteria were the predominant taxa at the phylum level for all six areas. In the Bohai Bay and the Loess Plateau region, which are the two largest apple producing areas, Proteobacteria and Actinobacteria had the highest relative abundance, respectively. Furthermore, soil bacterial diversity showed positive correlation with the mean annual temperature (MAT), soil organic matter (SOM) and pH. Excluding a direct effect on the apple fruit quality, MAT exerted an indirect influence through soil SOM and pH to alter the relative abundance of dominant taxa and shift the bacterial diversity, which affects the apple fruit titratable acids and soluble solids.

CONCLUSIONS

Geographic variables underlie apple orchard soil bacterial communities vary according to spatial scale. Environmental factors exert an indirect effect on apple fruit quality via shaping soil bacterial community.

SIGNIFICANCE AND IMPACT OF THE STUDY

This study provides a list of bacteria associated with environmental factors and the ecological attributes of their interactions in apple orchards, which will improve our ability to promote soil bacterial functional capabilities in order to reduce the fertilizer input and enhance the fruit quality.

Identification of two QTLs associated with high fruit acidity in apple using pooled genome sequencing analysis

Acidity is a critical component determining apple fruit quality. Previous studies reported two major acidity quantitative trait loci (QTLs) on linkage groups (LGs) 16 (Ma) and 8 (Ma3), respectively, and their homozygous genotypes mama and ma3ma3 usually confer low titratable acidity (TA) (<3.0 mg ml-1) to apple fruit. However, apples of genotypes Ma- (MaMa and Mama) or Ma3- (Ma3Ma3 and Ma3ma3) frequently show an acidity range spanning both regular (TA 3.0-10.0 mg ml-1) and high (TA > 10 mg ml-1) acidity levels. To date, the genetic control for high-acidity apples remains essentially unknown. In order to map QTLs associated with high acidity, two genomic DNA pools, one for high acidity and the other for regular acidity, were created in an interspecific F1 population Royal Gala (Malus domestica) × PI 613988 (M. sieversii) of 191 fruit-bearing progenies. By Illumina paired-end sequencing of the high and regular acidity pools, 1,261,640 single-nucleotide variants (SNVs) commonly present in both pools were detected. Using allele frequency directional difference and density (AFDDD) mapping approach, one region on chromosome 4 and another on chromosome 6 were identified to be putatively associated with high acidity, and were named Ma6 and Ma4, respectively. Trait association analysis of DNA markers independently developed from the Ma6 and Ma4 regions confirmed the mapping of Ma6 and Ma4. In the background of MaMa, 20.6% of acidity variation could be explained by Ma6, 28.5% by Ma4, and 50.7% by the combination of both. The effects of Ma6 and Ma4 in the background of Mama were also significant, but lower. These findings provide important genetic insight into high acidity in apple.

BTB-BACK-TAZ domain protein MdBT2-mediated MdMYB73 ubiquitination negatively regulates malate accumulation and vacuolar acidification in apple

As an important primary metabolite, malate plays a key role in regulating osmotic pressure, pH homeostasis, stress tolerance, and fruit quality of apple. The R2R3-MYB transcription factor (TF) MdMYB73 was identified as a protein that plays a critical role in determining malate accumulation and vacuolar acidification by directly regulating the transcription of aluminum-activated malate transporter 9 (MdALMT9), vacuolar ATPase subunit A (MdVHA-A), and vacuolar pyrophosphatase 1 (MdVHP1) in apple. In addition, the bHLH TF MdCIbHLH1 interacts with MdMYB73 and enhances the transcriptional activity of MdMYB73. Our previous studies demonstrated that the BTB-BACK-TAZ domain protein MdBT2 can degrade MdCIbHLH1 to influence malate accumulation and vacuolar acidification. However, the potential upstream regulators of MdMYB73 are currently unknown. In this study, we found that MdBT2 directly interacts with and degrades MdMYB73 through the ubiquitin/26S proteasome pathway to regulate malate accumulation and vacuolar acidification. A series of functional assays with apple calli and fruit showed that MdBT2 controls malate accumulation and vacuolar acidification in an MdMYB73-dependent manner. Overall, our findings shed light on the mechanism by which the BTB-BACK-TAZ domain protein MdBT2 regulates malate accumulation and vacuolar acidification by targeting MdMYB73 and MdCIbHLH1 for ubiquitination in apple. This information may help guide traditional breeding programs and fruit tree molecular breeding, and lead to improvements in fruit quality and stress tolerance.

Genomics-assisted prediction of salt and alkali tolerances and functional marker development in apple rootstocks

Background

Saline, alkaline, and saline-alkaline stress severely affect plant growth and development. The tolerance of plants to these stressors has long been important breeding objectives, especially for woody perennials like apple. The aims of this study were to identify quantitative trait loci (QTLs) and to develop genomics-assisted prediction models for salt, alkali, and salt-alkali tolerance in apple rootstock.

Results

A total of 3258 hybrids derived from the apple rootstock cultivars ‘Baleng Crab’ (Malus robusta Rehd., tolerant) × ‘M9’ (M. pumila Mill., sensitive) were used to identify 17, 13, and two QTLs for injury indices of salt, alkali, and salt–alkali stress via bulked segregant analysis. The genotype effects of single nucleotide polymorphism (SNP) markers designed on candidate genes in each QTL interval were estimated. The genomic predicted value of an individual hybrid was calculated by adding the sum of all marker genotype effects to the mean phenotype value of the population. The prediction accuracy was 0.6569, 0.6695, and 0.5834 for injury indices of salt, alkali, and salt–alkali stress, respectively. SNP182G on MdRGLG3, which changes a leucine to an arginine at the vWFA-domain, conferred tolerance to salt, alkali, and salt-alkali stress. SNP761A on MdKCAB, affecting the Kv_beta domain that cooperated with the linked allelic variation SNP11, contributed to salt, alkali, and salt–alkali tolerance in apple rootstock.

Conclusions

The genomics-assisted prediction models can potentially be used in breeding saline, alkaline, and saline-alkaline tolerant apple rootstocks. The QTLs and the functional markers may provide insight for future studies into the genetic variation of plant abiotic stress tolerance.

Anthocyanin Biosynthesis and Methylation of the MdMYB10 Promoter Are Associated with the Red Blushed-Skin Mutant in the Red Striped-Skin "Changfu 2" Apple

The color of apple skin, particularly anthocyanin-based coloration, is a key factor determining market acceptance. The mechanisms of anthocyanin accumulation in apples with different skin color patterns (i.e., striped and blushed) were analyzed. In total, 14 anthocyanins and 5 procyanidins were simultaneously assayed in red blushed-skin mutants (CF-B1 and CF-B2) and red striped-skin parents (CF-S1 and CF-S2), and 13 significant differences were revealed. Anthocyanin accumulation was significantly higher in the red blushed-skin apples than it was in the parents. The transcript levels of anthocyanin biosynthesis genes and regulatory factors (MdMYB10, MdbHLH3, and MdWD40) were associated with different skin color patterns during the coloring period at 4, 6, and 8 days after the fruits were debagged. The methylation levels of the MdMYB10 promoter regions -1203 to -779 bp, -1667 to -1180 bp, and -2295 to -1929 bp were associated with different skin color patterns, and there was more methylation in red striped-skin apples. These results improve our understanding of anthocyanin accumulation and its underlying molecular mechanism in apples with different skin color patterns, thereby providing valuable information for apple breeding.

Ethylene-responsive factor 4 is associated with the desirable rind hardness trait conferring cracking resistance in fresh fruits of watermelon

Fruit rind plays a pivotal role in alleviating water loss and disease and particularly in cracking resistance as well as the transportability, storability and shelf-life quality of the fruit. High susceptibility to cracking due to low rind hardness is largely responsible for severe annual yield losses of fresh fruits such as watermelon in the field and during the postharvest process. However, the candidate gene controlling the rind hardness phenotype remains unclear to date. Herein, we report, for the first time, an ethylene-responsive transcription factor 4 (ClERF4) associated with variation in rind hardness via a combinatory genetic map with bulk segregant analysis (BSA). Strikingly, our fine-mapping approach revealed an InDel of 11 bp and a neighbouring SNP in the ClERF4 gene on chromosome 10, conferring cracking resistance in F₂ populations with variable rind hardness. Furthermore, the concomitant kompetitive/competitive allele-specific PCR (KASP) genotyping data sets of 104 germplasm accessions strongly supported candidate ClERF4 as a causative gene associated with fruit rind hardness variability. In conclusion, our results provide new insight into the underlying mechanism controlling rind hardness, a desirable trait in fresh fruit. Moreover, the findings will further enable the molecular improvement of fruit cracking resistance in watermelon via precisely targeting the causative gene relevant to rind hardness, ClERF4.

Perspectives of CRISPR/Cas-mediated cis-engineering in horticulture: unlocking the neglected potential for crop improvement

Directed breeding of horticultural crops is essential for increasing yield, nutritional content, and consumer-valued characteristics such as shape and color of the produce. However, limited genetic diversity restricts the amount of crop improvement that can be achieved through conventional breeding approaches. Natural genetic changes in cis-regulatory regions of genes play important roles in shaping phenotypic diversity by altering their expression. Utilization of CRISPR/Cas editing in crop species can accelerate crop improvement through the introduction of genetic variation in a targeted manner. The advent of CRISPR/Cas-mediated cis-regulatory region engineering (cis-engineering) provides a more refined method for modulating gene expression and creating phenotypic diversity to benefit crop improvement. Here, we focus on the current applications of CRISPR/Cas-mediated cis-engineering in horticultural crops. We describe strategies and limitations for its use in crop improvement, including de novo cis-regulatory element (CRE) discovery, precise genome editing, and transgene-free genome editing. In addition, we discuss the challenges and prospects regarding current technologies and achievements. CRISPR/Cas-mediated cis-engineering is a critical tool for generating horticultural crops that are better able to adapt to climate change and providing food for an increasing world population.

Genetic variations of acidity in grape berries are controlled by the interplay between organic acids and potassium

In a grapevine segregating population, genomic regions governing berry pH were identified, paving the way for breeding new grapevine varieties best adapted to a warming climate. As a consequence of global warming, grapevine berry acidity is expected to dramatically decrease. Adapting grapevine (Vitis vinifera L.) varieties to the climatic conditions of the future requires a better understanding of the genetic architecture of acidity-related traits. For this purpose, we studied during five growing seasons 120 individuals from a grapevine biparental cross. Each offspring was genotyped by simple sequence repeats markers and by hybridization on a 20-K Grapevine Illumina^® SNP chip. Quantitative trait loci (QTLs) for pH colocalized with QTLs for the ratio between potassium and tartaric acid concentrations, on chromosomes 10, 11 and 13. Strong QTLs for malic acid concentration or for the malic acid-to-tartaric acid ratio, on chromosomes 6 and 8, were not associated with variations of pH but can be useful for controlling pH stability under high temperatures. Our study highlights the interdependency between acidity parameters and consequently the constraints and degrees of freedom for designing grapevine genotypes better adapted to the expected warmer climatic conditions. In particular, it is possible to create grapevine genotypes with a high berry acidity as the result of both high tartaric acid concentrations and low K⁺ accumulation capacities.

Parallel Bud Mutation Sequencing Reveals that Fruit Sugar and Acid Metabolism Potentially Influence Stress in Malus

Apple sugar and acid are the most important traits of apple fruit. Bud sport cultivars can provide abundant research materials for functional gene studies in apple. In this study, using bud sport materials with a rather different sugar and acid flavor, i.e., "Jonathan" and "Sweet Jonathan", we profiled the whole genome variations and transcriptional regulatory network during fruit developmental stages using whole genome sequencing and RNA-sequencing. Variation analysis identified 4,198,955 SNPs, 319,494 InDels, and 32,434 SVs between the two cultivars. In total, 4313 differentially expressed genes among all of the d 44,399 genes expressed were identified between the two cultivars during fruit development, and functional analysis revealed stress response and signal transduction related genes were enriched. Using 24,047 genes with a more variable expression value, we constructed 28 co-expression modules by weighted correlation network analysis. Deciphering of 14 co-expression modules associated with sugar or acid accumulation during fruit development revealed the hub genes associated with sugar and acid metabolism, e.g., MdDSP4, MdINVE, and MdSTP7. Furthermore, exploration of the intra network of the co-expression module indicated the close relationship between sugar and acid metabolism or sugar and stress. Motif-based sequence analysis of the 17 differentially expressed ATP-binding cassette transporter genes and Yeast one-hybrid assay identified and confirmed a transcription factor, MdBPC6, regulating the ATP-binding cassette (ABC) transporter genes and potentially participating in the apple fruit development or stress response. Collectively, all of the results demonstrated the use of parallel bud mutation sequencing and identified hub genes, and inferred regulatory relationships providing new information about apple fruit sugar and acid accumulation or stress response.

Fruit Salad in the Lab: Comparing Botanical Species to Help Deciphering Fruit Primary Metabolism

Although fleshy fruit species are economically important worldwide and crucial for human nutrition, the regulation of their fruit metabolism remains to be described finely. Fruit species differ in the origin of the tissue constituting the flesh, duration of fruit development, coordination of ripening changes (climacteric vs. non-climacteric type) and biochemical composition at ripeness is linked to sweetness and acidity. The main constituents of mature fruit result from different strategies of carbon transport and metabolism. Thus, the timing and nature of phloem loading and unloading can largely differ from one species to another. Furthermore, accumulations and transformations of major soluble sugars, organic acids, amino acids, starch and cell walls are very variable among fruit species. Comparing fruit species therefore appears as a valuable way to get a better understanding of metabolism. On the one hand, the comparison of results of studies about species of different botanical families allows pointing the drivers of sugar or organic acid accumulation but this kind of comparison is often hampered by heterogeneous analysis approaches applied in each study and incomplete dataset. On the other hand, cross-species studies remain rare but have brought new insights into key aspects of primary metabolism regulation. In addition, new tools for multi-species comparisons are currently emerging, including meta-analyses or re-use of shared metabolic or genomic data, and comparative metabolic flux or process-based modeling. All these approaches contribute to the identification of the metabolic factors that influence fruit growth and quality, in order to adjust their levels with breeding or cultural practices, with respect to improving fruit traits.

Marker-trait associations and genomic predictions of interspecific pear (Pyrus) fruit characteristics

Interspecific pear (Pyrus spp.) hybrid populations are often used to develop novel cultivars. Pear cultivar breeding is a lengthy process because of long juvenility and the subsequent time required for reliable fruit phenotyping. Molecular techniques such as genome-wide association (GWA) and genomic selection (GS) provide an opportunity to fast-forward the development of high-value cultivars. We evaluated the genetic architecture of 10 pear fruit phenotypes (including sensory traits) and the potential of GS using genotyping-by-sequencing of 550 hybrid seedlings from nine interrelated full-sib families. Results from GWA suggested a complex polygenic nature of all 10 traits as the maximum variance explained by each marker was less than 4% of the phenotypic variance. The effect-size of SNPs for each trait suggested many genes of small effect and few of moderate effect. Some genomic regions associated with pear sensory traits were similar to those reported for apple - possibly a result of high synteny between the apple and pear genomes. The average (across nine families) GS accuracy varied from 0.32 (for crispness) to 0.62 (for sweetness), with an across-trait average of 0.42. Further efforts are needed to develop larger genotype-phenotype datasets in order to predict fruit phenotypes of untested seedlings with sufficient efficiency.

Apple whole genome sequences: recent advances and new prospects

In 2010, a major scientific milestone was achieved for tree fruit crops: publication of the first draft whole genome sequence (WGS) for apple (Malus domestica). This WGS, v1.0, was valuable as the initial reference for sequence information, fine mapping, gene discovery, variant discovery, and tool development. A new, high quality apple WGS, GDDH13 v1.1, was released in 2017 and now serves as the reference genome for apple. Over the past decade, these apple WGSs have had an enormous impact on our understanding of apple biological functioning, trait physiology and inheritance, leading to practical applications for improving this highly valued crop. Causal gene identities for phenotypes of fundamental and practical interest can today be discovered much more rapidly. Genome-wide polymorphisms at high genetic resolution are screened efficiently over hundreds to thousands of individuals with new insights into genetic relationships and pedigrees. High-density genetic maps are constructed efficiently and quantitative trait loci for valuable traits are readily associated with positional candidate genes and/or converted into diagnostic tests for breeders. We understand the species, geographical, and genomic origins of domesticated apple more precisely, as well as its relationship to wild relatives. The WGS has turbo-charged application of these classical research steps to crop improvement and drives innovative methods to achieve more durable, environmentally sound, productive, and consumer-desirable apple production. This review includes examples of basic and practical breakthroughs and challenges in using the apple WGSs. Recommendations for "what's next" focus on necessary upgrades to the genome sequence data pool, as well as for use of the data, to reach new frontiers in genomics-based scientific understanding of apple.

Natural variation in cytokinin maintenance improves salt tolerance in apple rootstocks

Plants experiencing salt-induced stress often reduce cytokinin levels during the early phases of stress-response. Interestingly, we found that the cytokinin content in the apple rootstock "robusta" was maintained at a high level under salt stress. Through screening genes involved in cytokinin biosynthesis and catabolism, we found that the high expression levels of IPT5b in robusta roots were involved in maintaining the high cytokinin content. We identified a 42 bp deletion in the promoter region of IPT5b, which elevated IPT5b expression levels, and this deletion was linked to salt tolerance in robusta×M.9 segregating population. The 42 bp deletion resulted in the deletion of a Proline Response Element (ProRE), and our results suggest that ProRE negatively regulates IPT5b expression in response to proline. Under salt stress, the robusta cultivar maintains high cytokinin levels as IPT5b expression cannot be inhibited by proline due to the deletion of ProRE, leading to improve salt tolerance.