The OsMYC2-JA feedback loop regulates diurnal flower-opening time via cell wall loosening in rice.

Diurnal flower-opening time (DFOT), the time of spikelet opening during the day, is an important trait for hybrid rice (Oryza sativa L.) seed production. Hybrids between indica and japonica rice varieties have strong heterosis, but the parental lines usually have different, nonoverlapping DFOTs. This reduces the success of hybrid seed production in crosses between indica and japonica subspecies, thus hindering the utilization of indica and japonica inter-subspecies heterosis. However, little is known about the molecular mechanisms regulating DFOT in rice. Here, we obtained japonica rice lines with a DFOT 1.5 h earlier than the wild type by overexpressing OsMYC2, a gene encoding a key transcription factor in the jasmonate (JA) signaling pathway. OsMYC2 is activated by JA signaling and directly regulates the transcription of genes related to JA biosynthesis and cell wall metabolism. Overexpressing OsMYC2 led to significantly increased JA contents and decreased cellulose and hemicellulose contents in lodicule cells, as well as the softening of lodicule cell walls. This may facilitate the swelling of lodicules, resulting in early diurnal flower-opening. These results suggest that the OsMYC2-JA feedback loop regulates DFOT in rice via cell wall remodeling. These findings shed light on the understanding of regulatory mechanism of the DFOT of plants, which should promote the development of indica and japonica varieties suitable for hybrid rice breeding.

GWAS and Transcriptomic Analysis Identify OsRING315 as a New Candidate Gene Controlling Amylose Content and Gel Consistency in Rice.

Cooking quality is the main factor determining the market value of rice. Although several major genes and a certain number of QTLs controlling cooking quality have been identified, the genetic complexity and environmental susceptibility limit the further improvement for cooking quality by molecular breeding. This research conducted a genome-wide association study to elucidate the QTLs related to cooking quality including amylose content (AC), gel consistency (GC) and alkali spreading value (ASV) by using 450 rice accessions consisting of 300 indica and 150 japonica accessions in two distinct environments. A total of 54 QTLs were identified, including 25 QTLs for AC, 12 QTLs for GC and 17 QTLs for ASV. Among them, 10 QTLs were consistently observed by the same population in both environments. Six QTLs were co-localized with the reported QTLs or cloned genes. The Wx gene for AC and GC, and the ALK gene for ASV were identified in every population across the two environments. The qAC9-2 for AC and the qGC9-2 for GC were defined to the same interval. The OsRING315 gene, encoding an E3 ubiquitin ligase, was considered as the candidate gene for both qAC9-2 and qGC9-2. The higher expression of OsRING315 corresponded to the lower AC and higher GC. Three haplotypes of OsRING315 were identified. The Hap 1 mainly existed in the japonica accessions and had lower AC. The Hap 2 and Hap 3 were predominantly present in the indica accessions, associated with higher AC. Meanwhile, the GC of accessions harboring Hap 1 was higher than that of accessions harboring Hap 3. In addition, the distribution of the three haplotypes in several rice-growing regions was unbalanced. The three traits of cooking quality are controlled by both major and minor genes and susceptible to environmental factors. The expression level of OsRING315 is related to both AC and GC, and this gene can be a promising target in quality improvement by using the gene editing method. Moreover, the haplotypes of OsRING315 differentiate between indica and japonica, and reveal the differences in GC and AC between indica and japonica rice.

Endoplasmic reticulum stress response modulator OsbZIP39 regulates cadmium accumulation via activating the expression of defensin-like gene OsCAL2 in rice.

Accumulation of cadmium (Cd) in rice is not only harmful to the growth of plants but also poses a threat to human health. Exposure to Cd triggers unfolded protein response (UPR) within cells, a process that is still not completely understood. The study demonstrated that the lack of OsbZIP39, an essential endoplasmic reticulum (ER)-resident regulator of the UPR, resulted in decreased Cd intake and reduced Cd levels in the roots, stems, and grains of rice. Upon exposure to Cd stress, GFP-OsbZIP39 translocated from ER to nucleus, initiating UPR. Further investigation revealed that Cd treatment caused changes in sphingolipid levels in the membrane, influencing the localization and activation of OsbZIP39. Yeast one-hybrid and dual-LUC assays were conducted to validate the interaction between activated OsbZIP39 and the promoter of the defensin-like gene OsCAL2, resulting in an increase in its expression. Different variations were identified in the coding region of OsbZIP39, which may explain the varying levels of Cd accumulation observed in the indica and japonica subspecies. Under Cd treatment, OsbZIP39ind exhibited a more significant enhancement in the transcription of OsCAL2 compared to OsbZIP39jap. Our data suggest that OsbZIP39 positively regulates Cd uptake in rice, offering an encouraging objective for the cultivation of low-Cd rice.

Flooding Tolerance of Rice: Regulatory Pathways and Adaptive Mechanisms.

Rice is a major food crop for more than half of the world's population, while its production is seriously threatened by flooding, a common environmental stress worldwide. Flooding leads to oxygen deficiency, which is a major problem for submerged plants. Over the past three decades, significant progress has been made in understanding rice adaptation and molecular regulatory mechanisms in response to flooding. At the seed germination and seedling establishment stages, the CIPK15-SnRK1A-MYBS1 signaling cascade plays a central role in determining rice submergence tolerance. However, from seedlings to mature plants for harvesting, SUB1A- and SK1/SK2-regulated pathways represent two principal and opposite regulatory mechanisms in rice. In addition, phytohormones, especially gibberellins, induce adaptive responses to flooding throughout the rice growth period. This review summarizes the significant adaptive traits observed in flooded rice varieties and updates the molecular genetics and mechanisms of submergence tolerance in rice.

Pangenome-Wide Association Study and Transcriptome Analysis Reveal a Novel QTL and Candidate Genes Controlling both Panicle and Leaf Blast Resistance in Rice.

Cultivating rice varieties with robust blast resistance is the most effective and economical way to manage the rice blast disease. However, rice blast disease comprises leaf and panicle blast, which are different in terms of resistance mechanisms. While many blast resistant rice cultivars were bred using genes conferring resistance to only leaf or panicle blast, mining durable and effective quantitative trait loci (QTLs) for both panicle and leaf blast resistance is of paramount importance. In this study, we conducted a pangenome-wide association study (panGWAS) on 9 blast resistance related phenotypes using 414 international diverse rice accessions from an international rice panel. This approach led to the identification of 74 QTLs associated with rice blast resistance. One notable locus, qPBR1, validated in a F4:5 population and fine-mapped in a Heterogeneous Inbred Family (HIF), exhibited broad-spectrum, major and durable blast resistance throughout the growth period. Furthermore, we performed transcriptomic analysis of 3 resistant and 3 sensitive accessions at different time points after infection, revealing 3,311 differentially expressed genes (DEGs) potentially involved in blast resistance. Integration of the above results identified 6 candidate genes within the qPBR1 locus, with no significant negative effect on yield. The results of this study provide valuable germplasm resources, QTLs, blast response genes and candidate functional genes for developing rice varieties with enduring and broad-spectrum blast resistance. The qPBR1, in particular, holds significant potential for breeding new rice varieties with comprehensive and durable resistance throughout their growth period.

Technological Development and Advances for Constructing and Analyzing Plant Pangenomes.

A pangenome captures the genomic diversity for a species, derived from a collection of genetic sequences of diverse populations. Advances in sequencing technologies have given rise to three primary methods for pangenome construction and analysis: de novo assembly and comparison, reference genome-based iterative assembly, and graph-based pangenome construction. Each method presents advantages and challenges in processing varying amounts and structures of DNA sequencing data. With the emergence of high-quality genome assemblies and advanced bioinformatic tools, the graph-based pangenome is emerging as an advanced reference for exploring the biological and functional implications of genetic variations.

Haplotype-resolved genomes of octoploid species in Phyllanthaceae family reveal a critical role for polyploidization and hybridization in speciation.

The Phyllanthaceae family comprises a diverse range of plants with medicinal, edible, and ornamental value, extensively cultivated worldwide. Polyploid species commonly occur in Phyllanthaceae. Due to the rather complex genomes and evolutionary histories, their speciation process has been still lacking in research. In this study, we generated chromosome-scale haplotype-resolved genomes of two octoploid species (Phyllanthus emblica and Sauropus spatulifolius) in Phyllanthaceae family. Combined with our previously reported one tetraploid (Sauropus androgynus) and one diploid species (Phyllanthus cochinchinensis) from the same family, we explored their speciation history. The three polyploid species were all identified as allopolyploids with subgenome A/B. Each of their two distinct subgenome groups from various species was uncovered to independently share a common diploid ancestor (Ancestor-AA and Ancestor-BB). Via different evolutionary routes, comprising various scenarios of bifurcating divergence, allopolyploidization (hybrid polyploidization), and autopolyploidization, they finally evolved to the current tetraploid S. androgynus, and octoploid S. spatulifolius and P. emblica, respectively. We further discuss the variations in copy number of alleles and the potential impacts within the two octoploids. In addition, we also investigated the fluctuation of metabolites with medical values and identified the key factor in its biosynthesis process in octoploids species. Our study reconstructed the evolutionary history of these Phyllanthaceae species, highlighting the critical roles of polyploidization and hybridization in their speciation processes. The high-quality genomes of the two octoploid species provide valuable genomic resources for further research of evolution and functional genomics.

Controlling diurnal flower-opening time by manipulating the jasmonate pathway accelerates development of indica-japonica hybrid rice breeding.

Inter-subspecific indica-japonica hybrid rice (Oryza sativa) has the potential for increased yields over traditional indica intra-subspecies hybrid rice, but limited yield of F1 hybrid seed production (FHSP) hinders the development of indica-japonica hybrid rice breeding. Diurnal flower-opening time (DFOT) divergence between indica and japonica rice has been a major contributing factor to this issue, but few DFOT genes have been cloned. Here, we found that manipulating the expression of jasmonate (JA) pathway genes can effectively modulate DFOT to improve the yield of FHSP in rice. Treating japonica cultivar Zhonghua 11 (ZH11) with methyl jasmonate (MeJA) substantially advanced DFOT. Furthermore, overexpressing the JA biosynthesis gene OPDA REDUCTASE 7 (OsOPR7) and knocking out the JA inactivation gene CHILLING TOLERANCE 1 (OsHAN1) in ZH11 advanced DFOT by 1- and 2-h respectively; and knockout of the JA signal suppressor genes JASMONATE ZIM-DOMAIN PROTEIN 7 (OsJAZ7) and OsJAZ9 resulted in 50-min and 1.5-h earlier DFOT respectively. The yields of FHSP using japonica male-sterile lines GAZS with manipulated JA pathway genes were significantly higher than that of GAZS wildtype. Transcriptome analysis, cytological observations, measurements of elastic modulus and determination of cell wall components indicated that the JA pathway could affect the loosening of the lodicule cell walls by regulating their composition through controlling sugar metabolism, which in turn influences DFOT. This research has vital implications for breeding japonica rice cultivars with early DFOT to facilitate indica-japonica hybrid rice breeding.

The Function of SD1 on Shoot Length and its Pyramiding Effect on Shoot Length and Plant Height in Rice (Oryza sativa L.).

Strong seedling vigor is imperative to achieve stable seedling establishment and enhance the competitiveness against weeds in rice direct seeding. Shoot length (SL) is one of the important traits associated with seedling vigor in rice, but few genes for SL have been cloned so far. In the previous study, we identified two tightly linked and stably expressed QTLs for SL, qSL-1f and qSL-1d by genome-wide association study, and cloned the causal gene (LOC_Os01g68500) underlying qSL-1f. In the present study, we identify LOC_Os01g66100 (i.e. the semidwarf gene SD1), a well-known gene controlling plant height (PH) at the adult-plant stage, as the causal gene underlying qSL-1d through gene-based haplotype analysis and knockout transgenic verification. By measuring the phenotypes (SL and PH) of various haplotypes of the two genes and their knockout lines, we found SD1 and LOC_ Os01g68500 controlled both SL and PH, and worked in the same direction, which provided the directly genetic evidence for a positive correlation between SL and PH combined with the analysis of SL and PH in the diverse natural population. Moreover, the knockout transgenic experiments suggested that SD1 had a greater effect on PH compared with LOC_ Os01g68500, but no significant difference in the effect on SL. Further investigation of the pyramiding effects of SD1 and LOC_Os01g68500 based on their haplotype combinations suggested that SD1 may play a dominant role in controlling SL and PH when the two genes coexist. In this study, the effect of SD1 on SL at the seedling stage is validated. In total, two causal genes, SD1 and LOC_ Os01g68500, for SL are cloned in our studies, which controlled both SL and PH, and the suitable haplotypes of SD1 and LOC_ Os01g68500 are beneficial to achieve the desired SL and PH in different rice breeding objectives. These results provide a new clue to develop rice varieties for direct seeding and provide new genetic resources for molecular breeding of rice with suitable PH and strong seedling vigor.

Unravelling inversions: Technological advances, challenges, and potential impact on crop breeding.

Inversions, a type of chromosomal structural variation, significantly influence plant adaptation and gene functions by impacting gene expression and recombination rates. However, compared with other structural variations, their roles in functional biology and crop improvement remain largely unexplored. In this review, we highlight technological and methodological advancements that have allowed a comprehensive understanding of inversion variants through the pangenome framework and machine learning algorithms. Genome editing is an efficient method for inducing or reversing inversion mutations in plants, providing an effective mechanism to modify local recombination rates. Given the potential of inversions in crop breeding, we anticipate increasing attention on inversions from the scientific community in future research and breeding applications.

Sustenance Trial to Analyze the Effects of Black Soldier Fly Larvae Meal on the Reproductive Efficiency of Sows and the Hematological Properties of Suckling and Weaning Piglets.

The escalating demand for meat, driven by global population growth, necessitates sustainable solutions for animal feed production. This study investigated the effects of substituting conventional protein resources in sow and piglet dietary regimens with black soldier fly (BSF; Hermetia illucens) meal on reproductive efficiency, blood profile, piglet growth, and intestinal tissue morphology. The results indicate that substituting animal-derived and soy proteins with BSF meal does not compromise sow reproductive performance. Although no notable disparities were observed in piglet growth, the feed conversion ratio from the 28- to 35-day age marks were lower in the BSF-fed groups. This suggests that the animal protein-BSF substitution rate may require optimization, potentially involving chitin removal from BSF meal to enhance digestibility. Minor variations in the hematological composition and properties in piglets, with elevated high-density lipoprotein cholesterol levels in the high BSF group at the 28-day mark, were potentially attributable to the unique fatty acid composition of BSF meal. Moreover, this study potentiates future exploration into the efficacy of complete animal protein substitution with BSF meals on piglet nutrition and physiology, particularly in fattening pigs. The practical implementation of BSF meals in animal feed production holds promise for enhancing the sustainability of the swine industry.

Genome-wide association mapping combined with gene-based haplotype analysis identify a novel gene for shoot length in rice (Oryza sativa L.).

KEY MESSAGE: Genome-wide association mapping revealed a novel QTL for shoot length across multiple environments. Its causal gene, LOC_Os01g68500, was identified firstly through gene-based haplotype analysis, gene expression and knockout transgenic verification. Strong seedling vigor is an important breeding target for rice varieties used in direct seeding. Shoot length (SL) is one of the important traits associated with seedling vigor characterized by rapid growth of seedling, which enhance seedling emergence. Therefore, mining genes for SL and conducting molecular breeding help to develop varieties for direct seeding. However, few QTLs for SL have been fine mapped or cloned so far. In this study, a genome-wide association study of SL was performed in a diverse rice collection consisting of 391 accessions in two years, using phenotypes generated by different cultivation methods according to the production practice, and a total of twenty-four QTLs for SL were identified. Among them, the novel QTL qSL-1f on chromosome 1 could be stably detected across all three cultivation methods in the whole population and indica subpopulation. Through gene-based haplotype analysis of the annotated genes within the putative region of qSL-1f, and validated by gene expression and knockout transgenic experiments, LOC_Os01g68500 (i.e., Os01g0913100 in RAP-DB) was identified as the causal gene for SL, which has a single-base variation (C-to-A transversion) in its CDS region, resulting in the significant difference in SL of rice. LOC_Os01g68500 encodes a DUF538 (Domain of unknown function) containing protein, and the function of DUF538 protein gene on rice seedling growth is firstly reported in this study. These results provide a new clue for exploring the molecular mechanism regulating SL, and promising gene source for the molecular breeding in rice.

Targeting the Metabolic Enzyme PGAM2 Overcomes Enzalutamide Resistance in Castration-Resistant Prostate Cancer by Inhibiting BCL2 Signaling.

The next-generation androgen receptor (AR) inhibitor enzalutamide is the mainstay treatment for metastatic prostate cancer. Unfortunately, resistance occurs rapidly in most patients, and once resistance occurs, treatment options are limited. Therefore, there is an urgent need to identify effective targets to overcome enzalutamide resistance. Here, using a genome-wide CRISPR-Cas9 library screen, we found that targeting a glycolytic enzyme, phosphoglycerate mutase PGAM2, significantly enhanced the sensitivity of enzalutamide-resistant prostate cancer cells to enzalutamide both in vivo and in vitro. Inhibition of PGAM2 together with enzalutamide treatment triggered apoptosis by decreasing levels of the antiapoptotic protein BCL-xL and increasing activity of the proapoptotic protein BAD. Mechanistically, PGAM2 bound to 14-3-3ζ and promoted its interaction with phosphorylated BAD, resulting in activation of BCL-xL and subsequent resistance to enzalutamide-induced apoptosis. In addition, high PGAM2 expression, which is transcriptionally regulated by AR, was associated with shorter survival and rapid development of enzalutamide resistance in patients with prostate cancer. Together, these findings provide evidence of a nonmetabolic function of PGAM2 in promoting enzalutamide resistance and identify PGAM2 inhibition as a promising therapeutic strategy for enzalutamide-resistant prostate cancer. SIGNIFICANCE: PGAM2 promotes resistance to enzalutamide by activating antiapoptotic BCL-xL and suppressing apoptosis, indicating that PGAM2 is a potential target for overcoming enzalutamide resistance in prostate cancer.

Long-Term Dietary Fish Meal Substitution with the Black Soldier Fly Larval Meal Modifies the Caecal Microbiota and Microbial Pathway in Laying Hens.

Feeding laying hens with black soldier fly larval (BSFL) meal improves their performance. However, the beneficial mechanism of BSFL meals in improving the performance of laying hens remains unclear. This study investigated the effects of the BSFL diet on liver metabolism, gut physiology, and gut microbiota in laying hens. Eighty-seven Julia hens were randomly assigned to three groups based on their diets and fed maize grain-and soybean meal-based diets mixed with either 3% fish meal (control diet), 1.5% fish and 1.5% BSFL meals, or 3% BSFL meal for 52 weeks. No significant differences were observed in biochemical parameters, hepatic amino acid and saturated fatty acid contents, intestinal mucosal disaccharidase activity, and intestinal morphology between BSFL diet-fed and control diet-fed laying hens. However, the BSFL diet significantly increased the abundance of acetic and propionic acid-producing bacteria, caecal short-chain fatty acids, and modified the caecal microbial pathways that are associated with bile acid metabolism. These findings indicate that consuming a diet containing BSFL meal has minimal effects on plasma and liver nutritional metabolism in laying hens; however, it can alter the gut microbiota associated with short-chain fatty acid production as well as the microbial pathways involved in intestinal fat metabolism. In conclusion, this study provides evidence that BSFL can enhance enterocyte metabolism and gut homeostasis in laying hens.

The 14-3-3 protein OsGF14f interacts with OsbZIP23 and enhances its activity to confer osmotic stress tolerance in rice.

Drought, which can induce osmotic stress, is the leading environmental constraint on crop productivity. Plants in both agricultural and natural settings have developed various mechanisms to cope with drought stress. The identification of genes associated with drought stress tolerance and understanding the underlying regulatory mechanisms are prerequisites for developing molecular manipulation strategies to address this issue. Here, we reported that the G-BOX FACTOR 14-3-3f (14-3-3 protein OsGF14f) positively modulates osmotic stress tolerance in rice (Oryza sativa). OsGF14f transgenic lines had no obvious change in crucial agronomic traits including yield and plant height. OsGF14f is transcriptionally induced by PEG treatment, and in rice, overexpression or knockout of this gene leads to enhanced or weakened osmotic stress tolerance, respectively. Furthermore, OsGF14f positively regulates abscisic acid (ABA) responses by interacting with the core ABA-responsive transcription factor BASIC LEUCINE ZIPPER 23 (OsbZIP23) to enhance its transcriptional regulation activity toward downstream target genes. Further genetic analysis showed that OsGF14f is required for the full function of OsbZIP23 in rice osmotic response, and OsGF14f-mediated osmotic stress tolerance partially depends on OsbZIP23. Interestingly, OsGF14f is a direct target gene of OsbZIP23. Taken together, our findings reveal a genetic and molecular framework by which the OsGF14f-OsbZIP23 complex modulates rice osmotic response, providing targets for developing drought-tolerant crops.

Two complementary genes in a presence-absence variation contribute to indica-japonica reproductive isolation in rice.

Understanding the evolutionary forces in speciation is a central goal in evolutionary biology. Asian cultivated rice has two subspecies, indica and japonica, but the underlying mechanism of the partial reproductive isolation between them remains obscure. Here we show a presence-absence variation (PAV) at the Se locus functions as an indica-japonica reproductive barrier by causing hybrid sterility (HS) in indica-japonica crosses. The locus comprises two adjacent genes: ORF3 encodes a sporophytic pollen killer, whereas ORF4 protects pollen in a gametophytic manner. In F1 of indica-japonica crosses, pollen with the japonica haplotype, which lacks the sequence containing the protective ORF4, is aborted due to the pollen-killing effect of ORF3 from indica. Evolutionary analysis suggests ORF3 is a gene associated with the Asian cultivated rice species complex, and the PAV has contributed to the reproductive isolation between the two subspecies of Asian cultivated rice. Our analyses provide perspectives on rice inter-subspecies post-zygotic isolation, and will promote efforts to overcome reproductive barriers in indica-japonica hybrid rice breeding.

Progress in the study of functional genes related to direct seeding of rice.

Rice is a major food crop in the world. Owing to the shortage of rural labor and the development of agricultural mechanization, direct seeding has become the main method of rice cultivation. At present, the main problems faced by direct seeding of rice are low whole seedling rate, serious weeds, and easy lodging of rice in the middle and late stages of growth. Along with the rapid development of functional genomics, the functions of a large number of genes have been confirmed, including seed vigor, low-temperature tolerance germination, low oxygen tolerance growth, early seedling vigor, early root vigor, resistance to lodging, and other functional genes related to the direct seeding of rice. A review of the related functional genes has not yet been reported. In this study, the genes related to direct seeding of rice are summarized to comprehensively understand the genetic basis and mechanism of action in direct seeding of rice and to lay the foundation for further basic theoretical research and breeding application research in direct seeding of rice.

Genome-wide association mapping and gene expression analysis reveal candidate genes for grain chalkiness in rice.

Grain chalkiness is the main factor determining the market value of rice. Reducing chalkiness is an important breeding goal for genetic improvement of high quality rice. Identification of QTLs or genes controlling chalkiness is the prerequisite for molecular breeding in rice. Here, we conducted a genome-wide association study to identify QTLs associated with grain chalkiness including percentage of grains with chalkiness (PGWC) and degree of endosperm chalkiness (DEC) in 450 rice accessions consisting of 300 indica and 150 japonica rice in two environments. A total of 34 QTLs were identified, including 14 QTLs for PGWC and 20 QTLs for DEC. Among them, seven QTLs were commonly identified in two environments, and eight QTLs were simultaneously related to two traits. Based on the haplotype analysis, LD decay analysis, RNA-sequencing, qRT-PCR confirmation and haplotype comparisons, four genes (LOC_Os10g36170, LOC_Os10g36260, LOC_Os10g36340 and LOC_Os10g36610) were considered as the candidate genes for qDEC-10c1w,2wj , which could be identified in both environments and had the most significant p-value among the newly identified QTLs. These results provided new insight into the genetic basis of grain chalkiness and gene resources for improving quality by molecular breeding in rice.

Gap-free genome assembly and comparative analysis reveal the evolution and anthocyanin accumulation mechanism of Rhodomyrtus tomentosa.

Rhodomyrtus tomentosa is an important fleshy-fruited tree and a well-known medicinal plant of the Myrtaceae family that is widely cultivated in tropical and subtropical areas of the world. However, studies on the evolution and genomic breeding of R. tomentosa were hindered by the lack of a reference genome. Here, we presented a chromosome-level gap-free T2T genome assembly of R. tomentosa using PacBio and ONT long read sequencing. We assembled the genome with size of 470.35 Mb and contig N50 of ~43.80 Mb with 11 pseudochromosomes. A total of 33 382 genes and 239.31 Mb of repetitive sequences were annotated in this genome. Phylogenetic analysis elucidated the independent evolution of R. tomentosa starting from 14.37MYA and shared a recent WGD event with other Myrtaceae species. We identified four major compounds of anthocyanins and their synthetic pathways in R. tomentosa. Comparative genomic and gene expression analysis suggested the coloring and high anthocyanin accumulation in R. tomentosa tends to be determined by the activation of anthocyanin synthesis pathway. The positive selection and up-regulation of MYB transcription factors were the implicit factors in this process. The copy number increase of downstream anthocyanin transport-related OMT and GST gene were also detected in R. tomentosa. Expression analysis and pathway identification enriched the importance of starch degradation, response to stimuli, effect of hormones, and cell wall metabolism during the fleshy fruit development in Myrtaceae. Our genome assembly provided a foundation for investigating the origins and differentiation of Myrtaceae species and accelerated the genetic improvement of R. tomentosa.