- 1Laboratorio de Bioinformática Aplicada, Escuela de Ciencias Biológicas, Universidad Nacional de Costa Rica, Heredia, Costa Rica
- 2Laboratorio de Control Biológico, Escuela de Ciencias Agrarias, Universidad Nacional de Costa Rica, Heredia, Costa Rica
- 3Programa de Biotecnología Vegetal y Recursos Genéticos para el Fitomejoramiento (BIOVERFI), Escuela de Ciencias Agrarias, Universidad Nacional de Costa Rica, Heredia, Costa Rica
Background: The fungus Beauveria bassiana is widely used for agronomical applications, mainly in biological control. B. bassiana uses chitinase enzymes to degrade chitin, a major chemical component found in insect exoskeletons and fungal cell walls. However, until recently, genomic information on neotropical isolates, as well as their metabolic and biotechnological potential, has been limited.
Methods: Eight complete B. bassiana genomes of Neotropical origin and three references were studied to identify chitinase genes and its corresponding proteins, which were curated and characterized using manual curation and computational tools. We conducted a computational study to highlight functional differences and similarities for chitinase proteins in these Neotropical isolates.
Results: Eleven chitinase 1 genes were identified, categorized as chitinase 1.1 and chitinase 1.2. Five chitinase 2 genes were identified but presented a higher sequence conservation across all sequences. Interestingly, physicochemical parameters were more similar between chitinase 1.1 and chitinase 2 than between chitinase 1.1 and 1.2.
Conclusion: Chitinases 1 and 2 demonstrated variations, especially within chitinase 1, which presented a potential paralog. These differences were observed in their physical parameters. Additionally, CHIT2 completely lacks a signal peptide. This implies that CHIT1 might be associated with infection processes, while CHIT2 could be involved in morphogenesis and cellular growth. Therefore, our work highlights the importance of computational studies on local isolates, providing valuable resources for further experimental validation. Intrinsic changes within local species can significantly impact our understanding of complex pathogen-host interactions and offer practical applications, such as biological control.
1 Introduction
Beauveria bassiana (Ascomycota: Hypocreales) is a cosmopolitan fungus with the ability to survive as a saprophyte in soil and an as endophyte in plants; but can also act as an insect pathogen (Valero-Jiménez et al., 2016). The potential of this fungus to infect insects has been utilized in the agricultural industry to produce biopesticides (de Faria and Wraight, 2007). Unlike entomopathogenic viruses and bacteria, fungi do not infect insect hosts by ingestion. Instead, fungi initiate infection through penetration of the insect’s cuticle (Mascarin and Jaronski, 2016). To overcome the host’s cuticular layer, B. bassiana produces several enzymes capable of degrading cuticle constituents, among which proteases and chitinases stand out (Valero-Jiménez et al., 2016).
The conditions of the insect cuticle are considered adverse due to exposure to solar radiation, osmotic stress, and the insect’s defense responses. Therefore, rapid degradation of the cuticle would result in less exposure time to these conditions and thus a higher probability of infection success (Ortiz-Urquiza and Keyhani, 2016; Valero-Jiménez et al., 2016). For this reason, hydrolytic enzymes are considered important elements of pathogenesis. Chitinase is the most important enzyme in degrading the chitin polymer in the cuticular layer, and its activity has been associated with the virulence of entomopathogenic fungi (Pelizza et al., 2012; Dhawan and Joshi, 2017).
Chitinases are hydrolytic enzymes that degrade chitin into its oligo and monomeric components by hydrolyzing the β-1,4 N-acetyl-D-glucosaminide bonds (Le and Yang, 2019). According to their cleavage pattern, chitinases are divided into endo and exochitinases. Endochitinases can hydrolyze chitin at any point, generating products of variable size, while exochitinases do so from the non-reducing end of the chain, producing N,N′-diacetylchitobiose (Seidl, 2008). Additionally, based on the amino acid sequence, chitinases can be classified into families GH 18 and 19. These families show no sequence similarity to each other and differ in their catalytic mechanisms (Seidl, 2008).
Fungal chitinases not only play roles in pathogenesis, but also serve functions in morphogenesis, cell division, autolysis, among others (Le and Yang, 2019). The implication of chitinases in pathogenesis has prompted research development to enhance the understanding of these enzymes during the infection process. Despite numerous studies quantifying enzymatic production and gene expression encoding chitinases in B. bassiana (Fang et al., 2005; Pelizza et al., 2012; Al Khoury et al., 2019; Bhadani et al., 2021), there is still limited information regarding their structural and functional characteristics (Bhagwat et al., 2021). Several genomic resources are available for B. bassiana including whole genome sequences (Kim et al., 2016). More recently, eight assembled genomes for isolates of Neotropical origin have been published (Castro-Vásquez et al., 2022; Solano-González et al., 2023) and can serve as useful tools for comparison, as numerous organisms from this region are diverse.
The objective of this work was to compare, at an in silico level, gene sequences of two chitinases from Neotropical isolates of B. bassiana to determine their structural and functional characteristics to facilitate understanding of their biological functions. Furthermore, as building knowledge relying on previous genomic assemblies from these species (Solano-González et al., 2023), this study provides resources to further manipulate genetic sequences of chitinase through future experimental validation as a bottom line aiming to evaluate the potential of these fungi as biological control agents.
2 Materials and methods
2.1 Data obtention and chitinase coding sequences identification
Genomic sequences from six isolates of B. bassiana (Castro-Vásquez et al., 2022) belonging to the entomopathogenic fungi collection of the School of Agrarian Sciences (ECA) of the National University of Costa Rica (UNA) were used. Additionally, sequences from strains reported in the NCBI database were included as references (Supplementary Table 1).
For the identification of chitinase coding sequences, a local blastx was performed on the Laboratorio de Bioinformática Aplicada (LABAP) computational platform at UNA. All CHIT1 and CHIT2 protein sequences of B. bassiana reported in the Identical Protein Groups resource of NCBI (https://www.ncbi.nlm.nih.gov/ipg) were used as the database (Figure 1B). Matches were filtered based on identity percentage (≥90%). Subsequently, sequence fragments were retrieved using the Bedtools v.2.27.1 tool (Quinlan and Hall, 2010) and manually mapped within the genomes using Artemis v.18.0.0 (Rutherford et al., 2000) (Figures 1C, D) (Supplementary Table 2).
Figure 1. Schematic view of the methodology used to obtain the target sequences and their respective analyses. (A) Retrieval of chitinase protein sequences CHIT1 and CHIT2 from NCBI for database construction. (B) Identification of matches in the genomes of B. bassiana using Blastx. (C) Retrieval of sequence fragments using Bedtools v.2.27.1. (D) Manual identification of genes in the genomes using the Artemis tool v.18.0.0. (E) Verification of the identity of the complete gene sequences using an online Blastx. (F) Selection of genes based on coverage and identity criteria, annotation, and size. (G) Structural characterization using GSDS v2.0. (H) Translation of genes into proteins using the EMBOSS Transeq tool v.6.6.0. (I) Identification of conserved and non-conserved regions using Muscle v.3.8.31, Jalview v.2.11.6, and MEGA v.11. (J) In silico analysis of physicochemical and functional characteristics using SignalP, ProtParam InterPro, respectively; and analysis of secondary structures using Superfamily and SOPMA.
To verify the identity of the complete gene and its correspondence with chit1 and chit2 genes, an online blastx was performed (https://blast.ncbi.nlm.nih.gov/Blast.cgi) using default parameters (Figure 1E). Only the top three hits were considered as they exhibited the highest identity and coverage percentages (Supplementary Table 2). Subsequently, a final selection of these genes was carried out, to exclusively include genes relevant to the objectives of the study (those most similar to chit1 and chit2). The final selection was based on three criteria: 1) Genes whose identity and coverage percentage were ≥90 ± 3% with chit1 or chit2, 2) Genes already annotated as chitinase 1 or chitinase 2, however, those annotated differently from chit1 or chit2, were excluded, irrespective of whether they met the first criterion 3) Genes with a size greater than 3,000 bp were also discarded regardless of whether they satisfied the previous criteria (Figure 1F).
2.2 Chitinase structural in silico characterization
The structural annotation of the nucleotide sequences associated with chit1 and chit2 from the Neotropical fungal isolates was performed using GSDS v2.0 (http://gsds.gao-lab.org/). For this purpose, the complete sequences of each gene and their respective coding regions were used to visualize the location and size of introns and exons (Figure 1G).
2.3 Identification of conserved and non-conserved regions in chit1 and chit2
Translation of the peptide sequences corresponding to chit1 and chit2 genes from the studied genomes was carried out using the EMBOSS Transeq tool v.6.6.0 (https://www.ebi.ac.uk/Tools/st/emboss_transeq/) using default settings (Figure 1H). These sequences were used in conjunction with reference sequences to construct a multiple-sequence alignment using Muscle v.3.8.31. With the aim of visualizing sequence similarity, a multiple alignment of peptide sequences was performed. The alignments included sequences from Neotropical isolates and sequences from reference strains, totaling 27 sequences of CHIT1 and 11 sequences of CHIT2, respectively. The visualization was performed in Jalview v.2.11.6 to determine conserved and non-conserved regions among the different sequences (Figure 1I). A distance matrix was constructed using the aligned sequences of CHIT1 by replacing each amino acid with its numeric value based on relative entropy (Li et al., 2008). A hierarchical cluster dendrogram was drawn using the hclust function in RStudio (R Core Team, 2022).
On the other hand, for CHIT2 a Maximum likelihood tree was constructed from chit2 from B. bassiana draft isolates and references. This tree was constructed using Muscle alignments of 11 chit2 genes, implementing the JTT model and 1,000 replicates for bootstraps support, depicted in numbers. The tree shows Costa Rican isolates (BV-ECA 0, BV-ECA 26, and BV-ECA 31), Honduran isolate (BV-ECA 13), and the Puerto Rican isolate (BV-ECA 43) in addition to NCBI references.
2.4 In silico analysis of physicochemical and functional characteristics
SignalP6.0 (https://services.healthtech.dtu.dk/services/SignalP-6.0/) was used to predict the presence of potential signal peptides in all chitinases.
ExPASy ProtParam tool (https://web.expasy.org/protparam/) was used for the calculation and analysis of the physicochemical properties of the chitinases such as molecular weight, isoelectric point (pI), instability index (II), aliphatic index and the grand average of hydropathicity (GRAVY) were the evaluated.
The functional classification of the sequences corresponding to the Neotropical isolates was performed through a manual search and review of matches in protein databases associated with InterPro, for which its online version (https://www.ebi.ac.uk/interpro/search/sequence/) was used, implementing the default configuration (all selected databases). The SUPERFAMILY database (Pandurangan et al., 2019) was also used to determine the superfamily and subgroup associated with each of the chitinases (Figure 1J).
In addition, the prediction of secondary structures of chitinases was also carried out using the SOPMA web server (http://npsa-pbil.ibcp.fr/) (Geourjon and Deleage, 1995) maintaining the parameters established by default. Four conformational states were determined: helix, sheet, turn and coil.
A complete schematic view of the methodology is depicted in Figure 1.
3 Results
3.1 Data obtention and chitinase coding sequences identification
During data retrieval, a total of 50 genes similar to chit1 and 11 genes similar to chit2 were recovered (Supplementary Table 2). After the final selection process, 27 genes similar to chit1 met the established selection criteria, set as: identity and coverage percentage of ≥90 ± 3%, absence of gene annotations different from chit1 or chit2, as appropriate, and a size <3,000 bp (see Supplementary Tables 1, 2). Of these 27 genes, 11 belonged to the Neotropical isolates. As for the 11 chit2 genes, all met the selection criteria; and 5 belonged to the neotropical isolates. Selected gene sequences, as well as their respective protein sequences, were used in subsequent analyses.
3.2 Chitinase structural characterization
The structural characterization of the sequences corresponding to the ECA-UNA isolates shows that gene composition is very similar in terms of the quantity and distribution of introns. Regarding genes similar to chit1, only gene g5795 from the BV-ECA 27 isolate, did not present intronic regions. All other genes contained two introns of 44 and 48 bp or 53 and 59 bp, except for gene g1800 from BV-ECA 31, which had introns of 45 and 75 bp. Genes g5795 and g1800 showed the greatest difference in length, at 867 and 1,185 bp, respectively, while all others had lengths of 1,274 or 1,462 bp (Supplementary Figure 1). On the other hand, none of the sequences similar to chit2 from the neotropical isolates contain intronic regions. Furthermore, all these sequences present exactly a 1,044 bp size (Supplementary Figure 2).
3.3 Identification of conserved and non-conserved regions in CHIT1 and CHIT2
Most sequences similar to CHIT1 are grouped into two large clusters (Figures 2A, B). The first group, hereinafter referred to as CHIT1.1, encompassed 10 sequences, five of which were from ECA fungal accessions (g3951_BV-ECA43, g1021_BV-ECA26, g5905_BV-ECA, g1141_BV-ECA13, and g1800_BV-ECA31), while the remaining accessions were from reference strains. The second group, named CHIT1.2, was composed of ten sequences, five of these originated from the Neotropical isolates (g8786_BV-ECA31, g10088_BV-ECA43, g9108_BV-ECA26, g4716_BV-ECA0, and g899_BVECA13). Sequence alignment for the first group and the second group exhibited high conservation of residues throughout the entire peptide chain. On the other hand, the six remaining sequences that were not grouped show lesser similarity to the sequences of both groups (Figures 2A, B). Hence, these peptide sequences may also correspond to different chit genes. However, none of the six remaining sequences correspond to Neotropical isolates. Additionally, the clustering supported by the dendogram show CHIT1.1 and CHIT1.2 groups (Figure 2B).
Figure 2. Sequence alignment and phylogeny of CHIT1. (A) Partial alignment region of peptide sequences like chit1. Colored sections represent residues that are similar or identical between sequences. A higher intensity in color indicates a greater conservation of residues among the sequences. Areas with dashed lines represent gaps. On the left side of the image is the name of the sequence, consisting of the gene ID followed by the isolate or strain it comes from. In the figure, two groups of sequences are highlighted by labels, which show similarity in most of the residues. The first block is denoted as chit1.1 and the second block are denoted as chit1.2. (B) Hierarchical cluster dendrogram constructed based on a distance matrix from the aligned sequences of CHIT1 by replacing each amino acid with its numeric value based on relative entropy.
On the other hand, all CHIT2 sequences were highly conserved, with one or two amino acid substitutions. Therefore, both the peptide sequences from reference strains and the sequences from Neotropical isolates correspond to the same chit2 gene (Figures 3A, B).
Figure 3. Sequence alignment and phylogeny of CHIT2. (A) Partial alignment region of peptide sequences similar to chit2. The colored parts represent regions that are similar or identical between sequences. A higher intensity in color indicates a greater conservation of residues among the sequences. On the left side of the image is the name of the sequence, consisting of the gene ID followed by the isolate or strain. (B) Maximum likelihood tree constructed from chit2 genes from B. bassiana draft isolates and references. This tree was constructed using Muscle alignments of 11 chit2 genes, implementing the JTT model and 1,000 replicates for bootstraps support, depicted in numbers. The tree shows Costa Rican isolates (BV-ECA 0, BV-ECA 26, and BV-ECA 31), Honduran isolate (BV-ECA 13), and the Puerto Rican isolate (BV-ECA 43) in addition to NCBI references.
3.4 In silico analysis of physicochemical and functional characteristics
The presence of SEC/SPI signal peptides was predicted with a probability greater than 99% in the sequences of the CHIT1.1 group, except for g5795 BV-ECA 27 and g1800 BV-ECA 31. Additionally, the peptidase cleavage site of these chitinases was predicted between residues 21 and 22. In contrast with the CHIT1.2 group, the presence of signal peptides was predicted with a probability greater than 99% in all analyzed sequences, with cleavage site positions between residues 17 and 18. No signal peptides were detected in any of the sequences of the CHIT2 group.
The InterPro database was used to determine the classification and associated functions of the studied chitinases. All chitinases in the CHIT1.1 and CHIT1.2 clustered groups were classified as members of the glycosyl hydrolase family 18, except for chitinase g5795 from isolate BV-ECA 27, which was not assigned to any specific family. Additionally, members of the CHIT2 group were also not assigned to a particular family. However, the catalytic domain of glycosyl hydrolases was found in all three groups of chitinases.
In terms of functions, chitinases in the CHIT1.1 and CHIT1.2 groups were associated with the same GO terms: chitin binding and chitinase activity, except for g5795 BV-ECA 27, which again was not associated with any specific functional term. On the other hand, using the SUPERFAMILY database, all analyzed chitinases were grouped into the transglucosylases superfamily and classified as members of the type II chitinase group. Secondary structures of the chitinases were analyzed using the SOPMA tool (Supplementary Table 4). When comparing the ECA-UNA chitinases with the reference chitinases of each group, no significant variations were observed in the proportion of secondary structures, which coincided with the low standard deviation values. CHIT1.1 chitinases consistently displayed comparable proportions of alpha helix and random coil elements (37.82% and 39.76%, respectively), while beta sheets and beta turns were found in a smaller proportion (15.59% and 6.83%). On the other hand, CHIT1.2 and CHIT2 chitinases showed a greater difference between the proportion of random coil and alpha helix structures. In CHIT1.2, 46.88% of random coils and 34.62% of alpha helices were predicted, while in CHIT2, 44.81% of random coils and 30.23% of alpha helices were observed. Beta sheets were present in 13.02% in CHIT1.2% and 18.32% in CHIT2, while the proportion of beta turns was 5.49% and 6.63% in CHIT1.2 and CHIT2, respectively.
Regarding ExPASy ProtParam predictions, these showed chitinases from isolates BV-ECA 0, BV-ECA 13, BV-ECA 26 had the highest molecular weight (44000.79 Da) for the CHIT1.1 group, whereas the lowest weight was r found in isolate BV-ECA 27 (17344.35 Da), which did not group to either CHIT1.1 or CHIT1.2 groups. On the other hand, the molecular weight for chitinases from the CHIT1.2 group was very similar among all isolates, ranging from 49163.24 to 49198.42 Da. As for chitinase CHIT2, the estimated molecular weight for all sequences was 38501.57–38533.59 Da (Supplementary Table 3).
The isoelectric point (pI) was below 7 in all studied chitinases. For CHIT1.1, the pI value ranged from 4.72 to 4.94 among different isolates. For CHIT1.2, values between 5.16 and 5.42 were calculated, while the pI for chitinase CHIT2 was 4.64–4.68 (Supplementary Table 3).
Proteins with an instability index (II) below 40 are considered stable, while proteins with higher values are considered unstable (Ranjan et al., 2024). All studied chitinases had an instability index below 40, except for chitinase CHIT2 from isolate BV-ECA 31 (II 40.90), being the only chitinase classified as unstable. However, the estimated half-lives for all chitinases were greater than 20 h, according to the N-terminal rule and using the yeast model as a reference (Supplementary Table 3).
The aliphatic index (IA) of the CHIT1.1 chitinases varied between 75.43 and 85.37, while the value fluctuated between 61.60 and 61.83 for CHIT1.2 chitinases. CHIT2 chitinases exhibited IA values between 88.90 and 90.03. This suggests that all chitinases were considered thermally stable (Supplementary Table 3).
Both CHIT1 and CHIT2 chitinases showed negative values in the grand average of hydropathicity (GRAVY) parameter, indicating all chitinases are hydrophilic (Kyte and Doolittle, 1982; Uddin et al., 2020) (Supplementary Table 3).
4 Discussion
As a first approach toward its future exploitation, our work aimed to compare gene and protein sequences of two types of chitinases obtained from Neotropical isolates of B. bassiana. To determine their structural and functional characteristics and facilitate the understanding of their biological functions, we identified, manually curated and structurally annotated these sequences, as well as characterized their physical and biochemical characteristics. From the 16 chitinase genes identified in this study, we determined that chitinase 2 is more homogeneous than chitinase 1. This was evident based on the multiple sequence alignment and the dendrogram analysis of the protein sequences. We observed CHIT1.1 chitinases were similar to chitinases from isolates from South Korea (JEF-350: GenBank ID GCA_021365345.1) and Denmark (ARSEF8028: GenBank ID GCA_001682635.1); whereas CHIT1.2 chitinases grouped with isolates from South Korea (JEF-007: GenBank ID GCA_002871155.1) and China (D1-5: GenBank ID GCA_000770705.1). CHIT1.2 chitinases from Neotropical isolates grouped closer together to CHIT1.1 where isolate ECA-31 clustered with chitinases from other latitudes. This closer relationship of chitinase 1.1 from isolate ECA-31 suggests the genetic variations of B. bassiana may not necessarily be related to its geographic origin. In previous research, a low relationship between the geographic origin and genetic variation in Neotropical isolates of B. bassiana had been observed, including those analyzed in this study (Castro-Vásquez et al., 2021; Solano-González et al., 2023). It is possible to infer the presence of a potential paralog of CHIT1 (as seen in the CHIT.1 and CHIT1.2 gene groups). The presence of potential paralogous is relevant as higher number of copies with pathogenic function will impact the effectiveness of its infection processes (Gasmi et al., 2021), increasing the virulence of the fungus, therefore a relevant parameter to assess for future biological control applications. Chitinase g5795 BV-ECA is considered an outlier because, during the sequence retrieval, a frameshift was observed in the open reading frame, resulting in a protein with a start codon different from methionine. Therefore, to include this chitinase in the analysis, it was decided to define the start of the sequence at the nearest start codon, resulting in a chitinase with a significantly lower number of amino acids compared to the others.
Signal peptides are short sequences of amino acids that direct proteins toward a specific destination to fulfill a biological function (Owji et al., 2018; Teufel et al., 2022). The presence of the Secretory (Sec) Pathway and Signal Peptide (SEC/SPI peptide) in chitinases CHIT1.1 and CHIT1.2 indicated that these enzymes are likely secreted extracellularly. Although the presence of the peptide was not detected in all sequences of the chitinase 1.1 group, there is experimental evidence suggesting the involvement of chitinase CHIT1 in the infection process, indicating that this enzyme performs extracellular functions (Fang et al., 2005; Al Khoury et al., 2019). In contrast, none of the CHIT2 chitinases showed the presence of the signal peptide, suggesting that it may be an endogenous chitinase with functions related to morphogenesis and cell growth rather than being involved in the infection process. This could be associated with the high conservation of the CHIT2 sequences observed in the alignment, as being involved in basic functions for fungal development and might not require frequent environmental adaptations. We hypothesize that CHIT1 was exposed to greater evolutionary pressure due to interactions with the host and the environment (Vilcinskas, 2010), which in turn explains the variability among CHIT1 chitinases. Therefore, for the purposes of biological pest management, researchers and bioproduct developers should focus on the diverse extracellular chitinases, as these specifically interact with the host.
Functional analysis revealed that chitinases from all three groups (CHIT1.1, CHIT1.2, and CHIT2) possess the characteristic catalytic domain of glycosyl hydrolase family 18. Furthermore, they were classified within the same superfamily and subgroup of type II chitinases. Based on these findings, it is likely that CHIT1.1 and CHIT1.2 chitinases exhibit conserved functions and are products of a gene duplication event. However, the accumulation of amino acid sequence variations over time may introduce differences in physicochemical properties and secondary structures, leading to functionally similar proteins that could facilitate the adaptation of the fungus to diverse environmental conditions or different hosts.
Observed differences between the sequences of CHIT1.1 and CHIT1.2 are also reflected in their physicochemical characteristics. The latter is more water-soluble compared to the former, due to its lower GRAVY value (Kyte and Doolittle, 1982). Additionally, CHIT1.2 has a considerably lower aliphatic index, suggesting lower thermostability (Ikai, 1980). Differences were also identified in the isoelectric point (pI) and instability index; CHIT1.2 showed a less acidic pI and a lower instability index, which is associated with a longer in vivo half-life (Gasteiger et al., 2005). Interestingly, CHIT2 exhibited physicochemical parameters similar to CHIT1.1. The pI and GRAVY index values did not show significant differences between these two chitinases, unlike the comparisons between CHIT1.1 and CHIT1.2, or between CHIT1.2 and CHIT2. However, the latter showed the highest aliphatic index value, which can be associated with high thermal stability (Ikai, 1980). This characteristic may be of interest for biotechnological applications, as highly thermostable chitinases better withstand high temperatures and, therefore, have potential applications in fields such as agriculture, medicine, and the environment (Mathew et al., 2021). However, as it is likely an intracellular protein, its production may pose a greater challenge, as it requires genetic manipulation to direct its secretion extracellularly.
Liu et al. (2021) demonstrated that when exposed to different cuticle extracts, B. bassiana adjusted its infection strategy in response to new hosts. This adaptation involved differential regulation of chitinase expression within the same subgroup. In other words, the strain expressed different chitinases depending on the host being infected. Given B. bassiana’s broad host range (Oda et al., 2014), this plasticity may be linked to the fungus’s ability to vary its infection strategy. Moreover, this aligns with the results herein, regarding the abundance of GH18 family chitinases in this entomopathogen (Xiao et al., 2012).
Databases are invaluable resources for understanding the structural and functional characteristics of target proteins. However, the hierarchical nature of database annotations often results in non-specific or redundant information, hindering precise functional analysis, as consistently observed in our study. On the other hand, secondary structure analysis revealed a similar composition among chitinases from all three groups, suggesting significant structural conservation. Interestingly, chitinase CHIT1.2 exhibited a secondary structure composition more closely resembling CHIT2 than CHIT1.1. Previous studies by Bhagwat et al. (2021) on 15 B. bassiana chitinases yielded similar results, with a predominance of random coils and alpha helices. Given their shared ancestry, it is likely that the observed variations in secondary structure reflect the accumulation of changes over time, leading to functional adaptations.
There is abundant information related to the production of chitinases derived from fungi, however, reports detailing their purification and characteristics are scarce (Tupe et al., 2022). One of the study’s objectives was to provide resources to the scientific community that facilitate the development of future research in genetic and protein engineering. To serve as a foundation, we initially assessed the intrinsic physicochemical properties of chitinases knowing these parameters is significant for purification and isolation purposes (Ranjan et al., 2024), protein folding (led by hydrophobicity features) (Pinnamaneni et al., 2011; Stoykov et al., 2015), molecule interaction, functionality under variations of temperature (Islam et al., 2015) and pH conditions (Gomaa, 2021).
B. bassiana is a fungus known for its ability to affect a wide variety of hosts. However, different strains of this fungus show varied levels of virulence (Ortiz-Urquiza and Keyhani, 2016). This variability in virulence can pose a challenge for the commercialization of mycoinsecticides and their establishment as a strategy in integrated pest management. Gasmi et al. (2021) emphasized the importance of evaluating the genetic diversity of isolates and their effects on fungal performance. Similarly, Zhang et al. (2020) suggested that the study of genetic variants could be useful for improving the understanding of the different levels of virulence displayed by B. bassiana isolates towards their hosts. In the present study, we found possible genetic variants in the CHIT1 and CHIT2 chitinase sequences of Neotropical isolates. These variants could be the subject of further research to determine their possible effect on the phenotypic characteristics of the isolates, as well as the level of entomopathogenic potency among them.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author contributions
JS-V: Data curation, Formal Analysis, Investigation, Methodology, Software, Validation, Visualization, Writing–original draft, Writing–review and editing. AG-H: Writing–review and editing. RM-B: Writing–review and editing, Visualization. SS-G: Conceptualization, Funding acquisition, Methodology, Project administration, Resources, Software, Supervision, Validation, Writing–original draft, Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Universidad Nacional (https://ror.org/01t466c14) at Grant SIA-0717-19 and the Fondo Especial para la Educación Superior (FEES).
Acknowledgments
We want to thank the Centro Nacional de Computación Avanzada from CeNAT, CONARE for providing the cluster Kabré and the Centro de Gestión Tecnológica from Universidad Nacional for setting up the computational platform at the Laboratorio de Bioinformática Aplicada.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbinf.2024.1434442/full#supplementary-material
SUPPLEMENTARY FIGURE 1 | Genetic structure of chit1 in ECA-UNA isolates. The green color represents the exon while the dashed lines represent the introns.
SUPPLEMENTARY FIGURE 2 | Genetic structure of chit2 in ECA-UNA isolates. The green color represents the exon while the dashed lines represent the introns (absent in this case).
SUPPLEMENTARY TABLE 1 | Genes recovered from the search for matches between the ECA-UNA fungal genomes and the NCBI chitinases chit1 and chit2.
SUPPLEMENTARY TABLE 2 | Online blastx results for the verification of the identity of complete gene sequences.xlsx.
SUPPLEMENTARY TABLE 3 | In silico physicochemical properties of chitinase identified in Neotropical isolates of B. bassiana ECA-UNA.xlsx.
SUPPLEMENTARY TABLE 4 | Secondary structure analysis for CHIT1.1, CHIT1.2 and CHIT2.xlsx.
References
Al Khoury, C., Nemer, G., Guillot, J., Nour, A. A., and Nemer, N. (2019). Expression analysis of the genes involved in the virulence of Beauveria bassiana. Agri Gene 14, 100094. doi:10.1016/j.aggene.2019.100094
Bhadani, R. V., Gajera, H. P., Hirpara, D. G., Kachhadiya, H. J., and Dave, R. A. (2021). Metabolomics of extracellular compounds and parasitic enzymes of Beauveria bassiana associated with biological control of whiteflies (Bemisia tabaci). Pestic. Biochem. Physiol. 176, 104877. doi:10.1016/j.pestbp.2021.104877
Bhagwat, P., Amobonye, A., Singh, S., and Pillai, S. (2021). A comparative analysis of GH18 chitinases and their isoforms from Beauveria bassiana: an in-silico approach. Process Biochem. 100, 207–216. doi:10.1016/j.procbio.2020.10.012
Castro-Vásquez, R., Solano-González, S., Molina-Bravo, R., and Montero-Astúa, M. (2022). Draft genome sequences of eight isolates of Beauveria bassiana of neotropical origin. Microbiol. Resour. Announc 11, e0026722–22. doi:10.1128/mra.00267-22
Castro-Vásquez, R. M., Molina-Bravo, R., Hernández-Villalobos, S., Vargas-Martínez, A., González-Herrera, A., and Montero-Astúa, M. (2021). Identification and phylogenetic analysis of a collection of Beauveria spp. isolates from Central America and Puerto Rico. J. Invertebr. Pathol. 107642, 107642. doi:10.1016/j.jip.2021.107642
de Faria, M. R., and Wraight, S. P. (2007). Mycoinsecticides and Mycoacaricides: a comprehensive list with worldwide coverage and international classification of formulation types. Biol. Control 43, 237–256. doi:10.1016/j.biocontrol.2007.08.001
Dhawan, M., and Joshi, N. (2017). Enzymatic comparison and mortality of Beauveria bassiana against cabbage caterpillar Pieris brassicae LINN. Braz. J. Microbiol. 48, 522–529. doi:10.1016/j.bjm.2016.08.004
Fang, W., Leng, B., Xiao, Y., Jin, K., Ma, J., Fan, Y., et al. (2005). Cloning of Beauveria bassiana chitinase gene Bbchit1 and its application to improve fungal strain virulence. Appl. Environ. Microbiol. 71, 363–370. doi:10.1128/aem.71.1.363-370.2005
Gasmi, L., Baek, S., Kim, J. C., Kim, S., Lee, M. R., Park, S. E., et al. (2021). Gene diversity explains variation in biological features of insect killing fungus, Beauveria bassiana. Sci. Rep-uk 11, 91. doi:10.1038/s41598-020-78910-1
Gasteiger, E., Hoogland, C., Gattiker, A., Duvaud, S. E., Wilkins, M. R., Appel, R. D., et al. (2005). Protein identification and analysis tools on the ExPASy server. Methods Mol. Biol., 571–607. doi:10.1385/1-59259-890-0:571
Geourjon, C., and Deleage, G. (1995). SOPMA: significant improvements in protein secondary structure prediction by consensus prediction from multiple alignments. Bioinformatics 11 (6), 681–684. doi:10.1093/bioinformatics/11.6.681
Gomaa, E. Z. (2021). Microbial chitinases: properties, enhancement and potential applications. Protoplasma 258, 695–710. doi:10.1007/s00709-021-01612-6
Ikai, A. (1980). Thermostability and aliphatic index of globular proteins. J. Biochem. 88 (6), 1895–1898. doi:10.1093/oxfordjournals.jbchem.a133168
Islam, Md. S., Shahik, S.Md., Sohel, Md., Patwary, N. I. A., and Hasan, Md. A. (2015). In silico structural and functional annotation of hypothetical proteins of Vibrio cholerae O139. Genom. Inf. 13, 53–59. doi:10.5808/gi.2015.13.2.53
Kim, S., Lee, S. J., Nai, Y. S., Yu, J. S., Lee, M. R., Yang, Y. T., et al. (2016). Characterization of T-DNA insertion mutants with decreased virulence in the entomopathogenic fungus Beauveria bassiana JEF-007. Appl. Microbiol. Biotechnol. 100 (20), 8889–8900. doi:10.1007/s00253-016-7734-y
Kyte, J., and Doolittle, R. F. (1982). A simple method for displaying the hydropathic character of a protein. J. Mol. Biol. 157 (1), 105–132. doi:10.1016/0022-2836(82)90515-0
Le, B., and Yang, S. H. (2019). Microbial chitinases: properties, current state and biotechnological applications. World J. Microbiol. Biotechnol. 35, 144. doi:10.1007/s11274-019-2721-y
Li, C., Wang, J., Zhang, Y., and Wang, J. (2008). Similarity analysis of protein sequences based on the normalized relative entropy. Comb. Chem. High. Through Screen 11, 477–481. doi:10.2174/138620708784911500
Liu, J., Ling, Z., Wang, J., Xiang, T., Xu, L., Gu, C., et al. (2021). In vitro transcriptomes analysis identifies some special genes involved in pathogenicity difference of the Beauveria bassiana against different insect hosts. Microb. Pathog. 154, 104824. doi:10.1016/j.micpath.2021.104824
Mascarin, G. M., and Jaronski, S. T. (2016). The production and uses of Beauveria bassiana as a microbial insecticide. World J. Microbiol. Biotechnol. 32, 177. doi:10.1007/s11274-016-2131-3
Mathew, G. M., Madhavan, A., Arun, K. B., Sindhu, R., Binod, P., Singhania, R. R., et al. (2021). Thermophilic chitinases: structural, functional and engineering attributes for industrial applications. Appl. Biochem. Biotechnol. 193, 142–164. doi:10.1007/s12010-020-03416-5
Oda, H., Hatakeyama, Y., Yamamoto, Y., Enomoto, K., Shigano, T., and Iwano, H. (2014). Phylogenetic relationships among strains of the entomopathogenic fungus Beauveria bassiana (Hypocreales: clavicipitaceae) isolated from Japan. Appl. Entomology Zoology 49, 213–221. doi:10.1007/s13355-013-0239-6
Ortiz-Urquiza, A., and Keyhani, N. O. (2016). Molecular genetics of Beauveria bassiana infection of insects. Adv. Genet. 94, 165–249. doi:10.1016/bs.adgen.2015.11.003
Owji, H., Nezafat, N., Negahdaripour, M., Hajiebrahimi, A., and Ghasemi, Y. (2018). A comprehensive review of signal peptides: structure, roles, and applications. Eur. J. Cell Biol. 97 (6), 422–441. doi:10.1016/j.ejcb.2018.06.003
Pandurangan, A. P., Stahlhacke, J., Oates, M. E., Smithers, B., and Gough, J. (2019). The SUPERFAMILY 2.0 database: a significant proteome update and a new webserver. Nucleic Acids Res. 47 (D1), D490–D494. doi:10.1093/nar/gky1130
Pelizza, S. A., Elíades, L. A., Scorsetti, A. C., Cabello, M. N., and Lange, C. E. (2012). Entomopathogenic fungi from Argentina for the control of Schistocerca cancellata (Orthoptera: acrididae) nymphs: fungal pathogenicity and enzyme activity. Biocontrol Sci. Technol. 22, 1119–1129. doi:10.1080/09583157.2012.713910
Pinnamaneni, R., Kalidas, P., and Rao, K. R. S. S. (2011). Studies on the cloning and expression of Bbchit1 gene of Beauveria bassiana NCIM 1216. Indian J. Microbiol. 51, 396–402. doi:10.1007/s12088-011-0116-3
Quinlan, A. R., and Hall, I. M. (2010). BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics 26, 841–842. doi:10.1093/bioinformatics/btq033
Ranjan, P., Yadav, A., Behera, A. K., Singh, D. K., Singh, P., and Singh, G. P. (2024). Structural, functional and evolutionary analysis of wheat WRKY45 protein: a combined bioinformatics and MD simulation approach. Céréal. Res. Commun. 52, 489–499. doi:10.1007/s42976-023-00399-x
R Core Team (2022). R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing.
Rutherford, K., Parkhill, J., Crook, J., Horsnell, T., Rice, P., Rajandream, M.-A., et al. (2000). Artemis: sequence visualization and annotation. Bioinformatics 16, 944–945. doi:10.1093/bioinformatics/16.10.944
Seidl, V. (2008). Chitinases of filamentous fungi: a large group of diverse proteins with multiple physiological functions. Fungal Biol. Rev. 22, 36–42. doi:10.1016/j.fbr.2008.03.002
Solano-González, S., Castro-Vásquez, R., and Molina-Bravo, R. (2023). Genomic characterization and functional description of Beauveria bassiana isolates from Latin America. J. Fungi 9, 711. doi:10.3390/jof9070711
Stoykov, Y. M., Pavlov, A. I., and Krastanov, A. I. (2015). Chitinase biotechnology: production, purification, and application. Eng. Life Sci. 15, 30–38. doi:10.1002/elsc.201400173
Teufel, F., Almagro Armenteros, J. J., Johansen, A. R., Gíslason, M. H., Pihl, S. I., Tsirigos, K. D., et al. (2022). SignalP 6.0 predicts all five types of signal peptides using protein language models. Nat. Biotechnol. 40 (7), 1023–1025. doi:10.1038/s41587-021-01156-3
Tupe, S. G., Pathan, E. K., Ganger, S., Patil, S., and Deshpande, M. V. (2022). Fungal chitinolytic enzymes. Biol. Biotechnol. Appl., 185–201. doi:10.1007/978-981-16-3307-2_7
Uddin, A., Mazumder, T. H., Barbhuiya, P. A., and Chakraborty, S. (2020). Similarities and dissimilarities of codon usage in mitochondrial ATP genes among fishes, aves, and mammals. IUBMB life 72 (5), 899–914. doi:10.1002/iub.2231
Valero-Jiménez, C. A., Faino, L., Veld, D. S. in’t, Smit, S., Zwaan, B. J., and Kan, J. A. L. van (2016). Comparative genomics of Beauveria bassiana: uncovering signatures of virulence against mosquitoes. Bmc Genomics 17, 986. doi:10.1186/s12864-016-3339-1
Vilcinskas, A. (2010). Coevolution between pathogen-derived proteinases and proteinase inhibitors of host insects. Virulence 1 (3), 206–214. doi:10.4161/viru.1.3.12072
Xiao, G., Ying, S. H., Zheng, P., Wang, Z. L., Zhang, S., Xie, X. Q., et al. (2012). Genomic perspectives on the evolution of fungal entomopathogenicity in Beauveria bassiana. Sci. Rep. 2, 483. doi:10.1038/srep00483
Keywords: entomopathogen, bioinformatics, gene structure, fungi, physicochemical parameters, genomics, biocontrol
Citation: Segura-Vega J, González-Herrera A, Molina-Bravo R and Solano-González S (2024) Computational identification and characterization of chitinase 1 and chitinase 2 from neotropical isolates of Beauveria bassiana. Front. Bioinform. 4:1434442. doi: 10.3389/fbinf.2024.1434442
Received: 17 May 2024; Accepted: 30 September 2024;
Published: 18 October 2024.
Edited by:
David W. Ussery, University of Arkansas for Medical Sciences, United StatesReviewed by:
Leandro de Mattos Pereira, University of Porto, PortugalYaqian Li, Shanghai Jiao Tong University, China
Copyright © 2024 Segura-Vega, González-Herrera, Molina-Bravo and Solano-González. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Stefany Solano-González, stefany.solano.gonzalez@una.ac.cr