The Arabidopsis histone demethylase JMJ28 regulates CONSTANS by interacting with FBH transcription factors.

Arabidopsis thaliana CONSTANS (CO) is an essential transcription factor that promotes flowering by activating the expression of the floral integrator FLOWERING LOCUS T (FT). A number of histone modification enzymes involved in the regulation of flowering have been identified, but the involvement of epigenetic mechanisms in the regulation of the core flowering regulator CO remains unclear. Previous studies have indicated that the transcription factors, FLOWERING BHLH1 (FBH1), FBH2, FBH3, and FBH4, function redundantly to activate the expression of CO. In this study, we found that the KDM3 group H3K9 demethylase JMJ28 interacts with the FBH transcription factors to activate CO by removing the repressive mark H3K9me2. The occupancy of JMJ28 on the CO locus is decreased in the fbh quadruple mutant, suggesting that the binding of JMJ28 is dependent on FBHs. Furthermore, genome-wide occupancy profile analyses indicate that the binding of JMJ28 to the genome overlaps with that of FBH3, indicating a functional association of JMJ28 and FBH3. Together, these results indicate that Arabidopsis JMJ28 functions as a CO activator by interacting with the FBH transcription factors to remove H3K9me2 from the CO locus.

A nonsynonymous mutation in BhLS, encoding an acyl-CoA N-acyltransferase leads to fruit and seed size variation in wax gourd (Benincasa hispida).

Wax gourd (Benincasa hispida (Thunb.) Cogn., 2n = 2x = 24) is an economically important vegetable crop cultivated widely in many tropical and subtropical regions, including China, India, and Japan. Both fruit and seeds are prized agronomic attributes in wax gourd breeding and production. However, the genetic mechanisms underlying these traits remain largely unexplored. In this study, we observed a strong correlation between fruit size and seed size variation in our mapping population, indicating genetic control by a single gene, BhLS, with large size being dominant over small. Through bulk segregant analysis sequencing and fine mapping with a large F2 population, we precisely located the BhLS gene within a 47.098-kb physical interval on Chromosome 10. Within this interval, only one gene, Bhi10M000649, was identified, showing homology to Arabidopsis HOOKLESS1. A nonsynonymous mutation (G to C) in the second exon of Bhi10M000649 was found to be significantly associated with both fruit and seed size variation in wax gourd. These findings collectively highlight the pleiotropic effect of the BhLS gene in regulating fruit and seed size in wax gourd. Our results offer molecular insights into the variation of fruit and seed size in wax gourd and establish a fundamental framework for breeding wax gourd cultivars with desired traits.

MSI1 and HDA6 function interdependently to control flowering time via chromatin modifications.

MULTICOPY SUPPRESSOR OF IRA1 (MSI1) is a conserved subunit of Polycomb Repressive Complex 2 (PRC2), which mediates gene silencing by histone H3 lysine 27 trimethylation (H3K27Me3). Here, we demonstrated that MSI1 interacts with the RPD3-like histone deacetylase HDA6 both in vitro and in vivo. MSI1 and HDA6 are involved in flowering and repress the expression of FLC, MAF4, and MAF5 by removing H3K9 acetylation but adding H3K27Me3. Chromatin immunoprecipitation analysis showed that HDA6 and MSI1 interdependently bind to the chromatin of FLC, MAF4, and MAF5. Furthermore, H3K9 deacetylation mediated by HDA6 is dependent on MSI1, while H3K27Me3 mediated by PRC2 containing MSI1 is also dependent on HDA6. Taken together, these data indicate that MSI1 and HDA6 act interdependently to repress the expression of FLC, MAF4, and MAF5 through histone modifications. Our findings reveal that the HDA6-MSI1 module mediates the interaction between histone H3 deacetylation and H3K27Me3 to repress gene expression involved in flowering time control.

miR2105 and the kinase OsSAPK10 co-regulate OsbZIP86 to mediate drought-induced ABA biosynthesis in rice.

Mediating induced abscisic acid (ABA) biosynthesis is important for enhancing plant stress tolerance. Here, we found that rice (Oryza sativa L.) osa-miR2105 (miR2105) and the Stress/ABA-activated protein kinase (OsSAPK10) coordinately regulate the rice basic region-leucine zipper transcription factor (bZIP TF; OsbZIP86) at the posttranscriptional and posttranslational levels to control drought-induced ABA biosynthesis via modulation of rice 9-cis-epoxycarotenoid dioxygenase (OsNCED3) expression. OsbZIP86 expression is regulated by miR2105-directed cleavage of the OsbZIP86 mRNA. OsbZIP86 encodes a nuclear TF that binds to the promoter of the ABA biosynthetic gene OsNCED3. OsSAPK10 can phosphorylate and activate OsbZIP86 to enhance the expression of OsNCED3. Under normal growth conditions, altered expression of miR2105 and OsbZIP86 displayed no substantial effect on rice growth. However, under drought conditions, miR2105 knockdown or OsbZIP86 overexpression transgenic rice plants showed higher ABA content, enhanced tolerance to drought, lower rates of water loss, and more stomatal closure of seedlings, compared with wild-type rice Zhonghua 11; in contrast, miR2105 overexpression, OsbZIP86 downregulation, and OsbZIP86 knockout plants displayed opposite phenotypes. Collectively, our results show that the "miR2105-(OsSAPK10)-OsbZIP86-OsNCED3" module regulates the drought-induced ABA biosynthesis without penalty on rice growth under normal conditions, suggesting candidates for improving drought tolerance in rice.

The chromatin remodelling ATPase BRAHMA interacts with GATA-family transcription factor GNC to regulate flowering time in Arabidopsis.

BRAHMA (BRM) is the ATPase of the SWItch/Sucrose Non-Fermentable (SWI/SNF) chromatin remodelling complex, which is indispensable for transcriptional inhibition and activation, associated with vegetative and reproductive development in Arabidopsis thaliana. Here, we show that BRM directly binds to the chromatin of SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1), which integrates multiple flowering signals to regulate floral transition, leading to flowering. In addition, genetic and molecular analysis showed that BRM interacts with GNC (GATA, NITRATE-INDUCIBLE, CARBON METABOLISM INVOLVED), a GATA transcription factor that represses flowering by directly repressing SOC1 expression. Furthermore, BRM is recruited by GNC to directly bind to the chromatin of SOC1. The transcript level of SOC1 is elevated in brm-3, gnc, and brm-3/gnc mutants, which is associated with increased histone H3 lysine 4 tri-methylation (H3K4Me3) but decreased DNA methylation. Taken together, our results indicate that BRM associates with GNC to regulate SOC1 expression and flowering time.

Arabidopsis JMJ29 is involved in trichome development by regulating the core trichome initiation gene GLABRA3.

Plant trichomes are large single cells that are organized in a regular pattern and play multiple biological functions. In Arabidopsis, trichome development is mainly governed by the core trichome initiation regulators, including the R2R3 type MYB transcript factor GLABRA 1 (GL1), bHLH transcript factors GLABRA 3/ENHANCER OF GLABRA 3 (GL3/EGL3), and the WD-40 repeat protein TRANSPARENT TESTA GLABRA 1 (TTG1), as well as the downstream trichome regulator GLABRA 2 (GL2). GL1, GL3/EGL3, and TTG1 can form a trimeric activation complex to activate GL2, which is required for the trichome initiation and maintenance during cell differentiation. Arabidopsis JMJ29 is a JmjC domain-containing histone demethylase belonging to the JHDM2/KDM3 group. Members of the JHDM2/KDM3 group histone demethylases are mainly responsible for the H3K9me1/2 demethylation. In the present study, we found that the trichome density on leaves and inflorescence stems is significantly decreased in jmj29 mutants. The expression of the core trichome regulators GL1, GL2, and GL3 is decreased in jmj29 mutants as well. Furthermore, JMJ29 can directly target GL3 and remove H3K9me2 on the GL3 locus. Collectively, we found that Arabidopsis JMJ29 is involved in trichome development by directly regulating GL3 expression. These results provide further insights into the molecular mechanism of epigenetic regulation in Arabidopsis trichome development.

SWI3B and HDA6 interact and are required for transposon silencing in Arabidopsis.

Although the interplay of covalent histone acetylation/deacetylation and ATP-dependent chromatin remodelling is crucial for the regulation of chromatin structure and gene expression in eukaryotes, the underlying molecular mechanism in plants remains largely unclear. Here we show a direct interaction between Arabidopsis SWI3B, an essential subunit of the SWI/SNF chromatin-remodelling complex, and the RPD3/HDA1-type histone deacetylase HDA6 both in vitro and in vivo. Furthermore, SWI3B and HDA6 co-repress the transcription of a subset of transposons. Both SWI3B and HDA6 maintain transposon silencing by decreasing histone H3 lysine 9 acetylation, but increasing histone H3 lysine 9 di-methylation, DNA methylation and nucleosome occupancy. Our findings reveal that SWI3B and HDA6 may act in the same co-repressor complex to maintain transposon silencing in Arabidopsis.

HISTONE DEACETYLASE6 Acts in Concert with Histone Methyltransferases SUVH4, SUVH5, and SUVH6 to Regulate Transposon Silencing.

Histone deacetylases (HDACs) play important roles in regulating gene expression. In yeast and animals, HDACs act as components of multiprotein complexes that modulate transcription during various biological processes. However, little is known about the interacting proteins of plant HDACs. To identify the plant HDAC complexes and interacting proteins, we developed an optimized workflow using immunopurification coupled to mass spectrometry-based proteomics in Arabidopsis thaliana We found that the histone deacetylase HDA6 can interact with the histone methyltransferases SUVH4, SUVH5, and SUVH6 (SUVH4/5/6). Domain analysis revealed that the C-terminal regions of HDA6 and SUVH5 are important for their interaction. Furthermore, HDA6 interacts with SUVH4/5/6 and coregulates a subset of transposons through histone H3K9 methylation and H3 deacetylation. In addition, two phosphorylated serine residues, S427 and S429, were unambiguously identified in the C-terminal region of HDA6. Phosphomimetics (amino acid substitutions that mimic a phosphorylated protein) of HDA6 resulted in increased enzymatic activity, whereas the mutation of S427 to alanine in HDA6 abolished its interaction with SUVH5 and SUVH6, suggesting that the phosphorylation of HDA6 is important for its activity and function.

Identification and Characterization of Tomato SWI3-Like Proteins: Overexpression of SlSWIC Increases the Leaf Size in Transgenic Arabidopsis.

As the subunits of the SWI/SNF (mating-type switching (SWI) and sucrose nonfermenting (SNF)) chromatin-remodeling complexes (CRCs), Swi3-like proteins are crucial to chromatin remodeling in yeast and human. Growing evidence indicate that AtSWI3s are also essential for development and response to hormones in Arabidopsis. Nevertheless, the biological functions of Swi3-like proteins in tomato (Solanum lycopersicum) have not been investigated. Here we identified four Swi3-like proteins from tomato, namely SlSWI3A, SlSWI3B, SlSWI3C, and SlSWI3D. Subcellular localization analysis revealed that all SlSWI3s are localized in the nucleus. The expression patterns showed that all SlSWI3s are ubiquitously expressed in all tissues and organs, and SlSWI3A and SlSWI3B can be induced by cold treatment. In addition, we found that SlSWI3B can form homodimers with itself and heterodimers with SlSWI3A and SlSWI3C. SlSWI3B can also interact with SlRIN and SlCHR8, two proteins involved in tomato reproductive development. Overexpression of SlSWI3C increased the leaf size in transgenic Arabidopsis with increased expression of GROWTH REGULATING FACTORs, such as GRF3, GRF5, and GRF6. Taken together, our results indicate that SlSWI3s may play important roles in tomato growth and development.

Roles of the INO80 and SWR1 Chromatin Remodeling Complexes in Plants.

Eukaryotic genes are packed into a dynamic but stable nucleoprotein structure called chromatin. Chromatin-remodeling and modifying complexes generate a dynamic chromatin environment that ensures appropriate DNA processing and metabolism in various processes such as gene expression, as well as DNA replication, repair, and recombination. The INO80 and SWR1 chromatin remodeling complexes (INO80-c and SWR1-c) are ATP-dependent complexes that modulate the incorporation of the histone variant H2A.Z into nucleosomes, which is a critical step in eukaryotic gene regulation. Although SWR1-c has been identified in plants, plant INO80-c has not been successfully isolated and characterized. In this review, we will focus on the functions of the SWR1-c and putative INO80-c (SWR1/INO80-c) multi-subunits and multifunctional complexes in Arabidopsis thaliana. We will describe the subunit compositions of the SWR1/INO80-c and the recent findings from the standpoint of each subunit and discuss their involvement in regulating development and environmental responses in Arabidopsis.

A SUMO Ligase AtMMS21 Regulates the Stability of the Chromatin Remodeler BRAHMA in Root Development.

Chromatin remodeling is essential for gene expression regulation in plant development and response to stresses. Brahma (BRM) is a conserved ATPase in the SWI/SNF chromatin remodeling complex and is involved in various biological processes in plant cells, but the regulation mechanism on BRM protein remains unclear. Here, we report that BRM interacts with AtMMS21, a SUMO ligase in Arabidopsis (Arabidopsis thaliana). The interaction was confirmed in different approaches in vivo and in vitro. The mutants of BRM and AtMMS21 displayed a similar defect in root development. In the mms21-1 mutant, the protein level of BRM-GFP was significantly lower than that in wild type, but the RNA level of BRM did not change. Biochemical evidence indicated that BRM was modified by SUMO3, and the reaction was enhanced by AtMMS21. Furthermore, overexpression of wild-type AtMMS21 but not the mutated AtMMS21 without SUMO ligase activity was able to recover the stability of BRM in mms21-1 Overexpression of BRM in mms21-1 partially rescued the developmental defect of roots. Taken together, these results supported that AtMMS21 regulates the protein stability of BRM in root development.

The Arabidopsis SWI2/SNF2 Chromatin Remodeling ATPase BRAHMA Targets Directly to PINs and Is Required for Root Stem Cell Niche Maintenance.

BRAHMA (BRM), a SWI/SNF chromatin remodeling ATPase, is essential for the transcriptional reprogramming associated with development and cell differentiation in Arabidopsis thaliana. In this study, we show that loss-of-function mutations in BRM led to defective maintenance of the root stem cell niche, decreased meristematic activity, and stunted root growth. Mutations of BRM affected auxin distribution by reducing local expression of several PIN-FORMED (PIN) genes in the stem cells and impaired the expression of the stem cell transcription factor genes PLETHORA (PLT1) and PLT2. Chromatin immunoprecipitation assays showed that BRM could directly target to the chromatin of PIN1, PIN2, PIN3, PIN4, and PIN7. In addition, genetic interaction assays indicate that PLTs acted downstream of BRM, and overexpression of PLT2 partially rescued the stem cell niche defect of brm mutants. Taken together, these results support the idea that BRM acts in the PLT pathway to maintain the root stem cell niche by altering the expression of PINs.

Arabidopsis BREVIPEDICELLUS interacts with the SWI2/SNF2 chromatin remodeling ATPase BRAHMA to regulate KNAT2 and KNAT6 expression in control of inflorescence architecture.

BREVIPEDICELLUS (BP or KNAT1), a class-I KNOTTED1-like homeobox (KNOX) transcription factor in Arabidopsis thaliana, contributes to shaping the normal inflorescence architecture through negatively regulating other two class-I KNOX genes, KNAT2 and KNAT6. However, the molecular mechanism of BP-mediated transcription regulation remains unclear. In this study, we showed that BP directly interacts with the SWI2/SNF2 chromatin remodeling ATPase BRAHMA (BRM) both in vitro and in vivo. Loss-of-function BRM mutants displayed inflorescence architecture defects, with clustered inflorescences, horizontally orientated pedicels, and short pedicels and internodes, a phenotype similar to the bp mutants. Furthermore, the transcript levels of KNAT2 and KNAT6 were elevated in brm-3, bp-9 and brm-3 bp-9 double mutants. Increased histone H3 lysine 4 tri-methylation (H3K4me3) levels were detected in brm-3, bp-9 and brm-3 bp-9 double mutants. Moreover, BRM and BP co-target to KNAT2 and KNAT6 genes, and BP is required for the binding of BRM to KNAT2 and KNAT6. Taken together, our results indicate that BP interacts with the chromatin remodeling factor BRM to regulate the expression of KNAT2 and KNAT6 in control of inflorescence architecture.

The Arabidopsis SWI2/SNF2 chromatin Remodeler BRAHMA regulates polycomb function during vegetative development and directly activates the flowering repressor gene SVP.

The chromatin remodeler BRAHMA (BRM) is a Trithorax Group (TrxG) protein that antagonizes the functions of Polycomb Group (PcG) proteins in fly and mammals. Recent studies also implicate such a role for Arabidopsis (Arabidopsis thaliana) BRM but the molecular mechanisms underlying the antagonism are unclear. To understand the interplay between BRM and PcG during plant development, we performed a genome-wide analysis of trimethylated histone H3 lysine 27 (H3K27me3) in brm mutant seedlings by chromatin immunoprecipitation followed by next generation sequencing (ChIP-seq). Increased H3K27me3 deposition at several hundred genes was observed in brm mutants and this increase was partially supressed by removal of the H3K27 methyltransferase CURLY LEAF (CLF) or SWINGER (SWN). ChIP experiments demonstrated that BRM directly binds to a subset of the genes and prevents the inappropriate association and/or activity of PcG proteins at these loci. Together, these results indicate a crucial role of BRM in restricting the inappropriate activity of PcG during plant development. The key flowering repressor gene SHORT VEGETATIVE PHASE (SVP) is such a BRM target. In brm mutants, elevated PcG occupancy at SVP accompanies a dramatic increase in H3K27me3 levels at this locus and a concomitant reduction of SVP expression. Further, our gain- and loss-of-function genetic evidence establishes that BRM controls flowering time by directly activating SVP expression. This work reveals a genome-wide functional interplay between BRM and PcG and provides new insights into the impacts of these proteins in plant growth and development.

Regulation of flowering time by the histone deacetylase HDA5 in Arabidopsis.

The acetylation level of histones on lysine residues regulated by histone acetyltransferases and histone deacetylases plays an important but under-studied role in the control of gene expression in plants. With the aim of characterizing the Arabidopsis RPD3/HDA1 family histone deacetylase HDA5, we present evidence showing that HDA5 displays deacetylase activity. Mutants defective in the expression of HDA5 displayed a late-flowering phenotype. Expression of the flowering repressor genes FLC and MAF1 was up-regulated in hda5 mutants. Furthermore, the gene activation markers, histone H3 acetylation and H3K4 trimethylation on FLC and MAF1 chromatin were increased in hda5-1 mutants. Chromatin immunoprecipitation analysis showed that HDA5 binds to the chromatin of FLC and MAF1. Bimolecular fluorescence complementation assays and co-immunoprecipitation assays showed that HDA5 interacts with FVE, FLD and HDA6, indicating that these proteins are present in a protein complex involved in the regulation of flowering time. Comparing gene expression profiles of hda5 and hda6 mutants by RNA-seq revealed that HDA5 and HDA6 co-regulate gene expression in multiple development processes and pathways.

PHYTOCHROME INTERACTING FACTOR3 associates with the histone deacetylase HDA15 in repression of chlorophyll biosynthesis and photosynthesis in etiolated Arabidopsis seedlings.

PHYTOCHROME INTERACTING FACTOR3 (PIF3) is a key basic helix-loop-helix transcription factor of Arabidopsis thaliana that negatively regulates light responses, repressing chlorophyll biosynthesis, photosynthesis, and photomorphogenesis in the dark. However, the mechanism for the PIF3-mediated transcription regulation remains largely unknown. In this study, we found that the REDUCED POTASSIUM DEPENDENCY3/HISTONE DEACETYLASE1-type histone deacetylase HDA15 directly interacted with PIF3 in vivo and in vitro. Genome-wide transcriptome analysis revealed that HDA15 acts mainly as a transcriptional repressor and negatively regulates chlorophyll biosynthesis and photosynthesis gene expression in etiolated seedlings. HDA15 and PIF3 cotarget to the genes involved in chlorophyll biosynthesis and photosynthesis in the dark and repress gene expression by decreasing the acetylation levels and RNA Polymerase II-associated transcription. The binding of HDA15 to the target genes depends on the presence of PIF3. In addition, PIF3 and HDA15 are dissociated from the target genes upon exposure to red light. Taken together, our results indicate that PIF3 associates with HDA15 to repress chlorophyll biosynthetic and photosynthetic genes in etiolated seedlings.

HISTONE DEACETYLASE19 interacts with HSL1 and participates in the repression of seed maturation genes in Arabidopsis seedlings.

The seed maturation genes are specifically and highly expressed during late embryogenesis. In this work, yeast two-hybrid, bimolecular fluorescence complementation, and coimmunoprecipitation assays revealed that HISTONE DEACETYLASE19 (HDA19) interacted with the HIGH-LEVEL EXPRESSION OF SUGAR-INDUCIBLE GENE2-LIKE1 (HSL1), and the zinc-finger CW [conserved Cys (C) and Trp (W) residues] domain of HSL1 was responsible for the interaction. Furthermore, we found that mutations in HDA19 resulted in the ectopic expression of seed maturation genes in seedlings, which was associated with increased levels of gene activation marks, such as Histone H3 acetylation (H3ac), Histone H4 acetylation (H4ac), and Histone H3 Lys 4 tri-methylation (H3K4me3), but decreased levels of the gene repression mark Histone H3 Lys 27 tri-methylation (H3K27me3) in the promoter and/or coding regions. In addition, elevated transcription of certain seed maturation genes was also found in the hsl1 mutant seedlings, which was also accompanied by the enrichment of gene activation marks but decreased levels of the gene repression mark. Chromatin immunoprecipitation assays showed that HDA19 could directly bind to the chromatin of the seed maturation genes. These results suggest that HDA19 and HSL1 may act together to repress seed maturation gene expression during germination. Further genetic analyses revealed that the homozygous hsl1 hda19 double mutants are embryonic lethal, suggesting that HDA19 and HSL1 may play a vital role during embryogenesis.

Chromosome-level genome assembly and population genomics reveals crucial selection for subgynoecy development in chieh-qua.

Chieh-qua is an important cucurbit crop and very popular in South China and Southeast Asia. Despite its significance, its genetic basis and domestication history are unclear. In this study, we have successfully generated a chromosome-level reference genome assembly for the chieh-qua 'A36' using a hybrid assembly strategy that combines PacBio long reads and Illumina short reads. The assembled genome of chieh-qua is approximately 953.3 Mb in size and is organized into 12 chromosomes, with contig N50 of 6.9 Mb and scaffold N50 of 68.2 Mb. Notably, the chieh-qua genome is comparable in size to the wax gourd genome. Through gene prediction analysis, we have identified a total of 24 593 protein-coding genes in the A36 genome. Additionally, approximately 56.6% (539.3 Mb) of the chieh-qua genome consists of repetitive sequences. Comparative genome analysis revealed that chieh-qua and wax gourd are closely related, indicating a close evolutionary relationship between the two species. Population genomic analysis, employing 129 chieh-qua accessions and 146 wax gourd accessions, demonstrated that chieh-qua exhibits greater genetic diversity compared to wax gourd. We also employed the GWAS method to identify related QTLs associated with subgynoecy, an interested and important trait in chieh-qua. The MYB59 (BhiCQ0880026447) exhibited relatively high expression levels in the shoot apex of four subgynoecious varieties compared with monoecious varieties. Overall, this research provides insights into the domestication history of chieh-qua and offers valuable genomic resources for further molecular research.

Involvement of histone modifications in plant abiotic stress responses.

As sessile organisms, plants encounter various environmental stimuli including abiotic stresses during their lifecycle. To survive under adverse conditions, plants have evolved intricate mechanisms to perceive external signals and respond accordingly. Responses to various stresses largely depend on the plant capacity to modulate the transcriptome rapidly and specifically. A number of studies have shown that the molecular mechanisms driving the responses of plants to environmental stresses often depend on nucleosome histone post-translational modifications including histone acetylation, methylation, ubiquitination, and phosphorylation. The combined effects of these modifications play an essential role in the regulation of stress responsive gene expression. In this review, we highlight our current understanding of the epigenetic mechanisms of histone modifications and their roles in plant abiotic stress response.

Responses of differential metabolites and pathways to high temperature in cucumber anther.

Cucumber is one of the most important vegetable crops, which is widely planted all over the world. Cucumber always suffers from high-temperature stress in South China in summer. In this study, liquid chromatography-mass spectrometry (LC-MS) analysis was used to study the differential metabolites of cucumber anther between high-temperature (HT) stress and normal condition (CK). After HT, the pollen fertility was significantly reduced, and abnormal anther structures were observed by the paraffin section. In addition, the metabolomics analysis results showed that a total of 125 differential metabolites were identified after HT, consisting of 99 significantly upregulated and 26 significantly downregulated metabolites. Among these differential metabolites, a total of 26 related metabolic pathways were found, and four pathways showed significant differences, namely, porphyrin and chlorophyll metabolism; plant hormone signal transduction; amino sugar and nucleotide sugar metabolism; and glycine, serine, and threonine metabolism. In addition, pollen fertility was decreased by altering the metabolites of plant hormone signal transduction and amino acid and sugar metabolism pathway under HT. These results provide a comprehensive understanding of the metabolic changes in cucumber anther under HT.