Next Article in Journal
Circulating miRNA Signaling for Fatty Acid Metabolism in Response to a Maximum Endurance Test in Elite Long-Distance Runners
Previous Article in Journal
Carboxypeptidase Inhibitor LXN Expression in Endometrial Tissue Is Menstrual Cycle Phase-Dependent and Is Upregulated in Endometriotic Lesions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Identification of a Maize Chitinase Gene Family and the Induction of Its Expression by Fusarium verticillioides (Sacc.) Nirenberg (1976) Infection

by
Jesús Eduardo Cazares-Álvarez
1,
Paúl Alán Báez-Astorga
2,
Analilia Arroyo-Becerra
3 and
Ignacio Eduardo Maldonado-Mendoza
1,*
1
Departamento de Biotecnología Agrícola, Centro Interdisciplinario de Investigación para el Desarrollo Integral Regional (CIIDIR), Unidad Sinaloa, Instituto Politécnico Nacional, Guasave 81049, Sinaloa, Mexico
2
CONAHCYT—Departamento de Biotecnología Agrícola, Centro Interdisciplinario de Investigación para el Desarrollo Integral Regional (CIIDIR), Unidad Sinaloa, Instituto Politécnico Nacional, Guasave 81049, Sinaloa, Mexico
3
Laboratorio de Genómica Funcional y Biotecnología de Plantas, Centro de Investigación en Biotecnología Aplicada, Instituto Politécnico Nacional, Ex-Hacienda San Juan Molino Carretera Estatal Km 1.5, Santa Inés-Tecuexcomac-Tepetitla 90700, Tlaxcala, Mexico
*
Author to whom correspondence should be addressed.
Genes 2024, 15(8), 1087; https://doi.org/10.3390/genes15081087
Submission received: 24 July 2024 / Revised: 9 August 2024 / Accepted: 13 August 2024 / Published: 17 August 2024
(This article belongs to the Section Bioinformatics)

Abstract

:
Maize chitinases are involved in chitin hydrolysis. Chitinases are distributed across various organisms including animals, plants, and fungi and are grouped into different glycosyl hydrolase families and classes, depending on protein structure. However, many chitinase functions and their interactions with other plant proteins remain unknown. The economic importance of maize (Zea mays L.) makes it relevant for studying the function of plant chitinases and their biological roles. This work aims to identify chitinase genes in the maize genome to study their gene structure, family/class classification, cis-related elements, and gene expression under biotic stress, such as Fusarium verticillioides infection. Thirty-nine chitinase genes were identified and found to be distributed in three glycosyl hydrolase (GH) families (18, 19 and 20). Likewise, the conserved domains and motifs were identified in each GH family member. The identified cis-regulatory elements are involved in plant development, hormone response, defense, and abiotic stress response. Chitinase protein-interaction network analysis predicted that they interact mainly with cell wall proteins. qRT-PCR analysis confirmed in silico data showing that ten different maize chitinase genes are induced in the presence of F. verticillioides, and that they could have several roles in pathogen infection depending on chitinase structure and cell wall localization.

1. Introduction

Chitin is the second most abundant polysaccharide in nature and can be found in crustacean exoskeletons and fungal cell walls [1], in insects and animals [2]. A linear polymer of β-1,4-N-acetylglucosamine (GlcNAC), chitin occurs in three forms: α, β and γ. The most abundant form, α-chitin, can be found in the exoskeletons of insects and crustaceans; this form exhibits adjacent sheets along the c axis (crystallographic orientation) that are arranged in the same direction [3]. β-chitin has adjacent sheets along the c axis in the same direction, resulting in the parallel orientation of sheets [4]. β-chitin is found in squid pens [5] and insects [6]. γ-chitin has three chitin chains with alternating parallel and antiparallel aligned polymer chains, making this type of chitin the most complicated in structure [7].
Chitinases (EC 3.2.1.14) hydrolyze chitin glycosidic bonds, producing low-molecular-weight compounds. Chitinase genes are required for growth and development in chitin-synthesizing species. Chitinases are classified according to their amino acid sequences and belong to the GH18 (glycosyl hydrolase-18), GH19 and GH20 families, according to previously proposed criteria [8] as well as the CAZy database. Chitinases can be grouped into several classes: those that are within the GH19 family belong to classes I, II, IV, VI and VII, while GH18 members belong to classes III and V. The GH20 family consists of β-N-acetyl-D-hexosaminidases, which hydrolyze both GlcNAC and GalNAC [9].
Chitinases can also be classified according to the presence of specific modules (motifs) or conserved domains such as a pre-sequence, a catalytic domain, a linker domain, or the presence of a carbohydrate-binding domain [10]. Plant chitinases play an important role in plant defense against fungi containing chitin in their cell walls [11]. Chitinases belonging to classes I and IV possess a chitin-binding domain. In plants, this domain plays an important role in generating the chitin-derived signaling molecules (chito-oligosaccharides) that are typically released through binding of these enzymes to chitin in fungal cell walls, which is subsequently followed by chitin hydrolysis [12]. Chito-oligosaccharides (CHOs) are one of the main elicitors in plants that induce a series of defense mechanisms leading to the production of pathogenesis-related proteins, reactive oxygen species (ROS) and phytohormones [13,14,15].
Maize is an economically important crop worldwide, and it grows in a range of environments in many countries including Mexico, from where it originated. Maize is mainly used as human food, livestock feed, and in the production of bioethanol for energy [16]. This crop is susceptible to various fungal diseases, including Fusarium stalk, ear and root rot (SERR), a disease caused by F. verticillioides [17]. Different species of Fusarium (F. verticillioides, F. nygamai and F. thapsinum) are reported to coexist in mixed infections in a single maize plant [18]. Phytopathogenic fungi can induce the production of plant antifungal compounds including enzymes such as chitinases. In maize, 27 putative genes encoding endochitinases, including members of the GH18 (Class III) and GH19 families (classes I, II and IV), are induced by Trichoderma harzianum [19], Ustilago maydis [20] and Aspergillus flavus [21] colonization.
F. verticillioides is an aggressive phytopathogen that has evolved a mechanism to modify maize chitinases (ChitA and ChitB from the GH19 family) through the production of an effector chitinase-modifying protein (Cmp) [22]. This modification cuts the chitin-binding domain in the N-terminus upstream of a conserved polyglycine-rich site. Therefore, maize chitinases lacking the chitin-binding domain, although still active, cannot bind to the chitin substrate in the fungal cell wall. Cmp production has been confirmed in other fungal phytopathogens such as Bipolaris zeicola [23], Stenocarpella maydis [24] and Fusarium oxysporum [25]. The effect of Cmp activity is to limit fungal cell-wall chitin degradation and the production of the chito-oligomer elicitors, resulting in the reduction in plant defenses to allow these pathogens to infect maize cells.
Few studies have performed a genome-wide analysis of maize chitinases [19,21], and this could be explained by the fact that we know very little about the role of chitinases and their involvement in biotic and abiotic stresses, including maize–pathogen interactions. This is especially true in root tissues, such as during F. verticillioides infection. The availability of maize transcriptome data now makes it feasible to assess gene expression under different stress conditions. In the present study, we performed a genome-wide analysis to identify maize chitinase genes in order to characterize them according to their classification in GH families and classes. We examined the cis-elements present in their upstream promoter sequences, the presence of conserved motifs and domains in their encoded proteins, and their possible interaction with other maize proteins using co-expression network analysis. We also used transcriptome data to understand the putative involvement in the F. verticillioides-maize root interaction. Finally, this work provides insight and contributes to our comprehension of the genome-wide regulation of maize chitinase expression.

2. Materials and Methods

2.1. Identification of Putative Chitinase Gene

Nucleotide sequences of maize chitinase genes and their corresponding amino acid sequences were obtained from the Maize Genetics and Genomics Database (MaizeGDB) (https://www.maizegdb.org, accessed on 24 May 2023), using the B73 genome and “chitinase, chitin, endochitinase, and exochitinase” as discovery terms [26]. A total of 202 gene models were retrieved. From these, we removed all genes that were not annotated inside the Zm-B73-REFERENCE-NAM-5.0 genome database, resulting in 39 gene models. Thirty-nine amino acid sequences were analyzed using the Pfam database (https://www.ebi.ac.uk/interpro/) to confirm the presence of the glycoside hydrolase domain. Chromosomic localization, exon number, and transcript length were obtained from the Ensembl Plants database (https://www.plants.ensembl.org/). The molecular weight (kDa) and isoelectric point of chitinase amino acid sequences were obtained using the “Compute pI/MW” tool from the ExPAsy database (https://www.expasy.org/tools). The physical locations of chitinase genes were schematically drawn on their respective maize chromosomes using the online server MG2C (https://mg2c.iask.in/mg2c_v2.0/).

2.2. Gene Structure, Multiple Alignment, and Phylogenetic Analysis

Gene structure (exon–intron) was analyzed based on the coding sequence (CDS) of the chitinases at the Gene Structure Display Server (GSDS, https://gsds.gao-lab.org/) [27]. Evolutionary relationships among maize chitinase families were analyzed by performing an unrooted Neighbor-Joining phylogenetic tree constructed with the MEGA 11.0 software based on the Jones–Taylor–Thornton (JTT) model, using the MUSCLE alignment function with 1000 bootstrap replicates and default parameters for each chitinase family. The chitinase family tree was constructed for each chitinase class using orthologous genes from Sorghum bicolor, Oryza sativa and Triticum aestivum.

2.3. Prediction of cis-Acting Elements, Conserved Domains and Transcription Factor Binding Sites

To predict the cis-acting regulatory elements in the promoter region of the 39 Z. mays chitinase genes, a 2.0-kb sequence upstream of the translation initiation site (ATG) was analyzed, along with the 5′UTR region, using the PlantCARE website (http://bioinformatics.psb.ugent.be/webtools/plantcare/html) [28]. All amino acid sequences were analyzed using the NCBI-CDD database (https://www.ncbi.nlm.nih.gov/Structure/cdd/cdd.shtml) to find conserved domains, and the MEME website (https://meme-suite.org/index.html) to find conserved motifs [29]. The domains illustration was made using the Domain Illustrator software 2.0 version [30]. Subcellular localization prediction of each chitinase protein was obtained using WoLF PSORT (https://wolfpsort.hgc.jp/) [31], Plant-mSubP (https://bioinfo.usu.edu/Plant-mSubP/) [32] and MULocDeep (https://www.mu-loc.org/) [33], and signal peptide prediction was made by the SignalP 6.0 tool (https://services.healthtech.dtu.dk/services/SignalP-6.0/) [34]. To identify the binding sites of transcription factors in the promoter region of chitinases (2.0 kb upstream), the PlantRegMap database was used (http://plantregmap.gao-lab.org/) [35] with the following parameters: p-value ≤ 1 × 10−4 and a maximum number of 15.

2.4. In Silico Expression Analysis of Chitinase Genes Using RNA-Seq Data

In silico gene transcription analysis of the maize glycoside hydrolase family was performed using RNA-seq data in order to retrieve the expression data in response to abiotic stress. The fragment per kilobase of exon per million fragments mapped (FPKM) was obtained from the B73 V5 genome through a gene ID (Zm00001) search using qTeller in the MaizeGDB database [36]. Previous reports were used to study gene expression for abiotic stress using seedling [37] and leaf [38] samples. The extracted data were then log2 transformed and used to generate heatmaps via the TBtools package [39].

2.5. Quantitative Real-Time PCR Analysis

Previously, RNA-seq data were generated by our group from RNA samples of a rolled paper-towel assay using roots of 7-day-old maize seedlings (Garañón hybrid from Asgrow) inoculated with F. verticillioides P03 [40]. Differentially expressed genes (DEGs) were defined as those genes with adjusted p-value < 0.01 and Log2 fold change ±1. The differentially expressed genes (DEGs) obtained in the RNA-seq were analyzed and a list of chitinase DEGs with gene expression presented as fold-change value was generated. To confirm the induction of the maize root chitinase genes that respond to the Fv infection in the RNA-seq, gene expression analysis of the chitinase DEG genes was performed by RT-qPCR, using the SYBR Green Master Mix (Qiagen, Cat. No. 204074, Hilden, Germany) in a Rotor Gene-Q real-time PCR system instrument (Qiagen, Cat. No. 9001550, Hilden, Germany). The RNA used for RT-qPCR was the same as that used for RNA-seq. The RT-qPCR reactions were performed using two biological replicates with three technical replicates, and were run using 5 μL of SYBR Green Master Mix (2×), 10 μM of each primer and 10 ng of cDNA in a 10 μL final volume. The PCR program included a pre-heating step at 95 °C (5 min), followed by 40 cycles of 95 °C (30 s), 60 °C (30 s) and 72 °C (30 s). Relative quantification of chitinase genes was normalized to the cyclin-dependent kinase (CDK) gene [41], and the comparative threshold cycle method 2−∆∆Ct [42] was used to calculate the fold-change (FC) values in gene expression. The set of primers used in this work are listed in Supplementary Table S1.

2.6. Protein–Protein Interaction Network

To analyze protein–protein interactions, the STRING database (http://string-db.org) [43] was used to generate a co-expression network. All maize chitinase proteins were submitted to the STRING database. The minimum required score was set to medium confidence (0.400). The maximum number of interactors showed no more than 20 on the first shell, and no more than 10 on the second shell. The protein–protein network was obtained using text mining and experiments, as well as database and co-expression interaction sources.

3. Results

3.1. Genome-Wide Identification of Chitinase Genes in Z. mays

A total of 39 chitinase genes were obtained from the Zm-B73-REFERENCE-NAM-5.0 genome. Gene name and ID, chromosomal position, exon and transcript count, glycoside hydrolase family, molecular weight, and isoelectric point are presented in Table 1. The amino acid sequence length of the maize chitinase genes ranged from 96 (Zm00001eb270440) to 599 (Zm00001eb008880) amino acids. We found 17 chitinase genes belonging to the GH18 family (PF00704), 18 genes belonging to the GH19 family (PF00182), and 4 genes belonging to the GH20 (PF00728) family. Their predicted molecular weight (MW) ranged from 10.55 kDa (Zm00001eb270440) to 66.12 kDa (Zm00001eb008880), and their predicted isoelectric points ranged from 4.06 (Zm00001eb157820) to 10.35 (Zm00001eb272090). Most GH19 chitinases presented a basic isoelectric point, specifically Zm00001eb078730 and Zm00001eb425600 (Cta1 and Ctb1) from class IV, as well as nine GH18 chitinases and one from the GH20 family. An acidic isoelectric point was exhibited by GH19 family members of classes I, II, and IV; GH18 family members of class III; and three GH20 family members.

3.2. Phylogenetic Analysis and Multiple Sequence Alignment

An unrooted phylogenic tree was constructed to study the evolutionary relationships among the chitinase families in the maize genome. The chitinase protein phylogenetic tree was constructed using proteins from S. bicolor (36 sequences), T. aestivum (37 sequences) and O. sativa (32 sequences) (see Figure 1). The class classification revealed that chitinase members of the GH18 family belong to class III (containing hevamine, narbonin-like and SI-CLP protein sequences) and class V, and members of the GH19 family belong to classes I, II and IV. Although the Zm00001eb317090 gene belongs to family 19 and is grouped in the class IV clade, it does not contain a chitin-binding domain, and shares more protein characteristics that make it closer to class II proteins.

3.3. Chromosomal Location and Intron–Exon Architecture

All 39 chitinase genes were assigned to nine out of ten maize chromosomes (Figure 2). Chromosome number 6 contained the highest number of chitinase genes (seven), followed by chromosomes 3 and 8 (six chitinase genes in each chromosome). No chitinase genes were found in chromosome 9.
To study the diversity of maize chitinase genes, we next analyzed the exon–intron sequences. Among the thirty-nine chitinase genes in Z. mays, only four genes had more than two introns; four genes had two introns; eleven genes had one intron; and twenty genes did not have any introns. Zm000011eb266300 was the longest gene, containing 15 introns, followed by Zm000011eb365840 and Zm000011eb288150, with 14 and 13 introns, respectively, and Zm000011eb358410, with 7 introns (Figure 3).
Alternative splicing can occur in maize chitinase genes due to the presence of introns; Zm00001eb358410 (four splice variants), Zm00001eb002620 (three splice variants), Zm00001eb266300 (three splice variants), and Zm00001eb365840 (six splice variants) all presented more than one splicing transcript (Table S2). Since most chitinase genes had only one identified transcript, only one transcript (T001 from Table S2) from each gene was used to analyze each glycosyl hydrolase family.

3.4. Conserved Domains and Motifs

In order to better understand the structural differences between members in the chitinase gene families, we next identified the conserved domains and motifs (Figure 4). Motifs 6, 7, 8, 10 and 17 were exclusive to the GH18 gene family members, whereas motifs 1, 2, 3, 4, 5, 9, 13 and 16 were mainly found in the GH19 gene family. Members of the GH20 gene family contained conserved motifs 12, 14, 15, 18, 19 and 20 (Figure 4A, Table S3). As shown in Figure 4B, all protein sequences possessed domains typical of the glycosyl hydrolase families. According to their enzymatic activity, three different domains were present in members of class III in the GH18 gene family: hevamine, narbonin, and SI-CLP; other GH18 members had a domain typical of class V. However, GH18 family members did not contain any chitin-binding domain, which was only found in GH19 family members (ten out of eighteen). All GH20 family members presented the glycosyl hydrolase-20 and glycosyl hexaminidase HexB-like domains.
Additionally, almost all GH19 chitinase genes that presented chitin-binding domains were predicted to be extracellular proteins. However, Zm00001eb272050 did not contain a predicted signal peptide, but it had a chitin-binding domain and is predicted to be a cytoplasmic protein, and Zm00001eb078720 contained one signal peptide, and is predicted to be chloroplastic (Table S4).

3.5. Cis-Regulatory Elements and Transcription Factors in the Promoter Region

Cis-regulatory element analysis was performed to elucidate the possible regulatory mechanisms/pathways/signals of maize chitinase genes when elicited in response to hormones, abiotic stress, pathogen infections and plant development (Table S5). The hormone and abiotic stress cis-acting regulatory elements were widely present in chitinase genes. Defense-related elements had fewer repetitions among genes, and only a W-box element (involved in WRKY transcription factor recognition) was repeatedly found. The hormone-related elements ABRE (abscisic acid response), CGTCA-motif (methyl-jasmonate response), MYC (jasmonate signaling), as-1 and TGACG-motif (response to pathogens and methyl-jasmonate response) were the most abundant. We also found the abiotic stress cis-elements ARE (anaerobic induction), G-box (response to light) and STRE (general stress) in abundance. The AAGAA-motif (secondary xylem development) and CCGTCC-box (meristem activation) were the most frequently detected cis-regulatory elements related to plant developmental processes in chitinase upstream sequences. The methyl jasmonate response cis-elements were present in most chitinase promoters, followed by abiotic-stress and light-response elements. We also analyzed their gene expression (Supplementary Figure S1) and found that some chitinase genes may be transcriptionally induced under different stress conditions (cold, heat, salinity, UV light, and drought) (Table S6). The TF binding-site prediction revealed a total of 6281 binding sites for 18 TFs (Table S7) including AP2/ERF [44,45], B3 [46,47], BBR-BPC [48], bHLH [49], bZIP [50,51,52], C2H2 [53,54,55], Dof [56,57,58], G2-like [59,60,61], GATA [62], HD-ZIP [63,64,65], LBD [66,67], MIKC_MADS [68,69], MYB and MYB-related [70,71], NAC [72,73], TALE [74,75], TCP [76,77] and WRKY [78,79,80].

3.6. In Silico Analysis of Maize Chitinase Gene Expression in 7-Day-Old F. verticillioide Infected Roots

RNA-seq data from maize roots infected for 7 days with F. verticillioides (Fv) [40] were used to determine chitinase gene-expression patterns in maize roots (Table S8). Overall, RNA-seq data showed that 15 out of 39 chitinase genes were differentially expressed in maize roots in response to seven-day infection with Fv. All fifteen chitinase genes were upregulated. During RT-qPCR analysis, it was only possible to amplify 14 out of 15 differentially expressed chitinase genes, due to the sequence similarity among chitinase families. We were not able to amplify the Zm00001eb272050 gene by PCR, and therefore this gene was not included in the RT-qPCR analysis. Table 2 shows the RT-qPCR expression analysis of the 14 chitinase genes from maize roots identified in RNA-seq data under the Fv infection condition.
Ten out of the fourteen chitinase genes analyzed by RT-qPCR showed the same induced expression trend as observed in the RNA-seq data (Table 2). Three out of six GH18 gene family members, six out of seven GH19 members, and one GH20 member were induced by seven days of Fv infection.

3.7. Protein–Protein Interaction Network Analysis

Protein–protein interaction networks with the ten RT-qPCR validated chitinase genes from the maize–Fv interaction at seven days’ infection predicted 20 functional maize protein partners, exhibiting 100 interactions (Figure 5; Table S9). Zm00001eb317090 showed 30 interactions, followed by Zm00001eb008880 with 13 interactions and Zm00001eb167340 with 10 interactions. GH18 chitinases (Zm00001eb167340 and Zm00001eb168350) were predicted to interact with β-hexosaminidases, while GH19 (Zm00001eb272090) was the main maize chitinase partner to interact with cellulose synthase proteins, and GH20 (Zm00001eb008880) was related to β-galactosidase proteins (Table S9). These results show diverse patterns of maize chitinase co-expression with maize cell wall proteins.

4. Discussion

In plants, glycosyl hydrolases are present as large gene families that are involved in several biological and defense mechanisms, and they can be expressed in different cell compartments [81]. The GH18 and GH19 chitinase families are some of the most studied plant enzymes, due to their involvement in different metabolic processes, in addition to their synergistic function in plant development and defense responses [82]. GH18 chitinases have been reported in fungi [83], bacteria [84] and plants [85]. GH19 chitinases are predominantly found in plants (e.g., class I chitinases are found in Camellia sinensis [86], class II in soybean [87] and class IV in grape [88]), while GH18 class III chitinases have been reported in Rhododendron irroratum and Pteris ryukyuensis [89,90]. On the other hand, GH20 chitinase family members are β-hexosaminidases; these enzymes have been detected before maize seed germination and in scutellum in seedlings [91]. β-Hexosaminidases are proposed to participate in N-glycan processing, depending on their subcellular localization: N-glycan trimming occurs in the vacuole for vacuolar proteins, whereas N-glycans are processed on secretory glycoproteins for plasma membrane proteins [92].
Several maize chitinase genes described here have more than two introns, and some of these genes also presented alternative splicing (Table S2). Alternative splicing occurs when there is a deviation from constitutive splicing, in which the removal of certain introns is omitted, producing various forms of mature mRNA [93]. This process is involved in mediating diverse biological processes [94] and has an important role in plant abiotic-stress responses [95], as well as a direct effect on protein function, amino acid sequence, and enzymatic properties [96].
Protein domains are short amino acid sequences that can be considered as a functional site, or as structural and evolutionary units [97]. These units are considered important for protein function [98], and proteins can acquire several domains over time through gene fusion and exon recombination [99]. However, protein domains can suffer changes like mutations or deletions, and cause dramatic effects on gene function [100]. Plant chitinases usually have different domains such as catalytic, chitin-binding, and linker domains [10]. In the present work, all maize chitinases had at least one domain belonging to a glycosyl hydrolase family. GH18 sequences presented only one type of GH-conserved domain per gene, and no chitin-binding-related domains. Gene members of the GH20 family possess two domains: the glycosyl hydrolase-20 and the glycosyl hexaminidase HexB-like domain. Most GH19 chitinases presented either a glycosyl hydrolase-19 domain and/or a chitin-binding domain; the latter is able to bind and hydrolyze insoluble or soluble chitin forms, such as the GH19 chitinase from Hevea brasiliensis [101] or Streptomyces griseus. The chitin-binding domain of S. griseus differs slightly from those of plants, and it can also bind to cellulose [102]. For example, inserting the chitin-binding domain from Serratia marcescens into the structure of a chimeric chitinase from Trichoderma atroviride (Chit42) resulted in increased enzymatic activity and colloidal chitin binding, as well as higher antifungal activity against phytopathogenic fungi [103]. This demonstrates the importance of the chitin-binding domain during chitin degradation.
In plants, the chitin-binding domain is important due to its relatedness to mechanisms of avoidance of fungal phytopathogen infection, such as in the interactions between Fusarium solani and Nicotiana tabacum [104], Candida albicans and Moringa oleifera [105], and Trichoderma viride and Brassica juncea [106]. These studies suggest that chitinases with a chitin-binding domain could potentially be used as antifungal enzymes, since they can bind to chitin from the fungal cell wall and hydrolyze it to produce small chito-oligomers that elicit the plant defense [107]. These may occur in GH19 class I and IV chitinases, due to the presence of a chitin-binding domain [108].
Motifs are short amino acid sequences shared by protein family members that have a specific structural function [109]. In our study, maize chitinases presented the NYNG conserved motif (Motif 1, Table S3). This motif is generally located within a loop III structure close to the amino acid catalytic triad of the GH19 family (Glu221, Arg361 and Glu234) [110] that belongs to the lysozyme superfamily. GH19 maize family members, like most of the GH22, 23, 24 and 46 families, have two other main elements: an α-helix and a β-hairpin motif [111], which are important, since they contain glutamate residues acting as general acid catalysts that allow proteins to perform chitin catalysis. Thirty-two chitinases have signal peptides that localize them for their secretion. Likewise, the hevamine-type chitinase from H. brasiliensis also contains a signal peptide sequence at the N-terminus [112] for its vacuolar secretion. In maize, according to the WoLF PSORT (an extension of the PSORTII program) in silico subcellular prediction software, one chitinase is targeted to the vacuole (Zm00001eb167340), nine are targeted to the chloroplast, and eighteen are extracellular proteins (Table S4). GH18 proteins are predicted to have a different subcellular distribution as compared to GH19 proteins, where most chitinases are predicted to be extracellular. Plant chitinases are proposed to be targeted to the chloroplast, endoplasmic reticulum and vacuole, and these subcellular localizations have been tested in Arabidopsis thaliana [113].
Different cis-regulatory elements in maize chitinases have been found to be related to plant defense, development, response to abiotic stress, and hormones. The plant development elements AAGA-, CAT-box, CCGTCC, RY-element, HD-Zip1 and GCN-4 motifs were present in maize chitinase genes. The GCN-4 motif is fundamental to specific meristem gene expression [114]. Hormone-related elements such as ABRE, as-1, MYC, CGTCA-, and TGACG-motifs were also present in most chitinase genes. These cis-elements, such as ABRE (abscisic acid-response-related element), are related to several dehydration-response genes that trigger ABA production [115]. Likewise, Chi2 and Chi14 induction in Cucumis sativus L. was also observed in response to jasmonic acid, ethylene and salicylic acid treatments [116], suggesting that maize chitinases could be involved in hormone signaling in response to either abiotic or biotic stresses during development. ARE, G-box, STRE, Box4, and DRE-core were the most abundant cis-elements in abiotic stress-related categories such as general stress, as well as light and dehydration responses. These elements have also been reported in chitinase genes from Allium sativum L. [117], Thalassiosira pseudonana [26] and B. juncea [118]. The W-box in chitinase genes plays a role in the induction of plant defense genes against pathogenic fungi [119]. TC-rich repeat elements are also involved in the induction of defense genes, and they have been reported in C. sativus L. [120]. The GH18 and GH19 family members are induced under osmotic stress in Ammopiptanthus nanus [121] and Solanum lycopersicum [122]. According to our in silico maize chitinase expression analysis (Supplementary Figure S1), most GH18 genes are induced under cold, heat, UV, and salt stresses, such as a class III gene from A. thaliana that is induced in response to different environmental stress conditions [123]. Interestingly, most GH19 and GH20 genes are either repressed or do not show induction under different types of stress (drought and cold). This may be related to the activity of GH19 and GH20; GH19 genes are mainly induced in the presence of phytopathogenic organisms [124,125], while GH20 genes are involved in N-glycan processing [92].
Transcription factors (TFs) can be induced under abiotic or biotic stresses. TFs can be regulated by phosphorylation, inhibitory factors or de novo synthesis [126]. Maize chitinase genes exhibit target sites for several transcription factor families. These TFs are related to hormone response, plant development, and defense, including the ERE and WRKY TFs that are involved in the plant-pathogen defense response [127,128].
Previous studies have shown the induction of chitinases in maize plants infected with F. verticillioides and A. flavus [129,130], which increase the plant resistance to fungal pathogen infection. This effect may be related to the release of chito-oligosaccharides acting as pathogen-associated molecular patterns (PAMPs), which are derived from fungal cell wall degradation by chitinases [131]. Chitinases can accumulate in different plant organelles such as the chloroplast, vacuole and nucleus, or they can be secreted into the extracellular space [113]. The extracellular space, or apoplast, is considered to be a cellular compartment that plays important roles such as nutrient metabolism, plant growth, defense and signaling [132]. In the present study, we identified by RNA-seq 15 maize chitinase genes induced in response to a 7-day F. verticillioides infection; additionally, we confirmed by qRT-PCR the induction of ten genes. Five GH19 chitinase members (Zm00001eb317090, Zm00001eb272090, Zm00001eb346860, Zm00001eb078730 and Zm00001eb425600) are predicted to be extracellular proteins. Four of these genes are indicated in Figure 6: (a) Zm00001eb272090 [Chitinase 2], (b) Zm00001eb346860 [Chitinase 21], (c) Zm00001eb078730 [Chitinase A1], and (d) Zm00001eb425600 [Chitinase B1]). These genes possess a chitin-binding domain and motif 5, which are important for substrate recognition and chitin binding. ChitA and ChitB were previously reported to be targeted in a polyglycine-rich region by F. verticillioides chitinase modifying enzymes (Cmps) (Figure 6) [23]. This suggests that (a) Zm00001eb272090 and (b) Zm00001eb346860 may also be targets of Cmps, since they also possess a potential polyglycine site targeted by this effector protein. Taken together, these results suggest that these four extracellular chitinases could be targets of F. verticillioides Cmps, preventing them from (1) hydrolyzing the chitin present in the fungal cell wall, and (2) inducing the presence of oligosaccharides (PAMPs) that trigger the plant defense mechanisms. As a result, the maize root does not perceive the invasion by the fungus, allowing for the F. verticillioides infection to proceed. Zm00001eb317090 (Brittle stalk4; [g] in Figure 6) does not possess this Cmp region, and our co-expression network analysis predicted its interaction with other cell wall proteins (Table S9) that may be involved in cell wall strengthening [133], including two cytoplasmic/membrane GH18 members (Zm00001eb250900 [chitinase 28] and Zm00001eb168350 [chitinase 29]) (e and f in Figure 6, respectively).
The other induced chitinases may also participate in the response to pathogen attack, playing different roles in the vacuole/apoplast ([j] Zm00001eb167340, chitinase 27 [pI = 8.69]). In addition, other intracellularly localized induced chitinases ([i] Zm00001eb354540 [chitinase 23] and [h] Zm00001eb008880 [Exochitinase 2]; chloroplastic/extracellular) may play other functions in this particular interaction, such as BC15/OsCTL1 (chitinase-like 1), which is necessary for cellulose biosynthesis in rice [134]. Chloroplastic chitinases may also be regulated by the accumulation of ROS in response to fungus infection [135].
Co-expression network analysis showed interactions with other maize proteins according to text mining, experiments, databases, and co-expression. The most prominent interactions were between β-hexosaminidase and cellulose synthase, which play a major role in cell wall formation [136]. β-Hexosaminidase (HEXO) has been detected in different plant tissues and has a role during seed germination. HEXOs also show different types of catalytic activity and substrate affinity to N-glycans [137], and they are capable of hydrolyzing N-acetyl-hexosamine residues to produce low-molecular-weight compounds [138]. However, their physiological and biological functions are not fully understood, and it has been hypothesized that they are mainly involved in chitin degradation. On the other hand, cellulose synthase (generally found as an enzymatic complex) is responsible for cellulose synthesis at the plasma membrane [139]. Zm00001eb317090 (GH19, Brittle stalk4 [g]) is related to Zm00001eb250900 (GH18, Chitinase 28 [e]) and Zm00001eb168350 (GH18, Chitinase 29 [f]), which are predicted to be cytoplasmic/membrane proteins of class III narbonin chitinases. In addition, Zm00001eb317090 (g) possesses the GhCTL (chitinase-like proteins) consensus [140], which differs from the classical GH19 consensus; this difference in catalytic residues may allow this protein to participate in the synthesis of primary and secondary cell walls. However, although chitinases are known to bind to N-acetyl-D-glucosamine (GlcNAC) residues, plants do not produce chitin. Therefore, there must be a mechanism by which chitinases can be induced in plants, even if they are not under pathogen attack. For example, one previous study demonstrated that N-acetyl-D-glucosamine-containing arabinogalactan proteins are sensitive to endochitinases, and that they could be implicated in embryogenesis [141]. Since arabinogalactan proteins are typically found in membranes or cell walls, these may be some of the few proteins that contain GlcNAC, and they could explain why chitinases participate in cell wall biosynthesis

5. Conclusions

In summary, we analyzed maize chitinase genes in order to learn about their putative function and possible role in fungal infection. Our findings suggest that they may participate in different plant biological pathways, even though we addressed their involvement during fungal pathogen infection (e.g., F. verticillioides). We also found by qRT-PCR that maize chitinases can participate in Fv infection, according to their domain type and subcellular localization. This works opens new avenues to study novel biochemical mechanisms in which chitinases may play a major role.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes15081087/s1, Figure S1: Heat map of chitinase mRNA expression in Z. mays. (a) Maize chitinase gene expression in seedling samples; (b) maize chitinase gene expression in leaf samples. The color range indicates a gene expression gradient from repressed (blue) to induced (red) expression, relative to the untreated-plant control; Table S1: Chitinase primers used in this study; Table S2: Number of splice variants obtained from Ensembl Plants database; Table S3: Motif consensus; Table S4: Subcellular localization of maize chitinase; Table S5: cis-regulatory element in 2.0 kb upstream sequence of maize chitinase genes; Table S6: Abiotic-stress transcriptome data obtained from MaizeGDB database used in this study; Table S7: Transcription-factor site identification in maize chitinase-promoter regions; Table S8: Transcriptome data obtained from Baez-Astorga (2022) from maize-F. verticillioides interaction; Table S9: Functional protein-association-network analysis using STRING database.

Author Contributions

Conceptualization, J.E.C.-Á.; methodology, J.E.C.-Á.; validation, J.E.C.-Á.; formal analysis, J.E.C.-Á.; investigation, J.E.C.-Á.; resources, I.E.M.-M.; writing—original draft preparation, J.E.C.-Á.; writing—review and editing, P.A.B.-A., A.A.-B. and I.E.M.-M.; visualization, J.E.C.-Á. and P.A.B.-A.; supervision, A.A.-B. and I.E.M.-M.; project administration, I.E.M.-M.; funding acquisition, I.E.M.-M. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by CONAHCYT Fronteras de la Ciencia [grant number 2016-01-2510] and Instituto Politécnico Nacional (Grants SIP 20211507, SIP 20221357 and SIP 20232027). Jesús Eduardo Cazares Álvarez was granted a PhD fellowship [grant number 934735] from CONAHCYT and complementary fellowship support from IPN (BEIFI fellowship programme), and Paúl Alán Báez Astorga was granted a postdoctoral fellowship [grant number 70883] from CONAHCYT.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding author.

Acknowledgments

We thank Brandon Loveall for English proofreading of the manuscript.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Elieh-Ali-Komi, D.; Hamblin, M.R. Chitin and chitosan: Production and application of versatile biomedical nanomaterials. Am. Math. Mon. 2004, 111, 915. [Google Scholar]
  2. Rinaudo, M. Chitin and chitosan: Properties and applications. Prog. Polym. Sci. 2006, 31, 603–632. [Google Scholar] [CrossRef]
  3. Sahraee, S.; Milani, J.M. Chapter 8—Chitin and chitosan-based blends, composites, and nanocomposites for packaging applications. In Handbook of Chitin and Chitosan; Gopi, S., Thomas, S., Pius, A., Eds.; Elsevier: Amsterdam, The Netherlands, 2020; pp. 247–271. [Google Scholar]
  4. Aranaz, I.; Mengibar, M.; Harris, R.; Panos, I.; Miralles, B.; Acosta, N.; Galed, G.; Heras, A. Functional characterization of chitin and chitosan. Curr. Chem. Biol. 2009, 3, 203–230. [Google Scholar] [CrossRef]
  5. Lavall, R.L.; Assis, O.B.G.; Campana-Filho, S.P. β-Chitin from the pens of Loligo Sp.: Extraction and characterization. Bioresour. Technol. 2007, 98, 2465–2472. [Google Scholar] [CrossRef]
  6. Eisemann, C.H.; Binnington, K.C. The peritrophic membrane: Its formation, structure, chemical composition and permeability in relation to vaccination against ectoparasitic arthropods. Int. J. Parasitol. 1994, 24, 15–26. [Google Scholar] [CrossRef]
  7. Blackwell, J. Physical methods for the determination of chitin structure and conformation. Methods Enzymol. 1988, 161, 435–442. [Google Scholar] [CrossRef]
  8. Hamid, R.; Khan, M.A.; Ahmad, M.; Ahmad, M.M.; Abdin, M.Z.; Musarrat, J.; Javed, S. Chitinases: An update. J. Pharm. Bioallied Sci. 2013, 5, 21–29. [Google Scholar] [CrossRef] [PubMed]
  9. Henrissat, B.; Bairoch, A. New families in the classification of glycosyl hydrolases based on amino acid sequence similarities. Biochem. J. 1993, 293, 781–788. [Google Scholar] [CrossRef] [PubMed]
  10. Volpicella, M.; Leoni, C.; Fanizza, I.; Placido, A.; Pastorello, E.A.; Ceci, L.R. Overview of plant chitinases identified as food allergens. J. Agric. Food Chem. 2014, 62, 5734–5742. [Google Scholar] [CrossRef]
  11. Di Rosa, M.; Distefano, G.; Zorena, K.; Malaguarnera, L. Chitinases and immunity: Ancestral molecules with new functions. Immunobiology 2015, 221, 399–411. [Google Scholar] [CrossRef] [PubMed]
  12. Aboitiz, N.; Vila-Perelló, M.; Groves, P.; Asensio, J.L.; Andreu, D.; Cañada, F.J.; Jiménez-Barbero, J. NMR and modeling studies of protein-carbohydrate interactions: Synthesis, three-dimensional structure, and recognition properties of a minimum hevein domain with binding affinity for chitooligosaccharides. ChemBioChem. 2004, 5, 1245–1255. [Google Scholar] [CrossRef] [PubMed]
  13. Ali, S.; Ganai, B.A.; Kamili, A.N.; Bhat, A.A.; Mir, Z.A.; Bhat, J.A.; Tyagi, A.; Islam, S.T.; Mushtaq, M.; Yadav, P.; et al. Pathogenesis-related proteins and peptides as promising tools for engineering plants with multiple stress tolerance. Microbiol. Res. 2018, 212, 29–37. [Google Scholar] [CrossRef]
  14. Kumar, M.; Brar, A.; Yadav, M.; Chawade, A.; Vivekanand, V.; Pareek, N. Chitinases—Potential candidates for enhanced plant resistance towards fungal pathogens. Agriculture 2018, 8, 88. [Google Scholar] [CrossRef]
  15. Oyeleye, A.; Normi, Y.M. Chitinase: Diversity, limitations, and trends in engineering for suitable applications. Biosci. Rep. 2018, 38, BSR2018032300. [Google Scholar] [CrossRef] [PubMed]
  16. Ranum, P.; Peña-Rosas, J.P.; Garcia-Casal, M.N. Global maize production, utilization, and consumption. Ann. N. Y. Acad. Sci. 2014, 1312, 105–112. [Google Scholar] [CrossRef]
  17. Figueroa-Rivera, M.G.; Rodríguez-Guerra, R.; Guerrero-Aguilar, B.Z.; González-Chavira, M.M.; Pons-Hernández, J.L.; Jiménez-Bremont, J.F.; Ramírez-Pimentel, J.G.; Andrio-Enríquez, E.; Mendoza-Elos, M. Caracterización de especies de Fusarium asociadas a la pudrición de raíz de maíz en Guanajuato, México. Rev. Mex. Fitopatol. 2010, 28, 124–134. [Google Scholar]
  18. Leyva-Madrigal, K.Y.; Larralde-Corona, C.P.; Apodaca-Sánchez, M.A.; Quiroz-Figueroa, F.R.; Mexia-Bolaños, P.A.; Portillo-Valenzuela, S.; Ordaz-Ochoa, J.; Maldonado-Mendoza, I.E. Fusarium species from the Fusarium fujikuroi species complex involved in mixed infections of maize in northern Sinaloa, México. J. Phytopathol. 2015, 163, 486–497. [Google Scholar] [CrossRef]
  19. Shoresh, M.; Harman, G.E. Genome-wide identification, expression and chromosomal location of the genes encoding chitinolytic enzymes in Zea mays. Mol. Genet. Genom. 2008, 280, 173–185. [Google Scholar] [CrossRef]
  20. Ökmen, B.; Kemmerich, B.; Hilbig, D.; Wemhöner, R.; Aschenbroich, J.; Perrar, A.; Huesgen, P.F.; Schipper, K.; Doehlemann, G. Dual function of a secreted fungalysin metalloprotease in Ustilago maydis. New Phytol. 2018, 220, 249–261. [Google Scholar] [CrossRef] [PubMed]
  21. Hawkins, L.K.; Mylroie, J.E.; Oliveira, D.A.; Smith, J.S.; Ozkan, S.; Windham, G.L.; Williams, W.P.; Warburton, M.L. Characterization of the maize chitinase genes and their effect on Aspergillus flavus and aflatoxin accumulation resistance. PLoS ONE 2015, 10, e0126185. [Google Scholar] [CrossRef]
  22. Naumann, T.A.; Wicklow, D.T.; Price, N.P.J. Identification of a chitinase-modifying protein from Fusarium verticillioides: Truncation of a host resistance protein by a fungalysin metalloprotease. J. Biol. Chem. 2011, 286, 35358–35366. [Google Scholar] [CrossRef] [PubMed]
  23. Naumann, T.A.; Wicklow, D.T.; Kendra, D.F. Maize seed chitinase is modified by a protein secreted by Bipolaris zeicola. Physiol. Mol. Plant Pathol. 2009, 74, 134–141. [Google Scholar] [CrossRef]
  24. Naumann, T.A.; Price, N.P.J. Purification and in vitro activities of a chitinase-modifying protein from the corn ear rot pathogen Stenocarpella maydis. Physiol. Mol. Plant Pathol. 2019, 106, 74–80. [Google Scholar] [CrossRef]
  25. Karimi Jashni, M.; Dols, I.H.M.; Iida, Y.; Boeren, S.; Beenen, H.G.; Mehrabi, R. Synergistic action of a metalloprotease and a serine protease from Fusarium oxysporum cleaves chitin-binding tomato chitinases. Mol. Plant Microbes Interact. 2015, 28, 996–1008. [Google Scholar] [CrossRef] [PubMed]
  26. Cheng, H.; Shao, Z.; Lu, C.; Duan, D. Genome-wide identification of chitinase genes in Thalassiosira pseudonana and analysis of their expression under abiotic stresses. BMC Plant Biol. 2021, 21, 87. [Google Scholar] [CrossRef] [PubMed]
  27. Hu, B.; Jin, J.; Guo, A.Y.; Zhang, H.; Luo, J.; Gao, G. GSDS 2.0: An upgraded gene feature visualization server. Bioinformatics 2015, 31, 1296–1297. [Google Scholar] [CrossRef]
  28. Lescot, M.; Déhais, P.; Thijs, G.; Marchal, K.; Moreau, Y.; Van De Peer, Y.; Rouzé, P.; Rombauts, S. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002, 30, 325–327. [Google Scholar] [CrossRef]
  29. Bailey, T.L.; Elkan, C. Fitting a mixture model by expectation maximization to discover motifs in biopolymers. Proc. Int. Conf. Intell. Syst. Mol. Biol. 1994, 2, 28–36. [Google Scholar]
  30. Ren, J.; Wen, L.; Gao, X.; Jin, C.; Xue, Y.; Yao, X. DOG 1.0: Illustrator of protein domain structures. Cell Res. 2009, 19, 271–273. [Google Scholar] [CrossRef]
  31. Horton, P.; Park, K.J.; Obayashi, T.; Fujita, N.; Harada, H.; Adams-Collier, C.J.; Nakai, K. WoLF PSORT: Protein localization predictor. Nucleic Acids Res. 2007, 35, 585–587. [Google Scholar] [CrossRef]
  32. Sahu, S.S.; Loaiza, C.D.; Kaundal, R. Plant-MSubP: A computational framework for the prediction of single- and multi-target protein subcellular localization using integrated machine-learning approaches. AoB Plants 2020, 12, plz068. [Google Scholar] [CrossRef]
  33. Jiang, Y.; Wang, D.; Yao, Y.; Eubel, H.; Künzler, P.; Møller, I.M.; Xu, D. MULocDeep: A Deep-learning framework for protein subcellular and suborganellar localization prediction with residue-level interpretation. Comput. Struct. Biotechnol. J. 2021, 19, 4825–4839. [Google Scholar] [CrossRef] [PubMed]
  34. Teufel, F.; Almagro Armenteros, J.J.; Johansen, A.R.; Gíslason, M.H.; Pihl, S.I.; Tsirigos, K.D.; Winther, O.; Brunak, S.; von Heijne, G.; Nielsen, H. SignalP 6.0 predicts all five types of signal peptides using protein language models. Nat. Biotechnol. 2022, 40, 1023–1025. [Google Scholar] [CrossRef]
  35. Tian, F.; Yang, D.-C.; Meng, Y.-Q.; Jin, J.; Gao, G. PlantRegMap: Charting functional regulatory maps in plants. Nucleic Acids Res. 2019, 48, D1104–D1113. [Google Scholar] [CrossRef]
  36. Woodhouse, M.R.; Cannon, E.K.; Portwood, J.L.; Harper, L.C.; Gardiner, J.M.; Schaeffer, M.L.; Andorf, C.M. A pan-genomic approach to genome databases using maize as a model system. BMC Plant Biol. 2021, 21, 385. [Google Scholar] [CrossRef]
  37. Makarevitch, I.; Waters, A.J.; West, P.T.; Stitzer, M.; Hirsch, C.N.; Ross-Ibarra, J.; Springer, N.M. Transposable elements contribute to activation of maize genes in response to abiotic stress. PLoS Genet. 2015, 11, e1004915. [Google Scholar] [CrossRef]
  38. Forestan, C.; Aiese Cigliano, R.; Farinati, S.; Lunardon, A.; Sanseverino, W.; Varotto, S. Stress-induced and epigenetic-mediated maize transcriptome regulation study by means of transcriptome reannotation and differential expression analysis. Sci. Rep. 2016, 6, 30446. [Google Scholar] [CrossRef]
  39. Chen, C.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.; Xia, R. TBtools: An integrative toolkit developed for interactive analyses of big biological data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef]
  40. Báez-Astorga, P.A. Caracterización de los Mecanismos de Antagonismo de Bacillus Cereus B25 Contra Fusarium verticillioides y Análisis Transcriptómico de Raíces de Maíz Inoculadas con B25 y Fv. Ph.D. Thesis, Instituto Politécnico Nacional, CIIDIR-Sinaloa, Guasave, Sinaloa, Mexico, 2022. [Google Scholar]
  41. Lin, F.; Jiang, L.; Liu, Y.; Lv, Y.; Dai, H.; Zhao, H. Genome-wide identification of housekeeping genes in maize. Plant Mol. Biol. 2014, 86, 543–554. [Google Scholar] [CrossRef] [PubMed]
  42. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-delta delta c(t)) method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  43. Szklarczyk, D.; Franceschini, A.; Wyder, S.; Forslund, K.; Heller, D.; Huerta-Cepas, J.; Simonovic, M.; Roth, A.; Santos, A.; Tsafou, K.P.; et al. STRING V10: Protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015, 43, D447–D452. [Google Scholar] [CrossRef] [PubMed]
  44. Zang, Z.; Wang, Z.; Zhao, F.; Yang, W.; Ci, J.; Ren, X.; Jiang, L.; Yang, W. Maize ethylene response factor ZmERF061 is required for resistance to Exserohilum Turcicum. Front. Plant Sci. 2021, 12, 630413. [Google Scholar] [CrossRef] [PubMed]
  45. Zhu, Y.; Liu, Y.; Zhou, K.; Tian, C.; Aslam, M.; Zhang, B.; Liu, W.; Zou, H. Overexpression of ZmEREBP60 enhances drought tolerance in maize. J. Plant Physiol. 2022, 275, 153763. [Google Scholar] [CrossRef] [PubMed]
  46. Grant-Grant, S.; Schaffhauser, M.; Baeza-Gonzalez, P.; Gao, F.; Conéjéro, G.; Vidal, E.A.; Gaymard, F.; Dubos, C.; Curie, C.; Roschzttardtz, H. B3 transcription factors determine iron distribution and FERRITIN gene expression in embryo but do not control total seed iron content. Front. Plant Sci. 2022, 13, 870078. [Google Scholar] [CrossRef] [PubMed]
  47. Ahmad, B.; Zhang, S.; Yao, J.; Ur Rahman, M.; Hanif, M.; Zhu, Y.; Wang, X. Genomic organization of the B3-domain transcription factor family in grapevine (Vitis Vinifera L.) and expression during seed development in seedless and seeded cultivars. Int. J. Mol. Sci. 2019, 20, 4553. [Google Scholar] [CrossRef] [PubMed]
  48. Simonini, S.; Kater, M.M. Class I basic pentacysteine factors regulate homeobox genes involved in meristem size maintenance. J. Exp. Bot. 2014, 65, 1455–1465. [Google Scholar] [CrossRef]
  49. Zhang, T.; Lv, W.; Zhang, H.; Ma, L.; Li, P.; Ge, L.; Li, G. Genome-wide analysis of the Basic Helix-Loop-Helix (BHLH) transcription factor family in maize. BMC Plant Biol. 2018, 18, 235. [Google Scholar] [CrossRef] [PubMed]
  50. Wei, K.; Chen, J.; Wang, Y.; Chen, Y.; Chen, S.; Lin, Y.; Pan, S.; Zhong, X.; Xie, D. Genome-wide analysis of BZIP-encoding genes in maize. DNA Res. 2012, 19, 463–476. [Google Scholar] [CrossRef] [PubMed]
  51. Ying, S.; Zhang, D.F.; Fu, J.; Shi, Y.S.; Song, Y.C.; Wang, T.Y.; Li, Y. Cloning and characterization of a maize BZIP transcription factor, ZmbZIP72, confers drought and salt tolerance in transgenic Arabidopsis. Planta 2012, 235, 253–266. [Google Scholar] [CrossRef]
  52. Pautler, M.; Eveland, A.L.; Larue, T.; Yang, F.; Weeks, R.; Lunde, C.; Il Je, B.; Meeley, R.; Komatsu, M.; Vollbrecht, E.; et al. FASCIATED EAR4 encodes a BZIP transcription factor that regulates shoot meristem size in maize. Plant Cell 2015, 27, 104–120. [Google Scholar] [CrossRef]
  53. Liu, Q.; Wang, Z.; Xu, X.; Zhang, H.; Li, C. Genome-wide analysis of C2H2 Zinc-finger family transcription factors and their responses to abiotic stresses in Poplar (Populus Trichocarpa). PLoS ONE 2015, 10, e0134753. [Google Scholar] [CrossRef]
  54. Jiao, Z.; Wang, L.; Du, H.; Wang, Y.; Wang, W.; Liu, J.; Huang, J.; Huang, W.; Ge, L. Genome-wide study of C2H2 Zinc Finger gene family in Medicago Truncatula. BMC Plant Biol. 2020, 20, 401. [Google Scholar] [CrossRef] [PubMed]
  55. Ghosh, S.; Kant, R.; Pradhan, A.; Jha, G. RS_CRZ1, a C2H2-Type transcription factor is required for pathogenesis of Rhizoctonia Solani AG1-IA in tomato. Mol. Plant-Microbe Interact. 2021, 34, 26–38. [Google Scholar] [CrossRef] [PubMed]
  56. Chen, Y.; Cao, J. Comparative analysis of Dof transcription factor family in maize. Plant Mol. Biol. Rep. 2015, 33, 1245–1258. [Google Scholar] [CrossRef]
  57. Ruta, V.; Longo, C.; Lepri, A.; De Angelis, V.; Occhigrossi, S.; Costantino, P.; Vittorioso, P. The DOF transcription factors in seed and seedling development. Plants 2020, 9, 218. [Google Scholar] [CrossRef]
  58. Wu, J.; Chen, L.; Chen, M.; Zhou, W.; Dong, Q.; Jiang, H.; Cheng, B. The DOF-domain transcription factor ZmDOF36 positively regulates starch synthesis in transgenic maize. Front. Plant Sci. 2019, 10, 465. [Google Scholar] [CrossRef]
  59. Liu, F.; Xu, Y.; Han, G.; Zhou, L.; Ali, A.; Zhu, S.; Li, X. molecular evolution and genetic variation of G2-like transcription factor genes in maize. PLoS ONE 2016, 11, e0161763. [Google Scholar] [CrossRef]
  60. Qin, M.; Zhang, B.; Gu, G.; Yuan, J.; Yang, X.; Yang, J.; Xie, X. Genome-wide analysis of the G2-like transcription factor genes and their expression in different senescence stages of tobacco (Nicotiana Tabacum L.). Front. Genet. 2021, 12, 626352. [Google Scholar] [CrossRef] [PubMed]
  61. Murmu, J.; Wilton, M.; Allard, G.; Pandeya, R.; Desveaux, D.; Singh, J.; Subramaniam, R. Arabidopsis GOLDEN2-LIKE (GLK) transcription factors activate jasmonic acid (ja)-dependent disease susceptibility to the biotrophic pathogen Hyaloperonospora arabidopsidis, as well as ja-independent plant immunity against the necrotrophic pathogen Botrytis cinerea. Mol. Plant Pathol. 2014, 15, 174–184. [Google Scholar] [CrossRef]
  62. Behringer, C.; Schwechheimer, C. B-GATA transcription factors—Insights into their structure, regulation, and role in plant development. Front. Plant Sci. 2015, 6, 90. [Google Scholar] [CrossRef]
  63. Turchi, L.; Baima, S.; Morelli, G.; Ruberti, I. Interplay of HD-Zip II and III transcription factors in auxin-regulated plant development. J. Exp. Bot. 2015, 66, 5043–5053. [Google Scholar] [CrossRef] [PubMed]
  64. Shao, J.; Haider, I.; Xiong, L.; Zhu, X.; Hussain, R.M.F.; Övernäs, E.; Meijer, A.H.; Zhang, G.; Wang, M.; Bouwmeester, H.J.; et al. Functional analysis of the HD-Zip transcription factor genes Oshox12 and Oshox14 in rice. PLoS ONE 2018, 13, e0199248. [Google Scholar] [CrossRef] [PubMed]
  65. Tang, Y.; Bao, X.; Wang, S.; Liu, Y.; Tan, J.; Yang, M.; Zhang, M.; Dai, R.; Yu, X. A physic nut stress-responsive HD-Zip transcription factor, JcHDZ07, confers enhanced sensitivity to salinity stress in transgenic Arabidopsis. Front. Plant Sci. 2019, 10, 1–12. [Google Scholar] [CrossRef] [PubMed]
  66. Jia, R.; Li, C.; Wang, Y.; Qin, X.; Meng, L.; Sun, X. Genome-wide analysis of LBD transcription factor genes in Dendrobium catenatum. Int. J. Mol. Sci. 2022, 23, 2089. [Google Scholar] [CrossRef] [PubMed]
  67. Zhang, X.; He, Y.; He, W.; Su, H.; Wang, Y.; Hong, G.; Xu, P. Structural and functional insights into the LBD family involved in abiotic stress and flavonoid synthases in Camellia sinensis. Sci. Rep. 2019, 9, 15651. [Google Scholar] [CrossRef] [PubMed]
  68. Teo, Z.W.N.; Zhou, W.; Shen, L. Dissecting the function of MADS-Box transcription factors in orchid reproductive development. Front. Plant Sci. 2019, 10, 1474. [Google Scholar] [CrossRef] [PubMed]
  69. Liu, G.; Li, F.; Shi, G.; Wang, L.; Wang, L.; Fan, L. Identification of MADS-Box transcription factors in Iris Laevigata and functional assessment of IlSEP3 and IlSVP during flowering. Int. J. Mol. Sci. 2022, 23, 9950. [Google Scholar] [CrossRef] [PubMed]
  70. Dubos, C.; Stracke, R.; Grotewold, E.; Weisshaar, B.; Martin, C.; Lepiniec, L. MYB transcription factors in Arabidopsis. Trends Plant Sci. 2010, 15, 573–581. [Google Scholar] [CrossRef] [PubMed]
  71. Ambawat, S.; Sharma, P.; Yadav, N.R.; Yadav, R.C. MYB transcription factor genes as regulators for plant responses: An overview. Physiol. Mol. Biol. Plants 2013, 19, 307–321. [Google Scholar] [CrossRef]
  72. Wang, G.; Yuan, Z.; Zhang, P.; Liu, Z.; Wang, T.; Wei, L. Genome-wide analysis of NAC transcription factor family in maize under drought stress and rewatering. Physiol. Mol. Biol. Plants 2020, 26, 705–717. [Google Scholar] [CrossRef]
  73. Lu, M.; Ying, S.; Zhang, D.F.; Shi, Y.S.; Song, Y.C.; Wang, T.Y.; Li, Y. A maize stress-responsive NAC transcription factor, zmsnac1, confers enhanced tolerance to dehydration in transgenic Arabidopsis. Plant Cell Rep. 2012, 31, 1701–1711. [Google Scholar] [CrossRef]
  74. Peng, W.; Yang, Y.; Xu, J.; Peng, E.; Dai, S.; Dai, L.; Wang, Y.; Yi, T.; Wang, B.; Li, D.; et al. TALE transcription factors in sweet orange (Citrus Sinensis): Genome-wide identification, characterization, and expression in response to biotic and abiotic stresses. Front. Plant Sci. 2022, 12, 814252. [Google Scholar] [CrossRef]
  75. Wang, J.; Zhao, P.; Cheng, B.; Zhang, Y.; Shen, Y.; Wang, X.; Zhang, Q.; Lou, Q.; Zhang, S.; Wang, B.; et al. Identification of TALE transcription factor family and expression patterns related to fruit chloroplast development in tomato (Solanum lycopersicum L.). Int. J. Mol. Sci. 2022, 23, 4507. [Google Scholar] [CrossRef]
  76. Li, S. The Arabidopsis thaliana TCP transcription factors: A broadening horizon beyond development. Plant Signal. Behav. 2015, 10, e1044192. [Google Scholar] [CrossRef] [PubMed]
  77. Danisman, S. TCP transcription factors at the interface between environmental challenges and the plant’s growth responses. Front. Plant Sci. 2016, 7, 1930. [Google Scholar] [CrossRef] [PubMed]
  78. Li, H.; Gao, Y.; Xu, H.; Dai, Y.; Deng, D.; Chen, J. ZmWRKY33, a WRKY maize transcription factor conferring enhanced salt stress tolerances in Arabidopsis. Plant Growth Regul. 2013, 70, 207–216. [Google Scholar] [CrossRef]
  79. Cai, R.; Dai, W.; Zhang, C.; Wang, Y.; Wu, M.; Zhao, Y.; Ma, Q.; Xiang, Y.; Cheng, B. The maize WRKY transcription factor ZMWRKY17 negatively regulates salt stress tolerance in transgenic Arabidopsis plants. Planta 2017, 246, 1215–1231. [Google Scholar] [CrossRef] [PubMed]
  80. Wani, S.H.; Anand, S.; Singh, B.; Bohra, A.; Joshi, R. WRKY transcription factors and plant defense responses: Latest discoveries and prospects. Plant Cell Rep. 2021, 40, 1071–1085. [Google Scholar] [CrossRef]
  81. Minic, Z. Physiological roles of plant glycoside hydrolases. Planta 2008, 227, 723–740. [Google Scholar] [CrossRef]
  82. Rathore, A.S.; Gupta, R.D. Chitinases from bacteria to human: Properties, applications, and future perspectives. Enzyme Res. 2015, 2015, 791907. [Google Scholar] [CrossRef]
  83. Xia, G.; Jin, C.; Zhou, J.; Yang, S.; Zhang, S.; Jin, C. A novel chitinase having a unique mode of action from Aspergillus fumigatus YJ-407. Eur. J. Biochem. 2001, 268, 4079–4085. [Google Scholar] [CrossRef] [PubMed]
  84. Yang, S.; Fu, X.; Yan, Q.; Guo, Y.; Liu, Z.; Jiang, Z. Cloning, expression, purification and application of a novel chitinase from a thermophilic marine bacterium Paenibacillus barengoltzii. Food Chem. 2016, 192, 1041–1048. [Google Scholar] [CrossRef]
  85. Tanaka, J.; Takashima, T.; Abe, N.; Fukamizo, T.; Numata, T.; Ohnuma, T. Characterization of two rice GH18 chitinases belonging to family 8 of plant pathogenesis-related proteins. Plant Sci. 2023, 326, 111524. [Google Scholar] [CrossRef]
  86. Bordoloi, K.S.; Krishnatreya, D.B.; Baruah, P.M.; Borah, A.K.; Mondal, T.K.; Agarwala, N. Genome-wide identification and expression profiling of chitinase genes in tea Camellia sinensis (L.) O. Kuntze) under Biotic Stress Conditions. Physiol. Mol. Biol. Plants 2021, 27, 369–385. [Google Scholar] [CrossRef] [PubMed]
  87. Lv, P.; Zhang, C.; Xie, P.; Yang, X.; El-Sheikh, M.A.; Hefft, D.I.; Ahmad, P.; Zhao, T.; Bhat, J.A. Genome-wide identification and expression analyses of the chitinase gene family in response to white mold and drought stress in soybean (Glycine Max). Life 2022, 12, 1340. [Google Scholar] [CrossRef]
  88. Vincenzi, S.; Bierma, J.; Wickramasekara, S.I.; Curioni, A.; Gazzola, D.; Bakalinsky, A.T. Characterization of a grape class IV chitinase. J. Agric. Food Chem. 2014, 62, 5660–5668. [Google Scholar] [CrossRef] [PubMed]
  89. Zha, H.G.; Milne, R.I.; Zhou, H.X.; Chen, X.Y.; Sun, H. Identification and cloning of class ii and iii chitinases from alkaline floral nectar of Rhododendron irroratum, ericaceae. Planta 2016, 244, 805–818. [Google Scholar] [CrossRef]
  90. Kitaoku, Y.; Umemoto, N.; Ohnuma, T.; Numata, T.; Taira, T.; Sakuda, S.; Fukamizo, T. A Class III chitinase without disulfide bonds from the fern, Pteris ryukyuensis: Crystal structure and ligand-binding studies. Planta 2015, 242, 895–907. [Google Scholar] [CrossRef]
  91. Oikawa, A.; Itoh, E.; Ishihara, A.; Iwamura, H. Purification and Characterization of β-N-acetylhexosaminidase from maize seedlings. J. Plant Physiol. 2003, 160, 991–999. [Google Scholar] [CrossRef]
  92. Strasser, R.; Bondili, J.S.; Schoberer, J.; Svoboda, B.; Liebminger, E.; Glössl, J.; Altmann, F.; Steinkellner, H.; Mach, L. Enzymatic properties and subcellular localization of Arabidopsis β-N-acetylhexosaminidases. Plant Physiol. 2007, 145, 5–16. [Google Scholar] [CrossRef]
  93. Wang, Y.; Liu, J.; Huang, B.; Xu, Y.-M.; Li, J.; Huang, L.-F.; Lin, J.; Zhang, J.; Min, Q.-H.; Yang, W.-M.; et al. Mechanism of alternative splicing and its regulation. Biomed. Rep. 2015, 3, 152–158. [Google Scholar] [CrossRef] [PubMed]
  94. Irimia, M.; Penny, D.; Roy, S.W. Coevolution of genomic intron number and splice sites. Trends Genet. 2007, 23, 321–325. [Google Scholar] [CrossRef]
  95. Gao, B.; Chen, M.; Oliver, M.J. Alternative splicing: From abiotic stress tolerance to evolutionary genomics. Int. J. Mol. Sci. 2023, 24, 6708. [Google Scholar] [CrossRef]
  96. Kelemen, O.; Convertini, P.; Zhang, Z.; Wen, Y.; Shen, M.; Falaleeva, M.; Stamm, S. Function of alternative splicing. Gene 2013, 514, 1–30. [Google Scholar] [CrossRef]
  97. Wang, Y.; Zhang, H.; Zhong, H.; Xue, Z. Protein domain identification methods and online resources. Comput. Struct. Biotechnol. J. 2021, 19, 1145–1153. [Google Scholar] [CrossRef]
  98. Holm, L.; Sander, C. Parser for protein folding units. Proteins Struct. Funct. Genet. 1994, 268, 256–268. [Google Scholar] [CrossRef]
  99. Toll-riera, M.; Albà, M.M. Emergence of novel domains in proteins. BMC Evol. Biol. 2013, 13, 47. [Google Scholar] [CrossRef] [PubMed]
  100. Buljan, M.; Bateman, A. the evolution of protein domain families. Biochem. Soc. Trans. 2009, 37, 751–755. [Google Scholar] [CrossRef] [PubMed]
  101. Martínez-Caballero, S.; Cano-Sánchez, P.; Mares-Mejía, I.; Díaz-Sánchez, A.G.; Macías-Rubalcava, M.L.; Hermoso, J.A.; Rodríguez-Romero, A. Comparative study of two GH19 chitinase-like proteins from Hevea brasiliensis, one exhibiting a novel carbohydrate-binding domain. FEBS J. 2014, 281, 4535–4554. [Google Scholar] [CrossRef]
  102. Itoh, Y.; Kawase, T.; Nikaidou, N.; Fukada, H.; Mitsutomi, M.; Watanabe, T.; Itoh, Y. functional analysis of the chitin-binding domain of a family 19 chitinase from Streptomyces griseus HUT6037: Substrate-binding affinity and cis-dominant increase of antifungal function. Biosci. Biotechnol. Biochem. 2002, 66, 1084–1092. [Google Scholar] [CrossRef]
  103. Matroodi, S.; Motallebi, M.; Zamani, M.; Moradyar, M. Designing a new chitinase with more chitin binding and antifungal activity. World J. Microbiol. Biotechnol. 2013, 29, 1517–1523. [Google Scholar] [CrossRef]
  104. Badrhadad, A.; Nazarian, F.; Ahmad, F. Fusion of a chitin-binding domain to an antibacterial peptide to enhance resistance to Fusarium Solani in tobacco (Nicotiana tabacum). 3 Biotech 2018, 8, 391. [Google Scholar] [CrossRef]
  105. Neto, J.X.S.; Pereira, M.L.; Oliveira, J.T.A.; Rocha-Bezerra, L.C.B.; Lopes, T.D.P.; Costa, H.P.S.; Sousa, D.O.B.; Rocha, B.A.M.; Grangeiro, T.B.; Freire, J.E.C.; et al. A chitin-binding protein purified from Moringa oleifera seeds presents anticandidal activity by increasing cell membrane permeability and reactive oxygen species production. Front. Microbiol. 2017, 8, 980. [Google Scholar] [CrossRef]
  106. Fung, K.-L.; Zhao, K.-J.; He, Z.-M.; Chye, M.-L. Tobacco-expressed Brassica juncea chitinase BjCHI1 shows antifungal activity in vitro. Plant Mol. Biol. 2002, 50, 283–294. [Google Scholar] [CrossRef]
  107. Vaghela, B.; Vashi, R.; Rajput, K.; Joshi, R. Plant chitinases and their role in plant defense: A comprehensive review. Enzyme Microb. Technol. 2022, 159, 110055. [Google Scholar] [CrossRef]
  108. Taira, T.; Mahoe, Y.; Kawamoto, N.; Onaga, S.; Iwasaki, H.; Ohnuma, T.; Fukamizo, T. Cloning and characterization of a small family 19 chitinase from moss (Bryum coronatum). Glycobiology 2011, 21, 644–654. [Google Scholar] [CrossRef]
  109. Parras-Moltó, M.; Campos-Laborie, F.J.; García-Diéguez, J.; Rodríguez-Griñolo, M.R.; Pérez-Pulido, A.J. Classification of protein motifs based on subcellular localization uncovers evolutionary relationships at both sequence and functional levels. BMC Bioinform. 2013, 14, 229. [Google Scholar] [CrossRef] [PubMed]
  110. Ubhayasekera, W. Structure and function of chitinases from glycoside hydrolase family 19. Polym. Int. 2011, 60, 890–896. [Google Scholar] [CrossRef]
  111. Wohlkönig, A.; Huet, J.; Looze, Y.; Wintjens, R. Structural relationships in the lysozyme superfamily: Significant evidence for glycoside hydrolase signature motifs. PLoS ONE 2010, 5, e15388. [Google Scholar] [CrossRef]
  112. Bokma, E.; Spiering, M.; Chow, K.S.; Mulder, P.P.M.F.A.; Subroto, T.; Beintema, J.J. Determination of CDNA and genomic DNA sequences of hevamine, a chitinase from the rubber tree Hevea brasiliensis. Plant Physiol. Biochem. 2001, 39, 367–376. [Google Scholar] [CrossRef]
  113. Xu, F.; Fan, C.; He, Y. Chitinases in Oryza sativa ssp. Japonica and Arabidopsis thaliana. J. Genet. Genom. 2007, 34, 138–150. [Google Scholar] [CrossRef]
  114. Wu, C.Y.; Suzuki, A.; Washida, H.; Takaiwa, F. The GCN4 motif in a rice glutelin gene is essential for endosperm-specific gene expression and is activated by Opaque-2 in transgenic rice plants. Plant J. 1998, 14, 673–683. [Google Scholar] [CrossRef]
  115. Narusaka, Y.; Nakashima, K.; Shinwari, Z.K.; Sakuma, Y.; Furihata, T.; Abe, H.; Narusaka, M.; Shinozaki, K.; Yamaguchi-Shinozaki, K. Interaction between two cis-acting elements, ABRE and DRE, in ABA-dependent expression of Arabidopsis Rd29A gene in response to dehydration and high-salinity stresses. Plant J. 2003, 34, 137–148. [Google Scholar] [CrossRef]
  116. Xu, J.; Zhang, N.; Wang, K.; Xian, Q.; Dong, J.; Qi, X.; Chen, X. Chitinase Chi 2 positively regulates cucumber resistance against Fusarium oxysporum f. sp. cucumerinum. Genes 2022, 13, 62. [Google Scholar] [CrossRef]
  117. Filyushin, M.A.; Anisimova, O.K.; Kochieva, E.Z.; Shchennikova, A.V. Genome-wide identification and expression of chitinase class I genes in garlic (Allium sativum L.) cultivars resistant and susceptible to Fusarium proliferatum. Plants 2021, 10, 720. [Google Scholar] [CrossRef]
  118. Mir, Z.A.; Ali, S.; Shivaraj, S.M.; Bhat, J.A.; Singh, A.; Yadav, P.; Rawat, S.; Paplao, P.K.; Grover, A. Genome-wide identification and characterization of chitinase gene family in Brassica juncea and Camelina sativa in response to Alternaria brassicae. Genomics 2020, 112, 749–763. [Google Scholar] [CrossRef]
  119. Gao, Y.; Jia, S.; Wang, C.; Wang, F.; Wang, F.; Zhao, K. BjMYB1, a transcription factor implicated in plant defence through activating BjCHI1 chitinase expression by binding to a W-box-like element. J. Exp. Bot. 2016, 67, 4647–4658. [Google Scholar] [CrossRef]
  120. Bartholomew, E.S.; Black, K.; Feng, Z.; Liu, W.; Shan, N.; Zhang, X.; Wu, L.; Bailey, L.; Zhu, N.; Qi, C.; et al. Comprehensive analysis of the chitinase gene family in cucumber (Cucumis sativus L.): From gene identification and evolution to expression in response to Fusarium oxysporum. Int. J. Mol. Sci. 2019, 20, 5309. [Google Scholar] [CrossRef]
  121. Cao, S.; Wang, Y.; Li, Z.; Shi, W.; Gao, F.; Zhou, Y.; Zhang, G.; Feng, J. Genome-wide identification and expression analyses of the chitinases under cold and osmotic stress in Ammopiptanthus nanus. Genes 2019, 10, 472. [Google Scholar] [CrossRef]
  122. Cao, J.; Tan, X. Comprehensive analysis of the chitinase family genes in tomato (Solanum lycopersicum). Plants 2019, 8, 52. [Google Scholar] [CrossRef]
  123. Takenaka, Y.; Nakano, S.; Tamoi, M.; Sakuda, S.; Fukamizo, T. Chitinase gene expression in response to environmental stresses in Arabidopsis thaliana: Chitinase inhibitor allosamidin enhances stress tolerance. Biosci. Biotechnol. Biochem. 2009, 73, 1066–1071. [Google Scholar] [CrossRef]
  124. Sahai, A.S.; Manocha, M.S. Chitinases of fungi and plants: Their involvement in morphogenesis and host-parasite interaction. FEMS Microbiol. Rev. 1993, 11, 317–338. [Google Scholar] [CrossRef]
  125. Orlando, M.; Buchholz, P.C.F.; Lotti, M.; Pleiss, J. The GH19 engineering database: Sequence diversity, substrate scope, and evolution in glycoside hydrolase family 19. PLoS ONE 2021, 16, e0256817. [Google Scholar] [CrossRef]
  126. Armstrong, S.; Korcok, J.; Sims, S.M.; Dixon, S.J. Activation of transcription factors by extracellular nucleotides in immune and related cell types. Purinergic Signal. 2007, 3, 59–69. [Google Scholar] [CrossRef]
  127. Huo, T.; Wang, C.T.; Yu, T.F.; Wang, D.M.; Li, M.; Zhao, D.; Li, X.T.; Fu, J.D.; Xu, Z.S.; Song, X.Y. Overexpression of ZmWRKY65 transcription factor from maize confers stress resistances in transgenic Arabidopsis. Sci. Rep. 2021, 11, 4024. [Google Scholar] [CrossRef]
  128. Zang, Z.; Lv, Y.; Liu, S.; Yang, W.; Ci, J.; Ren, X.; Wang, Z.; Wu, H.; Ma, W.; Jiang, L.; et al. A novel ERF transcription factor, ZmERF105, positively regulates maize resistance to Exserohilum turcicum. Front. Plant Sci. 2020, 11, 850. [Google Scholar] [CrossRef]
  129. Lanubile, A.; Ferrarini, A.; Maschietto, V.; Delledonne, M.; Marocco, A.; Bellin, D. functional genomic analysis of constitutive and inducible defense responses to Fusarium verticillioides infection in maize genotypes with contrasting ear rot resistance. BMC Genom. 2014, 15, 710. [Google Scholar] [CrossRef]
  130. Shu, X.; Livingston, D.P.; Woloshuk, C.P.; Payne, G.A. Comparative histological and transcriptional analysis of maize kernels infected with Aspergillus flavus and Fusarium verticillioides. Front. Plant Sci. 2017, 8, 2075. [Google Scholar] [CrossRef]
  131. Bravo, J.M.; Campo, S.; Murillo, I.; Coca, M.; San Segundo, B. Fungus- and wound-induced accumulation of mRNA containing a class II chitinase of the pathogenesis-related protein 4 (PR-4) family of maize. Plant Mol. Biol. 2003, 52, 745–759. [Google Scholar] [CrossRef]
  132. Darino, M.; Kanyuka, K.; Hammond-Kosack, K.E. Apoplastic and vascular defenses. Essays Biochem. 2022, 66, 595–605. [Google Scholar] [CrossRef]
  133. Jiao, S.; Hazebroek, J.P.; Chamberlin, M.A.; Perkins, M.; Sandhu, A.S.; Gupta, R.; Simcox, K.D.; Yinghong, L.; Prall, A.; Heetland, L.; et al. Chitinase-like1 plays a role in stalk tensile strength in maize. Plant Physiol. 2019, 181, 1127–1147. [Google Scholar] [CrossRef]
  134. Wu, B.; Zhang, B.; Dai, Y.; Zhang, L.; Shang-Guan, K.; Peng, Y.; Zhou, Y.; Zhu, Z. Brittle culm15 encodes a membrane-associated chitinase-like protein required for cellulose biosynthesis in rice. Plant Physiol. 2012, 159, 1440–1452. [Google Scholar] [CrossRef]
  135. Rossi, F.R.; Krapp, A.R.; Bisaro, F.; Maiale, S.J.; Pieckenstain, F.L.; Carrillo, N. Reactive oxygen species generated in chloroplasts contribute to tobacco leaf infection by the necrotrophic fungus Botrytis cinerea. Plant J. 2017, 92, 761–773. [Google Scholar] [CrossRef]
  136. Lei, L.; Li, S.; Gu, Y. cellulose synthase complexes: Composition and regulation. Front. Plant Sci. 2012, 3, 75. [Google Scholar] [CrossRef]
  137. Alvisi, N.; van Noort, K.; Dwiani, S.; Geschiere, N.; Sukarta, O.; Varossieau, K.; Nguyen, D.-L.; Strasser, R.; Hokke, C.H.; Schots, A.; et al. β-hexosaminidases along the secretory pathway of Nicotiana benthamiana have distinct specificities toward engineered helminth n-glycans on recombinant glycoproteins. Front. Plant Sci. 2021, 12, 638454. [Google Scholar] [CrossRef]
  138. Hossain, M.A.; Roslan, H.A. Molecular phylogeny and predicted 3D structure of plant beta-D-N-Acetylhexosaminidase. Sci. World J. 2014, 2014, 186029. [Google Scholar] [CrossRef]
  139. Polko, J.K.; Kieber, J.J. The regulation of cellulose biosynthesis in plants. Plant Cell 2019, 31, 282–296. [Google Scholar] [CrossRef]
  140. Zhang, D.; Hrmova, M.; Wan, C.-H.; Wu, C.; Balzen, J.; Cai, W.; Wang, J.; Densmore, L.D.; Fincher, G.B.; Zhang, H.; et al. Members of a new group of chitinase-like genes are expressed preferentially in cotton cells with secondary walls. Plant Mol. Biol. 2004, 54, 353–372. [Google Scholar] [CrossRef]
  141. van Hengel, A.J.; Tadesse, Z.; Immerzeel, P.; Schols, H.; van Kammen, A.; de Vries, S.C. N-Acetylglucosamine and glucosamine-containing arabinogalactan proteins control somatic embryogenesis. Plant Physiol. 2001, 125, 1880–1890. [Google Scholar] [CrossRef]
Figure 1. Phylogenetic tree of Z. mays chitinase genes. GH18 includes family members of class III (hevamine-like domain; narbonin-like domain; SI-CLP-like domain); GH19 includes family members of class I, class II and class IV; and GH20 includes family members that are β-hexosaminidase-like proteins. The phylogenetic tree was built using the neighbor-joining (NJ) method and is based on the Jones–Taylor–Thornton (JTT) model, using the MUSCLE alignment function. The numbers at the tree nodes represent bootstrap percentage values. Bootstrap values are from 1000 replicates. Orthologous genes are labeled as Ta: T. aestivum; Sb: S. bicolor; Os: O. sativa.
Figure 1. Phylogenetic tree of Z. mays chitinase genes. GH18 includes family members of class III (hevamine-like domain; narbonin-like domain; SI-CLP-like domain); GH19 includes family members of class I, class II and class IV; and GH20 includes family members that are β-hexosaminidase-like proteins. The phylogenetic tree was built using the neighbor-joining (NJ) method and is based on the Jones–Taylor–Thornton (JTT) model, using the MUSCLE alignment function. The numbers at the tree nodes represent bootstrap percentage values. Bootstrap values are from 1000 replicates. Orthologous genes are labeled as Ta: T. aestivum; Sb: S. bicolor; Os: O. sativa.
Genes 15 01087 g001
Figure 2. Chitinase gene distribution in maize chromosomes. Gene names are indicated on the left and right sides of the chromosome. The Y-axis represents the distance (Mb) between genes. Their locations may be close in each chromosome, but they are not contiguous with each other.
Figure 2. Chitinase gene distribution in maize chromosomes. Gene names are indicated on the left and right sides of the chromosome. The Y-axis represents the distance (Mb) between genes. Their locations may be close in each chromosome, but they are not contiguous with each other.
Genes 15 01087 g002
Figure 3. The intron–exon structure of the chitinase-coding gene families in Z. mays. The blue bars indicate the 5′ and 3′ UTR untranslated regions, the yellow bars indicate coding sequences (CDS), and the gray lines indicate intron sequences.
Figure 3. The intron–exon structure of the chitinase-coding gene families in Z. mays. The blue bars indicate the 5′ and 3′ UTR untranslated regions, the yellow bars indicate coding sequences (CDS), and the gray lines indicate intron sequences.
Genes 15 01087 g003
Figure 4. Schematic structure of conserved maize chitinase motifs and domains. (A) Motifs. Colored boxes represent conserved motifs and gray lines represent non-conserved sequences; motif number sequences are shown in Table S3. (B) Domains. The conserved domains are represented by different colors.
Figure 4. Schematic structure of conserved maize chitinase motifs and domains. (A) Motifs. Colored boxes represent conserved motifs and gray lines represent non-conserved sequences; motif number sequences are shown in Table S3. (B) Domains. The conserved domains are represented by different colors.
Genes 15 01087 g004
Figure 5. Functional protein-association-network analysis of maize chitinases using the STRING database. Yellow lines: text mining network; black lines: co-expression; blue lines: from curated databases; violet lines: experimentally determined. Colored nodes represent the first shell of proteins interactors, and the white nodes represent the second shell of interactors.
Figure 5. Functional protein-association-network analysis of maize chitinases using the STRING database. Yellow lines: text mining network; black lines: co-expression; blue lines: from curated databases; violet lines: experimentally determined. Colored nodes represent the first shell of proteins interactors, and the white nodes represent the second shell of interactors.
Genes 15 01087 g005
Figure 6. A proposed model of maize chitinase induction in the presence of F. verticillioides. Four chitinases (a–d) possess the chitin-binding domain and Cmp site (fungalysin); these proteins may act synergistically to hydrolyze chitin in the fungal cell wall to produce chitin oligosaccharides (elicitors). These enzymes may also be the target of fungalysin, and their binding capacity could be eliminated; however, maize may produce other chitinases (i and j) that do not possess the chitin-binding domain/Cmp site to hydrolyze chitin. Due to the absence of a carbohydrate-binding domain, the chitin-binding capacity is low in these two chitinases, but their catalytic activity remains intact, and the plant could be inducing these proteins as an alternative strategy to degrade fungal chitin. On the other hand, proteins g and h, and e and f, may be induced in a coordinated mechanism to interact with cellulose synthases and to promote cell wall biosynthesis in response to Fv infection, in order to reinforce the plant cell wall. Different letters represent different chitinase proteins. Letters a–d: Zm00001eb272090, Zm00001eb346860, Zm00001eb078730, and Zm00001eb425600, respectively; letters e and f: Zm00001eb250900 and Zm00001eb168350; letters g and h: Zm00001eb317090, Zm00001eb008880; and letters i and j: Zm00001eb354540 and Zm00001eb167340.
Figure 6. A proposed model of maize chitinase induction in the presence of F. verticillioides. Four chitinases (a–d) possess the chitin-binding domain and Cmp site (fungalysin); these proteins may act synergistically to hydrolyze chitin in the fungal cell wall to produce chitin oligosaccharides (elicitors). These enzymes may also be the target of fungalysin, and their binding capacity could be eliminated; however, maize may produce other chitinases (i and j) that do not possess the chitin-binding domain/Cmp site to hydrolyze chitin. Due to the absence of a carbohydrate-binding domain, the chitin-binding capacity is low in these two chitinases, but their catalytic activity remains intact, and the plant could be inducing these proteins as an alternative strategy to degrade fungal chitin. On the other hand, proteins g and h, and e and f, may be induced in a coordinated mechanism to interact with cellulose synthases and to promote cell wall biosynthesis in response to Fv infection, in order to reinforce the plant cell wall. Different letters represent different chitinase proteins. Letters a–d: Zm00001eb272090, Zm00001eb346860, Zm00001eb078730, and Zm00001eb425600, respectively; letters e and f: Zm00001eb250900 and Zm00001eb168350; letters g and h: Zm00001eb317090, Zm00001eb008880; and letters i and j: Zm00001eb354540 and Zm00001eb167340.
Genes 15 01087 g006
Table 1. Information for the 39 Z. mays chitinase genes obtained from the Zm-B73-REFERENCE-NAM-5.0 genome.
Table 1. Information for the 39 Z. mays chitinase genes obtained from the Zm-B73-REFERENCE-NAM-5.0 genome.
Gene NameGene IDChr.Gene PositionTranscript Length (bp)Number of Amino AcidsGH FamilyIsoelectric Point (pI)
Chn1Zm00001eb1578203217635788–2176369891202295184.06
Chn7Zm00001eb301500710810239–108115841346312187.11
Chn12Zm00001eb1471203181373589–1813749701382311189.12
Chn13Zm00001eb1471403181400773–814021091269299185.92
Chn14Zm00001eb1471503181403132–1814043901155294184.93
Chn15Zm00001eb1471603181489354–181490403924307188.54
Chn16Zm00001eb301490710733581–107350521472307189.21
Chn17Zm00001eb0479401246418756–2464198321077308184.11
Chn19Zm00001eb42085010110156286–1101581621807484188.88
Chn25Zm00001eb2832806139958834–1399604201587368185.72
Chn26Zm00001eb2832606139888026–1398892021177315186.95
Chn27Zm00001eb16734047657537–76589881452301188.69
Chn28Zm00001eb2509005205495288–2054964931206286184.99
Chn29Zm00001eb168350413184755–131858561102286184.97
Chn30Zm00001d0414263120119356–120120333978325186.17
Chn31Zm00001eb3584108148327731–1483345251690430188.42
Chn34Zm00001eb169950422152266–221544111889399184.81
Bk4Zm00001eb3170907139523498–1395288521974328198.16
Chn2Zm00001eb272090692474902–9247606511642611910.35
Chn3Zm00001eb228510565159105–651603621258193199.46
Chn4Zm00001eb340820842175395–42176288807268198.09
Chn5Zm00001eb078720235057479–350587661160278197.86
Chn8Zm00001eb007850123078494–23079090597198195.92
Chn9Zm00001eb270440680869472–8086976229196199.1
Chn10Zm00001eb33235081145764–11469511188285196.97
Chn11Zm00001eb078740235087711–35095502926227197.57
Chn20Zm00001eb228500565072200–650733591160282195.1
Chn21Zm00001eb346860890582152–905835991350357198.04
Chn22Zm00001eb00262017344153–73488771087258199.1
Chn23Zm00001eb3545408132934929–1329361811253311195.97
Chn24Zm00001eb022500186493293–86494456940210194.79
Cta1Zm00001eb078730235084623–350859271202280198.44
Ctb1Zm00001eb42560010129884568–1298858881195281198.92
EPR4Zm00001eb2466405186238039–1862392561150271195.14
Prp10Zm00001eb272050692426037–924273051269379194.82
Exo1Zm00001eb266300642318401–423319032511545205.78
Exo2Zm00001eb008880127253118–272558162619599206.47
Exo3Zm00001eb2881506158303978–1583089282332529205.94
Exo4Zm00001eb3658408170849713–1708553992161525205.36
Table 2. Chitinase gene-expression data obtained by RNA-seq and validated by qRT-PCR.
Table 2. Chitinase gene-expression data obtained by RNA-seq and validated by qRT-PCR.
Gene IDSynonymsFC (Fold Change)
RNA-SeqRT-qPCR
GH18
Zm00001eb157820Chitinase chem 14.921.64 a
Zm00001eb047940Chitinase 175.140.70 a
Zm00001eb167340Chitinase 2713.413.61
Zm00001eb250900Chitinase 284.173.76
Zm00001eb168350Chitinase 292.713.49
Zm00001eb169950Chitinase 342.571.81 a
GH19
Zm00001eb317090Brittle stalk42.312.75
Zm00001eb078740Chitinase 119.580.67 a
Zm00001eb272090Chitinase 25.162.35
Zm00001eb346860Chitinase 217.272.51
Zm00001eb354540Chitinase 236.694.99
Zm00001eb078730Chitinase A14.442.03
Zm00001eb425600Chitinase B14.114.02
GH20
Zm00001eb008880Exochitinase 24.143.25
The cyclin-dependent kinase gene (Zm00001eb350890) was used for normalization. Genes with a Fold Change (FC) value ≥ 2 in RT-qPCR analyses were considered as a validated gene. RT-qPCR values are presented as the mean of the two biological replicates used for RNA-seq analysis with three technical replicates per each biological replicate. a Gene expression was not validated.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cazares-Álvarez, J.E.; Báez-Astorga, P.A.; Arroyo-Becerra, A.; Maldonado-Mendoza, I.E. Genome-Wide Identification of a Maize Chitinase Gene Family and the Induction of Its Expression by Fusarium verticillioides (Sacc.) Nirenberg (1976) Infection. Genes 2024, 15, 1087. https://doi.org/10.3390/genes15081087

AMA Style

Cazares-Álvarez JE, Báez-Astorga PA, Arroyo-Becerra A, Maldonado-Mendoza IE. Genome-Wide Identification of a Maize Chitinase Gene Family and the Induction of Its Expression by Fusarium verticillioides (Sacc.) Nirenberg (1976) Infection. Genes. 2024; 15(8):1087. https://doi.org/10.3390/genes15081087

Chicago/Turabian Style

Cazares-Álvarez, Jesús Eduardo, Paúl Alán Báez-Astorga, Analilia Arroyo-Becerra, and Ignacio Eduardo Maldonado-Mendoza. 2024. "Genome-Wide Identification of a Maize Chitinase Gene Family and the Induction of Its Expression by Fusarium verticillioides (Sacc.) Nirenberg (1976) Infection" Genes 15, no. 8: 1087. https://doi.org/10.3390/genes15081087

APA Style

Cazares-Álvarez, J. E., Báez-Astorga, P. A., Arroyo-Becerra, A., & Maldonado-Mendoza, I. E. (2024). Genome-Wide Identification of a Maize Chitinase Gene Family and the Induction of Its Expression by Fusarium verticillioides (Sacc.) Nirenberg (1976) Infection. Genes, 15(8), 1087. https://doi.org/10.3390/genes15081087

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop