Thaumatin-like Protein (TLP) Genes in Garlic (Allium sativum L.): Genome-Wide Identification, Characterization, and Expression in Response to Fusarium proliferatum Infection

Plant antifungal proteins include the pathogenesis-related (PR)-5 family of fungi- and other stress-responsive thaumatin-like proteins (TLPs). However, the information on the TLPs of garlic (Allium sativum L.), which is often infected with soil Fusarium fungi, is very limited. In the present study, we identified 32 TLP homologs in the A. sativum cv. Ershuizao genome, which may function in the defense against Fusarium attack. The promoters of A. sativum TLP (AsTLP) genes contained cis-acting elements associated with hormone signaling and response to various types of stress, including those caused by fungal pathogens and their elicitors. The expression of AsTLP genes in Fusarium-resistant and -susceptible garlic cultivars was differently regulated by F. proliferatum infection. Thus, in the roots the mRNA levels of AsTLP7–9 and 21 genes were increased in resistant and decreased in susceptible A. sativum cultivars, suggesting the involvement of these genes in the garlic response to F. proliferatum attack. Our results provide insights into the role of TLPs in garlic and may be useful for breeding programs to increase the resistance of Allium crops to Fusarium infections.


In Silico Genome-Wide Identification of TLP Genes in A. sativum cv. Ershuizao
A total of 32 complete TLP gene sequences were detected in the genome and transcriptome of A. sativum cv. Ershuizao (PRJNA606385, PRJNA607255) and denoted AsTLP1-32 (Table 1, Supplementary Table S1). Most of the genes (AsTLP3-20) are located on chromosome 2 and seven (AsTLP1, 2, and 21-25)-on chromosomes 1, 4, 6, 7, and 8 ( Figure 1), whereas the other seven (AsTLP26-32) have been found in scaffolds and do not match any chromosomes (assembly Garlic.V2.fa; Table 1). The sizes of the identified genes range from 615 to 2572 bp and those of the coding sequences (CDSs)-from 615 to 972 bp (Table 1).
In total, 15 conserved motifs were identified in putative AsTLPs, including motif corresponding to signal peptide. MEME analysis indicated that phylogenetically related

In Silico Analysis of AsTLP Transcription
The expression of AsTLP genes in various tissues (roots, bulbs, stems, leaves, buds, flowers, and sprouts) of A. sativum cv. Ershuizao was evaluated based on transcriptomics data (PRJNA607255). As the bulb is the part of garlic most susceptible to Fusarium infection and shows visible symptoms, gene expression analysis in this organ was more detailed and included eight developmental stages (from day 192 to day 227 after planting). Figure 4 shows mRNA expression levels of 16 AsTLP genes, for which Fragments Per Kilobase of transcript per Million mapped reads (FPKM) was ≥ 10 in at least one of the organs. AsTLP1-3, 5-11, and 16 had overall similar expression patterns (with the exception of AsTLP1 and 9 in the roots and AsTLP5-7, 9, 11, and 16 in the leaves), whereas those of AsTLP20, 21, and 23-25 were different.

AsTLP Promoter Analysis
Considering the stress-and hormone-responding role of TLPs [49], we searched for respective cis-acting elements in the 5'-UTR and promoter regions (1 kb upstream of the initiation codon) of AsTLP genes. The results indicated that the AsTLP regulatory regions contained 11 hormone-responsive and 10 stress-responsive elements ( Figure 5). Among the former, the most common were the ABA-responsive element (ABRE) [54], ethylene-responsive element (ERE) critical for ethylene-regulated transcription in plants [55], a highly conserved prolamin-box (P-box) recognized by the trans-activator P-box-binding factor (PBF) during endosperm development [56], and the CGTCA motif involved in MeJA and osmotic stress responsiveness [57]. Among the stress-responsive elements, the most common were anaerobic-responsive element (ARE) involved in the activation of anaerobic gene expression [58] and TC-rich repeats found in promoters of many plant disease-resistance genes [59]. The promoters most enriched in the regulatory elements were those of AsTLP11 (2 ABRE, 4 CGTCA motifs, 2 ERE, and 3 ARE), AsTLP17 (4 ERE and 3 ARE), AsTLP1, 24, and 25 (7, 4, and 3 ABRE, respectively), and AsTLP6 and 10 (4 and 5 ARE, respectively). The promoters of five genes (AsTLP8, 9, 23, 24, and 32) contained low-temperature-responsive (LTR) elements; among them, AsTLP23 had the highest number (6). The W-box, which responds to fungal elicitors and wounding, was About a third of the AsTLP genes (AsTLP1-3, 5-11, and 16) were strongly expressed in the stems and leaves. In the roots, the highest mRNA levels were observed for AsTLP1, 9, and 23 and the lowest-for AsTLP21, 24, and 25; the latter had the highest expression in the floral buds and stage 1 bulbs, which gradually decreased during bulb development. The expression of AsTLP20 and 23 in the bulbs showed bell-shape changes, with the maximum at stages 3 and 6, respectively. According to the expression profiles, the AsTLP genes could be divided into two groups: AsTLP1-3, 5-11, and 16 with the highest expression in the stems, leaves, and roots, and AsTLP20, 21, 24, and 25 with that in the developing bulbs, floral buds, and flowers. The first group includes all genes of phylogenetic clades V and VI, as well as one gene each from clades I and III; members of the second group are phylogenetically heterogeneous (Figure 2a). The expression pattern of the AsTLP23 gene had the features of both groups (Figure 4).

AsTLP Promoter Analysis
Considering the stress-and hormone-responding role of TLPs [49], we searched for respective cis-acting elements in the 5 -UTR and promoter regions (1 kb upstream of the initiation codon) of AsTLP genes. The results indicated that the AsTLP regulatory regions contained 11 hormone-responsive and 10 stress-responsive elements ( Figure 5). Among the former, the most common were the ABA-responsive element (ABRE) [54], ethyleneresponsive element (ERE) critical for ethylene-regulated transcription in plants [55], a highly conserved prolamin-box (P-box) recognized by the trans-activator P-box-binding factor (PBF) during endosperm development [56], and the CGTCA motif involved in MeJA and osmotic stress responsiveness [57]. Among the stress-responsive elements, the most common were anaerobic-responsive element (ARE) involved in the activation of anaerobic gene expression [58] and TC-rich repeats found in promoters of many plant disease-resistance genes [59]. The promoters most enriched in the regulatory elements were those of AsTLP11 (2 ABRE, 4 CGTCA motifs, 2 ERE, and 3 ARE), AsTLP17 (4 ERE and 3 ARE), AsTLP1, 24, and 25 (7, 4, and 3 ABRE, respectively), and AsTLP6 and 10 (4 and 5 ARE, respectively). The promoters of five genes (AsTLP8, 9, 23, 24, and 32) contained low-temperature-responsive (LTR) elements; among them, AsTLP23 had the highest number (6). The W-box, which responds to fungal elicitors and wounding, was found in only three genes (AsTLP1, 24, and 30). The rarest elements were ABA-responsive CARE (in AsTLP1 and 21) and auxin-responsive AuxRE (in AsTLP8 and 9) ( Figure 5). No pronounced similarity between the set of elements and the structure of genes was found, with the exception of some groups of genes, for example, AsTLP5 and 6 (clade V), 8 and 9 (clade VI), or 22, 28, and 29 (clade I).

AsTLP1-32 Expression in cv. Sarmat and Strelets Infected with F. proliferatum
To elucidate the role of the identified AsTLP1-32 genes in the garlic response to fungal infection, we performed AsTLP expression profiling in garlic cultivars resistant (cv. Sarmat) and susceptible (cv. Strelets) to FBR [60]. The AsTLP7 corresponds to clade V group of genes (AsTLP5-7 and 11), AsTLP8-to set of AsTLP3, 8, and 10 genes, AsTLP12-to set of AsTLP12, 15, 17, 22, and 26-29 genes, AsTLP19-to AsTLP13, 18, 19, and 30 genes, and AsTLP31-to AsTLP31 and 32 (clade II) genes (Supplementary Table S2). Gene expression was analyzed in the roots, stems (basal plates), and cloves of cv. Sarmat and cv. Strelets at 24 and 96 h post infection (hpi) with F. proliferatum and compared with uninfected control (Figures 6 and 7).  (Figures 6 and 7).  The results revealed that the AsTLP genes were transcribed in all analyzed organs, except for AsTLP24 not expressed in the cloves ( Figure 6). To determine the expression of individual genes in clades V and II, we performed sequencing of PCR-amplified products, which revealed that only AsTLP7 and AsTLP31 were transcribed in the tissues of both cultivars. For the other groups of highly homologous AsTLPs, the amplified PCR products were invariable within each group. The results revealed that the AsTLP genes were transcribed in all analyzed organs, except for AsTLP24 not expressed in the cloves ( Figure 6). To determine the expression of individual genes in clades V and II, we performed sequencing of PCR-amplified products, which revealed that only AsTLP7 and AsTLP31 were transcribed in the tissues of both cultivars. For the other groups of highly homologous AsTLPs, the amplified PCR products were invariable within each group.
In both uninfected cultivars, the highest expression was observed for AsTLP7, 8, and 31 ( Figure 7). Analysis of time-dependent changes revealed that most genes had similar expression dynamics in the two cultivars. Thus, the mRNA levels of AsTLP1, 7, 8, and 16 remained the same or increased from 24 to 96 hpi, whereas those of AsTLP12, 19-21, and 24 were either the same or decreased. However, the expression of AsTLP23 and 31 was increased from 24 to 96 hpi in the cloves of FBR-resistant cv. Sarmat but decreased in those of FBR-susceptible cv. Strelets (Figure 7).
Symptoms of F. proliferatum infection were observed only in FBR-sensitive cv. Strelets after 96 hpi: the roots were covered with white fungal mycelium, indicating the external growth of the pathogen. Infection changed the expression patterns of the AsTLP Figure 7. Time-dependent expression of AsTLP genes in the roots, stems (basal plates), and cloves of uninfected garlic cultivars resistant (cv. Sarmat) and susceptible (cv. Strelets) to FBR. The data were normalized to GAPDH and UBQ mRNA levels and presented as the mean ± SE (n = 3); * p < 0.01 indicates the difference between 24 and 96 h.
In both uninfected cultivars, the highest expression was observed for AsTLP7, 8, and 31 ( Figure 7). Analysis of time-dependent changes revealed that most genes had similar expression dynamics in the two cultivars. Thus, the mRNA levels of AsTLP1, 7, 8, and 16 remained the same or increased from 24 to 96 hpi, whereas those of AsTLP12, 19-21, and 24 were either the same or decreased. However, the expression of AsTLP23 and 31 was increased from 24 to 96 hpi in the cloves of FBR-resistant cv. Sarmat but decreased in those of FBR-susceptible cv. Strelets (Figure 7).
Symptoms of F. proliferatum infection were observed only in FBR-sensitive cv. Strelets after 96 hpi: the roots were covered with white fungal mycelium, indicating the external growth of the pathogen. Infection changed the expression patterns of the AsTLP genes compared to uninfected control. Thus, in FBR-resistant cv. Sarmat, the transcription of all genes was upregulated at 24 hpi, except for that of AsTLP1 in the roots and stems, which was upregulated at 96 hpi. These data suggest a correlation between the expression level of AsTLP1 and the appearance of disease symptoms. In the roots, the expression of most genes was increased from 24 to 96 hpi, except for that of AsTLP12 and 19 (unchanged) and AsTLP24 (decreased). Similar dynamics were observed in the cloves, except for AsTLP12 and 31, which were downregulated with time. In the stems, the expression of AsTLP7-9 and 31 decreased and that of AsTLP19 and 24 increased from 24 to 96 hpi, whereas that of AsTLP16 was unchanged ( Figure 6).
In FBR-susceptible cv. Strelets, all AsTLP genes were upregulated at 24 hpi compared to uninfected control, except for AsTLP9 in the roots and stem and AsTLP31 in the cloves, which were upregulated at 96 hpi ( Figure 6). The time-dependent AsTLP expression dynamics in cv. Strelets significantly differed from that in cv. Sarmat: the expression of AsTLP7-9, 20, 21, 23, and 24 was either unchanged or downregulated from 24 to 96 hpi, whereas tissue-dependent changes were observed for the other genes. Thus, AsTLP12 and 31 were downregulated with time in the roots and cloves and upregulated in the stems, AsTLP16 was downregulated in the roots and stems and upregulated in the cloves, and AsTLP19 was downregulated in the cloves and upregulated in the stems ( Figure 6).

Cloning and Characterization of CDSs of AsTLP Genes Differentially Expressed in FBR-Sensitive and -Resistant Cultivars
Considering the differential expression of some AsTLP genes in cv. Sarmat and cv. Strelets in response to F. proliferatum infection (Figure 6), we investigated the association of these genes with FBR resistance. AsTLP7, 16, and 20 were upregulated by the infection in the roots of cv. Sarmat, whereas AsTLP21 and 23 were upregulated in those of cv. Strelets. The CDSs of these genes were amplified, cloned, and sequenced, and the data deposited in NCBI GenBank (Table 2).  The cv. Ershuizao AsTLP sequences were used as references. The amino acid sequences of AsTLP7, 16, 20, 21, and 23 from cv. Sarmat and Strelets did not differ in size from those of cv. Ershuizao AsTLPs. Compared to the cv. Ershuizao AsTLP genes, the cv. Sarmat and Strelets genes contained 4-18 single nucleotide polymorphisms (SNPs). Among them, non-synonymous SNPs leading to amino acid substitutions were found in AsTLP7 (4), AsTLP16 (1, resulting in radical A159G), AsTLP20 (4, including radical G62R), AsTLP21 (1), and AsTLP23 (1) of both cultivars, whereas a cultivar-specific non-synonymous SNP (670T>C, leading to radical Y224H) was detected in AsTLP20 of cv. Strelets (Table 2).

Analysis of Regulatory Regions in the AsTLP Genes Differentially Expressed in cv. Sarmat and Strelets after F. proliferatum Infection
The regulatory sequences (promoter and 5 -UTR) of the AsTLP7, 16, 20, 21, and 23 genes, which were differently expressed in cv. Sarmat and Strelets in response to F. proliferatum infection, were amplified, sequenced, and analyzed for hormone-and stressresponsive elements (Figure 8). The regulatory region of only the AsTLP23 gene of cv.
Sarmat and Strelets differed from the AsTLP23 of cv. Ershuizao by the absence of the CGTCA-motif and one of the six LTR-elements (Figure 8). In general, the sequences of the regulatory regions in AsTLP7, 16,20,21, and 23 were the same in cv. Sarmat and Strelets, but differed from those in cv. Ershuizao: there were 3 indels (11, 4, and 2 bp) and 20 SNPs in AsTLP7, 8 SNPs in AsTLP20, 3 indels (2, 2, and 7 bp) and 7 SNPs in AsTLP21, and 5 SNPs in AsTLP23.

Analysis of Regulatory Regions in the AsTLP Genes Differentially Expressed in cv. Sarmat and Strelets after F. proliferatum infection
The regulatory sequences (promoter and 5′-UTR) of the AsTLP7, 16, 20, 21, and 23 genes, which were differently expressed in cv. Sarmat and Strelets in response to F. proliferatum infection, were amplified, sequenced, and analyzed for hormone-and stress-responsive elements (Figure 8). The regulatory region of only the AsTLP23 gene of cv. Sarmat and Strelets differed from the AsTLP23 of cv. Ershuizao by the absence of the CGTCA-motif and one of the six LTR-elements (Figure 8). In general, the sequences of the regulatory regions in AsTLP7, 16,20,21, and 23 were the same in cv. Sarmat and Strelets, but differed from those in cv. Ershuizao: there were 3 indels (11, 4, and 2 bp) and 20 SNPs in AsTLP7, 8 SNPs in AsTLP20, 3 indels (2, 2, and 7 bp) and 7 SNPs in AsTLP21, and 5 SNPs in AsTLP23.  The mutations found in the AsTLP7, 20, and 21 regulatory sequences did not affect cis-elements, whereas those in the AsTLP23 sequences did. Thus, two SNPs (at positions -674G>A and -676A>G) resulted in the loss of the CGTCA motif associated with MeJA response, another SNP (-521C>T) created the GARE motif (gibberellin-responsive element), and two SNPs (-640G>A and -641C>A) eliminated an LTR element (Figure 8).

Discussion
Garlic (A. sativum L.) belongs to the Allium genus, which comprises 971 species, thus being one of the largest genera of the Amaryllidaceae family widely distributed in the Northern Hemisphere from dry subtropics to boreal zones [61,62]. Significant garlic crop losses regularly occur worldwide, mostly because of diseases caused by soil fungi of Fusarium spp. [63,64], as well as by drought, osmotic, and cold stresses [35,65]. In many plants, TLPs play a critical role in the defense against both abiotic and biotic stresses [66], suggesting that they also perform a similar role in garlic and may potentially be used in Allium crop breeding programs.
In the present study, we identified and characterized 32 TLP genes in the A. sativum cv. Ershuizao genome [67] and amplified five of them (CDSs and 5 -regulatory regions) from garlic cv. Sarmat and Strelets differing in the susceptibility to FBR. All the genes encoded proteins homologous to the A. thaliana pathogenesis-related TLPs, including osmotin 34, implying the same functional activity. Among the identified AsTLP genes, eight were located in six chromosomes, where they probably have emerged because of local segmental duplications, and 17 were tandemly arranged in a single cluster in chromosome 2 (Figure 1), which may be due to tandem duplication through multiple episodes of unequal crossovers. This hypothesis is consistent with the clade distribution of AsTLP genes, when several tandemly duplicated genes were clustered together (Figure 2a), suggesting their origin from recent duplication events. Structure analysis revealed that most AsTLPs (26) had no introns or contained only a few (1-2) of them (Figure 2a). The intron-poor and intronless genes transcribed with less or without splicing are considered to have originated from intron-rich genes as a part of the adaptation strategy aimed to provide fast responses to different biotic and abiotic stresses, which is consistent with strong induction of intronless genes after stress [68].
All putative AsTLP1-32 proteins contained the conserved thaumatin domain and TLP family motif signatures (Figures 2b and 3), suggesting conservation of garlic TLPs in terms of regulation through stress/defense-related signaling and functional activity to hydrolyze fungal β-1,3-glucan. The presence of N-terminal signal peptides indicates targeting of mature AsTLPs to the apoplastic pathway, which is supported by GO analysis predicting AsTLP extracellular secretion and is consistent with a previous report on TLP localization to the intercellular space [69]. These data agree well with the function of TLPs in hydrolyzing the fungal cell wall and their GO-predicted role in defense response to fungal infection.
All AsTLP1-32 proteins contained fragments with significant homology to previously described thaumatin-like allergens from various plant species [31,[51][52][53]. Given the reports of IgE-mediated human hypersensitivity to raw garlic [70], it may be attributed, at least in part, to the AsTLPs expressed in the cloves ( Figure 6).
Seven AsTLPs (AsTLP5-11) were found to be homologous to A. thaliana osmotin 34 (Table 1), which functions in the initiation of ABA responses, including the regulation of ABA-induced proline synthesis [10], suggesting the involvement of AsTLP5-11 in the garlic response to osmotic stress.
It is known that fungi-mediated biotic stress activates plant immune system through perception of pathogen-associated molecular patterns [71] with following formation of reactive oxygen species, induction of salicylic acid and jasmonic acid signaling, and upregulation of PR genes [72]. The promoters of all PR5 family AsTLP genes contained cis-regulatory elements associated with the activation of immune mechanisms through responses to stresses (anaerobic conditions, dehydration, low temperature, salinization, heavy metals, and wounding), elicitors, and hormones, (Figure 6). A similar set of elements has also been found in the promoters of Gossypium barbadense and Rosa chinensis TLP genes [73,74] and of the other garlic PR1-5 genes [50,60]. Accordingly, the expression of AsTLP genes may be triggered by phytohormones (ethylene, salicylic acid, ABA, and MeJA), and stresses.
Ethylene and ABA homeostasis is known to be modulated in host plant in the response to fungal infection [60,75]. In accordance with this, the ERE and ABRE elements were found in the promoters of more than half of AsTLPs ( Figure 5), which suggests these genes activation by ethylene and ABA after Fusarium attack.
Ethylene is known to crosstalk with ABA in response to abiotic stresses [75], such as low temperature (freezing), drought, and salinity, which cause osmotic stress and, consequently, activation of SnRK2 kinases as key signal transducers in the ABA pathway [76]. Considering the set of discovered stress responsive elements, it can be assumed that almost all AsTLP genes are involved in the response to osmotic stress, except for AsTLP7, 9,16, and 20, promoters of which do not contain neither the CGTCA motif, nor the DRE1, LTR, and MBS elements. Only five genes (AsTLP8, 9, 23, 24, and 32) may be involved in the response to cold stress (especially AsTLP23, which contains 5 LTR elements) ( Figure 5).
Some AsTLP genes within the same clade (AsTLP5-7, 11 or AsTLP2, 3, 8-10) (Figure 2a) tended to be dissimilar in the promoter cis-regulatory motif patterns ( Figure 5), suggesting that their transcription in response to adverse conditions might be differentially regulated.
Analysis of TLP genes with the established function in plant defense has predicted two unique combinations of regulatory elements to distinguish between abiotic and biotic stress responses: genes regulated by abiotic stresses contain the ABRE element, whereas those responding to fungal attack do not [77]. These data suggest that AsTLP3, 4, 7, 14-17, 20-23, 26-29, and 31, whose promoters do not have ABRE ( Figures 5 and 8), may be involved in the antifungal defense of garlic.
Among the five AsTLP genes amplified from FBR-resistant cv. Sarmat and FBRsusceptible cv. Strelets, only one contained a cultivar-specific SNP leading to an amino acid substitution ( Table 2), suggesting that the antifungal resistance in garlic may not be attributed to AsTLP mutations, but rather to transcriptional and translational regulation of gene expression. This notion is supported by the differential expression of AsTLP genes in FBR-resistant and -susceptible cultivars infected with F. proliferatum. The most pronounced dissimilarity was the time-dependent upregulation of AsTLP7-9 (homologs of OSMOTIN 34) and AsTLP21 in the roots of cv. Sarmat and their downregulation in those of cv. Strelets (Figure 6), suggesting that the expression dynamics of these AsTLP genes may underlie the difference in antifungal resistance between the cultivars. Considering that the promoters of AsTLP8 and 9 contain the ABRE element, whereas those of AsTLP7 and 21 do not, the latter are more likely to have a role in the garlic defense responses against F. proliferatum.
It has been suggested that TLP genes can serve as molecular markers associated with resistance to fungal diseases [69]. Thus, our results could be useful for breeding programs aimed on increasing the resistance of garlic crops to Fusarium spp. by, for example, generating cultivars that overexpress the respective AsTLP genes, which may protect against fungal infections. In addition, it should be noted that, besides TLPs, the fungicidal effects of Allium plants are associated with chitinases and endo-1,3-β-glucanases, as well as miRNAs involved in positive (miR164a, miR168a, and miR393) and negative (miR394) regulation of resistance to Fusarium [48,50,60,78]. Moreover, Allium roots produce antifungal volatiles such as 2-methyl-2-pentenal and organosulfur compounds, as well as non-volatiles such as spirostanol, furostanol, and steroidal saponins, which inhibit Fusarium growth [60,79,80].

In Silico Identification and Structural Characterization of AsTLP Genes
The search for TLP genes was performed in the A. sativum cv. Ershuizao wholegenome (NCBI accession number: PRJNA606385, assembly Garlic.V2.fa) and transcriptome (PRJNA607255) sequences [68]. The thaumatin domain (http://pfam.xfam.org/family/ PF00314; accessed on 25 October 2021) was used as reference. All the selected sequences contained start and stop codons and full-length catalytic domains.
Multiple sequence alignment and structural analyses of the AsTLP genes and encoded proteins were conducted with MEGA 7.0.26 [81]. The phylogenetic dendrogram was constructed based on protein sequences using the MEGA 7.0.26 (ML method); confidence for tree topologies was estimated by bootstrap values of 1000 replicates.
To predict exon-intron composition, AsTLP genes and CDSs were analyzed with GSDS v2.0 [82]. Putative proteins were characterized by molecular weight, pI, and grand average

In Silico mRNA Expression Analysis
The expression of AsTLP genes in the roots, bulbs, stems (basal plates), leaves, buds, flowers, and sprouts was analyzed based on A. sativum cv. Ershuizao RNA-Seq data (FPKM; ID: PRJNA607255), normalized as FPKM [66], and visualized using Heatmapper [83]. Only transcripts with an average FPKM value ≥ 10 in at least one of the organs were used for heatmap construction.

Gene Identification
To amplify the AsTLP CDSs from garlic cultivars, gene-specific primers were designed based on A. sativum cv. Ershuizao transcriptomic data (NCBI project accession number: PRJNA607255) (Supplementary Table S2). CDNA (30 ng) from the roots of each cultivar accession was used as a template for PCR amplification at the following conditions: initial denaturation at 95 • C for 5 min, 35 cycles [denaturation at 95 • C for 30 s, primer annealing at 55 • C for 30 s, and extension at 72 • C for 2 min], and final extension at 72 • C for 5 min. PCR products of the expected size were purified by using the QIAEX ® II Gel Extraction kit (QIAGEN, Hilden, Germany), cloned in the pGEM ® -T Easy vector (Promega, Madison, WI, USA), and sequenced (3-5 clones for each accession) on ABI Prism 3730 DNA Sequencer (Applied Biosystems, Waltham, MA, USA) using the designed primers.

Plants, Fungi, and F. proliferatum Infection Assay
F. proliferatum was kindly provided by the Group of Experimental Mycology, Winogradsky Institute of Microbiology (Research Center of Biotechnology of the RAS, Moscow, Russia). The strain was originally isolated from the bulbs of garlic cv. Strelets; according to the pathogenicity test, the first signs of the disease appeared on the clove surface 5 days after infection [50].
Accessions of A. sativum cv. Sarmat and cv. Strelets (winter garlic of Russian breeding) resistant and susceptible to FBR, respectively, were kindly provided by the Federal Scientific Vegetable Center (Moscow region, Russia). The number of cloves used per biological replicate in the Fusarium infection assay was based on that of cloves in the bulb (5-7 for cv. Strelets and 7-11 for cv. Sarmat). In total, 12 cloves of each cultivar (6 infected and 6 uninfected) were processed (three biological replicates were used). Cloves were surfacesterilized in 70% ethanol for 3 min, rinsed with sterile water, placed in Petri dishes with wet filter paper, and incubated at +25 • C in the dark. After 72 h, active root growth was observed, and half of the cloves were infected by soaking in F. proliferatum conidial suspension (~10 6 conidia ml −1 ) for 5 min as previously described [60]), transferred to fresh Petri dishes, and incubated at +25 • C in the dark for 24 and 96 h (n = 3 cloves per each time point).The roots, stems, and cloves of the infected and uninfected samples were collected at each time point, frozen in liquid nitrogen, and stored at −80 • C. The time points were chosen according to the expression peaks of some PR genes, which were observed 1-3 days after inoculation with hemibiotrophic pathogens [84].

RNA Extraction and Quantitative Real-Time Reverse Transcription PCR (qRT-PCR)
Total RNA was extracted from individual roots, stems, and cloves (0.5 g of each tissue) using the RNeasy Plant Mini Kit (QIAGEN, Hilden, Germany), purified from genomic DNA (RNase free DNase set; QIAGEN), qualified by gel electrophoresis, and used for first-strand cDNA synthesis (GoScript Reverse Transcription System; Promega, Madison, USA) with an oligo-dT primer. RNA and cDNA concentrations were quantified by fluorimetry (Qubit ® Fluorometer, Thermo Fisher Scientific, Waltham, MA, USA) and qRT-PCR was performed in a CFX96 Real-Time PCR Detection System (Bio-Rad Laboratories, Hercules, CA, USA) with 3.0 ng cDNA, SYBR Green RT-PCR mixture (Syntol, Moscow, Russia), and specific primers (Supplementary Table S2). Because CDSs of some AsTLPs had a high degree of homology, universal primers were designed for such genes, which were grouped according to their homology (Supplementary Table S2). The following cycling conditions were used: initial denaturation at 95 • C for 5 min, 40 cycles of denaturation at 95 • C for 15 s, and annealing/extension at 60 • C for 40 s.
AsTLP gene expression was normalized using two reference garlic genes, GAPDH and UBQ [60], and the qRT-PCR results were statistically analyzed with Graph Pad Prism version 8 (GraphPad Software Inc., San Diego, CA, USA; https://www.graphpad.com/ scientific-software/prism/ (accessed on 30 October 2021)). The data were expressed as the mean ± standard deviation (SE) based on three technical replicates of three biological replicates for each combination of cDNA and primer pairs. The unequal variance (Welch's) t-test was applied to assess differences in gene expression; p < 0.01 was considered to indicate statistical significance.

Promoter and 5 -UTR Analysis
The search for specific cis-elements in the promoters and 5 -UTRs (1.0 kb regions upstream of the initiation codon) was performed using the PlantCARE database, which provides evaluation of cis-regulatory elements, enhancers, and repressors; (http://bioinformatics.psb. ugent.be/webtools/plantcare/html/; accessed on 25 November 2021).

Conclusions
We identified and characterized 32 genes encoding thaumatin-like proteins in A. sativum cv. Ershuizao genome. The AsTLP genes were distributed among six chromosomes and four scaffolds and might have been evolutionary originated from segmental or tandem duplications. AsTLP7, 16, 20, 21, and 23 homologs were amplified from garlic cultivars resistant and susceptible to Fusarium infection. The promoters of AsTLP genes contained distinct sets of cis-acting elements associated with hormone and stress reactivity, suggesting differential transcriptional regulation of garlic TLPs in response to pathogens and abiotic stresses, which was consistent with specific expression patterns of AsTLP genes in garlic cultivars infected with F. proliferatum. The transcription of AsTLP7-9, and 21 genes in the roots was downregulated in FBR-susceptible and upregulated in FBR-resistant cultivars, suggesting their particular involvement in the sensitivity of garlic to fungal infection. Our results provide the foundation for further functional characterization of the AsTLP genes using a reverse genetics strategy, and may contribute in the breeding of A. sativum cultivars with increased resistance to Fusarium infections as well as various abiotic stresses.

Funding:
The article was made with support of the Ministry of Science and Higher Education of the Russian Federation in accordance with agreement № 075-15-2020-907 date 16 November 2020 on providing a grant in the form of subsidies from the Federal budget of Russian Federation. The grant was provided for state support for the creation and development of a world-class scientific center "Agrotechnologies for the Future".
Institutional Review Board Statement: Not applicable.

Informed Consent Statement: Not applicable.
Data Availability Statement: AsTLP CDSs of A. sativum cv. Sarmat/cv. Strelets were deposited in NCBI (see Table 2).