Skip to main content

In silico and computational analysis of zinc finger motif-associated homeodomain (ZF-HD) family genes in chilli (Capsicum annuum L)

Abstract

Zinc finger-homeodomain (ZHD) proteins are mostly expressed in plants and are involved in proper growth and development and minimizing biotic and abiotic stress. A recent study identified and characterized the ZHD gene family in chilli (Capsicum annuum L.) to determine their probable molecular function. ZHD genes with various physicochemical characteristics were discovered on twelve chromosomes in chilli. We separated ZHD proteins into two major groups using sequence alignment and phylogenetic analysis. These groups differ in gene structure, motif distribution, and a conserved ZHD and micro-zinc finger ZF domain. The majority of the CaZHDs genes are preserved, early duplication occurred recently, and significant pure selection took place throughout evolution, according to evolutionary study. According to expression profiling, the genes were found to be equally expressed in tissues above the ground, contribute to plant growth and development and provide tolerance to biotic and abiotic stress. This in silico analysis, taken as a whole, hypothesized that these genes perform distinct roles in molecular and phytohormone signaling processes, which may serve as a foundation for subsequent research into the roles of these genes in other crops.

Peer Review reports

Introduction

Capsicum annuum L., commonly known as chilli pepper, is a highly diverse species cultivated and traded globally for centuries. In our recent genome-wide analysis, we performed an in silico characterization of the Zinc Finger-HD (ZF-HD) transcription factor gene family in Capsicum annuum L [1]. Chili peppers are known for their pungency, which is due to the presence of capsaicinoids, a group of secondary metabolites that provide protection against herbivores and pathogens [2]. Our analysis revealed that the ZF-HD transcription factor gene family in Capsicum annuum L. includes multiple members with varying degrees of conservation across different cultivars [3]. These findings have important implications for improving our understanding of the molecular mechanisms underlying the diverse phenotypic traits observed in chilli pepper cultivars, including heat, color, and shape. In addition to their culinary uses, chilli peppers have also been traditionally used for their anti-inflammatory, analgesic, and antioxidant properties [4]. Plant's physiological activities and expression regulation are strongly linked to Transcription Factors [5], which control growth, development, and stress response [6]. The DNA-binding domains of each TF family are distinct, forming their binding specificity. A group of transcription factors is generated through zinc finger-homeodomain (ZHD) genes which are found to contribute differently to plants' growth and development [7]. It has been discovered that ZHD genes are connected with photosynthesis, defense mechanism, and various stress in plants [4,5,6]. ZHD proteins are expressed in floral tissue in advance or, more specifically, indicating that they can regulate flowering in plants [8]. Different regulatory proteins contain zinc finger motifs. These motifs are composed of two pairs of conserved cysteine or histidine residues which provide consistency to the motifs into a finger-shaped loop vis coordinating with a single zinc ion [9].

Furthermore, a single protein can carry one or more zinc finger motifs, mainly contributing to DNA-binding protein–protein interaction [10]. The zinc finger-homeodomain (ZHD) gene family members were first discovered in the C4 plant species Flaveria [11]. Previously ZF-HD1 was identified in Arabidopsis to observe its function, where it was found that ZF-HD1 binds with the EARLY RESPONSE TO DEHYDRATION STRESS 1 (ERD1) promoter due to induction of salt, increasing transpiration rate as well as abscisic acid [12]. Moreover, when this ZF-HD1 overexpressed and joined with NAC genes, it assists in developing drought tolerance in plants [12]. In certain crops, such as rice the members of the ZF-HD group have been identified. Specifically, researchers have discovered fourteen zinc finger homeobox genes in rice that belong to the ZF-HD group [13]. After conducting in-depth evolutionary analysis, which included examining the phylogenetic tree, duplication events, and comparative analysis, researchers could identify ancestral relationships between the ZF-HD gene and other genes. This analysis revealed how the ZF-HD gene had evolved and diversified over time, providing insights into its functional and structural properties For instance, MIF genes (Macrophage Migration Inhibitory Factor) exhibit only zinc finger while ZF-HD has both zinc finger as well as homeodomain, which indicates that MIF genes might have been produced from ZF-HD through many evolutionary changes or ZF-HD originated from MIF through gaining homeodomain portion [14].

ZHD-regulated transcription factors have been identified and studied in hot pepper (Capsicum annuum L.), contributing to our understanding of their function in plant growth, development, and stress response. Therefore, this in silico paper described the evolutionary trend, protein–protein interaction, and expression pattern of ZHD TFs in chilli under optimal and adverse conditions and explained the possible mechanisms of ZHD transcription factors in chilli plants.

Materials and methods

Sequence retrieval from the database

Hidden Markov Model (HMM) based conserved domain of ZHD gene family (PF) firstly downloaded from the Pfam database (Pfam 35.0: http://pfam.xfam.org/). The conserved domain HMM profile was used for BLASTP search against Sol Genomic (https://solgenomics.net/organism/Capsicum_annuum/genome) and China National Gene Bank (https://db.cngb.org/) databases for the release 2.0 cv. Zunla-1 genome with an expected cut-off value of 0.01 [15]. SMART (http://smart.embl-heidelberg.de/) and NCBI CDD (Conserved Domain Database) (https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi) were used to analyze the simple molecular architecture of the retrieved predicted amino acid sequences to find out either the sequences possessed the PF domain or not. Those which did not contain the required conserved domain were excluded [12, 13].

Physiochemical properties determination

ProtParam tool was used to retrieve some physiochemical properties like protein length (amino acid residues), molecular weight, isoelectric point [16], and GRVY, while the subcellular localization was found in the Cell-Ploc-20 [17]. Gene IDs, chromosomal positions, direction, and sequences of proteins and CDS of the potential genes were collected from Sol Genomics, but genomic sequences were retrieved from the NCBI (https://www.ncbi.nlm.nih.gov/). All those predicted Chilli’s ZHD family genes were renamed following their chromosomal locations.

Conserved motif analysis, domain prediction, and exon–intron distribution

The CaZHD gene family proteins’ motifs analysis was done with the help of MEME suite (https://meme-suite.org/meme/) [18] with the default conditions included motifs 20, minimum width 6, and maximum width 50. The hit data from the NCBI CDD, Norwich phytogenic data from MEGA 11.0 and MEME suite data were used in TBtools to visualize the conserved domains and motifs (TB) [19]. The genomic and CDS sequences of the CaZHD gene family are used in the gene structure display server (GSDS) web tool to analyze their exon and intron distribution [20].

Multiple sequence analysis and phylogenetic tree analysis

Crops like Chilli (Capsicum annuum) 11 and Arabidopsis (15) amino acid sequences of ZHD family were aligned with the help of ClustalX v2.1 multiple sequence alignment tool [21], and then for colorful visualization, the aligned data exported to Genedoc (https://www.nrbcs.org/gfx/genedoc/ebinet.htm) [22]. ZHD protein sequences of Arabidopsis, Maize, Tomato and chilli were used to build phylogenetic tree using the neighbor-joining (NJ) technique with 1000 bootstrap replicates with the help of MEGA 11.0 software [23] and representation of the tree updated using iTOL (https://itol.embl.de/) [24].

Cis‑regulatory elements analysis and function determination

The 1500 upstream promotor regions were extracted from NCBI (Link) for cis-regulatory elements extraction as they bind to the transcription factors and regulate target gene functions [25]. The PlantCARE (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/) web tool retrieves 5 to 20 bp putative cis-elements for the promotor region [26]. The measured cis-regulatory elements results are visualized in a heatmap with the help of TBtools [19].

Gene duplication and synteny analysis

A NCBI BLASTP search was done among the ZHD gene protein sequences of the chilli using 80% sequence identity for determining gene duplication [27]. The synonymous substitution rate (Ks), nonsynonymous substitution rate (Ka), and Ka/Ks ratio among the duplicated gene pairs were calculated with the help of KaKs calculator 2.0 [28]. A well-established formula T = Ks/2λ (where λ = 6.5 × 10− 9) was used to measure the evolutionary divergence. Gene duplication events among the chilli-Arabidopsis, chilli-tomato, and chilli-maize were analyzed using the Multiple Collinearity Scan toolkit (MCScanX) [29]. To exhibit the syntenic relationship of the ZHD genes of Chilli, Micro Synteny view software in TBtools was used to construct a map [19].

Gene ontogeny analysis

The ShinyGO v0.75: Gene Ontology Enrichment Analysis + more (http://bioinformatics.sdstate.edu/go/) was used to obtain gene ontology (GO) annotation against Capsicum annuum where potential chilli candidate proteins and/or genes IDs were subjected. The p-value cut-off (FDR) at 0.01 is set to calculate GO enrichment.

Transcriptome analysis

The expression profile of Zulna-1 for different parts such as fruit, flower, leaf, meristem, root and stem were extracted from the NCBI GEO database to analyze the organ-specific expression profile of CaZHD gene family (https://www.ncbi.nlm.nih.gov/geo/) [30]. For expression profiling, RNA-seq data unit values Read Per Kilobases per Million mapped reads (RPKM) were log2 folded, transformed and hierarchically displayed the heatmap through Heatmap Illustrator in TBtools [19].

qRT-PCR analysis for Investigating CaZHD genes in root

The samples were subjected to RNA isolation using RNA-easy Isolation Reagent from Vazyme, Nanjing, China. The quality of the isolated RNA was assessed using 0.8% agarose gel electrophoresis, while the purity and concentration were determined using a NanoDrop 2000 Spectrophotometer. Samples with OD260/OD280 ratios ranging from 1.90 to 2.10 were considered suitable for further experiments. For the analysis of CaZHD genes, the RNA samples were reverse-transcribed into cDNA using Hifair™ II 1st Strand cDNA Synthesis SuperMix for qPCR. The subsequent qRT-PCR analyses were carried out on a Light Cycler 480 instrument with a 20-μL reaction mixture (Primers Table S1). The Livak method was employed to calculate the relative levels of gene transcripts. Each RT-qPCR analysis was performed using three biological replicates. Protein-Protein interaction analysis.

With a high confidence score of 0.7, Protein-Protein Interaction (PPI) analysis was performed using STRING v11.0 (https://string-db.org/) [31]. The interactome's functional enrichment analysis was done somewhere at 0.01 level. The PPI network was created using active interaction based on various sources, including text mining, experiments, gene fusion, databases, co-expression, and an interaction score > 0.4. The physical and functional roles of the major candidate genes implicated were determined using this interactome map.

Putative microRNA target site analysis

The first mature miRNA from the PmiREM server (https://www.pmiren.com/) was used to determine the target site of chilli's 11 CaZHD gene family. The CDS of the 11 genes was then compared to the mature miRNA using the PsRNA online server tool (https://www.zhaolab.org/psRNATarget/) using the default setting [32]. Using the Cytoscape program (https://www.omicshare.com/tools/), a connection between the predicted miRNA was created [29, 30].

Results

Identification and physicochemical property of CaZHD TF gene family

A total of 11 ZF-HD genes were identified from the Sol Genomics, consisting of two conserved domains, the ZF domain and Micro Zinc Finger ZF domain (Fig. 1). All the 11 ZHD genes of chilli were named CaZHD1-CaZHD10, and micro zinc finger CaMIF1 was designated based on their corresponding location on the chromosomes 1-12 top to bottom. The length of the genes ranged from 125 -362aa, where CaMIF1 and CaZHD7 have the lowest and highest values, respectively (Table 1). CaMIF cellular location was either the nucleus or chloroplast, while the remaining CaZHDs were localized in the nucleus (Table 1). The CaZHDs were found to be more conserved as only five among eleven have only one intron (Fig. 1).

Fig. 1
figure 1

Phylogenetic relationships and gene structures of the ZHD genes from chilli. A The phylogenetic tree was constructed based on the full-length sequences of CaZHD genes. B Intron–Exon structures of the CaZHD genes. Purple boxes indicate exons, red boxes UTR region, and black lines indicate introns

Table 1 Information of eleven (11) ZHD genes of chilli (Capsicum annuum)

Conserved domain and gene structure analysis of ZHD

The conserved domain analysis showed that all the ZHD genes of Arabidopsis and chilli had two completely conserved domains: ZF Domain of ZHD subfamily (Fig. 2) and the ZF Domain of MIF subfamily (Fig. 2). Twenty conserved motifs in the ZHD gene family were identified using MEME online tools in Arabidopsis and chilli (Fig. 3, Supplementary Table 1). The results showed that most of genes of the ZHD family had four typical motifs, motif 1 (ZF-HD_dimer), motif 2 (Homeo_ZF_HD), motif 3 (Homeo_ZF_HD superfamily), and motif 4 (ZF-HD_prot_N), indicating that these four motifs were relatively conserved among the ZHD gene family (Fig. 3). Motif 1 and motif 3 were present in most of the ZHD genes of Arabidopsis and chilli except CaZHD10, AtZHD14, respectively. Moreover, motif 2 was not present in CaZHD10, AtMIF1 and CaMIF1, while Arabidopsis and chilli MIF do not contain motif4.

Fig. 2
figure 2

The sequence logos are based on alignments of CaZHD domains (A). ZF domain of ZHD subfamily and MIF subfamily is highly conserved across all ZHD proteins in Chilli. Multiple alignment analysis of CaZHD domains was performed with MUSCLE (B)

Fig. 3
figure 3

The distribution of 20 motifs on ZHD proteins of Arabidopsis and chilli by using MEME version 4.9.0 and interlinking it with a phylogenetic tree to better understand their association. The bars represent motifs with different color codes for different types of motifs

Phylogenetic analysis and classification of the ZHD transcription factor family

The phylogenetic tree was constructed with ZHD members of chilli (11), Arabidopsis (15), Maize (21), and Solanum lycopersicon (22). The result showed that these 69 ZHD genes could be divided into two groups (Group I-II), where Group I was divided into two subgroups (Ia and Ib) and Group II was divided into three subgroups (IIa, IIb, and IIc) (Fig. 4). This classification was done based on the MIF group's presence or absence. Both MIF genes from chilli and Arabidopsis are present in Clade I, consisting of CaZHD1, CaZHD2, and CaZHD10, along with AtZHD5, AtZHD7, AtZHD13 and AtZHD14 (Supplementary Table 2). Subgroup IIc had no CaZHD gene, while IIa contained most of the CaZHD genes.

Fig. 4
figure 4

Phylogenetic Relationship among ZHD genes of Capsicum annuum, A. thaliana, S. lycopersicum, and Zea mays was studied. Capsicum annuum genes are marked with a red and blue triangle. The evolutionary history was inferred using the NJ method with 1000 Bootstrap. This analysis involved 69 ZHD genes. Evolutionary analyses were conducted in MEGA 11 [23]

Gene chromosome mapping and collinearity analysis

The physical location of the chilli CaZHDs was identified in the chromosomes, but only CaMIF was located in the scaffold region of the chilli genome region (Fig. 5a). We arranged them on a pseudo-chromosome, designated as 1-12, concatenated by the unplaced scaffolds. Chromosomes 2, 3, 4, 5, and 8 contained these 11 ZHD genes, while chromosome 2 contained the highest 3 genes (Fig. 5a). The inter-relation of the genes is shown with a blue line in Fig. 5a, where Fig. 5b clarifies the location of genes in the chromosomes and scaffold. To understand the evolutionary mechanism of the ZHD gene family of chilli, we analyzed the collinear relationship between chilli and Arabidopsis, chilli and tomato, and chilli and maize (Fig. 5b). The results show that 7 pairs of collinearity genes of ZHD were between chilli and tomato, followed by chilli and Arabidopsis (5 pairs), and the least was chilli and maize (1 pair).

Fig. 5
figure 5

a Distribution of CaZHDs on chilli chromosomes, lines predicting the possible gene duplication on different chromosomes; b Genome-wide synteny analysis of CaZHDs and chromosomal distribution and intrachromosomal relationships of ZHD genes between chilli- tomato, chilli- Arabidopsis and chilli- (monocot) Zea mays. Gray lines indicate all synteny blocks in the Capsicum annuum genome, and the red lines indicate duplicated ZHD gene pairs. The chromosome number is indicated on the top of each chromosome, showing the dominance of segmental duplication and the rare occurrence of tandem duplication

According to the result of the evolution analysis, the values of Ka, Ks, and Ka/Ks were obtained (Table 2). Nineteen gene pairs were identified using Tbtools. The value of Ka/Ks of each pair ranged from 0.23 to 0.29 (Ka/Ks < 1). This result indicated that all of them had undergone strong purifying selection. The CaMIF1 and CaZHD4 occurred between 213.50 MYA (Million years ago), while CaZHD6-CaZHD7 gene pair’s duplication happened most recently (around 42 MYA).

Table 2 Ka/Ks ratio duplicated gene pairs in chilli

Cis-acting element prediction of ZHD gene family

The results of the ZHD gene family of ciss-acting element prediction of the upstream 2000 bp sequence showed various environmental and stress response elements in the ZHD gene family of chilli (Fig. 6). All CaZHDs contained light response elements STREs, while 95% of the CaZHDs contained ERE, WUN-Motif, Box-4, MYB, MYCas-1 and W-box response elements. The most common cis-elements found in our analysis were GA-motif, P-box, TCA, GARE-Motif, TGA-element, and WRE3, which were present in relatively low numbers. The remaining cis-elements accounted for approximately 50-70% of the total identified. Most of the present cis-regulatory elements were phytohormone and abiotic stress-responsible elements. For instance, ABRE cis-element regulates abscisic acid, TGA-element regulates auxin phytohormone and TCA-element salicylic acid where light-responsive STRE, AT1-motif, and Box-4 and abiotic stress-responsive MYB, MYC, ERE, and ARE present in the different CaZHDs (Supplementary File S1).

Fig. 6
figure 6

Cis-regulatory elements in putative CaZHD promoter’s process are associated with different plant developmental processes

Expression pattern analysis of GRF gene family

Based on the transcriptome data of Capsicum annum, the tissue expression pattern of the ZHD gene family was analyzed, and fourteen samples, including Root, Stem, Leaf, Bud, Flower, F-Dev-1, F-Dev-2, F-Dev-3, F-Dev-4, F-Dev-5, F-Dev-6, F-Dev-7, F-Dev-8, and F-Dev-9 were selected for prediction. The result showed that the fourteen samples expressed the CaZHDs genes differently. The ZHD genes were expressed in Root, Stem, Leaf, Bud, Flower, F-Dev-1, F-Dev-2, F-Dev-3, F-Dev-4, F-Dev-5, F-Dev-6, F-Dev-7, F-Dev-8, and F-Dev-9, especially in the stem, leaf, flower bud, medium flower bud, and small fruit (Fig. 7a). Most the genes are expressed in leaf and fruit development stage 1 (F-Dev-1), while no gene expression in the root is sometimes expressed adversely. CaZHD2, 3, 5, 8 and 9 genes are responsible for the fruit maturity stage, while CaZHD6 is highly express in bud and flower. This CaZHD6 gene act against pathogens like nematodes resembling the AtZHD9 gene and CaZHD8, and CaZHD9 genes help ameliorate abiotic stress like drought and salinity (Table 3).

Fig. 7
figure 7

a The heat map shows the expression profile of the CaZHDs in different organs. b The heat map shows the expression profile of the CaZHDs in root organ under salt stress

Table 3 Gene ontology enrichment analysis of CaZHDs their GO functions, Gene’s expression, Orthologs in Arabidopsis, and their functions are presented in the table

Investigating of CaZHD genes in root organ

To investigate the biological roles of CaZHD genes family by examining their expression levels in the roots of chili plants (Capsicum spp.). Our aim was to understand how these genes respond to stress conditions, specifically a high sodium chloride (NaCl) concentration of 3000 ppm from 6 weeks after transplanting [46].

Using the quantitative real-time polymerase chain reaction (qPCR) technique, we analyzed the expression patterns of 11 CaZHD genes family in the chili roots (Fig. 7b). Our findings revealed a consistent down-regulation in the expression levels of most of the ZF-HD genes examined, such as CaMIF1, CaZHD1, CaZHD2, CaZHD4, CaZHD5, CaZHD6, CaZHD7, CaZHD8, CaZHD9 and CaZHD10 under the NaCl stress conditions. However, amidst this general down-regulation trend, we made an intriguing observation. CaZHD3 displayed an up-regulation in its expression specifically in response to the NaCl stress. This differential response of CaZHD3 sets it apart from the other ZF-HD genes suggesting that it may have a distinct role in the plant's response to salt stress. Our study focused on the roots of chili plants to better understand how this crucial plant organ responds at the molecular level to the stress caused by high sodium chloride (NaCl) concentrations. By employing qPCR, we were able to assess the expression levels of the ZHD genes and highlight the potential significance of CaZHD3 as a key player in the plant's adaptive response to high NaCl concentrations. Further investigations are warranted to unravel the specific functions and regulatory mechanisms of CaZHD3 within the root system of chili plants, with the aim of enhancing salt tolerance in this important crop species.

Annotation and ontology of the CaZHD gene family

To have a general knowledge of the genes of the CaZHD gene family, the transcripts of 11 CaZHDs were annotated and categorized with gene ontology (GO). Only two of the 11 transcripts were annotated and categorized into all three primary categories, Biological Process (BP), Molecular Function (MF), and Cellular Component (CC) (Fig. 8). These two annotated genes are highly responsible for the molecular function, including transcription cis-regulatory region binding, transcription regulatory region nucleic acid binding, sequence-specific double-stranded DNA binding, double-stranded DNA binding, sequence-specific DNA binding and transcription regulatory function. These GO results of the CaZHDs were in accordance with the transcription factors' functions. The network of the molecular function depicted how they interrelated with other molecular functions and how they work (Fig. 9).

Fig. 8
figure 8

Fold Enrichment chart representing the overlapping CaZHDs functions. Red color dot plots represent the more no. of genes involved in that process and vice versa for small blue sizes

Fig. 9
figure 9

Static Network Enrichment graph showing the network of CaZHDs functions. Darker nodes are more significantly enriched gene sets. Bigger nodes represent larger gene sets. Thicker edges represent more overlapped genes

MiRNA target site prediction and validation

MicroRNAs (miRNAs), a subclass of non-coding short RNAs with an average length of 21 nucleotides (nt), are essential for many biological processes in plants, including development, signal transduction, and responses to biotic and abiotic stress [47]. The chilli genome contains 29 conserved and 35 new miRNA families [48]. Only 9 miRNAs were found despite the prediction of 398 mature miRNA sequences against the CDSs of CaZHDs (Supplementary File S2, Table 5). Can-miR482a has ten target sites in CaZHD3, Can-miRN450a, Can-miRN450b, Can-miRN473a, Can-miRN473b, and Can-miRN473c has seven target sites in CaZHD2, and miRN4099 and miRN4186 have four target sites in CaZHD4 and CaZHD8, respectively (Table 4).

Table 4 CaZHDs targeting Putative miRNA functions and their targeted genes

Protein–protein interaction and predicted protein structure

We used the STRING database to predict potential interactions among the proteins (https://string-db.org/). Only 5 proteins among 11 CaZHD were correlated to each at the highest level (0.900) of confidence (Fig. 10). CaZHD2 is highly associated with other correlated genes and is speculated to play a central role in the expression and signaling. Furthermore, we illustrated protein structure with the help of online software to identify the protein structure.

Fig. 10
figure 10

Interaction network of the ZHD proteins in Capsicum annuum. Deep ash-colored lines indicate confidence level (We selected the highest level of confidence (0.900) as the measuring unit). Protein–protein interactions (PPIs) play a crucial role in cellular functions and biological processes, including cell–cell interactions and metabolic and developmental control in all organisms

Discussion

Zinc finger Homeodomain (ZHD) transcription factors are only found in plants and play vital roles in plant growth and development and protect against biotic and abiotic stress [49]. This gene family has so far been investigated in Arabidopsis, maize, tomato, rice, cotton and wheat but not in economically beneficial hot pepper or chilli (Capsicum annuum L.) [4, 5, 9, 50]. Sol genomic and BLASTp studies revealed that chilli has about 11 ZHD genes. These CaZHD genes have persevered as most of the genes have no introns, and the most recent duplication occurred around 42 MYA (Table 2). ZHD genes were exclusively found in terrestrial plants and were more prevalent during the evolution of angiosperms [50, 51]. For the structural analysis of CaZHD domains, numerous evolutionary and structural studies have been carried out. Various techniques, including investigations of synteny, gene duplication, gene structure, phylogenetic trees, and motif organization were used. The size and types of conserved sections essentially present in an acceptable sequence alignment dictate the function of a gene family. Sequence comparisons revealed that motifs 1 and 4, together known as the ZF dimer (Fig. 2), are consistent with those of other plant species [50, 51], suggesting that ZHD proteins have a structure similar to those of other plants. In our tree analysis, CaZHDs were divided into two groups with six subgroups (Fig. 4), which aligns with other phylogenetic studies on the crops [4, 5, 52]. In subgroup IIb, most of the proteins from Arabidopsis, Tomato, and chilli were present, whereas in subgroup IIc, the CaZHD protein was absent, which made up most of the monocot maize protein in the tree (Fig. 4). This difference indicated that proteins from both monocots and dicots have diverged. The evolutionary insights of gene families can be extracted from the architectures of the genes belonging to the gene family [53]. Our gene structure research has revealed that many of the ZF-HD/ZHD genes lack introns (Fig. 1), a feature that has been observed in other species as well [4, 52, 54]. This suggests that the absence of introns in these genes may contribute to their ability to withstand environmental stresses, as there are no mutations that could affect their function in response to stress Therefore, according to the phylogenetic study, we can speculate that CaZHD family is an old gene family that emerged after Angiosperm split from Gymnosperm and before Dicots split from Monocots. The gene family has significantly grown since then, with the most recent gene duplication occurring between 25 and 50 MYA. Even though the feature domains of the gene family have mostly remained preserved, the gene family members have substantially diversified in terms of their nucleotide sequences, locations in the genome, and associated functionalities. The CaZHDs gene family was divided into two groups along with total five subgroups or subfamilies based on phylogenetic analysis and the presence or lack of conserved domains; thus, each subfamily has a unique set of conserved domains and motifs [55].

Furthermore, the functional differentiation of the CaZHDs genes has also been revealed at the gene expression level in different plant tissues, across developmental stages, and cultivars. For instance, almost no expression of 11 CaZHDs gene transcripts was found in all plant roots (Fig. 6). CaMIF1, CaZHD2, 3, 4, 8 and 9 are highly expressed in the leaf and early fruit development stage, while CaZHD5 and CaZHD6 were highly expressed in buds and flowers (Fig. 6). Among all eleven genes, CaZHD5 and CaZHD6 help to ameliorate biotic stress, and CaZHD8 and CaZHD9 ameliorate abiotic stress (Table 3), which resembles to the recent study [56]. From all the results taken together, we speculated that these genes, which were highly expressed in the above-ground plant parts, helped plants overcome adverse environmental conditions. These genes were also co-expressed (Fig. 10), where CaZHD2 was centrally interrelated with others, and its main function was controlling the developmental-like phloem lineage [37, 39]. Therefore, it is clear that CaZHD2 plays a central role in developmental and biotic and abiotic stress tolerance by controlling the other gene expression. This result was acceptable as CaZHD2 [57,58,59] gene consisted of sequence-specific microRNAs, which helped to overcome biotic and abiotic stress (Table 5). During qtPCR analysis, our study found that most ZF-HD genes in the CaZHD gene family were down-regulated in the roots of chili plants under high NaCl stress, indicating their involvement in the plant's response to salt stress. However, CaZHD3 showed an intriguing up-regulation specifically in response to NaCl stress, suggesting it may have a distinct role in enhancing salt tolerance. Further investigations are needed to uncover the precise functions and regulatory mechanisms of CaZHD3 in the root system, contributing to a comprehensive understanding of salt stress adaptation in chili plants [55, 60,61,62]. Though some of the miRNA functions are still not found, we need to explore to know the more complex mechanisms. microRNAs are important plant regulators that regulate almost every biological process, from growth and development to combating pathogens and maintaining proper internal conditions [63, 64]. The molecular mechanism of the 11 CaZHDs was predicted through GO and network analysis. Only 2 genes showed molecular functions, and their molecular functions were highly correlated with each and made an intriguing network (Figs. 8 and 9). In addition, as the CaZHD genes' cis-regulatory elements control the phytohormonal signaling and abiotic stress tolerance mechanism, we can predict that these genes' microRNAs are also correlated with them as they are sequence-specific [62, 67, 68].

Table 5 CaZHDs targeting Putative miRNA functions along with their targeted genes

Conclusions

It has been determined that the CaZHD gene family in hot pepper (Capsicum annuum L.) consists of 11 genes. Before the Monocots and Dicots division and following the separation between Angiosperm and Gymnosperm, the CaZHD gene family most likely developed. The family was divided into two subfamilies, and these subfamilies differed substantially in terms of chromosomal position, nucleotide sequence, and GO annotation and categorization. Even while the expression of the genes in the family varied significantly in four-year-old plant tissues, four-year-old plant roots, and four-year-old plant roots from various cultivars, the functional relationships between the genes in the family persisted. Two CaZHDs were crucial in the chilli plant's biotic stress, and two were confirmed to control the plant's response to abiotic stress.

Availability of data and materials

All data generated or analysed during this study are included in this published article and its supplementary information files.

References

  1. Jang H-H, Lee J, Lee S-H, Lee Y-M. Effects of Capsicum annuum supplementation on the components of metabolic syndrome: a systematic review and meta-analysis. Sci Rep. 2020;10(1):20912.

    CAS  PubMed  PubMed Central  Google Scholar 

  2. de Silva DD, et al. Identification, prevalence and pathogenicity of Colletotrichum species causing anthracnose of Capsicum annuum in Asia. IMA Fungus. 2019;10(1):1–32.

    Google Scholar 

  3. Rasekh M, Karami H, Fuentes S, Kaveh M, Rusinek R, Gancarz M. Preliminary study non-destructive sorting techniques for pepper (Capsicum annuum L.) using odor parameter. LWT. 2022;164:113667.

    CAS  Google Scholar 

  4. Baruah PM, Krishnatreya DB, Bordoloi KS, Gill SS, Agarwala N. Genome wide identification and characterization of abiotic stress responsive lncRNAs in Capsicum annuum. Plant Physiol Biochem. 2021;162:221–36.

    CAS  PubMed  Google Scholar 

  5. Zhang H, et al. Comparative transcriptome analysis reveals molecular defensive mechanism of Arachis hypogaea in response to salt stress. Int J Genomics. 2020;2020:652–93.

  6. Gibson SI. Control of plant development and gene expression by sugar signaling. Curr Opin Plant Biol. 2005;8(1):93–102. https://doi.org/10.1016/j.pbi.2004.11.003.

    Article  CAS  PubMed  Google Scholar 

  7. Wang H, et al. Genome-Wide Identification, Evolution and Expression Analysis of the Grape (Vitis vinifera L.) Zinc Finger-Homeodomain Gene Family. Int J Mol Sci. 2014;15(4):5730–48. https://doi.org/10.3390/ijms15045730.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Tan QK-G, Irish VF. The Arabidopsis Zinc Finger-Homeodomain Genes Encode Proteins with Unique Biochemical Properties That Are Coordinately Expressed during Floral Development. Plant Physiol. 2006;140(3):1095–108. https://doi.org/10.1104/pp.105.070565.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Krishna SS. Structural classification of zinc fingers: survey and summary. Nucleic Acids Res. 2003;31(2):532–50. https://doi.org/10.1093/nar/gkg161.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Mackay JP, Crossley M. Zinc fingers are sticking together. Trends Biochem Sci. 1998;23(1):1–4. https://doi.org/10.1016/S0968-0004(97)01168-7.

    Article  CAS  PubMed  Google Scholar 

  11. Shalmani A, et al. Zinc Finger-Homeodomain Genes: Evolution, Functional Differentiation, and Expression Profiling Under Flowering-Related Treatments and Abiotic Stresses in Plants. Evol Bioinforma. 2019;15:117693431986793. https://doi.org/10.1177/1176934319867930.

    Article  Google Scholar 

  12. Tran L-SP, et al. Co-expression of the stress-inducible zinc finger homeodomain ZFHD1 and NAC transcription factors enhances expression of the ERD1 gene in Arabidopsis: ZFHD and NAC in drought stress response. Plant J. 2006;49(1):46–63. https://doi.org/10.1111/j.1365-313X.2006.02932.x.

    Article  CAS  Google Scholar 

  13. Jain M, Tyagi AK, Khurana JP. Genome-wide identification, classification, evolutionary expansion and expression analyses of homeobox genes in rice: Genomic analyses of homeobox genes in rice. FEBS J. 2008;275(11):2845–61. https://doi.org/10.1111/j.1742-4658.2008.06424.x.

    Article  CAS  PubMed  Google Scholar 

  14. Hu W, de Pamphilis CW, Ma H. Phylogenetic Analysis of the Plant-specific Zinc Finger-Homeobox and Mini Zinc Finger Gene Families. J Integr Plant Biol. 2008;50(8):1031–45. https://doi.org/10.1111/j.1744-7909.2008.00681.x.

    Article  CAS  PubMed  Google Scholar 

  15. Qin C, et al. Whole-genome sequencing of cultivated and wild peppers provides insights into Capsicum domestication and specialization. Proc Natl Acad Sci. 2014;111(14):5135–40.

    CAS  PubMed  PubMed Central  Google Scholar 

  16. Ahmed L, Rebaz O. Spectroscopic properties of Vitamin C: A theoretical work. Cumhur Sci J. 2020;41(4):916–28.

    Google Scholar 

  17. Chou K-C, Shen H-B. Cell-PLoc 2.0: an improved package of web-servers for predicting subcellular localization of proteins in various organisms. Nat Sci. 2010;2(10):1090.

    CAS  Google Scholar 

  18. Bailey TL, Johnson J, Grant CE, Noble WS. The MEME suite. Nucleic Acids Res. 2015;43(W1):W39–49.

    CAS  PubMed  PubMed Central  Google Scholar 

  19. Chen C, et al. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020;13(8):1194–202.

    CAS  PubMed  Google Scholar 

  20. Hu B, Jin J, Guo A-Y, Zhang H, Luo J, Gao G. GSDS 2.0: an upgraded gene feature visualization server. Bioinformatics. 2015;31(8):1296–7.

    PubMed  Google Scholar 

  21. Larkin MA, et al. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23(21):2947–8.

    CAS  PubMed  Google Scholar 

  22. Nicholas KB. GeneDoc: analysis and visualization of genetic variation. Embnew News. 1997;4:14.

    Google Scholar 

  23. Tamura K, Stecher G, Kumar S. MEGA11: molecular evolutionary genetics analysis version 11. Mol Biol Evol. 2021;38(7):3022–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  24. Letunic I, Bork P. Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021;49(W1):W293–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  25. Fang Y, You J, Xie K, Xie W, Xiong L. Systematic sequence analysis and identification of tissue-specific or stress-responsive genes of NAC transcription factor family in rice. Mol Genet Genomics. 2008;280(6):547–63.

    CAS  PubMed  Google Scholar 

  26. Lescot M, et al. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002;30(1):325–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  27. Song X, Li Y, Hou X. Genome-wide analysis of the AP2/ERF transcription factor superfamily in Chinese cabbage (Brassica rapa ssp. pekinensis). BMC Genomics. 2013;14(1):1–15.

    CAS  Google Scholar 

  28. Wang D-P, Wan H-L, Zhang S, Yu J. γ-MYN: a new algorithm for estimating Ka and Ks with consideration of variable substitution rates. Biol Direct. 2009;4(1):1–18.

    CAS  Google Scholar 

  29. Wang Y, et al. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012;40(7):e49–e49.

    CAS  PubMed  PubMed Central  Google Scholar 

  30. Barrett T, et al. NCBI GEO: archive for functional genomics data sets—update. Nucleic Acids Res. 2012;41(D1):D991–5.

    PubMed  PubMed Central  Google Scholar 

  31. Szklarczyk D, et al. The STRING database in 2017: quality-controlled protein–protein association networks, made broadly accessible. Nucleic Acids Res. 2017;45(D1):gkw937.

  32. Sridhar J, Sowmiya G, Sekar K, Rafi ZA. PsRNA: A Computing Engine for the comparative identification of putative small RNA locations within intergenic regions. Genom Proteomics Bioinform. 2010;8(2):127–34.

    CAS  Google Scholar 

  33. Pushpalatha C. et al. Antibacterial nanozymes: an emerging innovative approach to oral health management. Top Catal. 2010;8(2):1–12.

  34. Otasek D, Morris JH, Bouças J, Pico AR, Demchak B. Cytoscape automation: empowering workflow-based network analysis. Genome Biol. 2019;20(1):1–15.

    Google Scholar 

  35. Zapata L, Ding J, Willing EM, Hartwig B, Bezdan D, Jiao WB, Patel V, Velikkakam James G, Koornneef M, Ossowski S, Schneeberger K. Chromosome-level assembly of Arabidopsis thaliana L er reveals the extent of translocation and inversion polymorphisms. Proc Natl Acad Sci. 2016;113(28):E4052–60.

  36. Su G, Morris JH, Demchak B, Bader GD. Biological network exploration with Cytoscape 3. Curr Protoc Bioinforma. 2014;47(1):8–13.

    Google Scholar 

  37. Ascencio-Ibánez JT, et al. Global analysis of Arabidopsis gene expression uncovers a complex array of changes impacting pathogen response and cell cycle during geminivirus infection. Plant Physiol. 2008;148(1):436–54.

    PubMed  PubMed Central  Google Scholar 

  38. Bueso E, Serrano R, Pallás V, Sánchez-Navarro JA. Seed tolerance to deterioration in arabidopsis is affected by virus infection. Plant Physiol Biochem. 2017;116:1–8.

    CAS  PubMed  Google Scholar 

  39. Gaudinier A, et al. Transcriptional regulation of nitrogen-associated metabolism and growth. Nature. 2018;563(7730):259–64.

    CAS  PubMed  Google Scholar 

  40. Hu W, Feng B, Ma H. Ectopic expression of the Arabidopsis MINI ZINC FINGER1 and MIF3 genes induces shoot meristems on leaf margins. Plant Mol Biol. 2011;76(1):57–68.

    CAS  PubMed  Google Scholar 

  41. Altmann M, et al. Extensive signal integration by the phytohormone protein network. Nature. 2020;583(7815):271–6.

    CAS  PubMed  Google Scholar 

  42. Smit ME, et al. A PXY-mediated transcriptional network integrates signaling mechanisms to control vascular development in Arabidopsis. Plant Cell. 2020;32(2):319–35.

    CAS  PubMed  Google Scholar 

  43. Ballvora A, et al. Comparative sequence analysis of Solanum and Arabidopsis in a hot spot for pathogen resistance on potato chromosome V reveals a patchwork of conserved and rapidly evolving genome segments. BMC Genomics. 2007;8(1):1–13.

    Google Scholar 

  44. Perotti MF, Ribone PA, Cabello JV, Ariel FD, Chan RL. AtHB23 participates in the gene regulatory network controlling root branching, and reveals differences between secondary and tertiary roots. Plant J. 2019;100(6):1224–36.

    CAS  PubMed  Google Scholar 

  45. Barth O, Vogt S, Uhlemann R, Zschiesche W, Humbeck K. Stress induced and nuclear localized HIPP26 from Arabidopsis thaliana interacts via its heavy metal associated domain with the drought stress related zinc finger transcription factor ATHB29. Plant Mol Biol. 2009;69(1):213–26.

    CAS  PubMed  Google Scholar 

  46. Bollier N, et al. At-MINI ZINC FINGER2 and Sl-INHIBITOR OF MERISTEM ACTIVITY, a conserved missing link in the regulation of floral meristem termination in Arabidopsis and tomato. Plant Cell. 2018;30(1):83–100.

    CAS  PubMed  PubMed Central  Google Scholar 

  47. Hu J, Zhang H, Ding Y. Identification of conserved microRNAs and their targets in the model legume Lotus japonicus. J Biotechnol. 2013;164(4):520–4.

    CAS  PubMed  Google Scholar 

  48. Hwang D-G, et al. The hot pepper (Capsicum annuum) microRNA transcriptome reveals novel and conserved targets: a foundation for understanding microRNA functional roles in hot pepper. PLoS ONE. 2013;8(5):e64238.

    CAS  PubMed  PubMed Central  Google Scholar 

  49. Li Y, et al. Genome-wide identification and expression analysis of HD-ZIP I gene subfamily in Nicotiana tabacum. Genes. 2019;10(8):575.

    CAS  PubMed  PubMed Central  Google Scholar 

  50. Zhang J, et al. Gelsolin inhibits malignant phenotype of glioblastoma and is regulated by miR-654-5p and miR-450b-5p. Cancer Sci. 2020;111(7):2413–22.

    CAS  PubMed  PubMed Central  Google Scholar 

  51. Eyles RP. microRNA involvement in root organ formation and function in Medicago truncatula. 2014.

    Google Scholar 

  52. Shivaprasad PV, Chen H-M, Patel K, Bond DM, Santos BA, Baulcombe DC. A microRNA superfamily regulates nucleotide binding site–leucine-rich repeats and other mRNAs. Plant Cell. 2012;24(3):859–74.

    CAS  PubMed  PubMed Central  Google Scholar 

  53. Lai W, Zhu C, Hu Z, Liu S, Wu H, Zhou Y. Identification and transcriptional analysis of zinc finger-homeodomain (ZF-HD) family genes in cucumber. Biochem Genet. 2021;59(4):884–901.

    CAS  PubMed  Google Scholar 

  54. Zhou C, et al. Genome-wide analysis of zinc finger motif-associated homeodomain (ZF-HD) family genes and their expression profiles under abiotic stresses and phytohormones stimuli in tea plants (Camellia sinensis). Sci Hortic. 2021;281:109976.

    CAS  Google Scholar 

  55. Lu S, et al. CDD/SPARCLE: the conserved domain database in 2020. Nucleic Acids Res. 2020;48(D1):D265–8.

    CAS  PubMed  Google Scholar 

  56. Shalmani A, et al. Zinc finger-homeodomain genes: evolution, functional differentiation, and expression profiling under flowering-related treatments and abiotic stresses in plants. Evol Bioinforma. 2019;15:1176934319867930.

    Google Scholar 

  57. Islam MAU, et al. Genome-Wide Identification and In Silico Analysis of ZF-HD Transcription Factor Genes in Zea mays L. Genes. 2022;13(11):2112.

    CAS  PubMed  PubMed Central  Google Scholar 

  58. Khatun K, et al. Genome-wide analysis and expression profiling of zinc finger homeodomain (ZHD) family genes reveal likely roles in organ development and stress responses in tomato. BMC Genomics. 2017;18(1):1–16.

    Google Scholar 

  59. Nai G. et al. Genome-wide analysis and expression characterization of zinc finger homeodomain (ZHD) family genes responsed to different abiotic stresses and hormonal treatments in grape (Vitis vinifera L.),” In Review, preprint. 2020. https://doi.org/10.21203/rs.3.rs-29270/v1.

  60. Schultz J, Copley RR, Doerks T, Ponting CP, Bork P. SMART: a web-based tool for the study of genetically mobile domains. Nucleic Acids Res. 2000;28(1):231–4.

    CAS  PubMed  PubMed Central  Google Scholar 

  61. Taylor-Teeples M, et al. An Arabidopsis gene regulatory network for secondary cell wall synthesis. Nature. 2015;517(7536):571–5.

    CAS  PubMed  Google Scholar 

  62. Rustikawati R, Herison C, Sutrawati M, Umroh D. Assessment of salinity tolerance on chili pepper genotypes. E3S Web of Conferences. 2023.

    Google Scholar 

  63. Hu J, Gao Y, Zhao T, Li J, Yao M, Xu X. Genome-wide identification and expression pattern analysis of zinc-finger homeodomain transcription factors in tomato under abiotic stress. J Am Soc Hortic Sci. 2018;143(1):14–22.

    CAS  Google Scholar 

  64. Abdullah M, et al. Zinc finger-homeodomain transcriptional factors (ZHDs) in upland cotton (Gossypium hirsutum): Genome-wide identification and expression analysis in fiber development. Front Genet. 2018;9:357.

    CAS  PubMed  PubMed Central  Google Scholar 

  65. Zhu Q-H, Fan L, Liu Y, Xu H, Llewellyn D, Wilson I. miR482 regulation of NBS-LRR defense genes during fungal pathogen infection in cotton. PLoS ONE. 2013;8(12):e84390.

    PubMed  PubMed Central  Google Scholar 

  66. Muys BR, et al. miR-450a Acts as a Tumor Suppressor in Ovarian Cancer by Regulating Energy MetabolismmiR-450a Acts as a Tumor Suppressor in Ovarian Cancer. Cancer Res. 2019;79(13):3294–305.

    CAS  PubMed  PubMed Central  Google Scholar 

  67. Bueso E, et al. Arabidopsis thaliana homeobox25 uncovers a role for Gibberellins in seed longevity. Plant Physiol. 2014;164(2):999–1010.

    CAS  PubMed  Google Scholar 

  68. Wang W, Wu P, Li Y, Hou X. Genome-wide analysis and expression patterns of ZF-HD transcription factors under different developmental tissues and abiotic stresses in Chinese cabbage. Mol Genet Genomics. 2016;291(3):1451–64.

    CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

Not applicable.

Author information

Authors and Affiliations

Authors

Contributions

MAUI carried out research work and wrote initial draft of manuscript. JBN, and MS carried out the use of software and data analysis. QA, AS, and MAS edited final version of manuscript. All authors reviewed final version of manuscript and approved final version.

Corresponding authors

Correspondence to Muhammad Shafiq or Qurban Ali.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Additional file 2.

Additional file 3: Supplementary Table 1.

Information regarding Motifs.

Additional file 4: Supplementary Table 2.

Chili’s ZHD gene family distribution among groups based on phylogenetic analysis.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Islam, M.A.U., Nupur, J.A., Shafiq, M. et al. In silico and computational analysis of zinc finger motif-associated homeodomain (ZF-HD) family genes in chilli (Capsicum annuum L). BMC Genomics 24, 603 (2023). https://doi.org/10.1186/s12864-023-09682-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-023-09682-x

Keywords