The common wheat Chinese cultivar Sumai3 which carries Fhb1 and has consistently shown a major effect on resistance to FHB (Li et al., 2019; Su et al., 2019) was collected and maintained by Henan Institute of Science and Technology (HIST), Xinxiang China. FHB susceptible common wheat Jimai22 was collected from the Shandong Academy of Agricultural Sciences, and preserved by HIST. Bainong207 was developed and maintained by HIST. Sumai3, Jimai22, and Bainong207 were used for gene expression analysis and leaves from three different individuals of the same treatment were collected and mixed at each time point. Bainong207 was also used for the BSMV-VIGS assay. Tissue expression of TaSHMTs in root, stem, and leaf were examined at the adult stage of Sumai3. For the abiotic stress treatments, 14-day-old wheat seedlings of Bainong207 were treated with 20% PEG6000, 200 mmol NaCl, 4°C, 100 μmol methyl jasmonate (MeJA), 100 μmol H₂O₂, and 100 μmol ABA. Leaf samples were collected after 1, 12, and 24 h treatment. Bainong207 was used for expression analysis and grown in a climatic chamber at 23°C/18°C, with a 14 h light/10 h dark cycle, and 70% relative humidity.

Mixed races of Blumeria graminis f. sp. tritici (Bgt) were collected in fields in Xinxiang, China and preserved on seedlings of the high susceptible variety Sumai3 in the climatic chamber. For RNA extraction, total RNAs were extracted from seedling leaves of Bainong207 inoculated with Bgt at 0, 6, 24, 48, and 72 hpi (hours post-inoculation) using Trizol reagent (Vazyme, Nanjing, China) following the manufacturer’s protocol. The F. graminearum used in this study was a field isolate originating in Henan, China and preserved by HIST. The spikelets of Sumai 3 and Jimai 22 were inoculated at the early flowering stage with 20 μL fresh spores of F. graminearum in the middle of the heads by the single-floret inoculation method. The spore concentration was 10¹ conidia mL^–1. Three spikelets (from three inoculated spikes of different individuals) were collected at 0, 24, 36, 48, and 72 hpi. All these three materials of Bainong207, Sumai3, and Jimai22 were grown in a greenhouse at 23°C/18°C, with a 14 h light/10 h dark cycle, and 70% relative humidity.

The qRT-PCR procedure was performed as described by Hu et al. (2018). The first-strand cDNA was synthesized using the HiScript Q RT SuperMix for qRT-PCR Kit (Vazyme, Nanjing, China). The qRT-PCR was performed using the SYBR Green detection kit AceQ qPCR SYBR Green Master Mix (Vazyme, Nanjing, China) on LC 480II (Roche, German). The program used was as follows: 5 min at 95°C, followed by 40 cycles at 95°C for 10 s and 60°C for 20 s. The relative gene expression was calculated by comparative 2^–ΔΔCT method. The primers used in this study are listed in Supplementary Table 1 and the wheat TaTubulin gene was used as an internal control.

The genome-wide data for T. aestivum (Chinese Spring) from IWGSC² was downloaded (IWGSC, 2018). Data for Triticum urartu (Tu 2.0) were downloaded from the MBKBase database³ (Ling et al., 2018). Triticum dicoccoides (WEWSeq_v.1.0), Aegilops tauschii (Aet_v4.0), Hordeum vulgare (IBSC_v2), Arabidopsis thaliana (TAIR10), Solanum lycopersicum (SL3.0), Cucumis sativus (ASM407v2), and Vitis vinifera (12X) were downloaded from the Ensemble Plants⁴ to construct a local database. The typical SHMT domain (PF00464) was downloaded from the Pfam database as the search model⁵ (El-Gebali et al., 2019). As described in Xu et al. (2021), a new hidden Markov model (HMM) was established to ensure the reliability of search results. The high-quality protein set was obtained using the raw SHMT HMM (E-value < 1 × 10 ^–20 and manual validation of an intact SHMT domain), and then used to build a specific SHMT HMM by the hmmbuild from the HMMERv3 suite (Lozano et al., 2015). The specific HMM was used, and proteins with an E-value lower than 0.001 were retained. The longest transcript was included for the following analysis when a gene has multiple transcripts. Both The SMART (Simple Modular Architecture Research Tool)⁵ (Letunic et al., 2021) and Conserved Domains⁶ (Lu et al., 2020) were used to check the candidate SHMT protein sequences again. The proteins containing complete SHMT conserved domains were reserved for further analysis, and named sequentially to their species and location on the chromosomes; all gene names are listed in Supplementary Table 2.

Multiple sequence alignment of all these full-length SHMT proteins was performed with ClustalW using the default options in MEGA-X (Kumar et al., 2018). Phylogenetic trees were constructed using the Maximum likelihood method of MEGAX with 1000 bootstrap replicates (Kumar et al., 2018; Yu et al., 2020). The phylogenetic tree was visualized with EvolView⁷ (He et al., 2016). The exon-intron structure was generated using TBtools based on the full-length genome sequence and the corresponding coding sequences (Chen C. et al., 2020). Conserved motifs analysis was performed using the MEME program⁸. The parameters were as follows: the maximum number of motifs was set to 20 and the optimum width was 6–50 residues (Xie et al., 2018; Xu et al., 2021). The gene structure with motif composition was visualized by the TBtools (Chen C. et al., 2020).

A micro-collinearity analysis is of great value for understanding gene evolutionary history and TGT (Triticeae-Gene Tribe⁹) was used to trace the origin history of the target gene, and gene pair was also analyzed with TGT (Chen Y. et al., 2020). DIVERGE v3.0 software was used to analyze the functional diversification among the subgroups based on the selected protein sequences (Gu et al., 2013).

BSMV-VIGS was performed as described by Wang et al. (2010) with some modifications. The fragment of TaSHMT3A-1 with the length of 249 bp was amplified with primer pair VIGS-TaSHMT3A-1 (Supplementary Table 1), and the target fragment was inserted into the γ-strain of BSMV to construct BSMV: TaSHMT3A-1 vector. An in vitro transcription kit (mMESSAGEmMACHINE T7, Invitrogen, Waltham, MA, United States) was used to produce the virus RNA. The common wheat Bainong207 was used for the gene silencing assay. When the second leaves were fully extended, the leaves infected with the virus BSMV: TaSHMT3A-1, with BSMV:TaPDS- and BSMV:γ- infected leaves as controls. The fourth leaves fully expanded with clear virus symptoms were used for disease resistance evaluation, and the inoculation method was performed as described in Xiang et al. (2011) with some modifications. The fourth leaves were detached from Bainong207, wounded on the adaxial surface, and then placed on the PCR board to form arch bridges. For F. graminearum inoculation, the inoculum comprised 1.5 μL of conidial suspension with a concentration at 1 × 10⁶ conidia mL^–1. The conidial suspension of F. graminearum was applied to the fresh wound, and then the PCR board with leaves was placed in 25 mgL^–1 benzimidazole water, at the same time ensuring that the cut of the leaves contacted the water, and then sealed for moisturizing. The inoculated leaves were cultured in an incubator with 14 h light/8 h dark at 23°C. The lesion size was observed 3 and 5 days after inoculation. Target gene silencing efficiency was evaluated by qRT-PCR using the primer pair TaSHMT3A-1-Q (Supplementary Table 1).

In this study, a total of 64 genes with SHMT conserved domain (PF00464) were identified from 10 sequenced plant species, 37 from monocotyledon, including T. aestivum (12), T. urartu (3), Ae. Tauschii (4), T. dicoccoides (9), H. vulgare (4) and O. sativa (5), and 27 from dicotyledons, including Arabidopsis (7), S. lycopersicum (7), C. sativus (6), and V. vinifera (7). To investigate the evolutionary relationships of the SHMT proteins, all the above 64 proteins were used to construct a Maximum likelihood phylogenetic tree (Figure 1). According to the classification of Arabidopsis and soybean SHMT (Lakhssassi et al., 2019), SHMT is divided into two main classes (classes I and II) and four subclasses (classes I a–b, II a, and II b). Among them, Class II a subgroup only contains proteins from dicotyledons, and the other subgroups include proteins from monocotyledons and dicotyledons. Within the same subgroup, the SHMT members of monocotyledons and dicotyledons were clustered together respectively, indicating that the SHMT genes of monocotyledons and dicotyledons had experienced great differentiation in the process of evolution. SHMT members within the same subgroup of Triticeae species had high protein sequence similarity, and the evolutionary process in Triticeae species was relatively conservative.

The numbers of SHMT in each subgroup of all analyzed species are listed in Table 1. Compared with monocotyledons, dicotyledons have more SHMT genes, mainly reflected in class I b, class II a, and class II b subgroups (Table 1). In Triticeae species, T. dicoccoides and T. aestivum had 9 and 12 SHMTs, which was about two and three times that of diploid species, respectively; moreover, the number of SHMTs in each subgroup was almost 2–3 times that of the diploid species (Table 1). This indicated that the increased number of TaSHMTs in polyploid wheat was primarily due to genome polyploidization. Interestingly, there was no SHMT gene of T. urartu in class I b, and the SHMT genes of barley and Ae. Tauschii were distributed on chromosomes 1, 2, 3, and 4, while T. urartu lacked the SHMT gene on chromosome 3; moreover, the SHMT genes on chromosome 3 of Ae. tauschii and barley were classed in class I b (Figure 1 and Table 2). This may be due to either a poor reference genome sequence of T. urartu or gene loss events that occurred during the evolution of T. urartu. The proportion of SHMT gene and copy number in class I a of T. dicoccoides was higher than that of other Triticeae species and there was an SHMT gene, TdSHMT7B-1, on chromosome 7B, but there was no SHMT gene on chromosome 7 of other Triticeae species. From the perspective of evolutionary relationships, the number and evolution of SHMT genes have differentiated between monocotyledons and dicotyledons, and within Triticeae species. Therefore, gene structure and conserved motifs of SHMT in different species were further analyzed to help explore the evolutionary process.

To further understand SHMT functional divergence, the conserved motifs of these SHMT proteins were identified using MEME software. Twenty individual motifs were identified (Figure 2). Our results show that, among all the analyzed 60 SHMT, proteins contained motif 4, motif 7 (except for TdSHMT7B-1), motif 10, and motif 18. All the analyzed SHMTs in Class I contained motif 9 (except for TdSHMT7B-1) and motif 17 (except for TdSHMT7B-1 and CsSHMT4-2); however, motif 9 and motif 17 lacked in the protein sequence of class II a (Figures 2A,B). Exon–intron structure divergence plays an important role during the evolution of duplicate genes and different compositions of the motif are important for their functional diversity (Yu et al., 2020; Xu et al., 2021). The intron-exon structure was analyzed by the aligning the full-length cDNA with genomic DNA sequence (Figure 2). The structure analysis of SHMT genes indicated that closely related members had a similar exon-intron structure. In class I, the number of exons ranged from three to five; among the 38 analyzed SHMT genes, 35 genes had four exons, and only TdSHMT7B-1 and AtSHMT4 in class I a had three exons, TaSHMT3B-1 in class I b had five exons (Figures 2A,C). Compared with class I, there were great differences in the intron–exon structure between class I and class II, each SHMT gene in class II had multiple short exons. Class II a, which is only composed of dicotyledon genes, contained 11 exons. Class II b contained 17 analyzed genes in total, and the number of exons from 13 to 16, especially, TdSHMT4A-1, had 13 exons. TdSHMT4B-1 and AetSHMT4D-1 had 14 exons, CsSHMT2-2 had 16 exons, and all the other 13 genes had 15 exons (76.47%). These data suggest that the exon-intron structure and motifs in SHMT were highly correlated with phylogenetic relationships and special motifs in SHMT may play critical roles in specific functions.

The micro-collinearity analysis helps to understand the replication or loss events of specific genes in the process of evolution or domestication; provides an opportunity to explore the inheritance and variation of genes in a local region, and can trace the origin history of a gene (Chen Y. et al., 2020). TGT was used to trace the origin history of the target gene under the guidance of the best matched collinear region. To explore the origin of SHMT genes in Triticeae species, 4 SHMTs of subgenomes A of wheat were used as query genes for micro-collinearity analysis; the results showed that homologs of TaSHMT1A-1, TaSHMT2A-1, and TaSHMT4A-1 were found in the collinearity regions of T. urartu, Ae. Tauschii, subgenomes A and B of T. dicoccoides and subgenomes B and D of wheat (Supplementary Figure 1). This indicated that the evolution process of most SHMTs and their adjacent genes in the micro-collinearity regions of common wheat is relatively conservative, and it mainly enters common wheat through polyploidy. However, when TaSHMT3A-1 was used as a query gene, the micro-collinearity relationship showed that no homolog of TaSHMT3A-1 and its two neighboring genes was found in the collinearity region of T. urartu, and there were 27 genes inserted into the collinearity region of T. urartu (Figure 3A). The micro-collinearity relationship was further analyzed by removing the genome of T. urartu (Figure 3B). The results showed that TaSHMT3A-1 and its homologs showed high similarity in the micro-collinearity regions of different genomes, suggesting that the corresponding SHMT and its two neighboring genes in T. urartu may have been lost. Some duplication events in the local region also occurred during the evolution of T. urartu, or a poor reference genome sequence of T. urartu caused this phenomenon.

Based on the fact that there was an SHMT gene, TdSHMT7B-1, which was on chromosome 7B of T. dicoccoides, while other analyzed Triticeae species did not contain the SHMT gene in the homologous group 7, and the number of exons and conserved motifs of this gene was less than that of other genes in the same subclass, TdSHMT7B-1 was used as query gene for micro-collinearity analysis. However, there was no micro-collinearity of TdSHMT7B-1 and its neighboring genes with other analyzed Triticeae species; furthermore, TdSHMT7B-1 was a gene pair with TRITD7Bv1G144890 of Triticum turgidum. The above results suggest that SHMT in tetraploid wheat may have experienced different evolution processes from other SHMT genes.

DIVERGE v3.0 software was used to analyze the functional divergence of SHMT proteins in different subclasses (Table 3). The results showed that the θ value of Type -I ranged from 0.3337 to 0.6333, the degree of functional differentiation among subclasses fluctuated wildly, and the standard error was 0.06. The P-values of different subclasses were lower than 0.01, indicating that functional differentiation was caused by the change of evolution rate among subclasses (Gu et al., 2013; Xu et al., 2021). The Type -II analysis showed that θ ranged within 0.1584–0.3326 (Table 3), except class I a/I b. The differences among other subgroups reached a significant level (P < 0.01), indicating that there was functional differentiation caused by the constant evolution rate, but the change of corresponding amino acid characteristics among these subgroups meant that there were some changes of critical amino acid sites which caused the functional divergence among the above subclasses. The above results indicate that gene functional divergence of different classes (except class I a/I b) came from both changing of some critical amino acid sites and evolution rate.

The expression patterns of gene family members are helpful in predicting their potential biological functions. To elucidate the potential role of TaSHMTs, their expression patterns were studied by qRT-PCR. Expression patterns of TaSHMTs under two biotic stresses (powdery mildew pathogen and F. graminearum), abiotic stresses (PEG, NaCl, and 4°C), in different tissues (roots, stems, and leaves of the seedling stage) and two hormone treatments (Abscisic acid: ABA and H₂O₂) were analyzed. Because the sequence similarity of the three copy genes in different subgenomes of wheat is high, qRT-PCR primers could not effectively distinguish the three copy genes; therefore, the gene of subgenomes A was used to represent the relative expression of the three copy genes. The expression of TaSHMTs on chromosomes 4A, 4B, and 4D could not be detected in different tissues, so the relative expression data of TaSHMTs gene in homologous group 4 were not provided in this study.

The expression pattern of TaSHMT genes in different tissues, under various abiotic stresses and response to ABA and H₂O₂ treatments is shown in Figure 4. The expression of analyzed TaSHMTs in stem and leaf were higher than that in root (Figure 4A). With PEG and NaCl treatments, the relative expression of TaSHMT2A-1 and TaSHMT3A-1 was down-regulated, and the difference was significant, but the expression of TaSHMT1A-1 did not change significantly (Figures 4B,C). Under the cold treatment, the relative expression of all TaSHMT genes was up-regulated (Figure 4D). With ABA treatment, the expression patterns of the three TaSHMT genes were different, the relative expression of TaSHMT1A-1 did not change significantly at different time points (Figure 4E); TaSHMT2A-1 was down-regulated, TaSHMT3A-1 showed a rapid increase and reached peak levels at 12 h, then returned to the original level at 24 h (Figure 4E). The transcriptional responses of TaSHMT genes to H₂O₂ showed that TaSHMT1A-1 (39-fold), TaSHMT2A-1 (10-fold), and TaSHMT3A-1 (8-fold) were highly upregulated at 12 h and then returned to the original level at 24 h (Figure 4F).

After Bgt inoculation, the relative expression levels of TaSHMT1A-1 and TaSHMT3A-1 were down-regulated at 6 h, then reached the expression peak at 24 h; the relative expression levels of TaSHMT2A-1 were up-regulated at 6 h and then reached the expression peak at 24 h (Figure 5A); however, the absolute times of relative expression change were small. The expression patterns of TaSHMTs in the FHB resistant cultivar Sumai 3 and the susceptible cultivar Jimai 22 at different times after infection with F. graminearum were further analyzed by qRT-PCR. For Sumai 3 and Jimai 22, all the relative expressions of the three TaSHMT1A-1, TaSHMT2A-1, and TaSHMT3A-1 were up-regulated in both materials (Figures 5B,C). Remarkably, the transcript levels of TaSHMT3A-1 rapidly reached levels at 48 and 72 h were 46- and 101-fold higher than at 0 h in Sumai 3, respectively; furthermore, although the expression of TaSHMT3A-1 in the susceptible material Jimai22 was also up-regulated, the levels that at 48 and 72 h were threefold higher than at 0 h. Therefore, we deduced that TaSHMTs, especially the TaSHMT3A-1 might play an important role in the resistance response to FHB.

Previous studies have shown that the Jasmonic Acid signaling pathway is related to wheat FHB resistance (Xiang et al., 2011), and exogenous ABA treatment can increase wheat sensitivity to FHB (Qi et al., 2016). In this study, TaSHMT3A-1 gene quickly responded to the induction of ABA and MeJA (Figure 5A and Supplementary Figure 2), and the relative expression of TaSHMT3A-1 in the FHB resistant cultivar Sumai 3 was significantly increased by 101 times at 72 h. Therefore, TaSHMT3A-1 was selected to further analyze its potential role in wheat resistance to FHB. The BSMV-VIGS system was used to further characterize the function of the TaSHMT3A-1 gene in common wheat Bainong207. BSMV: SHMT3A-1, which carries a 249-bp TaSHMT3A-1 fragment was used to induce target silencing. The fourth fully expanded leaves of BSMV: TaSHMT3A-1-infected plants were detached, followed by inoculation with a fresh F. graminearum spore suspension and RNA extraction. Leaves of the same age from BSMV:γ-infected plants were collected and inoculated with fresh F. graminearum as controls. Three and five days after infection, BSMV: TaSHMT3A-1-infected leaves were more susceptible to FHB than those BSMV:γ-infected individuals (Figure 6A). The expression levels of TaSHMT3A-1 were significantly decreased, by 3–6-fold, as assessed by qRT-PCR (Figure 6B). Therefore, silencing the TaSHMT3A-1 gene could increase the susceptibility of Bainong207 to FHB.

Serine hydroxymethyltransferase is involved in the reversible interconversion of Ser and Gly, participates in amination and decarboxylation reactions, and plays an important role in the cell-carbon metabolic pathway (Bauwe and Kolukisaoglu, 2003; Voll et al., 2006). Previous studies mainly focused on the function of SHMT genes under abiotic stresses such as light intensity, salt, and drought (Voll et al., 2006; Cui et al., 2019; Mishra et al., 2019); however, studies on the role of the SHMT gene family in biotic stresses and its evolutionary relationship in important crops at the genomic level have been limited. Wheat is a heterologous hexaploid composed of three subgenomes A, B, and D. Studies reveal that hybridization of T. urartu (AA) and an Ae. speltoides-related grass (BB) produced a tetraploid ancestor T. turgidum (AABB). After the second hybridization of Ae. tauschii (DD) and T. turgidum, the hexaploid ancestor (AABBDD) of wheat was formed (Shewry, 2009; Allaby et al., 2017). With the related species of wheat being sequenced, the evolutionary history of wheat has been relatively clear (Avni et al., 2017; Luo et al., 2017; IWGSC, 2018; Ling et al., 2018), and it is an important goal of future research to infer when some wheat gene subfamilies expanded and occurred during the evolutionary process (Schilling et al., 2020).

In this study, bioinformatics methods were used to comprehensively analyze the SHMT gene family in wheat, and to investigate its potential role in development and biotic and abiotic stresses. A total of 64 genes containing the complete conserved domain of SHMT (PF00464) were screened from the genomic data of six monocotyledons and four dicotyledons and divided into four subclasses. The number of SHMT genes identified from the whole genome of H. vulgare (HH), T. urartu (AA), Ae. tauschii (DD), T. dicoccoides (BBAA) and wheat (BBAADD) were 4, 3, 4, 9, and 12, respectively. The proportion of SHMT number among species was consistent with the proportion of the corresponding genomic multiples; the number of genes without copy genes was almost the same (Table 1). Gene replication is the most common mechanism for gene family extension (Cannon et al., 2004; Yu et al., 2020) and results in functional differentiation, which is critical for environmental adaptation and speciation (Conant and Wolfe, 2008). This study showed no tandem duplication or segmental replication in wheat SHMT gene. Meanwhile, macro-collinear and micro-collinear analysis showed that the SHMT gene of wheat was derived from a polyploidization process rather than self-replication (Figure 3 and Supplementary Figures 1, 3). Interestingly, there was a TdSHMT7B-1 on chromosome 7B of T. dicoccoides and there was no gene pair in wheat and other related species analyzed in this study; however, there was a homologous gene TRITD7Bv1G144890 on chromosome 7B of heterotetraploid species T. turgidum. We extended our focus to the incomplete domain SHMT of Triticeae species. The chromosomes 7B and 4B of wheat; 7A of T. dicoccoides; 4D of Ae. tauschii, and 4A and 7A of T. urartu all had one incomplete SHMT gene (Supplementary Table 2). Compared with the SHMT gene in the same subgroup, TdSHMT7B-1 lacked motifs 7, 9, and 17 (Figure 2B). Structural and conserved motif analysis revealed that motif, intron/exon loss or gain may occur during the evolution of Triticeae SHMT and similar events also appeared in SHMT of soybean (Lakhssassi et al., 2019). Some studies suggest that, in monocots, a special gene loss event occurred in class II a, and a gene duplication event in class I a compensated for this loss (Lakhssassi et al., 2019). However, this study does not support the above viewpoint. Instead, we speculate that the loss of SHMT gene did occur in class II a of monocotyledons, but there was no particular replication in class I a to make up for it. On the contrary, the loss or regression of SHMT gene occurred in class I a of monocotyledons during the process of evolution. This difference in genetic evolution may be because the effect of the circadian clock on daily growth rhythms is different between monocotyledons and dicotyledons (Campoli et al., 2012). However, this speculation needs further research.

In the process of evolution, duplicated genes often undergo functional differentiation, leading to neofunctionalization, subfunctionalization or non-functionalization (Prince and Pickett, 2002). Compared with the traditional “paralog” and “ortholog,” “homolog” can more accurately and reasonably understand the relationship of genes in polyploid species (Chen Y. et al., 2020). To better understand the evolution of SHMT, Diverge3.0 was used to analyze the functional differentiation of the SHMT gene in each subclass. Type I differentiation represents the difference in gene evolution rate, and Type II differentiation represents the change in physical and chemical properties of amino acids (Yang et al., 2020). The results showed that the gene functional differences among the subgroups (except class I a/I b) came not only from changes of some key amino acid loci, but also from changes in evolution rate. The results of class I a/I b showed that the functional differentiation between the two subgroups was mainly Type I, and it was mainly determined by the sites conserved in one subclass and not conserved in the other subclass. Gene expansion provides conditions for the formation of new functions. The sequence of one gene remains relatively stable and the sequence of the other expanded gene changes, and functional differentiation occurs.

Since the sequence similarity of the three copy genes in different subgenomes of wheat is high (Pfeifer et al., 2014), we used the universal qRT-PCR primers to analyze the expression levels of corresponding three copy genes. The SHMT gene was shown to be closely related to photosynthesis, and with reduced activities of SHMT will usually show severe growth retardation (Heineke et al., 2001; Somerville, 2001). In the present study, the results showed that the expression level of TaSHMTs in stems and leaves were higher than in roots, which may be because SHMT plays an important role in plant one-carbon metabolism and photorespiration, and green tissues are important parts of photorespiration. The expressions of TaSHMT2A-1 and TaSHMT3A-1 were down-regulated under the stress of PEG and NaCl and TaSHMT1A-1 did not respond to these two treatments. These results are consistent with previous research results in diploid wheat Triticum monococcum (Bhuiyan et al., 2007), speculating that SHMTs also have similar functions as in T. monococcum. With ABA treatment, the expression levels of TaSHMT2A-1 were down-regulated, and TaSHMT3A-1 showed a rapid increase (Figure 4E). Furthermore, TaSHMT1A-1, TaSHMT2A-1, and TaSHMT2A-1 significantly responded to H₂O₂ (Figure 4F). Previous studies showed that AtSHMT1 could prevent cell death and reduce reactive oxygen species accumulation during salt stress (Zhou et al., 2012). Furthermore, AtSHMT1 play a crucial role in plant abiotic stress tolerance and ABA-induced stomatal movements (Liu et al., 2019), suggesting TaSHMT2A-1 and TaSHMT3A-1 may be involved in H₂O₂ - and ABA-induced tolerance to abiotic stress.

Wheat FHB is a serious disease of wheat, which seriously affects the human food and animal feed security (Li et al., 2019; Su et al., 2019; Wang et al., 2020). In this study, we found that the expression of TaSHMT3A-1 in the FHB resistant cultivar Sumai 3 was significantly increased by 101 times after infection by F. graminearum at 72 h (Figure 5B). When the TaSHMT3A-1 was silenced in Bainong207, the results showed that silencing the TaSHMT3A-1 gene could increase the susceptibility of Bainong207 to FHB. Given that the link of the transcriptional response of TaSHMT3A-1 to F. graminearum, ABA, MeJA and the result of BSMV-VIGS (Figures 4E,F, 6 and Supplementary Figure 2), it is reasonable for us to conclude that TaSHMT3A-1 can enhance wheat resistance to FHB. The previous study indicated that Arabidopsis shmt1-1 mutants were more susceptible than control plants infection with biotrophic and necrotrophic pathogens (Moreno et al., 2005). The soybean cyst nematode (SCN) resistance major genetic locus Rhg4 encodes an SHMT protein (Kandoth et al., 2017; Lakhssassi et al., 2019). Based on the above research, we speculate that the SHMT gene may play a vital role in plant disease and insect resistance. In this study, the function of TaSHMT3A-1 on FHB was verified through a leaf BSMV-VIGS assay. The effects and contribution of TaSHMT3A-1 gene on wheat FHB need to be further verified through stable genetic transformation plants. Although the function of the wheat TaSHMTs needs to be elucidated, this study revealed that TaSHMTs might have diverged new functions during the course of evolution.