Plant multiscale networks: charting plant connectivity by multi-level analysis and imaging techniques

Unlocking the Potential of the RUBY Reporter System: How to Address Its Challenges in Plant-Environment Interaction Research?

We summarized limitations of RUBY and proposed using a regulatable promoter system to control RUBY expression, thereby alleviating the metabolic stress. We aim to inspire scientists to optimize this reporter system, thereby applying it to study plant‐environment interactions in the future.

Large-volume fully automated cell reconstruction generates a cell atlas of plant tissues

The geometric shape and arrangement of individual cells play a role in shaping organ functions. However, analyzing multicellular features and exploring their connectomes in centimeter-scale plant organs remain challenging. Here, we established a set of frameworks named Large-Volume Fully Automated Cell Reconstruction (LVACR), enabling the exploration of three-dimensional (3D) cytological features and cellular connectivity in plant tissues. Through benchmark testing, our framework demonstrated superior efficiency in cell segmentation and aggregation, successfully addressing the inherent challenges posed by light sheet fluorescence microscopy (LSFM) imaging. Using LVACR, we successfully established a cell atlas of different plant tissues. Cellular morphology analysis revealed differences of cell clusters and shapes in between different poplar (P. simonii Carr. and P. canadensis Moench.) seeds, whereas topological analysis revealed that they maintained conserved cellular connectivity. Furthermore, LVACR spatiotemporally demonstrated an initial burst of cell proliferation, accompanied by morphological transformations at an early stage in developing the shoot apical meristem. During subsequent development, cell differentiation produced anisotropic features, thereby resulting in various cell shapes. Overall, our findings provided valuable insights into the precise spatial arrangement and cellular behavior of multicellular organisms, thus enhancing our understanding of the complex processes underlying plant growth and differentiation.

Three-dimensional reconstruction and multiomics analysis reveal a unique pattern of embryogenesis in Ginkgo biloba

Ginkgo (Ginkgo biloba L.) is one of the earliest extant species in seed plant phylogeny. Embryo development patterns can provide fundamental evidence for the origin, evolution, and adaptation of seeds. However, the architectural and morphological dynamics during embryogenesis in G. biloba remain elusive. Herein, we obtained over 2,200 visual slices from 3 stages of embryo development using micro-computed tomography imaging with improved staining methods. Based on 3-dimensional (3D) spatiotemporal pattern analysis, we found that a shoot apical meristem with 7 highly differentiated leaf primordia, including apical and axillary leaf buds, is present in mature Ginkgo embryos. 3D rendering from the front, top, and side views showed 2 separate transport systems of tracheids located in the hypocotyl and cotyledon, representing a unique pattern of embryogenesis. Furthermore, the morphological dynamic analysis of secretory cavities indicated their strong association with cotyledons during development. In addition, we identified genes GbLBD25a (lateral organ boundaries domain 25a), GbCESA2a (cellulose synthase 2a), GbMYB74c (myeloblastosis 74c), GbPIN2 (PIN-FORMED 2) associated with vascular development regulation, and GbWRKY1 (WRKYGOK 1), GbbHLH12a (basic helix-loop-helix 12a), and GbJAZ4 (jasmonate zim-domain 4) potentially involved in the formation of secretory cavities. Moreover, we found that flavonoid accumulation in mature embryos could enhance postgerminative growth and seedling establishment in harsh environments. Our 3D spatial reconstruction technique combined with multiomics analysis opens avenues for investigating developmental architecture and molecular mechanisms during embryogenesis and lays the foundation for evolutionary studies of embryo development and maturation.

Cross-regulation of cytoskeleton and calcium signaling at plant-pathogen interface

During plant-pathogen interactions, cytoskeleton and calcium signaling work independently as well as in coordination with each other for developing preformed and induced defense responses. A cell wall (CW) - plasma membrane (PM) - cytoskeleton (CS) continuum is maintained by coordination of cytoskeleton and calcium signaling. The current review is focused on the current knowledge of cytoskeleton‑calcium cross-regulation during plant-pathogen interactions. Implications of recent technological developments in the existing toolkit that can address the outstanding questions of cytoskeleton‑calcium coordination plant immunity are also discussed.

An overview of aptamer: Design strategy, prominent applications, and potential challenge in plants

Aptamers, serving as highly efficient molecular recognition and biotechnology tools, have garnered increasing interest in the realm of plant science in recent years. Aptamers are synthetic single-stranded short nucleotides or peptides, that bind targets with high specificity and affinity, triggering precise biological responses. As an alternative to antibodies, aptamers present promising avenues for advancement in biological researches. Aptamers function in a range of fields, encompassing cell signaling, drug development, biosensor technology, as well as botany, agricultural and forestry sciences. In this review, we introduce classifications and screening methods of aptamers, as well as aptamer-based technologies, highlighting their significant contributions to recent advancements. With their powerful functionality and ability to bind targets with high specificity and affinity, aptamers offer promising opportunities for breakthroughs in plant research.

A century journey of organelles research in the plant endomembrane system

We are entering an exciting century in the study of the plant organelles in the endomembrane system. Over the past century, especially within the past 50 years, tremendous advancements have been made in the complex plant cell to generate a much clearer and informative picture of plant organelles, including the molecular/morphological features, dynamic/spatial behavior, and physiological functions. Importantly, all these discoveries and achievements in the identification and characterization of organelles in the endomembrane system would not have been possible without: (1) the innovations and timely applications of various state-of-art cell biology tools and technologies for organelle biology research; (2) the continuous efforts in developing and characterizing new organelle markers by the plant biology community; and (3) the landmark studies on the identification and characterization of the elusive organelles. While molecular aspects and results for individual organelles have been extensively reviewed, the development of the techniques for organelle research in plant cell biology is less appreciated. As one of the ASPB Centennial Reviews on "organelle biology," here we aim to take a journey across a century of organelle biology research in plants by highlighting the important tools (or landmark technologies) and key scientists that contributed to visualize organelles. We then highlight the landmark studies leading to the identification and characterization of individual organelles in the plant endomembrane systems.

SCAB1 coordinates sequential Ca2+ and ABA signals during osmotic stress induced stomatal closure in Arabidopsis

Hyperosmotic stress caused by drought is a detrimental threat to plant growth and agricultural productivity due to limited water availability. Stomata are gateways of transpiration and gas exchange, the swift adjustment of stomatal aperture has a strong influence on plant drought resistance. Despite intensive investigations of stomatal closure during drought stress in past decades, little is known about how sequential signals are integrated during complete processes. Here, we discovered that the rapid Ca2+ signaling and subsequent abscisic acid (ABA) signaling contribute to the kinetics of both F-actin reorganizations and stomatal closure in Arabidopsis thaliana, while STOMATAL CLOSURE-RELATED ACTIN BINDING PROTEIN1 (SCAB1) is the molecular switch for this entire process. During the early stage of osmotic shock responses, swift elevated calcium signaling promotes SCAB1 phosphorylation through calcium sensors CALCIUM DEPENDENT PROTEIN KINASE3 (CPK3) and CPK6. The phosphorylation restrained the microfilament binding affinity of SCAB1, which bring about the F-actin disassembly and stomatal closure initiation. As the osmotic stress signal continued, both the kinase activity of CPK3 and the phosphorylation level of SCAB1 attenuated significantly. We further found that ABA signaling is indispensable for these attenuations, which presumably contributed to the actin filament reassembly process as well as completion of stomatal closure. Notably, the dynamic changes of SCAB1 phosphorylation status are crucial for the kinetics of stomatal closure. Taken together, our results support a model in which SCAB1 works as a molecular switch, and directs the microfilament rearrangement through integrating the sequentially generated Ca2+ and ABA signals during osmotic stress induced stomatal closure.

Membrane microdomains: Structural and signaling platforms for establishing membrane polarity

Cell polarity results from the asymmetric distribution of cellular structures, molecules, and functions. Polarity is a fundamental cellular trait that can determine the orientation of cell division, the formation of particular cell shapes, and ultimately the development of a multicellular body. To maintain the distinct asymmetric distribution of proteins and lipids in cellular membranes, plant cells have developed complex trafficking and regulatory mechanisms. Major advances have been made in our understanding of how membrane microdomains influence the asymmetric distribution of proteins and lipids. In this review, we first give an overview of cell polarity. Next, we discuss current knowledge concerning membrane microdomains and their roles as structural and signaling platforms to establish and maintain membrane polarity, with a special focus on the asymmetric distribution of proteins and lipids, and advanced microscopy techniques to observe and characterize membrane microdomains. Finally, we review recent advances regarding membrane trafficking in cell polarity establishment and how the balance between exocytosis and endocytosis affects membrane polarity.

VA-TIRFM-based SM kymograph analysis for dwell time and colocalization of plasma membrane protein in plant cells

BACKGROUND

The plasma membrane (PM) proteins function in a highly dynamic state, including protein trafficking and protein homeostasis, to regulate various biological processes. The dwell time and colocalization of PM proteins are considered to be two important dynamic features determining endocytosis and protein interactions, respectively. Dwell-time and colocalization detected using traditional fluorescence microscope techniques are often misestimated due to bulk measurement. In particular, analyzing these two features of PM proteins at the single-molecule level with spatiotemporal continuity in plant cells remains greatly challenging.

RESULTS

We developed a single molecular (SM) kymograph method, which is based on variable angle-total internal reflection fluorescence microscopy (VA-TIRFM) observation and single-particle (co-)tracking (SPT) analysis, to accurately analyze the dwell time and colocalization of PM proteins in a spatial and temporal manner. Furthermore, we selected two PM proteins with distinct dynamic behaviors, including AtRGS1 (Arabidopsis regulator of G protein signaling 1) and AtREM1.3 (Arabidopsis remorin 1.3), to analyze their dwell time and colocalization upon jasmonate (JA) treatment by SM kymography. First, we established new 3D (2D+t) images to view all trajectories of the interest protein by rotating these images, and then we chose the appropriate point without changing the trajectory for further analysis. Upon JA treatment, the path lines of AtRGS1-YFP appeared curved and short, while the horizontal lines of mCherry-AtREM1.3 demonstrated limited changes, indicating that JA might initiate the endocytosis of AtRGS1. Analysis of transgenic seedlings coexpressing AtRGS1-YFP/mCherry-AtREM1.3 revealed that JA induces a change in the trajectory of AtRGS1-YFP, which then merges into the kymography line of mCherry-AtREM1.3, implying that JA increases the colocalization degree between AtRGS1 and AtREM1.3 on the PM. These results illustrate that different types of PM proteins exhibit specific dynamic features in line with their corresponding functions.

CONCLUSIONS

The SM-kymograph method provides new insight into quantitively analyzing the dwell time and correlation degree of PM proteins at the single-molecule level in living plant cells.

Spatial organization and connectivity of wood rays in Pinus massoniana xylem based on high-resolution μCT-assisted network analysis

MAIN CONCLUSION

Spatial organization and connectivity of wood rays in Pinus massoniana was comprehensively viewed and regarded as anatomical adaptions to ensure the properties of rays in xylem. Spatial organization and connectivity of wood rays are essential for understanding the wood hierarchical architecture, but the spatial information is ambiguous due to small cell size. Herein, 3D visualization of rays in Pinus massoniana was performed using high-resolution μCT. We found brick-shaped rays were 6.5% in volume fractions, nearly twice the area fractions estimated by 2D levels. Uniseriate rays became taller and wider during the transition from earlywood to latewood, which was mainly contributed from the height increment of ray tracheids and widened ray parenchyma cells. Furthermore, both volume and surface area of ray parenchyma cells were larger than ray tracheids, so ray parenchyma took a higher proportion in rays. Moreover, three different types of pits for connectivity were segmented and revealed. Pits in both axial tracheids and ray tracheids were bordered, but the pit volume and pit aperture of earlywood axial tracheids were almost tenfold and over fourfold larger than ray tracheids. Contrarily, cross-field pits between ray parenchyma and axial tracheids were window-like with the principal axis of 31.0 μm, but its pit volume was approximately one-third of axial tracheids. Additionally, spatial organization of rays and axial resin canal was analyzed by a curved surface reformation tool, providing the first evidence of rays close to epithelial cells inward through the resin canal. Epithelial cells had various morphologies and large variations in cell size. Our results give new insights into the organization of radial system of xylem, especially the connectivity of rays with adjacent cells.

PagUNE12 encodes a basic helix-loop-helix transcription factor that regulates the development of secondary vascular tissue in poplar

Secondary growth in woody plants generates new cells and tissues via the activity of the vascular cambium and drives the radial expansion of stems and roots. It is regulated by a series of endogenous factors, especially transcription factors. Here, we cloned the basic helix-loop-helix (bHLH) transcription factor gene UNFERTILIZED EMBRYO SAC12 (UNE12) from poplar (Populus alba × Populus glandulosa Uyeki) and used biochemical, molecular, and cytological assays to investigate the biological functions and regulatory mechanism of PagUNE12. PagUNE12 mainly localized in the nucleus and possessed transcriptional activation activity. It was widely expressed in vascular tissues, including primary phloem and xylem and secondary phloem and xylem. Poplar plants overexpressing PagUNE12 showed significantly reduced plant height, shorter internodes, and curled leaves compared with wild-type plants. Optical microscopy and transmission electron microscopy revealed that overexpressing PagUNE12 promoted secondary xylem development, with thicker secondary cell walls than wild-type poplar. Fourier transform infrared spectroscopy, confocal Raman microscopy, and 2D Heteronuclear Single Quantum Correlation analysis indicated that these plants also had increased lignin contents, with a lower relative abundance of syringyl lignin units and a higher relative abundance of guaiacyl lignin units. Therefore, overexpressing PagUNE12 promoted secondary xylem development and increased the lignin contents of secondary xylem in poplar, suggesting that this gene could be used to improve wood quality in the future.

Endoplasmic reticulum membrane contact sites: cross-talk between membrane-bound organelles in plant cells

Eukaryotic cells contain organelles surrounded by monolayer or bilayer membranes. Organelles take part in highly dynamic and organized interactions at membrane contact sites, which play vital roles during development and response to stress. The endoplasmic reticulum extends throughout the cell and acts as an architectural scaffold to maintain the spatial distribution of other membrane-bound organelles. In this review, we highlight the structural organization, dynamics, and physiological functions of membrane contact sites between the endoplasmic reticulum and various membrane-bound organelles, especially recent advances in plants. We briefly introduce how the combined use of dynamic and static imaging techniques can enable monitoring of the cross-talk between organelles via membrane contact sites. Finally, we discuss future directions for research fields related to membrane contact.

Dense Sampling Approaches for Psychiatry Research: Combining Scanners and Smartphones

Together, data from brain scanners and smartphones have sufficient coverage of biology, psychology, and environment to articulate between-person differences in the interplay within and across biological, psychological, and environmental systems thought to underlie psychopathology. An important next step is to develop frameworks that combine these two modalities in ways that leverage their coverage across layers of human experience to have maximum impact on our understanding and treatment of psychopathology. We review literature published in the last 3 years highlighting how scanners and smartphones have been combined to date, outline and discuss the strengths and weaknesses of existing approaches, and sketch a network science framework heretofore underrepresented in work combining scanners and smartphones that can push forward our understanding of health and disease.

Membrane Dynamics Regulated by Cytoskeleton in Plant Immunity

The plasma membrane (PM), which is composed of a lipid layer implanted with proteins, has diverse functions in plant responses to environmental triggers. The heterogenous dynamics of lipids and proteins in the plasma membrane play important roles in regulating cellular activities with an intricate pathway that orchestrates reception, signal transduction and appropriate response in the plant immune system. In the process of the plasma membrane participating in defense responses, the cytoskeletal elements have important functions in a variety of ways, including regulation of protein and lipid dynamics as well as vesicle trafficking. In this review, we summarized how the plasma membrane contributed to plant immunity and focused on the dynamic process of cytoskeleton regulation of endocytosis and exocytosis and propose future research directions.

SMPD1 expression profile and mutation landscape help decipher genotype-phenotype association and precision diagnosis for acid sphingomyelinase deficiency

BACKGROUND

Acid sphingomyelinase deficiency (ASMD) disorder, also known as Niemann-Pick disease (NPD) is a rare genetic disease caused by mutations in SMPD1 gene, which encodes sphingomyelin phosphodiesterase (ASM). Except for liver and spleen enlargement and lung disease, two subtypes (Type A and B) of NDP have different onset times, survival times, ASM activities, and neurological abnormalities. To comprehensively explore NPD's genotype-phenotype association and pathophysiological characteristics, we collected 144 NPD cases with strict quality control through literature mining.

RESULTS

The difference in ASM activity can differentiate NPD type A from other subtypes, with the ratio of ASM activity to the reference values being lower in type A (threshold 0.045 (4.45%)). Severe variations, such as deletion and insertion, can cause complete loss of ASM function, leading to type A, whereas relatively mild missense mutations generally result in type B. Among reported mutations, the p.Arg3AlafsX76 mutation is highly prevalent in the Chinese population, and the p.R608del mutation is common in Mediterranean countries. The expression profiles of SMPD1 from GTEx and single-cell RNA sequencing data of multiple fetal tissues showed that high expressions of SMPD1 can be observed in the liver, spleen, and brain tissues of adults and hepatoblasts, hematopoietic stem cells, STC2_TLX1-positive cells, mesothelial cells of the spleen, vascular endothelial cells of the cerebellum and the cerebrum of fetuses, indicating that SMPD1 dysfunction is highly likely to have a significant effect on the function of those cell types during development and the clinicians need pay attention to these organs or tissues as well during diagnosis. In addition, we also predicted 21 new pathogenic mutations in the SMPD1 gene that potentially cause the NPD, signifying that more rare cases will be detected with those mutations in SMPD1. Finally, we also analysed the function of the NPD type A cells following the extracellular milieu.

CONCLUSIONS

Our study is the first to elucidate the effects of SMPD1 mutation on cell types and at the tissue level, which provides new insights into the genotype-phenotype association and can help in the precise diagnosis of NPD.

Multiscale microscopy to decipher plant cell structure and dynamics

New imaging methodologies with high contrast and molecular specificity allow researchers to analyze dynamic processes in plant cells at multiple scales, from single protein and RNA molecules to organelles and cells, to whole organs and tissues. These techniques produce informative images and quantitative data on molecular dynamics to address questions that cannot be answered by conventional biochemical assays. Here, we review selected microscopy techniques, focusing on their basic principles and applications in plant science, discussing the pros and cons of each technique, and introducing methods for quantitative analysis. This review thus provides guidance for plant scientists in selecting the most appropriate techniques to decipher structures and dynamic processes at different levels, from protein dynamics to morphogenesis.

The Cytoskeleton in Plant Immunity: Dynamics, Regulation, and Function

The plant cytoskeleton, consisting of actin filaments and microtubules, is a highly dynamic filamentous framework involved in plant growth, development, and stress responses. Recently, research has demonstrated that the plant cytoskeleton undergoes rapid remodeling upon sensing pathogen attacks, coordinating the formation of microdomain immune complexes, the dynamic and turnover of pattern-recognizing receptors (PRRs), the movement and aggregation of organelles, and the transportation of defense compounds, thus serving as an important platform for responding to pathogen infections. Meanwhile, pathogens produce effectors targeting the cytoskeleton to achieve pathogenicity. Recent findings have uncovered several cytoskeleton-associated proteins mediating cytoskeletal remodeling and defense signaling. Furthermore, the reorganization of the actin cytoskeleton is revealed to further feedback-regulate reactive oxygen species (ROS) production and trigger salicylic acid (SA) signaling, suggesting an extremely complex role of the cytoskeleton in plant immunity. Here, we describe recent advances in understanding the host cytoskeleton dynamics upon sensing pathogens and summarize the effectors that target the cytoskeleton. We highlight advances in the regulation of cytoskeletal remodeling associated with the defense response and assess the important function of the rearrangement of the cytoskeleton in the immune response. Finally, we propose suggestions for future research in this area.

Arabidopsis ORP2A mediates ER-autophagosomal membrane contact sites and regulates PI3P in plant autophagy

Autophagy is an intracellular degradation system for cytoplasmic constituents which is mediated by the formation of a double-membrane organelle termed the autophagosome and its subsequent fusion with the lysosome/vacuole. The formation of the autophagosome requires membrane from the endoplasmic reticulum (ER) and is tightly regulated by a series of autophagy-related (ATG) proteins and lipids. However, how the ER contacts autophagosomes and regulates autophagy remain elusive in plants. In this study, we identified and demonstrated the roles of Arabidopsis oxysterol-binding protein-related protein 2A (ORP2A) in mediating ER-autophagosomal membrane contacts and autophagosome biogenesis. We showed that ORP2A localizes to both ER-plasma membrane contact sites (EPCSs) and autophagosomes, and that ORP2A interacts with both the ER-localized VAMP-associated protein (VAP) 27-1 and ATG8e on the autophagosomes to mediate the membrane contact sites (MCSs). In ORP2A artificial microRNA knockdown (KD) plants, seedlings display retarded growth and impaired autophagy levels. Both ATG1a and ATG8e accumulated and associated with the ER membrane in ORP2A KD lines. Moreover, ORP2A binds multiple phospholipids and shows colocalization with phosphatidylinositol 3-phosphate (PI3P) in vivo. Taken together, ORP2A mediates ER-autophagosomal MCSs and regulates autophagy through PI3P redistribution.

ADP-ribosylation factor D1 modulates Golgi morphology, cell plate formation, and plant growth in Arabidopsis

ADP-ribosylation factor (ARF) family proteins, one type of small guanine-nucleotide-binding (G) proteins, play a central role in regulating vesicular traffic and organelle structures in eukaryotes. The Arabidopsis (Arabidopsis thaliana) genome contains more than 21 ARF proteins, but relatively little is known about the functional heterogeneity of ARF homologs in plants. Here, we characterized the function of a unique ARF protein, ARFD1B, in Arabidopsis. ARFD1B exhibited both cytosol and punctate localization patterns, colocalizing with a Golgi marker in protoplasts and transgenic plants. Distinct from other ARF1 homologs, overexpression of a dominant-negative mutant form of ARFD1B did not alter the localization of the Golgi marker mannosidase I (ManI)-RFP in Arabidopsis cells. Interestingly, the ARFD1 artificial microRNA knockdown mutant arfd1 displayed a deleterious growth phenotype, while this phenotype was restored in complemented plants. Further, confocal imaging and transmission electron microscopy analyses of the arfd1 mutant revealed defective cell plate formation and abnormal Golgi morphology. Pull-down and liquid chromatography-tandem mass spectrometry analyses identified Coat Protein I (COPI) components as interacting partners of ARFD1B, and subsequent bimolecular fluorescence complementation, yeast (Saccharomyces cerevisiae) two-hybrid, and co-immunoprecipitation assays further confirmed these interactions. These results demonstrate that ARFD1 is required for cell plate formation, maintenance of Golgi morphology, and plant growth in Arabidopsis.

Molecular studies of cellulose synthase supercomplex from cotton fiber reveal its unique biochemical properties

Cotton fiber is a highly elongated and thickened single cell that produces large quantities of cellulose, which is synthesized and assembled into cell wall microfibrils by the cellulose synthase complex (CSC). In this study, we report that in cotton (Gossypium hirsutum) fibers harvested during secondary cell wall (SCW) synthesis, GhCesA 4, 7, and 8 assembled into heteromers in a previously uncharacterized 36-mer-like cellulose synthase supercomplex (CSS). This super CSC was observed in samples prepared using cotton fiber cells harvested during the SCW synthesis period but not from cotton stem tissue or any samples obtained from Arabidopsis. Knock-out of any of GhCesA 4, 7, and 8 resulted in the disappearance of the CSS and the production of fiber cells with no SCW thickening. Cotton fiber CSS showed significantly higher enzyme activity than samples prepared from knock-out cotton lines. We found that the microfibrils from the SCW of wild-type cotton fibers may contain 72 glucan chains in a bundle, unlike other plant materials studied. GhCesA4, 7, and 8 restored both the dwarf and reduced vascular bundle phenotypes of their orthologous Arabidopsis mutants, potentially by reforming the CSC hexamers. Genetic complementation was not observed when non-orthologous CesA genes were used, indicating that each of the three subunits is indispensable for CSC formation and for full cellulose synthase function. Characterization of cotton CSS will increase our understanding of the regulation of SCW biosynthesis.

Cytology, transcriptomics, and mass spectrometry imaging reveal changes in late-maturation elm (Ulmus pumila) seeds

During seed maturation, the seed deposits storage compounds (starches, oils, and proteins), synthesizes defense compounds, produces a seed coat, initiates embryo dormancy, and becomes desiccated. During the late-maturation stage, seed storage compound contents and compositions change dramatically. Although maturation has been extensively studied in model species and crops, it remains less well characterized in woody perennial plants. In this study, we conducted morphological and cytological observations, transcriptome profiling, and chemical constituent analysis of elm (Ulmus pumila L.) seeds during the late-maturation stage. Light and electron microscopy revealed that closely packed yet discrete lipid bodies frequently surrounded the densely stained protein bodies, and the protein bodies became irregular or even partially disintegrated at the end of seed development. RNA-seq detected substantial transcriptome changes during the late-maturation stage, and pathway enrichment analysis showed that the differentially expressed genes were associated with phenylpropanoid biosynthesis, starch and sucrose metabolism, plant-pathogen interactions, and hormone signal transduction. Furthermore, we used mass spectrometry imaging to detect the relative intensity and spatial distribution of fatty acids, phospholipids, and waxes in elm seeds. Our findings provide a framework for understanding the changes in cytological features and chemical composition during the final stage of elm seed development, and a detailed reference for seed development in woody plants.