Epitranscriptomic mark N6-methyladenosine (m6A) is the most predominant internal modification in RNAs, which plays pivotal roles in response to diverse stresses. Multiple environmental stresses caused by withering process can greatly influence the accumulation of specialized metabolites and the formation of tea flavor. However, little is known about the effects of m6A-mediated regulatory mechanism on flavor-related metabolisms in tea leaves. Here, we explored m6A-mediated regulatory mechanism and its impacts on flavonoid and terpenoid metabolisms under solar-withering using integrated RNA methylome and transcriptome. Dynamic changes in global m6A levels of tea leaves are mainly controlled by two m6A erasers (CsALKBH4A and CsALKBH4B) under solar-withering. Differentially methylated peak (DMP)-associated genes under different shading rates of solar-withering were identified and found to be enriched in terpenoid biosynthesis and spliceosome pathways. Further analyses indicated that CsALKBH4-driven RNA demethylation can not only directly affect the accumulation of volatile terpenoids by mediating the stability and abundance of terpenoid biosynthesis-related genes, but also indirectly influence the contents of flavonoids, catechins, and theaflavins via triggering the alternative splicing (AS)-mediated regulation. Our findings underscored a novel layer of epitranscriptomic gene regulation in tea flavor-related metabolic pathways and established a compelling link between m6A-mediated regulatory mechanism and the formation of high-quality flavor in tea leaves under solar-withering.

Recessive genic male sterility (RGMS) provides an effective approach for the commercial exploitation of heterosis, especially in Brassicacrops. Although some artificial RGMS mutants have been reported in B. rapa, no causal genes derived from these natural mutants have been identified so far. In this study, a spontaneous RGMS mutant Bcajh97-01A derived from the ‘Aijiaohuang’ line traced back to the 1980s was identified. Genetic analysis revealed that the RGMS trait was controlled by a single locus in the Bcajh97-01A/Bsystem. Bulk segregant analysis (BSA) in combination with linkage analysis was employed to delimit the causal gene to an approximate 129 kb interval on chromosome A02. The integrated information of transcriptional levels and the predicted genes in the target region indicated that the Brmmd1 (BraA02g017420) encoding a PHD-containing nuclear protein was the most likely candidate gene. In the RGMS line, a 374 bp miniature inverted-repeat transposable element (MITE) was inserted into the first exon to prematurely stop the Brmmd1 gene translation, thus blocking the normal expression of this gene at the tetrad stage in the Bcajh97-01A. Additionally, a co-segregating structure variation (SV) marker was developed to rapidly screen the RGMS progenies from Bcajh97-01A/B system. Our findings reveal that BraA02g017420 is the causal gene responsible for the RGMS trait. This study lays a foundation for marker-assisted selection and further molecular mechanism exploration of pollen development in B. rapa.