Jasmonate-ZIM domain or JAZ genes control jasmonic acid signaling, which in turn controls a range of stress reactions and developmental activities in plants (Delgado et al. 2021). In this study, we have executed a genome-wide analysis of the JAZ gene family in cotton, and our classification corroborates with previous studies. With the advent of NGS technologies, the first genome to be sequenced belonged to A. thaliana, a plant of the Brassicaceae family (Initiative 2000). Other plant species, including G. arboreum, G. barbadense, G. hirsutum, G. herbaceum, G. raimondii, O. sativa, S. tuberosum, S. bicolor, B. rapa, and S. lycopersicum, then had their whole genome sequences completed (Murat et al. 2016). The availability of these genomes enhanced our comprehension of the evolutionary links among the various plant species. It provided a strong basis for the discovery of genome-wide genes and functional studies of JAZ genes (Koenig and Weigel 2015, Warwick et al. 2010). The JAZ subfamily has been the subject of earlier research that has concentrated on gene identification and expression response to different hormones or stress treatments in a variety of model plants, including Arabidopsis, and crops, including rice, maize, wheat, and soybean (Ebel et al. 2018, Huang et al. 2017, Ye et al. 2009, Zhang et al. 2015, Zhu et al. 2013). The investigations on the TIFY family have been described in B. rapa, B. oleracea, and B. napus exclusive to plants. In the oldest terrestrial plant moss, two TIFY genes have been found, whereas the TIFY gene family has been seen in all terrestrial plants (Saha et al. 2016, Swanson et al. 2004). Hence, it is postulated that the emergence of the TIFY gene family occurred during the evolutionary transition of plants from aqueous to land habitats. However, comprehensive characterization and comparative evolutionary information between all cotton species including G. herbaceum with a focus on the JAZ subfamily is still limited.
The JAZ gene family is known for its significant involvement in regulating several physiological processes related to plant growth and stress responses. These functions are mediated via the signaling pathways of jasmonic acid (JA), including seed germination, flower development, and the plant's reaction to environmental factors such as salt, drought, high temperature, wounds, and diseases (He et al. 2018, Ju et al. 2019, Liu, Zhang, Li and Xia 2019, Yu et al. 2018, Zhao et al. 2020). These proteins are actively involved in the cellular response to stress and serve as key mediators in the signaling pathways of phytohormones, facilitating the activation of genes responsible for defensive mechanisms. In addition to this, JAZ proteins play a crucial role in maintaining a delicate equilibrium between growth and defense mechanisms, hence preventing the adverse consequences associated with excessive immune responses triggered by the jasmonic acid (JA) defense pathway (Guo et al. 2018).
Even among closely related species, there were significant changes in genome size, shape, and copy number (Zhang et al. 2018). These variations were confined to repetitive sequences and impacted certain gene families engaged in various physiological processes in plants (Lysak et al. 2009). In this study, 13 JAZ genes identified A. thaliana, 14 in G. arboreum, 24 in G. barbadense, 24 in G. hirsutum, 14 in G. herbaceum, 15 in G. raimondii, 15 in O. sativa, 13 in S. tuberosum, 16 in S. bicolor, 26 in B. rapa, and 13 in S. lycopersicum. All JAZ genes were divided into 7 groups from A to G. Based on their cotyledons, all JAZ genes arranged themselves into each group. Clades A, B, C, D, and F contained protein members (Delgado, Mora-Poblete, Ahmar, Chen and Figueroa 2021) from monocots and dicots. Clade E has only monocots S. bicolor and O. sativa, while Clade G has only dicots. Genes in the same clade might suggest their similar or complementary physiological functions in plants.
The precise interaction between cis-acting elements (CARE) and transcription factors (TFs) inside the promoter region plays a crucial role in the regulation of gene expression. This mechanism serves as a fundamental process for biological signal transmission and facilitates gene cooperation with other genes. In Gossypium JAZ genes, a comprehensive analysis revealed the presence of 67 probable cis-acting regulatory elements. These elements are associated with 17 or more distinct activities. The regulatory elements identified in this study include the TC-rich repeat element, which is associated with defense and stress responsiveness (Lv et al. 2023), Additionally, the LTR element is involved in low-temperature responsiveness (Du et al. 2022), MYB-binding sites (MBS) are implicated in drought-inducibility (Heidari et al. 2021), while the MYB binding site (MRE) and ACE are involved in light responsiveness (Hartmann et al. 2005), the ABRE associated with ABA signaling pathway (Nakashima and Yamaguchi-Shinozaki 2013), AuxRR-core and TGA elements related to Auxin-responsiveness (Chen and Liu 2019), TCA-elements involved in salicylic acid-responsive elements (Chen et al. 2005), CGTCA and TGACG motifs involved in MeJA responsiveness (Wang et al. 2011), P-box, TATC-box, and GARE-motif involved in gibberellin-responsive. In addition, several cis-acting elements are linked to hormones, such as ABREs, P-boxes, TCA elements, and TGA elements. The presence of these cis-acting areas might possibly have a substantial impact on the involvement of JAZ genes in the disease resistance pathway regulated by jasmonic acid (JA).
In order to get a deeper understanding of the structural organization of JAZ genes in G. herbaceum, their conserved motifs were found. Two conserved motifs, namely JAZ and TIFY, were found in G. herbaceum. The first motif, identified by its Interpro ID as PF09425, is recognized as a Jas motif. This motif can engage in a binding interaction with jasmonate, adopting a partly unwound helical conformation. Furthermore, it can engage in interactions with Myc proteins, which serve as important mediators of jasmonate signaling in plants (Zhang et al. 2015). The second motif, identified as Interpro ID: IPR010399, corresponds to a Tify domain of a sequence of 36 amino acids. This domain is found solely in Embryophyta, a taxonomic group that encompasses terrestrial plants. The naming of this particular domain is attributed to its presence of the amino acid pattern (TIF[F/Y]XG), which is known for its remarkable conservation. However, it was previously denoted as the Zim domain (NISHII et al. 2000, Shikata et al. 2004, Vanholme et al. 2007, White 2006).
The examination of synteny among JAZ genes in the four cotton species revealed a robust collinearity pattern despite the presence of chromosomal rearrangements or gene duplication events. Gene duplication events are a significant mechanism by which plants undergo evolutionary processes and facilitate the expansion of gene families (Magadum et al. 2013). Research findings indicate that a significant proportion of angiosperms, ranging from 70% to 80%, have undergone gene duplication or polyploidization events (Qiao et al. 2019). From a biological evolutionary standpoint, fragment replication, tandem duplication, and translocation are mechanisms that are very successful in producing novel genes and facilitating the development of resistance against external threats (Panchy et al. 2016).
The present study highlights the ubiquitous significance of GhJAZ genes in stress physiology and highlights the crucial role these genes play in reacting to different environmental stressors. This pattern is duplicated in other plant species (Hu et al. 2013, Thines et al. 2007). JAZ01 exhibits observable differential regulation under salinity, heat, and drought. This is consistent with its function in ion balance under salinity stress (Qi et al. 2011), triggering heat shock responses to counteract thermal stress (Clarke et al. 2004) and jasmonate signaling pathways for drought stress management via ABA-dependent mechanisms (Seo et al. 2009). Because of its broad-spectrum stress reactivity, JAZ01 may be used as a target in crop breeding initiatives to increase stress tolerance.
Another important insight we obtained from this study was that there is a significant alternation of expression of other stress-related genes in treated plants. JAZ1, and associated genes including ACO, MPK3, WRKY40, MPK6, ERF1, and HSP40 cummulatively play a major role in plant stress tolerance [20], therefore, we accessed their response when it comes to reacting to different abiotic stresses. JAZ01, a component of the jasmonate signaling pathway, plays a vital role in regulating defensive responses in a number of different stress circumstances [20]. ACO, a crucial component in the production of ethylene, plays a pivotal role in regulating responses to drought and salt conditions by modulating ethylene levels. The members of the MAP kinase gene family having a role in signaling, MPK3 and MPK6, enhance the plant's resistance to environmental stress (Pitzschke and Hirt 2010, Sinha et al. 2011). WRKY40, a transcription factor, modulates the expression of genes that respond to abiotic signals by guiding the response to unfavorable conditions (Jiang and Deyholos 2009). ERF1 basically transmits the ethylene signals to plants, regulating the transcript level of stress-related genes (Müller and Munné-Bosch 2015). Additionally, HSP40, a unique heat shock protein, has a role in protecting plants from high temperatures by aiding the process of folding of proteins and their degradation (Sarkar et al. 2013).
The upregulation of JAZ01, ACO, MPK3, and WRKY40 genes and the specific regulation of MPK6, ERF1, and HSP40 genes in cotton plants under heat, salt, and drought conditions underlines their crucial function in the defense mechanisms against abiotic stress. The involvement of JAZ01 in the transmission of jasmonate signals and ACO in the generation of ethylene suggests a synchronized hormonal control that enhances the ability to withstand stress (Van de Poel et al. 2015, Wasternack and Hause 2013). The constant MPK3 enhanced expression across all stresses indicates its participation in widespread stress response. In contrast, MPK6 and ERF1 show stress-specific responses. Additionally, they interact with other stress-responsive pathways, such as abscisic acid (Müller and Munné-Bosch 2015, Rodriguez et al. 2010). The expression pattern of HSP40 suggests its role in enhancing protein stability under heat stress. Downregulation under exposure to salt and drought conditions suggests a change in chaperones to deal with (Wang et al. 2004). Increased expression of WRKY40 in all stress conditions gives it a wide-ranging regulatory role in defense pathways, such as osmotic balance, cellular safeguarding, and water conservation (Rushton et al. 2012). Collectively, JAZ01 and associated genes play important roles in a complex adaptative response, augmenting a plant's ability to endure unfavorable environmental conditions.