Genome-wide identification of GhRLCK-VII subfamily genes in Gossypium hirsutum and investigation of their functions in resistance to Verticillium wilt

Background The receptor-like cytoplasmic kinases subfamily VII (RLCK-VII) is critical in regulating plant growth, development, and pattern-triggered immunity. However, a comprehensive exploration of these genes in the allotetraploid Gossypium hirsutum is still lacking. This study aimed to identify RLCK-VII genes in G. hirsutum and investigate their evolutionary history, structural features, expression patterns, and role in plant defense. Results Seventy-two RLCK-VII genes in the G. hirsutum genome were unveiled and classified into nine groups following their phylogenetic analysis with Arabidopsis thaliana. Group VII-1 was the largest, accounting for 28%, while Groups VII-2 and VII-3 had only one member each. The analysis using MCScanX revealed that these 72 genes formed 166 collinear gene pairs and were resided on 26 chromosomes of G. hirsutum, suggesting that they were derived from whole genome segmental duplication events. Their calculated Ka/Ks values were below one, implying the occurrence of purification selection during the evolution and inhibition of gene function differentiation/loss. All members of the RLCK-VII subfamily possessed two conserved domains, PKinase-Tyr and PKinase, and several conserved PBS1 kinase subdomains, individually included in one of the ten motifs identified using MEME. The RNA-Seq results showed that RLCK-VII genes exhibited different spatiotemporal expression, indicating their involvement in cotton growth, development, and defense responses to Verticillium dahliae. The transcription patterns of RLCK-VII genes found by RNA-Seq were further validated using qRT-PCR assays after inoculating “20B12” (cotton cultivar) with “V991” (V. dahliae). The virus-induced gene silencing (VIGS) assays uncovered that two RLCK-VII genes (Gohir.A13G227248 and Gohir.A10G219900) were essential to G. hirsutum resistance to Verticillium wilt. Conclusions These observations offer valuable insight into the attributes and roles of RLCK-VII genes in G. hirsutum, potentially enable the breeding of new cotton cultivars with enhanced resistance to Verticillium wilt. Supplementary Information The online version contains supplementary material available at 10.1186/s12870-023-04435-0.


Background
Verticillium wilt (VW) is a destructive soil-borne plant illness caused by Verticillium dahliae infection, mainly damaging their vascular bundle system.The infection can cause cotton leaves to turn yellow, wilt, and even die, decreasing lint output and fiber quality and causing significant economic losses worldwide [1].Developing new cultivars resistant to V. dahliae is the most productive and earth-friendly mean to manage the disease and support sustainable cotton production.However, no stable and durable cotton cultivars resistant to VW have yet been successfully developed.Therefore, identifying and characterizing genes associated with VW resistance using bioinformatic analysis and expression pattern analysis followed by function validation and genetic transformation will facilitate the creation of new cotton germplasms and cultivars with disease resistance.
Receptor-like kinase (RLK) is an important transmembrane protein in plants that made up of an extracellular receptor domain for signal perception, a single transmembrane region, an intracellular near-membrane domain, and cytoplasmic kinase domains, i.e., Pkinase (pfam00069) and Pkinase-Tyr (pfam07714), for signal transduction [2,3].The receptor-like cytoplasmic kinase (RLCK) is a special protein kinase family in the RLK superfamily and loses an extracellular signal sequence (ectodomain) and a transmembrane domain.Most RLCK proteins contain Ser/Thr intracellular kinase domains, and a minority of RLCK proteins also include other domains, such as U-box, leucine-rich repeat (LRR), pentatricopeptide repeat (PPR), WD40, and epidermal growth factor (EGF) [4].Arabidopsis and rice have 147 AtRLCKs and 379 OsRLCKs members identified, respectively [4].RLCKs in Arabidopsis are classified into 19 subfamilies based on phylogenetic relationships, referring to RLCK-II, RLCK-IV-RLCK-XVI, etc. [5].Earlier explorations have uncovered that RLCK family genes are essential to plant growth, development, signal transduction, and response to biological and abiotic stresses [2,4,6,7].
Remarkably, Arabidopsis RLCK subfamily VII genes have been comprehensively studied.A systematical identification indicated that Arabidopsis RLCK-VII subfamily has 46 members, including Botrytis-induced kinase 1 (BIK1), AvrPphB SUSCEPTIBLE 1 (PBS1), PBS1-Like1 (PBL1)-PBL43, and constitutive differential growth 1 (CDG1).Many RLCK-VII members play important roles in pattern-triggered immune (PTI) signaling [8].Rao et al. [8] uncovered their shared and distinct functions in plant development and receptor kinases-mediated immunity.For example, BIK1, a member of the RLCK-VII subfamily of A. thaliana, effectively improves the immune response against Botrytis cinerea infection [9].Latest studies have unveiled that BIK1 contributes to the pattern-triggered immunity (PTI) signaling pathway and co-activates PTI and effector-triggered immunity (ETI) signaling pathways to enhance plant immune responses.BIK1 phosphorylates RESPIRATORY BURST OXIDASE HOMOLOGUE (Rboh) D, a key molecular node that links PTI and ETI, leading to the formation of reactive oxygen species (ROS) that trigger plant immune responses, ultimately enhancing the plant immune response [2,10,11].PBS1 and PBL1, which also belong to the RLCK-VII subfamily, can phosphorylate Rboh to produce ROS after the exogenous application of FLG22.These findings indicate that PBS1, PBL1, and BIK1 have redundant functions in the FLG22-mediated PTI immune response [12].Arabidopsis receptor-interacting protein kinase (RIPK), i.e., AtPBL14, a member of the RLCK-VII subfamily, is critical to plant immune responses.AtRIPK activates AtRPMI receptor protein and phosphorylates AtRIN4 protein to improve plants' resistance to Pseudomonas syringae pv.tomato DC3000 (Pst DC3000).The ripk mutants in A. thaliana reduce RIN4 phosphorylation and weaken the AtRIN4-mediated immune defense response [13].While initially implicated in plant immune responses, RIPK mutants display a notable phenotype of reduced and shortened root hair.Studies have demonstrated that RIPK affects plant root hair development [12].The CaPIK1 gene, a member of the RLCK-VII subfamily with close phylogenetic relationship to AtPBL5 or AtPBL6, is highly expressed in pepper, promoting ROS accumulation and enhancing immune ability against Xanthomonas campestris pv.vesicatoria.PIK1 gene silencing reduces resistance to pathogens [14].
Increasing reports on RLCK-VII subfamily genes in A. thaliana, rice, and other plants have shown that RLCK-VII genes may also participate in the response of G. hirsutum to V. dahliae.However, only a few studies have reported on RLK superfamily's role in cotton, such as GbRLK, which is essential in regulating the drought and high salinity stresses pathway [15].Dou et al. [16] identified 29 wall-associated kinases (WAKs) in G. hirsutum (GhWAK) genes, a group of the RLKs, uncovering the potential significance of GhWAK genes in cotton fiber cell development.Zhang et al. [17] identified 58, 66, and 99 WAK/WAKL genes in G. arboreum, G. raimondii, and G. hirsutum, respectively.Expression profiling analyses have shown that abiotic stresses induce GhWAK/ WAKL expression.Furthermore, Feng et al. [18] unraveled a WAK-like kinase (GhWAKL) from upland cotton (G.hirsutum) and noted an elevated GhWAKL expression following V. dahliae infection and in response to salicylic acid (SA), implying that GhWAKL could serve as a promising molecular target for enhancing VW resistance in cotton.Other members, such as cysteine-rich receptor-like-kinases (CRKs), lectin receptor-like kinase (GhLecRK-2), and leucine-rich repeat-RLKs (LRR-RLKs), also contribute to the defense response [19][20][21].Unfortunately, there is no reports specifically addressing RLCK-VII subfamily genes in G. hirsutum.
While studies of the expression patterns and regulation of RLCK-VII subfamily genes are crucial for genetic improvement and developing resistant cotton cultivars, a systematic investigation of RLCK-VII members in G. hirsutum has yet to be reported.Therefore, this study identified and characterized RLCK-VII subfamily members based on genome assembly and annotation combined with the transcriptome sequencing (RNA-Seq) dataset of G. hirsutum.Subsequently, we validated the expression profile using qRT-PCR and verified gene functions using virus-induced gene silencing (VIGS) assays.These findings offer novel insights into the molecular regulatory system of cotton defense responses against V. dahliae and may facilitate the breeding of resistant cotton cultivars.

Identification and classification of RLCK-VII genes in G. hirsutum based on the phylogenetic tree
We identified genes encoding 4,513 proteins with the two conserved pfam00069 and pfam07714 domains in the G. hirsutum genome by using HMMER v3.3.1 software.Additionally, we found 949 proteins with an e-value < 1e − 5 by employing local BLASTP with the 46 amino acid sequences of A. thaliana RLCK-VII proteins as the queries.Because all these proteins were also included in the 4,513 proteins identified using HMMER software, they were further checked using SMART, NCBI CDD databases, DeepTMHMM, and local inter-ProScan (Table S2).We classified proteins containing kinase domains but lacking transmembrane domains and ectodomains as members of the RLCK-VII subfamily, ultimately identifying 72 RLCK-VII genes in the G. hirsutum genome (Table 1).
To classify the RLCK-VII genes in G. hirsutum, we generated a phylogenetic tree based on the amino acid sequences of 72 GhRLCK-VII members and 46 AtRLCK-VII members (Table S3, Fig. 1).According to the classification of AtRLCK-VII members, we classified the 72 GhRLCK-VII members into nine groups, named RLCK-VII-1 to RLCK-VII-9 [8].We further divided the RLCK-VII-1 group into RLCK-VII-1-A and RLCK-VII-1-B.The quantity of members within each group varied, with RLCK-VII-1 being the largest (including 20 genes, accounting for 28%) and VII-2 and VII-3 being the smallest (each with only one member) groups.There were 12 genes in the RLCK-VII-4 group, 15 genes in the RLCK-VII-6 group, 2 genes in the RLCK-VII-7 group, 6 genes in each of the RLCK-VII-5 group and VII-8 group, and 4 genes in the RLCK-VII-9 group.In addition, 5 genes were not classified into any of the groups.The protein sequence distance (diversity) of GhRLCK-VII and AtRLCK-VII in each group ranged from 0.2945 to 0.6059.The classification of GhRLCK-VII and sequences similarity of GhRLCK-VII and AtRLCK-VII (in a model plant) provided clues for investigating the roles of GhRLCK-VII members.
The physicochemical parameters of GhRLCK-VII proteins were analyzed using ProtParam software (Table 1).Their amino acid residues were in the range from 355 (Gohir.A08G028600) to 515 (Gohir.A10G161100), with an average of 426.Forty-nine (68%) GhRLCK-VII proteins had 400−500 amino acid residues, 20 (27.8%) had fewer than 400 amino acid residues, and 3 (4.2%)with more than 500 amino acid residues.The number of amino acid residues was significantly different between GhRLCK-VII-5 and GhRLCK-VII-3 (p < 0.05).Their molecular weight ranged from 39,666.44 to 56,753.74Da, with 51 (70.8%) proteins ranging from 41,388.09 to 49,683.9Da.Most proteins in the GhRLCK-VII subfamily were basic proteins, with 86% of proteins having pI values > 7.0 and only 14% of proteins having pI values < 7.0.The protein instability index varied from 23.98 to 47.81, with a mean value of 37.1.The aliphatic index of GhRLCK-VII proteins was less than 100, in the range of 66.42-92.52,indicating that all 72 proteins had a wide range of temperature stability.The grand average of hydropathy (GRAVY) values were all negative, implying that all 72 GhRLCK-VII proteins were hydrophilic.Analysis using the SignalP server 6.0 revealed that GhRLCK-VII proteins had no signaling peptides.Subcellular localization analysis unveiled that all GhRLCK-VII proteins were intracellular proteins, with 80.5% of proteins located in the cytoplasm and 19.4% of proteins localized in the chloroplasts.
To further understand the GhRLCK-VII subfamily genes, their "exon-intron" structures were investigated, and the diagrams showed similar gene structures (Fig. 3).The number of exons was 1 to 8, with 69 (96%) genes containing 4-7 exons.Gohir.D03G125400 of the RLCK-VII-8 group contained the most (8) exons, whole Gohir.A08G141200 and Gohir.D08G162300, which were not classified into any group, each had only two exons.The number, length, and location of exons/introns were comparable among genes in the same group, indicating that the 72 genes kept roughly conserved structures during evolution, with some divergences among different groups.

Gene duplication and collinearity analysis of RLCK-VII subfamily
To further understand the evolution of GhRLCK-VII subfamily, we investigated tandem duplication and segmental duplication events of GhRLCK-VII genes using MCScanX (Fig. 4).Our findings suggested that segmental duplication was the primary driver of subfamily's expansion during evolution, with only six GhRLCK-VII genes involved in tandem duplication (Table S6).We found a total of 266 pairwise collinear relationships involving one gene (100 pairwise collinearity between 7 RLCK-VII genes and others) or two genes (166 pairwise collinearity among 65 genes) in the GhRLCK-VII subfamily.Of the 7 genes of the GhRLCK-VII subfamily involved in the 100 pairwise collinearity, none had GhRLCK-VII counterparts, suggesting that they might have been lost or undergone prolonged positive (adaptive) selection.We also observed 26, 24, and 28 pairwise collinear relationships among chromosomal regions containing genes in the RLCK-VII-1 group, in the RLCK-VII-4 group, and in the RLCK-VII-6 group, respectively.Additionally, we found 21 pairwise collinear relationships in subgenome A t and 17 pairwise colinear relationships in subgenome D t , respectively, in addition to the collinear relationships between subgenome A t and D t .Furthermore, we observed the collinearity of single genes corresponding to multiple genes in the GhRLCK-VII subfamily.These results indicated that the GhRLCK-VII subfamily underwent segmental duplication events during evolution, resulting in either one or several times of  duplication events of certain chromosomal segments.
Our collinearity analysis indicated that chromosomal segmental duplication contributed significantly to the expansion of RLCK-VII subfamily genes in G. hirsutum during evolution.
To investigate the evolutionary selection pressure on the GhRLCK-VII subfamily after chromosomal segmental duplication, the Ka/Ks ratios of collinear gene pairs was calculated using KaKs_Calculator.The results uncovered that 166 gene pairs involved in segmental duplication had Ka/Ks ratios less than 1, implying the occurrence of strict purification selection during evolution and severe inhibition of functional differentiation (Table S6).We inferred that the genes arising from segmental duplication were relatively conservative in evolution, preserving their previous protein structures and functions without losing biological activity, ultimately enhancing specific biological functions (such as stress tolerance) in plants.

Identification of cis-acting regulatory elements in GhRLCK-VII gene promoters
The gene promoter region often possesses cis-acting regulatory elements linked to its function.We found 48 cis-acting elements using the program PlantCARE in the upstream region of GhRLCK-VII genes (Table S7).We focused on 11 important elements for visualization based on the characteristics of the RLCK-VII gene subfamily, excluding general transcription regulatory elements and unknown-functional elements (Fig. 5).We observed CAAT-box, TATA-box, and light response elements in the promoter region of 72 GhRLCK-VII genes.In addition, we found TC-rich repeats, low-temperature responsiveness region, and MYB binding site (MBS) in 28 genes, salicylic acid responsiveness element and TCA element in 41 genes, ABA-responsive element in 48 genes, anaerobic induction element ARE in 65 genes, methyl jasmonate (MeJA) reaction element in 46 genes, and gibberellin induction element in 37 genes.Most genes contained multiple cis-acting elements.These findings implied the involvement of GhRLCK-VII subfamily in plant growth, development, hormone response, and resistance to

Expression profile analysis based on RNA-Seq datasets and qRT-PCR assays
To get insights into the expression patterns of GhRLCK-VII subfamily, we analyzed RNA-Seq datasets SRP166405 from various tissues of G. hirsutum and SRP192537 and SRP328396 after inoculation with V. dahliae.Our results demonstrated distinct expression patterns of GhRLCK-VII genes in different tissues and developmental stages of G. hirsutum (Fig. 6, Table S8), with Gohir.A08G161700 being highly expressed in petals, filaments, and anthers, Gohir.A09G148400 being expressed significantly higher than other genes in the bract, and Gohir.D03G125400 being generally higher expressed in all tested tissues.Moreover, Gohir.D10G212300 and Gohir.D11G073300 genes were prominently expressed in petals and anthers, signaling their possible involvement in the growth and development of corresponding tissues, respectively.Moreover, some GhRLCK-VII genes, including Gohir.A04G072800, were down-regulated during ovule and fiber development, while others, including Gohir.We also explored GhRLCK-VII gene expression in cotton inoculated with V. dahliae based on the RNA-Seq datasets SRP192537 and SRP328396.Our finding demonstrated that 72 genes were significantly up-regulated in V. dahliae-tolerant cultivar Zhongzhimian-2 by 58-fold at 24 h, 72 h, and 120 h of post-inoculation with V. dahliae (strain 2 and strain 3), but only changed 25-fold in V. dahliae-susceptible cultivar Xinluzao-36.Several genes, including Gohir.D11G079000, were differentially expressed only in Zhongzhimian-2 after V. dahliae infection, while other genes, including Gohir.D13G232700, were differentially expressed in both cultivars at some time-points.The up-regulation of Gohir.D02G100700, Gohir.A13G227264, Gohir.A13G227248, Gohir.A10G219900 and Gohir.A02G085000 and Fig. 4 Collinear relationships of RLCK-VII genes, revealing segmental duplications of RLCK-VII genes and gene-pairs involved in segmental duplication.The red-colored genes represent the collinearity between RLCK-VII and other family genes, and the unmarked red genes represent the collinearity within the RLCK-VII family downregulation of Gohir.A05G034100, Gohir.D05G035500, and Gohir.D11G079000 in Zhongzhimian-2 after V. dahliae inoculation suggested that the increased transcription of more GhRLCK-VII genes in Zhongzhimian-2 than in Xinluzao-36 contributed to the tolerance of Zhongzhimian-2 to V. dahliae (Fig. 6).
To verify the involvement of RLCK-VII subfamily genes in response to V. dahliae V991, we analyzed 9 genes from G. hirsutum "20B12" using qRT-PCR (Fig. 7).Our findings showed differential regulation of RLCK-VII subfamily genes at different time-points of post V. dahliae inoculation, with some genes being up-regulated while others being down-regulated.For example, Fig. 5 Predicted cis-elements in RLCK-VII promoters the expression levels of Gohir.A05G034100, Gohir.D05G035500, and Gohir.D11G079000 were down-regulated after inoculation, and the expression levels of Gohir.D02G100700, Gohir.A13G227264, Gohir.A13G227248, Gohir.A10G219900, and Gohir.A02G085000 genes were up-regulated after inoculation.These results were in agreement with RNA-Seq data, further backing the hypothesis that GhRLCK-VII subfamily genes are responsible for cotton resistance to V. dahliae.

Function characterization of GhRLCK-VII genes in cotton respondence to V. dahliae through VIGS
We utilized three different perspectives, including DNA (many cis-acting elements (respondence to stresses) in the 2000 bp sequences upstream), mRNA (significantly increased transcript abundance after inoculation with V. dahliae), and amino acid (phylogenetic classification and conserved domain), to infer that Gohir.A13G227248 and Gohir.A10G219900 may positively contribute to cotton resistance to V. dahliae.To further understand their function, we designed primers to amplify 342-bp and 345-bp fragments of these genes, respectively, and constructed TRV2::Gohir.A13G227248 and TRV2::Gohir.A10G219900 vectors by inserting them into the pTRV2 vector.After confirming their sequences, we transformed the constructs individually into Agrobacterium.After 14 days of agroinoculation, cotton seedlings infiltrated with TRV2::PDS exhibited consistently uniform bleaching Fig. 7 Expression pattern of RLCK VII genes in G. hirsutum "20B12" treated with V. dahliae "V991" measured using qRT-PCR.The line charts represent the gene expression trends in the RNA-Seq datasets, while the bar charts depict the gene expression level in qRT-PCR assays (an albino phenotype) in newly sprouted leaves (Fig. 8a), while plants infiltrated with TRV2::Gohir.A13G227248 and TRV2::Gohir.A10G219900 exhibited reduced much lower expression of target genes compared with in the control plants infiltrated with TRV2::GUS (P < 0.01; Fig. 8b), indicating successful silencing of the target genes.
Subsequently, we inoculated cotton seedlings with V. dahliae and observed that the chlorosis spots arose in the leaf margins and veins after 7 days of inoculation.Compared with the control, the agroinoculated plants exhibited a higher number of chlorosis spots early.By about 40 days after inoculation, the plants agroinoculated with TRV2::Gohir.A13G227248 and TRV2::Gohir.A10G219900 showed more severe symptoms, including smaller but extensive chlorosis, gradual browning, wilting, and necrosis (Fig. 8c).By 60 days after inoculation, the leaves of agroinoculated plants developed more areas of dead brown tissues surrounded by larger areas of yellowing and showed a defoliating pathotype, and the stem's vascular bundles were seriously browning.We calculated the disease index and observed that the cotton seedlings infiltrated with TRV2::Gohir.A13G227248 and TRV2::Gohir.A10G219900 were more susceptible to V. dahliae, as they showed a higher disease index than the control (P < 0.05, Fig. 8d).Additionally, the cotton seedling agroinoculated with TRV2::Gohir.A13G227248 exhibited more serious disease symptoms than the cotton seedling agroinoculated with TRV2::Gohir.A10G219900, with the disease index being 87.5% higher.Based on the combined results of qRT-PCR and the observed phenotype of VW, we inferred that Gohir.A13G227248 and Gohir.A10G219900 were essential for cotton resistance to V. dahliae.

Discussion
Accumulating evidence has uncovered that RLCK-VII subfamily are actively involved in plant immune signaling pathways and essential in plant disease resistance, tolerance, growth, development, and other biological processes.In Arabidopsis, there are 46 RLCK-VII proteins that function in pattern-triggered immune signaling and confer plant resistance to pathogens [8].However, a comprehensive research is still missing on RLCK-VII subfamily genes in G. hirsutum, which often encounters pathogen attacks.
In this study, we identified 72 GhRLCK-VII subfamily genes in G. hirsutum through bioinformatics analysis.As reported previously, most plant species have undergone large-scale whole-genome duplication events, including segmental duplication and tandem duplication, leading to a larger genome size.Allotetraploid G. hirsutum possesses much more and larger GhRLCK-VII genes than diploid A. thaliana (46 genes), indicating a positive correlation between genome and RLCK-VII subfamily size [3,5].
Chromosomal segment duplication and tandem duplication provided raw genetic materials for new gene generation (gene gains) and contributed to phenotypic changes, facilitating species evolution [22].Our collinearity analysis indicated that most duplication events resulted from chromosomal segmental duplication and all GhRLCK-VII genes were implicated in segment duplication.These results showed that following segmental duplication, the function of the newly paralogous genes remained relatively conserved during evolution, with preserved protein structures and functions.Additionally, our collinearity analysis found gene loss, pseudogene formation, or dissimilation.Dissimilation of duplicated genes may have facilitated the co-evolution of plants and pathogens, while gene loss or pseudogene formation slightly reduced gene redundancy.Overall, segmental duplication improved plants' biological function and adaptability to diverse stresses.
Ka/Ks values were computed for pairwise collinear genes after MCScanX analysis to estimate selective pressure.The values indicated that the RLCK-VII genes in G. hirsutum underwent strict purification selection after expansion and retained primitive gene function.However, for a minority of GhRLCK-VII genes that lost their collinear counterpart due to gene loss, pseudogenization (defunction), or dissimilation, the related Ka/Ks values could not be computed.
To classify the GhRLCK-VII genes, we used Arabidopsis RLCK-VII subfamily, which has been systematically investigated [8], as the reference and grouped them into nine groups according to the phylogenetic tree (Fig. 1).Our classification revealed that 5 genes did not fall into any group.Compared to Arabidopsis, the sizes of groups VII-1, VII-4, VII-5, and VII-6 were expanded 2−3 times in G. hirsutum.For example, AtPBS1 in the group VII-1 was highly homologous to Gohir.A01G021700, Gohir.A01G021800, and Gohir.D01G020200 in terms of phylogenetic relationship.In terms of amino acid sequences, motif arrangements, and conserved subdomains (Table S9), RLCK-VII members within a group (intragroup) were similar to each other but different (intergroup) from those in other groups.Moreover, all genes within a phylogenetic group possess a similar exon-intron structure, which may be evidence of their derivation from a common ancestor [23].
Our identification using local HMMER search and online SMART analyses (http://smart.embl.de)indicated that 72 members of the RLCK-VII subfamily had two conserved domains, PKinase-Tyr (pfam07714) and PKinase (pfam00069), which overlapped each other.However, local InterProScan analyses showed that every member contained only one domain, with 44 and 28 members containing pfam07714 and pfam00069, respectively.This discrepancy may be due to differences in the algorithms implemented by each program, highlighting the need to integrate results from different programs.Additionally, local InterProScan analyses showed that prosite PS50011 (protein kinase domain) was present in 72 members, PS00107 (Protein kinases ATP-binding region signature) and PS00108 (Serine/Threonine protein kinases active-site signature) were present in 63 and 62 members, respectively, and 53 members contained three prosites.Furthermore, local InterProScan analyses also showed 19 members contained a SMART00220 (catalytic domain).These abundant sites implied the diverse roles of members of the RLCK-VII subfamily.
The conserved kinase active sites create a stable milieu for RLCKs to execute their function in the cytoplasm.Previous studies indicated that mojority RLCK kinase active sites are aspartic acid (D).Fan et al. [6] showed that the kinase activity sites for RLCK-VII subfamily are 'IYRDIKASNIL' , and every subfamily in RLCK harbors specific kinase activity sites.However, we found 17 kinds of activity sites (excluding 'IYRDIKASNIL') in the GhRLCK-VII subfamily, among which the predominant sites are 'IYRDFKTSNILL' and 'IYRDFKASNILL' , which are not present in other subfamilies in the RLCK family [6].The sequence pattern of activity sites in GhRLCK-VII is 'IYRD[FL]K[TSA]SNIL' present in the aforementioned motif-3.
The GDK motif is a cleavage motif, where cleavage occurs during plant immune response, i.e., AthPBS1 is cleaved by the effector AvrPphB.We found that 36 (50%) GhRLCK-VII members harbor the 'GDK' motif.This cleavage exposes the SEMPH motif in the C-terminal loop, that is recognized and bound by the NLR protein family member RPS5, activating RPS5 and subsequently inducing a hypersensitive response (HR) [24].We found that Gohir.A01G021700 and Gohir.D01G020200 harbor an STRPH motif rather than the SEMPH recognition motif, consistent with TaPBS1 in wheat [24].Sun et al. [24] found that recognition of PBS1 by RPS5 requires a negatively charged amino acid residue at the "E" position of the SEMPH motif.Therefore, our results support the highly diverse resistance recognition motifs in PBS1 orthologs.
The examination of PBS1 orthologs in various plant species revealed that all the orthologs contained several preserved subdomains of PBS1 kinase [25].Analysis using the program PS_scan revealed the presence of several subdomains in GhRLCK-VII members, such as GxGxxG, K, E, DxxxxN, DFG, APE, DxxxxG and R, showing similar sequence patterns to 'GEGGFG', ' [ S9).This indicates the conservation and/or diversity between GhRLCK-VII members and PBS1 orthologs.Moreover, MEME analyses indicated that motif-1 contained subdomains DFG, GDK, APE, and DxxxxG; motif-2 contained subdomains K and E; motif-4 contained subdomain GxGxxG [25].These conserved kinase subdomains and motifs may be critical to directing immune response of GhRLCK-VII members, allowing them to interact with a various downstream proteins and giving rise to diverse responses depending on the effector's nature [25].
The RNA-Seq and qRT-PCR assays demonstrated that GhRLCK-VII genes exhibited diverse tissue-specific expression patterns.For instance, Gohir.A12G166000 was overexpressed in anther, while Gohir.A06G147500 was highly expressed in 25 DPA fiber.The expression of Gohir.A11G062700 was high in most tissues but decreased after inoculation with V. dahliae, implying that RLCK-VII genes may contribute to specific growth and development stages, in line with several previous investigations [2,8].
Genes in the same phylogenetic group may contribute redundancies that originated from whole genome duplication events.However, whether they differ in spatiotemporal expression patterns has not been completely investigated.We found that Gohir.A10G204300, Gohir.D10G212300, and Gohir.D13G232700 of the RLCK-VII-4 group were highly expressed in some tissues (torus, bracts, sepals, petals, filaments, anthers, ovules, and fibers), but lowly expressed in others (roots, stems, and leaves) and leaves after V. dahliae infection.Reversely, both Gohir.A05G095400 and Gohir.A13G227264 in the RLCK-VII-8 group were highly expressed in anther and 10 DPA fiber, respectively.Therefore, redundant and phylogenetically related members of the RLCK-VII subfamily may have either similar expression patterns, resulting in biological promotion, or different and specific expression patterns, reducing redundancy.
Investigating tissue-specific expression profiles offers valuable information regarding the contribution of GhRLCK-VII genes in plant growth and development.Furthermore, our understanding of the potential involvement of GhRLCK-VII genes in plant responses to biotic and abiotic stresses remains limited.Yang et al. [26] showed that Gbvdr6, an RLK gene of G. barbadense, enhanced the resistance to V. dahliae in transgenic Arabidopsis and cotton plants.GhWAK7A, a wall-associated kinases (WAKs) gene belonging to RLKs, also enforced cotton defense against V. dahliae infections [27].RLCK-VII members are generally regulated by various pathogen effectors [8].In this study, analyses of the RNA-Seq datasets (SRP192537 and SRP328396) from cotton inoculated by V. dahliae revealed that some genes were up-regulated while others were down-regulated.Our qRT-PCR assays provided convincing supplements and validation, followed by gene knockdown through VIGS and symptoms phenotyping.The expression profile after V. dahliae inoculation and VIGS assays of Gohir.A13G227248 (ortholog of AtPBL37) and Gohir.A10G219900 (ortholog of AtBIK1) indicated that GhRLCK-VII genes were involved and might play certain regulatory roles in responding to V. dahliae.Veronese et al. [9] reported that BIK1 modulated various cellular factors necessary to defense responses against pathogen infection and normal root hair growth, connecting defense response regulation with growth and development.Our results are in line with the reported role of BIK1 in Arabidopsis.

Conclusion
Our comprehensive identification revealed 72 GhRLCK-VII genes in G. hirsutum.Characterization of their DNA and amino acid sequences provides insights into their location, composition, evolution, conservation, and diversity.RNA-Seq datasets and qRT-PCR unveiled diverse spatiotemporal expression patterns, and VIGS assays further indicated the potential impacts of GhRLCK-VII genes on cotton defense responses against pathogens and on plant growth and development.Our findings open novel avenues for further studying disease response pathways and enhance our comprehension of GhRLCK-VII functions.
To uncover members of the RLCK-VII subfamily in G. hirsutum, the protein sequences of AtRLCK-VII members were utilized as queries and aligned with all protein sequences of G. hirsutum by running BLASTP locally.The HMMER v3.3.1 software was used to build the HMM profile based on the two conserved domains, and the proteins containing both conserved domains were identified using hmmsearch.The candidate protein sequences of RLCK-VII subfamily were further validated by screening the intersection of BLASTP subjects and HMMsearch results and submitted to the SMART database (http://smart.embl.de/),local InterProScan, Phobius (https://phobius.sbc.su.se/), and DeepTMHMM (https://dtu.biolib.com/DeepTMHMM) to pinpoint proteins lacking transmembrane domains and ectodomains.Subsequently, subcellular localization was projected using Plant-mPLoc (http://www.csbio.sjtu.edu.cn/bioinf/plant-multi/) [29].Proteins containing kinase domains lacking transmembrane domains and ectodomains were ultimately identified as members of the G. hirsutum RLCK-VII subfamily.
To analyze the relationship between G. hirsutum RLCK-VII members and 46 members of A. thaliana, their amino acid sequences were analyzed using the ClustalX program for multiple sequence alignment.Subsequently, a phylogenetic tree was generated using the neighborjoining (NJ) method with a bootstrap value of 1000 via MEGA software [30].Finally, based on the phylogenetic tree, G. hirsutum RLCK-VII members were categorized into several groups.

Chromosome localization, physicochemical properties, and signal peptide of GhRLCK-VII members
The chromosome length and physical location information of genes on G. hirsutum chromosomes was extracted from the GFF3 annotation file.The positions of GhRLCK-VII genes on chromosomes were graphically displayed using MapChart.The amino acids count, theoretical isoelectric points (pI), molecular weight, protein instability index, aliphatic index, and grand average of hydrophobicity (GRAVY) of GhRLCK-VII proteins were calculated using the ProtParam (http://web.expasy.org/ProtParam/)[31].Signal peptides were projected using SignalP-6.0(http://www.cbs.dtu.dk/services/SignalP)[32].

Collinearity and evolutionary relationship among GhRLCK-VII genes
To investigate gene duplication during evolution within species, the collinearity between genes and genes' physical location on chromosomes of GhRLCK-VII subfamily members was analyzed using MCScanX and visualized using Circos [35,36].Non-synonymous mutation rate (Ka), synonymous mutation rate (Ks), and Ka/Ks value of selection pressure were obtained using the KaKs_Calculator [37].

Inoculation with V. dahliae and qRT-PCR validation
Nine genes were selected to verify the expression profiles of GhRLCK-VII subfamily genes using qRT-PCR assays with specific primers (Table S1) designed using Oligo 7.60 and synthesized by the Xian Qingke Biotechnology Company.In detail, highly virulent V. dahliae strain "V991" from our laboratory was grown on potato dextrose agar (PDA) medium in the dark at 25 °C for about 10 days.After rinsed with sterile distilled water, V. dahliae was adjusted 1 × 10 6 conidia/mL.Meanwhile, V. dahliae tolerant G. hirsutum cv."20B12" was grown in an incubator.At the 3-leaf stage, the cotton seedlings were infected with V. dahliae using the dipping root method, transplanted to an incubator with new soil, and cultured for different durations [40].
The seedling roots were gathered at 0 h (CK), 24 h, 72 h, and 120 h post-inoculation with 3 or more independent plants.Total RNAs were extracted using the RNA Prep Pure Plant Kit (Tiangen Biochemical Technology (Beijing) Co., Ltd.) as per the manufacturer's instructions.
RNA concentration and integrity were measured using NanoDrop and 1.2% agarose gel electrophoresis, respectively.cDNAs were generated using the HiFi Script gDNA Removal cDNA Synthesis Kit (Kang Wei Century Company) and subjected to qRT-PCR assay in a 20-µL reaction system on a QuantStudio 7 Flex Real-time Fluorescence Quantitative System with UBQ7 as the internal control.The reaction was proceeded at pre-denaturation at 95 °C for 10 min prior to 40 cycles of 95 °C for 15 s, 60 °C for 1 min, and 72 °C for 15 s.The relative gene expression was calculated using the 2 −ΔΔCt method.

Virus-induced gene silencing (VIGS) assay and symptom phenotyping
RNAs were isolated from "20B12" leaves, converted to single-strand cDNA, and used to amplify Gohir.A13G227248 and Gohir.A10G219900 using primers with KpnI and XbaI restriction sites and overhangs designed using Oligo 7.60 software (Table S1).After treated with with KpnI and XbaI, the generated fragments were cloned into TRV2 vector.After sequencing verification, the favorable plasmid was prepared and transformed into Agrobacterium GV3101.The agrobacterium strains carrying TRV2::Gohir.A13G227248, TRV2::Gohir.A10G219900, TRV2::PDS, TRV2::GUS, and TRV1 were grown in YEB liquid medium supplemented with kanamycin, rifampicin and gentamicin at 28 °C at 200 rpm.After culturing for 10-12 h when OD 600 value reached 1.0, bacteria were precipitated, resuspended in an Agrobacterium induction buffer with 20 g/L sucrose, 5 g/L MS salts, 1.95 g/L MES, pH 5.6, and placed at 22 °C for 2 h.The TRV1 and TRV2::gene cultures were combined at a 1:1 (vol/vol) ratio and inoculated into the abaxial side of the cotyledons of 14-day-old leaf-free G. hirsutum cv."20B12" seedlings.The inoculated plants were incubated for 24 h in the dark, followed by a '16 h light/8 h dark' photoperiod.When the cotton seedlings infiltrated with TRV2::PDS showed a bleached leaf phenotype, leaf RNAs were extracted and used for qRT-PCR assay to evaluate the silencing efficiency.Meanwhile, cotton seedlings' roots were dip-inoculated with V. dahliae "V991" conidia suspension.After dip-inoculation, the cotton seedlings were transplanted and grown in an incubator with new soil to observe symptoms.At 35 days of post-dip-inoculation, the number of diseased plants and the degree of disease were recorded to calculate the disease index.

Fig. 1
Fig. 1 Phylogenetic tree depicting proteins encoded by 72 GhRLCK-VII genes in G. hirsutum genome and 46 AtRLCK-VII genes in A. thaliana genome

Fig. 2
Fig. 2 Localization of GhRLCK-VII subfamily genes on the 26 chromosomes in G. hirsutum genome

Fig. 3
Fig. 3 Phylogenetic relationship, motif architecture, gene structure (exon-intron organization), and conserved domains of GhRLCK-VII genes.(A) Phylogenetic tree.(B) Motif composition and distribution.The colored boxes indicate different conserved motifs identified using MEME.(C) Gene structure (exon-intron organization).The sequences of UTR, CDS, and PKc_like superfamily are represented in different colors

Fig. 6
Fig. 6 Expression profile of RLCK-VII genes based on three RNA-Seq datasets

Table 1
Physicochemical properties of RLCK-VII subfamily proteins in G. hirsutumS4 and Table