Search tips
Search criteria 


Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
Sci Rep. 2016; 6: 29022.
Published online 2016 June 29. doi:  10.1038/srep29022
PMCID: PMC4926273

Discovery and identification of candidate genes from the chitinase gene family for Verticillium dahliae resistance in cotton


Verticillium dahliae, a destructive and soil-borne fungal pathogen, causes massive losses in cotton yields. However, the resistance mechanism to V. dahilae in cotton is still poorly understood. Accumulating evidence indicates that chitinases are crucial hydrolytic enzymes, which attack fungal pathogens by catalyzing the fungal cell wall degradation. As a large gene family, to date, the chitinase genes (Chis) have not been systematically analyzed and effectively utilized in cotton. Here, we identified 47, 49, 92, and 116 Chis from four sequenced cotton species, diploid Gossypium raimondii (D5), G. arboreum (A2), tetraploid G. hirsutum acc. TM-1 (AD1), and G. barbadense acc. 3–79 (AD2), respectively. The orthologous genes were not one-to-one correspondence in the diploid and tetraploid cotton species, implying changes in the number of Chis in different cotton species during the evolution of Gossypium. Phylogenetic classification indicated that these Chis could be classified into six groups, with distinguishable structural characteristics. The expression patterns of Chis indicated their various expressions in different organs and tissues, and in the V. dahliae response. Silencing of Chi23, Chi32, or Chi47 in cotton significantly impaired the resistance to V. dahliae, suggesting these genes might act as positive regulators in disease resistance to V. dahliae.

Chitinases (EC, which are found in a variety of organisms, catalyze the hydrolysis of the β-1-4-linkage in the N-acetyl-D-glucosamine polymer of chitin, a major component of fungal cell walls and arthropod exoskeletons1. Hence, they play a major role in defense by directly attacking invading fungal pathogens. Analysis of the amino acid sequences of these enzymes reveals that chitinases have multiple domains; for example, a typical plant chitinase possesses a signal region at the N-terminal and a catalytic domain, and only vacuolar chitinases have C-terminal regions. In addition, some chitinases also have one or more chitin-binding domains (CBDs) behind the N-terminal signal region, and a variable linking region following the chitin-binding domains2,3,4,5.

Of the glycoside hydrolase (GH) gene families, families 18 and 19 are considered as chitinases due to the sequence similarity of their catalytic domains2. Based on the amino acid sequence homology, three-dimensional (3D) structures, substrate specificity, mechanisms of catalytic reactions, sensitivity to inhibitors and other characteristics, chitinases have been phylogenetically categorized into five classes (class I-V)3,4,5,6,7,8. Of them, GH family 18 enzymes in classes III and V are more widespread and distributed in a greater diversity of organisms such as fungi, bacteria, viruses, animals, and higher plants than other enzymes. However, members of GH family 19, which contains class I, II, and IV chitinases, are found mainly in higher plants and some bacteria8. Structurally, class I chitinases have an N-terminal cysteine-rich CBD and a catalytic domain (CatD) at the C-terminal. These two domains are linked by rich proline and glycine, and the length and composition of the linkers are variable9. Class II chitinases have high sequence homology with class I chitinases at the CatD but they lack the CBDs and linker regions. Class III chitinases possess lysozyme activity, and do not have any sequence homology with either class I or II chitinases. However, they do have a high sequence similarity to chitinases in the same class in different plants10. Class IV chitinases, which have CBDs and CatDs, are similar to the class I chitinases, but both domains are smaller than those in class I due to one deletion in the CBD and three deletions in the CatD. Class V chitinases possess two CBDs in tandem11. The multiple domains and the complex structures in chitinases are indicative of the variety of functions of these enzymes.

On the basis of their physicochemical properties and enzymatic activity, plant chitinases are considered as pathogenesis-related (PR) proteins and have been categorized as PR-3, PR-4, PR-8, and PR-114. In practice, plant chitinases are more propitious to inhibit the growth of fungi by producing hypersensitive reactions and inducing defense responses12. These enzymes also generate or degrade some signaling molecules that function by regulating the cell polarity, cell movement and cell division in plant growth and development13,14. Moreover, the expression patterns or production of the chitinases can be regulated by various abiotic stresses such as wounds, cold, drought, heavy metals, ozone, excessive salinity, and UV light15,16,17,18, and a variety of plant phytohormones including jasmonic acid (JA), salicylic acid (SA), ethylene (ET), cytokinin (CK), and indole acetic acid (IAA)19.

Cotton is the most important source of natural textile fibers and oilseed in the world, but the production of cotton is constrained by various biotic and abiotic stresses. Among the biotic stresses, fungi are the most severe pathogenic microorganisms. They deliver effector molecules into the plant cell to promote virulence and cause disease20. V. dahliae, a soil-borne fungal pathogen, is considered as the most serious factor in cotton yield and quality loss21,22. The most effective method of developing pathogen-resistant varieties is to mine key genes. Members of the chitinase family are pathogenesis-related proteins, and have been identified in almost all plants. Many of these genes have been proven to inhibit fungal growth in vitro23,24, and the overexpression of chitinase genes from some plants or the fungus Trichoderma species had shown to enhance resistance to fungal pathogens in a variety of transgenic plants25,26,27,28. However, to date, the chitinase genes in cotton have not been systematically analyzed. Although several chitinases from other organisms have been transformed into cotton, the resistance of these transgenic lines to pathogens is not obvious29.

The availability of data on the whole-genome of Gossypium in different cotton species, such as G. raimondii (D5)30, G. arboreum (A2)31, G. hirsutum acc. TM-1 (AD1)32, and G. barbadense acc. 3–79 (AD2)33, made it possible to identify chitinase genes on a genome-wide scale, and to enrich our understanding of the molecular mechanism of resistance to V. dahliae in cotton. Here, we surveyed chitinase members in four sequenced cotton species, and analyzed their phylogenetic relationships, gene structure, and expression patterns in different tissues and in response to pathogens and hormone treatments. Furthermore, we verified the functional roles of three chitinase genes that are significantly induced in cotton in response to V. dahliae by virus-induced gene silencing (VIGS) analysis. Our studies might provide effective gene resources for improving cotton V. dahliae resistance in future cotton-breeding programs.


Genome-wide identification of chitinase genes and their chromosomal distribution

The whole genome sequence scaffolds of four sequenced cotton species, G. raimondii30, G. arboreum31, G. hirsutum acc. TM-132, and G. barbadense acc. 3–7933 were used for the genome-wide exploration of chitinase genes in Gossypium. We used the HMMER 3.0 and the Pfam protein family databases with the Glyco_hydro_18 (PF00704) and Glyco_hydro_19 (PF00182) domains to search for chitinases in the four cotton genomes. SMART34 and INTERPROSCAN35 programs were used to verify the predicted genes. As a result, we identified 47, 49, 92, and 116 Chis in G. raimondii, G. arboreum, G. hirsutum acc. TM-1, and G. barbadense acc. 3–79, respectively (Supplementary Table S1). From the phylogenetic view, one member in the diploid G. raimondii would correspond to one homologous gene in G. arboreum and two homeologs from the A and D subgenomes in tetraploid G. hirsutum acc. TM-1 and G. barbadense acc. 3–79. Here, 38 members showed such a one-to-one correspondence. Other inconsistencies might result from the different sequencing methods, assembly error in partial chromosomal regions, or duplication in the evolutionary process of Gossypium. In particular, in G. barbadense acc. 3–79, more tandem Chis were found compared with the other three cotton species and need to be further confirmed.

The nomenclature of Chis in Gossypium followed their chromosome orders. As the G. raimondii chromosomes had been integrated with the D-subgenome in G. hirsutum36, the Chis in G. raimondii were named preferentially from GrChi1 to GrChi47, based on the reordered G. raimondii chromosomes (Fig. 1). The corresponding orthologs in G. arboreum, G. hirsutum acc. TM-1, and G. barbadense acc. 3–79 were named GaChis, GhChis, and GbChis, respectively, with the same number for orthologs as in G. raimondii (Supplementary Table S1).

Figure 1
Chromosomal distribution of chitinase genes in G. raimondii.

Chis were matched to 12 scaffolds in G. raimondii except in the D4 chromosome. In detail, chitinase members were distributed among the Gossypium chromosomes unevenly, with D3, D7, and D9 chromosomes observed to have a relatively high density of Chis and other chromosomes containing sparse Chis. To further ascertain the expansion mechanism of the Chis in G. raimondii, we analyzed their tandem and segmental duplication events on the 12 chromosomes. Eleven tandem duplications involved in 30 Chis were detected in the G. raimondii genome, indicating that tandem duplications were the major factor in the expansion of more than 60% of members of the gene family. Moreover, 10 segmental duplication events in 47 Chis were identified. The results were quite different from those in Arabidopsis thaliana and Populus trichocarpa, with 24 Chis found in each chromosome in A. thaliana, but 37 Chis were distributed in only eight of 19 chromosomes in P. trichocarpa37.

Phylogenetic classification and structural analysis of chitinase genes in G. raimondii

With GrChis as example, the phylogenetic relationships were further analyzed by constructing an unrooted tree using the maximum likelihood method. The results indicated that the 47 GrChis were clustered into six well-supported groups; namely groups A, B, C, D, E, and F (Fig. 2a). These classifications were consistent with the presence of conserved domains and catalytic domains (Fig. 2b), and exon-intron organizations (Fig. 2c).

Figure 2
Phylogenetic classification and structural analysis of chitinase genes in G. raimondii.

Thirty chitinases belonging to the GH 18 family were divided into groups A, B, and C, and seventeen members of the GH 19 family were classified into groups D, E, and F. The genes in group A matched those in class III of the previous classifications3,4,5,6,7,8. The chitinase genes of groups B and C corresponded to class V. The class II genes were divided in groups D and E. The genes in class I were classified into group E, and the chitinase genes in group F corresponded to class IV (Fig. 2a).

The motif locations of chitinases are displayed schematically in Fig. 2b. Totally, ten motif components in the 47 GrChis were analyzed. Of them, GrChis in GH 18 gene family were involved in the motifs 1, 3, 4, 7, 8, 9 and 10, and those in GH 19 gene family were possessed of the motifs 2, 5 and 6, respectively. Moreover, we also analyzed the chitin-binding conserved domains in the 47 GrChis, while only eight GrChis belonging to GH 19 gene family was detected, implying that the chitin-binding conserved domain in most GrChis had lost in the evolutionary process. The type, order, number, and motif location of these chitinases were similar within each group, but significantly different to the other groups. In detail, 90% (9 out of 10) of the members in group A displayed the same five motif components (motifs 1, 3, 4, 7, and 8) and were ranked in an identical order, but had many obvious differences to the other five groups. The previously reported class V was further divided into groups B and C in the present study, with many differences in motif locations. Most of the GrChis in group B shared motifs 9 and 10, but GrChis in the group C were mainly motif 9. GrChis in groups D, E, and F contained the same three motifs (motifs 2, 5, and 6), which were distributed in the same order. This led to these 3 groups having a closer evolutionary and phylogenetic relationship. These findings indicate that the motif compositions of each group in chitinase gene family are relatively conserved (Supplemental data 1).

To gain further insights into the structural evolution and phylogenetic classification of the genes, the exon-intron organization of GrChis was also analyzed by comparing the corresponding genomic sequences with their coding sequences (CDS) through the software Gene Structure Display Server (GSDS) 2.0 ( (Fig. 2c). As a result, most genes in the same group generally shared exon-intron organizations with the same number of exons or introns. For example, seven of the genes in group A only had one exon, but the other three genes had two exons. However, the genes of the groups B and C, which were categorized into class V in previous studies, had distinct structural divergences. Most genes in group B only possessed one exon, while more than two exons were found in group C. The majority of genes in groups D and E had three exons, while most genes in group F had two exons.

Overall, the motif locations, motif components, and exon-intron structures of the chitinases in G. raimondii supported the phylogenetic classification well. Moreover, the highly conserved sequences of chitinase family members in the same group indicated these genes were subject to duplications during their evolution.

Expression profiles of chitinase genes in G. hirsutum acc. TM-1

To investigate the roles of chitinase genes in tetraploid cotton, transcriptome data from G. hirsutum acc. TM-1 vegetative tissues (root, stem, and leaf), floral tissues (petal and anther), ovule tissues (−3, 0, and 3 DPA), and fiber tissues (10, 20, and 25 DPA) at different developmental stages were used to investigate the gene expression patterns32. With FPKM > 1.0, a total of 59 GhChis were found to be expressed, and showed diverse developmental and spatial regulation in the various tissues (Fig. 3).

Figure 3
Transcriptional profiling of chitinase genes in different tissues and organs in G. hirsutum acc. TM-1.

Most of the chitinase genes in group A and group B were predominantly expressed in vegetative organs. GhChi47A in group C was ubiquitously expressed in both vegetative and reproductive organs, GhChi14A/D was undetectable in the ovule, GhChi40A/D was undetectable in the petal and stamen, and GhChi46D was expressed highly in the leaf, petal and the fiber tissues, while other members were only highly expressed in vegetative organs. GhChi32A/D in group D were expressed at high levels in both vegetative and reproductive organs, yet GhChi27A/D and GhChi37A/D were most highly expressed in the stem and fiber tissues. Moreover, eight genes in group E were predominantly expressed in vegetative organs, while GhChi17A/D had higher expression levels at 25DPA fibers. In addition, GhChi3A/D and GhChi4A/D were predominantly expressed in the ovule at different developmental stages. A total of ten genes in group F were detected. Of these, GhChi28A/D, GhChi30D, GhChi31A/D, GhChi34A/D had high expression levels in roots and leaves. Furthermore, the genes including GhChi28A/D, GhChi31A/D, and GhChi34A/D were also highly expressed in the reproductive organs. These findings indicate that chitinase genes share diverse but overlapping expression patterns in various tissues and organs, suggesting that chitinase genes have conserved structures but diverse functions.

Expression analysis of chitinase genes in response to V. dahlia

The expression patterns of chitinase genes have been studied in a range of species, and many genes have been demonstrated to function in resistance to biotic and abiotic stresses14,15,16. To confirm the induction of chitinase genes after V. dahliae inoculation, the expression patterns of 23 chitinase genes with detectable expression levels in root tissues following transcriptome analysis were investigated. Of them, eight genes showed high expression levels at various treatment time points, showing that they were significantly induced by V. dahliae, and quickly reached a peak at different time points (Fig. 4). In detail, the expression of six genes, Chi2, Chi14, Chi17, Chi23, Chi28, and Chi32, significantly increased after V. dahliae inoculation, reaching a peak at 96 h. The other two genes, Chi25 and Chi47, were significantly induced and quickly reached a peak at 48 h after V. dahliae treatment. Therefore, the eight genes might be responsible for V. dahliae resistance in cotton.

Figure 4
Real-time qRT-PCR analysis of the chitinase genes in response to Verticillium dahliae in G. barbadense cv. Hai7124.

Functional characterization of candidate genes in cotton responsive to V. dahliae through virus-induced gene silencing (VIGS)

By integrating the phylogenetic classification, the expression patterns in various organs and tissues, and in response to V. dahilae, the functions of chitinase genes are multiple and complex, with the important roles in disease resistance. To further investigate the function of Chis in V. dahliae resistance, we selected three candidate genes, Chi23, Chi32, and Chi47, with different expression levels that were significantly induced after V. dahliae inoculation, for functional identification through VIGS analysis.

We constructed TRV: Chi23, TRV: Chi32, and TRV: Chi47 vectors to silence endogenous genes in Hai7124, with TRV: 00 considered as the mock treatment. To validate the viability of VIGS in cotton, CLOROPLASTOS ALTERADOS 1 (CLA1), which encodes 1-deoxy-D-xylulose-5- phosphate synthase, was used as an indicator gene and was silenced to produce plants with a bleached leaf phenotype38. All constructs were infiltrated into at least 20 Hai7124 seedlings at 8 days post-emergence, and untreated plants were also cultured in the same environment. Two weeks later, all individuals infiltrated with TRV: CLA1 showed highly uniform bleaching in newly emerged leaves (Fig. S1). Meanwhile, the cotton seedlings confirmed to be infiltrated with all constructs and those with the mock treatment were selected for RNA isolation and quantitative real-time PCR (qPCR) analysis. The untreated and mock-treated plants showed the same or higher expression levels of Chi23, Chi32, and Chi47. However, the transcripts of these three genes showed that they were significantly silenced in plants infiltrated with TRV: Chi23, TRV: Chi32, and TRV: Chi47 (P < 0.01) (Fig. 5a).

Figure 5
Silencing of Chi23, Chi32, and Chi47 significantly impaired the resistance to V. dahliae in cotton.

After the target genes were silenced, we inoculated the cotton seedlings with V. dahliae isolate V991 by dip-infection with a final concentration of 1 × 107 spores per milliliter. In parallel, two cotton cultivars, Hai7124 and Junmian 1, were used as controls resistant to and susceptible to V. dahliae, respectively, and these were also inoculated with V991. Ten days after inoculation, the seedlings of the Junmian 1 plants showed obvious cotyledon wilting. However, a small number of Hai7124 plants with the leaf wilting phenotype appeared at least 15 days after inoculation with V. dahliae. As a result, Hai7124 and Junmian 1 plants were confirmed to be resistant and susceptible to V. dahliae, respectively. The phenotypes of the two cotton cultivars at 20 and 25 days after inoculation with V. dahliae are shown in Fig. S2.

About three weeks after inoculation, the cotton seedlings of gene-silenced plants, particularly those containing TRV: Chi32 and TRV: Chi47, exhibited more wilting and etiolated leaves than the control TRV: 00 seedlings (Fig. 5b). Similarly, all the true leaves in the Junmian 1 seedlings were defoliated 25 days after inoculation (Fig. S2). The defoliated phenotype was also observed in TRV: Chi32 and TRV: Chi47 seedlings 30 days after inoculation. In summary, silencing Chi23, Chi32, and Chi47 obviously increased susceptibility to V. dahliae in Hai7124 seedlings. In addition, we used at least 20 plants per treatment to calculate the rate of diseased plants. As a result, the TRV: 00 control plants had few wilted leaves, similar to the number seen in Hai7124 seedlings without injection, and the rate of average diseased leaves to healthy leaves was approximately 55% 35 days after inoculation. However, about 95% of the Chi32- silenced and Chi47-silenced plants showed leaves wilting or exfliation that were similar to the Junmian 1 susceptible control plants, of which 100% were diseased 35 days after inoculation (Supplemental data 2). Furthermore, about 80% of the Chi23-silenced plants were diseased: this was higher than the percentage of control plants, but lower than the percentage of Chi32-silenced, Chi47-silenced, and Junmian 1 plants (Fig. 5c). By comparing the rate of disease in treated and control plants, the silencing of Chi23, Chi32, and Chi47 were found to significantly confer the resistance to V. dahliae (P < 0.01). These results indicate that Chi32 and Chi47 are important genes in resistance to V. dahliae infection in cotton, and have potential uses in cotton disease-resistant breeding.

The expression patterns of chitinase genes are also regulated by various plant phytohormones, including JA, SA, ET, CK, and IAA19. To further investigate the molecular mechanisms of Chi23, Chi32 and Chi47, we analyzed the expression of four stress-related signaling compounds, ET, JA, SA, and the hydrogen peroxide (H2O2). Chi32 and Chi47 were found to be induced by one or more of these compounds. The transcript levels of Chi32 significantly increased and accumulated 4–24 h after JA treatment. Chi47 was also significantly upregulated and reached a peak 12 h after JA or H2O2.treatment (Fig. 6). These results imply that Chi32 and Chi47 participate in JA and other signaling pathways to protect plants from various stresses.

Figure 6
Expression patterns of two chitinase genes treated with stress-related signals in cotton.


Chitinases are found in a variety of species including plants, animals, fungi, and bacteria1. They not only play a major role in defense against fungal pathogens16, but also have an important function in the regulation of growth and development in plants13,14. As a larger gene family, it is important to elucidate the multiple functions of chitinase genes and to determine their corresponding functional specificity for utilization in genetic analyses and breeding. To date, systematic genome-wide investigation of chitinase genes has been reported in many species; however, knowledge of the chitinase genes in cotton is limited, and their systematic investigation has not been reported.

In Arabidopsis, 24 chitinase genes have been found, distributed between all five of the chromosomes39. In rice, 37 chitinase genes have been found and these are positioned on all chromosomes except chromosome 739. A model woody plant, P. trichocarpa, contains 37 chitinase genes with complete ORFs and these are located on eight of the 19 chromosomes37. A total of 26 chitinase genes were found in the diploid banana (AA) genome40. Here, we systematically identified 47, 49, 92, and 116 candidate chitinase genes in four sequenced cotton species; the diploid cottons G. raimondii and G. arboreum, and the tetraploid cottons, G. hirsutum acc. TM-1 and G. barbadense acc. 3–79, respectively. It is noted that more tandem Chis were found in G. barbadense acc. 3–79 than the other three cotton species, which might result from the different sequencing methods used, assembly errors in partial chromosomal regions, or evolutionary expansion, and this should be further confirmed. In the diploid cotton G. raimondii, the chitinase genes were anchored onto 12 chromosomes and could be divided into six subgroups. Characteristics of the genes, such as signal peptides (SPs), GPI-anchor site, and transmembrane domains, were investigated, and the results indicated multiple complex physiological functions in the defense against various stresses and the enhancement of growth and development in cotton.

In accordance with previous reports, the transcription patterns of plant chitinase genes showed both constitutive and induced expression in the present study14,15,16,26. Constitutively expressed chitinases mainly participate in physiological functions such as flower development, symbiotic interactions, seed development, somatic embryogenesis, cell division, and programmed cell death14,41. However, induced expression of chitinases have been recorded in response to numerous abiotic stresses, such as wounding, drought, salinity, cold, heavy metals, ozone, frost, and UV light15,16,17,18, as well as biotic stresses, such as infection with fungi, bacteria, viruses, viroids and nematodes24,25,26,27,28,41. Furthermore, the expression of chitinases is regulated by various phytohormones, such as JA, SA, ET, IAA, and CK19.

The variation in the expression patterns of the chitinase genes is related to their functional diversity. For example, a rice chitinase III protein, drought-induced protein 3 (DIP3), was significantly induced by drought, salt, and low temperature, and functions by regulating the plant stress response42. In grapevine, the expression levels of Chi1b and CH3 rapidly increased after UV-C exposure, and the chitinase activity rose between 4- to 13-fold in bunchstems43. ScChi, an acidic class III chitinase in sugarcane, showed different expression levels when treated with H2O2, methyl jasmonate (MeJA), SA, abscisic acid (ABA), polyethylene glycol (PEG), NaCl, CuCl2, and low temperature (4 °C), suggesting that it confers positive responses to biotic and abiotic stresses44. In this study, we systematically analyzed the expression patterns of chitinase family genes in cotton. Using the transcriptome data of G. hirsutum acc. TM-1 vegetative tissues (root, stem, and leaf), floral tissues (petal and anther), ovule tissues (−3, 0, and 3 DPA), and fiber tissues (10, 20, and 25 DPA) at different developmental stages, a total of 59 GhChis were found to have diverse developmental and spatial regulation patterns (Fig. 3). Further, we investigated the expression of 23 GhChis with detectable expression in root tissues after V. dahliae inoculation. Eight genes were found to have significantly higher expression levels following inoculation and were therefore thought to be induced by V. dahliae (Fig. 4). These findings indicate that chitinase genes have diverse functions. It is hoped that these data will help mine key genes in various breeding targets such as V. dahliae resistance in cotton.

Chitinases have two different mechanisms for the hydrolysis of chitins: substrate-assisted catalysis and acid catalysis45,46. Although plants have no chitin in their cell walls, chitinase genes play key roles in multiple physiological processes, including growth and development and non-specific stress responses. Plants produce many forms of chitinases, including vacuolar and apoplastic chitinases, in the natural environment or in response to stresses. The production of chitinases had been shown to relate to some biotic or abiotic factors15,47. In addition, many chitinase genes sourced from plants and fungi have been successfully over-expressed in other plants to combat stresses. Transgenic tobacco (Nicotiana tabacum) plants that overexpress CHIT33 and CHIT42 from the mycoparasitic fungus Trichoderma harzianum, show broad resistance to biotic stresses (fungal and bacterial pathogens) and abiotic stresses (salinity and heavy metals)48. Cht-2 and Cht-3, class-I chitinase genes, are constitutively expressed in japonica rice to enhance resistance to blast (Magnaporthe grisea)26. Transgenic banana (Musa acuminata ‘Gros Michel’) with either of two rice chitinase genes (rcc2 or rcg3) showed enhanced resistance to black leaf streak disease caused by the fungus Mycosphaerella fijiensis27. A rice chitinase gene (Cht-2; RCC2) was introduced into the calli of Italian ryegrass (Lolium multiflorum Lam.), and dramatically increased its resistance to crown rust disease49. A transgenic peanut (Arachis hypogaea L.) expressing tobacco chitinase was resistant to leaf spot disease50. An endochitinase (ech42) was cloned from T. harzianum and introduced into apple to increase resistance to apple scab51. Ech42 was also isolated from the biocontrol fungus Trichoderma virens and successfully expressed in tomato and tobacco to enhance their tolerance to fungal pathogens28. In addition, the transient co-expression of class IV chitinase in pepper (CaChitIV) and pepper receptor-like cytoplasmic protein kinase (CaPIK1) in pepper plants (Capsicum annuum) enhanced the CaPIK1-triggered cell death response, the abundance of reactive oxygen species (ROS), and the number of nitric oxide (NO) bursts. In contrast, co-silencing of both CaChitIV and CaPIK1 reduced induction of the cell death response, ROS and NO bursts, and defense response genes, and hence, increased susceptibility to Xanthomonas campestris pv. vesicatoria (Xcv) infection52. Co-expression of a rice basic chitinase gene (RCH10) and a modified maize ribosome-inactivating protein (MOD1) in transgenic rice significantly enhanced resistance to sheath blight25. Interestingly, transgenic plants expressing chitinase genes in combination with other PR proteins, such as a modified maize ribosome-inactivating protein or β-1, 3-glucanases, showed enhanced resistance to fungal pathogens compared to those where a single chitinase gene was introduced25. Together, these results suggest that chitinase genes play important roles in growth, development, and non-specific stress responses, and play a particular role in the defense against fungal pathogens.

In this study, based on genome-wide identification, phylogenetic classification, and expression patterns, we further validated the roles of several chitinase genes in cotton V. dahliae resistance. We silenced the Chi23, Chi32, and Chi47 genes in G. barbadense cv. Hai7124 and found through VIGS analysis that it led to Verticillium resistance. Two weeks after infiltration with TRV: Chi23, TRV: Chi32, and TRV: Chi47, seedlings exhibited obvious silencing compared with the controls, and the silencing of these genes significantly enhanced cotton susceptibility to V. dahliae (Fig. 5). Statistical analysis suggested that silencing Chi23, Chi32, or Chi47 increased the susceptibility of Hai7124 to V. dahliae. Interestingly, the Chi32-silenced and Chi47-silenced plants showed particularly serious Verticillium wilt that was similar to that seen in the susceptible plant, Junmian 1. These findings imply that Chi32 and Chi47 are important defense-related genes and play crucial roles in resistance to V. dahliae in cotton. These systematic analyses of the chitinase family in cotton have provided crucial information on the molecular mechanisms of disease resistance, and provide a solid foundation for disease-resistant cotton breeding.

From the evolutionary view, G. hirsutum and G. barbadense probably originated from a single hybridization event between A- and D- diploid species, however, the two have very different agronomic and fiber quality characteristics53. Modern G. hirsutum cultivars have high-yield properties and dominate more than 90% of worldwide cotton production, while G. barbadense, characterized by its elite fiber qualities and resistance to V. dahliae, accounts for less than 10%32,33. Even the high similarity of orthologous genes in the two cotton species, the differences involved in the single nucleotide polymorphisms (SNP), expression level, proteins structures, and the variance in epigenetics significantly affected their phenotype characteristics54,55. Here, we found that Chi23, Chi32 and Chi47 were significantly upregulated at different expression levels in cotton roots after inoculation with V. dahlia, with higher expression in G. barbadense cv. Hai7124 than G. hirsutum cv. Junmian 1 (Fig. S3). Further function analysis showed that silencing of these three genes in Hai7124, especially Chi32 or Chi47, significantly increased the susceptibility of the cotton to V. dahliae. Taken together, our results suggested that resistant genes could be effectively mined from resistant materials, and these Chis might act as positive regulators in disease resistance to V. dahliae. Their function in disease resistance will be further elucidated by transgenic expression analysis in susceptible cotton cultivars.

The chitinase genes exhibit multiple complex functions in development and stress responses, and seem to be mainly involved in active or passive defense against pathogens. However, the production of chitinases is also regulated by the generation or degradation of signal molecules such as JA, ET, SA, CK, H2O2, and auxin19. In the bark and sapwood of Norway spruce (Picea abies L.), the expression levels of a class IV chitinase gene (PaChi4) was obviously increased when treated with the necrotropic pathogen Heterobasidion parviporum and MeJA56. In tomato (Solanum lycopersicon L.), SlERF1, an ethylene response factor (ERF), was engineered for Rhizopus nigricans resistance by accumulating chitinase (CHI) and phenylalanine ammonia lyase (PAL)57. In Casuarina equisetifolia, the expression of CeChi1 were induced by pathogen infection or SA treatment58. Ltchi7, a class III chitinase, was significantly induced in roots of Lotus tenuis and L. japonicas when treated with salt, drought, H2O2, and ABA59. Among these endogenous signaling molecules, the SA-mediated defense responses were effective against hemi-biotrophs and biotrophs, and were critical for systemic acquired resistance, whereas JA tended to act with ET to induce resistance to necrotrophic pathogens60.

Many signaling molecules such as the JA, ET, SA, and H2O2 positively contributed to the disease resistance to V. dahliae in cotton38,61. Accumulating evidence indicated that V. dahliae-responsive genes such as GbSSI238, GhSSN62, GbERF63, and GbWARKY164, which were involved in regulating the biosynthesis of signaling molecules, played crucial roles in defense against Verticillium wilt. In this study, the expression of Chi23, Chi32 and Chi47 was further investigated following treatment with the stress-related signals ET, JA, SA, and H2O2. As a result, the expression levels of Chi32 and Chi47 were significantly increased after JA treatment, and the expression of Chi47 was also upregulated following H2O2 treatment. Expression analysis showed that Chi32 and Chi47 may be involved in resistance to stress through regulation of JA or H2O2. This, along with VIGS analysis, suggests that JA acts as a key factor in engineering the chitinase gene-mediated resistance to V. dahliae in cotton; however the detailed functional mechanisms remains to be investigated.

Materials and Methods

Identification of the chitinase gene family in four sequenced cotton species

The available genomic database of four sequenced cotton species, G. raimondii, G. arboreum, G. hirsutum acc. TM-1, and G. barbadense acc. 3–79, were downloaded from,,, and, respectively. We used the Pfam protein family databases with the Glyco_hydro_18 (PF00704) and Glyco_hydro_19 (PF00182)65 domains and the HMMER software version 3.066 to identify the chitinase genes in the four cotton species. Further, the chitinase genes were isolated from the four databases using SeqHunter software version 1.0. In addition, the sequences were further confirmed the existence of the conserved domains of Glyco_hydro_18 or Glyco_hydro_19 using SMART34 and INTERPROSCAN35.

Mapping, phylogenetic tree construction, and structural analysis

The distribution of chitinase genes in Gossypium was determined by MapInspect (, where the chromosome orders from D1 to D13 and A1 to A13 referred to the map published by Wang et al.36. G. raimondii chromosomes were integrated with the D-subgenome in G. hirsutum36 and reordered from D1 to D13. The nomenclature of chitinase genes (GrChis) in G. raimondii was based on the order of the chromosomes. The corresponding orthologs in G. arboreum, G. hirsutum acc. TM-1, and G. barbadense acc. 3–79 were named GaChis, GhChis, and GbChis, respectively, following their orthologs in G. raimondii.

Multiple sequence alignments of chitinase proteins were carried out using ClustalX (ver.1.83). A phylogenetic tree was built by MEGA 5.1 software ( using the maximum likelihood method. The exon/intron structures were analyzed by aligning the genomic DNA sequences with their corresponding coding sequences using the online Gene Structure Display Server (GSDS) program ( Motifs of chitinase proteins were investigated statistically using MEME (, which set the maximum number of motifs at 10.

Plant materials and treatments

Seeds of G. barbadense cv. Hai7124, G. hirsutum cv. Junmian 1, and G. hirsutum cv. Jinmian 19 were used for expression analysis under different stress treatments. They were grown in the different controlled environment chambers under the same conditions of a 16 h light/8 h dark cycle at 28 °C for 2 weeks.

The seedlings of Hai7124 and Junmian 1, which exhibited resistance and susceptibility to V. dahliae, respectively, were inoculated with the fungal pathogen (V. dahliae) with the dip-inoculation method67. V991, a highly aggressive defoliating strain of V. dahliae, was cultured on potato dextrose agar medium at 24 °C for 4–5 d, and then incubated in Czapek’s medium (NaNO3, 0.3% w/v; MgSO4, 0.1% w/v; KH2PO4, 0.1% w/v; FeSO4, 0.0002% w/v; KCl, 0.1% w/v; Sucrose, 3% w/v; pH6.0) at 25 °C for 5 d. We then adjusted the concentration to 107 conidia per ml using deionized water for inoculation of the seedlings38. The seedling roots were harvested with three repeats at different time points (0, 24, 48, 96, and 144 hours) after V991 treatment, and then quick-frozen in liquid nitrogen and stored at −80 °C for use.

G. hirsutum cv. Jinmian 19, which exhibits high tolerance to abiotic stress, was used for the stress signal treatments. The leaves of three-week-old cotton seedlings were sprayed with 100 μM JA, 1 mM ET, 100 mM SA, or 10 mM H2O2 (ddH2O as a solvent control), and the leaves were collected at different time points (0, 1, 4, 12, and 24 hours). Three repeats were harvested at each time point and stored at −80 °C after quick-freezing in liquid nitrogen.

RNA isolation and expression pattern analysis

Total RNA was extracted from cotton seedling leaves or roots using the CTAB-acidic phenol extraction method68. RNA was then treated with DNase I (Invitrogen, to remove genomic DNA, and 2 μg of total RNA was used for first-strand cDNA synthesis using the Superscript first-strand synthesis system (Invitrogen, Foster City, CA).

The primer pairs used for qRT- PCR analysis were designed by Beacon Designer 7.0 according to the chitinase gene sequences in G. raimondii (Supplementary Table S2). The amplified fragment lengths were set between 80 base pair (bp) and 200 bp, and the annealing temperature was controlled at 60 °C. The cotton histone3 (AF024716) gene was used as the reference gene for qRT- PCR.

The real-time PCR amplification reactions were performed on an ABI 7500 Real Time PCR System (Applied Biosystems, USA) using SYBR Green (Vazyme, Nanjing, China) with three replicates. The amplification parameters were: denaturation at 95 °C for 10 min, 40 cycles of denaturation at 95 °C for 15 s, annealing at 60 °C for 15 s, and extension at 72 °C for 15 s. Expression data from three biologically independent experiments were analyzed and presented as means ± S.D.

To analyze the expression of chitinase genes, we obtained G. hirsutum acc. TM-1 high-throughput RNA-sequencing data from Zhang et al.32 from vegetative tissues (root, stem, and leaf), floral tissues (petal and stamen), ovule, and fiber tissues at −3, 0, 3, 10, 20, and 25 days post anthesis (DPA), respectively. Finally, Log2 (FPKM) was used to calculate GhChis expression levels, where FPKM referred to fragments per kilobase of exon model per million mapped reads with Cufflinks software (

Cloning of the chitinase genes in G. barbadense cv. Hai7124

Based on the known sequences of four cotton species, we designed gene-specific primers to amplify the homologous genes of chitinases with complete open reading frames (ORFs) in G. barbadense cv. Hai7124, a species resistant to V. dahliae (Supplementary Table S2). Standard PCR reactions were performed using High-fidelity ExTaq DNA Polymerase (TaKaRa Biotechnology [Dalian] Co. Ltd., China). All the PCR products were cloned into pMD19-T cloning Vector (TaKaRa, Dalian, China) and transformed into E. coli DH5α. At least six clones per gene were randomly selected and sequenced.

Vector construction and functional characterization of candidate genes via virus-induced gene silencing (VIGS) analysis

The pTRV1 and pTRV vectors for VIGS analysis were generously provided by Dr. Libo Shan of Texas A & M University (College Station, TX, USA). pTRV: GhCLA1, where GhCLA1 (Cloroplastos alterados 1) encodes 1-deoxy-D-xylulose-5-phosphate synthase, was used as a control69. The TRV: Chi23 construct contained a 444-bp fragment of the Chi23 cDNA, the TRV: Chi32 construct included a 454-bp fragment of Chi32, and the TRV: Chi47 construct included a 387-bp fragment of Chi47. These fragments were cloned from Hai7124 by PCR using primers with EcoRI/SacI enzyme sites for TRV: Chi23 and TRV: Chi32, and EcoRI/BamHI enzyme sites for TRV: Chi47, for insertion into TRV. The primer pairs used for constructing the VIGS vectors are displayed in Table S2.

We transformed the vectors, TRV1, TRV: 00, TRV: GhCLA1, TRV: Chi23, TRV: Chi32, and TRV: Chi47 into Agrobacterium tumefaciens strain GV3101, and cultured these Agrobacterium at 28 °C in LB medium containing the antibiotics 50 μg/mL kanamycin and 25 μg/mL rifampicin for 2 days. We then cultured the Agrobacterium in liquid LB medium (50 μg/mL kanamycin, 25 μg/mL rifampicin) on a shaking table with 160 rpm per minute at 28 °C for 16 hours. Finally, we inoculated the Agrobacterium at a concentration of 1:100 into the liquid LB medium that contained the antibiotics rifampicin (25 μg/mL) and kanamycin (50 μg/ mL), shook the mixture, and allowed the Agrobacterium to grow to an OD600 of 0.5 at 28 °C. The Agrobacterium was pelleted using centrifugation at 4000 rpm at room temperature for 10 min, and the OD600 value was adjusted to 2.0 with infiltration media (10 mM MgCl2, 10 mM 2-Morpholinoethanesulfonic acid, and 200 μM acetosyringone).

The cell suspensions carrying TRV1, TRV: 00, TRV: Chi23, TRV: Chi32, and TRV: Chi47 were incubated at room temperature for 3 h, and then infiltrated into the eight-day-old Hai7124 cotton seedlings with two fully expanded cotyledons using a needleless 1 mL syringe at a 1:1 ratio. For mock treatment and the technical control, the same plants were infiltrated with a 1:1 mixture of Agrobacterium carrying TRV1 and TRV: 00 or TRV1 and TRV: CLA1, respectively. All the plants were grown in the same growth chamber at 23/21 °C (day/night), with a 16 h light/8 h dark cycle. About 8 days after Agrobacterium infiltration, the seedlings infiltrating with TRV: CLA1 showed a bleached leaf phenotype. Two weeks after Agrobacterium infiltration, the TRV: CLA1 plants exhibited highly uniform bleaching in newly emerged leaves, and then we selected three seedlings per treatment and harvested the leaves of the VIGS plants for the RNA isolation. Meanwhile, the cotton seedlings containing the susceptible control Junmian 1, the resistant control Hai7124, and all the VIGS plants were removed from the soil and dip-inoculated with V. dahliae conidia suspension (1 × 107conidia/mL) as previously described67. All plants were grown in the same growth chamber at 25/23 °C (day/night), with a 16 h light/8 h dark cycle for seven weeks to investigate the rate of the diseased leaves to healthy leaves. The VIGS experiments were repeated at least three times and each treatment was applied to more than 20 plants to increase the reliability of the results.

Additional Information

How to cite this article: Xu, J. et al. Discovery and identification of candidate genes from the chitinase gene family for Verticillium dahliae resistance in cotton. Sci. Rep. 6, 29022; doi: 10.1038/srep29022 (2016).

Supplementary Material

Supplementary Information:
Supplementary Table S1:
Supplementary Table S2:


This program was financially supported in part by National Natural Science Foundation in China (31171590), the National Transgenic Program (2011ZX08005-004), Key R&D program in Jiangsu Province (BE2015360), Six talent peaks project in Jiangsu province (2015-NY-002), the Priority Academic Program Development of Jiangsu Higher Education Institutions (010-809001), and JCIC-MCP project (No. 10).


Author Contributions Experiments were designed by W.G. Experiments were performed by J.X., X.X., L.T., G.L.W., X.Z. and X.W. J.X., X.W. and W.G. drafted the manuscript, W.G. and X.W. revised the manuscript. All authors read and approved the final manuscript.


  • Lin W. et al. . Genetic engineering of rice for resistance to sheath blight. Nat. Biotech. 13, 686–691 (1995).
  • Li H. & Greene L. H. Sequence and structural analysis of the chitinase insertion domain reveals two conserved motifs involved in chitin-binding. PLoS One, 5, e8654 (2010). [PMC free article] [PubMed]
  • Shinshi H., Neuhaus J. M., Ryals J. & Meins F. Jr. Structure of a tobacco endochitinase gene: evidence that different chitinase genes can arise by transposition of sequences encoding a cysteine-rich domain. Plant Mol. Biol. 14, 357–368 (1990). [PubMed]
  • Collinge D. B. et al. . Plant chitinases. Plant J. 3, 31–40 (1993). [PubMed]
  • Beintema J. J. Structural features of plant chitinases and chitin binding proteins. FEBS Lett. 350, 159–163 (1994). [PubMed]
  • Melchers L. S. et al. . A new class of tobacco chitinase homologous to bacterial exo-chitinases displays antifungal activity. Plant J. 5, 469–480 (1994). [PubMed]
  • Hamel F., Boivin R., Tremblay C. & Bellemar G. Structural and evolutionary relationships among chitinases of flowering plants. J. Mol. Evol. 44, 614–624 (1997). [PubMed]
  • Passarinho P. & de Vries S. C. Arabidopsis chitinases: a genomic survey. Arabidopsis book 1, e0023 (2002). [PMC free article] [PubMed]
  • Renner T. & Specht C. D. Molecular and functional evolution of class I chitinases for plant carnivory in the Caryophyllales. Mol. Biol. Evol. 29, 2971–2985 (2012). [PubMed]
  • Jekel P. A., Hartmann B. H. & Beintema J. J. The primary structure of hevamine, an enzyme with lysozyme/chitinase activity from Hevea brasiliensis latex. Eur. J. Biochem. 200, 123–130 (1991). [PubMed]
  • Zhao K. J. & Chye M. L. Methyl jasmonate induces expression of a novel Brassica juncea chitinase with two chitin-binding domains. Plant Mol. Biol. 40, 1009–1018 (1999). [PubMed]
  • Punja Z. K. & Zhang Y. Y. Plant chitinases and their roles in resistance to fungal diseases. J. Nematol. 25, 526–540 (1993). [PMC free article] [PubMed]
  • van der Holst P. P., Schlaman H. R. & Spaink H. P. Proteins involved in the production and perception of oligosaccharides in relation to plant and animal development, Curr. Opin. Struct. Biol. 11, 608–616 (2001). [PubMed]
  • Passarinho P. A., Van Hengel A. J., Fransz P. F. & de Vries S. C. Expression pattern of the Arabidopsis thaliana AtEP3/AtchitIV endochitinase gene. Planta 212, 556–567 (2001). [PubMed]
  • Hamel F. & Bellemare G. Characterization of a class I chitinase gene and of wound-inducible, root and flower-specific chitinase expression in Brassica napus. Biochim. Biophys. Acta. 1263, 212–220 (1995). [PubMed]
  • Cohen-Kupiec R. & Chet I. The molecular biology of chitin digestion. Curr. Opin. Biotechnol. 9, 270–277 (1998). [PubMed]
  • Hanfrey C., Fife M. & Buchanan-Wollaston V. Leaf senescence in Brassica naupus: expression of genes encoding pathogenesis related proteins. Plant Mol. Biol. 30, 595–609 (1996).
  • Bishop J. G., Dean A. M. & Mitchell-Olds T. Rapid evolution in plant chitinases: molecular targets of selection in plant-pathogen coevolution. Proc. Natl. Acad. Sci. USA 97, 5322–5327 (2000). [PubMed]
  • Kasprzewska A. Plant chitinases-regulation and function. Cell. Mol. Biol. Lett. 8, 809–824 (2003). [PubMed]
  • Dangl J. L. & Jones J. D. Plant pathogens and integrated defence responses to infection. Nature 411, 826–833 (2001). [PubMed]
  • Klosterman S. J., Atallah Z. K., Vallad G. E. & Subbarao K. V. Diversity, pathogenicity, and management of verticillium species. Annu. Rev. Phytopathol. 47, 39–62 (2009). [PubMed]
  • Bhat R. G. & Subbarao K. V. Host range specificity in Verticillium dahliae. Phytopathology 89, 1218–1225 (1999). [PubMed]
  • Leah R., Tommerup H., Svendsen I. & Mundy J. Biochemical and molecular characterization of three barley seed proteins with antifungal properties. J. Biol. Chem. 266, 1564–1573 (1991). [PubMed]
  • Cao R. et al. . Mycoparasitism of endophytic fungi isolated from reed on soilborne phytopathogenic fungi and production of cell wall-degrading enzymes in vitro. Curr. Microbiol. 59, 584–592 (2009). [PubMed]
  • Kim J. K. et al. . Co-expression of a modified maize ribosome-inactivating protein and a rice basic chitinase gene in transgenic rice plants confers enhanced resistance to sheath blight. Transgenic Res. 12, 475–484 (2003). [PubMed]
  • Nishizawa Y. et al. . Enhanced resistance to blast (Magnaporthe grisea) in transgenic Japonica rice by constitutive expression of rice chitinase. Theor. Appl. Genet. 99, 383–390 (1999). [PubMed]
  • Kovács G. et al. . Expression of a rice chitinase gene in transgenic banana (‘Gros Michel’, AAA genome group) confers resistance to black leaf streak disease, Transgenic Res. 22, 117–130 (2013). [PMC free article] [PubMed]
  • Shah M. R., Mukherjee P. K. & Eapen S. Expression of a fungal endochitinase gene in transgenic tomato and tobacco results in enhanced tolerance to fungal pathogens. Physiol. Mol. Biol. Plants 16, 39–51 (2010). [PMC free article] [PubMed]
  • Emani C. et al. . Enhanced fungal resistance in transgenic cotton expressing an endochitinase gene from Trichoderma virens. Plant Biotechnol. J. 1, 321–336 (2003). [PubMed]
  • Wang K. et al. . The draft genome of a diploid cotton Gossypium raimondii. Nat. Genet. 44, 1098–1103 (2012). [PubMed]
  • Li F. et al. . Genome sequence of the cultivated cotton Gossypium arboreum. Nat. Genet. 46, 567–572 (2014). [PubMed]
  • Zhang T. et al. . Sequencing of allotetraploid cotton (Gossypium hirsutum L. acc. TM-1) provides a resource for fiber improvement. Nat. Biotech. 33, 531–537 (2015). [PubMed]
  • Yuan D. et al. . The genome sequence of Sea-Island cotton (Gossypium barbadense) provides insights into the allopolyploidization and development of superior spinnable fibres. Sci. Rep. 5, 17662 (2015). [PMC free article] [PubMed]
  • Letunic I., Doerks T. & Bork P. SMART 7: recent updates to the protein domain annotation resource. Nucleic. Acids. Res. 40, D302–305 (2012). [PMC free article] [PubMed]
  • Hunter S. et al. . InterPro in 2011: new developments in the family and domain prediction database. Nucleic. Acids. Res. 40, D306–312 (2012). [PMC free article] [PubMed]
  • Wang S. et al. . Sequence-based ultra-dense genetic and physical maps reveal structural variations of allopolyploid cotton genomes. Genome Biol. 16, 108 (2015). [PMC free article] [PubMed]
  • Jiang C., Huang R. F., Song J. L., Huang M. R. & Xu L. A. Genome-wide analysis of the chitinase gene family in Populus trichocarpa. J. Genet. 92, 121–125 (2013). [PubMed]
  • Gao W. 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 12, 3690–3703 (2013). [PMC free article] [PubMed]
  • Xu F., Fan C. & He Y. Chitinases in Oryza sativassp. Japonica and Arabidopsis thaliana. J. Genet. Genomics 34, 138–150 (2007). [PubMed]
  • Backiyarani S. et al. . Genome-wide analysis and differential expression of chitinases in banana against root lesion nematode (Pratylenchus coffeae) and eumusa leaf spot (Mycosphaerella eumusae) pathogens. Appl. Biochem. Biotechnol. 175, 3585–3598 (2015). [PubMed]
  • Cletus J., Balasubramanian V., Vashisht D. & Sakthivel N. Transgenic expression of plant chitinases to enhance disease resistance. Biotechnol. Lett. 35, 1719–1732 (2013). [PubMed]
  • Guo X. L., Bai L. R., Su C. Q., Shi L. R. & Wang D. W. Molecular cloning and expression of drought-induced protein 3 (DIP3) encoding a class III chitinase in upland rice. Genet. Mol. Res. 12, 6860–6870 (2013). [PubMed]
  • Petit A. N. et al. . Low responsiveness of grapevine flowers and berries at fruit set to UV-C irradiation. J. Exp. Bot. 60, 1155–1162 (2009). [PubMed]
  • Su Y. et al. . ScChi, encoding an acidic class III chitinase of sugarcane, confers positive responses to biotic and abiotic stresses in sugarcane. Int. J. Mol. Sci. 15, 2738–2760 (2014). [PMC free article] [PubMed]
  • van Aalten D. M. et al. . Structural insights into the catalytic mechanism of a family 18 exo-chitinase. Proc. Natl. Acad. Sci. USA 98, 8979–8984 (2001). [PubMed]
  • Brameld K. A. & Goddard W. A. 3rd. The role of enzyme distortion in the single displacement mechanism of family 19 chitinases. Proc. Natl. Acad. Sci. USA 95, 4276–4281 (1998). [PubMed]
  • Chen X. L. et al. . Cloning and functional analysis of the chitinase gene promoter in peanut. Genet. Mol. Res. 14, 12710–22 (2015). [PubMed]
  • de las Mercedes Dana M., Pintor-Toro J. A. & Cubero B. Transgenic tobacco plants overexpressing chitinases of fungal origin show enhanced resistance to biotic and abiotic stress agents. Plant Physiol. 142, 722–730 (2006). [PubMed]
  • Takahashi W. et al. . Increased resistance to crown rust disease in transgenic Italian ryegrass (Lolium multiflorum Lam.) expressing the rice chitinase gene. Plant Cell Rep. 23, 811–818 (2005). [PubMed]
  • Rohini V. K. & Sankara Rao K. Transformation of peanut (Arachis hypogaea L.) with tobacco chitinase gene: variable response of transformants to leaf spot disease. Plant Sci. 160, 889–898 (2001). [PubMed]
  • Bolar J. P. et al. . Expression of endochitinase from Trichoderma harzianum in transgenic apple increases resistance to apple scab and reduces vigor. Phytopathology 90, 72–77 (2000). [PubMed]
  • Kim D. S., Kim N. H. & Hwang B. K. The Capsicum annuum class IV chitinase ChitIV interacts with receptor-like cytoplasmic protein kinase PIK1 to accelerate PIK1-triggered cell death and defence responses. J. Exp. Bot. 66, 1987–1999 (2015). [PMC free article] [PubMed]
  • Brubaker C. L., Bourland F. M. & Wendel J. F.: The origin and domestication of cotton. In Cotton: origin, history, technology and production. Volume 1.1 edition. Edited by: Smith W. C., Cothren T. J.New York, John Wiley and Sons, 3–31 (1999).
  • Zhu H. Y. et al. . Structure, expression differentiation and evolution of duplicated fiber developmental genes in Gossypium barbadense and G. hirsutum. BMC Plant Biol. 11, 40 (2011). [PMC free article] [PubMed]
  • Chen Z. J. Genomic and epigenetic insights into the molecular bases of heterosis. Nat. Rev. Genet. 14, 471–82 (2013). [PubMed]
  • Yaqoob N. et al. . Defence-related gene expression in bark and sapwood of Norway spruce in response to Heterobasidion parviporum and methyl jasmonate. Physiol. Mol. Plant Pathol. 77, 10–16 (2012).
  • Pan X. Q., Fu D. Q., Zhu B. Z., Lu C. W. & Luo Y. B. Overexpression of the ethylene response factor SlERF1 gene enhances resistance of tomato fruit to Rhizopus nigricans. Postharv. Biol. Technol. 75, 28–36 (2013).
  • Veluthakkal R. & Dasgupta M. G. Isolation and characterization of pathogen defense-related class I chitinase from the actinorhizal tree Casuarina equisetifolia. For. Path. 42, 467–480 (2012).
  • Tapia G., Morales-Quintana L., Inostroza L. & Acuna H. Molecular characterisation of Ltchi7, a gene encoding a Class III endochitinase induced by drought stress in Lotus spp. Plant Biol. 13, 69–77 (2011). [PubMed]
  • Dempsey D. A. & Klessig D. F. SOS - too many signals for systemic acquired resistance? Trends Plant Sci. 17, 538–545 (2012). [PubMed]
  • Bari R. & Jones J. D. Role of plant hormones in plant defence responses. Plant Mol. Bio. 69, 473–488 (2009). [PubMed]
  • Sun L. et al. . Cotton cytochrome P450 CYP82D regulates systemic cell death by modulating the octadecanoid pathway. Nat. Commun. 5, 5372 (2014). [PMC free article] [PubMed]
  • Yang C. L. et al. . Cotton major latex protein 28 functions as a positive regulator of the ethylene responsive factor 6 in defense against Verticillium dahliae. Mol. Plant 8, 399–411 (2015). [PubMed]
  • Li C. et al. . Cotton WRKY1 mediates the plant defense-to-development transition during infection of cotton by Verticillium dahliae by activating JASMONATE ZIM-DOMAIN1 expression. Plant Physiol. 166, 2179–2194 (2014). [PubMed]
  • Finn R. D. et al. . Pfam: the protein families database. Nucleic Acids Res. 42, D222–230 (2014). [PMC free article] [PubMed]
  • Eddy S. R. Accelerated profile HMM searches. PLoS Comput. Biol. 7, e1002195 (2011). [PMC free article] [PubMed]
  • Wang F. X. et al. . Proteomic analysis of the sea-island cotton roots infected by wilt pathogen Verticillium dahliae. Proteomics 11, 4296–4309 (2011). [PubMed]
  • Jiang J. X. & Zhang T. Z. Extraction of total RNA in cotton tissues with CTAB-acidic phenolic method. Cotton Sci. 15, 166–167 (2003).
  • Wang X. Y., Lv K., Cai C. P., Xu J. & Guo W. Z. Establishment and application of TRV-mediated virus-induced gene silencing in cotton. Acta Agronom. Sin. 40, 1356–1363 (2014).

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group