Analysis of functional divergence within two structurally related glycoside hydrolase families

Silencing gram-negative bacteria binding protein 1 decreases the immunity of Tribolium castaneum against bacteria

Gram-negative bacteria binding proteins (GNBPs) have the ability to recognize molecular patterns associated with microbial pathogens (PAMPs), leading to the activation of immune responses downstream. In the genome of Tribolium castaneum, three GNBP genes have been identified; however, their immunological roles remain unexplored. In our study, a GNBP1, designated as TcGNBP1, were identified from the cDNA library of T. castaneum. The coding sequence of TcGNBP1 consisted of 1137 bps and resulted in the synthesis of a protein comprising 378 amino acids. This protein encompasses a signal peptide, a low-complexity region, and a glycoside hydrolase 16 domain. TcGNBP1 was strongly expressed in early adult stages, and mainly distributed in hemolymph and gut. Upon being challenged with Escherichia coli or Staphylococcus aureus, the transcript levels of TcGNBP1 were significantly changed at different time points. Through molecular docking and ELISA analysis, it was observed that TcGNBP1 has the ability to interact with lipopolysaccharides, peptidoglycan, and β-1, 3-glucan. Based on these findings, it was further discovered that recombinant TcGNBP1 can directly bind to five different bacteria in a Ca2+-dependent manner. After knockdown of TcGNBP1 with RNA interference, expression of antimicrobial peptide genes and prophenoloxidase (proPO) activity were suppressed, the susceptibility of T. castaneum to E. coli or S. aureus infection was enhanced, leading to low survival rate. These results suggest a regulatory mechanism of TcGNBP1 in innate immunity of T. castaneum and provide a potential molecular target for dsRNA-based insect pest management.

Phylogenetic systematics of Butyrivibrio and Pseudobutyrivibrio genomes illustrate vast taxonomic diversity, open genomes and an abundance of carbohydrate-active enzyme family isoforms

Butyrivibrio and Pseudobutyrivibrio dominate in anaerobic gastrointestinal microbiomes, particularly the rumen, where they play a key role in harvesting dietary energy. Within these genera, five rumen species have been classified (Butyrivibrio fibrisolvens, Butyrivibrio hungatei, Butyrivibrio proteoclasticus, Pseudobutyrivibrio ruminis and Pseudobutyrivibrio xylanivorans) and more recently an additional Butyrivibrio sp. group was added. Given the recent increase in available genomes, we re-investigated the phylogenetic systematics and evolution of Butyrivibrio and Pseudobutyrivibrio. Across 71 genomes, we show using 16S rDNA and 40 gene marker phylogenetic trees that the current six species designations (P. ruminis, P. xylanivorans, B. fibrisolvens, Butyrivibrio sp., B. hungatei and B. proteclasticus) are found. However, pangenome analysis showed vast genomic variation and a high abundance of accessory genes (91.50–99.34 %), compared with core genes (0.66–8.50 %), within these six taxonomic groups, suggesting incorrectly assigned taxonomy. Subsequent pangenome accessory genomes under varying core gene cut-offs (%) and average nucleotide identity (ANI) analysis suggest the existence of 42 species within 32 genera. Pangenome analysis of those that still group within B. fibrisolvens, B. hungatei and P. ruminis, based on revised ANI phylogeny, also showed possession of very open genomes, illustrating the diversity that exists even within these groups. All strains of both Butyrivibrio and Pseudobutyrivibrio also shared a broad range of clusters of orthologous genes (COGs) (870), indicating recent evolution from a common ancestor. We also demonstrate that the carbohydrate-active enzymes (CAZymes) predominantly belong to glycosyl hydrolase (GH)2, 3, 5, 13 and 43, with numerous within family isoforms apparent, likely facilitating metabolic plasticity and resilience under dietary perturbations. This study provides a major advancement in our functional and evolutionary understanding of these important anaerobic bacteria.

Glucanases and Chitinases

In many yeast and fungi, β-(1,3)-glucan and chitin are essential components of the cell wall, an important structure that surrounds cells and which is responsible for their mechanical protection and necessary for maintaining the cellular shape. In addition, the cell wall is a dynamic structure that needs to be remodelled along with the different phases of the fungal life cycle or in response to extracellular stimuli. Since β-(1,3)-glucan and chitin perform a central structural role in the assembly of the cell wall, it has been postulated that β-(1,3)-glucanases and chitinases should perform an important function in cell wall softening and remodelling. This review focusses on fungal glucanases and chitinases and their role during fungal morphogenesis.

Systematic deletions in the cellobiohydrolase (CBH) Cel7A from the fungus Trichoderma reesei reveal flexible loops critical for CBH activity

Glycoside hydrolase family 7 (GH7) cellulases are some of the most efficient degraders of cellulose, making them particularly relevant for industries seeking to produce renewable fuels from lignocellulosic biomass. The secretome of the cellulolytic model fungus Trichoderma reesei contains two GH7s, termed TrCel7A and TrCel7B. Despite having high structural and sequence similarities, the two enzymes are functionally quite different. TrCel7A is an exolytic, processive cellobiohydrolase (CBH), with high activity on crystalline cellulose, whereas TrCel7B is an endoglucanase (EG) with a preference for more amorphous cellulose. At the structural level, these functional differences are usually ascribed to the flexible loops that cover the substrate-binding areas. TrCel7A has an extensive tunnel created by eight peripheral loops, and the absence of four of these loops in TrCel7B makes its catalytic domain a more open cleft. To investigate the structure-function relationships of these loops, here we produced and kinetically characterized several variants in which four loops unique to TrCel7A were individually deleted to resemble the arrangement in the TrCel7B structure. Analysis of a range of kinetic parameters consistently indicated that the B2 loop, covering the substrate-binding subsites -3 and -4 in TrCel7A, was a key determinant for the difference in CBH- or EG-like behavior between TrCel7A and TrCel7B. Conversely, the B3 and B4 loops, located closer to the catalytic site in TrCel7A, were less important for these activities. We surmise that these results could be useful both in further mechanistic investigations and for guiding engineering efforts of this industrially important enzyme family.

GH16 and GH81 family β-(1,3)-glucanases in Aspergillus fumigatus are essential for conidial cell wall morphogenesis

The fungal cell wall is a rigid structure because of fibrillar and branched β-(1,3)-glucan linked to chitin. Softening of the cell wall is an essential phenomenon during fungal morphogenesis, wherein rigid cell wall structures are cleaved by glycosylhydrolases. During the search for glycosylhydrolases acting on β-(1,3)-glucan, we identified seven genes in the Aspergillus fumigatus genome coding for potential endo-β-(1,3)-glucanase. ENG1 (previously characterized and named ENGL1, Mouyna et al., ), belongs to the Glycoside-Hydrolase 81 (GH81) family, while ENG2 to ENG7, to GH16 family. ENG1 and four GH16 genes (ENG2-5) were expressed in the resting conidia as well as during germination, suggesting an essential role during A. fumigatus morphogenesis. Here, we report the effect of sequential deletion of AfENG2-5 (GH16) followed by AfENG1 (GH81) deletion in the Δeng2,3,4,5 mutant. The Δeng1,2,3,4,5 mutant showed conidial defects, with linear chains of conidia unable to separate while the germination rate was not affected. These results show, for the first time in a filamentous fungus, that endo β-(1,3)-glucanases are essential for proper conidial cell wall assembly and thus segregation of conidia during conidiation.

GNBP domain of Anopheles darlingi: are polymorphic inversions and gene variation related to adaptive evolution?

Anopheles darlingi is the main malaria vector in humans in South America. In the Amazon basin, it lives along the banks of rivers and lakes, which responds to the annual hydrological cycle (dry season and rainy season). In these breeding sites, the larvae of this mosquito feed on decomposing organic and microorganisms, which can be pathogenic and trigger the activation of innate immune system pathways, such as proteins Gram-negative binding protein (GNBP). Such environmental changes affect the occurrence of polymorphic inversions especially at the heterozygote frequency, which confer adaptative advantage compared to homozygous inversions. We mapped the GNBP probe to the An. darlingi 2Rd inversion by fluorescent in situ hybridization (FISH), which was a good indicator of the GNBP immune response related to the chromosomal polymorphic inversions and adaptative evolution. To better understand the evolutionary relations and time of divergence of the GNBP of An. darlingi, we compared it with nine other mosquito GNBPs. The results of the phylogenetic analysis of the GNBP sequence between the species of mosquitoes demonstrated three clades. Clade I and II included the GNBPB5 sequence, and clade III the sequence of GNBPB1. Most of these sequences of GNBP analyzed were homologous with that of subfamily B, including that of An. gambiae (87 %), therefore suggesting that GNBP of An. darling belongs to subfamily B. This work helps us understand the role of inversion polymorphism in evolution of An. darlingi.

Genomics review of holocellulose deconstruction by aspergilli

SUMMARY

Biomass is constructed of dense recalcitrant polymeric materials: proteins, lignin, and holocellulose, a fraction constituting fibrous cellulose wrapped in hemicellulose-pectin. Bacteria and fungi are abundant in soil and forest floors, actively recycling biomass mainly by extracting sugars from holocellulose degradation. Here we review the genome-wide contents of seven Aspergillus species and unravel hundreds of gene models encoding holocellulose-degrading enzymes. Numerous apparent gene duplications followed functional evolution, grouping similar genes into smaller coherent functional families according to specialized structural features, domain organization, biochemical activity, and genus genome distribution. Aspergilli contain about 37 cellulase gene models, clustered in two mechanistic categories: 27 hydrolyze and 10 oxidize glycosidic bonds. Within the oxidative enzymes, we found two cellobiose dehydrogenases that produce oxygen radicals utilized by eight lytic polysaccharide monooxygenases that oxidize glycosidic linkages, breaking crystalline cellulose chains and making them accessible to hydrolytic enzymes. Among the hydrolases, six cellobiohydrolases with a tunnel-like structural fold embrace single crystalline cellulose chains and cooperate at nonreducing or reducing end termini, splitting off cellobiose. Five endoglucanases group into four structural families and interact randomly and internally with cellulose through an open cleft catalytic domain, and finally, seven extracellular β-glucosidases cleave cellobiose and related oligomers into glucose. Aspergilli contain, on average, 30 hemicellulase and 7 accessory gene models, distributed among 9 distinct functional categories: the backbone-attacking enzymes xylanase, mannosidase, arabinase, and xyloglucanase, the short-side-chain-removing enzymes xylan α-1,2-glucuronidase, arabinofuranosidase, and xylosidase, and the accessory enzymes acetyl xylan and feruloyl esterases.

Bioinformatic characterisation of genes encoding cell wall degrading enzymes in the Phytophthora parasitica genome

Background

A critical aspect of plant infection by the majority of pathogens is penetration of the plant cell wall. This process requires the production and secretion of a broad spectrum of pathogen enzymes that target and degrade the many complex polysaccharides in the plant cell wall. As a necessary framework for a study of the expression of cell wall degrading enzymes (CWDEs) produced by the broad host range phytopathogen, Phytophthora parasitica, we have conducted an in-depth bioinformatics analysis of the entire complement of genes encoding CWDEs in this pathogen’s genome.

Results

Our bioinformatic analysis indicates that 431 (2%) of the 20,825 predicted proteins encoded by the P. parasitica genome, are carbohydrate-active enzymes (CAZymes) involved in the degradation of cell wall polysaccharides. Of the 431 proteins, 337 contain classical N-terminal secretion signals and 67 are predicted to be targeted to the non-classical secretion pathway. Identification of CAZyme catalytic activity based on primary protein sequence is difficult, nevertheless, detailed comparisons with previously characterized enzymes has allowed us to determine likely enzyme activities and targeted substrates for many of the P. parasitica CWDEs. Some proteins (12%) contain more than one CAZyme module but, in most cases, multiple modules are from the same CAZyme family. Only 12 P. parasitica CWDEs contain both catalytically-active (glycosyl hydrolase) and non-catalytic (carbohydrate binding) modules, a situation that contrasts with that in fungal phytopathogens. Other striking differences between the complements of CWDEs in P. parasitica and fungal phytopathogens are seen in the CAZyme families that target cellulose, pectins or β-1,3-glucans (e.g. callose). About 25% of P. parasitica CAZymes are solely directed towards pectin degradation, with the majority coming from pectin lyase or carbohydrate esterase families. Fungal phytopathogens typically contain less than half the numbers of these CAZymes. The P. parasitica genome, like that of other Oomycetes, is rich in CAZymes that target β-1,3-glucans.

Conclusions

This detailed analysis of the full complement of P. parasitica cell wall degrading enzymes provides a framework for an in-depth study of patterns of expression of these pathogen genes during plant infection and the induction or repression of expression by selected substrates.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-785) contains supplementary material, which is available to authorized users.

An update of DIVERGE software for functional divergence analysis of protein family

DIVERGE is a software system for phylogeny-based analyses of protein family evolution and functional divergence. It provides a suite of statistical tools for selection and prioritization of the amino acid sites that are responsible for the functional divergence of a gene family. The synergistic efforts of DIVERGE and other methods have convincingly demonstrated that the pattern of rate change at a particular amino acid site may contain insightful information about the underlying functional divergence following gene duplication. These predicted sites may be used as candidates for further experiments. We are now releasing an updated version of DIVERGE with the following improvements: 1) a feasible approach to examining functional divergence in nearly complete sequences by including deletions and insertions (indels); 2) the calculation of the false discovery rate of functionally diverging sites; 3) estimation of the effective number of functional divergence-related sites that is reliable and insensitive to cutoffs; 4) a statistical test for asymmetric functional divergence; and 5) a new method to infer functional divergence specific to a given duplicate cluster. In addition, we have made efforts to improve software design and produce a well-written software manual for the general user.