101 Dothideomycetes genomes: A test case for predicting lifestyles and emergence of pathogens.

Dothideomycetes is the largest class of kingdom Fungi and comprises an incredible diversity of lifestyles, many of which have evolved multiple times. Plant pathogens represent a major ecological niche of the class Dothideomycetes and they are known to infect most major food crops and feedstocks for biomass and biofuel production. Studying the ecology and evolution of Dothideomycetes has significant implications for our fundamental understanding of fungal evolution, their adaptation to stress and host specificity, and practical implications with regard to the effects of climate change and on the food, feed, and livestock elements of the agro-economy. In this study, we present the first large-scale, whole-genome comparison of 101 Dothideomycetes introducing 55 newly sequenced species. The availability of whole-genome data produced a high-confidence phylogeny leading to reclassification of 25 organisms, provided a clearer picture of the relationships among the various families, and indicated that pathogenicity evolved multiple times within this class. We also identified gene family expansions and contractions across the Dothideomycetes phylogeny linked to ecological niches providing insights into genome evolution and adaptation across this group. Using machine-learning methods we classified fungi into lifestyle classes with >95 % accuracy and identified a small number of gene families that positively correlated with these distinctions. This can become a valuable tool for genome-based prediction of species lifestyle, especially for rarely seen and poorly studied species.

Comparative genomics of Mortierella elongata and its bacterial endosymbiont Mycoavidus cysteinexigens.

Endosymbiosis of bacteria by eukaryotes is a defining feature of cellular evolution. In addition to well-known bacterial origins for mitochondria and chloroplasts, multiple origins of bacterial endosymbiosis are known within the cells of diverse animals, plants and fungi. Early-diverging lineages of terrestrial fungi harbor endosymbiotic bacteria belonging to the Burkholderiaceae. We sequenced the metagenome of the soil-inhabiting fungus Mortierella elongata and assembled the complete circular chromosome of its endosymbiont, Mycoavidus cysteinexigens, which we place within a lineage of endofungal symbionts that are sister clade to Burkholderia. The genome of M. elongata strain AG77 features a core set of primary metabolic pathways for degradation of simple carbohydrates and lipid biosynthesis, while the M. cysteinexigens (AG77) genome is reduced in size and function. Experiments using antibiotics to cure the endobacterium from the host demonstrate that the fungal host metabolism is highly modulated by presence/absence of M. cysteinexigens. Independent comparative phylogenomic analyses of fungal and bacterial genomes are consistent with an ancient origin for M. elongata - M. cysteinexigens symbiosis, most likely over 350 million years ago and concomitant with the terrestrialization of Earth and diversification of land fungi and plants.

Exploring the genomic diversity of black yeasts and relatives (Chaetothyriales, Ascomycota).

The order Chaetothyriales (Pezizomycotina, Ascomycetes) harbours obligatorily melanised fungi and includes numerous etiologic agents of chromoblastomycosis, phaeohyphomycosis and other diseases of vertebrate hosts. Diseases range from mild cutaneous to fatal cerebral or disseminated infections and affect humans and cold-blooded animals globally. In addition, Chaetothyriales comprise species with aquatic, rock-inhabiting, ant-associated, and mycoparasitic life-styles, as well as species that tolerate toxic compounds, suggesting a high degree of versatile extremotolerance. To understand their biology and divergent niche occupation, we sequenced and annotated a set of 23 genomes of main the human opportunists within the Chaetothyriales as well as related environmental species. Our analyses included fungi with diverse life-styles, namely opportunistic pathogens and closely related saprobes, to identify genomic adaptations related to pathogenesis. Furthermore, ecological preferences of Chaetothyriales were analysed, in conjuncture with the order-level phylogeny based on conserved ribosomal genes. General characteristics, phylogenomic relationships, transposable elements, sex-related genes, protein family evolution, genes related to protein degradation (MEROPS), carbohydrate-active enzymes (CAZymes), melanin synthesis and secondary metabolism were investigated and compared between species. Genome assemblies varied from 25.81 Mb (Capronia coronata) to 43.03 Mb (Cladophialophora immunda). The bantiana-clade contained the highest number of predicted genes (12 817 on average) as well as larger genomes. We found a low content of mobile elements, with DNA transposons from Tc1/Mariner superfamily being the most abundant across analysed species. Additionally, we identified a reduction of carbohydrate degrading enzymes, specifically many of the Glycosyl Hydrolase (GH) class, while most of the Pectin Lyase (PL) genes were lost in etiological agents of chromoblastomycosis and phaeohyphomycosis. An expansion was found in protein degrading peptidase enzyme families S12 (serine-type D-Ala-D-Ala carboxypeptidases) and M38 (isoaspartyl dipeptidases). Based on genomic information, a wide range of abilities of melanin biosynthesis was revealed; genes related to metabolically distinct DHN, DOPA and pyomelanin pathways were identified. The MAT (MAting Type) locus and other sex-related genes were recognized in all 23 black fungi. Members of the asexual genera Fonsecaea and Cladophialophora appear to be heterothallic with a single copy of either MAT-1-1 or MAT-1-2 in each individual. All Capronia species are homothallic as both MAT1-1 and MAT1-2 genes were found in each single genome. The genomic synteny of the MAT-locus flanking genes (SLA2-APN2-COX13) is not conserved in black fungi as is commonly observed in Eurotiomycetes, indicating a unique genomic context for MAT in those species. The heterokaryon (het) genes expansion associated with the low selective pressure at the MAT-locus suggests that a parasexual cycle may play an important role in generating diversity among those fungi.

Large-scale genomic analyses with machine learning uncover predictive patterns associated with fungal phytopathogenic lifestyles and traits.

Invasive plant pathogenic fungi have a global impact, with devastating economic and environmental effects on crops and forests. Biosurveillance, a critical component of threat mitigation, requires risk prediction based on fungal lifestyles and traits. Recent studies have revealed distinct genomic patterns associated with specific groups of plant pathogenic fungi. We sought to establish whether these phytopathogenic genomic patterns hold across diverse taxonomic and ecological groups from the Ascomycota and Basidiomycota, and furthermore, if those patterns can be used in a predictive capacity for biosurveillance. Using a supervised machine learning approach that integrates phylogenetic and genomic data, we analyzed 387 fungal genomes to test a proof-of-concept for the use of genomic signatures in predicting fungal phytopathogenic lifestyles and traits during biosurveillance activities. Our machine learning feature sets were derived from genome annotation data of carbohydrate-active enzymes (CAZymes), peptidases, secondary metabolite clusters (SMCs), transporters, and transcription factors. We found that machine learning could successfully predict fungal lifestyles and traits across taxonomic groups, with the best predictive performance coming from feature sets comprising CAZyme, peptidase, and SMC data. While phylogeny was an important component in most predictions, the inclusion of genomic data improved prediction performance for every lifestyle and trait tested. Plant pathogenicity was one of the best-predicted traits, showing the promise of predictive genomics for biosurveillance applications. Furthermore, our machine learning approach revealed expansions in the number of genes from specific CAZyme and peptidase families in the genomes of plant pathogens compared to non-phytopathogenic genomes (saprotrophs, endo- and ectomycorrhizal fungi). Such genomic feature profiles give insight into the evolution of fungal phytopathogenicity and could be useful to predict the risks of unknown fungi in future biosurveillance activities.

Functional analysis of family GH36 α-galactosidases from Ruminococcus gnavus E1: insights into the metabolism of a plant oligosaccharide by a human gut symbiont.

Ruminococcus gnavus belongs to the 57 most common species present in 90% of individuals. Previously, we identified an α-galactosidase (Aga1) belonging to glycoside hydrolase (GH) family 36 from R. gnavus E1 (M. Aguilera, H. Rakotoarivonina, A. Brutus, T. Giardina, G. Simon, and M. Fons, Res. Microbiol. 163:14-21, 2012). Here, we identified a novel GH36-encoding gene from the same strain and termed it aga2. Although aga1 showed a very simple genetic organization, aga2 is part of an operon of unique structure, including genes putatively encoding a regulator, a GH13, two phosphotransferase system (PTS) sequences, and a GH32, probably involved in extracellular and intracellular sucrose assimilation. The 727-amino-acid (aa) deduced Aga2 protein shares approximately 45% identity with Aga1. Both Aga1 and Aga2 expressed in Escherichia coli showed strict specificity for α-linked galactose. Both enzymes were active on natural substrates such as melibiose, raffinose, and stachyose. Aga1 and Aga2 occurred as homotetramers in solution, as shown by analytical ultracentrifugation. Modeling of Aga1 and Aga2 identified key amino acids which may be involved in substrate specificity and stabilization of the α-linked galactoside substrates within the active site. Furthermore, Aga1 and Aga2 were both able to perform transglycosylation reactions with α-(1,6) regioselectivity, leading to the formation of product structures up to [Hex](12) and [Hex](8), respectively. We suggest that Aga1 and Aga2 play essential roles in the metabolism of dietary oligosaccharides and could be used for the design of galacto-oligosaccharide (GOS) prebiotics, known to selectively modulate the beneficial gut microbiota.

Carbohydrate-active enzymes revealed in Coptotermes formosanus (Isoptera: Rhinotermitidae) transcriptome.

Coptotermes formosanus is one of the most destructive wood-feeding termites. To understand the molecular mechanisms that regulate the development of the termite, a normalized C. formosanus cDNA library was constructed using mixed RNA isolated from workers, soldiers, nymphs and alates of both sexes. The sequencing of this library generated 131 636 expressed sequence tags (ESTs) and 25 939 assembled unigenes. The carbohydrate-active enzymes (CAZymes) revealed in this library were analysed in the present report. A total of 509 putative CAZymes were identified. Diverse cellulolytic enzymes were uncovered from both the host termite and from symbionts harboured by the termite, which were possibly the result of the high efficiency of cellulose utilization. CAZymes associated with trehalose biosynthetic and metabolic pathways were also identified, which are potential regulators of the physiological activities of trehalose, an important insect blood sugar. Representative CAZyme coding genes in glycoside hydrolase family 1 (GH1) were quantitatively analysed. The results showed that the five GH1 ß-glucosidase genes were expressed differentially among different castes and one of them was female alate-specific. Overall, the normalized EST library provides a comprehensive genetic resource of C. formosanus and will serve a diverse range of research areas. The CAZymes represent one of the repositories of enzymes useful for physiological studies and applications in sugar-based biofuel production.

Glycoside hydrolases and glycosyltransferases: families and functional modules.

The past year has witnessed the expected increase in the number of solved structures of glycoside hydrolases and glycosyltransferases, and their constitutive modules. These structures show that, while glycoside hydrolases display an extraordinary variety of folds, glycosyltransferases and carbohydrate-binding modules appear to belong to a much smaller number of folding families.

Structural and sequence-based classification of glycoside hydrolases.

The diversity of oligo- and polysaccharides provides an abundance of biological roles for these carbohydrates. The enzymes hydrolysing these compounds, the glycoside hydrolases, therefore mediate a wealth of biological functions. Glycoside hydrolases fall into a number of sequence-based families. The recent analysis of these families, coupled with the burgeoning number of 3D structures, provides a detailed insight into the structure, function and catalytic mechanism of these enzymes.

Structural and electronic determinants of lytic polysaccharide monooxygenase reactivity on polysaccharide substrates.

Lytic polysaccharide monooxygenases (LPMOs) are industrially important copper-dependent enzymes that oxidatively cleave polysaccharides. Here we present a functional and structural characterization of two closely related AA9-family LPMOs from Lentinus similis (LsAA9A) and Collariella virescens (CvAA9A). LsAA9A and CvAA9A cleave a range of polysaccharides, including cellulose, xyloglucan, mixed-linkage glucan and glucomannan. LsAA9A additionally cleaves isolated xylan substrates. The structures of CvAA9A and of LsAA9A bound to cellulosic and non-cellulosic oligosaccharides provide insight into the molecular determinants of their specificity. Spectroscopic measurements reveal differences in copper co-ordination upon the binding of xylan and glucans. LsAA9A activity is less sensitive to the reducing agent potential when cleaving xylan, suggesting that distinct catalytic mechanisms exist for xylan and glucan cleavage. Overall, these data show that AA9 LPMOs can display different apparent substrate specificities dependent upon both productive protein-carbohydrate interactions across a binding surface and also electronic considerations at the copper active site.

Structures and mechanisms of glycosyl hydrolases.

The wealth of information provided by the recent structure determinations of many different glycosyl hydrolases shows that the substrate specificity and the mode of action of these enzymes are governed by exquisite details of their three-dimensional structures rather than by their global fold.

Influenza B virus neuraminidase can synthesize its own inhibitor.

BACKGROUND: Neuraminidase, one of the two surface glycoproteins of influenza virus, cleaves terminal sialic acid residues from glycolipids or glycoproteins. Its crystal structure is known at high resolution, but the mechanism of glycosyl hydrolysis remains unclear. RESULTS: We have determined the crystal structure at 1.8 A resolution of two complexes of influenza B/Beijing neuraminidase containing either the reaction product, sialic acid, or the transition state analogue inhibitor, 2,3-dehydro-2-deoxy-N-acetylneuraminic acid (DANA). The sialic acid is bound in a distorted 'boat' conformation closely resembling that of bound DANA, stabilized by a conserved tyrosine residue (Tyr408). This distortion also gives rise to a suicidal side reaction that converts sialic acid to DANA at a low rate. CONCLUSIONS: The mechanism of neuraminidase action is distinct from that of other known glycosyl hydrolases. Substrate distortion appears to be the driving force in glycosyl bond hydrolysis and the proton required for catalysis can probably be donated by water, rather than by residues in the active site, thus allowing the enzyme to operate at high pH. The side reaction converting sialic acid to DANA appears reasonably favourable, and it is unclear how this is minimized by the enzyme.

The crystal structures of Sinapis alba myrosinase and a covalent glycosyl-enzyme intermediate provide insights into the substrate recognition and active-site machinery of an S-glycosidase.

BACKGROUND: Myrosinase is the enzyme responsible for the hydrolysis of a variety of plant anionic 1-thio-beta-D-glucosides called glucosinolates. Myrosinase and glucosinolates, which are stored in different tissues of the plant, are mixed during mastication generating toxic by-products that are believed to play a role in the plant defence system. Whilst O-glycosidases are extremely widespread in nature, myrosinase is the only known S-glycosidase. This intriguing enzyme, which shows sequence similarities with O-glycosidases, offers the opportunity to analyze the similarities and differences between enzymes hydrolyzing S- and O-glycosidic bonds. RESULTS: The structures of native myrosinase from white mustard seed (Sinapis alba) and of a stable glycosyl-enzyme intermediate have been solved at 1.6 A resolution. The protein folds into a (beta/alpha)8-barrel structure, very similar to that of the cyanogenic beta-glucosidase from white clover. The enzyme forms a dimer stabilized by a Zn2+ ion and is heavily glycosylated. At one glycosylation site the complete structure of a plant-specific heptasaccharide is observed. The myrosinase structure reveals a hydrophobic pocket, ideally situated for the binding of the hydrophobic sidechain of glucosinolates, and two arginine residues positioned for interaction with the sulphate group of the substrate. With the exception of the replacement of the general acid/base glutamate by a glutamine residue, the catalytic machinery of myrosinase is identical to that of the cyanogenic beta-glucosidase. The structure of the glycosyl-enzyme intermediate shows that the sugar ring is bound via an alpha-glycosidic linkage to Glu409, the catalytic nucleophile of myrosinase. CONCLUSIONS: The structure of myrosinase shows features which illustrate the adaptation of the plant enzyme to the dehydrated environment of the seed. The catalytic mechanism of myrosinase is explained by the excellent leaving group properties of the substrate aglycons, which do not require the assistance of an enzymatic acid catalyst. The replacement of the general acid/base glutamate of O-glycosidases by a glutamine residue in myrosinase suggests that for hydrolysis of the glycosyl-enzyme, the role of this residue is to ensure a precise positioning of a water molecule rather than to provide general base assistance.

The kappa-carrageenase of P. carrageenovora features a tunnel-shaped active site: a novel insight in the evolution of Clan-B glycoside hydrolases.

BACKGROUND: kappa-carrageenans are gel-forming, sulfated 1,3-alpha-1,4-beta-galactans from the cell walls of marine red algae. The kappa-carrageenase from the marine, gram-negative bacterium Pseudoalteromonas carrageenovora degrades kappa-carrageenan both in solution and in solid state by an endoprocessive mechanism. This beta-galactanase belongs to the clan-B of glycoside hydrolases. RESULTS: The structure of P. carrageenovora kappa-carrageenase has been solved to 1.54 A resolution by the multiwavelength anomalous diffraction (MAD) method, using a seleno-methionine-substituted form of the enzyme. The enzyme folds into a curved beta sandwich, with a tunnel-like active site cavity. Another remarkable characteristic is the presence of an arginine residue at subsite -1. CONCLUSIONS: The crystal structure of P. carrageenovora kappa-carrageenase is the first three-dimensional structure of a carrageenase. Its tunnel-shaped active site, the first to be reported for enzymes other than cellulases, suggests that such tunnels are associated with the degradation of solid polysaccharides. Clan-B glycoside hydrolases fall into two subgroups, one with catalytic machinery held by an ancestral beta bulge, and the other in which it is held by a regular beta strand. At subsite -1, all of these hydrolases exhibit an aromatic amino acid that interacts with the hexopyranose ring of the monosaccharide undergoing catalysis. In addition, in kappa-carrageenases, an arginine residue recognizes the sulfate-ester substituents of the beta-linked kappa-carrageenan monomers. It also appears that, in addition to the nucleophile and acid/base catalysts, two other amino acids are involved with the catalytic cycle, accelerating the deglycosylation step.

Family A cellulases: two essential tryptophan residues in endoglucanase III from Trichoderma reesei.

Three tryptophan residues are readily oxidised by N-bromosuccinimide in endoglucanase III from Trichoderma reesei. Evidence was obtained that the residue first modified is situated in the cellulose-binding domain and the second in the enzyme's catalytic site. The latter influences the binding and hydrolysis of soluble substrates. The modification of a third residue does not further affect the catalytic properties. The present results complement published data concerning other identified catalytic residues, and help to clarify the active site structure of family A cellulases.

Crystal structure of a cohesin module from Clostridium cellulolyticum: implications for dockerin recognition.

In the assembly of the Clostridium cellulolyticum cellulosome, the multiple cohesin modules of the scaffolding protein CipC serve as receptors for cellulolytic enzymes which bear a dockerin module. The X-ray structure of a type I C. cellulolyticum cohesin module (Cc-cohesin) has been solved using molecular replacement, and refined at 2.0 A resolution. Despite a rather low sequence identity of 32 %, this module has a fold close to those of the two Clostridium thermocellum cohesin (Ct-cohesin) modules whose 3D structures have been determined previously. Cc-cohesin forms a dimer in the crystal, as do the two Ct-cohesins. We show here that the dimer exists in solution and that addition of dockerin-containing proteins dissociates the dimer. This suggests that the dimerization interface and the cohesin/dockerin interface may overlap. The nature of the overall surface and of the dimer interface of Cc-cohesin differ notably from those of the Ct-cohesin modules, being much less polar, and this may explain the species specificity observed in the cohesin/dockerin interaction of C. cellulolyticum and C. thermocellum. We have produced a topology model of a C. cellulolyticum dockerin and of a Cc-cohesin/dockerin complex using homology modeling and available biochemical data. Our model suggests that a special residue pair, already identified in dockerin sequences, is located at the center of the cohesin surface putatively interacting with the dockerin.

Solution structure of the module X2 1 of unknown function of the cellulosomal scaffolding protein CipC of Clostridium cellulolyticum.

Multidimensional, homo- and heteronuclear magnetic resonance spectroscopy combined with dynamical annealing has been used to determine the structure of a 94 residue module (X2 1) of the scaffolding protein CipC from the anaerobic bacterium Clostridium cellulolyticum. An experimental data set comprising 1647 nuclear Overhauser effect-derived restraints, 105 hydrogen bond restraints and 66 phi torsion angle restraints was used to calculate 20 converging final solutions. The calculated structures have an average rmsd about the mean structure of 0.55(+/-0.11) A for backbone atoms and 1.40(+/-0.11) A for all heavy atoms when fitted over the secondary structural elements. The X2 1 module has an immunoglobulin-like fold with two beta-sheets packed against each other. One sheet contains three strands, the second contains four strands. An additional strand is intercalated between the beta-sandwich, as well as two turns of a 3(.10) helix. X2 1 has a surprising conformational stability and may act as a conformational linker and solubility enhancer within the scaffolding protein. The fold of X2 1 is very similar to that of telokin, titin Ig domain, hemolin D2 domain, twitchin immunoglobulin domain and the first four domains of the IgSF portion of transmembrane cell adhesion molecule. As a consequence, the X2 1 module is the first prokaryotic member assigned to the I set of the immunoglobulin superfamily even though no sequence similarity with any member of this superfamily could be detected.

Specificity mapping of cellulolytic enzymes: classification into families of structurally related proteins confirmed by biochemical analysis.

The specificities of 15 cellulolytic enzymes have been examined using chromophoric glycosides derived from D-glucose, cellobiose, higher cellooligosaccharides, lactose, D-xylose, and beta-(1,4)-xylobiose. Coinciding with a classification based on hydrophobic cluster analysis of amino acid sequences, six families each showing a characteristic specificity pattern were observed. Furthermore, in these cases where the anomeric forms of reaction products were determined, results seem to indicate conservation of intrinsic reaction mechanism (single or double displacement) within each family. On the other hand, the low molecular weight substrates do not discriminate exo- from endocellulases. This functional differentiation is speculated to originate from the presence, in exoenzymes, of a tunnel-shaped active site formed by extra loops in their structure.

Cloning and characterization of a xylanase gene from corn strains of Erwinia chrysanthemi.

The gene encoding a 42-kDa endoxylanase was cloned from Erwinia chrysanthemi strain D1. Sequencing of this gene, called xynA, showed that it encoded a primary protein product of 413 amino acids with an unusual and long (31 amino acid) leader peptide that was cleaved during secretion to the bacterial periplasm. This protein is distinct from xylanases in glycohydrolase families 10 and 11 and, instead, appears to be intermediate between families 5 and 30. The xynA gene is located downstream from a gene with high homology to ATP-dependent RNA helicases and the Escherichia coli recD gene. Large amounts of the mature xylanase were produced by E. coli cells carrying a T7 expression plasmid construct and the protein was isolated from the bacterial periplasmic fraction by chromatography on a CM Bio-gel column. Marker exchange mutagenesis of the xynA gene eliminated the ability of strain D1 to produce detectable extracellular xylanase activity but did not affect virulence on corn leaves.

Structure-based multiple alignment of extracellular pectate lyase sequences.

Pectate lyases are secreted virulence factors which degrade the pectate component of plant cell walls. The evolutionary-based multiple alignment of extracellular pectate lyases has been corrected using three-dimensional structural information derived from Erwinia chyrsanthemi pectate lyases C and E. The new multiple alignment reveals invariant amino acids likely to be involved in two different enzymatic functions.

Hidden domains and active site residues in beta-glycanase-encoding gene sequences?

Reading-frame corrective shifts in the nucleotide sequence upstream, within, or downstream from the putative coding region of several beta-glycanase-encoding genes reported in the literature reveal hidden active-site residues or even additional domains, including a cellulose-binding domain on a beta-mannanase-encoding gene. These findings also help in assigning, to cellulase family A, two enzymes previously found to lack sequence similarity with known cellulase families.