Identification of cell wall-associated kinases as important regulators involved in Gossypium hirsutum resistance to Verticillium dahliae
BMC Plant Biology volume 21, Article number: 220 (2021)
Verticillium wilt, caused by the soil borne fungus Verticillium dahliae, is a major threat to cotton production worldwide. An increasing number of findings indicate that WAK genes participate in plant−pathogen interactions, but their roles in cotton resistance to V. dahliae remain largely unclear.
Here, we carried out a genome-wide analysis of WAK gene family in Gossypium hirsutum that resulted in the identification of 81 putative GhWAKs, which were all predicated to be localized on plasma membrane. In which, GhWAK77 as a representative was further located in tobacco epidermal cells using transient expression of fluorescent fusion proteins. All GhWAKs could be classified into seven groups according to their diverse protein domains, indicating that they might sense different outside signals to trigger intracellular signaling pathways that were response to various environmental stresses. A lot of cis-regulatory elements were predicted in the upstream region of GhWAKs and classified into four main groups including hormones, biotic, abiotic and light. As many as 28 GhWAKs, playing a potential role in the interaction between cotton and V. dahliae, were screened out by RNA-seq and qRT-PCR. To further study the function of GhWAKs in cotton resistance to V. dahliae, VIGS technology was used to silence GhWAKs. At 20 dpi, VIGSed plants exhibited more chlorosis and wilting than the control plants. The disease indices of VIGSed plants were also significantly higher than those of the control. Furthermore, silencing of GhWAKs significantly affected the expression of JA- and SA-related marker genes, increased the spread of V. dahliae in the cotton stems, dramatically compromised V. dahliae-induced accumulation of lignin, H2O2 and NO, but enhanced POD activity.
Our study presents a comprehensive analysis on cotton WAK gene family for the first time. Expression analysis and VIGS assay provided direct evidences on GhWAKs participation in the cotton resistance to V. dahliae.
Tetraploid Gossypium hirsutum is the most widely cultivated cotton species in the world and represents an important source of natural fiber and oilseed. Verticillium wilt, caused by the soil borne fungus Verticillium dahliae, is a major threat to cotton production . Identification and characterization of genes associated with resistance is an important basis for potential understanding on the interaction between cotton and V. dahliae, which is necessary for the development of novel disease management methods and new varieties resistant to Verticillium wilt.
Plants live in a complex environment crowded with biotic stresses mainly caused by various phytopathogens and pests, and expose to abiotic stresses including cold, hot, drought and salinity. To overcome these stress challenges, plants have evolved a complex and efficient defense signaling network, which includes monitoring systems to perceive different stress-derived signals triggering specific defense responses . Cell wall, a dynamic structure surrounding plant cell, has emerged as an essential monitoring system [3, 4]. Some receptor-like kinases (RLKs) have been identified as cell wall integrity sensors that are responsible for the communication between the cell wall and cytoplasm. Typically, RLKs contain a signal peptide (SP), transmembrane (TM) domain, and cytoplasmic kinase domain. They can be classified into more than 21 subfamilies according to their diverse extracellular domains . Of which, wall-associated kinases (WAKs) are distinguished from the other RLKs by the presence of their unique extracellular epidermal growth factor (EGF)-like domains [5, 6].
In Arabidopsis thaliana, WAKs are encoded by 5 WAKs and 22 WAKLs (WAK-like genes) . So far, WAK gene family was also identified in other plants, including Oryza sativa , Brassica rapa  and Populus trichocarpa . It has been demonstrated that some WAKs are involved in plant development, abiotic and biotic stress responsiveness. Notably, most of WAKs were characterized from Arabidopsis and rice. Arabidopsis AtWAK1, the first identified WAK gene in plant, was shown to contribute to the immune response [11, 12]. A rice WAK gene, OsDEES1 (DEFECT IN EARLY EMBRYO SAC1), played a role in the regulation of early embryo sac development . OsiWAK1 (O. sativa indica WAK-1) and HvWAK1 (Hordeum vulgare WAK-1) were involved in plant root development [14, 15]. Xa4, encoding a WAK in rice, conferred race-specific durable resistance against Xanthomonas oryzae pv. oryzae by reinforcing the cell wall and increasing the production of jasmonate-isoleucine and phytoalexins . OsWAK1 (O. sativa WAK) and OsWAK25 were up-regulated by wounding and salicylic acid (SA), and their overexpression led to higher resistance in transgenic rice lines against Magnaporthe oryzae [17, 18]. The other four rice WAK genes, including OsWAK14, OsWAK91, OsWAK92 and OsWAK112d, were also suggested to be required for resistance to M. oryzae by loss-of-function mutants . Beyond rice and Arabidopsis, WAKs have been characterized in response to pathogens as well in other plants, such as tomato SlWAK1 (conferring resistance to Pseudomonas syringae) , maize ZmWAK (conferring resistance to Sporisorium reilianum)  and ZmWAK-RLK1 (conferring resistance to Setosphaeria turcica) .
An increasing number of findings indicate that WAK genes participate in plant−pathogen interactions. Therefore, in our study, we used the latest G. hirsutum genome sequence data (HAU version 1.1  to explore the WAK gene family, representing the first genome-wide identification of GhWAKs. Moreover, two GhWAKs were functionally characterized in response to V. dahliae infection using VIGS (virus induced gene silencing).
GhWAKs identification and localization
In total, 81 GhWAKs as candidates were identified and named according to their chromosomal locations. These GhWAKs were marked on the physical map of 18 chromosomes (Fig. 1a) and one scaffold664 (GhWAK65). A total of 34 and 46 GhWAKs were distributed in the A and D sub-genomes, respectively. Chromosome D02 harbored the largest number of GhWAKs with 20 genes. Six pairs of tandem duplication events were found, including GhWAK16/17, GhWAK36/37, GhWAK43/44–49, GhWAK50/52, GhWAK61/62 and GhWAK69/70/71. These results revealed that the evolution and expansion of GhWAKs happened in G. hirsutum, especially on chromosome D02. The detailed information about GhWAKs, including gene ID, open reading frame (ORF) length, amino acid length, protein molecular weight and isoelectric point, instability index and subcellular localization, was listed in Table 1.
All GhWAKs were predicated to be localized on plasma membrane (PM) (Table 1). In which, GhWAK77 as a representative was further located in tobacco epidermal cells using transient expression of fluorescent fusion proteins. The images clearly showed that fluorescent signal corresponding to the sole gfp (green fluorescent protein) gene was observed in PM, cytoplasm and nucleus. However, the fluorescent signal corresponding to GhWAK77-gfp was solely shown in PM (Fig. 1b). These suggested that GhWAKs might be a potential connector responsible for communication between inside and outside of the cell.
GhWAKs have conservative kinase domains and diverse extracellular domains
The majority of GhWAKs have 3–4 introns and show similar exon-intron structure (Figure S1). A total of six conserved protein domains were identified in GhWAKs, including GUB_WAK_bind (wall-associated receptor kinase galacturonan-binding, PF13947), WAK (wall-associated kinase, PF08488), WAK assoc. (wall-associated receptor kinase C-terminal, PF14380), EGF (EGF, PF00008; cEGF, PF12662; hEGF, PF12661; EGF_CA, PF07645; EGF_3, PF12947), DUF1199 (domain of unknown function, PF06712) and protein kinase domain (pkinase, PF00069; pkinase_Tyr, PF07714; kinase-like, PF14531; protein-kinase domain of FAM69, PF12260) (Fig. 2a). Cytoplasmic, extracellular and TM regions were predicated in the majority of GhWAKs, further indicating that they were PM proteins. Typical WAK encodes a transmembrane protein with a cytoplasmic kinase domain and an extracellular region. However, several proteins showed uncommon structural characteristics, such as the kinase domain in extracellular region, double TMs and kinase domains. All GhWAKs were classified into seven groups according to their protein domain analysis (Fig. 2b). The members in Group I, Group II and Group III were typical WAKs that contain EGF domain in extracellular region. The other four Groups, including IV, V, VI and VII, do not contain EGF. GhWAKs in Group I and IV contain both WAK and GUB domain. Inversely, GhWAKs in Group III neither contain WAK nor GUB domain. GhWAKs in Group II, VI and VII only contain GUB domain. However, II and VI are one-GUB-domain groups, and VII are two-GUB-domain group. GhWAKs in Group V only contain WAK domain. Additionally, DUF1199 domain was found in GhWAK31 and GhWAK77. Different types and numbers of extracellular domains were present in GhWAKs, indicating that they might sense or bind different outside signaling to trigger intracellular signaling pathways that control plant development and response to various environmental stresses.
Prediction of putative cis-regulatory elements in GhWAK promoters
The 2-kb region upstream of the translation start site of all GhWAKs were considered the promoter and analyzed for the potential roles of cis-regulatory elements (Fig. 3). These cis-regulatory elements were classified into four main groups including hormones, biotic, abiotic and light. Twelve hormone-responsive regulatory elements associated with abscisic acid (ABA) (ABRE, ABRE4 and AT-ABRE), auxin (IAA) (AuxRR-core, TGA-box and TGA-element), methyl jasmonate (MeJA) (CGTCA-motif), gibberellin (GA) (GARE-motif, P-box and TATC-box), SA (TCA-element) and ethylene (ET) (ERE), were identified. Of which, ABRE-motif, CGTCA-motif and ERE were enriched in the most of GhWAK promoters, indicating that they might be widely induced by ABA, JA and ET. The biotic stress-related regulatory elements, such as AT-rich, TC-rich repeats, W-box, WUN-motif, WRE3, JERE and box S, were involved in elicitor-mediated activation, wounding and pathogen responsiveness. In addition, eight abiotic-responsive regulatory elements, associated with anaerobic induction (ARE and GC-motif), low-temperature responsiveness (LTR), drought-inducibility (MBS, DRE core and DRE1), heat shock, osmotic stress, low pH, nutrient starvation (STRE) and stress-related (TCA), were identified in the GhWAK promoter regions. Moreover, various light-responsive elements were present in the promoters of GhWAKs. Especially, Box 4 and G-Box were widely harbored. These results indicated that GhWAKs might play vital roles in the response to various stresses, hormones and light.
GhWAKs were significantly induced by V. dahliae infection
To identify GhWAKs that were related to V. dahliae infection, two-fold changes were applied in transcript expression profiles from RNA-seq as minimum cutoffs. As a result, 26 GhWAKs were screened out, including 17 up-regulated and 9 down-regulated genes (Fig. 4a). Of which, 11 GhWAKs, including GhWAK5, GhWAK9, GhWAK77, GhWAK10, GhWAK45, GhWAK47, GhWAK78, GhWAK48, GhWAK31, GhWAK26 and GhWAK72, were significantly up-regulated in at least three time points, suggesting that they continuously responded to V. dahliae infection. Their expression profiles were further verified through real-time quantitative reverse transcription PCR (qRT-PCR). The expression results of them in response to the V. dahliae infection from qRT-PCR were consistent with those found in RNA-seq data (Fig. 4b). Due to the high degree of sequence similarity in GhWAKs family, it was difficult to design specific primers for four gene pairs, including GhWAK4/GhWAK45, GhWAK5/GhWAK49, GhWAK10/GhWAK55, and GhWAK31/GhWAK77. The results of qRT-PCR indicated that these four pairs of GhWAKs were dramatically up-regulated. According to RNA-seq data, GhWAK4, GhWAK49 and GhWAK55 did not show to be up-regulated. Thus, the expression changes found using qRT-PCR probably more represent the responses of GhWAK45, GhWAK5 and GhWAK10 to V. dahliae infection. In addition, the other 45 GhWAKs that did not show differential expression in RNA-seq data were further detected through qRT-PCR. As a result, GhWAK1 and GhWAK69 showing higher transcription levels in cotton seedlings inoculated with V. dahliae than that in control was screened out complementally (Fig. 4c). Finally, a total of 28 GhWAKs were found to play a potential role in the interaction between cotton and V. dahliae.
Silencing GhWAKs compromised cotton resistance to Verticillium wilt
GhWAK26 and GhWAK77 showed obviously and persistently up-regulated expression to the infection from V. dahliae (Fig. 4). In addition, they contain cis-elements in their promoters associated with MeJA and SA, which play key roles in cotton resistance to V. dahliae. Thus, to further reveal the function of GhWAKs in cotton resistance to V. dahliae, GhWAK26 and GhWAK77 were prioritized for study as representatives here using tobacco rattle virus (TRV) based VIGS system. At approximately two weeks post-infiltration with a mixture of Agrobacterium cultures containing pTRV1 and pTRV2-CLA1, a strong photobleaching phenotype was shown on the newly emerging true leaves (Fig. 5a), indicating that VIGS system worked well. Then, the expression of GhWAK26 and GhWAK77 was detected in the leaves infiltrated with pTRV2-GhWAK26 and pTRV2-GhWAK77, respectively. As shown in Fig. 5b, the expression of GhWAK26 and GhWAK77 was reduced by about 80%, suggesting VIGS triggered their silencing in cotton plants. At 20 days post inoculation (dpi), VIGSed plants (Fig. 5d and e) exhibited more chlorosis and wilting than the control plants infiltrated with Agrobacterium cultures containing empty vector pTRV1 and pTRV2 (Fig. 5c). The disease indices of VIGSed plants were also significantly higher than those of the control at 15 dpi and 20 dpi (Fig. 5f). Therefore, the results of VIGS assays suggested that GhWAK26 and GhWAK77 were important participants in cotton resistance to V. dahliae infection.
Silencing GhWAKs increased the spread of V. dahliae in cotton stems
After inoculation, V. dahliae in cotton stem was detected by PCR. No specific amplification products from V. dahliae were shown in CK at 5 dpi and 7 dpi, indicating that V. dahliae had not yet invaded the stems or multiplied in large quantities (Fig. 6a and Figure S2). However, at 5 dpi, few specific products from V. dahliae were amplified in GhWAK26-silenced and GhWAK77-silenced plants, representing a small amount of pathogen invasion. Further, at 7 dpi, the bright bands amplified from GhWAK26-silenced and GhWAK77-silenced plant stems appeared on agarose gels, indicating that V. dahliae had invaded largely. In addition, pathogen isolation on potato dextrose agar (PDA) showed that a large number of V. dahliae grew out from the stems of GhWAK26-silenced and GhWAK77-silenced cotton plants, while no mycelium was shown from the control (Fig. 6b). Both PCR detection and PDA culture results suggested that silencing GhWAKs significantly increased the spread of V. dahliae in the cotton stems.
Lignin is considered to play an important role in preventing cotton from the infection of V. dahliae. Therefore, we further compared the changes of lignin content in GhWAK-silenced cotton stems with CK. The results showed that the lignin content in GhWAK-silenced plants was significantly lower than that in CK (Fig. 6c), which might affect the stem structure and then reduce the prevention of cotton from V. dahliae infection.
Silencing GhWAKs dramatically affected V. dahliae-induced H2O2, nitric oxide (NO) and peroxidase (POD)
The content of H2O2 and NO, and POD activity in GhWAK-silenced plants inoculated with V. dahliae were further measured. GhWAKs silencing caused lower levels of H2O2 at 6 h post inoculation (hpi), 12 hpi and 24 hpi (Fig. 7a and b). Both GhWAK26- and GhWAK77-silenced plants accumulated greatly depressed levels of NO comparing with CK (Fig. 7c and d). However, the activity of POD significantly elevated in GhWAK26- and GhWAK77-silenced plants at 6 hpi, 24 hpi and 48 hpi, except at 12 hpi (Fig. 7e and f).
Silencing GhWAKs significantly affected the expression of JA and SA-related marker genes
Further, the expression of several JA and SA-related marker genes involved in plant defense signaling pathways was detected. The expression of JAZ1 (jasmonate-zim-domain protein), JAZ3, JAZ6, LOX1 (lipoxygenase) (JA-related marker genes), PR3 (pathogenesis related protein) and NPR1 (nonexpresser of PR protein) (SA-related marker genes) were significantly down-regulated after silencing GhWAK26 in cotton (Fig. 8a). In GhWAK77-silenced plants, JAZ6 and three important genes involved in the SA signaling pathway, including ICS1 (isochorismate synthase), NPR1 and EDS1 (enhanced disease susceptibility), were down-regulated comparing with control. On the contrary, the expression of JAZ1 and LOX1 were significantly up-regulated due to the silencing of GhWAK77 (Fig. 8b). These results indicated that GhWAK26 and GhWAK77 might involve in cotton resistance to V. dahliae through SA and JA signaling pathways.
WAK gene family has been analyzed in several plant species, such as A. thaliana , O. sativa , P. trichocarpa  and B. rapa . Some WAKs have been implicated in the response to pathogen infection. Examples are Arabidopsis Wak1 , maize ZmWAK-RLK1 (Htn1) and ZmWAK (qHSR1) [21, 22], wheat Stb6 and TaWAK6 [24, 25], rice Xa4, OsWAK1 and OsWAK91 [16, 17, 26], and orange CsWAKL08 , conferring host plant disease resistance. In the present work, a total of 81 GhWAKs were systematically identified and analyzed for the first time from a high-quality G. hirsutum genome (Table 1) . Of which, 28 GhWAKs were potentially involved into the interaction between cotton and V. dahliae (Fig. 4). Especially, silencing of GhWAK26 or GhWAK77 dramatically reduced the resistance of cotton plants to V. dahliae infection (Fig. 5), suggesting that WAKs were important resistance genes during cotton–pathogen interactions.
At the PM, RLKs as cell-surface receptors can perceive and process extracellular danger signals to trigger plant defense responses . WAK belongs to RLK subfamily. All GhWAKs contain a typical eukaryotic kinase domain that is mostly present in intracellular region and relatively well conserved (Fig. 2a). In addition, GhWAKs locate on PMs in all probability (Table 1, Fig. 1), suggesting that GhWAKs have potential roles in communicating between inside and outside of the cell. In order to penetrate plant roots to gain access to the xylem and to spread in the vascular system, V. dahliae usually secretes various toxins and carbohydrate active enzymes, including glycoproteins and cell wall-degrading enzymes [29, 30]. Therefore, it is conceivable that V. dahliae infection affects plant cell wall integrity (CWI) and generates some degradation products, which are important defense signals . In the extracellular region, GhWAKs contain five different domains (Fig. 2a), which may sense CWI or interact with different components of these extracellular matrix, such as glycine-rich protein, pectin and oligogalacturonides (OGs) [32,33,34].
At present, the molecular mechanism of WAK-mediated resistance remains largely unknown. However, some defence responses associated with WAKs have been reported, including cell wall reinforcement , pathogenesis-related genes activation , SA or JA accumulation , POD and superoxide dismutase activities , and reactive oxygen species (ROS) homeostasis . Here, silencing GhWAKs resulted in the up- or down-regulation of several genes (Fig. 8) and depressed cotton resistance to V. dahliae. Among them, JAZ and LOX are associated with JA-mediated defense responses . NPR1, ICS1 and EDS1 are associated with SA-mediated defense responses . The two phytohormones, JA and SA, have been known to be involved into the regulation of plant resistance against V. dahliae [37, 38]. In addition, some hormone-responsive and biotic stress-related regulatory elements were enriched in the promoters of GhWAKs (Fig. 3). Thus, these findings suggest that GhWAK function as a mediator to active intracellular SA and JA signaling pathways to regulate cotton resistance.
V. dahliae is a vascular pathogen that penetrates the host roots and then extends to other overground parts of plant through the process of transpiration [29, 37]. The improvement of physical, chemical and structural barriers, such as ROS, NO, cell wall, lignin, callose and POD, contributes to preventing expansion and reducing colonization of V. dahliae in cotton tissues [37, 39,40,41]. In this study, more V. dahliae was detected in GhWAK26-silenced or GhWAK77-silenced plants with lower lignin contents than in CK (Fig. 6). Moreover, silencing of GhWAKs in cotton plants dramatically compromised V. dahliae-induced accumulation of H2O2 and NO, but enhanced POD activity (Fig. 7). These findings demonstrate that GhWAKs play roles in preventing pathogen spreading at least in part by regulating the accumulation of lignin, H2O2 and NO, and the activity of POD. Overall, these results augment our knowledge about cotton WAK gene family, and particularly promote the understanding on their function in disease resistance.
In this study, we carried out a genome-wide analysis of WAK gene family in G. hirsutum with the identification of 81 putative GhWAKs, which might sense different outside signals to trigger intracellular signaling pathways that response to various environment-stresses. Of which, 28 GhWAKs with potential roles in the interaction between cotton and V. dahliae were screened out. Silencing GhWAKs could significantly affect the expression of JA- and SA-related marker genes, increased the spread of V. dahliae in the cotton stems, dramatically compromised V. dahliae-induced accumulation of lignin, H2O2 and NO, but enhanced POD activity. These results provided direct evidences that GhWAKs participate in the cotton resistance to V. dahliae. Finally, a model for how GhWAKs were involved in cotton resistance to V. dahliae was proposed (Fig. 9).
Identification and bioinformatics analysis of GhWAKs
The amino acid and nucleotide sequences of WAKs from Arabidopsis accessed from TAIR website (https://www.Arabidopsis.org/) were queried against G. hirsutum genome database (HAU) in CottonFGD (https://cottonfgd.org/) using BLAST program (E-value < 0.01) [7, 23]. The obtained putative GhWAKs were further identified by HMMER software (HMM Database = Pfam; Significance E-values < 0.01) (https://www.ebi.ac.uk/Tools/hmmer/search/hmmscan) to confirm the presence of conserved protein domains.
Functional sites and transmembrane topology for all putative GhWAKs were analyzed through PROSITE database (https://prosite.expasy.org/) and Phobius database (http://phobius.sbc.su.se/), respectively. The number of amino acids, molecular weight, theoretical isoelectric point and instability index of proteins were analyzed using ExPASy program (http://www.expasy.org/). Prediction of protein subcellular localization was performed using CELLO v2.5 (http://cello.life.nctu.edu.tw/) and ProtComp 9.0 (http://www.softberry.com/berry.phtml?topic=protcomppl&group=programs&subgroup=proloc). Signal peptides were predicted using SignalP 5.0 (http://www.cbs.dtu.dk/services/SignalP/).
Analysis of chromosomal location, genes structure and cis-elements
The information about physical chromosomal locations and gene structures of GhWAKs was extracted from the gene annotations in gene feature format (GFF) files, which were downloaded from the CottonFGD website and analyzed by TBtools software . The potential promoter sequences, 2 kb upstream of GhWAKs, were also extracted from G. hirsutum genome database. The cis-elements in the potential promoters were predicted using PlantCARE databases (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/).
Plant materials and V. dahliae inoculation
The seeds of Nicotiana benthamiana and G. hirsutum cv. Nongda 601 (ND601) were preserved at the State Key Laboratory of North China Crop Improvement and Regulation, Hebei Agricultural University, China. N. benthamiana was grown in the greenhouse about 5 weeks at 21 °C with 14/10 h (light/dark) photoperiod. ND601 were grown in the greenhouse at 25 °C under a 14-h light/10-h dark cycle with relative humidity about 70%. Cotton seedlings inoculation with V. dahliae strain Linxi 2–1 (107 spores ml− 1) was performed as previously described .
Proteins subcellular localization
The ORF of GhWAK77 (without the stop codon) was amplified by PCR with primers gWAK77-F and gWAK77-R (Table S1), and then introduced into entry vector pDONR™207 by attB/attP recombination reaction, as described by the manufacturer (Invitrogen). The GhWAK77 fragment was transferred from the entry clone to expression vector pEarlyGate103  with attL/attR recombinant reaction, as described by the manufacturer (Invitrogen). The recombinant expression vector was introduced into Agrobacterium tumefaciens GV3101, cultured and infiltrated into four-week-old tobacco leaves via the method described by . After 2 days, GFP signal in the tobacco leaf epidermal cells was examined using a laser scanning microscope (FluoView FV1000; Olympus).
RNA-seq data and qRT-PCR analysis
The transcription patterns of GhWAKs in cotton roots after inoculation with V. dahliae were analyzed using high-through RNA-seq data published previously . Log2Fold change were calculated from FPKM (fragments per kilobase of exon model per million mapped) and used for the heat map of hierarchical clustering with the TBtools v0.67 software . Total RNA was extracted using EASYspin Plant RNA kit (Aidlab, Beijing, China) according to the manufacturer’s instructions. The quality and concentration of RNA were detected by 1.5% agarose gel electrophoresis and NanoDrop™ 1000 spectrophotometer (Thermo Fisher Scientific), respectively. cDNA was synthesized with a reverse transcription kit (ReverTra Ace® qPCR RT Master Mix with gDNA Remover, TaKaRa, Dalian, China). qRT-PCR was performed using 7500 Real Time PCR System (Applied Biosystems, USA) with THUNDERBIRD®SYBR® qPCR Mix (TaKaRa, Dalian, China). The 2-ΔΔCt method was used to calculate the relative expression of genes. GhHis 3 was used as internal reference. Three biological repeats were taken for each treatment.
VIGS assays in cotton
The vectors for VIGS, pTRV1 and pTRV2, were kindly provided by Professor Liu Yule of Tsinghua University . The fragments from GhWAKs were amplificated by PCR and inserted into the pTRV2 vector between EcoR I and Kpn I. The constructed vectors were separately transferred into A. tumefaciens strain GV3101 by freeze-thaw method . VIGS in cotton was performed as described previously . At least 30 plants were used per treatment, and each treatment was repeated three times. Plant resistance to V. dahliae was assayed by analyzing disease index .
Detection and isolation of V. dahliae in cotton stems
At 5 dpi and 7 dpi, 1 cm and 0.5 cm of samples excised at a height of 0.5 cm stem above ground were used for detection and isolation of V. dahliae, respectively. V. dahliae detection by PCR was performed using primers P1 and P2 . V. dahliae isolation from cotton stems was carried out according to the previous method . Twenty-four individual plants were sampled for each treatment and repeated three times.
Measurements of NO, H2O2 and POD activity
The first true leaves of cotton seedlings were powdered in the mortar with liquid nitrogen and homogenized using 50 mM sodium phosphate buffer (pH 7.0). After centrifugation (14,000 g, 20 min), the supernatants were used for the determination of NO, H2O2 and POD activity with commercialized assay kits (Nanjing Jiancheng Bioengineering Institute, China), following the manuals. The total protein concentration of the supernatants was measured using Pierce™ BCA Protein Assay Kit (Thermo Scientific).
Primers and statistical analysis
All primers used in this study were listed in Table S1. Differences between measured values were analyzed using software GraphPad Prism® 8 (GraphPad, San Diego, CA, USA). A two-way ANOVA with multiple comparisons (Sidak’s test) was used to compare gene expression in cotton roots between inoculated with V. dahliae and inoculated with water (CK) at the same hpi, disease indices, H2O2 and NO content, and POD activity between GhWAK-silenced plants and CK. A one-way ANOVA with Dennett’s multiple-comparisons test was used to compare lignin content between GhWAK-silenced plants and CK. The P-value less than 0.05 was assumed to be statistically significant.
Availability of data and materials
All GhWAKs sequence information is available in the Cotton Functional Genomics Database (CottonFGD) (https://cottonfgd.org/about/download.html). The data generated or analyzed during the current study are included in this published article and its supplemental data files and available from the corresponding author on reasonable request.
Cell wall integrity
- EDS :
Enhanced disease susceptibility
Epidermal growth factor
Fragments per kilobase of exon model per million mapped
Hours post inoculation
- ICS :
- JAZ :
- LOX :
- NPR :
Nonexpresser of PR protein
Pathogenesis related protein
Reactive oxygen species
Virus-induced gene silencing
Cai YF, He XH, Mo JC, Sun Q, Yang JP, Liu JG. Molecular research and genetic engineering of resistance to Verticillium wilt in cotton: a review. Afr J Biotechnol. 2009;8(25):7363–72.
Verhage A, van Wees SC, Pieterse CM. Plant immunity: it's the hormones talking, but what do they say? Plant Physiol. 2010;154(2):536–40. https://doi.org/10.1104/pp.110.161570.
Bacete L, Mélida H, Miedes E, Molina A. Plant cell wall-mediated immunity: cell wall changes trigger disease resistance responses. Plant J. 2018;93(4):614–36. https://doi.org/10.1111/tpj.13807.
Novaković L, Guo TT, Bacic A, Sampathkumar A, Johnson K. Hitting the wall-sensing and signaling pathways involved in plant cell wall remodeling in response to abiotic stress. Plants. 2018;7(4):89. https://doi.org/10.3390/plants7040089.
Shiu SH, Bleecker AB. Receptor-like kinases from Arabidopsis form a monophyletic gene family related to animal receptor kinases. Proc Natl Acad Sci U S A. 2001;98(19):10763–8. https://doi.org/10.1073/pnas.181141598.
Si K, Tanizawa K. Involvement of epidermal growth factor-like domain of NELL proteins in the novel protein–protein interaction with protein kinase C. Biochem Biophys Res Commun. 1999;265(3):752–7.
Verica JA, He ZH. The cell wall-associated kinase (WAK) and WAK-like kinase gene family. Plant Physiol. 2002;129(2):455–9. https://doi.org/10.1104/pp.011028.
de Oliveira LFV, Christoff AP, de Lima JC, de Ross BCF, Sachetto-Martins G, Margis-Pinheiro M, et al. The wall-associated kinase gene family in rice genomes. Plant Sci. 2014;229:181–92. https://doi.org/10.1016/j.plantsci.2014.09.007.
Zhang B, Li P, Su TB, Li PR, Xin XY, Wang WH, et al. Comprehensive analysis of wall-associated kinase genes and their expression under abiotic and biotic stress in Chinese cabbage (Brassica rapa ssp. pekinensis). J Plant Growth Regul. 2020;39:72–86. https://doi.org/10.1007/s00344-019-09964-3.
Tocquard K, Lafon-Placette C, Auguin D, Muries B, Bronner G, Lopez D, et al. In silico study of wall-associated kinase family reveals large-scale genomic expansion potentially connected with functional diversification in Populus. Tree Genet Genom. 2014;10(5):1135–47. https://doi.org/10.1007/s11295-014-0748-7.
He ZH, Fujiki M, Kohorn BD. A cell wall-associated, receptor-like protein kinase. J Biol Chem. 1996;271(33):19789–93. https://doi.org/10.1074/jbc.271.33.19789.
He ZH, He D, Kohorn BD. Requirement for the induced expression of a cell wall associated receptor kinase for survival during the pathogen response. Plant J. 1998;14(1):55–63. https://doi.org/10.1046/j.1365-313X.1998.00092.x.
Wang N, Huang HJ, Ren ST, Li JJ, Sun Y, Sun DY, et al. The rice wall-associated receptor-like kinase gene OsDEES1 plays a role in female gametophyte development. Plant Physiol. 2012;160(2):696–707. https://doi.org/10.1104/pp.112.203943.
Kanneganti V, Gupta AK. RNAi mediated silencing of a wall associated kinase, OsiWAK1 in Oryza sativa results in impaired root development and sterility due to anther indehiscence. Physiol Mol Biol Plants. 2011;17(1):65–77. https://doi.org/10.1007/s12298-011-0050-1.
Kaur R, Singh K, Singh J. A root-specific wall-associated kinase gene, HvWAK1, regulates root growth and is highly divergent in barley and other cereals. Funct Integr Genomic. 2013;13(2):167–77. https://doi.org/10.1007/s10142-013-0310-y.
Hu K, Cao JB, Zhang J, Xia F, Ke YG, Zhang HT, et al. Improvement of multiple agronomic traits by a disease resistance gene via cell wall reinforcement. Nat Plants. 2017;3:17009. https://doi.org/10.1038/nplants.2017.9.
Li H, Zhou SY, Zhao WS, Su SC, Peng YL. A novel wall-associated receptor-like protein kinase gene, OsWAK1, plays important roles in rice blast disease resistance. Plant Mol Biol. 2009;69(3):337–46. https://doi.org/10.1007/s11103-008-9430-5.
Harkenrider M, Sharma R, De Vleesschauwer D, Tsao L, Zhang XT, Chern M, et al. Overexpression of rice wall-associated kinase 25 (OsWAK25) alters resistance to bacterial and fungal pathogens. PLoS One. 2016;11(1):e0147310. https://doi.org/10.1371/journal.pone.0147310.
Delteil A, Gobbato E, Cayrol B, Estevan J, Michel-Romiti C, Dievart A, et al. Several wall-associated kinases participate positively and negatively in basal defense against rice blast fungus. BMC Plant Biol. 2016;16(1):17. https://doi.org/10.1186/s12870-016-0711-x.
Rosli HG, Zheng Y, Pombo MA, Zhong SL, Bombarely A, Fei ZJ, et al. Transcriptomics-based screen for genes induced by flagellin and repressed by pathogen effectors identifies a cell wall-associated kinase involved in plant immunity. Genome Biol. 2013;14(12):R139. https://doi.org/10.1186/gb-2013-14-12-r139.
Zuo WL, Chao Q, Zhang N, Ye JR, Tan GQ, Li BL, et al. A maize wall-associated kinase confers quantitative resistance to head smut. Nat Genet. 2015;47(2):151–7. https://doi.org/10.1038/ng.3170.
Hurni S, Scheuermann D, Krattinger SG, Kessel B, Wicker T, Herren G, et al. The maize disease resistance gene Htn1 against northern corn leaf blight encodes a wall-associated receptor-like kinase. Proc Natl Acad Sci U S A. 2015;112(28):8780–5. https://doi.org/10.1073/pnas.1502522112.
Wang M, Tu L, Yuan D, Zhu D, Shen C, Li JY, et al. Reference genome sequences of two cultivated allotetraploid cottons, Gossypium hirsutum and Gossypium barbadense. Nat Genet. 2019;51(2):224–9. https://doi.org/10.1038/s41588-018-0282-x.
Saintenac C, Lee WS, Cambon F, Rudd JJ, King RC, Marande W, et al. Wheat receptor-kinase-like protein Stb6 controls gene-for-gene resistance to fungal pathogen Zymoseptoria tritici. Nat Genet. 2018;50(3):368–74. https://doi.org/10.1038/s41588-018-0051-x.
Marta DB, Yuliya K, Andrzej Z, Przemysław W, Anna NO, Wojciech MK, et al. TaWAK6 encoding wall-associated kinase is involved in wheat resistance to leaf rust similar to adult plant resistance. PLOS One. 2020;15(1):e0227713. https://doi.org/10.1371/journal.pone.0227713.
Al-Bader N, Meier A, Geniza M, Gongora YS, Oard J, Jaiswal P. Loss of premature stop codon in the Wall-Associated Kinase 91 (OsWAK91) gene confers sheath blight disease resistance in rice. bioRxiv. 2019;625509. https://doi.org/10.1101/625509.
Li Q, Hu AH, Qi JJ, Dou WF, Qin XJ, Zou XP, et al. CsWAKL08, a pathogen-induced wall-associated receptor-like kinase in sweet orange, confers resistance to citrus bacterial canker via ROS control and JA signaling. Hortic Res. 2020;7(1):42. https://doi.org/10.1038/s41438-020-0263-y.
van der Burgh AM, Joosten MHAJ. Plant immunity: thinking outside and inside the box. Trends Plant Sci. 2019;24(7):587–601. https://doi.org/10.1016/j.tplants.2019.04.009.
Fradin EF, Thomma BP. Physiology and molecular aspects of Verticillium wilt diseases caused by V. dahliae and V. albo-atrum. Mol Plant Pathol. 2006;7(2):71–86. https://doi.org/10.1111/j.1364-3703.2006.00323.x.
Klosterman SJ, Subbarao KV, Kang S, Veronese P, Gold SE, Thomma BP, et al. Comparative genomics yields insights into niche adaptation of plant vascular wilt pathogens. PLoS Pathog. 2011;7(7):e1002137. https://doi.org/10.1371/journal.ppat.1002137.
Rui Y, Dinneny JR. A wall with integrity: surveillance and maintenance of the plant cell wall under stress. New Phytol. 2020;225(4):1428–39. https://doi.org/10.1111/nph.16166.
Giarola V, Krey S, von den Driesch B, Bartels D. The Craterostigma plantagineum glycine-rich protein CpGRP1 interacts with a cell wall-associated protein kinase 1 (CpWAK1) and accumulates in leaf cell walls during dehydration. New Phytol. 2016;210(2):535–50. https://doi.org/10.1111/nph.13766.
Kohorn BD, Kohorn SL. The cell wall-associated kinases, WAKs, as pectin receptors. Front Plant Sci. 2012;3:88. https://doi.org/10.3389/fpls.2012.00088.
Brutus A, Sicilia F, Macone A, Cervone F, De Lorenzo G. A domain swap approach reveals a role of the plant wall-associated kinase 1 (WAK1) as a receptor of oligogalacturonides. Proc Natl Acad Sci U S A. 2010;107(20):9452–7. https://doi.org/10.1073/pnas.1000675107.
Wasternack C, Hause B. Jasmonates: biosynthesis, perception, signal transduction and action in plant stress response, growth and development. An update to the 2007 review in Annals of Botany. Ann Bot. 2013;111(6):1021–58. https://doi.org/10.1093/aob/mct067.
Zheng XY, Zhou M, Yoo H, Prunedapaz JL, Spivey NW, Kay SA, et al. Spatial and temporal regulation of biosynthesis of the plant immune signal salicylic acid. Proc Natl Acad Sci U S A. 2015;112(30):9166–73. https://doi.org/10.1073/pnas.1511182112.
Zhang Y, Wang X, Rong W, Yang J, Li ZK, Wu LQ, et al. Histochemical analyses reveal that stronger intrinsic defenses in Gossypium barbadense than in G. hirsutum are associated with resistance to Verticillium dahliae. Mol Plant-Microbe Interact. 2017;30(12):984–96. https://doi.org/10.1094/MPMI-03-17-0067-R.
Gao W, Long L, Zhu L, Xu L, Gao WH, Sun LQ, et al. Proteomic and virus-induced gene silencing (VIGS) analyses reveal that gossypol, brassinosteroids, and jasmonic acid contribute to the resistance of cotton to Verticillium dahliae. Mol Cell Proteomics. 2013;12(12):3690–703. https://doi.org/10.1074/mcp.M113.031013.
Yang J, Zhang Y, Wang XF, Wang WQ, Li ZK, Wu JH, et al. HyPRP1 performs a role in negatively regulating cotton resistance to V. dahliae via the thickening of cell walls and ROS accumulation. BMC Plant Biol. 2018;18(1):339. https://doi.org/10.1186/s12870-018-1565-1.
Xu L, Zhu L, Tu L, Liu LL, Yuan DJ, Jin L, et al. Lignin metabolism has a central role in the resistance of cotton to the wilt fungus Verticillium dahliae as revealed by RNA-Seq-dependent transcriptional analysis and histochemistry. J Exp Bot. 2011;62(15):5607–21. https://doi.org/10.1093/jxb/err245.
Zhang J, Hu HL, Wang XN, Yang YH, Zhang CJ, Zhu HQ, et al. Dynamic infection of Verticillium dahliae in upland cotton. Plant Biol. 2020;22(1):90–105. https://doi.org/10.1111/plb.13037.
Chen CJ, Chen H, Zhang Y, Thomas HR, Frank MH, He YH, et al. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020;13(8):1194–202. https://doi.org/10.1016/j.molp.2020.06.009.
Earley KW, Haag JR, Pontes O, Opper K, Juehne T, Song K, et al. Gateway-compatible vectors for plant functional genomics and proteomics. Plant J. 2006;45(4):616–29. https://doi.org/10.1111/j.1365-313X.2005.02617.x.
Sparkes IA, Runions J, Kearns A, Hawes C. Rapid, transient expression of fluorescent fusion proteins in tobacco plants and generation of stably transformed plants. Nat Protoc. 2006;1(4):2019–25. https://doi.org/10.1038/nprot.2006.286.
Liu Y, Schiff M, Marathe R, Dinesh-Kumar SP. Tobacco Rar1, EDS1 and NPR1/NIM1 like genes are required for N-mediated resistance to tobacco mosaic virus. Plant J. 2002;30(4):415–29. https://doi.org/10.1046/j.1365-313X.2002.01297.x.
Chen H, Nelson RS, Sherwood JL. Enhanced recovery of transformants of Agrobacterium tumefaciens after freeze-thaw transformation and drug selection. BioTechniques. 1994;16(4):664–8 670.
Gao XQ, Britt RC Jr, Shan LB, He P. Agrobacterium-mediated virus-induced gene silencing assay in cotton. J Vis Exp. 2011;54:e2938. https://doi.org/10.3791/2938.
Zhang BL, Yang YW, Chen TZ, Yu WG, Liu TL, Li HJ, et al. Island cotton Gbve1 gene encoding a receptor-like protein confers resistance to both defoliating and non-defoliating isolates of Verticillium dahliae. PLoS One. 2012;7(12):e51091. https://doi.org/10.1371/journal.pone.0051091.
Zhu YY, Wang YY, Lyon BR. PCR detection of Verticillium dahliae from diseased plants. Acta Phytopathol Sinica. 1999;29(3):250–5.
Fradin EF, Zhang Z, Juarez Ayala JC, Castroverde CD, Nazar RN, Robb J, et al. Genetic dissection of Verticillium wilt resistance mediated by tomato Ve1. Plant Physiol. 2009;150(1):320–32. https://doi.org/10.1104/pp.109.136762.
All authors are grateful to the laboratory members for help, advice and discussion.
This work was supported by grants from the China Agricultural Research System (CARS-15-03), the Science and Technology Support Program of Hebei Province (16226307D) and the Outstanding Youth Found of Hebei Province (C2019204365). The funding bodies provided the financial support to this research, including experimental design and implementation, sampling and data analysis. No funder played the role in data collection and analysis and writing the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Primers list.
. Gene structures of GhWAKs.
. Detection of V. dahliae in cotton stems by PCR. M, marker. DNA templates were from V. dahliae spores as positive control (V), water as negative control (W), and cotton seedling stems (lane 1–4 for 5 dpi and lane 5–8 for 7 dpi).
About this article
Cite this article
Yang, J., Xie, M., Wang, X. et al. Identification of cell wall-associated kinases as important regulators involved in Gossypium hirsutum resistance to Verticillium dahliae. BMC Plant Biol 21, 220 (2021). https://doi.org/10.1186/s12870-021-02992-w