Genome-Wide Analysis of Q-Type C2H2 ZFP Genes in Response to Biotic and Abiotic Stresses in Sugar Beet

Simple Summary A plant’s C2H2-type zinc finger proteins (C2H2-ZFPs) play crucial roles in the process of plant growth and development, as well as various stress responses. The Q-type ZFP family, which contains a conserved “QALGGH”, has been reported in many plants. Sugar beet is an important crop for sugar production. Salt stress and viral infection significantly reduce both sugar yield and processing quality of sugar beet. So far, the genome-wide analysis of Q-type C2H2 ZPFs and their expression pattern in sugar beet have not been analyzed yet. This study analyzed 35 Q-type ZFPs in sugar beet and their expression patterns under salt stress and virus. These results will provide theoretical evidence for understanding the functions of Q-type ZFPs. Abstract A plant’s Q-type C2H2-type ZFP plays key roles in plant growth and development and responses to biotic and abiotic stresses. Sugar beet (Beta vulgaris L.) is an important crop for sugar production. Salt stress and viral infection significantly reduce the root yield and sugar content of sugar beet. However, there is a lack of comprehensive genome-wide analyses of Q-type C2H2 ZFPs and their expression patterns in sugar beet under stress. In this study, 35 sugar beet Q-type C2H2 ZFPs (BvZFPs) containing at least one conserved “QALGGH” motif were identified via bioinformatics techniques using TBtools software. According to their evolutionary relationship, the BvZFPs were classified into five subclasses. Within each subclass, the physicochemical properties and motif compositions showed strong similarities. A Ka/Ks analysis indicated that the BvZFPs were conserved during evolution. Promoter cis-element analysis revealed that most BvZFPs are associated with elements related to phytohormone, biotic or abiotic stress, and plant development. The expression data showed that the BvZFPs in sugar beet are predominantly expressed in the root. In addition, BvZFPs are involved in the response to abiotic and biotic stresses, including salt stress and viral infection. Overall, these results will extend our understanding of the Q-type C2H2 gene family and provide valuable information for the biological breeding of sugar beet against abiotic and biotic stresses in the future.


Introduction
Transcription factors (TFs) function as crucial molecular regulators for gene expression in all organisms [1].In plants, a number of TF families have been identified, including bZIP [2], WRKY [3], NAC [4], MYB [5], and zinc finger proteins (ZFPs) [6], which play important roles in many biological processes, such as growth, development, reproduction, and stress responses.Among these, ZFPs represent one of the largest TF families in plants [6].ZFPs contain a variable number of zinc finger domains, each of which consists of cysteine (Cys) and histidine (His) residues combined with a zinc ion to form a threedimensional finger-type structure.According to the number and arrangement of Cys and His residues, ZFPs have been classified into C2H2, C2HC, C2HC5, CCCH, C3HC4, C4, C4HC3, C6, and C8 [6].
C2H2-type ZFPs, also known as TFIIIA-type zinc fingers, are among the most extensively studied and abundant ZFPs in eukaryotes [7].Recently, the in silico genome-wide identification and functional characterization of plant C2H2-type ZFPs were well analyzed in many species, including 321 members in Glycine max [8], 301 in Brassica rapa [9], 218 in Medicago truncatula [10], 204 in Triticum aestivum [11], 189 in Oryza sativa [12], 176 in Arabidopsis thaliana [13], 150 in Zea mays [14], 145 in Sorghum bicolor [15], 129 in Cucumis sativus [16], 109 in Populus trichocarpa [17], and 104 in Solanum lycopersicum [18].These findings indicate that C2H2 ZFPs are ubiquitous in the plant kingdom, playing important regulatory roles in various biologic processes, including development and organogenesis, as well as responses to stresses and defense.The C2H2-type ZFPs in A. thaliana are categorized into three sets (A, B, and C), with each set further divided into several different subsets, such as C1, C2, and C3 [19].To date, the plant C1 subset has been the most extensively investigated, and its members have been further classified into five subclasses (C1-1i to C1-5i) based on the number of zinc finger domains [20].Moreover, the majority of plant C2H2 ZFPs contain a highly conserved QALGGH sequence within their zinc finger domain (CX 2-4 CX 3 FX 3 QALGGHX 3-5 H), identifying them as the plant-specific Q-type subfamily of C2H2 ZFP.
The first Q-type C2H2 ZFP was discovered in petunia, and a total of 21 Q-type ZFPs were identified based on their specific structures [21,22].The conserved QALGGH motif was found to be essential for DNA binding activity [23].In addition to the conserved zinc finger domain, the N-terminal of some plant C2H2 ZFPs contains a B-box and an L-box motif, while the C-terminal has an ethylene-responsive element binding-factor-associated amphiphilic repression (EAR) motif.The B-box region acts as a nuclear localization signal (NLS).The L-box motif, usually consisting of a core sequence of EXEXXAXCLXXL, is thought to relate to protein-protein interactions.The EAR motif, also known as the DLNbox, has been identified as playing a role in transcriptional repression [24].
The plant Q-type C2H2 ZFP is involved in plant development and various abiotic stress responses, such as drought, salt, osmotic, low temperature, and oxidative stresses [25].For example, the expression of AZF1, AZF2, AZF3, STZ/ZAT10, ZAT11, and ZAT18 is strongly induced by drought, salt, cold stresses, or abscisic acid treatment [26], while the expression levels of ZAT7 and ZAT12 are upregulated by oxidative stress, heat shock, or wounding in A. thaliana [27].ZAT7 has been shown to be involved in plant growth suppression and increased tolerance to salinity stress depending on its EAR motif [28].Overexpressed Zat7 was more tolerant to salinity stress than seedlings of wild-type plants in Arabidopsis [28].Transgenic plants overexpressing Zat12 could tolerate oxidative stress in Arabidopsis [27].More interestingly, ZAT18 was initially identified as a positive regulator of drought stress tolerance in Arabidopsis [29].A recent study indicated that Pseudomonas syringae induces ZAT18 expression to repress the transcription of EDS1 for bacterial infection [30].ZAT18 overexpressing plants were more susceptible to Pst DC3000 compared to Col-0, while ZAT18-KO plants displayed enhanced resistance in Arabidopsis [30].In T. aestivum, 47 Q-type C2H2-ZFPs were identified, and the expression of the majority of TaZFP genes was responsive to drought stress in either leaf or root [31].In O. sativa, OsZFP179 was characterized as a salt-responsive gene, and it was found to enhance salt tolerance in transgenic rice plants [32].Until now, more Q-type C2H2 ZFPs have been identified in A. thaliana, T. aestivum, Brassica oleracea, S. lycopersicum, and Medicago.sativa via genomewide identification [8,20,31,33,34].
Sugar beet (Beta vulgaris), a member of the Amaranthaceae family [35], is an important crop in temperate climates zone, accounting for 20-30% of the world's sugar production.Sugar beet also provides essential raw materials for bioethanol, animal feed, pulp, pectin extract, and functional-food-related industries [36].Sugar beet is frequently subjected to various biotic and abiotic challenges that reduce both sugar yield and processing quality [37].Salinity is a major abiotic stress that limits plant growth and development [38].
Although sugar beet is a salt-tolerant crop, prolonged exposure to salt stress can result in a significant yield loss in beet production [39].In addition, rhizomania caused by beet necrotic yellow vein virus (BNYVV) stands as one of the most severe biotic threats to sustainable beet production globally [40].Susceptible varieties infected by BNYVV exhibit the pronounced lateral rootlet proliferation of taproot and yellow veins on systemically infected leaves [41].This leads to significant losses in root yield as well as a decline in sugar content.Given these challenges, the identification and analysis of sugar beet genes involved in abiotic and biotic stress is crucial, as it offers genetic resources for molecular breeding.Recently, the genome of industrial diploid sugar beet (2n = 18 chromosomes) was sequenced, accelerating sugar beet breeding for tolerances and resistances against abiotic or biotic stresses [35].So far, the structure and function of the BvbZIP and BvWRKY family, as well as their expression pattern under salt stress, have been genome-wide analyzed for sugar beet [42,43].However, information about Q-type C2H2-ZFPs in sugar beet remains unknown.Since Q-type ZFP genes play vital roles in plant development and stress responses, it is essential to identify and analyze the Q-type ZFP gene family in the sugar beet genome.
In this study, we identified 35 Q-type ZFP genes (BvZFPs) in sugar beet comprising different numbers of zinc finger domains.The phylogenetic relationships, genomic location, gene structure, chromosome distribution, gene duplication, and cis-regulatory elements are also explored.In addition, we determined their mRNA expression profiles in both leaf and root tissues under salt stress and BNYVV infection.This study aims to provide a comprehensive understanding of the Q-type C2H2 gene family and shed light on the roles of the BvZFP family in sugar beet, especially the functional characterization of plant development processes and stress responses.This research study can help improve plant quality, stress resistance, and enhanced crop production via genetic modifications in Q-type C2H2 gene family members in sugar beet.

Chromosomal Distribution, Protein Characterization, and Amino Acid Properties
The starting and ending positions of all BvZFP genes on each chromosome were retrieved from the sugar beet gene annotation database, and the results were visualized using the Gene Location Visualize program of Tbtools.Out of 35 genes, 34 BvZFP genes were located on the nine chromosomes, while the remaining 1 may be located on unannotated intergenic regions of the genome.The members of sugar beet BvZFP were renamed according to the location on the chromosome.Protein properties, including the length of the amino acid (aa), molecular masses (MW), theoretical isoelectric point (pI), instability index, and subcellular localization, were predicted using ExPASy Server (https://web.expasy.org/protparam/(accessed on 27 March 2023)).

Multiple Sequence Alignment and Phylogenetic Tree Construction
The full-length sequences of Arabidopsis C2H2 C1-2i zinc finger proteins (Q-type ZFPs) were downloaded from the TAIR database (http://www.arabidopsis.org/(accessed on 16 March 2023)).The multiple sequence alignment of Q-type zinc finger protein members in sugar beet and Arabidopsis was performed using MEGA software with the default parameters.Full-length amino acid sequences were aligned using ClustalX, and the phylogenetic tree was constructed using the neighbor-joining method with the following parameters: bootstrap method with 1000 replicates and partial deletion.

Analysis of Members of the C2H2 ZFP Gene Family
The full-length protein sequences of each subfamily were analyzed for conserved motifs using Multiple Expectation Maximization for Motif Elicitation (MEME, https:// meme-suite.org/meme/(accessed on 22 March 2023)) [46].We used the classic discovery mode and adjusted the parameters as follows: distribution of motifs = any number of repetitions, the maximum number of motifs = 10, and the optimum motif width range from 6 to 50 (inclusive).The gene structure and conserved motifs of C2H2 genes were visualized using Gene Structure View tools in Tbtools software.

Promoter Analysis of the BvZFP Genes in Sugar Beet
The 2000 bp promoter sequences upstream of the 35 BvZFP start codons were extracted from the sugar beet EL10_1.0 genome.The predicted cis-acting elements on promoters with their positional information were identified using the PlantCARE online tool (http: //bioinformatics.psb.ugent.be/webtools/plantcare/html/(accessed on 23 March 2023)).The obtained results were used to predict putative stress-and hormone-responsive cisacting elements and visualized by Simple BioSequence Viewer in TBtools.

Collinearity Analysis of Arabidopsis and Sugar Beet BvZFP Genes
Arabidopsis was selected for collinearity analysis with sugar beet.The genome sequences and annotation files of Arabidopsis were downloaded from NCBI databases, and the chromosome length and location information for the ZFP genes on the genome of the two species were extracted.Using the One Step MCScanX tool of TBtools, we investigated gene replication events and collinearity relationships for gene pairs between two species.All data were visualized via the use of the Advanced Circos program of TBtools software.Ka and Ks substitution between gene pairs were also calculated using the Simple Ka/Ks Calculator tool.

BvZFP Genes Expression Network Analysis
The transcriptome data of sugar beet under salt stress and BNYVV were downloaded from the NCBI SRA database (accession number: PRJNA666117, https://www.ncbi.nlm.nih.gov/bioproject/PRJNA666117/ (accessed on 27 March 2023)), which could be used to analyze the expression pattern of BvZFP genes under 300 mM salt stress.In addition, the database of E-MTAB-8187 from NCBI was used to analyze the expression pattern of BvZFP genes under BNYVV infection.GraphPad Prism was used to map gene expressions.
To identify BvZFP genes between the two groups, the expression level of each transcript was calculated based on the number of fragments per kilobase of exons per million mapped reads (FPKM).RSEM (http://deweylab.biostat.wisc.edu/rsem/(accessed on 1 April 2023) was used to determine gene abundance.The R statistical package software edgeR (empirical analysis of digital gene expression in R, http://www.bioconductor.org/packages/2.12/bioc/html/edgeR.html (accessed on 1 April 2023)) was used for differential expression analysis.

Identification and Chromosomal Localization of the C2H2 Q-Type ZFP Subclass in Sugar Beet
To identify C2H2 ZFP in the B. vulgaris genome, we utilized hidden Markov model (HMM) files (PF13912, PF12756, and PF00096) to conduct a genome-wide HMM search.This search yielded 104 putative non-redundant sugar beet C2H2 ZFP proteins (BvZFPs).This number was greater than that present in PlantTFDB, where 64 BvZFPs have been deposited for sugar beet (http://planttfdb.gao-lab.org/family.php?sp=Bvu&fam=C2H2 (accessed on 17 March 2023)).Subsequently, 35 Q-type ZFPs, characterized by the specific CX 2 CX 3 FX 3 QALGGHX 3-5 HX domain, were manually selected (Supplementary File S1: Tables S1 and S2).According to the sugar beet genome database [47], we generated a map that detailed the physical positions of the Q-type BvZFPs (Figure 1).These Q-type BvZFPs were then renamed from BvZFP1 to BvZFP35 based on their physical positions on sugar beet chromosomes (Chr).All BvZFP genes were distributed widely and unevenly on nine Chrs, except for BvZFP35, which was located on unmapped scaffolds (Figure 1).Chr6 contained the largest number of BvZFP members and contained 8 BvZFPs, followed by Chr2, Chr3, Chr9, Chr1, Chr5, and Chr7, each of which contained 3 to 5 BvZFPs.Chr4 and Chr8 contained relatively fewer BvZFP members, with only two and one gene, respectively.

Identification and Chromosomal Localization of the C2H2 Q−Type ZFP Subclass in Sugar Beet
To identify C2H2 ZFP in the B. vulgaris genome, we utilized hidden Markov model (HMM) files (PF13912, PF12756, and PF00096) to conduct a genome-wide HMM search.This search yielded 104 putative non-redundant sugar beet C2H2 ZFP proteins (BvZFPs).This number was greater than that present in PlantTFDB, where 64 BvZFPs have been deposited for sugar beet (http://planttfdb.gao-lab.org/family.php?sp=Bvu&fam=C2H2 (accessed on 17 March 2023)).Subsequently, 35 Q−type ZFPs, characterized by the specific CX2CX3FX3QALGGHX3-5HX domain, were manually selected (Additional file 1: Tables S1  and S2).According to the sugar beet genome database [47], we generated a map that detailed the physical positions of the Q−type BvZFPs (Figure 1).These Q−type BvZFPs were then renamed from BvZFP1 to BvZFP35 based on their physical positions on sugar beet chromosomes (Chr).All BvZFP genes were distributed widely and unevenly on nine Chrs, except for BvZFP35, which was located on unmapped scaffolds (Figure 1).Chr6 contained the largest number of BvZFP members and contained 8 BvZFPs, followed by Chr2, Chr3, Chr9, Chr1, Chr5, and Chr7, each of which contained 3 to 5 BvZFPs.Chr4 and Chr8 contained relatively fewer BvZFP members, with only two and one gene, respectively.

Characterization of the Sugar Beet C2H2 Q−Type ZFP Subclass
Based on the number of zinc finger domains and the spacing between the two His residues, the 35 BvZFP genes were divided into four groups, including 22 members in the C1-1i group (one zinc finger domain), 6 members in the C1-2i group, 6 members in the C1-3i group, and 1 member in the C1-4i group (Additional file 1: Table S3).Interestingly, several zinc finger domains with certain modifications of the "QALGGH" motif were observed for all 3i and 4i members (Additional file 1: Table S3).According to a previous study [16], the modified zinc finger domains were classified as the M-type.Furthermore, the protein properties of these genes were predicted using the ExPASy Server, and the results are shown in Table 1 and Additional file 1: Table S4.The protein lengths of the 1i group, 2i group, and 3i group ranged from 162 to 323, 186 to 456, and 237 to 572 amino acids for 1i group, 2i group, and 3i group, respectively.The molecular weights of the 1i

Characterization of the Sugar Beet C2H2 Q-Type ZFP Subclass
Based on the number of zinc finger domains and the spacing between the two His residues, the 35 BvZFP genes were divided into four groups, including 22 members in the C1-1i group (one zinc finger domain), 6 members in the C1-2i group, 6 members in the C1-3i group, and 1 member in the C1-4i group (Supplementary File S1: Table S3).Interestingly, several zinc finger domains with certain modifications of the "QALGGH" motif were observed for all 3i and 4i members (Supplementary File S1: Table S3).According to a previous study [16], the modified zinc finger domains were classified as the M-type.Furthermore, the protein properties of these genes were predicted using the ExPASy Server, and the results are shown in Table 1 and Supplementary File S1: Table S4.The protein lengths of the 1i group, 2i group, and 3i group ranged from 162 to 323, 186 to 456, and 237 to 572 amino acids for 1i group, 2i group, and 3i group, respectively.The molecular weights of the 1i group's proteins ranged from 18,233 to 35,834 Da, with an average weight of 25,264 Da.For the 2i group, the molecular protein weights ranged from 20,667 to 50,063 Da, with an average weight of 32,905 Da.The average molecular weight of 3i group proteins is 50,602 Da, with individual weights ranging from 25,331 to 62,896 Da.For all BvZFPs, the theoretical isoelectric point (pI) fell between 5.69 and 9.22, and the instability index varied from 37.2 to 75.14.The GRAVY values, ranging from −1.27 to −0.358, revealed that sugar beet Q-type BvZFPs were hydrophilic proteins.In addition, subcellular localization prediction suggested that all BvZFPs are located in the nucleus.

Phylogenetic Analysis of Q-Type BvZFP Genes
To explore the evolutionary relationship between C2H2 Q-type genes in sugar beet and A. thaliana (AtZFP), we constructed a phylogenetic tree using MEGA10 based on the alignment of 93 Q-type ZFPs amino acid sequences at the whole protein level, which included 35 from sugar beet and 58 from Arabidopsis.The resulting tree classified the members into six major clades, including C1-1Q-A, C1-1Q-B, C1-2Q, C1-QM-A, C1-QM-B, and C1-QM-C (Figure 2).Remarkably, ZFP proteins that possessed the same types and numbers of zinc finger domains were clustered into the same clade.In the 1i group, 15 BvZFPs together with 16 AtZFPs were grouped in the C1-1Q-A clade, while 7 BvZFPs with 11 AtZFPs were grouped into C1-1Q-B clades.Both C1-1QA and C1-1QB clades contained one conserved "QALGGH" motif in their protein sequences.For the members of the 2i group, five BvZFPs along with 11 AtZFPs belonged to the C1-2Q clade.Following the same nomenclature as the C1-1Q clade, the members in the C1-2Q clade contained two conserved "QALGGH" motifs.The rest of the BvZFP members, which presented a variety of combinations with different numbers of Q-and M-type zinc finger domains, were defined as the C1-QM clades.For 3i and 4i group members, four BvZFPs were clustered into C1-QM-A clade together with nine AtZFPs, while four BvZFPs and seven AtZFPs were grouped into the clade C1-QM-B (Figure 2).Moreover, an unrooted phylogenetic tree comprising 35 BvZFPs was constructed.Based on sequence similarity and tree topology, these BvZFPs were further classified into three clades, C1-1Q, C1-2Q, and C1-QM, and five subclades, including C1-1Q-A, C1-1Q-B, C1-2Q, C1-QM-A, and C1-QM-B (Figure 3A), according to the classification in Figure 2.

Gene Structure and Conservative Motif Analysis of C2H2 Q−Type BvZFPs
We analyzed the number of exon and intron structures of all 35 Q−type BvZFP As shown in Figure 3, 30 BvZFPs were intronless, 4 BvZFPs contained a single intro only 1 BvZFP had four introns.All members with introns except for BvZFP17 belon the 1i group.To better understand the characteristic regions of BvZFP proteins, w the MEME online website to identify predicted conserved motifs.A total of 10 con motifs were identified, with their details presented in Additional file 1: Table S5.M and 2 were recognized as Q−type zinc finger domains (Figure 3D,F), while motif 4 sented the M−type zinc finger domains (Supplement Figure S3A).Motif 1 was wid tributed in all 35 BvZFP proteins (Figure 3C), and motif 2 was found in members C1-2Q, C1-QM-A, and C1-QM-B clades (Figure 3D).Motif 4 was specific to the C clades.Motif 3, known as the EAR motif (Figure 3E), had conserved amino acid sign "LxLxL" or "DLNx(1-2)P".Among these BvZFPs, 18 BvZFPs contained the "LxLxL of EAR motif, while 14 BvZFPs contained the "DLNx(1-2)P" type of EAR motif.I tion, BvZFPs EAR motifs were predominantly located at the C-terminus region, exc BvZFP21 and BvZFP14.Motif 5, termed the L-BOX motif, was characterized by sequence of EXEXXAXCLXXL (Figure 3G).All 14 members in the C1-2Q and C

Gene Structure and Conservative Motif Analysis of C2H2 Q-Type BvZFPs
We analyzed the number of exon and intron structures of all 35 Q-type BvZFP genes.As shown in Figure 3, 30 BvZFPs were intronless, 4 BvZFPs contained a single intron, and only 1 BvZFP had four introns.All members with introns except for BvZFP17 belonged to the 1i group.To better understand the characteristic regions of BvZFP proteins, we used the MEME online website to identify predicted conserved motifs.A total of 10 conserved motifs were identified, with their details presented in Supplementary File S1: Table S5.Motifs 1 and 2 were recognized as Q-type zinc finger domains (Figure 3D,F), while motif 4 represented the M-type zinc finger domains (Supplement Figure S1A).Motif 1 was widely distributed in all 35 BvZFP proteins (Figure 3C), and motif 2 was found in members of the C1-2Q, C1-QM-A, and C1-QM-B clades (Figure 3D).Motif 4 was specific to the C1-QM clades.Motif 3, known as the EAR motif (Figure 3E), had conserved amino acid signatures "LxLxL" or "DLNx(1-2)P".Among these BvZFPs, 18 BvZFPs contained the "LxLxL" type of EAR motif, while 14 BvZFPs contained the "DLNx(1-2)P" type of EAR motif.In addition, BvZFPs EAR motifs were predominantly located at the C-terminus region, except for BvZFP21 and BvZFP14.Motif 5, termed the L-BOX motif, was characterized by a core sequence of EXEXXAXCLXXL (Figure 3G).All 14 members in the C1-2Q and C1-QM clades contained Motif 5. Notably, some motifs were clade-specific.For example, motifs 10 and 6 mainly existed in C1-2Q.Motif 9 appeared in members of C1-2Q, C1-QM, and C1-1Q-B.Motif 8 was unique to C1-1Q-A (Figure 3).In summary, it was found that the BvZFPs clustered in the same subclades shared a similar motif composition, suggesting functional similarities among these evolutionarily conserved BvZFPs within the same clades.

Genomic Collinearity Analysis of C2H2 Q-Type BvZFPs between Sugar Beet and Arabidopsis
To investigate the potential evolution processes of C2H2 Q-type BvZFP genes, we analyzed the synteny relationship of 34 BvZFPs using TBtools.Intraspecific collinearity analysis showed that three pairs of BvZFPs exhibited collinearity within sugar beet.The ratio of the nonsynonymous substitution rate to the synonymous substitution rate (Ka/Ks) is widely used for evaluating the selective pressure of duplication events [34].Generally, Ka/Ks > 1 indicates positive selection, Ka/Ks = 1 indicates neutral selection, and Ka/Ks < 1 indicates purification [48].The Ka/Ks value of BvZAT7 and BvZAT32 is 0.1419 (Figure 4, Supplementary File S1: Table S10-1).Moreover, the phylogenetic tree analysis showed that these two genes were on the same branch.The results indicated that these synteny genes have undergone purifying selection in their evolutionary history.To further explore the evolutionary relationship of C2H2 Q-type BvZFPs, interspecific synteny comparisons between B. vulgaris and A. thaliana were also made.In our study, we identified 24 pairs of ZFP genes in sugar beet and Arabidopsis (Figure 5).The Ka/Ks ratios of these gene pairs ranged from 0.1 to 0.4, with an average of 0.25 (Supplementary File S1: Table S10-2), suggesting that the C2H2 Q-type ZFP gene remained relatively conserved in different species.
functional similarities among these evolutionarily conserved BvZFPs within the same clades.

Genomic Collinearity Analysis of C2H2 Q−Type BvZFPs between Sugar Beet and Arabidopsis
To investigate the potential evolution processes of C2H2 Q−type BvZFP genes, we analyzed the synteny relationship of 34 BvZFPs using TBtools.Intraspecific collinearity analysis showed that three pairs of BvZFPs exhibited collinearity within sugar beet.The ratio of the nonsynonymous substitution rate to the synonymous substitution rate (Ka/Ks) is widely used for evaluating the selective pressure of duplication events [34].Generally, Ka/Ks > 1 indicates positive selection, Ka/Ks = 1 indicates neutral selection, and Ka/Ks < 1 indicates purification [48].The Ka/Ks value of BvZAT7 and BvZAT32 is 0.1419 (Figure 4, Additional file 1: Table S10-1).Moreover, the phylogenetic tree analysis showed that these two genes were on the same branch.The results indicated that these synteny genes have undergone purifying selection in their evolutionary history.To further explore the evolutionary relationship of C2H2 Q−type BvZFPs, interspecific synteny comparisons between B. vulgaris and A. thaliana were also made.In our study, we identified 24 pairs of ZFP genes in sugar beet and Arabidopsis (Figure 5).The Ka/Ks ratios of these gene pairs ranged from 0.1 to 0.4, with an average of 0.25 (Additional file 1: Table S10-2), suggesting that the C2H2 Q−type ZFP gene remained relatively conserved in different species.

Promoter Analysis of the Q−Type BVZFP Genes in Sugar Beet
Cis-acting regulatory elements play essential roles in modulating the plant response to biotic and abiotic stresses.Therefore, the 2000 bp promoter regions located upstream of the BvZFP genes were extracted from the sugar beet EL10_1.0 genome database [44].These sequences were then analyzed for cis-acting elements using the PlantCARE website.The results showed that 385 cis-regulatory elements belonging to 26 categories were obtained (Additional file 1: Table S6).The basal promoter elements, such as TATA-box and CAATbox, and unannotated elements were excluded from this count.Most cis-elements were related to hormonal response and stress signal responsiveness (Additional file 1: Table S7) Among them, 197 elements were found to be related to hormonal responses, including 63 abscisic-acid-responsive elements (ABRE), 62 MeJA-responsive elements (TGACG-motif) 30 salicylic-acid-responsive elements (TCA-element), 29 gibberellin-responsive elements (TATC-box), and 13 auxin-responsive elements (TGA-element).Meanwhile, 46 elements associated with plant development, 47 elements that respond to abiotic stress, and 22 elements related to biotic stress were identified.In addition, the abiotic-stress-responsive elements were related to drought and low-temperature responses.As shown in Figure 6, 22 BvZFPs contain more than 10 functional elements from different categories.We also found that BvZFP8, BvZFP34, BvZFP24, and BvZFP17 had a higher number of cis-regulatory elements, while BvZAT2 and BvZAT26 contained fewer cis-elements (Figure 6).The presence of these cis-regulatory elements in the BvZFP promoter regions suggested their involvement in the regulation of various plant pathways.

Promoter Analysis of the Q-Type BVZFP Genes in Sugar Beet
Cis-acting regulatory elements play essential roles in modulating the plant response to biotic and abiotic stresses.Therefore, the 2000 bp promoter regions located upstream of the BvZFP genes were extracted from the sugar beet EL10_1.0 genome database [44].These sequences were then analyzed for cis-acting elements using the PlantCARE website.The results showed that 385 cis-regulatory elements belonging to 26 categories were obtained (Supplementary File S1: Table S6).The basal promoter elements, such as TATAbox and CAAT-box, and unannotated elements were excluded from this count.Most cis-elements were related to hormonal response and stress signal responsiveness (Supplementary File S1: Table S7).Among them, 197 elements were found to be related to hormonal responses, including 63 abscisic-acid-responsive elements (ABRE), 62 MeJAresponsive elements (TGACG-motif), 30 salicylic-acid-responsive elements (TCA-element), 29 gibberellin-responsive elements (TATC-box), and 13 auxin-responsive elements (TGAelement).Meanwhile, 46 elements associated with plant development, 47 elements that respond to abiotic stress, and 22 elements related to biotic stress were identified.In addition, the abiotic-stress-responsive elements were related to drought and low-temperature responses.As shown in Figure 6, 22 BvZFPs contain more than 10 functional elements from different categories.We also found that BvZFP8, BvZFP34, BvZFP24, and BvZFP17 had a higher number of cis-regulatory elements, while BvZAT2 and BvZAT26 contained fewer cis-elements (Figure 6).The presence of these cis-regulatory elements in the BvZFP promoter regions suggested their involvement in the regulation of various plant pathways.The phylogenetic tree was developed on MEGA 10 using neighbor-joining phylogenetic method analysis.Both the bootstrap test and the approximate likelihood ratio test were set to 1000 times.

Expression Profiles Analysis of Q−Type BvZFPs in Different Tissue
To investigate the expression patterns of Q−type BvZFP genes in different tissues, we analyzed sugar beet transcriptome data for leaf and root tissues from three-pair-euphyllastage sugar beet, which were taken from the SRA database [49].A heatmap was constructed to visualize the expression patterns via TBtools software (Figure 7, Additional file 1: Table S8).In the control group, some BvZFPs displayed a tissue-specific expression pattern.Specifically, BvZFP3, BvZFP8, BvZFP9, BvZFP16, BvZFP18, BvZFP19, BvZFP32, and BvZFP35 were only expressed in roots.Meanwhile, BvZFP2, BvZFP6, BvZFP7, BvZFP14, BvZFP17, BvZFP20, BvZFP27, BvZFP28, BvZFP30, and BvZFP34 were expressed in both leaf and root tissues, while 17 BvZFPs did not show expression in either tissues.Interestingly, the expression levels of all BvZFPs in root tissues were generally higher than that in leaf tissue, except for BvZFP14.The phylogenetic tree was developed on MEGA 10 using neighbor-joining phylogenetic method analysis.Both the bootstrap test and the approximate likelihood ratio test were set to 1000 times.

Expression Profiles Analysis of Q-Type BvZFPs in Different Tissue
To investigate the expression patterns of Q-type BvZFP genes in different tissues, we analyzed sugar beet transcriptome data for leaf and root tissues from three-pair-euphyllastage sugar beet, which were taken from the SRA database [49].A heatmap was constructed to visualize the expression patterns via TBtools software (Figure 7, Supplementary File S1: Table S8).In the control group, some BvZFPs displayed a tissue-specific expression pattern.Specifically, BvZFP3, BvZFP8, BvZFP9, BvZFP16, BvZFP18, BvZFP19, BvZFP32, and BvZFP35 were only expressed in roots.Meanwhile, BvZFP2, BvZFP6, BvZFP7, BvZFP14, BvZFP17, BvZFP20, BvZFP27, BvZFP28, BvZFP30, and BvZFP34 were expressed in both leaf and root tissues, while 17 BvZFPs did not show expression in either tissues.Interestingly, the expression levels of all BvZFPs in root tissues were generally higher than that in leaf tissue, except for BvZFP14.

Responses of Q-Type BvZFP Genes under Salt Treatment and Viral Infection
To further provide insight into the response of BvZFPs to abiotic and biotic stress, transcriptome data from the SRA database [49] were used to investigate the expression patterns of BvZFPs under salt stress or viral infection (Supplementary File S1: Tables S8 and S9).Under a 300 mM NaCl treatment, the expression of most BvZFPs was changed.At 72 h after treatment, BvZFP2, BvZFP6, BvZFP17, BvZFP30, and BvZFP34 were significantly up-regulated in leaves, while BvZFP2, BvZFP14, BvZFP16, and BvZFP34 were significantly up-regulated in roots (Figure 7).Comparing the expression patterns between root and leaf tissues, we found that BvZFP2 and BvZFP34 were up-regulated in both tissues, suggesting their important role in response to salt stress in these tissues.

Responses of Q−Type BvZFP Genes under Salt Treatment and Viral Infection
To further provide insight into the response of BvZFPs to abiotic and biotic stress, transcriptome data from the SRA database [49] were used to investigate the expression patterns of BvZFPs under salt stress or viral infection (Additional file 1: Table S8 and S9).Under a 300 mM NaCl treatment, the expression of most BvZFPs was changed.At 72 h after treatment, BvZFP2, BvZFP6, BvZFP17, BvZFP30, and BvZFP34 were significantly upregulated in leaves, while BvZFP2, BvZFP14, BvZFP16, and BvZFP34 were significantly up-regulated in roots (Figure 7).Comparing the expression patterns between root and leaf tissues, we found that BvZFP2 and BvZFP34 were up-regulated in both tissues, suggesting their important role in response to salt stress in these tissues.
Transcription factors in plants also play a pivotal role in the response to pathogen infections.Using beet necrotic yellow vein virus-infected sugar beet transcriptome data from the database [50], It was revealed that nearly half of the BvZFP genes exhibited differential expression at varying degrees.Among these genes, BvZFP6, BvZFP8, BvZFP28, and BvZFP34 were significantly up-regulated, while BvZFP2, BvZFP16, BvZFP17, and BvZFP30 were moderately induced.Conversely, the expression levels of BvZFP13, BvZFP25, BvZFP26, and BvZFP33 were down-regulated (Figure 8).The above results indicated that certain BvZFP genes may be involved in either promoting or inhibiting viral infection by regulating downstream target genes.Transcription factors in plants also play a pivotal role in the response to pathogen infections.Using beet necrotic yellow vein virus-infected sugar beet transcriptome data from the database [50], It was revealed that nearly half of the BvZFP genes exhibited differential expression at varying degrees.Among these genes, BvZFP6, BvZFP8, BvZFP28, and BvZFP34 were significantly up-regulated, while BvZFP2, BvZFP16, BvZFP17, and BvZFP30 were moderately induced.Conversely, the expression levels of BvZFP13, BvZFP25, BvZFP26, and BvZFP33 were down-regulated (Figure 8).The above results indicated that certain BvZFP genes may be involved in either promoting or inhibiting viral infection by regulating downstream target genes.

Discussion
C2H2 ZFPs are one of the most extensively studied transcription factors that play crucial roles in many biological processes in eukaryotic organisms [6,21].The Q−type ZFP, a plant-specific subfamily of C2H2−ZFPs, was found to be involved in plant development, as well as various stress responses [15,34,51].So far, Q−type ZFPs have been extensively studied in many species, such as A. thaliana, S. lycopersicum, Solanum tuberosum, T. aestivum, O. sativa, and P. trichocarpa (Table 2).However, this subfamily has yet to be explored in B. vulgaris, one of the most important crops for sugar production.In this study, we performed a comprehensive genome-wide investigation of Q−type C2H2 ZFPs in sugar beet.

Discussion
C2H2 ZFPs are one of the most extensively studied transcription factors that play crucial roles in many biological processes in eukaryotic organisms [6,21].The Q-type ZFP, a plant-specific subfamily of C2H2-ZFPs, was found to be involved in plant development, as well as various stress responses [15,34,51].So far, Q-type ZFPs have been extensively studied in many species, such as A. thaliana, S. lycopersicum, Solanum tuberosum, T. aestivum, O. sativa, and P. trichocarpa (Table 2).However, this subfamily has yet to be explored in B. vulgaris, one of the most important crops for sugar production.In this study, we performed a comprehensive genome-wide investigation of Q-type C2H2 ZFPs in sugar beet.A total of 104 C2H2 BvZFPs were first identified from the current B. vulgaris genome using bioinformatic analysis.This number of BvZFP was lower than most previously studied species, except for Vitis vinifera and C. annuum (Table 2).Among these genes, 35 Q-type C2H2 BvZFPs, with each containing at least one conserved "QALGGH"-type zinc finger domain, were selected for further investigation.As expected, sugar beet also has the fewest Q-type ZFP genes compared with other sequenced plant genomes (Table 2).These Q-type ZFP genes were renamed from BvZFP1 to BvZFP35 based on their physical positions on the sugar beet chromosomes.Thirty-four members were unevenly distributed in all nine chromosomes, except for BvZFP35 (LOC104893654) (Figure 1).Based on the number of zinc finger domains and the spacing between the two His residues, 63% (22 out of 35) members belonged to C1-1i group, 34% of members existed in both C1-2i and C1-3i groups, and only one gene BvZFP22 was categorized in the C1-4i group.We did not identify any Q-type ZFPs in sugar beet that possessed five zinc finger domains, which are present in Arabidopsis and rice genomes [12,20].Interestingly, we found several zinc finger domains with certain modifications to the "QALGGH" motif (M-type) for all 3i and 4i members.The lengths of the BvZFP proteins differed among groups, ranging from 162 to 323, 186 to 456, and 237 to 572 amino acids for 1i group, 2i group, and 3i group, respectively.All Q-type BvZFPs were predicted to be localized within the nucleus, indicating that these proteins indeed function as transcription factors in the nucleus.Phylogenetic analyses of 35 BvZFPs, using 58 Arabidopsis ZFPs as templates, classified 35 Q-type BvZFP genes into five major clades according to the number of "KS/KA/RS/RA/QA-LGGH" motifs they contained, including C1-1Q-A, C1-1Q-B, C1-2Q, C1-QM-A, and C1-QM-B (Figure 2).Proteins clustered into the same clade may have similar functions under various stresses and close evolutionary relationships.Among them, C1-2Q clusters with Arabidopsis C1-2i, which is highly conserved during evolution.It has been reported that mutations in the QALGGH sequence greatly affect DNA-binding activity [59].Cis-acting elements are involved in the regulation of gene activity and serve as fundamental molecular switches during transcriptional regulation [60].A previous study confirmed that many of these elements, such as ABREs and DREs, have been reported to widely participate in abiotic stress responses in Artemisia annua, A. thaliana, and T. aestivum [31].In plants, hormones, such as auxin, abscisic acid, and gibberellin, play an important role in the growth and development of plants and the response to adversity stress [61,62].However, no study has investigated these key regulatory elements in sugar beet.Among the BvZFP gene family members, 385 cis-regulatory elements were identified, which contain elements related to hormone response (abscisic acid, gibberellin, ethylene, and auxin), stress response elements (low temperature and drought), and growth-and development-related response elements.These results spotlight potential candidate genes for anti-abiotic stress, although their specific functions need to be confirmed by further investigation.Previous studies have shown that genes with fewer introns are more prone to be activated in response to stress [63].In this study, approximately 86% (30/35) of BvZFPs have no introns, except for five BvZFPs, which indicates that BvZFPs could rapidly react to external stimuli.Collinearity analysis, an essential analytical strategy in comparative genomics, illuminates both large-scale and small-scale molecular evolutionary events across species.In this study, about 70% of the 34 mapped BvZFP genes are collinear with respect to four C2H2-ZFP genes on Arabidopsis chromosomes.The Ka/Ks ratios showed the conservation of BvZFPs throughout evolution, aligning with the previous views that regarded the C2H2 family as evolutionarily stable [64,65].All this detailed information helps us better understand and screen for appropriate C2H2 gene family members in sugar beet.
Besides the conserved zinc finger domains, many previously identified plant ZFP proteins have EAR motifs.These motifs, characterized by conserved amino acid signatures "LxLxL" or "DLNxxP" in their respective C-terminal regions [66], play a key role in transcriptional repression [67,68].It has been reported that 26% of wheat TaZFP members contained potential EAR motifs [31], while approximately 55% of grape VviZFP members possessed this motif [56].Genome-wide analysis suggests that the EAR motif is conserved across all plant species and is involved in developmental and stress-related processes [69].For instance, an EAR repressor named NIMIN1 negatively regulates the expression of the PR1 defense gene in Arabidopsis [70].AtZP1, which contains an EAR motif, negatively regulates Arabidopsis root hair initiation and elongation [71].Most strikingly, our results reveal that 32 BvZFPs have at least one EAR motif, accounting for 91% of BvZFP members.As a key repression domain, alterations to any individual residue within the EAR can reduce or completely eliminate its transcriptional repression capability [72].
Previous studies demonstrated that many plant Q-type ZFP genes exhibited different expression patterns across various tissues [73][74][75][76].For example, the Q-type genes of cucumber could be clustered into four groups according to their expression levels; CsZFP genes in group 1 showed high expression levels, while CsZFPs in group 4 exhibited little to no expression in all tissues [16].In potatoes, some Q-type StZFP genes displayed a tissuespecific expression pattern.Specifically, PG0005486 and PG0030311 were predominantly expressed in leaves and roots, respectively [58].In strawberries, FaZAT10 was highly expressed in roots, followed by leaves and stems [74].Moreover, RNA-seq data analysis for Q-type TaZFPs in different wheat tissues showed that 75% of TaZFPs were predominantly expressed in roots, implying their potential role in regulating wheat root development [31].
In the current study, we also found that many BvZFPs displayed a specific expression pattern between leaves and roots (Figure 6).With respect to the three-pair-euphylla-stage seedlings of sugar beet, 46% BvZFPs (16 out of 35) were not observed to be expressed in either leaf or root tissue, suggesting that these genes might not influence development at this growth stage.In contrast, 23% BvZFPs (8 out of 35 genes) were expressed preferentially in roots while 29% BvZFPs (10 out of 35) were expressed in both leaf and root tissues to variable degrees.Among the genes expressed in both tissues, root expression levels for BvZFPs were relatively higher than in leaf tissue except for one gene.In addition, BvZFP6, BvZFP27, and BvZFP30, each containing two "QALGGH" domains, were significantly expressed in roots.Therefore, the results suggested that many Q-type BvZFP genes are involved in sugar beet development under normal growth conditions, especially root development.
Plants are frequently affected by various abiotic and biotic.To adapt to these stresses, plants have evolved sophisticated mechanisms.Numerous Q-type ZFP genes have been found in response to diverse stresses, such as drought, salt, osmotic, low temperature, oxidative stress, and pathogen infection in many plant species [16,30,32].Studies have shown that the overexpression of MdZAT5 promotes the expression of anthocyanin-biosynthesisrelated genes NHX1 and ABI1 to actively regulate anthocyanin synthesis and increases sensitivity to salt stress in apple calli and Arabidopsis [75].The ZFP gene Csa6G303740 was significantly up-regulated in response to drought, heat, and salt stresses, indicating it could be involved in the regulation of abiotic stresses in cucumbers [16].Potato StZFP1, which responds to salt, dehydration, and infection by Phytophthora infestans, has the potential to improve salt tolerance in transgenic tobacco [76].
In this work, we analyzed publicly available transcriptome data and found that five BvZFP genes, BvZFP2, BvZFP6, BvZFP17, BvZFP30, and BvZFP34, were up-regulated in leaves or/and roots under salt treatment, indicating that these genes may play roles in salt response in sugar beet.These genes belong to the C1-2i and C1-3i groups, and all contain two Q-type zinc finger domains.In addition, BvZFP2, BvZFP17, and BvZFP 30 are grouped into the same C1-QM clade, while BvZFP6 and BvZFP34 are clustered into the C1-2Q clade.For biotic stress, twelve BvZFP genes were differentially expressed by BNYVV infection.Interestingly, BvZFP2, BvZFP6, BvZFP17, BvZFP30, and BvZFP34 were induced either by salt treatment or viral infection.On the other hand, some BvZFP genes only responded to viral infection; for example, BvZFP8 and BvZFP28 were significantly up-regulated in the BNYVV-infected samples, and BvZFP13, BvZFP25, BvZFP26, and BvZFP33 were downregulated exclusively due to viral infection.These four down-regulated genes all belong to the C1-1i group and are clustered into the C1-1Q-A clade.BvZFP8 and BvZFP28 are in the C1-1Q-B and C1-2Q clades, respectively.
Generally, members in the same cluster may have conserved functions.Previous studies have shown that Arabidopsis ZAT6 is strongly induced by various stresses, including salt, cold, dehydration, and pathogen infection [77].Plants overexpressing AtZAT6 improved resistance to salt, drought, freezing stresses, and pathogen infection.In our study, the BvZFP30 clustered with AtZAT6 in C1-QM was differentially expressed under salt in sugar beet, indicating that BvZFP30 may have functions that are similar to AtZAT6 under abiotic stress.ZAT18 was initially shown to be a positive regulator of drought stress tolerance in Arabidopsis [29].A recent study has indicated that Pseudomonas syringae induces ZAT18 expression to repress EDS1 transcription during bacterial infection [30].Furthermore, BvZFP6, BvZFP34, and BvZFP28, which are clustered with ZAT18, are induced by BNYVV infection.Given that both BvZFPs and ZAT18 are highly expressed under the attack of pathogens, it is plausible that these genes may have similar functions with respect to responding to biotic stress.Taken together, these results indicate that BvZFP genes are involved in plant development and stress resistance in sugar beet, particularly in response to both abiotic and biotic stresses.

Conclusions
In summary, we performed a genome-wide exploration of Q-type ZFPs in sugar beet, and a total of 35 Q-type BvZFP genes that contain at least one conserved "QALGGH" motif were identified.We analyzed the physiochemical properties, genomic location, gene duplication event, gene structure, conserved motif compositions, and cis-regulatory elements of these BvZFP genes.Phylogenetic analysis revealed that 35 Q-type BvZFP genes were classified into five subclades.Moreover, the BvZFPs clustered within the same clade generally shared a similar motif composition.In addition, the expression profiles of BvZFPs in leaf and root tissues, as well as their responses to salt stress and viral infection, were analyzed.Q-type BvZFPs are predominantly expressed in roots and are enriched with members that are responsive to abiotic and biotic stresses.These results will enhance our understanding of the BvZFP gene family and provide valuable information for the further functional analysis of Q-type BvZFPs relative to the abiotic and biotic stress tolerance of sugar beet.Meanwhile, this study provides a theoretical basis for the biological breeding of sugar beet against salt stress and viral infection in the future.

Figure 1 .
Figure 1.Chromosomal localization of Q−type BvZFP genes on sugar beet chromosomes.The 34 members are distributed over 9 chromosomes.

Figure 1 .
Figure 1.Chromosomal localization of Q-type BvZFP genes on sugar beet chromosomes.The 34 members are distributed over 9 chromosomes.

Biology 2023 ,
12, x FOR PEER REVIEW 7 of 20 numbers of zinc finger domains were clustered into the same clade.In the 1i group, 15 BvZFPs together with 16 AtZFPs were grouped in the C1-1Q-A clade, while 7 BvZFPs with 11 AtZFPs were grouped into C1-1Q-B clades.Both C1-1QA and C1-1QB clades contained one conserved "QALGGH" motif in their protein sequences.For the members of the 2i group, five BvZFPs along with 11 AtZFPs belonged to the C1-2Q clade.Following the same nomenclature as the C1-1Q clade, the members in the C1-2Q clade contained two conserved "QALGGH" motifs.The rest of the BvZFP members, which presented a variety of combinations with different numbers of Q-and M-type zinc finger domains, were defined as the C1-QM clades.For 3i and 4i group members, four BvZFPs were clustered into C1-QM-A clade together with nine AtZFPs, while four BvZFPs and seven AtZFPs were grouped into the clade C1-QM-B (Figure 2).

Figure 2 .
Figure 2. Phylogenetic relationships of Q−type C2H2−ZFPs between sugar beet and Arabidopsis.The unrooted phylogenetic tree was constructed using MEGA10 via the neighbor-joining method with 1000 bootstrap replicates.The tree was divided into six phylogenetic clusters.The red circles represent Arabidopsis ZFPs, the blue circles represent sugar beet ZFPs.

Figure 3 .
Figure 3. Motif distributions and gene-structure analysis of Q−type BvZFP genes.(A) The p netic tree was built using the NJ method with a bootstrap value of 1000.(B) Exons, introns, untranslated region (UTR) are represented by yellow rectangles, gray lines, and green rec respectively.(C) The conserved motifs in BvZFP proteins (1-10) are shown in different col gray lines represent relative protein lengths.(D-G) The common conserved motifs in BvZ teins.

Figure 3 .
Figure 3. Motif distributions and gene-structure analysis of Q-type BvZFP genes.(A) The phylogenetic tree was built using the NJ method with a bootstrap value of 1000.(B) Exons, introns, and the untranslated region (UTR) are represented by yellow rectangles, gray lines, and green rectangles, respectively.(C) The conserved motifs in BvZFP proteins (1-10) are shown in different colors.The gray lines represent relative protein lengths.(D-G) The common conserved motifs in BvZFP proteins.

Figure 4 .
Figure 4. Distribution and collinearity of Q−type BvZFP genes in sugar beet.The outer circle represents the location of BvZFP genes on the chromosomes, and the inner circle histogram and heat map represent gene density.The grey lines represent all collinear genes on the sugar beet genome, and the colored lines connect the collinear BvZFP gene pairs.

Figure 4 .
Figure 4. Distribution and collinearity of Q-type BvZFP genes in sugar beet.The outer circle represents the location of BvZFP genes on the chromosomes, and the inner circle histogram and heat map represent gene density.The grey lines represent all collinear genes on the sugar beet genome, and the colored lines connect the collinear BvZFP gene pairs.

Figure 5 .
Figure 5. Collinearity analysis of Q−type C2H2−ZFPs in sugar beet and Arabidopsis genomes.The outer circle represents the different chromosomes of the two species, and the inner circle histogram and heat map represent the gene density of each species.The grey curves indicate the collinear gene regions within the genomes of the two species, and the colored curves emphasize the specific collinearity relationships among Q−type BvZFP genes between the two species.

Figure 5 .
Figure 5. Collinearity analysis of Q-type C2H2-ZFPs in sugar beet and Arabidopsis genomes.The outer circle represents the different chromosomes of the two species, and the inner circle histogram and heat map represent the gene density of each species.The grey curves indicate the collinear gene regions within the genomes of the two species, and the colored curves emphasize the specific collinearity relationships among Q-type BvZFP genes between the two species.

Figure 6 .
Figure 6.Distribution of cis-acting elements in the promoter of Q−type BvZFP members in sugar beet.The phylogenetic tree was developed on MEGA 10 using neighbor-joining phylogenetic method analysis.Both the bootstrap test and the approximate likelihood ratio test were set to 1000 times.

Figure 6 .
Figure 6.Distribution of cis-acting elements in the promoter of Q-type BvZFP members in sugar beet.The phylogenetic tree was developed on MEGA 10 using neighbor-joining phylogenetic method analysis.Both the bootstrap test and the approximate likelihood ratio test were set to 1000 times.

Table 1 .
The physical and chemical properties of Q-type BvZFPs in sugar beet.

Table 2 .
The number of C2H2 and Q-type zinc finger proteins in different species.