Novel phylogenetic methods are needed for understanding gene function in the era of mega-scale genome sequencing.

Ongoing large-scale genome sequencing projects are forecasting a data deluge that will almost certainly overwhelm current analytical capabilities of evolutionary genomics. In contrast to population genomics, there are no standardized methods in evolutionary genomics for extracting evolutionary and functional (e.g. gene-trait association) signal from genomic data. Here, we examine how current practices of multi-species comparative genomics perform in this aspect and point out that many genomic datasets are under-utilized due to the lack of powerful methodologies. As a result, many current analyses emphasize gene families for which some functional data is already available, resulting in a growing gap between functionally well-characterized genes/organisms and the universe of unknowns. This leaves unknown genes on the 'dark side' of genomes, a problem that will not be mitigated by sequencing more and more genomes, unless we develop tools to infer functional hypotheses for unknown genes in a systematic manner. We provide an inventory of recently developed methods capable of predicting gene-gene and gene-trait associations based on comparative data, then argue that realizing the full potential of whole genome datasets requires the integration of phylogenetic comparative methods into genomics, a rich but underutilized toolbox for looking into the past.

Molecular characterization of ß-endoglucanase from antagonistic Trichoderma saturnisporum isolate GITX-Panog (C) induced under mycoparasitic conditions.

The endoglucanase belonging to glycoside hydrolase family 61 are little studied. In present study, a ß-endoglucanase of ~37 kDa induced on autoclaved mycelium of Fusarium oxysporum was cloned and characterized. The molecular characterization of ß-endoglucanase encoding gene revealed presence of a single intron and an open reading frame of 1044-bp which encoded a protein of 347 amino acid residues. The phylogenetic analysis of Eglu revealed its similarity to endo-ß-glucanases of other Trichoderma spp. The catalytic site of ß-endoglucanase contained Asp, Asn, His and Tyr residues. The cDNA encoding ß-glucanase was cloned into E. coli and Pichia pastoris using pQUA-30 and pPIC9K vector system, respectively. The comparison of structure revealed that most similar structure to Eglu is Hypocrea jecorina template 5o2w.1.A of glycoside hydrolase family 61.The biochemical characterization of ß-endoglucanase purified from T. saturnisporum isolate and the recombinant protein expressed in E. coli and P. pastoris was active under acidic conditions with a pH optima of 5 and temperature optima of 60 °C. The purified and expressed enzyme preparation was able to inhibit growth of F.oxysporum at 1 × 105 spores/mL which clearly revealed its significance in plant pathogen suppression.

Transcriptome-Guided Identification of Carbohydrate Active Enzymes (CAZy) from the Christmas Island Red Crab, Gecarcoidea natalis and a Vote for the Inclusion of Transcriptome-Derived Crustacean CAZys in Comparative Studies.

The Christmas Island red crab, Gecarcoidea natalis, is an herbivorous land crab that consumes mostly fallen leaf litter. In order to subsist, G. natalis would need to have developed specialised digestive enzymes capable of supplying significant amounts of metabolisable sugars from this diet. To gain insights into the carbohydrate metabolism of G. natalis, a transcriptome assembly was performed, with a specific focus on identifying transcripts coding for carbohydrate active enzyme (CAZy) using in silico approaches. Transcriptome sequencing of the midgut gland identified 70 CAZy-coding transcripts with varying expression values. At least three newly discovered putative GH9 endo-ß-1,4-glucanase ("classic cellulase") transcripts were highly expressed in the midgut gland in addition to the previously characterised GH9 and GH16 (ß-1,3-glucanase) transcripts, and underscoring the utility of whole transcriptome in uncovering new CAZy-coding transcripts. A highly expressed transcript coding for GH5_10 previously missed by conventional screening of cellulase activity was inferred to be a novel endo-ß-1,4-mannase in G. natalis with in silico support from homology modelling and amino acid alignment with other functionally validated GH5_10 proteins. Maximum likelihood tree reconstruction of the GH5_10 proteins demonstrates the phylogenetic affiliation of the G. natalis GH5_10 transcript to that of other decapods, supporting endogenous expression. Surprisingly, crustacean-derived GH5_10 transcripts were near absent in the current CAZy database and yet mining of the transcriptome shotgun assembly (TSA) recovered more than 100 crustacean GH5_10s in addition to several other biotechnological relevant CAZys, underscoring the unappreciated potential of the TSA database as a valuable resource for crustacean CAZys.

A New Salt-Tolerant Thermostable Cellulase from a Marine Bacillus sp. Strain.

A salt-tolerant cellulase secreted by a marine Bacillus sp. SR22 strain with wide resistance to temperature and pH was purified and characterized. Its approximate mass was 37 kDa. The endoglucanase, named as Bc22Cel, was purified by ammonium sulfate precipitation, gel filtration chromatography, and extraction from the gel after non-reducing sodium dodecyl sufate-polyacrylamide gel electrophoresis. The optimal pH value and temperature of Bc22Cel were 6.5 and 60°C, respectively. The purified Bc22Cel showed a considerable halophilic property, being able to maintain more than 70% of residual activity even when pre-incubated with 1.5 M NaCl for 1 h. Kinetic analysis of the purified enzyme showed the Km and Vmax to be 0.704 mg/ml and 29.85 µmol·ml-1·min-1, respectively. Taken together, the present data indicate Bc22Cel as a potential and useful candidate for industrial applications, such as the bioconversion of sugarcane bagasse to its derivatives.

Origin, evolution, and divergence of plant class C GH9 endoglucanases.

BACKGROUND: Glycoside hydrolases of the GH9 family encode cellulases that predominantly function as endoglucanases and have wide applications in the food, paper, pharmaceutical, and biofuel industries. The partitioning of plant GH9 endoglucanases, into classes A, B, and C, is based on the differential presence of transmembrane, signal peptide, and the carbohydrate binding module (CBM49). There is considerable debate on the distribution and the functions of these enzymes which may vary in different organisms. In light of these findings we examined the origin, emergence, and subsequent divergence of plant GH9 endoglucanases, with an emphasis on elucidating the role of CBM49 in the digestion of crystalline cellulose by class C members. RESULTS: Since, the digestion of crystalline cellulose mandates the presence of a well-defined set of aromatic and polar amino acids and/or an attributable domain that can mediate this conversion, we hypothesize a vertical mode of transfer of genes that could favour the emergence of class C like GH9 endoglucanase activity in land plants from potentially ancestral non plant taxa. We demonstrated the concomitant occurrence of a GH9 domain with CBM49 and other homologous carbohydrate binding modules, in putative endoglucanase sequences from several non-plant taxa. In the absence of comparable full length CBMs, we have characterized several low strength patterns that could approximate the CBM49, thereby, extending support for digestion of crystalline cellulose to other segments of the protein. We also provide data suggestive of the ancestral role of putative class C GH9 endoglucanases in land plants, which includes detailed phylogenetics and the presence and subsequent loss of CBM49, transmembrane, and signal peptide regions in certain populations of early land plants. These findings suggest that classes A and B of modern vascular land plants may have emerged by diverging directly from CBM49 encompassing putative class C enzymes. CONCLUSION: Our detailed phylogenetic and bioinformatics analysis of putative GH9 endoglucanase sequences across major taxa suggests that plant class C enzymes, despite their recent discovery, could function as the last common ancestor of classes A and B. Additionally, research into their ability to digest or inter-convert crystalline and amorphous forms of cellulose could make them lucrative candidates for engineering biofuel feedstock.

Functional diversity for biomass deconstruction in family 5 subfamily 5 (GH5_5) of fungal endo-ß1,4-glucanases.

Endo-ß1,4-glucanases in glycosyl hydrolase family 5 (GH5) are ubiquitous enzymes in the multicellular fungi and are common components of enzyme cocktails for biomass conversion. We recently showed that an endo-glucanase of subfamily 5 of GH5 (GH5_5) from Sporotrichum thermophile (StCel5A) was more effective at releasing glucose from pretreated corn stover, when part of an eight-component synthetic enzyme mixture, compared to its closely related counterpart from Trichoderma reesei, TrCel5A. StCel5A and TrCel5A belong to different clades of GH5_5 (GH5_5_1 and GH5_5_2, respectively). To test whether the superior activity of StCel5A was a general property of all enzymes in the GH5_5_2 clade, StCel5A, TrCel5A, and two additional members of each subfamily were expressed in a common host that had been engineered to suppress its native cellulases (T. reesei Δxyr1) and compared against each other alone on pure substrates, in synthetic mixtures on pure substrates, and against each other in synthetic mixtures on real biomass. The results indicated that superiority is a unique property of StCel5A and not of GH5_5_2 generally. The six Cel5A enzymes had significant differences in relative activities on different substrates, in specific activities, and in sensitivities to mannan inhibition. Importantly, the behavior of the six endo-glucanases on pure cellulose substrates did not predict their behavior in combination with other cellulolytic enzymes on a real lignocellulosic biomass substrate.

Characterization of a novel manganese dependent endoglucanase belongs in GH family 5 from Phanerochaete chrysosporium.

The cDNA encoding a putative glycoside hydrolase family 5, which has been predicted to be an endoglucanase (PcEg5A), was cloned from Phanerochaete chrysosporium and expressed in Pichia pastoris. PcEg5A contains a carbohydrate-binding domain and two important amino acids, E209 and E319, playing as proton donor and nucleophile in substrate catalytic domain. SDS-PAGE analysis indicated that the recombinant endoglucanase 5A (rPcEg5A) has a molecular size of 43 kDa which corresponds with the theoretical calculation. Optimum pH and temperature were found to be 4.5-6.0, and 50°C-60°C, respectively. Moreover, rPcEg5A exhibited maximal activity in the pH range of 3.0-8.0, whereas over 50% of activity still remained at 20°C and 80°C. rPcEg5A was stable at 60°C for 12 h incubation, indicating that rPcEg5A is a thermostable enzyme. Manganese ion enhanced the enzyme activity by 77%, indicating that rPcEg5A is a metal dependent enzyme. The addition of rPcEg5A to cellobiase (cellobiohydrolase and ß-glucosidase) resulted in a 53% increasing saccharification of NaOH-pretreated barley straw, whereas the glucose release was 47% higher than that cellobiase treatment alone. Our study suggested that rPcEg5A is an enzyme with great potential for biomass saccharification.

A sequence-anchored linkage map of the plant-parasitic nematode Meloidogyne hapla reveals exceptionally high genome-wide recombination.

Root-knot nematodes (Meloidogyne spp.) cause major yield losses to many of the world's crops, but efforts to understand how these pests recognize and interact with their hosts have been hampered by a lack of genetic resources. Starting with progeny of a cross between inbred strains (VW8 and VW9) of Meloidogyne hapla that differed in host range and behavioral traits, we exploited the novel, facultative meiotic parthenogenic reproductive mode of this species to produce a genetic linkage map. Molecular markers were derived from SNPs identified between the sequenced and annotated VW9 genome and de novo sequence of VW8. Genotypes were assessed in 183 F2 lines. The colinearity of the genetic and physical maps supported the veracity of both. Analysis of local crossover intervals revealed that the average recombination rate is exceptionally high compared with that in other metazoans. In addition, F2 lines are largely homozygous for markers flanking crossover points, and thus resemble recombinant inbred lines. We suggest that the unusually high recombination rate may be an adaptation to generate within-population genetic diversity in this organism. This work presents the most comprehensive linkage map of a parasitic nematode to date and, together with genomic and transcript sequence resources, empowers M. hapla as a tractable model. Alongside the molecular map, these progeny lines can be used for analyses of genome organization and the inheritance of phenotypic traits that have key functions in modulating parasitism, behavior, and survival and for the eventual identification of the responsible genes.

Molecular cloning of glycoside hydrolase family 45 cellulase genes from brackish water clam Corbicula japonica.

We previously reported endogenous Glycoside Hydrolase Family (GHF) 9 beta-1,4-glucanase gene, CjCel9A, from common Japanese freshwater clam Corbicula japonica. Here we identified another endogenous beta-1,4-glucanase genes which belong to GHF45 (CjCel45A, CjCel45B). Both genes encode ORF of 627 bp corresponding to 208 amino acids. CjCel45A and CjCel45B are different in 5' and 3'-untranslated regions and six nucleotides in the ORF. CjCEL45 has only one GHF45 catalytic domain without any carbohydrate binding modules as is the case with other molluskan GHF45 enzymes. Phylogenetic analysis and genomic structure of CjCel45 gene implies that this gene is likely to be acquired from fungi by common ancestor of mollusks. Reverse transcription (RT)-PCR analysis and in situ hybridization revealed that CjCel45A is likely to be expressed in the secretory cells in the digestive gland, suggesting that this cellulase is produced in the same site as CjCEL9A. CjCEL45A was successfully expressed in E. coli cells and zymographic analysis of the recombinant CjCEL45A showed that CjCEL45A is a functional beta-1,4-glucanase. The finding of multiple cellulase genes in C. japonica strongly supports our hypothesis that this species function as a cellulose decomposer in estuarine environments.

Enzyme diversity of the cellulolytic system produced by Clostridium cellulolyticum explored by two-dimensional analysis: identification of seven genes encoding new dockerin-containing proteins.

The enzyme diversity of the cellulolytic system produced by Clostridium cellulolyticum grown on crystalline cellulose as a sole carbon and energy source was explored by two-dimensional electrophoresis. The cellulolytic system of C. cellulolyticum is composed of at least 30 dockerin-containing proteins (designated cellulosomal proteins) and 30 noncellulosomal components. Most of the known cellulosomal proteins, including CipC, Cel48F, Cel8C, Cel9G, Cel9E, Man5K, Cel9M, and Cel5A, were identified by using two-dimensional Western blot analysis with specific antibodies, whereas Cel5N, Cel9J, and Cel44O were identified by using N-terminal sequencing. Unknown enzymes having carboxymethyl cellulase or xylanase activities were detected by zymogram analysis of two-dimensional gels. Some of these enzymes were identified by N-terminal sequencing as homologs of proteins listed in the NCBI database. Using Trap-Dock PCR and DNA walking, seven genes encoding new dockerin-containing proteins were cloned and sequenced. Some of these genes are clustered. Enzymes encoded by these genes belong to glycoside hydrolase families GH2, GH9, GH10, GH26, GH27, and GH59. Except for members of family GH9, which contains only cellulases, the new modular glycoside hydrolases discovered in this work could be involved in the degradation of different hemicellulosic substrates, such as xylan or galactomannan.

Cloning and identification of novel cellulase genes from uncultured microorganisms in rabbit cecum and characterization of the expressed cellulases.

A metagenomic cosmid library was prepared in Escherichia coli from DNA extracted from the contents of rabbit cecum and screened for cellulase activities. Eleven independent clones expressing cellulase activities (four endo-beta-1,4-glucanases and seven beta-glucosidases) were isolated. Subcloning and sequencing analysis of these clones identified 11 cellulase genes; the encoded products of which shared less than 50% identities and 70% similarities to cellulases in the databases. All four endo-beta-1,4-glucanases and all seven beta-glucosidases, respectively, belonged to glycosyl hydrolase family 5 (GHF 5) and family 3 (GHF 3) and formed two separate branches in the phylogenetic tree. Ten of the 11 cloned cellulases exhibited highest activities at pH 5.5 approximately 7.0 and 40 approximately 55 degrees C, a condition similar to that in the rabbit cecum. All the four endo-beta-1,4-glucanases could hydrolyze a wide range of beta-1,4-, beta-1,4/beta-1,3- or beta-1,3/beta-1,6-linked polysaccharides. One endo-beta-1, 4-glucanase gene, umcel5G, was overexpressed in E. coli, and the purified recombinant enzyme was characterized in detail. The enzymes cloned in this work represented at least some of the cellulases operating efficiently in the rabbit cecum. This work provides the first snapshot on the cellulases produced by bacteria in rabbit cecum.

Evaluation of cellulase preparations for hydrolysis of hardwood substrates.

Seven cellulase preparations from Penicillium and Trichoderma spp. were evaluated for their ability to hydrolyze the cellulose fraction of hardwoods (yellow poplar and red maple) pretreated by organosolv extraction, as well as model cellulosic substrates such as filter paper. There was no significant correlation among hydrolytic performance on pretreated hardwood, based on glucose release, and filter paper activity. However, performance on pretreated hardwood showed significant correlations to the levels of endogenous beta-glucosidase and xylanase activities in the cellulase preparation. Accordingly, differences in performance were reduced or eliminated following supplementation with a crude beta-glucosidase preparation containing both activities. These results complement a previous investigation using softwoods pretreated by either organosolv extraction or steam explosion. Cellulase preparations that performed best on hardwood also showed superior performance on the softwood substrates.

Specific characteristics of family 45 endoglucanases from Mucorales in the use of textiles and laundry.

We examined the characteristics of family 45 endoglucanases (glycoside hydrolases family 45; GH45) from Mucorales belonging to Zygomycota in the use of textiles and laundry. The defibrillation activities on lyocell fabric of family 45 endoglucanases from Mucorales, such as RCE1 and RCE2 from Rhizopus oryzae, MCE1 and MCE2 from Mucor circinelloides, and PCE1 from Phycomyces nitens, were much higher than those of the other family 45 endoglucanases. By contrast, family 45 endoglucanases from Mucorales were less resistant to anionic surfactant and oxidizing agent, main components in detergents, than the other family 45 endoglucanases. RCE1 consists of two distinct modules, a catalytic module and a carbohydrate-binding module family 1 (CBM1), and these common specific characteristics were considered to due to the catalytic module, but not to the CBM1.

Phylogenetic analysis of the plant endo-beta-1,4-glucanase gene family.

Phylogenetic analysis of the endo-beta-1,4-glucanase gene family of Arabidopsis and other plants revealed a clear distinction in three subfamilies (alpha, beta, and gamma). The alpha- and beta-subfamily contains proteins believed to be involved in a number of physiological roles such as elongation, ripening, and abscission. The gamma-subfamily is composed of proteins that are predicted to have a membrane-spanning domain and to be localized at the plasma membrane. Some of these proteins have been linked to cellulose biosynthesis by serving to hydrolyze a lipid-linked intermediate that acts as a primer for the elongation of beta-glucan chains during cellulose synthesis at the plasma membrane. Similar glucanases are important in cellulose biosynthesis in bacteria. Searches in the genomes of unrelated organisms that make cellulose, such as Ciona intestinalis and Dictyostelium discoideum, revealed the presence of membrane-linked endo-beta-1,4-glucanases and it is suggested that these might also have a role in cellulose synthesis.

Studies of Thermobifida fusca plant cell wall degrading enzymes.

I have been studying the Thermobifida fusca cellulose degrading proteins for the past 25 years. In this period, we have purified and characterized the six extracellular cellulases and an intracellular beta- glucosidase used by T. fusca for cellulose degradation, cloned and sequenced the structural genes encoding these enzymes, and helped to determine the 3-dimensional structures of two of the cellulase catalytic domains. This research determined the mechanism of a novel class of cellulase, family 9 processive endoglucanases, and helped to show that there were two types of exocellulases, ones that attacked the non-reducing ends of cellulose and ones that attacked the reducing ends. It also led to the sequencing of the T. fusca genome by the DOE Joint Genome Institute. We have studied the mechanisms that regulate T. fusca cellulases and have shown that cellobiose is the inducer and that cellulase synthesis is repressed by any good carbon source. A regulatory protein (CelR) that functions in the induction control has been purified, characterized, and its structural gene cloned and expressed in E. coli. I have also carried out research on two rumen bacteria, Prevotella ruminicola and Fibrobacter succinogenes, in collaboration with Professor James Russell, helping to arrange for the genomes of these two organisms to be sequenced by TIGR, funded by a USDA grant to the North American Consortium for Genomics of Fibrolytic Ruminal Biology.

Proteinaceous inhibitors of endo-beta-glucanases.

Both plants and filamentous phytopathogens secrete proteins that inhibit endo-beta-glucanases. The first endo-beta-glucanase inhibitor proteins to be discovered are XEGIP, a tomato protein that inhibits fungal xyloglucan-specific endo-beta-1,4-glucanases, and GIP1, an oomycete protein that inhibits endo-beta-1,3-glucanases produced by the plant host. These inhibitor proteins act by forming high-affinity complexes with their endoglucanase ligands. A family of XEGIP-like proteins has been identified. At least one member of this family (extracellular dermal glycoprotein, EDGP) has been shown to have endoglucanase-inhibitor activity, while other members have sequence similarity to a xylanase inhibitor from wheat (TAXI-1). The oomycete inhibitor GIP1 is a catalytically inactive serine protease homolog (SPH) whose structure is unrelated to XEGIP. Both types of inhibitor proteins are likely to affect the interactions of plants with filamentous phytopathogens, and a basic model describing their roles in pathogenesis is proposed.

Multiple forms of endo-1,4-beta-glucanases in the endosperm of Euphorbia heterophylla L.

Germinating seeds of Euphorbia heterophylla L. contain endo-1,4-beta-glucanases which degrade carboxymethylcellulose (CMC). The activity decreased approximately 66% in extracts of endosperm containing isopropanol or ethanol. The endoglucanases were isolated from endosperm extracts using ammonium sulphate fractionation followed by Sephacryl S-100-HR chromatography resulting in two main peaks: I and II. Peak I endoglucanase was further purified about 15-fold on DEAE-Sephadex A50 and then by affinity chromatography (CF11-cellulose). Peak II endoglucanases were further purified 10-fold on CM-cellulose chromatography. The results indicated the occurrence of a 66 kDa endoglucanase (fractionated by SDS-PAGE and visualized by activity staining using Congo Red). Several acidic (pI 3.0 to 5.7) and basic (pI 8.5 to 10.0) forms from both peaks which differed in their capacities for degrading CMC or xyloglucans from Copaifera langsdorffii or Hymenaea courbaril were detected.

Cellulase genes from the parabasalian symbiont Pseudotrichonympha grassii in the hindgut of the wood-feeding termite Coptotermes formosanus.

Cellulase genes of Pseudotrichonympha grassii (Hypermastigida: Eucomonymphidae), the symbiotic flagellate in the hindgut of the wood-feeding termite Coptotermes formosanus, were isolated and characterized. The nucleotide sequences of the major cellulase component in the hindgut of C. formosanus were determined based on its N-terminal amino acid sequence. The five isolated nucleotide sequences (PgCBH-homos) had an open reading frame of 1350 bp showing similarity to catalytic domains of glycoside hydrolase family (GHF) 7 members, and primary structure comparison with GHF7 members whose tertiary structures are well-characterized revealed the overall similarity between PgCBH-homo and the catalytic domain of a processive cellulase Cel7A (formerly CBHI) from the aerobic fungus Trichoderma reesei. Functional expression of PgCBH-homos in Escherichia coli, using the carboxymethylcellulose-Congo red assay, demonstrated the actual cellulolytic activity of PgCBH-homo. RT-PCR showed that PgCBH-homos were expressed, from the three flagellates in the hindgut, specifically in P. grassii.

Clostridium thermocellum cellulase CelT, a family 9 endoglucanase without an Ig-like domain or family 3c carbohydrate-binding module.

The celT gene of Clostridium thermocellum strain F1 was found downstream of the mannanase gene man26B [Kurokawa J et al. (2001) Biosci Biotechnol Biochem 65:548-554] in pKS305. The open reading frame of celT consists of 1,833 nucleotides encoding a protein of 611 amino acids with a predicted molecular weight of 68,510. The mature form of CelT consists of a family 9 cellulase domain and a dockerin domain responsible for cellulosome assembly, but lacks a family 3c carbohydrate-binding module (CBM) and an immunoglobulin (Ig)-like domain, which are often found with family 9 catalytic domains. CelT devoid of the dockerin domain (CelTDeltadoc) was constructed and purified from a recombinant Escherichia coli, and its enzyme properties were examined. CelTDeltadoc showed strong activity toward carboxymethylcellulose (CMC) and barley beta-glucan, and low activity toward xylan. The V(max) and K(m) values were 137 micro mol min(-1) mg(-1) and 16.7 mg/ml, respectively, for CMC. Immunological analysis indicated that CelT is a catalytic component of the C. thermocellum F1 cellulosome. This is the first report describing the characterization of a family 9 cellulase without an Ig-like domain or family 3c CBM.

Cloning and relational analysis of 15 novel fungal endoglucanases from family 12 glycosyl hydrolase.

Cellulases belong to the large family of glycosyl hydrolases (GHs) and are produced by a variety of bacteria and fungi. These extracellular enzymes act as endoglucanases (EGs), cellobiohydrolases or beta-glucosidases. In this paper, we describe molecular screening for EGs from the GH family 12. Using three homologous sequence boxes deduced from five previously known members of the family, we analysed 22 cellulase-producing fungal strains obtained from a diverse area of the fungal kingdom. Polymerase chain reactions using degenerate primers designed to the homologous protein boxes were used to identify the family 12 homologues. Several fungi showed the presence of multiple versions of the gene, while amino acid sequence analysis showed diversity in 15 novel members of the family, ranging from 26% to 96% similarity. Our sequence analysis shows that the phylogenetic tree of family 12 EGs can be divided into four subfamilies: 12-1 (fungal group I), 12-2 (fungal group II), 12-3 ( Streptomyces group in which Rhodothermus marinus fits) and 12-4 ( Thermophiles group). Erwinia carotovora may form a new subgroup.