Plant Communications

Arabidopsis (Arabidopsis thaliana) is a vital model organism in plant biology and genetics. In the last two decades, researchers have made significant progresses in characterizing the chromatin conformation and epigenomic information within the Arabidopsis genome. This information includes but is not limited to the higher structure of chromosomes, histone modification, DNA methylation, and chromatin accessibility. The results of these studies have provided an additional layer of information that complements the DNA sequence data. However, utilizing such knowledge poses a challenge for certain groups that lack bioinformatics analysts or adequate computing resources. A user-friendly and reproducible platform for accessing this information is urgently needed. In this study, we have developed a comprehensive epigenomic database for Arabidopsis (AraENCODE http://glab.hzau.edu.cn/AraENCODE), which comprises a total of 4,511 data libraries, including published chromatin conformation capture datasets (Hi-C, HiChIP), epigenomic datasets (ChIP-Seq, ATAC-Seq, FAIRE-Seq, MNase-Seq, DNase-Seq, BS-seq), and transcriptome data (RNA-Seq, miRNA-Seq). Furthermore, we have incorporated various existing resources, such as single nucleotide polymorphisms (SNPs), cis-regulatory modules, and multi-omics associations. We aim to provide a novel platform for investigating the regulation of epigenetic and chromatin interactions in Arabidopsis in relation to biological processes.

Cassava is the most important starch sources, a tropical model crop. We constructed nearly T2T genomes of cultivar AM560, wild ancestors FLA4047 and W14, pan-genome of 24 representatives and a clarified evolutionary tree with 486 accessions. Comparison of SVs and SNVs between the ancestors and cultivated cassavas revealed predominant expansion, contraction of genes and gene families. Significantly selective sweeping occurred in the cassava genomes in 122 footprints with 1,519 candidate domestication genes. We identified selective mutations in MeCSK and MeFNR3 promoting photoreaction associated with MeNADP-ME of C4 assimilation in modern cassava. Co-evolved retardation of floral primordia and initiation of storage roots arose from MeCOL5 mutants with altered bindings to MeFT1, MeFT2 and MeTFL2. MebHLHs evolved to regulate the biosynthesis, transport and endogenous remobilization of cyanogenic glucosides, with new functionalities of MeMATE1, MeGTR in selected sweet cassava. These findings enhanced comprehensive knowledge and database on the evolution and breeding of cassava. HIGHLIGHTSO_LIThree nearly T2T cassava genomes of cultivar AM560 and its wild ancestors FLA4047 and W14. C_LIO_LIA species-level cassava panSV haplotype map across 346,322 structural variations over 31,362 gene families and 96,032,008 SNPs and InDels variations globally and a clarified evolutionary tree with 486 accessions. C_LIO_LISelective mutations in MeCSK and MeFNR3 promoted photoreaction associated with MeNADP-ME of C4 assimilation shaped the C3-C4 intermediate photosynthesis of modern cassava. C_LIO_LICoevolution of floral primordia contrary to initiated storage root is pivotal for the domestication of cassava, and arose from MeCOL5 mutants altered the binding with MeFT1, MeFT2 (SP6A), and MeTFL2. C_LI

Protein-protein interactions (PPIs) play a crucial role in enabling proteins to carry out their functions within various biological processes (Hui et al., 2003). Since the introduction of the yeast two-hybrid (Y2H) method for PPI detection in 1989 (Fields and Song, 1989), the identification of PPIs has become a significant focus in modern biological research. PPI goes beyond examining individual proteins, allowing researchers to establish a comprehensive network that regulates biological processes. Rice, as a key model organism in plant biological studies, has been at the forefront of PPI research. In 2008, prominent rice scientists in China called for concerted efforts to define a comprehensive protein-protein interaction network experimentally, which aimed to facilitate the prediction of the functional mechanisms operating throughout a plants lifecycle (Zhang et al., 2008). With efforts for 2 decades, the experimentally identified rice PPIs have reached over ten thousand. Several public databases have been established to systematically collate and store PPIs, including STRING (Szklarczyk et al., 2019), BioGRID (Oughtred et al., 2020), IntAct (del Toro et al., 2022), PRIN (Gu et al., 2011), RicePPINet (Liu et al., 2017) and RiceNet v2 (Lee et al., 2015). However, most PPI datasets in rice stem from computational predictions, while experiment-based rice PPI datasets are fragmented due to the lack of systematic profiling at the rice PPIome level, which largely hinders information sharing in the rice research community. To bridge this gap, we constructed the Port of Protein-Protein Interactomes (POPPIN; https://riceome.hzau.edu.cn/poppin/), an integrated database dedicated to sharing experimentally verified PPIs and functional clues in rice. Empowered by high-throughput PPIome profiling technologies and text mining assisted by a large language model (Huang et al., 2025; Liu et al., 2025), POPPIN currently has deposited over 150,451 pieces of rice PPI-related information. Additionally, POPPIN provides detailed protein information, including GO annotations, subcellular localizations, domains, trait ontology (TO) information, and hyperlinks to external biological databases. Through offering a user-friendly web interface for search and dynamic network visualization, POPPIN serves as the first large-scale, experiment-based database for searchable PPIs in rice, and has the potential to be extended to other species under this structural framework.

Dear Editor, Since the first plant receptor-like kinase (RLK) gene ZmPK1 was cloned from Zea mays in 1990 (Walker & Zhang, 1990), this large gene family has been extensively studied and shown to play crucial roles in growth, development, and immunity (Tang et al., 2017). RLKs are widespread in the plant kingdom, while the biological functions of most RLKs remain largely elusive (Dievart et al., 2020). Given RLKs share a conserved monophyletic RLK/Pelle kinase domain, RLKs in several model plants are classified into distinct families by extracellular domains (ECDs) (Shiu & Bleecker, 2001). However, independent domain shuffling in specific lineages drives the origin of novel families, which raises a question: how about the landscape of RLKs in the whole plant kingdom? Previously, sequence homology-based methods have been widely used for RLK identification and classification, which probably will miss the distantly related proteins but with similar structures and potential novel families unmentioned in the literature. The academic community urgently requires a dedicated database for a systematic overview of the RLK gene family, providing data support for in-depth research on RLK genes. Here, we used a topology-based method to accurately isolate the RLKomes from proteomes. The obtained RLKomes were further classified into (sub)families based on ECD domains. We constructed a comprehensively curated plant RLK database (https://biotec.njau.edu.cn/rlkdb/), which contains valuable resources for investigating the origin and evolution of the RLK family and multiple online tools for personalized analysis.

In agriculture, propiconazole (PCZ) controls excessive growth in flowering Chinese cabbage but poses dietary safety risks due to residue accumulation. Therefore, identifying novel PCZ targets and breeding PCZ-free cultivars is critical for the safe production of flowering Chinese cabbage. Here, we identified three P4-ATPase flippase homologs aminophospholipid ATPase 3 (BraALA3a/b/c) in flowering Chinese cabbage that function as sensitive targets for PCZ. These proteins exhibit high binding affinity for PCZ, which directly inhibits their ATPase activity. Overexpression of the BraALA3 homologs enhanced plant growth and increased sensitivity to PCZ, whereas knockdown led to dwarfism and reduced sensitivity. Based on these findings, we identified editable active sites via protoplast-based screening. Genetic transformation of one such site yielded BraALA3a/braala3aK200T mutant lines, which displayed a dwarf and compact architecture. These findings provide a precise molecular target for developing PCZ-free germplasm in flowering Chinese cabbage through gene editing.

The increasing availability of omics data from non-model plant species has created a pressing need for centralized, user-friendly platforms that maximize the utility of these datasets in a FAIR manner. Here, we introduce PlantaeViz, a web-based tool designed for the integration, visualization, and analysis of multi-omics data across a wide range of plant species. PlantaeViz offers advanced functionalities, including gene catalogues built from curated literature, transcriptomic meta-analyses presented as gene expression atlases, gene co-expression and regulatory networks with on-the-fly ontology analyses, cistrome visualization and metabolomics-transcriptomics integration, among other tools, providing a robust framework for hypothesis generation and biological interpretation. Gene Cards applications, tailored to each plant species, provide both community-curated and automatically generated functional annotation information. One of the platforms core features is its big data approach: over 58,000 publicly available SRA transcriptomic samples have been processed and visualized to date. Significant efforts have been made in orthology assessment using multiple layers of evidence, as well as in the automatic classification and standardization of omics metadata through regular expressions and data mining. As a result, around 90% of transcriptomic runs have been successfully classified according to sample tissue. These data have been used to construct gene networks via computationally intensive methods based on diverse algorithms. We present a study case to illustrate the platforms integration and exploratory capabilities. PlantaeViz bridges genomics and functional knowledge between model and non-model plant species and aims to expand its species catalogue of species in the future, democratizing access to large-scale plant omics data. Further developments will include the incorporation of additional data types, and the implementation of new tools to further support plant research across diverse biological contexts.

Momilactone A, an important plant labdane-related diterpenoid, functions as a phytoalexin against pathogens and an allelochemical against neighboring plants. The genes involved in biosynthesis of momilactone A are found in clusters, i.e., MABGCs (Momilactone A biosynthetic gene clusters), in the rice and barnyardgrass genomes. How MABGCs originate and evolve is still not clear. Here, we integrated results from comprehensive phylogeny and comparative genomic analyses of the core genes of MABGC-like clusters and MABGCs in 40 monocot plant genomes, providing convincing evidence for the birth and evolution of MABGCs in grass species. The MABGCs found in the PACMAD clade of the core grass lineage (including Panicoideae and Chloridoideae) originated from a MABGC-like cluster in Triticeae (BOP clade) via horizontal gene transfer (HGT) and followed by recruitment of MAS and CYP76L1 genes. The MABGCs in Oryzoideae originated from PACMAD through another HGT event and lost CYP76L1 afterwards. The Oryza MABGC and another Oryza diterpenoid cluster c2BGC are two distinct clusters, with the latter being originated from gene duplication and relocation within Oryzoideae. Further comparison of the expression patterns of the MABGC genes between rice and barnyardgrass in response to pathogen infection and allelopathy provides novel insights into the functional innovation of MABGCs in plants. Our results demonstrate HGT-mediated origination of MABGCs in grass and shed lights into the evolutionary innovation and optimization of plant biosynthetic pathways.

Sapindaceae is a family of flowering plants, also known as the soapberry family, comprising 141 genera and about 1900 species (Pedro et al., 2010). Most of them are distributed in tropical and subtropical regions, including trees, shrubs, also woody or herbaceous vines. Some are dioecious, while others are monoecious. Many Sapindaceae species possess great economic value; some furnish delicious fruits, like lychee (Litchi chinensis), longan (Dimocarpus longan), rambutan (Nephelium lappaceum); and ackee (Blighia sapida) - the national fruit of Jamaica; some produce abundance secondary metabolites, like saponin from soapberry (Sapindus mukorossi), and seed oil from yellowhorn (Xanthoceras sorbifolium); some yield valuable timber including maple (Acer spp.) and buckeye (Aesculus glabra); and some are of great herbal medicinal value, like balloon-vine (Cardiospermum halicacabum). In the last decade, with the rocketing of next generation sequencing (NGS) and genomic technologies, the full genome sequences of several Sapindaceae plants have been resolved (Lin et al., 2017; Liang et al., 2019; Yang et al., 2019; Zhang et al., 2021; Hu et al., 2022; Xue et al., 2022). Among them, our recent publication of the lychee genome attracted broad attention (Edger, 2022; Hu et al., 2022; Lyu, 2022). Now the post-genome era arrives for Sapindaceae, however, there is no public genomic database available for any Sapindaceae species, let alone an integrative database for the whole Sapindaceae family. A unified data platform is in urgent need to collect, manage and share relevant data resources. Therefore, we integrated our home-brew NGS data with all publicly available data for seven Sapindaceae plants and constructed the Sapinaceae Genomic DataBase, named SapBase (www.sapindaceae.com), in order to provide genomic resources and an online powerful analytic platform for scientific research on Sapinaceae species and comparative studies with other plants.

Analysis of spatiotemporal patterns of gene expression is crucial to decode biological systems responses. High-throughput sequencing allows in-depth transcriptome analyses and experimental designs, providing valuable reference expression atlases. Specifically, testing overrepresentations of tissue-specific genes based on these atlases can provide valuable insights; however, such an approach is not accessible to inexperienced users. Here, we introduce TOTEM (TOol for Tissue-EnrichMent), a web tool designed to calculate enrichment values per tissue by identifying tissue-specific genes from an organ/organism of interest given a user gene list. Results are visually represented, and users gene classified. The utility of TOTEM is manifest when using integrated single cell expression atlases, enabling the study of complicated tissues, with the maximum possible resolution. Its effectiveness is validated by the study of BRL3 role in stress specifically from the vascular tissues. Finally, TOTEMs modular design allows for continual integration of new experiments. TOTEM can be freely accessed at: https://totemwebtool.com.

Gene duplications have long been recognized as a driving force in the evolution of genes, giving rise to novel functions. The soybean genome is characterized by a large extent of duplicated genes. However, the extent and mechanisms of functional divergence among these duplicated genes in soybean remain poorly understood. In this study, we revealed that tandem duplication of MYB genes, which occurred specifically in the Phaseoleae lineage, exhibited a stronger purifying selection in soybean compared to common bean. To gain insights into the diverse functions of these MYB genes in anthocyanin biosynthesis, we examined the expression, transcriptional activity, metabolite, and evolutionary history of four MYB genes (GmMYBA5, GmMYBA2, GmMYBA1 and Glyma.09g235000), which were presumably generated by tandem duplication in soybean. Our data revealed that Glyma.09g235000 had become a pseudogene, while the remaining three MYB genes exhibited strong transcriptional activation activity and promoted anthocyanin biosynthesis in different soybean tissues. Furthermore, GmMYBA5 produced distinct compounds in Nicotiana benthamiana leaves compared to GmMYBA2 and GmMYBA1 due to variations in their DNA binding domains. The lower expression of anthocyanin related genes in GmMYBA5 resulted in lower levels of anthocyanins compared to GmMYBA2 and GmMYBA1. Metabolomics analysis further demonstrated the diverse and differential downstream metabolites, suggesting their functional divergence in metabolites following gene duplication. Together, our data provided evidence of functional divergence within the MYB gene cluster following tandem duplication, which shed light on the potential evolutionary direction of gene duplications during legume evolution.

Tinospora sagittata (Oliv.) Gagnep. is an important medicinal tetraploid plant in the Menispermaceae family. Its tuber, namely "Radix Tinosporae" used in Traditional Chinese Medicine, is rich in medicinal terpenoids and benzylisoquinoline alkaloids (BIAs), To enhance understanding the biosynthesis of medicinal compounds, we, herein, report the assembly of a high quality chromosome-scale genome with both PacBio HiFi and Illumina sequencing technologies. The size of assembled genome was 2.33 Gb consisting of 4070 scaffolds (N50=42.06Mb), of which 92.05% were assigned to 26 pseudochromosomes in A and B sub-genomes. A phylogenetic analysis with the T. sagittata and other 16 plant genomes estimated the evolutionary placement of T. sagittata and its divergence time in Ranunculales. Further genome evolution analysis characterized one round tandem duplication about 1.5 million years ago (MYA) and one whole-genome duplication (WGD) about 86.9 MYA. WGD contributed to the duplication of clade-specific cytochrome P450 gene family in Ranunculales. Moreover, sequencing mining obtained genome-wide genes involved in the biosynthesis of alkaloids and terpenoids. TsA02G014550, one candidate, was functionally characterized to catalyze the formation of (S)-canadine in the jatrorrhizine biosynthetic pathway. Taken together, the assembled genome of T. sagittata provides useful sequences to understand the biosynthesis of jatrorrhizine and other BIAs in plants.

With the rapid accumulation of diverse omics datasets, achieving efficient management and integrative analysis of plant multi-omics data remains a major challenge. Conventional solutions rely on constructing web-based databases, which often demand substantial programming expertise and long-term financial support. To address these limitations, we developed the Plant Multi-omics Database Construction System (PlantMDCS)-a locally deployable, user-friendly, and collaborative platform that unifies database construction and downstream multi-omics analysis within a graphical environment. PlantMDCS adopts a decoupled front-end/back-end architecture. The back end serves as the core engine for data management and computation, and is responsible for the storage, preprocessing, integration, and hierarchical association of multi-omics data. Once initialized, the front end supports the complete research workflow, including data import, querying, integrative analysis and visualization. All operations can be performed without programming, while local resource usage is dominated by disk storage required for user-provided datasets rather than sustained computational overhead. Benchmarking across plant species ranging from Arabidopsis to hexaploid wheat demonstrated that database construction can be completed within minutes, independent of genome size or data complexity. PlantMDCS is designed for local deployment to ensure data security, while allowing multi-user collaboration within local networks and supporting controlled remote access for teams distributed across different regions. Overall, PlantMDCS offers a secure and sustainable framework that integrates data management and analysis within a unified system. This design shifts multi-omics research away from fragmented file-based processing toward persistent, database-driven exploration, thereby enhancing analytical efficiency and reproducibility.

Senna tora (S. tora) is a commonly used Chinese medicinal plant due to the presence of the bioactive compounds anthraquinones in its mature seeds. Seed size is an important factor that affect S. tora yield quantity and quality. However, the mechanism regulating seed size and aurantio-obtusin biosynthesis in S.tora is still unclear. In this study, we identified the ubiquitin receptor StDA1 in S.tora that served as a negative regulator to seed formation and limited seed enlargement. Antisense overexpression of StDA1 led to larger seeds in S. tora and promoted the accumulation of aurantio-obtusin. In contrast, overexpression of StDA1 in S.tora resulted in a relative decrease in aurantio-obtusin accumulation. Moreover, StDA1 can directly bind to StHDR1and regulate its degradation through the 26S proteasome to regulate seed size and aurantio-obtusin accumulation. We also found that the StDA1-StHDR1 module is responsive to the MeJA via StTCP4.1, which in turn affects the accumulation of aurantio-obtusin. Overall, we have identified a protein complex that regulates the accumulation of aurantio-obtusin, StTCP4.1-StDA1-StHDR1, as a potential target for improving S.tora yield quantity and quality.

Bulk RNA sequencing and single-cell RNA sequencing provide complementary information on tissue and cell-type-specific gene expression. Bulk RNA sequencing enables the construction of gene association networks that identify co-expressed genes involved in shared pathways, whereas single-cell RNA sequencing maps their expression to cell types. However, most platforms provide access to either bulk RNA sequencing or single-cell RNA sequencing analysis, making it difficult to connect tissue-level co-expression with cell-type-specific expression. PlantNetX was developed as a web-based platform that integrates both data types. Although PlantNetX currently focuses on rice (Oryza sativa) and includes 70 quality-controlled RNA sequencing datasets comprising 1,198 sequencing libraries, together with nine single-cell RNA sequencing datasets containing more than 580,000 cells, including recently released datasets not consistently represented in existing platforms, it was designed to incorporate additional plant species, datasets, and analytical tools. PlantNetX provides Mutual Rank-based co-expression analysis, global and tissue-specific gene association networks, interactive visualization, and cell-type-specific expression summaries. The platform was validated with published examples of plant cell-wall biosynthesis and root-hair growth and retrieved gene association and expression patterns. Under standardized testing conditions, PlantNetX had a shorter mean response time than the other databases assessed. PlantNetX will support research in plant cell-wall biosynthesis, pathway discovery, functional genomics, and crop improvement. HighlightsPlantNetX integrates gene co-expression networks with cell-type expression, enabling fast identification and biological interpretation of candidate genes across tissues and individual cells. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=120 SRC="FIGDIR/small/741827v1_ufig1.gif" ALT="Figure 1"> View larger version (60K): org.highwire.dtl.DTLVardef@1207e06org.highwire.dtl.DTLVardef@31c3acorg.highwire.dtl.DTLVardef@12560b0org.highwire.dtl.DTLVardef@eeee7c_HPS_FORMAT_FIGEXP M_FIG C_FIG

With the development of high-throughput experimental technologies, large-scale RNA sequencing (RNA-Seq) data have been and continue to be produced, but have led to challenges in extracting relevant biological knowledge hidden in the produced high-dimensional gene expression matrices. Here, we present easyMF, a user-friendly web platform that aims to facilitate biological discovery from large-scale transcriptome data through matrix factorization (MF). The easyMF platform enables users with little bioinformatics experience to streamline transcriptome analysis from raw reads to gene expression and to decompose expression matrix from thousands of genes to a handful of metagenes. easyMF also offers a series of functional modules for metagene-based exploratory analysis with an emphasis on functional gene discovery. As a modular, containerized and open-source platform, easyMF can be customized to satisfy users specific demands and deployed as a web server for broad applications. easyMF is freely available at https://github.com/cma2015/easyMF. We demonstrated the application of easyMF with four case studies using 940 RNA sequencing datasets from maize (Zea mays L.).

Grain yield and salt tolerance are critical for crop production. However, the genetic and biochemical basis underlying the trade-off of these characters remain poorly described in crops. We show here that SiPLATZ12 transcription factor positively regulates multiple elite yield traits at the expense of salt tolerance in foxtail millet. SiPLATZ12 overexpression increases seed size, panicle length, and stem diameter, while reduces plant height and salt tolerance of foxtail millet. A 9-bp insertion in the SiPLATZ12 promoter has significant effects on the different expression of SiPLATZ12, multiple yield traits, and salt tolerance between foxtail millet and its wild ancestor, green foxtail. Moreover, SiPLATZ12 upregulates the expression of genes involved in seed development, but repressing the transcription of most NHX, SOS, and CBL genes to regulate Na+, K+ and pH homeostasis. Therefore, our results uncover a domesticated site that could be used to improve grain yield and salt tolerance in foxtail millet.

The development of robust biomarkers enables accurate prediction of complex phenotypes. However, the dynamic nature of biomarkers is often underestimated since their quantitative changes during development are directly connected to phenotypic transformations, influencing both crop agronomic traits and human diseases. Here, we performed network analysis of untargeted metabolite profiles to investigate height heterosis in rice, which is dynamic that varies during development and is a key determinant of yield heterosis. We found that the levels of pyruvaldehyde were predictive of height heterosis specific at the seedling stage, while 4-hydroxycinnamic acid positively correlated with height heterosis across four developmental stages. We identified metabolic pathways associated with height heterosis and found that metabolomic changes during the elongation stage had a greater impact than those in other stages. Finally, 11 heterosis-associated pathways were developed into metabolomic biomarkers through random forest analysis, successfully predicting height heterosis in an independent population under different growth conditions. This study elucidates the metabolomic landscape of dynamic height heterosis in rice and develops pathway biomarkers for complex phenotypes, demonstrating robustness across diverse populations, environments, and developmental stages.

The genus Citrus is characterized by a reticulate evolutionary history with frequent hybridization, making it an intriguing subject for genome evolution investigation. Citrus changshan Y. B. Chang (Huyou) is a unique landrace first discovered in Zhejiang Province, China with premium fruit quality. The evolutionary origin of Huyou has puzzled local botanists and growers. Here, we sequenced a 120-years-old "ancestral tree" of Huyou using PacBio long read and Hi-C sequencing and assembled 2 high-quality haplotype-resolved genomes HY1 and HY2. Huyou displayed a genome heterozygosity level at 3.07%, among the highest in published citrus genomes. Using a k-mer-based tracing approach, we explicitly resolved that HY1 genome contained 87.8% mandarin, 7.3% pummelo, 0.2% citron origin, whereas HY2 had 85.0% pummelo, 2.9% mandarin, 0.3% citron, implying a hybridization event between mandarin and pummelo. Phylogeny dating showed that HY1 (2.0 Mya) and HY2 (2.18 Mya) had diverged earlier than the split of Citrus clementina and Citrus reticulata, and the split of Citrus grandis and Citrus maxima, respectively. We observed clear chromosomal recombination on chr8 and chr9 in HY1, which may have occurred after the ancestral hybridization. Further transcriptome analyses in 6 tissues revealed a strong allelic dominance of HY2 over HY1 in root tissue and moderately in stem, leaf, flower, and fruits. KEGG enrichment analyses revealed that genes related to antioxidants biosynthesis and lipid metabolisms were most significantly affected by allelic imbalance. This first report of allelic imbalance in citrus species support Huyou as an interesting model to investigate genome evolution following distant hybridization.

Achieving high-throughput and precise phenotypic quantification and imaging modalities of stomatal and epidermal cells across diverse species remains a primary bottleneck in elucidating the mechanisms of stomatal dynamics, epidermal patterning, and environmental adaptation of plants. Here, we developed EpiReasoner, an artificial intelligence framework comprising a vision module, EpiVision, and a knowledge-based reasoning module, EpiBrain, for the quantitative phenotypic analysis and domain-specific knowledge reasoning of stomatal complexes and pavement cells in plants. Operating across bright-field, scanning electron microscopy, and differential interference contrast modalities, EpiVision achieves precise instance segmentation in various monocotyledonous, dicotyledonous, and fern species. Its performance significantly surpasses current state-of-the-art models. Moreover, we defined 23 quantitative indices describing stomatal cell morphology and spatial distribution. For domain-specific tasks such as phenotype prediction, genotype deduction, and molecular mechanism reasoning, EpiBrain demonstrates a human preference rate significantly higher than that of general-purpose large language models, including GPT-5 and Claude Sonnet 4. The application of EpiReasoner to phenotypic data of stomatal density derived from a tomato natural population of 170 accessions successfully identified a major quantitative trait locus on chromosome 8. The candidate gene, SKP1-interaction partner 19L (SKIP19L), encoding an F-box family protein, exhibited severe allele frequency drift during tomato domestication, which is highly consistent with the adaptive trend of reduced stomatal density under artificial selection. EpiReasoner provides a novel paradigm that unifies visual phenomics and knowledge-driven reasoning for the biology of stomata and pavement cells, thereby significantly accelerating scientific discovery in plant science.

Pandanus amaryllifolius is a valuable aromatic crop widely cultivated in tropical forest understories. Here, A telomere-to-telomere, gap-free genome assembly of the Chinese cultivar Taishan Banlan was generated, and the molecular basis underlying shading-mediated regulation of leaf aroma biosynthesis was elucidated for the first time through integrated genomic, transcriptomic, and metabolomic analyses. Volatile metabolomics identified 1599 volatile organic compounds in leaf tissue, with terpenoids as the primary contributors to leaf aroma. Shading modulated the aroma quality of P. amaryllifolius primarily via terpenoid biosynthesis, yielding aroma profiles dominated by sweet, fruity, green, and woody notes. Integrated transcriptomic profiling revealed that terpenoid biosynthesis is co-regulated by structural genes in the MVA/MEP pathways, TPS family members, and key transcription factors, which collectively modulate the biosynthesis and accumulation of the three core differential terpenoids. Our findings provide valuable resources for the targeted genetic improvement and efficient cultivation optimization of P. amaryllifolius.