Chinese bayberry (Myrica rubra or Morella rubra; 2n = 16) produces fruit with a distinctive flavor, high nutritional, and economic value. However, previous versions of the bayberry genome lack sequence continuity. Moreover, to date, no large-scale germplasm resource association analysis has examined the allelic and genetic variations determining fruit quality traits. Therefore, in this study, we assembled a telomere-to-telomere (T2T) gap-free reference genome for the cultivar ‘Zaojia’ using PacBio HiFi long reads. The resulting 292.60 Mb T2T genome, revealed 8 centromeric regions, 15 telomeres, and 28 345 genes. This represents a substantial improvement in the genome continuity and integrity of Chinese bayberry. Subsequently, we re-sequenced 173 accessions, identifying 6 649 674 single nucleotide polymorphisms (SNPs). Further, the phenotypic analyses of 29 fruit quality-related traits enabled a genome-wide association study (GWAS), which identified 1937 SNPs and 1039 genes significantly associated with 28 traits. An SNP cluster pertinent to fruit color was identified on Chr6: 3407532 to 5 153 151 bp region, harboring two MYB genes (MrChr6G07650 and MrChr6G07660), exhibiting differential expression in extreme phenotype transcriptomes, linked to anthocyanin synthesis. An adjacent, closely linked gene, MrChr6G07670 (MLP-like protein), harbored an exonic missense variant and was shown to increase anthocyanin production in tobacco leaves tenfold. This SNP cluster, potentially a quantitative trait locus (QTL), collectively regulates bayberry fruit color. In conclusion, our study presented a complete reference genome, uncovered a suite of allelic variations related to fruit-quality traits, and identified functional genes that could be harnessed to enhance fruit quality and breeding efficiency of bayberries.

Twig blight disease is the primary disease that affects the production of Myrica rubra in China. It was reported that exogenous brassinolide (BL) can improve disease resistance in plants. Here, we examined the effects of exogenous BL on disease resistance, chlorophyll contents, antioxidant enzyme activity, ROS accumulation, and key gene expression of M. rubra to analyze the mechanism of BR-induced resistance of twig blight disease in M. rubra. The results demonstrated that 2.0 mg·L⁻¹ of BL could significantly lessen the severity of twig blight disease in M. rubra. Exogenous BL increased the contents of chlorophyll a, chlorophyll b, carotenoids, and total chlorophyll. Moreover, exogenous BL also significantly enhanced the activity of antioxidant enzymes such as superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT), and decreased malondialdehyde (MDA) content and reactive oxygen species (ROS) accumulation in leaves, such as H₂O₂ and O₂^·−. Additionally, exogenous BL dramatically up-regulated the expression of pathogenesis-related (PR) genes such as MrPR1, MrPR2, and MrPR10, as well as important genes such as MrBAK1, MrBRI1, and MrBZR1 involved in brassinosteroid (BR) signaling pathway. The transcriptome analysis revealed that a total of 730 common differentially expressed genes (DEGs) under BL treatment were found, and these DEGs were primarily enriched in four Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Based on these findings, nine important candidate genes related to the resistance of twig blight disease under BL treatment were further identified. In this study, we elucidated the effects of exogenous BL on enhancing the resistance of M. rubra to twig blight disease and preliminary analyzed the potential mechanism of resistance induction, which will provide a crucial foundation for the management and prevention of twig blight disease in M. rubra.

The ripening process of Chinese bayberries (Myrica rubra) is intricate, involving a multitude of molecular interactions. Here, we integrated transcriptomic and metabolomic analysis across three developmental stages of the Myrica rubra (M. rubra) to elucidate these processes. A differential gene expression analysis categorized the genes into four distinct groups based on their expression patterns. Gene ontology and pathway analyses highlighted processes such as cellular and metabolic processes, including protein and sucrose metabolism. A metabolomic analysis revealed significant variations in metabolite profiles, underscoring the dynamic interplay between genes and metabolites during ripening. Flavonoid biosynthesis and starch and sucrose metabolism were identified as key pathways, with specific genes and metabolites playing crucial roles. Our findings provide insights into the molecular mechanisms governing fruit ripening in M. rubra and offer potential targets for breeding strategies aimed at enhancing fruit quality.

Fruit development and softening play pivotal roles in determining fruit quality and post-harvest shelf life in Chinese bayberry (Myrica rubra). However, the specific role of beta (β)-galactosidase, particularly β-galactosidase of M. rubra (MrBGAL), in facilitating fruit softening remains unclear. In this study, we aimed to address this gap by investigating the involvement of MrBGALs genes in fruit softening. We identified all 15 MrBGALs and conducted a comprehensive analysis, including phylogenetic relationships, gene structure, protein motifs, co-linearity, and expression patterns. Using phylogenetic analysis, we classified all MrBGALs into five distinct groups. Additionally, cis-element prediction and comparative genome analysis provided insightful clues about the functionality of MrBGALs. Transcriptome data revealed unique expression patterns of MrBGALs throughout various fruit development stages. These findings introduce valuable candidate genes that can contribute to unraveling the functions and molecular mechanisms governing fruit development and softening in Chinese bayberry.