Plant somatic embryogenesis (SE) is an in vitro biological process wherein bipolar structures are induced to form somatic cells and regenerate into whole plants. MicroRNA (miRNA) is an essential player in plant SE. However, the mechanism of microRNA408 (miR408) in SE remains elusive. Here, we used stable transgenic technology in longan (Dimocarpus longan) embryogenic calli to verify the mechanism by which miR408 promotes cell division and differentiation of longan early SE. dlo-miR408-3p regulated riboflavin biosynthesis by targeting nudix hydrolase 23 (DlNUDT23), a previously unidentified gene mediating N6-methyladenosine (m⁶A) modification and influencing RNA homeostasis and cell cycle gene expression during longan early SE. We showed that DlMIR408 overexpression (DlMIR408-OE) promoted 21-nt miRNA biosynthesis. In DlMIR408-OE cell lines, dlo-miR408-3p targeted and downregulated DlNUDT23, promoted riboflavin biosynthesis, decreased flavin mononucleotide (FMN) accumulation, promoted m⁶A level, and influenced miRNA homeostasis. DNA replication, glycosylphosphatidylinositol (GPI)-anchor biosynthesis, the pentose phosphate pathway, and taurine and hypotaurine metabolism were also closely associated with riboflavin metabolism. In a riboflavin feeding assay, dlo-miR408-3p and pre-miR408 were upregulated and DlNUDT23 was downregulated, increasing the m⁶A level and cell division and differentiation in longan globular embryos. When riboflavin biosynthesis was inhibited, dlo-miR408-3p was downregulated and DlNUDT23 was upregulated, which decreased m⁶A modification and inhibited cell division but did not inhibit cell differentiation. FMN artificial demethylated m⁶A modification affected the homeostasis of precursor miRNA and miRNA. Our results revealed a mechanism underlying dlo-miR408-3p-activated riboflavin biosynthesis in which DlNUDT23 is targeted, m⁶A modification is dynamically mediated, and cell division is affected, promoting early SE in plants.

DNA methylation plays essential roles in gene regulation, chromatin structure stability, gene imprinting, X chromosome inactivation and embryonic development. However, the dynamics and functions of DNA methylation during the early stage of longan (Dimocarpus longan) somatic embryogenesis (SE) are still unclear. In this study, we carried out whole genome bisulphite sequencing and transcriptome sequencing analyses for embryogenic callus (EC), incomplete compact pro-embryogenic cultures (ICpEC) and globular embryos (GE) in an early SE system. At a global level, the DNA 5-methylcytosine content in EC, ICpEC and GE was 24.59, 19.65 and 19.74%, respectively, suggesting a global decrease in DNA methylation from EC to ICpEC and then a slight increase from ICpEC to GE. Differentially methylated region (DMR) analysis showed that hypomethylation mainly occurred in CHH contexts. Gene ontology and Kyoto encyclopedia of genes and genomes analysis of hypomethylated-CHH-DMR-associated genes revealed that zein biosynthesis, fatty acid biosynthesis, circadian rhythm and mitophagy pathways were involved in longan early SE. Expression patterns of DNA methyltransferase and demethylase genes during longan early SE suggested that the decrease in DNA methylation was probably regulated by DNA methyltransferase genes and the DNA demethylase gene REPRESSOR OF SILENCING 1 (ROS1). The correlation between DNA hypomethylation and gene expression revealed that decreased DNA methylation did not cause extensive changes in gene expression during early longan SE and that gene expression may be affected by methylation changes in gene and downstream regions. Inhibiting DNA methylation with 5-azacytidine treatment in EC promoted the formation of GE and enhanced the capability of longan SE. Our results suggest that DNA demethylation has important roles in longan SE development.

Somatic embryogenesis (SE), like zygotic embryo development, is a progressive process. Early SE is the beginning of a switch from a somatic to an embryogenic state and is an important stage for initiating chromatin reprogramming of SE. Previous studies suggest that changes in chromatin accessibility occur during early SE, although information on the 3D structure of chromatin is not yet available. Here, we present a chromosome-level genome assembly of longan (Dimocarpus longan) using PacBio combined with Hi-C scaffolding, which resulted in a 446 Mb genome assembly anchored onto 15 scaffolds. During early SE, chromatin was concentrated and then decondensed, and a large number of long terminal repeat retrotransposons (LTR-RTs) were enriched in the local chromatin interaction region, suggesting LTR-RTs were involved in chromatin reorganization. Early SE was accompanied by the transformation from A to B compartments, and the interactions between B compartments were enhanced. Results from chromatin accessibility, monomethylation of histone H3 at lysine 4 (H3K4me1) modification, and transcription analyses further revealed a gene regulatory network for cell wall thickening during SE. Particularly, we found that the H3K4me1 differential peak binding motif showed abnormal activation of ethylene response factor (ERF) transcription factors and participation in SE. The chromosome-level genomic and multi-omics analyses revealed the 3D conformation of chromatin during early SE, providing insight into the molecular mechanisms underlying cell wall thickening and the potential regulatory networks of TFs during early SE in D. longan. These results provide additional clues for revealing the molecular mechanisms of plant SE.

Circular RNAs (circRNAs) are widely involved in plant growth and development. However, the function of circRNAs in plant somatic embryogenesis (SE) remains elusive. Here, by using high-throughput sequencing, a total of 5029 circRNAs were identified in the three stages of longan (Dimocarpus longan Lour.) early SE. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that differentially expressed (DE) circRNA host genes were enriched in the ‘non-homologous end-joining’ (NHEJ) and ‘butanoate metabolism’ pathways. In addition, the reactive oxygen species (ROS) content during longan early SE was determined. The results indicated that ROS-induced DNA double-strand breaks may not depend on the NHEJ repair pathway. Correlation analyses of the levels of related metabolites (glutamate, γ-aminobutyrate and pyruvate) and the expression levels of circRNAs and their host genes involved in butanoate metabolism were performed. The results suggested that circRNAs may act as regulators of the expression of cognate mRNAs, thereby affecting the accumulation of related compounds. A competing endogenous RNA (ceRNA) network of DE circRNAs, DE mRNAs, DE long noncoding RNAs (lncRNAs) and DE microRNAs (miRNAs) was constructed. The results showed that the putative targets of the noncoding RNA (ncRNAs) were significantly enriched in the KEGG pathways ‘mitogen-activated protein kinase signaling’ and ‘nitrogen metabolism’. Furthermore, the expression patterns of the candidate circRNAs, lncRNAs, miRNAs and mRNAs confirmed the negative correlation between miRNAs and ceRNAs. In addition, two circRNA overexpression vectors were constructed to further verify the ceRNA network correlations in longan early SE. Our study revealed the potential role of circRNAs in longan early SE, providing new insights into the intricate regulatory mechanism underlying plant SE.

Anthracnose, a prevalent fungal disease in tea plantations, cause substantial economic losses in tea production. Identifying resistance-associated genes in tea plants is crucial for developing anthracnose-resistant cultivars. This study used eight tea samples with differential anthracnose resistance for phenotypic evaluation, weighted gene co-expression network analysis (WGCNA) of RNA-seq data, WGCNA- QTL co-localization to identify resistance gene, and qRT-PCR validation of candidate genes. in vitro pathogen inoculation assay revealed that the lesion diameters of the eight samples ranged from 1.45 mm to 4.5 mm (ANOVA p = 4.4 ×10−10 ). Using the ‘Longjing 43’ reference genome, transcriptome assembly achieved 93.9% gene detection rate (31,509/33,557 genes). WGCNA categorized expressed genes into 30 modules with the purple module (containing 907 genes) showing positive trait correlation and the yellow-green module (containing 781 genes) exhibiting negative correlation. Integration of WGCNA and QTL mapping identified two high-confidence candidate genes within LG08 QTL intervals. Both genes exhibited significant upregulation (t-test p < 0.01) in tea plant leaves following Colletotrichum spore inoculation. These findings provide actionable genetic targets for marker-assisted breeding of anthracnose-resistant tea cultivars.

Polyamines (PAs) play an important regulatory role in many basic cellular processes and physiological and biochemical processes. However, there are few studies on the identification of PA biosynthesis and metabolism family members and the role of PAs in the transition of plant embryogenic calli (EC) into globular embryos (GE), especially in perennial woody plants. We identified 20 genes involved in PA biosynthesis and metabolism from the third-generation genome of longan (Dimocarpus longan Lour.). There were no significant differences between longan and other species regarding the number of members, and they had high similarity with Citrus sinensis. Light, plant hormones and a variety of stress cis-acting elements were found in these family members. The biosynthesis and metabolism of PAs in longan were mainly completed by DlADC2, DlSAMDC2, DlSAMDC3, DlSPDS1A, DlSPMS, DlCuAOB, DlCuAO3A, DlPAO2 and DlPAO4B. In addition, 0.01 mmol∙L⁻¹ 1-aminocyclopropane-1-carboxylic acid (ACC), putrescine (Put) and spermine (Spm), could promote the transformation of EC into GE, and Spm treatment had the best effect, while 0.01 mmol∙L⁻¹ D-arginine (D-arg) treatment inhibited the process. The period between the 9th and 11th days was key for the transformation of EC into GE in longan. There were higher levels of gibberellin (GA), salicylic acid (SA) and abscisic acid (ABA) and lower levels of indole-3-acetic acid (IAA), ethylene and hydrogen peroxide (H₂O₂) in this key period. The expression levels in this period of DlADC2, DlODC, DlSPDS1A, DlCuAOB and DlPAO4B were upregulated, while those of DlSAMDC2 and DlSPMS were downregulated. These results showed that the exogenous ACC, D-arg and PAs could regulate the transformation of EC into GE in longan by changing the content of endogenous hormones and the expression levels of PA biosynthesis and metabolism genes. This study provided a foundation for further determining the physicochemical properties and molecular evolution characteristics of the PA biosynthesis and metabolism gene families, and explored the mechanism of PAs and ethylene for regulating the transformation of plant EC into GE.