The utilization of metabolomics approaches to explore the metabolic mechanisms underlying plant fitness and adaptation to dynamic environments is growing, highlighting the need for an efficient and user-friendly toolkit tailored for analyzing the extensive datasets generated by metabolomics studies. Current protocols for metabolome data analysis often struggle with handling large-scale datasets or require programming skills. To address this, we present MetMiner (https://github.com/ShawnWx2019/MetMiner), a user-friendly, full-functionality pipeline specifically designed for plant metabolomics data analysis. Built on R shiny, MetMiner can be deployed on servers to utilize additional computational resources for processing large-scale datasets. MetMiner ensures transparency, traceability, and reproducibility throughout the analytical process. Its intuitive interface provides robust data interaction and graphical capabilities, enabling users without prior programming skills to engage deeply in data analysis. Additionally, we constructed and integrated a plant-specific mass spectrometry database into the MetMiner pipeline to optimize metabolite annotation. We have also developed MDAtoolkits, which include a complete set of tools for statistical analysis, metabolite classification, and enrichment analysis, to facilitate the mining of biological meaning from the datasets. Moreover, we propose an iterative weighted gene co-expression network analysis strategy for efficient biomarker metabolite screening in large-scale metabolomics data mining. In two case studies, we validated MetMiner's efficiency in data mining and robustness in metabolite annotation. Together, the MetMiner pipeline represents a promising solution for plant metabolomics analysis, providing a valuable tool for the scientific community to use with ease.

Temperature sensitivity and tolerance play a key role in plant survival and production. Perennial ryegrass (Lolium perenne L.), widely cultivated in cool-season for forage supply and turfgrass, is extremely susceptible to high temperatures, therefore serving as an excellent grass for dissecting the genomic and genetic basis of high-temperature adaptation. In this study, expression analysis revealed that LpHsfA2, an important gene associated with high-temperature tolerance in perennial ryegrass, is rapidly and substantially induced under heat stress. Additionally, heat-tolerant varieties consistently display elevated expression levels of LpHsfA2 compared with heat-sensitive ones. Comparative haplotype analysis of the LpHsfA2 promoter indicated an uneven distribution of two haplotypes (HsfA2^Hap1 and HsfA2^Hap2) across varieties with differing heat tolerance. Specifically, the HsfA2^Hap1 allele is predominantly present in heat-tolerant varieties, while the HsfA2^Hap2 allele exhibits the opposite pattern. Overexpression of LpHsfA2 confers enhanced thermotolerance, whereas silencing of LpHsfA2 compromises heat tolerance. Furthermore, LpHsfA2 orchestrates its protective effects by directly binding to the promoters of LpHSP18.2 and LpAPX1 to activate their expression, preventing the non-specific misfolding of intracellular protein and the accumulation of reactive oxygen species in cells. Additionally, LpHsfA4 and LpHsfA5 were shown to engage directly with the promoter of LpHsfA2, upregulating its expression as well as the expression of LpHSP18.2 and LpAPX1, thus contributing to enhanced heat tolerance. Markedly, LpHsfA2 possesses autoregulatory ability by directly binding to its own promoter to modulate the self-transcription. Based on these findings, we propose a model for modulating the thermotolerance of perennial ryegrass by precisely regulating the expression of LpHsfA2. Collectively, these findings provide a scientific basis for the development of thermotolerant perennial ryegrass cultivars.

It has been proposed that cortical fine actin filaments are needed for the morphogenesis of pavement cells (PCs). However, the precise role and regulation mechanisms of actin filaments in PC morphogenesis are not well understood. Here, we found that Arabidopsis thaliana ACTIN DEPOLYMERIZING FACTOR9 (ADF9) is required for the morphogenesis of PC, which is negatively regulated by the R2R3 MYELOBLASTOSIS (MYB) transcription factor MYB52. In adf9 mutants, the lobe number of cotyledon PCs was significantly reduced, while the average lobe length did not differ significantly compared to that of wild type (Col-0), except for the variations in cell area and circularity, whereas the PC shapes in ADF9 overexpression seedlings showed different results. ADF9 decorated actin filaments, and colocalized with plasma membrane. The extent of filament bundling and actin filament bundling activity in adf9 mutant decreased. In addition, MYB52 directly targeted the promoter of ADF9 and negatively regulated its expression. The myb52-2 mutant showed increased lobe number and cell area, reduced cell circularity of PCs, and the PC phenotypes were suppressed when ADF9 was knocked out. Taken together, our data demonstrate that actin filaments play an important role in the morphogenesis of PC and reveal a transcriptional mechanism underlying MYB52 regulation of ADF9-mediated actin filament bundling in PC morphogenesis.

Leaves play a crucial role in the growth and development of rice (Oryza sativa) as sites for the production of photosynthesis. Early leaf senescence leads to substantial drops in rice yields. Whether and how DNA methylation regulates gene expression and affects leaf senescence remains elusive. Here, we demonstrate that mutations in rice ARGONAUTE 2 (OsAGO2) lead to premature leaf senescence, with chloroplasts in Osago2 having lower chlorophyll content and an abnormal thylakoid structure compared with those from wild-type plants. We show that OsAGO2 associates with a 24-nt microRNA and binds to the promoter region of OsNAC300, which causes DNA methylation and suppressed expression of OsNAC300. Overexpressing OsNAC300 causes the similar premature leaf senescence as Osago2 mutants and knocking out OsNAC300 in the Osago2 mutant background suppresses the early senescence of Osago2 mutants. Based on yeast one-hybrid, dual-luciferase, and electrophoresis mobility shift assays, we propose that OsNAC300 directly regulates transcription of the key rice aging gene NAC-like, activated by APETALA3/PISTILLATA (OsNAP) to control leaf senescence. Our results unravel a previously unknown epigenetic regulatory mechanism underlying leaf senescence in which OsAGO2–OsNAC300–OsNAP acts as a key regulatory module of leaf senescence to maintain leaf function.

Cotton fiber length is basically determined by well-coordinated gene expression and phosphatidylinositol phosphates (PIPs) accumulation during fiber elongation but the regulatory mechanism governing PIPs transport remains unknown. Here, we report a MYB transcription factor GhMYB30D04 in Gossypium hirsutum that promotes fiber elongation through modulating the expression of PIP transporter gene GhLTPG1. Knockout of GhMYB30D04 gene in cotton (KO) results in a reduction of GhLTPG1 transcripts with lower accumulation of PIPs, leading to shorter fibers and lower fiber yield. Conversely, GhMYB30D04 overexpression (GhMYB30D04-OE) causes richer PIPs and longer cotton fibers, mimicking the effects of exogenously applying PIPs on the ovules of GhMYB30D04-KO and wild type. Furthermore, GhMYB30D04 interacts with GhHD1, the crucial transcription factor of fiber initiation, to form an activation complex stabilized by PIPs, both of which upregulate GhLTPG1 expression. Comparative omics-analysis revealed that higher and extended expressions of LTPG1 in fiber elongation mainly correlate with the variations of the GhMYB30D04 gene between two cotton allotetraploids, contributing to longer fiber in G. babardense. Our work clarifies a mechanism by which GhHD1–GhMYB30D04 form a regulatory module of fiber elongation to tightly control PIP accumulation. Our work still has an implication that GhMYB30D04–GhHD1 associates with development transition from fiber initiation to elongation.

Anchorene, identified as an endogenous bioactive carotenoid-derived dialdehyde and diapocarotenoid, affects root development by modulating auxin homeostasis. However, the precise interaction between anchorene and auxin, as well as the mechanisms by which anchorene modulates auxin levels, remain largely elusive. In this study, we conducted a comparative analysis of anchorene's bioactivities alongside auxin and observed that anchorene induces multifaceted auxin-like effects. Through genetic and pharmacological examinations, we revealed that anchorene's auxin-like activities depend on the indole-3-pyruvate-dependent auxin biosynthesis pathway, as well as the auxin inactivation pathway mediated by Group II Gretchen Hagen 3 (GH3) proteins that mainly facilitate the conjugation of indole-3-acetic acid (IAA) to amino acids, leading to the formation of inactivated storage forms. Our measurements indicated that anchorene treatment elevates IAA levels while reducing the quantities of inactivated IAA–amino acid conjugates and oxIAA. RNA sequencing further revealed that anchorene triggers the expression of numerous auxin-responsive genes in a manner reliant on Group II GH3s. Additionally, our in vitro enzymatic assays and biolayer interferometry (BLI) assay demonstrated anchorene's robust suppression of GH3.17-mediated IAA conjugation with glutamate. Collectively, our findings highlight the significant role of carotenoid-derived metabolite anchorene in modulating auxin homeostasis, primarily through the repression of GH3-mediated IAA conjugation and inactivation pathways, offering novel insights into the regulatory mechanisms of plant bioactive apocarotenoids.

Eleocharis vivipara, an amphibious sedge in the Cyperaceae family, has several remarkable properties, most notably its alternate use of C₃ photosynthesis underwater and C₄ photosynthesis on land. However, the absence of genomic data has hindered its utility for evolutionary and genetic research. Here, we present a high-quality genome for E. vivipara, representing the first chromosome-level genome for the Eleocharis genus, with an approximate size of 965.22 Mb mainly distributed across 10 chromosomes. Its Hi–C pattern, chromosome clustering results, and one-to-one genome synteny across two subgroups indicates a tetraploid structure with chromosome count 2n = 4x = 20. Phylogenetic analysis suggests that E. vivipara diverged from Cyperus esculentus approximately 32.96 million years ago (Mya), and underwent a whole-genome duplication (WGD) about 3.5 Mya. Numerous fusion and fission events were identified between the chromosomes of E. vivipara and its close relatives. We demonstrate that E. vivipara has holocentromeres, a chromosomal feature which can maintain the stability of such chromosomal rearrangements. Experimental transplantation and cross-section studies showed its terrestrial culms developed C₄ Kranz anatomy with increased number of chloroplasts in the bundle sheath (BS) cells. Gene expression and weighted gene co-expression network analysis (WGCNA) showed overall elevated expression of core genes associated with the C₄ pathway, and significant enrichment of genes related to modified culm anatomy and photosynthesis efficiency. We found evidence of mixed nicotinamide adenine dinucleotide - malic enzyme and phosphoenolpyruvate carboxykinase type C₄ photosynthesis in E. vivipara, and hypothesize that the evolution of C₄ photosynthesis predates the WGD event. The mixed type is dominated by subgenome A and supplemented by subgenome B. Collectively, our findings not only shed light on the evolution of E. vivipara and karyotype within the Cyperaceae family, but also provide valuable insights into the transition between C₃ and C₄ photosynthesis, offering promising avenues for crop improvement and breeding.

During their co-evolution with herbivorous insects, plants have developed multiple defense strategies that resist pests, such as releasing a blend of herbivory-induced plant volatiles (HIPVs) that repel pests or recruit their natural enemies. However, the responses of insects to HIPVs in maize (Zea mays L.) are not well understood. Here, we demonstrate that the Asian corn borer (ACB, Ostrinia furnacalis), a major insect pest of maize, shows a preference for maize pre-infested with ACB larvae rather than being repelled by these plants. Through combined transcriptomic and metabolomics analysis of ACB-infested maize seedlings, we identified two substances that explain this behavior: (E)-4,8-dimethylnona-1,3,7-triene (DMNT) and (3E,7E)-4,8,12-trimethyltrideca-1,3,7,11-tetraene (TMTT). DMNT and TMTT attracted ACB larvae, and knocking out the maize genes responsible for their biosynthesis via gene editing impaired this attraction. External supplementation with DMNT/TMTT hampered the larvae's ability to locate pre-infested maize. These findings uncover a novel role for DMNT and TMTT in driving the behavior of ACB. Genetic modification of maize to make it less detectable by ACB might be an effective strategy for developing maize germplasm resistant to ACB and for managing this pest effectively in the field.

Disease resistance is often associated with compromised plant growth and yield due to defense-growth tradeoffs. However, key components and mechanisms underlying the defense-growth tradeoffs are rarely explored in maize. In this study, we find that ZmSKI3, a putative subunit of the SUPERKILLER (SKI) complex that mediates the 3′-5′ degradation of RNA, regulates both plant development and disease resistance in maize. The Zmski3 mutants showed retarded plant growth and constitutively activated defense responses, while the ZmSKI3 overexpression lines are more susceptible to Curvularia lunata and Bipolaris maydis. Consistently, the expression of defense-related genes was generally up-regulated, while expressions of growth-related genes were mostly down-regulated in leaves of the Zmski3-1 mutant compared to that of wild type. In addition, 223 differentially expressed genes that are up-regulated in Zmski3-1 mutant but down-regulated in the ZmSKI3 overexpression line are identified as potential target genes of ZmSKI3. Moreover, small interfering RNAs targeting the transcripts of the defense- and growth-related genes are differentially accumulated, likely to combat the increase of defense-related transcripts but decrease of growth-related transcripts in Zmski3-1 mutant. Taken together, our study indicates that plant growth and immunity could be regulated by both ZmSKI3-mediated RNA decay and post-transcriptional gene silencing in maize.