Plant regeneration in the new era: from molecular mechanisms to biotechnology applications.

Plants or tissues can be regenerated through various pathways. Like animal regeneration, cell totipotency and pluripotency are the molecular basis of plant regeneration. Detailed systematic studies on Arabidopsis thaliana gradually unravel the fundamental mechanisms and principles underlying plant regeneration. Specifically, plant hormones, cell division, epigenetic remodeling, and transcription factors play crucial roles in reprogramming somatic cells and reestablishing meristematic cells. Recent research on basal non-vascular plants and monocot crops has revealed that plant regeneration differs among species, with various plant species using distinct mechanisms and displaying significant differences in regenerative capacity. Conducting multi-omics studies at the single-cell level, tracking plant regeneration processes in real-time, and deciphering the natural variation in regenerative capacity will ultimately help understand the essence of plant regeneration, improve crop regeneration efficiency, and contribute to future crop design.

Dynamic root microbiome sustains soybean productivity under unbalanced fertilization.

Root-associated microbiomes contribute to plant growth and health, and are dynamically affected by plant development and changes in the soil environment. However, how different fertilizer regimes affect quantitative changes in microbial assembly to effect plant growth remains obscure. Here, we explore the temporal dynamics of the root-associated bacteria of soybean using quantitative microbiome profiling (QMP) to examine its response to unbalanced fertilizer treatments (i.e., lacking either N, P or K) and its role in sustaining plant growth after four decades of unbalanced fertilization. We show that the root-associated bacteria exhibit strong succession during plant development, and bacterial loads largely increase at later stages, particularly for Bacteroidetes. Unbalanced fertilization has a significant effect on the assembly of the soybean rhizosphere bacteria, and in the absence of N fertilizer the bacterial community diverges from that of fertilized plants, while lacking P fertilizer impedes the total load and turnover of rhizosphere bacteria. Importantly, a SynCom derived from the low-nitrogen-enriched cluster is capable of stimulating plant growth, corresponding with the stabilized soybean productivity in the absence of N fertilizer. These findings provide new insights in the quantitative dynamics of the root-associated microbiome and highlight a key ecological cluster with prospects for sustainable agricultural management.

Clinical PDT dose dosimetry for pleural Photofrin-mediated photodynamic therapy.

Significance: Photodynamic therapy (PDT) is an established cancer treatment utilizing light-activated photosensitizers (PS). Effective treatment hinges on the PDT dose-dependent on PS concentration and light fluence-delivered over time. We introduce an innovative eight-channel PDT dose dosimetry system capable of concurrently measuring light fluence and PS concentration during treatment. Aim: We aim to develop and evaluate an eight-channel PDT dose dosimetry system for simultaneous measurement of light fluence and PS concentration. By addressing uncertainties due to tissue variations, the system enhances accurate PDT dosimetry for improved treatment outcomes. Approach: The study positions eight isotropic detectors strategically within the pleural cavity before PDT. These detectors are linked to bifurcated fibers, distributing signals to both a photodiode and a spectrometer. Calibration techniques are applied to counter tissue-related variations and improve measurement accuracy. The fluorescence signal is normalized using the measured light fluence, compensating for variations in tissue properties. Measurements were taken in 78 sites in the pleural cavities of 20 patients. Results: Observations reveal minimal Photofrin concentration variation during PDT at each site, juxtaposed with significant intra- and inter-patient heterogeneities. Across 78 treated sites in 20 patients, the average Photofrin concentration for all 78 sites is 4.98 µM, with a median concentration of 4.47 µM. The average PDT dose for all 78 sites is 493.17 µMJ/cm2, with a median dose of 442.79 µMJ/cm2. A significant variation in PDT doses is observed, with a maximum difference of 3.1 times among all sites within one patient and a maximum difference of 9.8 times across all patients. Conclusions: The introduced eight-channel PDT dose dosimetry system serves as a valuable real-time monitoring tool for light fluence and PS concentration during PDT. Its ability to mitigate uncertainties arising from tissue properties enhances dosimetry accuracy, thus optimizing treatment outcomes and bolstering the effectiveness of PDT in cancer therapy.

E3 ubiquitin ligase TaSDIR1-4A activates membrane-bound transcription factor TaWRKY29 to positively regulate drought resistance.

Drought is a deleterious abiotic stress factor that constrains crop growth and development. Post-translational modification of proteins mediated by the ubiquitin-proteasome system is an effective strategy for directing plant responses to stress, but the regulatory mechanisms in wheat remain unclear. In this study, we showed that TaSDIR1-4A is a positive modulator of the drought response. Overexpression of TaSDIR1-4A increased the hypersensitivity of stomata, root length and endogenous abscisic acid (ABA) content under drought conditions. TaSDIR1-4A encodes a C3H2C3-type RING finger protein with E3 ligase activity. Amino acid mutation in its conserved domain led to loss of activity and altered the subcellular localization. The membrane-bound transcription factor TaWRKY29 was identified by yeast two-hybrid screening, and it was confirmed as interacting with TaSDIR1-4A both in vivo and in vitro. TaSDIR1-4A mediated the polyubiquitination and proteolysis of the C-terminal amino acid of TaWRKY29, and its translocation from the plasma membrane to the nucleus. Activated TaWRKY29 bound to the TaABI5 promoter to stimulate its expression, thereby positively regulating the ABA signalling pathway and drought response. Our findings demonstrate the positive role of TaSDIR1-4A in drought tolerance and provide new insights into the involvement of UPS in the wheat stress response.

Fractionated Photofrin-Mediated Photodynamic Therapy Significantly Improves Long-Term Survival.

This study investigates the effect of fractionated (two-part) PDT on the long-term local control rate (LCR) using the concentration of reactive oxygen species ([ROS]rx) as a dosimetry quantity. Groups with different fractionation schemes are examined, including a 2 h interval between light delivery sessions to cumulative fluences of 135, 180, and 225 J/cm2. While the total treatment time remains constant within each group, the division of treatment time between the first and second fractionations are explored to assess the impact on long-term survival at 90 days. In all preclinical studies, Photofrin is intravenously administered to mice at a concentration of 5 mg/kg, with an incubation period between 18 and 24 h before the first light delivery session. Fluence rate is fixed at 75 mW/cm2. Treatment ensues via a collimated laser beam, 1 cm in diameter, emitting light at 630 nm. Dosimetric quantities are assessed for all groups along with long-term (90 days) treatment outcomes. This study demonstrated a significant improvement in long-term survival after fractionated treatment schemes compared to single-fraction treatment, with the optimal 90-day survival increasing to 63%, 86%, and 100% vs. 20%, 25%, and 50%, respectively, for the three cumulative fluences. The threshold [ROS]rx for the optimal scheme of fractionated Photofrin-mediated PDT, set at 0.78 mM, is significantly lower than that for the single-fraction PDT, at 1.08 mM.

GPRC5A regulates proliferation and oxidative stress by inhibiting the STAT3/Socs3/c-MYC pathway in hepatocellular carcinoma.

The G protein-coupled receptor, class C, group 5, member A (GPRC5A) plays a key role in various diseases, but its effect on hepatocellular carcinoma (HCC) and the potential underlying mechanisms remains unclear. In the present study, we explored the effect of GPRC5A on the progression of HCC and further explored its mechanism of action. The results revealed that the expression of GPRC5A was lower in HCC tissues and cells. Overexpression of GPRC5A suppressed the proliferation and epithelial-mesenchymal transition (EMT) of HCC cells. In addition, overexpression of GPRC5A induced oxidative stress and apoptosis. Further study showed that overexpression of GPRC5A inhibited the expression of STAT3/Socs3/c-MYC related-protein and the NLRP3 inflammasome. Moreover, the STAT3/Socs3/c-MYC and NLRP3 inflammasome was involved in the effect of GPRC5A on HCC cells. These results suggest that GPRC5A suppresses proliferation and EMT, induces oxidative stress and leads to apoptosis of HCC cells, potentially by regulating STAT3/Socs3/c-MYC signalling and the NLRP3 inflammasome. These findings suggest that GPRC5A has an anti-tumor effect in the formation of HCC, and the molecular therapy of GPRC5A provides a theoretical basis for treating HCC.

Targeting CXCL9/10/11–CXCR3 axis: an important component of tumor-promoting and antitumor immunity

Chemokines are chemotactic-competent molecules composed of a family of small cytokines, playing a key role in regulating tumor progression. The roles of chemokines in antitumor immune responses are of great interest. CXCL9, CXCL10, and CXCL11 are important members of chemokines. It has been widely investigated that these three chemokines can bind to their common receptor CXCR3 and regulate the differentiation, migration, and tumor infiltration of immune cells, directly or indirectly affecting tumor growth and metastasis. Here, we summarize the mechanism of how the CXCL9/10/11–CXCR3 axis affects the tumor microenvironment, and list the latest researches to find out how this axis predicts the prognosis of different cancers. In addition, immunotherapy improves the survival of tumor patients, but some patients show drug resistance. Studies have found that the regulation of CXCL9/10/11–CXCR3 on the tumor microenvironment is involved in the process of changing immunotherapy resistance. Here we also describe new approaches to restoring sensitivity to immune checkpoint inhibitors through the CXCL9/10/11–CXCR3 axis (AU)

Superior gluten structure and more small starch granules synergistically confer dough quality for high amylose wheat varieties.

High amylose wheat (HAW) has potential health benefits but its dough structure is usually inferior. Wheat dough is a complex mixture and its structure is influenced by the physicochemical properties of gluten and starch. In this study, we investigated the starch granule development, gluten structure, starch properties, pasting, and thermal properties of flour, as well as the rheological properties of dough in wheat variety Xinong 836 with high amylose content (33.57%) and its parents. The results showed that Xinong 836 wheat starch contained more small starch granules, which was consistent with the microstructural results of starch granules in grain filling stage. Moreover, Xinong 836 wheat starch showed highest swelling power and water solubility. Importantly, the flour of Xinong 836 wheat had the highest protein content and wet gluten content and Xinong 836 wheat gluten showed highest ß-sheets content and disulfide bond content than its parents Zhengmai 7698 and Xinong 979, which conferring to more compact microscopic networks of dough, thereby contributing to the higher peak viscosity (PV), final viscosity (FV), and setback viscosity (SB) in the flour of Xinong 836. Our finding elucidated that the stability of gluten and properties of starch synergistically affected the pasting and thermal properties of the flour paste, and the presence of more small starch granules contributed to dough with a rather dense structure in HAW Xinong 836. Thus, superior gluten structure and more small starch granules have synergistic effects on enhancing the gluten-starch interaction, thereby contributing to better dough quality.

Genome-wide analysis and identification of TaRING-H2 gene family and TaSDIR1 positively regulates salt stress tolerance in wheat.

Salt stress is an abiotic stress factor that limits high yields, and thus identifying salt tolerance genes is very important for improving the tolerance of salt in wheat. In this study we identified 274 TaRING-H2 family members and analyzed their gene positions, gene structures, conserved structural domains, promoter cis-acting elements and covariance relationships. And we investigated TaRING-H2-120 (TaSDIR1) in salt stress. Transgenic lines exhibited higher salt tolerance in the germination and seedling stages. Compared with the wild type, overexpression of TaSDIR1 upregulated the expression of genes encoding enzymes related to the control of reactive oxygen species (ROS), thereby reducing the accumulation of ROS, as well as increased the expression of ion transport-related genes to limit the inward flow of Na+ in vivo and maintain a higher K+/Na+ ratio. The expression levels of these genes were opposite in lines where TaSDIR1 was silenced by BSMV-VIGS, and the silenced wheat exhibited higher salt sensitivity. Arabidopsis mutants and heterologous TaSDIR1 overexpressing lines had similar salt stress tolerance phenotypes. We also demonstrated that TaSDIR1 interacted with TaSDIR1P2 in vivo and in vitro. A sequence of 80-100 amino acids in TaSDIR1P2 encoded a coiled coil domain that was important for the activity of E3 ubiquitin ligase, and it was also the core region for the interaction between TaSDIR1 and TaSDIR1P2. Overall, our results suggest that TaSDIR1 positively regulates salt stress tolerance in wheat.

Genome-wide analysis and identification of light-harvesting chlorophyll a/b binding (LHC) gene family and BSMV-VIGS silencing TaLHC86 reduced salt tolerance in wheat.

The discovery and identification of gene families by using wide-genome and public databases is an effective way to gain initial insight into gene function, which also is one of the current hot spots of research. Chlorophyll ab-binding proteins (LHC) are important for photosynthesis and widely involved in plant adversity stress. However, the study in wheat has not been reported. In this study, we identified 127 TaLHC members from common wheat which were unevenly distributed on all chromosomes except 3B and 3D. All members divided into three subfamilies, LHC a, LHC b and the LHC t which was only discovered in wheat. All of them had maximum expression in leaves and contained multiple light-responsive cis-acting element, which were evidence of the extensive involvement of LHC families in photosynthesis. In addition, we also analyzed their collinear relationship, targeting relationship with miRNA and their responses under different stresses. Based on these analyses, it was found that TaLHC86 was an excellent candidate gene for stress resistance. The full-length ORF of TaLHC86 was 792 bp and was localized on the chloroplasts. The salt tolerance of wheat was reduced when BSMV-VIGS silenced TaLHC86, and the photosynthetic rate and electron transport were also seriously affected. This study made a comprehensive analysis of the TaLHC family and found that TaLHC86 was a good gene for salt tolerance.

A pathogen effector FOLD diversified in symbiotic fungi.

Pathogenic fungi use secreted effector proteins to suppress immunity and support their infection, but effectors have also been reported from fungi that engage in nutritional symbioses with plants. Sequence-based effector comparisons between pathogens and symbiotic arbuscular mycorrhizal (AM) fungi are hampered by the huge diversity of effector sequences even within closely related microbes. To find sequence-divergent but structurally similar effectors shared between symbiotic and pathogenic fungi, we compared secreted protein structure models of the AM fungus Rhizophagus irregularis to known pathogen effectors. We identified proteins with structural similarity to known Fusarium oxysporum f. sp. lycopersici dual domain (FOLD) effectors, which occur in low numbers in several fungal pathogens. Contrastingly, FOLD genes from AM fungi (MycFOLDs) are found in enlarged and diversified gene families with higher levels of positive selection in their C-terminal domains. Our structure model comparison suggests that MycFOLDs are similar to carbohydrate-binding motifs. Different MycFOLD genes are expressed during colonisation of different hosts and MycFOLD-17 transcripts accumulate in plant intracellular arbuscules. The exclusive presence of MycFOLDs across unrelated plant-colonising fungi, their inducible expression, lineage-specific sequence diversification and transcripts in arbuscules suggest that FOLD proteins act as effectors during plant colonisation of symbiotic and pathogenic fungi.

Lutein inhibits tumor progression through the ATR/Chk1/p53 signaling pathway in non-small cell lung cancer.

Lung cancer is the leading cause of cancer-related death. In particular, non-small cell lung cancer (NSCLC) accounts for approximately 85% of all lung cancer cases. Due to tumor resistance and the toxicity of chemotherapeutic agents, it is increasingly critical to discover novel, potent antitumorigenic drugs for treating NSCLC. Lutein, a carotenoid, has been reported to exert toxic effects on cells in several tumor types. However, the detailed functions and underlying mechanisms of lutein in NSCLC remain elusive. The present study showed that lutein significantly and dose-dependently inhibited cell proliferation, arrested the cell cycle at the G0/G1 phase, and induced apoptosis in NSCLC cells. RNA-sequencing analysis revealed that the p53 signaling pathway was the most significantly upregulated in lutein-treated A549 cells. Mechanistically, lutein exerted antitumorigenic effects by inducing DNA damage and subsequently activating the ATR/Chk1/p53 signaling pathway in A549 cells. In vivo, lutein impeded tumor growth in mice and prolonged their survival. In conclusion, our findings demonstrate the antitumorigenic potential of lutein and reveal its molecular mechanism of action, suggesting that lutein is a promising candidate for clinical NSCLC treatment.

Targeting CXCL9/10/11-CXCR3 axis: an important component of tumor-promoting and antitumor immunity.

Chemokines are chemotactic-competent molecules composed of a family of small cytokines, playing a key role in regulating tumor progression. The roles of chemokines in antitumor immune responses are of great interest. CXCL9, CXCL10, and CXCL11 are important members of chemokines. It has been widely investigated that these three chemokines can bind to their common receptor CXCR3 and regulate the differentiation, migration, and tumor infiltration of immune cells, directly or indirectly affecting tumor growth and metastasis. Here, we summarize the mechanism of how the CXCL9/10/11-CXCR3 axis affects the tumor microenvironment, and list the latest researches to find out how this axis predicts the prognosis of different cancers. In addition, immunotherapy improves the survival of tumor patients, but some patients show drug resistance. Studies have found that the regulation of CXCL9/10/11-CXCR3 on the tumor microenvironment is involved in the process of changing immunotherapy resistance. Here we also describe new approaches to restoring sensitivity to immune checkpoint inhibitors through the CXCL9/10/11-CXCR3 axis.

Multispectral singlet oxygen luminescent dosimetry (MSOLD) for Photofrin-mediated Photodynamic Therapy.

Direct detection of singlet-state oxygen ([1O2]) constitutes the holy grail dosimetric method for type II PDT, a goal that can be quantified using multispectral singlet oxygen dosimetry (MSOLD). However, the short lifetime and extremely weak nature of the singlet oxygen signal produced has given rise to a need to improve MSOLD signal-to-noise ratio. This study examines methods for optimizing MSOLD signal acquisition, specifically employing an orthogonal arrangement between detection and PDT treatment light, consisting of two fiber optics - connected to a 632-nm laser and an InGaAs detector respectively. Light collected by the InGaAs detector is then passed through a filter wheel, where spectral emission measurements are taken at 1200 nm, 1240 nm, 1250 nm, 1270 nm, and 1300 nm. The data, after fitting to the fluorescence background and a gaussian-fit for the singlet oxygen peak, is established for the background-subtracted singlet oxygen emission signal. The MSOLD signal is then compared with the singlet oxygen explicit dosimetry (SOED) results, based on direct measurements of in-vivo light fluence (rate), in-vivo Photofrin concentration, and tissue oxygenation concentration. This study focuses on validating the sensitivity and minimum detectability of MSOLD signal in various in-vitro conditions. Finally, the MSOLD device will be tested in Photofrin-mediated PDT for mice bearing Radiation-Induced Fibrosarcoma (RIF) tumors.

Transcriptional regulation mechanism of wheat varieties with different nitrogen use efficiencies in response to nitrogen deficiency stress.

BACKGROUND: As one of the microelements, nitrogen play essential roles in cereal production. Although the use of chemical fertilizers has significantly improved the yield of wheat, it has also caused increasingly adverse environmental pollution. Revealing the molecular mechanism manipulating wheat nitrogen use efficiency (NUE), and cultivating wheat germplasms with high nitrogen use efficiency has become important goals for wheat researchers. In this study, we investigated the physiological and transcriptional differences of three wheat cultivars with different NUE under low nitrogen stress. RESULTS: The results showed that, under low nitrogen conditions, the activities of nitrogen metabolism-related enzymes (GS, NR, GDH), antioxidant enzymes (SOD, POD, CAT) and soluble protein contents of ZM366 (high NUE cultivar) were higher than those of JD8 (low NUE cultivar). The hybrid cultivar of ZM366 and JD8 showed mid-parent or over-parent heterosis. Transcriptome analysis revealed that 'alanine, aspartate and glutamate metabolism', 'terpenoid backbone biosynthesis' and 'vitamin B6 metabolism' pathways play key roles in nitrogen use efficiency in wheat. The significant enhancement of the 'Calvin cycle' and 'photorespiration' in ZM366 contributed to its higher level of carbon metabolism under low nitrogen stress, which is an important attribute differs from the other two varieties. In addition, the activation of ABA signal transduction and biosynthesis pathways also helps to maintain NUE under low- nitrogen conditions. Moreover, bHLH transcription factors were also found to play a positive role in wheat NUE. CONCLUSIONS: In conclusion, these results enriched our knowledge of the mechanism of wheat NUE, and provided a theoretical basis for improving wheat NUE and breeding new cultivars.

Golgi-localized putative S-adenosyl methionine transporters required for plant cell wall polysaccharide methylation.

Polysaccharide methylation, especially that of pectin, is a common and important feature of land plant cell walls. Polysaccharide methylation takes place in the Golgi apparatus and therefore relies on the import of S-adenosyl methionine (SAM) from the cytosol into the Golgi. However, so far, no Golgi SAM transporter has been identified in plants. Here we studied major facilitator superfamily members in Arabidopsis that we identified as putative Golgi SAM transporters (GoSAMTs). Knockout of the two most highly expressed GoSAMTs led to a strong reduction in Golgi-synthesized polysaccharide methylation. Furthermore, solid-state NMR experiments revealed that reduced methylation changed cell wall polysaccharide conformations, interactions and mobilities. Notably, NMR revealed the existence of pectin 'egg-box' structures in intact cell walls and showed that their formation is enhanced by reduced methyl esterification. These changes in wall architecture were linked to substantial growth and developmental phenotypes. In particular, anisotropic growth was strongly impaired in the double mutant. The identification of putative transporters involved in import of SAM into the Golgi lumen in plants provides new insights into the paramount importance of polysaccharide methylation for plant cell wall structure and function.

TaNBR1, a Novel Wheat NBR1-like Domain Gene Negatively Regulates Drought Stress Tolerance in Transgenic Arabidopsis.

Drought stress is an important factor that severely affects crop yield and quality. Autophagy has a crucial role in the responses to abiotic stresses. In this study, we explore TaNBR1 in response to drought stress. Expression of the TaNBR1 gene was strongly induced by NaCl, PEG, and abscisic acid treatments. The TaNBR1 protein is localized in the Golgi apparatus and autophagosome. Transgenic Arabidopsis plants overexpressing TaNBR1 exhibited reduced drought tolerance. When subjected to drought stress, compared to the wild-type (WT) lines, the transgenic overexpressing TaNBR1 plants had a lower seed germination rate, relative water content, proline content, and reduced accumulation of antioxidant enzymes, i.e., superoxide dismutase, peroxidase, and catalase, as well as higher chlorophyll losses, malondialdehyde contents, and water loss. The transgenic plants overexpressing TaNBR1 produced much shorter roots in response to mannitol stress, in comparison to the WT plants, and they exhibited greater sensitivity to abscisic acid treatment. The expression levels of the genes related to stress in the transgenic plants were affected in response to drought stress. Our results indicate that TaNBR1 negatively regulates drought stress responses by affecting the expression of stress-related genes in Arabidopsis.

Membrane-bound transcription factor TaNTL1 positively regulates drought stress tolerance in transgenic Arabidopsis.

Drought negatively affects plant growth and development to cause major yield losses in crops. Transcription factors (TFs) play important roles in abiotic stress response signaling in plant. However, the biological functions of membrane-bound transcription factors (MTFs) in abiotic stress have rarely been studied in wheat. In this study, we identified a homologue of the maize ZmNTL1 gene in wheat, which was designated as TaNTL1. TaNTL1 is a NAC family MTF (NTM1-like, NTL proteins) encoding 481 amino acid residues with a transmembrane motif at the C-terminal. Quantitative results and expression profile analysis showed that TaNTL1 could respond to drought. We demonstrated the transcriptional activity of TaNTL1 and that it could specifically bind to NAC recognition cis-acting elements (NACBS). The full-length TaNTL1 protein localized in the plasma membrane and TaNTL1 lacking the transmembrane motif (TaNTL1-ΔTM) localized in the nucleus. TaNTL1 was proteolytically activated by PEG6000 and abscisic acid (ABA). Phenotypic and physiological analyses showed that overexpression transgenic Arabidopsis exhibited enhanced drought resistance, which was greater with TaNTL1-ΔTM than TaNTL1. Transient silencing of TaNTL1 significantly reduced the resistance to drought stress in wheat. Germination by the TaNTL1 and TaNTL1-ΔTM transgenic Arabidopsis seeds was also hypersensitive to ABA. Most of the stress-related genes in transgenic plants were upregulated under drought conditions. These results suggest that MTF TaNTL1 is a positive regulator of drought and it may function by entering the nucleus through cleavage.

Comparative transcriptome and DNA methylation analysis in temperature-sensitive genic male sterile wheat BS366.

BACKGROUND: Known as the prerequisite component for the heterosis breeding system, the male sterile line determines the hybrid yield and seed purity. Therefore, a deep understanding of the mechanism and gene network that leads to male sterility is crucial. BS366, a temperature-sensitive genic male sterile (TGMS) line, is male sterile under cold conditions (12 °C with 12 h of daylight) but fertile under normal temperature (20 °C with 12 h of daylight). RESULTS: During meiosis, BS366 was defective in forming tetrads and dyads due to the abnormal cell plate. During pollen development, unusual vacuolated pollen that could not accumulate starch grains at the binucleate stage was also observed. Transcriptome analysis revealed that genes involved in the meiotic process, such as sister chromatid segregation and microtubule-based movement, were repressed, while genes involved in DNA and histone methylation were induced in BS366 under cold conditions. MethylRAD was used for reduced DNA methylation sequencing of BS366 spikes under both cold and control conditions. The differentially methylated sites (DMSs) located in the gene region were mainly involved in carbohydrate and fatty acid metabolism, lipid metabolism, and transport. Differentially expressed and methylated genes were mainly involved in cell division. CONCLUSIONS: These results indicated that the methylation of genes involved in carbon metabolism or fatty acid metabolism might contribute to male sterility in BS366 spikes, providing novel insight into the molecular mechanism of wheat male sterility.