Asima Rasool1
Farrukh Azeem1*
Mahmood Ur-Rahman1
Muhammad Rizwan2Muhammad Hussnain Siddique1Daniyah Habiballah Bay3
Najat Binothman4Najla Amin T. Al Kashgry5
Sameer H. Qari6*- 1Department of Bioinformatics and Biotechnology, Government College University Faisalabad, Faisalabad, Pakistan
- 2Department of Environmental Sciences, Government College University Faisalabad, Faisalabad, Pakistan
- 3Department of Biology, Faculty of Applied Science, Umm Al-Qura University, Makkah, Saudi Arabia
- 4Department of Chemistry, College of Sciences & Arts, King Abdulaziz University, Rabigh, Saudi Arabia
- 5Department of Biology, College of Science, Taif University, Taif, Saudi Arabia
- 6Department of Biology, A1-Jumum University College, Umm A1-Qura University, Makkah, Saudi Arabia
The two-component system (TCS) genes are involved in a wide range of physiological processes in prokaryotes and eukaryotes. In plants, the TCS elements help in a variety of functions, including cell proliferation, response to abiotic and biotic stresses, leaf senescence, nutritional signaling, and division of chloroplasts. Three different kinds of proteins make up the TCS system in plants. These are known as HKs (histidine kinases), HPs (histidine phosphotransfer), and RRs (response regulators). We investigated the genome of Gossypium raimondii and discovered a total of 59 GrTCS candidates, which include 23 members of the HK family, 8 members of the HP family, and 28 members of the RR family. RR candidates are further classified as type-A (6 members), type-B (11 members), type-C (2 members), and pseudo-RRs (9 members). The GrTCS genes were analyzed in comparison with the TCS components of other plant species such as Arabidopsis thaliana, Cicer arietinum, Sorghum bicolor, Glycine max, and Oryza sativa. This analysis revealed both conservation and changes in their structures. We identified 5 pairs of GrTCS syntenic homologs in the G. raimondii genome. All 59 TCS genes in G. raimondii are located on all thirteen chromosomes. The GrTCS promoter regions have several cis-regulatory elements, which function as switches and respond to a wide variety of abiotic stresses. RNA-seq and real-time qPCR analysis showed that the majority of GrTCS genes are differentially regulated in response to salt and cold stress. 3D structures of GrTCS proteins were predicted to reveal the specific function. GrTCSs were docked with abscisic acid to assess their binding interactions. This research establishes the groundwork for future functional studies of TCS elements in G. raimondii, which will further focus on stress resistance and overall development.
Introduction
The two-component system (TCS) was discovered as the crucial signaling transduction system performing key processes both in eukaryotes and prokaryotes, including plants (Mizuno, 1997; Gao and Stock, 2009). These physiological processes involve cell division, stress tolerance, chloroplast division, nutrient signaling, and leaf senescence (Hwang and Sheen, 2001; Urao et al., 2001). A number of genes assist in these processes but three types of proteins have been distinguished so far. It was first studied in bacteria. This system has two main divisions, one is the histidine kinase sensor (HK), which is membrane associated and the second is a cytoplasmic response regulator (RR) with a receiver domain called Rec domain (West and Stock, 2001). The function of HK is the direct phosphorylation of its related RR when a certain signal is received from the environment. In response, the activity of RR is regulated by this phosphorylation. The sensor HK has two domains called input and transmitter domain. The input domains collect signals and in response, the transmitter domain regulates its HK activity after the autophosphorylation of a histidine residue. Autophosphorylation is a process in which a protein kinase catalyzes the transfer of a phosphate group from ATP to a hydroxyl group on one of its amino acid residues. This process results in the activation of the kinase itself and may also enable the recruitment and activation of downstream signaling molecules (Dabravolski and Isayenkov, 2023). The receiver domain present on RR transfers the phosphate group to conserved aspartic residue. Consequently, the output domain found in RR behaves as transcription factor (Stock et al., 2000). Over time, a protein family known as histidine phosphotransfer proteins (HPs) has developed a multi-step phosphorylation process in eukaryotes. It is assumed to be a link between the histidine kinase and the response regulator for phosphoryl group transfer because of the existence of an XHQXKGSSXS motif (Schaller et al., 2008; Tiwari et al., 2021). A multifaceted TCS transduction signaling system has been characterized in eukaryotic species mainly in higher plants. In plants, this system consists of three signaling elements: histidine phosphotransfer proteins (HPs), hybrid HKs, and RRs (Urao et al., 2000; Hwang et al., 2002).
Whole genome sequencing has helped to identify TCS genes in different plant species, i.e. A. thaliana, T. aestivum, O. sativa, and G. max. The role of TCS genes in these plants has also been studied, which respond to stress. The A. thaliana genome contains 8 HKs, 6 HPs, and 33 RRs (Hutchison and Kieber, 2002). Ethylene- (ERS), phytochrome- (PHY), and cytokinin-receptors are all HK subgroups. Additionally, there are three HKs (CKl2/HK5, HK1s, and CKI1) that are not known to belong to any particular group in A. thaliana. There are numerous transmembrane domains, an input and a REC domain, as well as an H residue (conserved) that includes a transmitter domain. All these makeup HK’s general structure. Few Ethylene receptors and phytochromes lacked conserved residues and motifs such as EIN4, ETR2, and ERS2 are called divergent of HKs. This is because they are unable to execute HK activation. The ERS subgroup comprises three different domains: a HisKA, a transmembrane domain that is involved in interacting with ethylene, and a GAF that participates in the interactions of protein-protein. There are two subgroups of ERS family on the basis of similarities in the sequences of their proteins such as AETR1 and AERS1. Each of the 5 candidates of the ERS family discovered in A. thaliana has a C2H4 domain that facilitates the binding of ethylene. These potential candidates have been given the names ERS1, ERS2, EIN4, ETR1, and ETR2 (Ahmad et al., 2020; Chen et al., 2020b). Phytochromes are receptors that respond to light stimulation and also regulate the growth and development of plants. The phytochrome subfamily consists of five members as PHYE, PHYD, PHYA, PHYB, and PHYC. Each of the 5 receptors has an amino- and a carboxyl-terminal domain, which are the two structural regions that make up each receptor. To enhance light absorbance and photo-reversible action, a tetrapyrrole chromophore has been covalently bonded to the domain at amino terminus. Carboxyl terminal has two domains, which are referred to as PAS (2) and HisKA (1). The transmission of signals is facilitated by these domains. The three transmembrane HKs work as cytokinin receptors known as AHK2, AHK3, and AHK4. They are detected due to the presence of CHASE domain in their structure (Hutchison and Kieber, 2002; Paik and Huq, 2019).
The HP family contains a domain known as Hpt, which is essential for transferring phosphate groups from HK to the RR family. This domain must contain the conserved pattern XHQXKGSSXS (Gupta et al., 2020). The five TCS members of A. thaliana are known as AHP5, AHP4, AHP3, AHP1, and AHP2 (all have Hpt domain). On the other hand, AHP6’s conserved motif does not contain H residue. So, it is referred to as a pseudo HP. AHP6 was referred to as a negative regulator of the cytokinin signaling pathway because it is incapable of acting as an HP. Gene structure, conserved motifs, and domains have been used to categorize members of the RR family into 3 distinct types such as type C, type B, and type A. Type A response regulators possess REC domain and function as proteins which are responsive to cytokinin. Type B response regulators have REC domain at N-terminus and a domain at C-terminus. Cytokinin is essential for the activation of type C response regulators and they possess domain architectures identical to type-B. Pseudo response regulators, also known as PRRs are characterized by the absence of a conserved residue D. They are unable to be phosphorylated. In cytokinin signaling, when A. thaliana is exposed to a stressor, the cytokinin receptors make a negative response (Grefen and Harter, 2004). AHKs collect signals and are autophosphorylated in the transmitter domain at conserved histidine residues. Afterward, the conserved aspartate residues are phosphorylated with phosphoryl groups. In the end, the transmitter domain transmits the signal to the receiver domain in the RR B-type and then act as a transcription factor. In this mechanism type-B RRs activate the type-A RRs in the nucleus. In addition, type-B RRs can also bind to different cis-elements of the promoter genes, like stress related genes and MAPKs responding to stress (Ishida et al., 2009).
In T. aestivum, 62 genes of TCS have been characterized of which 7 candidates from the HK family, 10 from the HP family, and 45 from the RR family (Gahlaut et al., 2014). In O. sativa, 37 genes of TCS have been characterized. Both families HK and HP have the same number of members (5), however, RR family has 27 candidates (Pareek et al., 2006). TCS signaling system is the crucial pathway in the stress signaling transduction, such as salt, drought, cold, temperature, and heat (Gupta et al., 2020; West and Stock, 2001; Zwack and Rashotte, 2015; Rasul et al., 2017). In A. thaliana, all HKs and HPs respond either negatively or positively to salinity, cold, and drought stress. A-type RRs are involved in the osmosis-stress response by working with abscisic acid and the response may be negative or positive (Miyata et al., 1998; Zwack and Rashotte, 2015). In soybean, the dehydration-sensitive TCS genes have been characterized (Mochida et al., 2010). The role of TCS genes in tomatoes has also been studied in response to stress (He et al., 2016a).
Cotton, an economically significant crop, serves as a model for various studies and is expected to be a common source of organic fiber (Li et al., 2014). G. raimondii (target plant) is diploid cotton species having a D-genome (Kirungu et al., 2018; Khan et al., 2020). It is used to produce fabric, industrial moisturizers, eatable oil, and chemicals that are essential for various industries and applications (Khadi et al., 2010). It is very sensitive to drought, pathogens, and inappropriate temperature, which reduces production and fiber reliability. TCS signaling transduction system plays a key role in the stress signaling of plants (Sharif et al., 2019; Azeem et al., 2022). The present research’s goal was to uncover TCS genes in G. raimondii since the detailed study and identification of stress-responding TCS genes can help to combat the susceptibility and mechanism of tolerance. A comparative phylogenetic tree was generated using the protein sequences of the G. raimondii and five other species (A. thaliana, G. max, C. arietinum, S. bicolor, and O. sativa). The study employed several methods, including motif and domain analysis, cis-element prediction, RNA-seq data expression profiling, real-time qPCR, and molecular docking analysis to better understand the role TCS genes in G. Raimondii stress response.
Material and methods
Identification of two-component system genes in Gossypium raimondii
Protein sequences of already known two-component system genes in A. thaliana were taken from the Ensemble Plants database (https://plants.ensembl.org/index.html). These sequences were used as queries to perform the NCBI-BLASTp program against the target species (G. raimondii) to identify TCS genes (Mahram and Herbordt, 2010). All identified genes were analyzed by using different databases i.e. SMART (http://smart.embl.de/) (Schultz et al., 2000), CDD (https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi) (Marchler-Bauer et al., 2015) and Pfam (https://pfam.xfam.org/) (Finn et al., 2009). The purpose of this analysis was to validate the identified sequences that must contain TCS-specific domains which are needed to perform normal TCS functions. Genomic information retrieved from NCBI for additional analysis, including protein length, chromosomal location and the number of exons and introns. Furthermore, the ExPASY ProtParam program (https://web.expasy.org/protparam/) was used to determine physiochemical characteristics such as molecular weight, isoelectric point, instability index, aliphatic index, and GRAVY (Gasteiger et al., 2005). TMHMM 2.0 (https://services.healthtech.dtu.dk/service.php?TMHMM-2.0) (Chen et al., 2003) and CELLO v.2.5 (http://cello.life.nctu.edu.tw/) (Yu et al., 2006) were used to identify transmembrane domains and subcellular localizations respectively.
Prediction of gene structure and conserved motifs
Gene Structure Display Server (GSDS) (http://gsds.gao-lab.org/) (Hu et al., 2015) was used to display the intron-exon patterns. Genomic and coding sequences of TCS genes from the target species were used for this objective. MEME program (https://meme-suite.org/meme/tools/meme) was employed to identify the particular conserved motifs of identified sequences. All default settings were used except the number of motifs which was fixed at 20 (Bailey et al., 2006).
MSA, phylogenetic, and putative promoter region analysis
To comprehend the ancestral relationship of GrTCS genes, the identified sequences were aligned with reference sequences using ClustalW (https://www.genome.jp/tools-bin/clustalw), a multiple sequence alignment program (Thompson et al., 2003). The phylogenetic study was carried out using the MEGA7 (Kumar et al., 2016) and IQ-TREE (http://iqtree.cibiv.univie.ac.at/) programs (Nguyen et al., 2015). The Neighbor-joining and Maximum Likelihood approaches have been utilized to generate the phylogenetic tree, which had 1000 bootstrap values. The iTOL (https://itol.embl.de/) (Letunic and Bork, 2007) web application was used to create an interactive representation of the phylogenetic tree. NCBI database was used to extract upstream genomic sequences (1,000 bp) of G. raimondii’s TCS genes. These sequences were input into plantCARE online program (https://bioinformatics.psb.ugent.be/webtools/plantcare/html/) to find out what possible cis-regulatory elements they had (Lescot et al., 2002).
Chromosomal localization, gene duplications, and syntenic analysis
GrTCS chromosomes map of genetic relationships was created employing TBtool (advanced Circos) (Chen et al., 2022). In addition, the NCBI database was utilized to determine the chromosomal locations of all identified genes in G. raimondii. Duplication occurrences were found by using the offline DNAsp tool. Ka and Ks were calculated to estimate the selection pressure on the duplicated GrTCS genes. We computed the divergence time using the following equation: Time = Ks/2x (x = 6.56 109–9) years (Mushtaq et al., 2021).
Differential expression analysis of TCS genes in Gossypium raimondii
To investigate the expression profiles of putative TCS genes in G. raimondii under salt and cold stress conditions, we used the NCBI-SRA database (https://www.ncbi.nlm.nih.gov/sra) (Leinonen et al., 2011) to download two transcriptomes for comparative analysis: PRJNA601953 (BioProject’s accession of salt stress) and PRJNA554555 (BioProject’s accession of salt and cold stress) (https://www.ncbi.nlm.nih.gov/assembly/GCF000003195.3/) was used to get the annotated genome with an extension of.gtf and.fna. The G. raimondii genome sequence indices were generated using Bowtie2, and the resulting pair-end clean reads were mapped to the G. raimondii genome(Langmead and Salzberg, 2012). Cufflinks tool was used to determine the expression levels of the annotated genes in the target genome (Ghosh and Chan, 2016). Each GrTCS’s normalized FPKM value was calculated. The heatmap for differently expressed genes was created using TBtool (Chen et al., 2020a).
Growth conditions, plant materials, and real-time qPCR
For the confirmation of RNA-seq based expression pattern, G. raimondii plants were grown in a climate-controlled growth chamber (16h light & 8h dark) at 30 °C and 26°C in light and dark, respectively (Azeem et al., 2022). Plants were treated with 250mM NaCl to impose salt stress (After the expansion of the first true leaf). For cold stress treatment, plants were kept at 4 °C. Leaf samples (three biological replicates) were collected at 12h after stress. The gene-specific primers were designed by using the online tool “Oligo Calculator” (http://mcb.berkeley.edu/labs/krantz/tools/oligocalc.html) and primer specificity was verified by NCBI PrimerBLAST program (https://www.ncbi.nlm.nih.gov/tools/primer-blast/. An internal control gene (GrActin) was used for normalization of expression (Sun et al., 2015). SYBR Green (iTaq Universal Super Mix) and CFX96 Touch™ Real-Time PCR Detection System were used for qRT-PCR analysis using the 2−ΔΔCT method (Livak and Schmittgen, 2001).
Structure prediction of differentially expressed genes
Structure prediction is essential to reveal the specific function of identified elements. Swiss-model (homology modeling server) (https://swissmodel.expasy.org/interactive) (Schwede et al., 2003) and I-TASSER (threading approach-based server) (https://zhanggroup.org/I-TASSER/) (Yang et al., 2014) were used in order to make predictions about the 3D structures of significantly up-regulated, down-regulated, and zero expression GrTCS genes in response to salt and cold stresses. After then, these structures validate by using another online tool Saves (https://saves.mbi.ucla.edu/). These structures were further refined by utilizing the online server Galaxy refine located on Galaxy-web (https://galaxy.seoklab.org/cgi-bin/submit.cgi?type=REFINE).
Docking of differentially expressed genes
We used Moe (Vilar et al., 2008) to find the top interacting GrTCS genes. Abscisic acid (ABA) was used as a ligand to dock against the structures of these genes. It is a hormone that aids in plant growth and development. Its structure was retrieved from PubChem online database (https://pubchem.ncbi.nlm.nih.gov/) (Kim et al., 2016).
Results
Identification of TCS elements in Gossypium raimondii
The primary goal of this research was to evaluate the TCS in G. raimondii. We found 59 genes associated with G. raimondii’s TCS using whole-genome sequencing. For this purpose, We used the BLASTp program to search for TCS genes in G. raimondii by using the known TCS elements of A. thaliana (47 members) as query sequences. All found genes were given names based on the names of A. thaliana TCS genes, and then they were categorized into three subfamilies: Histidine Kinase (23 members), Histidine Phosphotransfer (8 members), and Response Regulator proteins (28 members) (Table 1).
Table 1 A list of the TCS genes from various plants that have been identified.
HK protein family in Gossypium raimondii
G. raimondii has 23 members of HK family in its genome (Table 1). GrHKs are much more numerous than other HKs reported in non-leguminous plants: T. aestivum has 7 HKs, A. thaliana and O. sativa have 8 HKs, and S. bicolor has 13 HKs. In contrast, other leguminous plants, such as G.max (21), C. arietinum (18), and C. melo L. (17) contain a comparable amount of HK elements, suggesting their value in the plants (Table 1, Figure 1A). The three subgroups of the HK family are AHK, Ethylene receptor, and Phytochrome. AHK is composed of 11 members from the HK family: GrHK1, GrHK2.1, GrHK2.2, GrHK3.1, GrHK3.2, GrHK4.1, GrHK4.2, GrHK5.1, GrHK5.2, GrCKI1.1, and GrCKI1.2 that are larger than the model plant A. thaliana (6). They all have a HisKA domain with site of His phosphorylation, a HATPase_c domain, and a photoreceptor-like REC domain. The Cytokinin receptors GrHK2.1, GrHK2.2, GrHK3.1, GrHK3.2, GrHK4.1, and GrHK4.2 have an extra TM domain and the CHASE domain which serve as binding of Cytokinin (Figure 2).
Figure 1 Phylogenetic evaluation of HK (A), HP (B) and RR (C) proteins from A. thaliana, S. bicolor, O. sativa, C. arietinum, G. raimondii and G. max. In MEGA X, neighbor-joining methodology of 1000 bootstraps was utilized to create the tree. Particular colors are used to symbolize each species and subgroups. Candidates from GrTCS are denoted by a dashed-lines and dark blue color.
Figure 2 Symbolic representation of GrTCS protein domains.
Ethylene is a plant hormone that aids in aging, growth and development. Ethylene receptor elements are required for a variety of functions such as leaf abscission, propagation of root and shoot, fruit ripeness and softness. The ethylene response is mediated by 5 ethylene receptors in the model plant A. thaliana. G. raimondii, on the other hand, possesses seven ethylene receptor elements: GrETR1.1, GrETR1.2, GrETR2, GrERS1.1, GrERS1.2, GrERS2, and GrEIN4. G. raimondii’s ethylene receptors contain a HisKa domain, GAF (protein-protein interaction) domain, a HATPase_c domain (excluding GrEIN4), and a REC domain (except GrERS1.1 and GrERS1.2) (Figure 2).
Plants use phytochromes as photoreceptors to respond to photons and promote plant growth and development. In A. thaliana, candidates of phytochrome family have been identified as AtPHYA, AtPHYB, AtPHYC, AtPHYD, and AtPHYE. Phytochromes have a PHY domain at the N-terminus that implicated in the absorbing of light, a HisKa domain, two PAS domains that potentially engaged in transmission of signal, and a GAF domain. The structure of the sensor H proteins is comparable to that of soluble proteins like phytochromes. At the C-terminal, they have a HisKa domain for transmission of signal and at the N-terminal, there is a sensor domain. This family is termed as divergent HKs because it lacked the necessary 5 conserved motifs. They possess additional Ser/Thr kinase activity rather than HK activity. In G. raimondii, five phytochrome candidates were identified (GrPHYA1, GrPHYA2, GrPHYB, GrPHYC and GrPHYE). These members have HisKa, PHY, GAF, PAS and HATPase c domains required for light and transmission of signals, indicating that they were real photoreceptors (Figure 2).
HP protein family in Gossypium raimondii
The Histidine phospho-transfer protein family has six members in A. thaliana (AHP1-AHP6). Due to the existence of H conserved residue, the first five candidates are real Histidine phospho-transfer proteins. AHP6 is referred to as pseudo HP because it lacks the H residue required for phosphate acquisition from donor protein (Bleecker, 1999). Each AHPs candidate (excluding AHP6) has a phosphorylation pattern that is conserved (XHQXKGSSXS). The N residue in AHP6 replaces the H residue in the phosphorylation motif. GrHP1.1, GrHP1.2, GrHP3, GrHP4.1, GrHP4.2, GrHP5, GrHP6.1, and GrHP6.2 were identified as HP candidates with an HPt conserved domain in G. raimondii (Figure 1B). First six members (GrHP1.1-GrHP5) have phosphorylation motif with H residue. So, they were considered real HPs. Remaining two candidates (GrHP6.1, GrHP6.2) were labelled as pseudo HPs since their H residue replaces with N residue in phosphorylation motif.
RR protein family in Gossypium raimondii
The ultimate response to many environmental stimuli is regulated by response regulators (RR family). In G. raimondii, twenty-eight candidates of response regulator family both typical and pseudo were found. A. thaliana and O. sativa are model plants, having 33 and 22 RRs, respectively (Table 1) (Du et al., 2007). In the TCS signalling cascade, response regulators function as terminal components and they are switches which regulate the phosphorylation. Due to this key function, they catalysis the transfer of phosphoryl group to the conserved domain’s Asp residue. In A. thaliana, RRs have conserved residues such as aspartic acid (D) and lysine (K). The Rec domain contains these residues.
There are three sub-types of RR family (type-A, type-B, and type-C) on the basis of conserved domains. Type-A of response regulators possess REC domain (conserved residue D) and also contain a C-terminal extension. Type-B response regulators have two essential domains: REC domain and Myb (DNA binding) domain. Type-C response regulators exhibit a pattern comparable to type-A but lack a C-terminal extension. Some response regulators are considered as PRRs (pseudo-RRs) because they do not have D conserved residue (replaces with E residue) in REC domain and contain a CCT motif at the C-terminus. In G. raimondii, from 28 RRs six, eleven, two, and nine belong to type-A, type-B, type-C and PRRs respectively (Figure 1C).The GrRR5, GrRR6, GrRR8, GrRR9, GrRR16 and GrRR17 were discovered as type-A RR candidates in G. raimondii. These type-A GrRRs, like their counterparts (A. thaliana RRs), have a REC conserved domain (Figure 2). The type-B RR’s most typical candidates are nuclear proteins. They vary from type-A RRs in that they have a Myb (DNA-binding) domain. They have been shown to serve as transcriptional factors. G. raimondii has eleven candidates of type-B response regulators: GrRR1, GrRR2.1, GrRR2.2, GrRR10, GrRR11, GrRR12, GrRR14, GrRR18.1, GrRR18.2, GrRR21 and GrRR23. This is greater than their counterparts O. sativa L. (7) Z. mays (9) but lower than model plant A. thaliana (12). Myb and REC both domains are present in all type-B GrRR members excluding GrRR2.1 and GrRR23 (Figure 2). G. raimondii has two candidates for the type-C response regulator family (GrRR22 and GrRR24), which are similar to their counterpart A. thaliana. They contain a REC domain similar to type-A response regulators but have a short C-terminus (Figure 2). Phylogenetic study revealed that type-C RRs are not closely linked to type-A candidates and do not respond to cytokinin receptors.
In many plants, there is an additional subtype of RRs referred to as PRRs (diverging RRs). Model plant A. thaliana contains 9 PRRs with REC domain but they don’t have the DDK (conserved motif). Similarly, G. raimondii possesses nine PRRs named GrPRR1, GrPRR2.1, GrPRR2.2, GrPRR3, GrPRR5.1, GrPRR5.2, GrPRR7.1, GrPRR7.2, and GrPRR9. On the basis of their C-terminus extension, divergent RRs are classified into two types such as type-B PRRs contain two candidates (GrPRR2.1, 2.2) with Myb domain and clock PRRs possesses seven members: GrPRR1, GrPRR3, GrPRR5.1, GrPRR5.2, GrPRR7.1, GrPRR7.2, and GrPRR9. All clock PRRs have a CCT motif.
Genomic and physiochemical parameters of the GrTCS elements
All genomic and physiochemical features of fifty-nine identified GrTCS elements are exhibited in Table 2. These genes are located at one to thirteen chromosomes. There are a large number of GrTCS genes (12) on chromosome 7. They have 2-17 and 1-16 exons and introns, respectively. Protein length of GrTCS genes varied from 134 to 1275 amino acids. Physical and chemical parameters of GrTCS proteins such as molecular weight, pI, instability index, aliphatic index, and GRAVY are ranging from 15029.42 to 142735.62, 4.67 to 9.25, 28.10 to 71.06, 63.68 to 108.69 and -0.777 to 0.185, respectively. GrTCS elements are present in the nuclear, mitochondrial, and plasma membranes (Table 2).
Table 2 Nomenclating, physicochemical characteristics, and genomic features of GrTCS proteins.
Phylogenetic analysis
The purpose of this work was to investigate the evolutionary and phylogenetic relationships of G. raimondii’s TCS genes with their counterparts; a comparative tree was generated of target specie (G. raimondii), A. thaliana, O. sativa, C. arietinum, S. bicolor and G. max (Supplementary Figure 1). TCS elements were classified as HKs, HPs and RRs. HKs (Histidine kinases) were categorized into three groups: AHKs, ethylene receptors and phytochromes. AHK has six subfamilies named AHK2, AHK3, AHK4, AHK1, AHK5 and CKl1. First three subfamilies are also known as cytokinin receptors because they have additional cytokinin binding CHASE domain. These cytokinin receptors (GrHK2.1, GrHK2.2, GrHK3.1, GrHK3.2, GrHK4.1, and GrHK4.2) were discovered in G. raimondii which are similar to the AHK4, AHK3, and AHK2 found in A. thaliana. On the basis of the phylogenetic tree, these receptors were orthologues of A. thaliana’HKs. In A. thaliana, cytokinin receptors responsible for variety of functions including stress responsive, propagation of root and shoot, fruit ripeness and softness, and it is possible that G. raimondii’s cytokinin receptors have similar features.
AHK1 is a TM protein that participates in osmo-sensing and is highly exhibited in A. thaliana roots when exposed to salinity stress. In G. raimondii, GrHK5.1 and GrHK5.2 are orthologous of AHK5 (also known as CKI2). GrHK1, GrCKI1.1 and GrCKI1.2 are orthologous of AHK1 and CKI1 respectively. Phylogenetic study revealed that identified ethylene receptors GrETR1.1, GrETR1.2, GrERS1.1 and GrERS1.2 are exact homologous to AETR1 and AERS1. GrETR2 and GrERS2 are paralogues, they have orthologous relationship to AtETR2 and AtERS2. On the other hand, GrEIN4 is orthologous of A. thaliana EIN4.
Plant growth and development are controlled by phytochromes. These are light responsive elements (also referred to as photoreceptors). These receptors are responsible for signal transmission in model plant A. thaliana. In target specie, GrPHYA1, GrPHYA2, GrPHYB, GrPHYC and GrPHYE perform the same functions as in A. thaliana. G. raimondii has eight HPs, all of which are closely related to A. thaliana and O. sativa in terms of phylogeny. The phylogenetic relations between these HPs and their counterparts (A. thaliana) were used to group them. GrHP1.1 and GrHP1.2 were classified as AHP1-like. GrHP3 and GrHP5 were named AHP3-like and AHP5-like. GrHP4.1 and GrHP4.2 were classified as AHP4-like. GrHP6.1 and GrHP6.1 were grouped as AHP6-like HPs.
For the phylogenetic study of response regulators, protein sequences from different species such as A. thaliana, G. max, S. bicolor, C. arietinum, O. sativa and G. raimondii were employed. There are 4 sub-types of RR family (type-A, type-B, type-C and PRRs) on the basis of conserved domains. These individuals control the final reactions to environmental challenges. Divergent RRs are not regarded as real RRs because they don’t have the DDK (conserved motif). Because GrRRs are closely related to A. thaliana and O. sativa in terms of phylogeny. So, it was determined that they are real RRs. In plants, type-A RR is regarded as new RR family because it is not present in unicellular algae. They have been proposed to conduct several novel activities in plants. The REC domain as well as the DDK motif are present in RRs which are required for phosphate group acceptance. Cytokinin is the primary stimulator for type-A RRs, and type-B RRs play a role in cytokinin induction.
Analysis of gene structure and conserved motifs
The intron-exon pattern of genes is an important characteristic that offers knowledge about the gene’s various functions. So, we further investigated the intron-exon pattern of the GrTCS genes. The findings indicate that cytokinin receptors (GrHK1-GrHK5.2) have 11-17 exons and 10-16 introns (Figures 3A, B), which are comparable to A. thaliana cytokinin receptors (exons from 11 to 14 and introns from 10 to 13). AtCKl1 has nine introns and nine exons, whereas GrCKI1.1 and GrCKI1.2 each have six introns and seven exons. Ethylene responsive elements GrETR1.1, GrERS1.1 and GrERS1.2 were found to have five introns and six exon while GrETR2, GrERS2 and GrEIN4 have two introns and three exons. GrETR1.2 was shown to have eight exons and seven introns. This intron-exon pattern is identical to the intron-exon organisation of A. thaliana ethylene sensors. GrPHYA1, GrPHYA2, GrPHYB, GrPHYE were discovered to contain four introns and five exons among G. raimondii phytochrome family. GrPHYC was shown to contain four exons and three exons. The HP candidates GrHP1.1 and GrHP3 have 7 exons and 6 introns, GrHP1.2, GrHP4.1, GrHP4.2 and GrHP5 have 6 exons and 5 introns while GrHP6.1 and GrHP6.2 have 5 exons and 4 introns.
Figure 3 Structural features of GrTCS. (A) Phylogeny of GrTCS genes. (B) Shows the GrTCS gene structure assessment. In this section, introns, exons and UTRs are denoted by black lines, yellow and blue boxes, respectively. (C) The GrTCS conserved motifs are shown in this final section. Each specific color box represents a specific motif.
The number of introns in the G. raimondii response regulators ranged from one to fourteen while the exons number ranged from two to fifteen. GrPRR2.2 has greatest number of introns-exons (14 and 15 respectively) which is more than the introns-exons number in A. thaliana RRs. These findings revealed that identified members have a similar intron-exon pattern with their phylogenetic ties. The MEME tool was then used to estimate the twenty conserved motifs in identified TCS elements (Figure 3C). G. raimondii cytokinin receptors have 18, 17, 15, 8, and 2 conserved motifs, which are likewise seen in A. thaliana. They also have 4.5, 7, 11 and 12 additional conserved motifs. GrCKI1.1 has 2, 7, 11, 12 and 17 conserved motifs. The same motifs found in GrCKI1.2 (except 2). From these conserved motifs 2 and 17 similar to AtCKI1. Ethylene receptors have conserved motifs that are similar to cytokinin receptors and with the 1 and 14 additional motifs. G. raimondii phytochromes have four conserved motifs (3, 11, 14 and 19). Except for motif 3, these members share three motifs with A. thaliana: 11, 14, and 19. There were only two conserved motifs (11 and 16) in GrHPs.
All type-A response regulators have 3, 2 and 1 conserved motifs; GrRR16 and GrRR17 do not have 3 conserved motifs. Type-B response regulators have similar motifs to type-A, but they include two additional conserved motifs (13 and 6). Two conserved motifs (1 and 2) are found in Type-C RRs (GrRR22 and GrRR24). The PRRs have similar motifs to type-A response regulators but with the 6, 13 and 19 additional motifs. Similar motifs present in same families signify that their functions or sequences are not more diverse. Consequently, two component system is conserved according to evolutionary studies.
Genomic distribution, gene duplication, and syntenic analysis of the GrTCS elements
A syntenic study was used to identify the chromosomal gene locations and duplicated occurrences of the GrTCS elements in order to investigate their genomic distribution. All 59 detected TCS elements in G. raimondii are located at one to thirteen chromosomes. There are a large number of GrTCS genes (12) on chromosome 7, whereas chromosomes 1 and 2 have just two genes. G. raimondii HK proteins are dispersed on all thirteen chromosomes except chromosomes 4, 10 and 12. The HP proteins are found on chromosomes 6, 7, 9, 10 and 11. The RR proteins are present on all thirteen chromosomes (Supplementary Figure 2).
Duplication events of G. raimondii TCS gene families were investigated. Duplication events offer essential material for development and assist in genome analysis to investigate the evolution of novel genes. Many plants have an increase in the frequency of duplication events (tandem or segmental). We identified syntenic homologues (5 pairs) in the G. raimondii genome for duplication events. It revealed that GrHK3.1/GrHK3.2, GrETR1.1/GrETR1.2, GrERS1.1/GrERS1.2, GrPHYA1/GrPHYA2 and GrHP6.1/GrHP6.2 are genes with the segmental duplication. Expansion of G. raimondii’s TCS proteins occurs mostly by segmental duplications. This pattern is similar to several plants like A. thaliana and G. max. We computed both synonymous and non synonymous substitutions (Ks and Ka) of GrTCS proteins duplications, as well as their ratio (Ka/Ks). The Ks values have been employed to estimate the time required for duplication. GrTCSs Segment duplicates had Ks ranging 1.483-3.968. As a result, the diverging period spanned 113.0182927-302.4085366 MYA (Table 3).
Table 3 G. raimondii duplicated gene pairs displaying the computing Ka/Ks, Ka and Ks values and diverging time.
Putative promoter analysis of GrTCS genes
The Promoter region of G. raimondii’s TCS elements was studied to anticipate the cis-elements and obtain a better knowledge of their functions and transcriptional control. Numerous abiotic responsive and hormonal related cis-elements were discovered, including the TATA box and CAAT box which were found in nearly all GrTCS genes. In 26 GrTCS, TGACG and CGTCA motifs (MeJa responsive elements) have been discovered. ABRE element is involved in the abscisic acid response and was discovered in 36 of the 59 GrTCS genes. The ATCT, ATC, TCCC, AT1, GATA, GT1, TCT, CAG, and ACA motifs, as well as the I-box, Box 4, G-Box, G-box, Box-II, AE-box, MRE, chs-cMA1a, and chs-cMA2a cis-elements, were all connected to light responsiveness in GrTCS. The O2-site was found in 11 GrTCS genes and was associated with zein metabolism control.
GrTCS promoter regions also contained a TCA-element, which was associated with salicylic acid responsiveness. The GARE motif, P, and TATC boxes were found in 17 genes and were determined to be responsible for gibberellin responsiveness. In 39 GrTCS genes, ARE cis-element was present and revealed to be responsible for anaerobic induction. GrTCS found two cis-elements (TGA element and AuxRR-core) linked to auxin responsiveness. LTR and MBS elements were discovered in 12 GrTCS genes that respond to low temperature and drought respectively (Supplementary Figure 3). These findings suggest that GrTCS proteins serve key role in signaling pathways during activities of plant growth and development, as well as they were responding to different abiotic stresses.
Differential expression pattern of the TCS genes in Gossypium raimondii
The G. raimondii is known to be more tolerant to salinity stress and cold stress than the cultivated cotton species G. hirsutum. Therefore, understanding the mechanisms of stress response in G. raimondii can provide insights into improving stress tolerance in cultivated cotton. Transcriptome analysis can identify genes and pathways that are differentially regulated under combined stress conditions, and provide a better understanding of how plants respond to stress at the molecular level. This information can be used to identify key genes and pathways involved in stress response and to develop molecular markers for selecting stress-tolerant cotton cultivars. Furthermore, G. raimondii is a valuable model system for studying the molecular basis of stress tolerance in plants. As a diploid species with a smaller genome size than G. hirsutum, it is easier to study the functional genomics of G. raimondii. Transcriptome analysis can provide a foundation for further functional analysis of genes and pathways involved in stress tolerance in cotton.
The NCBI server’s Sequence Read Archive (SRA) Database was utilized to obtain publically accessible RNA sequencing data in order to investigate the differential expression of 59 identified GrTCS genes in response to cold and salt stresses (Supplementary Table 2). It was observed that GrPRR2.1, GrPRR2.2, GrPRR9, GrPHYE, and GrHK1 were highly upregulated in reponse to both stresses (salt and cold). On the contrary, GrRR1, GrRR5, GrRR8, and GrRR10 were highly downregulated in response to both stresses. Meanwhile, GrPRR5.1, GrPRR5.2, and GrPRR7.2 were upregulated in response to cold stress and downregulated in response to salt stress. Similarly, GrHP4.2, GrHK5.2, and GrRR11 were upregulated in response to salt sress and downregulated in response to cold stress (Figure 4). These findings suggest that GrTCS genes may serve key role in responding to different abiotic stresses.
Figure 4 RNA-seq based expression profile of GrTCS genes. The heatmap is graphical representation of GrTCSs response against salt and cold stresses. Dark red, blue and light yellow color indicates the up, down regulation and non-responsive expression, respectively.
A relative real-time qPCR analysis was carried out to confirm the expression profiles of highly differentially expressed genes obtained from RNA-seq data (Figure 5). The expression profiles correlated with RNA-seq data with some variations. The GrHK1, GrPHYE, GrERS1.1, GrPRR2.1 and GrPRR2.2 were significantly upregulated in response to both stresses. Similarly, GrRR1 and GrRR5 were downregulated in response to both stresses. The transcripts of GrPRR9, GrPRR5.2 were upregulated and GrHP4.2, GrRR10, GrRR11, GrPRR7.2 were downregulated in response to cold stress and vice versa for salt stress. Interestingly, no variation of expression was observed for GrHK5.2, GrRR8 and GrPRR5.1 for both stresses (Figure 6).
Figure 5 Relative Real-time qPCR based expression profiling of GrTCs members.
Figure 6 Docking of seven GrTCS proteins with abcisic acid (plant-hormone). The protein structures are indicated by dark blue color, while the ligand is indicated by dark red color (abcisic acid). The conserved residues are likewise labeled in blue in this figure.
3D structure prediction of differentially expressed GrTCS genes
Proteins are tough to explain because these are complicated chemical structures that include a vast range of winding topologies as well as atoms that may take on a variety of different configurations. During this investigation, some accurate predictions about the 3D-structures of 7 different TCS proteins including GrHK4.1, GrHP4.2, GrPHYA2, GrPRR5.1, GrPRR5.2, GrPRR7.1 and GrPRR9 were made (Supplementary Figure 4). These genes were selected for further analysis based on whether they showed high upregulation, downregulation, or no expression at all in response to cold or salt stress. GrHK4.1 expression highly decreased in response to salt stress and increased to some extent against cold stress. The expression of GrHP4.2 was shown to be significantly upregulated in reaction to salt stress, whereas it was found to be almost decreased in response to cold stress. GrPHYA2 expression was found to be significantly reduced in response to cold stress, and it was also shown to be reduced, although to a lesser level, in response to salt stress. In reaction to salt stress, GrPRR5.1 expression was found to be significantly reduced, but it was shown to be highly elevated in response to cold stress. GrPRR5.2 expression highly increased in response to cold and decreased to some extent against salt stress. The expression of GrPRR7.1 was shown to significantly decrease in reaction to cold and has almost zero expression in response to salt stress. In reaction to both salt and cold stresses, GrPRR9 expression was found to be significantly elevated. The final structures exhibited favorable stereochemical properties. GrHK4.1, GrHP4.2, GrPHYA2, GrPRR5.1, GrPRR5.2, GrPRR7.1, and GRPRR9 genes exhibited strong overall quality. In the predicted models, helices are shown in yellow, whereas sheets and strands are shown in green. Proteins that are members of the same family often possess a structural symmetry that is quite identical. The structures of the proteins belonging to the RR family are almost identical.
Docking of differentially expressed GrTCS genes
The decision to use abscisic acid (ABA) as the ligand for docking studies with TCS genes in G. raimondii was based on the known role of ABA in plant responses to abiotic stress, particularly salinity and cold stress. ABA is a key phytohormone involved in regulating plant responses to environmental stresses, including water deficit, salt stress, and cold stress. ABA signaling pathways are known to interact with two-component systems in plants, making ABA a relevant ligand for studying these interactions. Furthermore, the use of ABA as the ligand in docking studies was motivated by the goal of understanding the specific role of ABA in regulating the responses of TCS genes to salinity and cold stress. By using ABA as the ligand, the study aimed to identify specific interactions between ABA and TCS genes that are important for mediating stress responses. It is also important to note that the selection of ABA as the ligand does not preclude the possibility of studying other plant hormones and their interactions with TCS genes in future research. However, by focusing on ABA in this study, we were able to gain a more detailed understanding of the specific role of ABA in regulating stress responses in G. raimondii.
Molecular Docking (MD) of abscisic acid was conducted with seven different GrTCS proteins to investigate the binding interactions (Figure 6). The findings of the Molecular Docking demonstrated that all of the GrTCSs had a strong capacity to bind. Docking complexes possess the best score and RMSD value possible, which is a range from -13.3339338 (GrPRR5.1) to -6.74946165 (GrPRR9) for their S-score and from 1.04909098 (GrPRR5.1) to 2.34853554 (GrPRR9) for their RMSD. GrHK4.1’s conserved residues include MET 928, ILE 925, GLN 899, LEU 900, TRP 950, HIS 902, GLY 948, THR 947, GLU 929, and ALA 932. The bond lengths of these residues vary from 1.686 Å (THR 947) to 3.294 Å (GLN 899). GrHP4.2’s conserved residues were ARG 14, GLU 31, LEU 33, GLN 34, ASN 38, THR 37, and PRO 39 and bond lengths of these residues ranged from 1.604 Å to 3.201 Å. GrPHYA2 possessed 12 conserved residues including ARG 411, THR 412, ASP 437, VAL 439, MET 436, LYS 440, LEU 959, VAL 957, LEU 438, ASN 408, GLY 955, and ASP 958 and their bond lengths ranged from 1.482 Å to 2.863 Å. GrPRR5.1 possessed conserved residues including GLN 475, LEU 464, LEU 478, PHE 479, PRO 482, VAL 462, THR 461, MET 504, THR 467 and GLN 466. The bond lengths of these residues vary from 1.480 Å to 2.758 Å. GrPRR5.2 possessed conserved residues including GLN 550, THR 584, VAL 565, ARG 587, SER 560, ASN 547, GLU 551, ALA 542 and TYR 537. Their bond lengths ranged from 1.260 Å to 2.951 Å. GrPRR7.1’s conserved residues include SER 250, SER 249, ARG 321, GLU 231, THR 247, SER 245 and ASN 232. The bond lengths of these residues vary from 1.306 Å to 3.318 Å. GrPRR9 contained 12 conserved residues such as TYR 490, GLY 491, SER 527, GLN 526, GLN 503, ASP 524, LEU 532, ALA 563, HIS 539, THR 559, TYR 525, and GLN 506. Their bond lengths ranged from 1.502 Å to 2.669 Å (Supplementary Table 1).
Discussion
The growth and production of plants can be hindered due to their susceptibility to various biotic and abiotic stressors. Plants with their immobility are not able to resist these circumstances. As a result, they have developed several signalling systems to survive (Gupta et al., 2020; Rasul et al., 2017). The two component system (TCS) is essential for signalling, and contributes to the development of plants with enhanced characteristics like stress tolerance (Gupta et al., 2020). Members of TCS have been studied in several plant species, but in target plant (G. raimondii), the diversification of TCS gene family is never revealed. We used multiple analyses to identify and validate TCS genes in the current study.
TCS genes have been found in several plant species, such as T. aestivum, A. thaliana, O. sativa, S. bicolor, G.max, C. arietinum, and C. melo L. (Liu et al., 2020) G. raimondii has 59 genes of TCS in its genome that is more than TCS members of A. thaliana, O. sativa, and S. bicolor (47, 37 and 37, respectively) but less than G. max, T. aestivum, C. arietinum and C. melo L. (98, 62, 51 and 51, respectively). TCS genes are scattered on various chromosomes in A. thaliana, Brassica rapa, Zizania latifolia, and Cucumis melo L., and these genes are duplicated segmentally and tandemly. The majority of segmentally duplicated genes have been detected in B. rapa (Liu et al., 2014). There is one tandem incidence in Cucumis sativus L (He et al., 2016b). Citrullus lanatus, on the other hand, had one tandem and two segmental occurrences. There are 19 duplicated genes in Z. latifolia (He et al., 2020).
The five gene pairs (GrHK3.1/GrHK3.2, GrETR1.1/GrETR1.2, GrERS1.1/GrERS1.2, GrPHYA1/GrPHYA2, and GrHP6.1/GrHP6.2) represent segmental duplication in G. raimondii, leading to an enlargement of the GrTCS gene family. Similar segmental duplications have also been observed in other plants such as B. rapa (Liu et al., 2014), C. arietinum (Ahmad et al., 2020), and A. thaliana (Du et al., 2007). In B. rapa, 61 out of 85 TCS candidate genes were found to be segmentally duplicated. A. thaliana and C. arietinum have segmental duplication rates of approximately 35.7% and 55.6%, respectively. Solanum lycopersicum has both segmental and tandem duplications, with synonymous substitutions of 0.6-0.79 and 46-60 for the two types, respectively (He et al., 2016a). The divergence time for segmentally duplicated genes in G. raimondii ranges from 113.01 to 302.40 million years ago, whereas the divergence time for tandemly duplicated genes in S. lycopersicum ranges from 5.96 to 26.55 million years ago. Tandem duplication appears to be more common than segmental duplication, as it may confer a greater ability to respond to different stressors (Firon et al., 2012).
According to phylogenetic analysis, identified genes in the target plant were grouped into separate families and subfamilies, as in other plants with TCS systems (Arabidopsis, Rice, Sorghum bicolor, and Tomato). Candidates belonging to the three subgroups HK, RR, and HPs are found in all of the plants listed above. The functional domains that are preserved in G. raimondii and A. thaliana were used to classify the subgroups. AtCKI1 is referred to as a Hybrid HK (L) because of its role in cytokinin transduction (Ishida et al., 2009). GrCKI1.1 and GrCKI1.2 performed a similar function. The CKI2 subgroup includes two genuine HKs, GrHK5.1 and GrHK5.2. A. thaliana response regulators are characterized by the presence of a REC domain. Response regulators have a comparable domain in the target plant G. raimondii. The presence of a conserved domain (CHASE) places C. arietinum cytokinin receptors in a distinct clade in phylogenetic analysis. Likewise, G. raimondii cytokinin receptors have an extra TM domain and the CHASE domain.
The Promoter region of the G. raimondii’s TCS elements was studied to identify the cis-elements and aided in the identification of switches engaged in the regulation of transcriptional activity of downstream candidates. The existence of many light, salicylic acid, abscisic acid, drought, hormone, stress and defense-related, low-temperature, and auxin-responsive switches were shown by our recent research findings(Ahmad et al., 2020; Zameer et al., 2021). In many other plants, such switches have already been discovered in their promoter regions. Like light, ABRE, gibberellin, and hormone responsiveness switches present in histidine kinase family of bananas (Dhar et al., 2019), while response regulators have many stress related switches. A considerable number of cis-elements such as drought and abscisic acid responsiveness elements were discovered in C. sativus L. and C. lanatus. Aside from abiotic stresses, In B. rapa, response regulators of type-A possess an additional binding site known as GARP, which may contribute to an increase in available binding sites for type-B response regulators and lead to transcriptional boosting (Liu et al., 2014). This may contribute to the availability of binding sites for type-B response regulators, leading to transcriptional boosting. Type-A response regulators were discovered to be somewhat dependent on type-B response regulators. As a result, The TCS gene family has been found to respond to both biotic and abiotic stresses, including hormonal stress.
Plant growth and development may be influenced by abiotic stressors including heat, water, and salinity. In plants, TCS genes have a function in managing the response against several abiotic challenges. The expression study of this gene family may lead to a better understanding of its role in environmental adaptation. AtHK1 functions as a strong activator against dehydration and salinity stressors in A. thaliana (Tran et al., 2007). In the present investigation, it was found that the GrTCS gene family has different expressions in different tissues. GrPHYE, GrHP4.2, GrPRR2.1, and GrPRR9 were shown to be increased in leaves against salt stress. Numerous RR candidates and only one HK candidate such as GrPRR2.1, GrPRR2.2, GrPRR5.1, GrPRR5.2, GrPRR7.2, GrPRR9, and GrPHYE were increased in response to cold stress. However, GrHK4.1, GrHK4.2, GrPRR5.1, and GrRR8 all had lower expression levels against salinity. Several HK and RR candidates, including GrHK2.2, GrHK3.2, GrPHYA2, GrRR5, GrRR14, and GrPRR7.1 were discovered to be decreased in response to cold stress. Most HP candidates showed no changes in expression when treated to cold stress. Chickpea members such as CarHK2, 3 and 4 have been exhibited in most tissues. Two potential Chickpea HK1 and HK5 genes have been found in the plant’s pods and shoots. The majority of expression of genes belonging to the CarRR family was found in floral buds. In a similar manner, candidates for rice’s HK seemed to be detected in the leaf and roots whereas HP candidates have been detected in the leaves. Response regulators have been discovered in spikelets, leaf, and stems. The melons with the greatest level of RR expressions provide evidence that these elements play a significant role in the root cytokinin signalling process (Liu et al., 2020). The Chinese cabbage had the same findings as the previous study (Liu et al., 2014).
One of the most significant challenges facing contemporary agriculture is the impact of abiotic stressors, including salinity and water deficit. A growing body of data suggests that TCS elements play a crucial role in responding to these stressors. During a recent study, 59 GrTCS elements were identified, and it was revealed that certain genes were up-regulated or down-regulated in response to cold and salt stress. Thirty-six of the 59 GrTCS elements were affected by cold stress, while 26 were affected by salt stress. The expression profile of TCS elements has a significant impact on their response to abiotic stressors.
Conclusion
In current study, we identified 59 candidates for the TCS, which were categorized into three subfamilies. We conducted in-depth analysis of protein classifications, predictions of domains, evolutionary connections, gene architectures, dispersion of genes on different chromosomes, and occurrences of duplicate genes. Our results indicate that TCS has highly conserved sequences as well as specialized domains that can play important roles in executing various functions. We found that the phylogenetic link between GrTCS elements and TCS genes from other plant species is more closely related. Notably, we observed that GrHK candidates were duplicated segmentally, which led to the duplication of genes and ultimately to their extension. Differential expression in response to cold and salt stresses revealed that GrPRR2.1, GrPRR2.2, GrPRR9, GrPHYE, and GrHK1 were highly upregulated, while GrRR1, GrRR5, GrRR8, and GrRR10 were highly downregulated in response to both stresses. Furthermore, GrPRR5.1, GrPRR5.2, and GrPRR7.2 were upregulated in response to cold stress and downregulated in response to salt stress, whereas GrHP4.2, GrHK5.2, and GrRR11 were upregulated in response to salt stress and downregulated in response to cold stress. Real-time qPCR analysis confirmed the RNA-seq data, indicating the potential role of GrTCS genes in responding to abiotic stresses. These findings provide a better understanding of the signaling cascades and the means by which G. raimondii might improve its capacity to endure the effects of different biotic and abiotic stressors.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
FA, AR, SQ: investigation, data analysis, writing-original draft preparation, writing-review & editing. MU-R: methodology. MR, MS, DH: participated in the experiment, and data analysis. visualization. NA: Methodology. FA: conceptualization. All authors contributed to the article and approved the submitted version.
Funding
Deanship of Scientific Research at Umm A1-Qura University (Grant Code: 22UQU4281560DSR10).
Acknowledgments
The authors would like to thank the Deanship of Scientific Research at Umm A1-Qura University for supporting this work by Grant Code (22UQU4281560DSR10).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2023.1138048/full#supplementary-material
Supplementary Figure 1 | Cummulative phylogenetic tree of GrTCS. proteins from A. thaliana, S. bicolor, O. sativa, C. arietinum, G. raimondii and G. max. In MEGA X, neighbor-joining methodology of 1000 bootstraps was utilized to create the tree. Particular colors are used to symbolize each species and subgroups. Candidates from GrTCS are denoted by a dashed-lines and dark blue color.
Supplementary Figure 2 | The synteny of GrTCS members. The chromosomal positions are indicated by the genes on various bar blocks in circle. The five duplicated pairs are represented by color arcs in red, green, yellow, blue, and pink.
Supplementary Figure 3 | GrTCS’ promoter evaluation. cis-regulatory components are graphically shown. These elements indicated by distinct colors, while the numbers in the bars show the frequency of elements.
Supplementary Figure 4 | Final 3D structures of seven GrTCS proteins. Different colours were used in these models’ predictions to represent various sheets, helical structures and domains.
References
Ahmad, B., Azeem, F., Ali, M. A., Nawaz, M. A., Nadeem, H., Abbas, A., et al. (2020). Genome-wide identification and expression analysis of two component system genes in cicer arietinum. Genomics 112, 1371–1383. doi: 10.1016/j.ygeno.2019.08.006
Azeem, F., Zameer, R., Rehman Rashid, M. A., Rasul, I., Ul-Allah, S., Siddique, M. H., et al. (2022). Genome-wide analysis of potassium transport genes in gossypium raimondii suggest a role of GrHAK/KUP/KT8, GrAKT2.1 and GrAKT1.1 in response to abiotic stress. Plant Physiol. Biochem. 170, 110–122. doi: 10.1016/J.PLAPHY.2021.11.038
Bailey, T. L., Williams, N., Misleh, C., Li, W. W. (2006). MEME: discovering and analyzing DNA and protein sequence motifs. Nucleic Acids Res. 34, W369–W373. doi: 10.1093/NAR/GKL198
Bleecker, A. B. (1999). Ethylene perception and signaling: An evolutionary perspective. Trends Plant Sci. 4, 269–274. doi: 10.1016/S1360-1385(99)01427-2
Chen, C., Chen, H., Zhang, Y., Thomas, H. R., Frank, M. H., He, Y., et al. (2020a). TBtools: An integrative toolkit developed for interactive analyses of big biological data. Mol. Plant 13, 1194–1202. doi: 10.1016/J.MOLP.2020.06.009
Chen, C., Wu, Y., Xia, R. (2022). A painless way to customize circos plot: From data preparation to visualization using TBtools. iMeta 1, e35. doi: 10.1002/IMT2.35
Chen, Y., Hu, G., Rodriguez, C., Liu, M., Binder, B. M., Chervin, C. (2020b). Roles of SlETR7, a newly discovered ethylene receptor, in tomato plant and fruit development. Hortic. Res. 7, 17. doi: 10.1038/S41438-020-0239-Y/41989365/41438_2020_ARTICLE_239.PDF
Chen, Y., Yu, P., Luo, J., Jiang, Y. (2003). Secreted protein prediction system combining CJ-SPHMM, TMHMM, and PSORT. Mamm. Genome 14, 859–865. doi: 10.1007/s00335-003-2296-6
Dabravolski, S. A., Isayenkov, S. V. (2023). The role of histidine kinase signalling in response to salt stress. Plant Soil., 1–15. doi: 10.1007/s11104-023-05881-3
Dhar, Y. V., Lakhwani, D., Pandey, A., Singh, S., Trivedi, P. K., Asif, M. H. (2019). Genome-wide identification and interactome analysis of members of two-component system in banana. BMC Genomics 20, 674. doi: 10.1186/s12864-019-6050-1
Du, L., Jiao, F., Chu, J., Jin, G., Chen, M., Wu, P. (2007). The two-component signal system in rice (Oryza sativa l.): A genome-wide study of cytokinin signal perception and transduction. Genomics 89, 697–707. doi: 10.1016/j.ygeno.2007.02.001
Finn, R. D., Mistry, J., Tate, J., Coggill, P., Heger, A., Pollington, J. E., et al. (2009). The pfam protein families database. Nucleic Acids Res. 38, D211–D222. doi: 10.1093/NAR/GKP985
Firon, N., Pressman, E., Meir, S., Khoury, R., Altahan, L. (2012). Ethylene is involved in maintaining tomato (Solanum lycopersicum) pollen quality under heat-stress conditions. AoB PLANTS 2012, 1–9. doi: 10.1093/aobpla/pls024
Gahlaut, V., Mathur, S., Dhariwal, R., Khurana, J. P., Tyagi, A. K., Balyan, H. S., et al. (2014). A multi-step phosphorelay two-component system impacts on tolerance against dehydration stress in common wheat. Funct. Integr. Genomics 14, 707–716. doi: 10.1007/S10142-014-0398-8/FIGURES/3
Gao, R., Stock, A. M. (2009). Biological insights from structures of two-component proteins. Annu. Rev. Microbiol. 63, 133. doi: 10.1146/ANNUREV.MICRO.091208.073214
Gasteiger, E., Hoogland, C., Gattiker, A., Wilkins, M. R., Appel, R. D., Bairoch, A. (2005). Protein identification and analysis tools on the ExPASy server. Proteomics Protoc. Handb. 571–607. doi: 10.1385/1-59259-890-0:571
Ghosh, S., Chan, C. K. K. (2016). Analysis of RNA-seq data using TopHat and cufflinks. Methods Mol. Biol. 1374, 339–361. doi: 10.1007/978-1-4939-3167-5_18
Grefen, C., Harter, K. (2004). Plant two-component systems: Principles, functions, complexity and cross talk. Planta 219, 733–742. doi: 10.1007/s00425-004-1316-4
Gupta, P., Yadav, C., Singh, D., Nongpiur, R. C., Singla-Pareek, S. L., Pareek, A. (2020). Two component mediated stress signaling in plants. Protein Kinases Stress Signal. Plants 649, 1–19. doi: 10.1002/9781119541578.CH1
He, Y., Liu, X., Ye, L., Pan, C., Chen, L., Zou, T., et al. (2016a). Genome-wide identification and expression analysis of two-component system genes in tomato. Int. J. Mol. Sci. 17. doi: 10.3390/ijms17081204
He, Y., Liu, X., Zou, T., Pan, C., Qin, L., Chen, L., et al. (2016b). Genome-wide identification of two-component system genes in cucurbitaceae crops and expression profiling analyses in cucumber. Front. Plant Sci. 7. doi: 10.3389/fpls.2016.00899
He, L., Zhang, F., Wu, X., Hu, Y., Dong, L., Dewitte, W., et al. (2020). Genome-wide characterization and expression of two-component system genes in cytokinin-regulated gall formation in zizania latifolia. Plants 9, 1–26. doi: 10.3390/plants9111409
Hu, B., Jin, J., Guo, A.-Y., Zhang, H., Luo, J., Gao, G. (2015). GSDS 2.0: an upgraded gene feature visualization server. Bioinformatics 31, 1296–1297. doi: 10.1093/bioinformatics/btu817
Hutchison, C. E., Kieber, J. J. (2002). Cytokinin signaling in arabidopsis. Plant Cell 14 Suppl, S47–S59. doi: 10.1105/tpc.010444
Hwang, I., Sheen, J. (2001). Two-component circuitry in arabidopsis cytokinin signal transduction. Nature 413, 383–389. doi: 10.1038/35096500
Hwang, I., Chen, H.-C., Sheen, J. (2002). Two-component signal transduction pathways in arabidopsis. Plant Physiol. 129, 500–515. doi: 10.1104/pp.005504
Ishida, K., Niwa, Y., Yamashino, T., Mizuno, T. (2009). A genome-wide compilation of the two-component systems in lotus japonicus. DNA Res. 16, 237–247. doi: 10.1093/dnares/dsp012
Khadi, B. M., Santhy, V., Yadav, M. S. (2010). Cotton: An introduction. Biotechnol. Agric. For. 65, 1–14. doi: 10.1007/978-3-642-04796-1_1/COVER
Khan, M. A., Wahid, A., Ahmad, M., Tahir, M. T., Ahmed, M., Ahmad, S., et al. (2020). “World cotton production and consumption: An overview,” in Cotton production and uses. Eds. Ahmad, S., Hasanuzzaman, M. (Singapore: Springer Singapore), 1–7. doi: 10.1007/978-981-15-1472-2_1
Kim, S., Thiessen, P. A., Bolton, E. E., Chen, J., Fu, G., Gindulyte, A., et al. (2016). PubChem substance and compound databases. Nucleic Acids Res. 44, D1202–D1213. doi: 10.1093/nar/gkv951
Kirungu, J. N., Deng, Y., Cai, X., Magwanga, R. O., Zhou, Z., Wang, X., et al. (2018). Simple sequence repeat (SSR) genetic linkage map of d genome diploid cotton derived from an interspecific cross between gossypium davidsonii and gossypium klotzschianum. Int. J. Mol. Sci. 19, 204. doi: 10.3390/IJMS19010204
Kumar, S., Stecher, G., Tamura, K. (2016). MEGA7 : Molecular evolutionary genetics analysis version 7 . 0 for bigger datasets brief communication. Mol. Biol. Evol. 33, 1870–1874. doi: 10.1093/molbev/msw054
Langmead, B., Salzberg, S. L. (2012). Fast gapped-read alignment with bowtie 2. Nat. Methods 94 (9), 357–359. doi: 10.1038/nmeth.1923
Leinonen, R., Sugawara, H., Shumway, M., Collaboration, on behalf of the I. N. S. D (2011). The sequence read archive. Nucleic Acids Res. 39, D19–D21. doi: 10.1093/NAR/GKQ1019
Lescot, M., Déhais, P., Thijs, G., Marchal, K., Moreau, Y., Van de Peer, Y., et al. (2002). PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 30, 325–327. doi: 10.1093/nar/30.1.325
Letunic, I., Bork, P. (2007). Interactive tree of life (iTOL): An online tool for phylogenetic tree display and annotation. Bioinformatics 23, 127–128. doi: 10.1093/BIOINFORMATICS/BTL529
Li, Q., Xiao, G., Zhu, Y.-X. (2014). Single-nucleotide resolution mapping of the gossypium raimondii transcriptome reveals a new mechanism for alternative splicing of introns. Mol. Plant 7, 829–840. doi: 10.1093/mp/sst175
Liu, P., Wang, S., Wang, X., Yang, X., Li, Q., Wang, C., et al. (2020). Genome-wide characterization of two-component system (TCS) genes in melon (Cucumis melo l.). Plant Physiol. Biochem. 151, 197–213. doi: 10.1016/j.plaphy.2020.03.017
Liu, Z., Zhang, M., Kong, L., Lv, Y., Zou, M., Lu, G., et al. (2014). Genome-wide identification, phylogeny, duplication, and expression analyses of two-component system genes in Chinese cabbage (Brassica rapa ssp. pekinensis). DNA Res. 21, 379–396. doi: 10.1093/dnares/dsu004
Livak, K. J., Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2–ΔΔCT method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262
Mahram, A., Herbordt, M. C. (2010). Fast and accurate NCBI BLASTP: Acceleration with multiphase FPGA-based prefiltering. Proc. Int. Conf. Supercomput. 73–82. doi: 10.1145/1810085.1810099
Marchler-Bauer, A., Derbyshire, M. K., Gonzales, N. R., Lu, S., Chitsaz, F., Geer, L. Y., et al. (2015). CDD: NCBI’s conserved domain database. Nucleic Acids Res. 43, D222–D226. doi: 10.1093/NAR/GKU1221
Miyata, S. I., Urao, T., Yamaguchi-Shinozaki, K., Shinozaki, K. (1998). Characterization of genes for two-component phosphorelay mediators with a single HPt domain in arabidopsis thaliana. FEBS Lett. 437, 11–14. doi: 10.1016/S0014-5793(98)01188-0
Mizuno, T. (1997). Compilation of all genes encoding two-component phosphotransfer signal transducers in the genome of escherichia coli. DNA Res. 4, 161–168. doi: 10.1093/DNARES/4.2.161
Mochida, K., Yoshida, T., Sakurai, T., Yamaguchi-Shinozaki, K., Shinozaki, K., Tran, L. S. P. (2010). Genome-wide analysis of two-component systems and prediction of stress-responsive two-component system members in soybean. DNA Res. 17, 303–324. doi: 10.1093/dnares/dsq021
Mushtaq, N., Munir, F., Gul, A., Amir, R., Paracha, R. Z. (2021). Genome-wide analysis, identification, evolution and genomic organization of dehydration responsive element-binding (DREB) gene family in solanum tuberosum. PeerJ 9, e11647. doi: 10.7717/PEERJ.11647/SUPP-11
Nguyen, L.-T., Schmidt, H. A., Von Haeseler, A., Minh, B. Q. (2015). IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 32, 268–274. doi: 10.1093/molbev/msu300
Paik, I., Huq, E. (2019). Plant photoreceptors: Multi-functional sensory proteins and their signaling networks. Semin. Cell Dev. Biol. 92, 114–121. doi: 10.1016/J.SEMCDB.2019.03.007
Pareek, A., Singh, A., Kumar, M., Kushwaha, H. R., Lynn, A. M., Singla-Pareek, S. L. (2006). Whole-genome analysis of oryza sativa reveals similar architecture of two-component signaling machinery with arabidopsis. Plant Physiol. 142, 380–397. doi: 10.1104/PP.106.086371
Rasul, I., Nadeem, H., Siddique, M. H., Atif, R. M., Ali, M. A., Umer, A., et al. (2017). Plants sensory-response mechanisms for salinity and heat stress. J. Anim. Plant Sci. 27, 490–502.
Schaller, G. E., Kieber, J. J., Shiu, S.-H. (2008). Two-component signaling elements and histidyl-aspartyl phosphorelays †. Arab. B. 6, e0112. doi: 10.1199/tab.0112
Schultz, J., Copley, R. R., Doerks, T., Ponting, C. P., Bork, P. (2000). SMART: A web-based tool for the study of genetically mobile domains. Nucleic Acids Res. 28, 231–234. doi: 10.1093/nar/28.1.231
Schwede, T., Kopp, J., Guex, N., Peitsch, M. C. (2003). SWISS-MODEL: an automated protein homology-modeling server. Nucleic Acids Res. 31, 3381–3385. doi: 10.1093/NAR/GKG520
Sharif, I., Aleem, S., Farooq, J., Rizwan, M., Younas, A., Sarwar, G., et al. (2019). Salinity stress in cotton: Effects, mechanism of tolerance and its management strategies. Physiol. Mol. Biol. Plants 25, 807–820. doi: 10.1007/s12298-019-00676-2
Stock, A. M., Robinson, V. L., Goudreau, P. N. (2000). Two-component signal transduction. Annu. Rev. Biochem. 69, 183–215. doi: 10.1146/annurev.biochem.69.1.183
Sun, R., He, Q., Zhang, B., Wang, Q. (2015). Selection and validation of reliable reference genes in Gossypium raimondii. Biotechnol. Lett. 37, 1483–1493. doi: 10.1007/S10529-015-1810-8/TABLES/5
Thompson, J. D., Gibson, T. J., Higgins, D. G. (2003). Multiple sequence alignment using ClustalW and ClustalX. Curr. Protoc. Bioinforma. 00, 2.3.1–2.3.22. doi: 10.1002/0471250953.BI0203S00
Tiwari, M., Yadav, M., Singh, B., Pandey, V., Nawaz, K., Bhatia, S. (2021). Evolutionary and functional analysis of two-component system in chickpea reveals CaRR13, a TypeB RR, as positive regulator of symbiosis. Plant Biotechnol. J. 19, 2415–2427. doi: 10.1111/PBI.13649
Tran, L.-S. P., Urao, T., Qin, F., Maruyama, K., Kakimoto, T., Shinozaki, K., et al. (2007). Functional analysis of AHK1/ATHK1 and cytokinin receptor histidine kinases in response to abscisic acid, drought, and salt stress in arabidopsis. Proc. Natl. Acad. Sci. U. S. A. 104, 20623–20628. doi: 10.1073/pnas.0706547105
Urao, T., Yamaguchi-Shinozaki, K., Shinozaki, K. (2000). Two-component systems in plant signal transduction. Trends Plant Sci. 5, 67–74. doi: 10.1016/S1360-1385(99)01542-3
Urao, T., Yamaguchi-Shinozaki, K., Shinozaki, K. (2001). Plant histidine kinases: An emerging picture of two-component signal transduction in hormone and environmental responses. Sci. STKE 2001, 18. doi: 10.1126/STKE.2001.109.RE18
Vilar, S., Cozza, G., Moro, S. (2008). Medicinal chemistry and the molecular operating environment (MOE): Application of QSAR and molecular docking to drug discovery. Curr. Top. Med. Chem. 8, 1555–1572. doi: 10.2174/156802608786786624
West, A. H., Stock, A. M. (2001). Histidine kinases and response regulator proteins in two-component signaling systems. Trends Biochem. Sci. 26, 369–376. doi: 10.1016/S0968-0004(01)01852-7
Yang, J., Yan, R., Roy, A., Xu, D., Poisson, J., Zhang, Y. (2014). The I-TASSER suite: protein structure and function prediction. Nat. Methods 121 (12), 7–8. doi: 10.1038/nmeth.3213
Yu, C. S., Chen, Y. C., Lu, C. H., Hwang, J. K. (2006). Prediction of protein subcellular localization. Proteins Struct. Funct. Bioinforma. 64, 643–651. doi: 10.1002/PROT.21018
Zameer, R., Sadaqat, M., Fatima, K., Fiaz, S., Rasul, S., Zafar, H., et al. (2021). Two-component system genes in sorghum bicolor: Genome-wide identification and expression profiling in response to environmental stresses. Front. Genet. 12. doi: 10.3389/fgene.2021.794305
Keywords: two-component system, HKS, HPS, RRs, differential expression, 3D structures, molecular-docking, ABA
Citation: Rasool A, Azeem F, Ur-Rahman M, Rizwan M, Hussnain Siddique M, Bay DH, Binothman N, Al Kashgry NAT and Qari SH (2023) Omics-assisted characterization of two-component system genes from Gossypium Raimondii in response to salinity and molecular interaction with abscisic acid. Front. Plant Sci. 14:1138048. doi: 10.3389/fpls.2023.1138048
Received: 05 January 2023; Accepted: 09 March 2023;
Published: 31 March 2023.
Edited by:
Sajid Masood, Bahauddin Zakariya University, PakistanReviewed by:
Vidya Niranjan, R.V. College of Engineering (RVCE), IndiaMuhammad Tahir, University of Minnesota Twin Cities, United States
Copyright © 2023 Rasool, Azeem, Ur-Rahman, Rizwan, Hussnain Siddique, Bay, Binothman, Al Kashgry and Qari. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Farrukh Azeem, YXpldWFmQGhvdG1haWwuY29t; Sameer H. Qari, c2hxYXJpQHVxdS5lZHUuc2E=