Abstract

PlantCAZyme is a database built upon dbCAN (database for automated carbohydrate active enzyme annotation), aiming to provide pre-computed sequence and annotation data of carbohydrate active enzymes (CAZymes) to plant carbohydrate and bioenergy research communities. The current version contains data of 43 790 CAZymes of 159 protein families from 35 plants (including angiosperms, gymnosperms, lycophyte and bryophyte mosses) and chlorophyte algae with fully sequenced genomes. Useful features of the database include: (i) a BLAST server and a HMMER server that allow users to search against our pre-computed sequence data for annotation purpose, (ii) a download page to allow batch downloading data of a specific CAZyme family or species and (iii) protein browse pages to provide an easy access to the most comprehensive sequence and annotation data.

Database URL: http://cys.bios.niu.edu/plantcazyme/

Introduction

Lignocellulosic biofuels have received great attentions in the past decade for obvious economic and environmental reasons [1]. Other than using starch-based plant materials as the feedstock, lignocellulosic biofuels use inedible plant biomass materials, which however are very recalcitrant to be degraded to release fermentable sugars. The bioenergy research community thus has major interests in genetically modifying plants in order to develop low-cost biofuels [2]. To achieve this goal, researchers need to know which genes should be modified to acquire the desired plants with lower recalcitrance to enzymatic degradation. Therefore biomass-related enzyme databases are highly needed to promote the development of transgenic biofuel crops [3]. Carbohydrate-Active enzymes (CAZymes) are enzymes responsible for the synthesis, degradation and modification of storage and structural biomass polysaccharides [4] and thus are the most important enzymes for bioenergy research. CAZymes are not only found in plants and bacteria, but also in fungi and animals, responsible for the synthesis, degradation and modification of all the glycoconjugates in nature including glycoproteins and glycolipids. Therefore they are also fundamentally important for general carbohydrate and glycobiology research [4].

CAZymes are present in all life kingdoms and particularly abundant in plants [5]. Since 1998, the CAZyme database, known as CAZy, has started to collect experimentally (biochemically, genetically and structurally) characterized CAZyme proteins and classify them into protein families and so far has created 330 families (as of May 2013) of six classes based on sequence homology: GHs (glycoside hydrolases), GTs (glycosyltransferases), CEs (carbohydrate esterases), PLs (polysaccharide lyases), AAs (auxiliary activities) and CBMs (carbohydrate binding modules) [6]. It then populated each family by including homologs from GenBank, UniProt and PDB databases using both BLAST and protein domain/motif search strategies as well as expert manual inspection of sequence alignment [4, 7]. CAZy is an extremely useful resource for its most original classification scheme and high-quality manual curation, and thus has been widely accepted by the carbohydrate research community.

A great demand of an automated CAZyme annotation emerged in the past few years due to the production of thousands of completed plant and microbial genomes and metagenomes. However CAZy database does not provide automated CAZyme annotation. In view of this need, in 2012 we have developed a web server named dbCAN, to allow users to submit the newly sequenced genomes for an automated CAZyme annotation [8]. Behind the web server are hidden Markov models (HMMs) of the 330 CAZyme families; each HMM represents the sequence alignment of conserved signature domains of each family, which were retrieved from annotated CAZyme protein sequences of the CAZy database. dbCAN website has received over thousands of visits from many countries after publication, demonstrating its impact on the research of CAZymes.

The availability of the 330 CAZyme HMMs has also made it possible to build a dedicated database for plant CAZymes. With regard to similar resources, the CAZy database covers only two (Arabidopsis thaliana and Oryzasativa) out of over 40 sequenced plant and algal genomes; all sequenced bioenergy crops (e.g. poplar, switchgrass, soghum) and evolutionarily important organisms (e.g. moss, spike moss, algae) were not included. Two other databases, pDAWG [9] and Rice GT [10], are limited to a small number of CAZyme families and genomes. There are also a few other databases such as the Cell Wall Genomics database [11] and the Cell Wall Navigator database [12], which only contain a very small number of CAZyme families. Therefore, the development of PlantCAZyme is a timely and highly significant addition to the toolbox of plant carbohydrate and bioenergy research.

Construction and Content

Collection of CAZyme sequences

Over 40 plant and algal genomes are completed and most of them are available in the Phytozome database [13]. To collect the plant CAZyme protein sequences, we used 330 dbCAN HMMs as query and scanned 35 genomes (Table 1), including 34 Phytozome genomes of 23 dicots, six monocots, one moss, one spike moss, two chlorophyte algae, as well as one gymnosperm genome [14] that is not available in Phytozome, using the HMMER 3.0 package as the homology search tool [15] with default parameters (E-value < 10 and output in parseable table of per-domain hits). The HMMER output was further processed to keep the significant hits as described in below.

Table 1.

Thirty-five plant and algal genomes that are included in the PlantCAZyme database

SpeciesCladeSource# of genes# of CAZyme genes% of CAZyme genes
Volvox carteriChlorophytePhytozome14 9711981.32
Chlamydomonas reinhardtiiChlorophytePhytozome20 4972851.39
Physcomitrella patensBryophytaPhytozome21 1738574.05
Selaginella moellendorffiiLycophytaPhytozome22 2859194.12
Picea abiesGymnospermCongenie71 15818432.59
Aquilegia coeruleaDicotPhytozome24 82310994.43
Arabidopsis lyrataDicotPhytozome32 67012323.77
Arabidopsis thalianaDicotPhytozome27 41612244.46
Brassica rapaDicotPhytozome40 90518124.43
Capsella rubellaDicotPhytozome26 52112114.57
Carica papayaDicotPhytozome27 7698453.04
Citrus clementinaDicotPhytozome24 55310984.47
Citrus sinensisDicotPhytozome25 37910834.27
Cucumis sativusDicotPhytozome21 50310084.69
Eucalyptus grandisDicotPhytozome36 37617114.70
Fragaria vescaDicotPhytozome65 66211051.68
Glycine maxDicotPhytozome54 17523544.35
Gossypium raimondiiDicotPhytozome37 50516484.39
Linum usitatissimumDicotPhytozome43 47120184.64
Malus domesticaDicotPhytozome63 51422203.50
Manihot esculentaDicotPhytozome30 66614424.70
Medicago truncatulaDicotPhytozome44 13511732.66
Mimulus guttatusDicotPhytozome26 71812714.76
Phaseolus vulgarisDicotPhytozome27 19713514.97
Populus trichocarpaDicotPhytozome41 33517514.24
Prunus persicaDicotPhytozome27 86412884.62
Ricinus communisDicotPhytozome31 22111353.64
Thellungiella halophilaDicotPhytozome26 35111324.30
Vitis viniferaDicotPhytozome26 34610964.16
Brachypodium distachyonMonocotPhytozome26 55212434.68
Oryza sativaMonocotPhytozome39 23413633.47
Panicum virgatumMonocotPhytozome65 87826243.98
Setaria italicaMonocotPhytozome35 47114874.19
Sorghum bicolorMonocotPhytozome27 60813344.83
Zea maysMonocotPhytozome39 65614753.72
SpeciesCladeSource# of genes# of CAZyme genes% of CAZyme genes
Volvox carteriChlorophytePhytozome14 9711981.32
Chlamydomonas reinhardtiiChlorophytePhytozome20 4972851.39
Physcomitrella patensBryophytaPhytozome21 1738574.05
Selaginella moellendorffiiLycophytaPhytozome22 2859194.12
Picea abiesGymnospermCongenie71 15818432.59
Aquilegia coeruleaDicotPhytozome24 82310994.43
Arabidopsis lyrataDicotPhytozome32 67012323.77
Arabidopsis thalianaDicotPhytozome27 41612244.46
Brassica rapaDicotPhytozome40 90518124.43
Capsella rubellaDicotPhytozome26 52112114.57
Carica papayaDicotPhytozome27 7698453.04
Citrus clementinaDicotPhytozome24 55310984.47
Citrus sinensisDicotPhytozome25 37910834.27
Cucumis sativusDicotPhytozome21 50310084.69
Eucalyptus grandisDicotPhytozome36 37617114.70
Fragaria vescaDicotPhytozome65 66211051.68
Glycine maxDicotPhytozome54 17523544.35
Gossypium raimondiiDicotPhytozome37 50516484.39
Linum usitatissimumDicotPhytozome43 47120184.64
Malus domesticaDicotPhytozome63 51422203.50
Manihot esculentaDicotPhytozome30 66614424.70
Medicago truncatulaDicotPhytozome44 13511732.66
Mimulus guttatusDicotPhytozome26 71812714.76
Phaseolus vulgarisDicotPhytozome27 19713514.97
Populus trichocarpaDicotPhytozome41 33517514.24
Prunus persicaDicotPhytozome27 86412884.62
Ricinus communisDicotPhytozome31 22111353.64
Thellungiella halophilaDicotPhytozome26 35111324.30
Vitis viniferaDicotPhytozome26 34610964.16
Brachypodium distachyonMonocotPhytozome26 55212434.68
Oryza sativaMonocotPhytozome39 23413633.47
Panicum virgatumMonocotPhytozome65 87826243.98
Setaria italicaMonocotPhytozome35 47114874.19
Sorghum bicolorMonocotPhytozome27 60813344.83
Zea maysMonocotPhytozome39 65614753.72
Table 1.

Thirty-five plant and algal genomes that are included in the PlantCAZyme database

SpeciesCladeSource# of genes# of CAZyme genes% of CAZyme genes
Volvox carteriChlorophytePhytozome14 9711981.32
Chlamydomonas reinhardtiiChlorophytePhytozome20 4972851.39
Physcomitrella patensBryophytaPhytozome21 1738574.05
Selaginella moellendorffiiLycophytaPhytozome22 2859194.12
Picea abiesGymnospermCongenie71 15818432.59
Aquilegia coeruleaDicotPhytozome24 82310994.43
Arabidopsis lyrataDicotPhytozome32 67012323.77
Arabidopsis thalianaDicotPhytozome27 41612244.46
Brassica rapaDicotPhytozome40 90518124.43
Capsella rubellaDicotPhytozome26 52112114.57
Carica papayaDicotPhytozome27 7698453.04
Citrus clementinaDicotPhytozome24 55310984.47
Citrus sinensisDicotPhytozome25 37910834.27
Cucumis sativusDicotPhytozome21 50310084.69
Eucalyptus grandisDicotPhytozome36 37617114.70
Fragaria vescaDicotPhytozome65 66211051.68
Glycine maxDicotPhytozome54 17523544.35
Gossypium raimondiiDicotPhytozome37 50516484.39
Linum usitatissimumDicotPhytozome43 47120184.64
Malus domesticaDicotPhytozome63 51422203.50
Manihot esculentaDicotPhytozome30 66614424.70
Medicago truncatulaDicotPhytozome44 13511732.66
Mimulus guttatusDicotPhytozome26 71812714.76
Phaseolus vulgarisDicotPhytozome27 19713514.97
Populus trichocarpaDicotPhytozome41 33517514.24
Prunus persicaDicotPhytozome27 86412884.62
Ricinus communisDicotPhytozome31 22111353.64
Thellungiella halophilaDicotPhytozome26 35111324.30
Vitis viniferaDicotPhytozome26 34610964.16
Brachypodium distachyonMonocotPhytozome26 55212434.68
Oryza sativaMonocotPhytozome39 23413633.47
Panicum virgatumMonocotPhytozome65 87826243.98
Setaria italicaMonocotPhytozome35 47114874.19
Sorghum bicolorMonocotPhytozome27 60813344.83
Zea maysMonocotPhytozome39 65614753.72
SpeciesCladeSource# of genes# of CAZyme genes% of CAZyme genes
Volvox carteriChlorophytePhytozome14 9711981.32
Chlamydomonas reinhardtiiChlorophytePhytozome20 4972851.39
Physcomitrella patensBryophytaPhytozome21 1738574.05
Selaginella moellendorffiiLycophytaPhytozome22 2859194.12
Picea abiesGymnospermCongenie71 15818432.59
Aquilegia coeruleaDicotPhytozome24 82310994.43
Arabidopsis lyrataDicotPhytozome32 67012323.77
Arabidopsis thalianaDicotPhytozome27 41612244.46
Brassica rapaDicotPhytozome40 90518124.43
Capsella rubellaDicotPhytozome26 52112114.57
Carica papayaDicotPhytozome27 7698453.04
Citrus clementinaDicotPhytozome24 55310984.47
Citrus sinensisDicotPhytozome25 37910834.27
Cucumis sativusDicotPhytozome21 50310084.69
Eucalyptus grandisDicotPhytozome36 37617114.70
Fragaria vescaDicotPhytozome65 66211051.68
Glycine maxDicotPhytozome54 17523544.35
Gossypium raimondiiDicotPhytozome37 50516484.39
Linum usitatissimumDicotPhytozome43 47120184.64
Malus domesticaDicotPhytozome63 51422203.50
Manihot esculentaDicotPhytozome30 66614424.70
Medicago truncatulaDicotPhytozome44 13511732.66
Mimulus guttatusDicotPhytozome26 71812714.76
Phaseolus vulgarisDicotPhytozome27 19713514.97
Populus trichocarpaDicotPhytozome41 33517514.24
Prunus persicaDicotPhytozome27 86412884.62
Ricinus communisDicotPhytozome31 22111353.64
Thellungiella halophilaDicotPhytozome26 35111324.30
Vitis viniferaDicotPhytozome26 34610964.16
Brachypodium distachyonMonocotPhytozome26 55212434.68
Oryza sativaMonocotPhytozome39 23413633.47
Panicum virgatumMonocotPhytozome65 87826243.98
Setaria italicaMonocotPhytozome35 47114874.19
Sorghum bicolorMonocotPhytozome27 60813344.83
Zea maysMonocotPhytozome39 65614753.72

Selection of golden standard datasets for accuracy benchmark

Since the CAZymes of Arabidopsis and rice have been annotated in the CAZy database, we have used these two genomes to calculate the sensitivity (or recall) and positive predictive value (or precision) of our CAZyme data. It is worth mentioning that the ‘annotated’ CAZymes of CAZy include not only experimentally characterized proteins, but also proteins that are deemed to be true homologs of the characterized proteins. For example, there are only three Arabidopsis proteins experimentally characterized to be GH17 enzymes (http://www.cazy.org/GH17_characterized.html); however 51 Arabidopsis proteins are listed as GH17 enzymes (http://www.cazy.org/GH17_eukaryota.html). The reason is that CAZy database annotates CAZymes from the GenBank database, including those from Arabidopsis and rice, by combining homology search and expert curation (e.g. manual inspection of sequence alignment for characteristic amino acid motifs [7]). Most of the Arabidopsis CAZymes including those experimentally uncharacterized have been manually curated by CAZy developers and published in 2001 [16]. The similar approach has also been applied to the annotation of poplar CAZymes in 2006 [17]. Due to its high-quality manual curation and rich functional annotation, CAZy was used as a golden standard dataset to assess automated CAZyme annotation by the CAZymes Analysis Toolkit (CAT) [18] and the dbCAN database [8].

There are also other protein family and function classification databases such as Pfam [19], KOG (eukaryotic orthologous groups) [20], KEGG Orthology (KO) [21], SUPERFAMILY [22], PANTHER [23], Gene Ontology (GO) [24] and many others. Each database has its own strength and focus (e.g. on protein domain or evolution or pathway or structure) and has much redundancy among each other (i.e. one protein family is described in multiple databases). Therefore integration efforts such as InterPro database [25] and CDD database [26] attempted to integrate all these different protein family databases into one framework to remove redundancy. Many of these resources are extremely useful for genome annotation purpose. For example, in the plant genomics community Phytozome [13], Gramene [27] and PLAZA [28] used the above resources to construct and compare protein families across different plants. In addition, ENZYME database [29] created the nomenclature system (i.e. the Enzyme Commission/EC numbers) of all characterized enzymes and associated biochemical reactions. Other databases such as Priam [30], CatFam [31], EFICAz [32] and PlantCyc [33] employed the EC classification system to either define enzyme family models or reconstruct metabolic pathways.

However, unlike CAZy, dbCAN and PlantCAZyme, all the above resources are not specifically designed for CAZymes but rather are general protein family/classification databases. As their mission is to cover all protein families in nature as broadly as possible, they do not have a focus and often miss some families of certain protein class, which is one of the reasons for the need of many specialized databases for individual protein families/classes such as [6, 34–37] (see more at http://www.oxfordjournals.org/nar/database/subcat/3/10). For example, Pfam only covers 142 out of 330 CAZyme families [8]. As a matter of fact, most of these 142 families were initially defined and annotated (from literature curation) by CAZy database and then were included into Pfam as HMMs, which makes Pfam not an ideal resource for CAZyme annotation. In addition, it is well known that one single CAZyme family could contain proteins with different biochemical activities and one biochemical activity could be carried by multiple CAZyme families [4]. For example, the CAZyme GH5 family contains characterized proteins with 20 different EC numbers (manually curated at http://www.cazy.org/GH5.html) and the cellulase (EC 3.2.1.4) activity is found in more than 10 GH families [38]. This makes it impossible to compare dbCAN HMM-based search and EC-based databases (e.g. Priam and CatFam) in terms of CAZyme assignment. Therefore, one cannot evaluate the CAZyme family assignment by comparing to the general protein family/classification databases. Since we aim to assess if we have retrieved all CAZyme homologs using the HMMs built from CAZy annotated proteins, CAZy database is naturally selected as the gold standard dataset to evaluate our performance.

Accuracy benchmark with Arabidopsis and rice data

As discussed in our dbCAN article [8], two criteria significantly impact the sensitivity and precision of our automated CAZyme annotation. One is E-value and the other is coverage, which is defined to measure the fraction of CAZyme domains covered in the alignment. We have tested the performance of dbCAN-based search on all of the CAZyme families as a whole (denoted as All) using different combinations of E-values and overage cutoffs. Figure 1 shows the F-measure values of different parameter combinations for the All sets of Arabidopsis (Figure 1A) and rice (Figure 1B), where F-measure = 2 × (Sensitivity × Precision) / (Sensitivity + Precision). We then selected the combination that gave the highest F-measure value and presented them in Tables 2 and 3. The more detailed information about how to calculate Sensitivity and Precision is provided in the Supplementary Tables S1–S12.

Evaluation of the impact of E-value and coverage parameters to the accuracy of pre-computed PlantCAZyme sequence data for Arabidopsis and rice; x-axis (horizontal): E-value, y-axis (vertical): F-measure, Z-axis: coverage. For both species, E-value < 1e–23 and coverage > 0.2 gave the highest F-measure. The detailed calculations are provided in Supplementary Table S1 and Supplementary Data.
Figure 1.

Evaluation of the impact of E-value and coverage parameters to the accuracy of pre-computed PlantCAZyme sequence data for Arabidopsis and rice; x-axis (horizontal): E-value, y-axis (vertical): F-measure, Z-axis: coverage. For both species, E-value < 1e–23 and coverage > 0.2 gave the highest F-measure. The detailed calculations are provided in Supplementary Table S1 and Supplementary Data.

Tables 2 and 3 show that the coverage >0.2 and E-value < 1e-23 combination gave the best F-measure for both Arabidopsis (F-measure = 0.91, sensitivity = 0.89 and precision = 0.92) and rice (F-measure = 0.85, sensitivity = 0.84 and precision = 0.85). We have also performed evaluation for the five CAZyme classes separately, which suggests that the best F-measure varies for different CAZyme classes (Tables 2 and 3). Overall the largest two classes GT and GH (81% of CAZyme families) in both plants have higher F-measures than the three smaller classes CE, PL and CBM. It also suggests that: (i) to annotate GH proteins, one should use a very relax coverage cutoff or the sensitivity will be low (Supplementary Tables S4 and Supplementary Data); (ii) to annotate CE families a very stringent E-value cutoff and coverage cutoff should be used; otherwise the precision will be very low due to a very high false positive rate (Supplementary Tables S5 and Supplementary Data). Although it would work best to use different parameter combinations for different CAZyme classes and for different plants, we decided to use coverage > 0.2 and E-value < 1e-23 as the universal threshold, as this setting agrees in both dicots and monocots and makes the parsing process less complicated and easy to reproduce by others.

Table 2.

The E-value and Coverage cutoffs that lead to the best F-measure in Arabidopsis

Arabidopsis# of CAZyme familiesE-valueCoverageF-measureSensitivityPrecision
All981.00E-230.20.9092367620.8940719140.924924925
GT431.00E-110.250.9376344090.9478260870.927659574
GH361.00E-160.050.9748110830.9699248120.979746835
CE51.00E-290.950.9457411340.9176470590.975609756
PL21.00E-300.250.9705882350.9705882350.970588235
CBM101.00E-120.750.796137730.8214285710.772357724
Arabidopsis# of CAZyme familiesE-valueCoverageF-measureSensitivityPrecision
All981.00E-230.20.9092367620.8940719140.924924925
GT431.00E-110.250.9376344090.9478260870.927659574
GH361.00E-160.050.9748110830.9699248120.979746835
CE51.00E-290.950.9457411340.9176470590.975609756
PL21.00E-300.250.9705882350.9705882350.970588235
CBM101.00E-120.750.796137730.8214285710.772357724
Table 2.

The E-value and Coverage cutoffs that lead to the best F-measure in Arabidopsis

Arabidopsis# of CAZyme familiesE-valueCoverageF-measureSensitivityPrecision
All981.00E-230.20.9092367620.8940719140.924924925
GT431.00E-110.250.9376344090.9478260870.927659574
GH361.00E-160.050.9748110830.9699248120.979746835
CE51.00E-290.950.9457411340.9176470590.975609756
PL21.00E-300.250.9705882350.9705882350.970588235
CBM101.00E-120.750.796137730.8214285710.772357724
Arabidopsis# of CAZyme familiesE-valueCoverageF-measureSensitivityPrecision
All981.00E-230.20.9092367620.8940719140.924924925
GT431.00E-110.250.9376344090.9478260870.927659574
GH361.00E-160.050.9748110830.9699248120.979746835
CE51.00E-290.950.9457411340.9176470590.975609756
PL21.00E-300.250.9705882350.9705882350.970588235
CBM101.00E-120.750.796137730.8214285710.772357724
Table 3.

The E-value and coverage cutoffs that lead to the best F-measure in Rice

Rice# of CAZyme familiesE-valueCoverageF-measureSensitivityPrecision
All971.00E-230.20.8451696810.8406193080.849769585
GT441.00E-100.350.9063817930.9089316990.903846154
GH351.00E-130.10.924153310.917452830.930952381
CE51.00E-280.950.9135452520.9056603770.921568627
PL21.00E-300.70.8275862070.750.923076923
CBM91.00E-160.450.7160316320.8571428570.614814815
Rice# of CAZyme familiesE-valueCoverageF-measureSensitivityPrecision
All971.00E-230.20.8451696810.8406193080.849769585
GT441.00E-100.350.9063817930.9089316990.903846154
GH351.00E-130.10.924153310.917452830.930952381
CE51.00E-280.950.9135452520.9056603770.921568627
PL21.00E-300.70.8275862070.750.923076923
CBM91.00E-160.450.7160316320.8571428570.614814815
Table 3.

The E-value and coverage cutoffs that lead to the best F-measure in Rice

Rice# of CAZyme familiesE-valueCoverageF-measureSensitivityPrecision
All971.00E-230.20.8451696810.8406193080.849769585
GT441.00E-100.350.9063817930.9089316990.903846154
GH351.00E-130.10.924153310.917452830.930952381
CE51.00E-280.950.9135452520.9056603770.921568627
PL21.00E-300.70.8275862070.750.923076923
CBM91.00E-160.450.7160316320.8571428570.614814815
Rice# of CAZyme familiesE-valueCoverageF-measureSensitivityPrecision
All971.00E-230.20.8451696810.8406193080.849769585
GT441.00E-100.350.9063817930.9089316990.903846154
GH351.00E-130.10.924153310.917452830.930952381
CE51.00E-280.950.9135452520.9056603770.921568627
PL21.00E-300.70.8275862070.750.923076923
CBM91.00E-160.450.7160316320.8571428570.614814815

Annotation data

We have further generated extensive bioinformatics annotation data for the plant CAZyme sequences by running various bioinformatics tools against different databases. As shown in Figure 2, these data include functional annotation (conserved functional domains, Gene Ontology annotation, top matches in the non-redundant protein database [NCBI-nr] and expressed sequence tag (EST) database), structural annotation [top matches in the Protein Data Bank (PDB), predicted transmembrane domains, signal peptides, coiled regions, hydropathy plot], phylogenetic annotation (orthologous groups of the CAZyme domains, multiple sequence alignment, phylogenetic tree) and miscellaneous data (nucleotide coding sequences, CAZyme signature domain sequences, genomic location, external links, publications, etc.).

A schematic architecture of the PlantCAZyme database
Figure 2.

A schematic architecture of the PlantCAZyme database

Utility and Discussion

Implementation and user interface

All the data were integrated and presented through a web interface powered by MySQL+PHP+JavaScript. As shown in Figure 2, the protein centric display page is used to present the sequence and annotation of each CAZyme protein. The website has a download page that allows users to download CAZyme sequences of a particular species or a particular CAZyme family. Both the CAZyme signature domain sequences and the full-length sequences are available for any species or any family.

A BLAST page and a HMMER (annotate) page were included to allow users to submit their own sequences for annotation, which are very useful to annotate sequences that are not included in our database. For BLAST search, users can submit both protein and nucleotide sequences and the databases for BLAST search can be chosen from: (i) the CAZy database that contains full-length GenBank protein sequences annotated in the CAZy database, (ii) the plant CAZyme domain sequences (not the full length) that are compiled in our PlantCAZyme database containing the CAZyme signature domains identified by dbCAN search. The results are returned as a webpage with a tabular output of the BLAST program.

For HMMER page, users must submit protein sequences as query and the database is the dbCAN’s HMMs. Since HMMs are built for each CAZyme family to represent the signature domain, this type of search is a better way than BLAST search to annotate new protein sequences with the modular CAZyme domain architecture.

In addition to sequence search, the keyword search function was also implemented. The top-right corner of each webpage has a search box, where users can search the database with a keyword. There are two options for keyword search: unformatted searching and formatted searching. For unformatted searching you enter a query with no formatting. This will run the query only against the following fields: (i) ID, e.g. AT2G46570.1, (ii) Family, e.g. CBM10, (iii) Species, e.g. A.thaliana and (iv) Domain, e.g. Cellulose_synt. Formatted searching allows users to be more specific and search through more fields. Formatted searches are done by indicating formatting with the use of brackets []. For example, if users want to search for the species A.thaliana, they can search ‘Arabidopsis thaliana[Species]’, which will bring up anything with a species containing ‘Arabidopsis’ or ‘thaliana’. Users can write more than one specifier in a query. So if users only wanted the AA1 family, they could write the query as ‘Arabidopsis[Species] thaliana[Species] AA1[Family]’. These specifiers are all strung together in an AND fashion, so a result will only appear if it matches all of the criteria users have given. Currently the keyword search only allows exact match and does not allow partial match and wildcard, which will be considered in the future.

A help page is designed to provide all necessary information for browsing, querying, downloading and searching the website and the database.

Use cases

If users want to retrieve all CAZyme proteins of A. thaliana, there will be three options. (i) Users can go to the download page, browse by species and locate the species to download the FASTA format sequences of full-length proteins or just the CAZyme domains. (ii) They can also go to the homepage, browse by species, click on the species and link to the family browse page of A. thaliana. There they can view which CAZyme families are in A. thaliana and how many genes are in each family, as well as a clickable genomic location plot. This Arabidopsis thaliana browse page also has a link to the complete HMMER output, where hits that did not pass our filters (coverage > 0.3 and E-value < 1e-5) can also be retrieved. Clicking on each family will present a new page with the list of proteins of that family, and further clicking on the ID will open the protein browse page. (iii) The last way is to perform a keyword search in the following format: (Arabidopsis thaliana)[species] or Arabidopsis[Species] thaliana[Species], which will return a table with all the Arabidopsis thaliana CAZyme IDs.

Similarly, if users want to retrieve CAZyme proteins of a specific family, say GT8, they will have the three options too: (i) download all GT8 proteins at the download page, (ii) browse by family at the homepage and (iii) use the keyword search function: GT8[family].

If users have a dataset (e.g. a newly sequenced genome) to be annotated for CAZymes, they can upload the FASTA sequences to our computing server through the BLAST page or the annotate (HMMER) page. The job will be run and the result will be returned with the CAZyme match information. If a huge dataset (>5000 sequences) needs to be processed, we recommend that users download the BLAST databases (CAZyDB or PlantCAZyme) or the HMM database (dbCAN) at our download page and run the searches on their local computers.

Future work

We plan to update the database at least once a year. We plan to include more species in the future, particularly selected plants and algae that do not have completed genomes. We will use transcriptomes of species such as ferns, liverworts, charophytic green algae (CGA), basal angiosperms, as they are important for the evolutionary study of CAZymes in plants and algae. The automatic collection of CAZyme sequences will also be further improved, e.g. by considering applying different parsing thresholds for different plant clades and by supplementing the HMMER search with BLAST search. We will also develop new web applications to display duplicated genes and orthologous genes of CAZymes on the chromosomes to allow comparative and evolutionary study of CAZymes.

PlantCAZyme is the first web resource dedicated to provide pre-computed CAZyme sequence and annotation data for all sequenced plants and algae. We expect it will be a highly useful tool to the plant cell wall and bioenergy research communities.

Acknowledgements

A.E. computed most of the data and implemented the database and website. R.T. contributed to the data collection. N.M. developed an early version of the database. Y.Y. conceived the database, supervised the entire project, computed some of the data and wrote the paper. The authors acknowledge the Department of Computer Science of NIU for providing free access to the Linux computing cluster Gaea and our lab members for helpful discussions. The authors thank all the reviewers for their good suggestions to improve this article.

Funding

Y.Y. was funded by the Research & Artistry Award and the startup package from Northern Illinois University. A.E. and N.M. were supported by Undergraduate Research Assistantships. Funding for open access charge: Northern Illinois University Libraries.

Conflict of interest. None declared.

References

1

Rubin
E.M.
(
2008
)
Genomics of cellulosic biofuels
.
Nature
,
454
,
841
845
.

2

Himmel
M.E.
Ding
S.Y.
Johnson
D.K.
et al. . (
2007
)
Biomass recalcitrance: engineering plants and enzymes for biofuels production
.
Science
,
315
,
804
807
.

3

Yin
Y.
(
2014
) In:
Gupta
V.
Tuohy
M.
Kubicek
C.
Saddler
J.
Xu
F.
(eds).
Bioenergy Research: Advances and Applications
.
Elsevier BV
,
The Netherlands
, pp.
95
107
.

4

Cantarel
B.L.
Coutinho
P.M.
Rancurel
C.
et al. . (
2009
)
The Carbohydrate-Active EnZymes database (CAZy): an expert resource for Glycogenomics
.
Nucleic Acids Res.
,
37
,
D233D238.

5

Coutinho
P.M.
Stam
M.
Blanc
E.
et al. . (
2003
)
Why are there so many carbohydrate-active enzyme-related genes in plants?
Trend Plant Sci.
,
8
,
563565.

6

Lombard
V.
Golaconda Ramulu
H.
Drula
E.
et al. . (
2014
)
The carbohydrate-active enzymes database (CAZy) in 2013
.
Nucleic Acids Res.
,
42
,
D490
D495
.

7

Coutinho
P.M.
Henrissat
B.
(
2010
),
Annual Plant Reviews
.
New Jersey, United States
:
Wiley-Blackwell
, pp.
93
107
.

8

Yin
Y.B.
Mao
X.Z.
Yang
J.C.
et al. . (
2012
)
dbCAN: a web resource for automated carbohydrate-active enzyme annotation
.
Nucleic Acids Res.
,
40
,
W445
W451
.

9

Mao
F.L.
Yin
Y.B.
Zhou
F.F.
et al. . (
2009
)
pDAWG: An integrated database for plant cell wall genes
.
Bioenerg Res.
,
2
,
209
216
.

10

Cao
P.J.
Bartley
L.E.
Jung
K.H.
et al. . (
2008
)
Construction of a rice glycosyltransferase phylogenomic database and identification of rice-diverged glycosyltransferases
.
Mol. Plant
,
1
,
858
877
.

11

Yong
W.
Link
B.
O'Malley
R.
et al. . (
2005
)
Genomics of plant cell wall biogenesis
.
Planta
,
221
,
747
751
.

12

Girke
T.
Lauricha
J.
Tran
H.
et al. . (
2004
)
The cell wall navigator database. A systems-based approach to organism-unrestricted mining of protein families involved in cell wall metabolism
.
Plant Physiol.
,
136
,
3003
3008
.

13

Goodstein
D.M.
Shu
S.
Howson
R.
et al. . (
2012
)
Phytozome: a comparative platform for green plant genomics
.
Nucleic Acids Res.
,
40
,
D1178
D1186
.

14

Nystedt
B.
Street
N.R.
Wetterbom
A.
et al. . (
2013
)
The Norway spruce genome sequence and conifer genome evolution
.
Nature
,
497
,
579
584
.

15

Eddy
S.R.
(
2011
)
Accelerated profile HMM searches
.
PLoS Computat. Biol.
,
7
,
e1002195.

16

Henrissat
B.
Coutinho
P.M.
Davies
G.J.
(
2001
)
A census of carbohydrate-active enzymes in the genome of Arabidopsis thaliana
.
Plant Mol. Biol.
,
47
,
55
72
.

17

Geisler-Lee
J.
Geisler
M.
Coutinho
P.M.
et al. . (
2006
)
Poplar carbohydrate-active enzymes. Gene identification and expression analyses
.
Plant Physiol.
,
140
,
946
962
.

18

Park
B.H.
Karpinets
T.V.
Syed
M.H.
et al. . (
2010
)
CAZymes Analysis Toolkit (CAT): Web service for searching and analyzing carbohydrate-active enzymes in a newly sequenced organism using CAZy database
.
Glycobiology
,
20
,
1574
1584
.

19

Finn
R.D.
Bateman
A.
Clements
J.
et al. . (
2014
)
Pfam: the protein families database
.
Nucleic Acids Res.
,
42
,
D222
D230
.

20

Tatusov
R.L.
Fedorova
N.D.
Jackson
J.D.
et al. . (
2003
)
The COG database: an updated version includes eukaryotes
.
BMC Bioinformatics
,
4
,
41
.

21

Kanehisa
M.
Goto
S.
Sato
Y.
et al. . (
2012
)
KEGG for integration and interpretation of large-scale molecular data sets
.
Nucleic Acids Res.
,
40
,
D109
D114
.

22

Gough
J.
Chothia
C.
(
2002
)
SUPERFAMILY: HMMs representing all proteins of known structure. SCOP sequence searches, alignments and genome assignments
.
Nucleic Acids Res.
,
30
,
268
272
.

23

Mi
H.Y.
Muruganujan
A.
Thomas
P.D.
(
2013
)
PANTHER in 2013: modeling the evolution of gene function, and other gene attributes, in the context of phylogenetic trees
.
Nucleic Acids Res.
,
41
,
D377
D386
.

24

Ashburner
M.
Ball
C.A.
Blake
J.A.
et al. . (
2000
)
Gene ontology: tool for the unification of biology
.
Nat Genet
,
25
,
25
29
.

25

Hunter
S.
Apweiler
R.
Attwood
T.K.
et al. . (
2009
)
InterPro: the integrative protein signature database
.
Nucleic Acids Res.
,
37
,
D211
D215
.

26

Marchler-Bauer
A.
Lu
S.N.
Anderson
J.B.
et al. . (
2011
)
CDD: a conserved domain database for the functional annotation of proteins
.
Nucleic Acids Res.
,
39
,
D225
D229
.

27

Ware
D.
Jaiswal
P.
Ni
J.J.
et al. . (
2002
)
Gramene: a resource for comparative grass genomics
.
Nucleic Acids Res.
,
30
,
103
105
.

28

Van Bel
M.
Proost
S.
Wischnitzki
E.
et al. . (
2012
)
Dissecting plant genomes with the PLAZA comparative genomics platform
.
Plant Physiol.
,
158
,
590
600
.

29

Bairoch
A.
(
2000
)
The ENZYME database in 2000
.
Nucleic Acids Res.
,
28
,
304
305
.

30

Claudel-Renard
C.
Chevalet
C.
Faraut
T.
et al. . (
2003
)
Enzyme-specific profiles for genome annotation: PRIAM
.
Nucleic Acids Res.
,
31
,
6633
6639
.

31

Yu
C.G.
Zavaijevski
N.
Desai
V.
et al. . (
2009
)
Genome-wide enzyme annotation with precision control: Catalytic families [CatFam] databases
.
Proteins
,
74
,
449
460
.

32

Tian
W.D.
Arakaki
A.K.
Skolnick
J.
(
2004
)
EFICAz: a comprehensive approach for accurate genome-scale enzyme function inference
.
Nucleic Acids Res.
,
32
,
6226
6239
.

33

Mueller
L.A.
Zhang
P.F.
Rhee
S.Y.
(
2003
)
AraCyc: a biochemical pathway database for Arabidopsis
.
Plant Physiol.
,
132
,
453
460
.

34

Guo
A.Y.
Chen
X.
Gao
G.
et al. . (
2008
)
PlantTFDB: a comprehensive plant transcription factor database
.
Nucleic Acids Res.
,
36
,
D966
D969
.

35

Fawal
N.
Li
Q.
Savelli
B.
et al. . (
2013
)
PeroxiBase: a database for large-scale evolutionary analysis of peroxidases
.
Nucleic Acids Res.
,
41
,
D441
4
.

36

Saier
M.H.
Jr.
Tran
C.V.
Barabote
R.D.
(
2006
)
TCDB: the transporter classification database for membrane transport protein analyses and information
.
Nucleic Acids Res.
,
34
,
D181
D186
.

37

Rawlings
N.D.
Barrett
A.J.
Bateman
A.
(
2012
)
MEROPS: the database of proteolytic enzymes, their substrates and inhibitors
.
Nucleic Acids Res.
,
40
,
D343
D350
.

38

Sukharnikov
L.O.
Cantwell
B.J.
Podar
M.
et al. . (
2011
)
Cellulases: ambiguous nonhomologous enzymes in a genomic perspective
.
Trends Biotechnol.
,
29
,
473
479
.

Author notes

Citation details: Ekstrom,A., Taujale,R., McGinn,N. et al. PlantCAZyme: a database for plant carbohydrate-active enzymes. Database (2014) Vol. 2014: article ID bau079; doi:10.1093/database/bau079

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.

Supplementary data