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.

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.

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.

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.

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.

Comparative genomics analysis of Lactobacillus species associated with weight gain or weight protection.

BACKGROUND: Some Lactobacillus species are associated with obesity and weight gain while others are associated with weight loss. Lactobacillus spp. and bifidobacteria represent a major bacterial population of the small intestine where lipids and simple carbohydrates are absorbed, particularly in the duodenum and jejunum. The objective of this study was to identify Lactobacillus spp. proteins involved in carbohydrate and lipid metabolism associated with weight modifications. METHODS: We examined a total of 13 complete genomes belonging to seven different Lactobacillus spp. previously associated with weight gain or weight protection. We combined the data obtained from the Rapid Annotation using Subsystem Technology, Batch CD-Search and Gene Ontology to classify gene function in each genome. RESULTS: We observed major differences between the two groups of genomes. Weight gain-associated Lactobacillus spp. appear to lack enzymes involved in the catabolism of fructose, defense against oxidative stress and the synthesis of dextrin, L-rhamnose and acetate. Weight protection-associated Lactobacillus spp. encoded a significant gene amount of glucose permease. Regarding lipid metabolism, thiolases were only encoded in the genome of weight gain-associated Lactobacillus spp. In addition, we identified 18 different types of bacteriocins in the studied genomes, and weight gain-associated Lactobacillus spp. encoded more bacteriocins than weight protection-associated Lactobacillus spp. CONCLUSIONS: The results of this study revealed that weight protection-associated Lactobacillus spp. have developed defense mechanisms for enhanced glycolysis and defense against oxidative stress. Weight gain-associated Lactobacillus spp. possess a limited ability to breakdown fructose or glucose and might reduce ileal brake effects.

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.

Glycosyltransferases: structures, functions, and mechanisms.

Glycosyltransferases catalyze glycosidic bond formation using sugar donors containing a nucleoside phosphate or a lipid phosphate leaving group. Only two structural folds, GT-A and GT-B, have been identified for the nucleotide sugar-dependent enzymes, but other folds are now appearing for the soluble domains of lipid phosphosugar-dependent glycosyl transferases. Structural and kinetic studies have provided new insights. Inverting glycosyltransferases utilize a direct displacement S(N)2-like mechanism involving an enzymatic base catalyst. Leaving group departure in GT-A fold enzymes is typically facilitated via a coordinated divalent cation, whereas GT-B fold enzymes instead use positively charged side chains and/or hydroxyls and helix dipoles. The mechanism of retaining glycosyltransferases is less clear. The expected two-step double-displacement mechanism is rendered less likely by the lack of conserved architecture in the region where a catalytic nucleophile would be expected. A mechanism involving a short-lived oxocarbenium ion intermediate now seems the most likely, with the leaving phosphate serving as the base.

Structural enzymology of carbohydrate-active enzymes: implications for the post-genomic era.

Simple and complex carbohydrates have been described as "the last frontier of molecular and cell biology". The enzymes that are required for the synthesis and degradation of these compounds provide an enormous challenge in the post-genomic era. This reflects both the extreme chemical and functional diversity of sugars and the difficulties in characterizing both the substrates and the enzymes themselves. The vast myriad of enzymes involved in the synthesis, modification and degradation of oligosaccharides and polysaccharides is only just being unveiled by genomic sequencing. These so-called "carbohydrate-active enzymes" lend themselves to classification by sensitive sequence similarity detection methods. The modularity, often extremely complex, of these enzymes must first be dissected and annotated before high throughput characterization or "structural genomics" approaches may be employed. Once achieved, modular analysis also permits collation of a detailed "census" of carbohydrate-active enzymes for a whole organism or throughout an ecosystem. At the structural level, improvements in X-ray crystallography have opened up a three-dimensional understanding of the way these enzymes work. The mechanisms of many of the glycoside hydrolase families are becoming clearer, yet glycosyltransferases are only slowly revealing their secrets. What is clear from the genomic and structural data is that if we are to harness the latent power of glycogenomics, scientists must consider distant sequence relatives revealed by the sequence families or other sensitive detection methods.

The location of the ligand-binding site of carbohydrate-binding modules that have evolved from a common sequence is not conserved.

Polysaccharide-degrading enzymes are generally modular proteins that contain non-catalytic carbohydrate-binding modules (CBMs), which potentiate the activity of the catalytic module. CBMs have been grouped into sequence-based families, and three-dimensional structural data are available for half of these families. Clostridium thermocellum xylanase 11A is a modular enzyme that contains a CBM from family 6 (CBM6), for which no structural data are available. We have determined the crystal structure of this module to a resolution of 2.1 A. The protein is a beta-sandwich that contains two potential ligand-binding clefts designated cleft A and B. The CBM interacts primarily with xylan, and NMR spectroscopy coupled with site-directed mutagenesis identified cleft A, containing Trp-92, Tyr-34, and Asn-120, as the ligand-binding site. The overall fold of CBM6 is similar to proteins in CBM families 4 and 22, although surprisingly the ligand-binding site in CBM4 and CBM22 is equivalent to cleft B in CBM6. These structural data define a superfamily of CBMs, comprising CBM4, CBM6, and CBM22, and demonstrate that, although CBMs have evolved from a relatively small number of ancestors, the structural elements involved in ligand recognition have been assembled at different locations on the ancestral scaffold.

Glycosidase active site mutations in human alpha-L-iduronidase.

Mucopolysaccharidosis type I (MPS I; McKusick 25280) results from a deficiency in alpha-L-iduronidase activity. Using a bioinformatics approach, we have previously predicted the putative acid/base catalyst and nucleophile residues in the active site of this human lysosomal glycosidase to be Glu182 and Glu299, respectively. To obtain experimental evidence supporting these predictions, wild-type alpha-L-iduronidase and site-directed mutants E182A and E299A were individually expressed in Chinese hamster ovary-K1 cell lines. We have compared the synthesis, processing, and catalytic properties of the two mutant proteins with wild-type human alpha-L-iduronidase. Both E182A and E299A transfected cells produced catalytically inactive human alpha-L-iduronidase protein at levels comparable to the wild-type control. The E182A protein was synthesized, processed, targeted to the lysosome, and secreted in a similar fashion to wild-type alpha-L-iduronidase. The E299A mutant protein was also synthesized and secreted similarly to the wild-type enzyme, but there were alterations in its rate of traffic and proteolytic processing. These data indicate that the enzymatic inactivity of the E182A and E299A mutants is not due to problems of synthesis/folding, but to the removal of key catalytic residues. In addition, we have identified a MPS I patient with an E182K mutant allele. The E182K mutant protein was expressed in CHO-K1 cells and also found to be enzymatically inactive. Together, these results support the predicted role of E182 and E299 in the catalytic mechanism of alpha-L-iduronidase and we propose that the mutation of either of these residues would contribute to a very severe clinical phenotype in a MPS I patient.

A census of carbohydrate-active enzymes in the genome of Arabidopsis thaliana.

The synthesis, modification, and breakdown of carbohydrates is one of the most fundamentally important reactions in nature. The structural and functional diversity of glycosides is mirrored by a vast array of enzymes involved in their synthesis (glycosyltransferases), modification (carbohydrate esterases) and breakdown (glycoside hydrolases and polysaccharide lyases). The importance of these processes is reflected in the dedication of 1-2% of an organism's genes to glycoside hydrolases and glycosyltransferases alone. In plants, these processes are of particular importance for cell-wall synthesis and expansion. starch metabolism, defence against pathogens, symbiosis and signalling. Here we present an analysis of over 730 open reading frames representing the two main classes of carbohydrate-active enzymes, glycoside hydrolases and glycosyltransferases, in the genome of Arabidopsis thaliana. The vast importance of these enzymes in cell-wall formation and degradation is revealed along with the unexpected dominance of pectin degradation in Arabidopsis, with at least 170 open-reading frames dedicated solely to this task.

The iota-carrageenase of Alteromonas fortis. A beta-helix fold-containing enzyme for the degradation of a highly polyanionic polysaccharide.

Carrageenans are gel-forming hydrocolloids extracted from the cell walls of marine red algae. They consist of d-galactose residues bound by alternate alpha(1-->3) and beta(1-->4) linkages and substituted by one (kappa-carrageenan), two (iota-carrageenan), or three (lambda-carrageenan) sulfate-ester groups per disaccharide repeating unit. Both the kappa- and iota-carrageenan chains adopt ordered conformations leading to the formation of highly ordered aggregates of double-stranded helices. Several kappa-carrageenases and iota-carrageenases have been cloned from marine bacteria. Kappa-carrageenases belong to family 16 of the glycoside hydrolases, which essentially encompasses polysaccharidases specialized in the hydrolysis of the neutral polysaccharides such as agarose, laminarin, lichenan, and xyloglucan. In contrast, iota-carrageenases constitute a novel glycoside hydrolase structural family. We report here the crystal structure of Alteromonas fortis iota-carrageenase at 1.6 A resolution. The enzyme folds into a right-handed parallel beta-helix of 10 complete turns with two additional C-terminal domains. Glu(245), Asp(247), or Glu(310), in the cleft of the enzyme, are proposed as candidate catalytic residues. The protein contains one sodium and one chloride binding site and three calcium binding sites shown to be involved in stabilizing the enzyme structure.

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.

Enzyme-catalyzed condensation reaction in a mammalian alpha-amylase. High-resolution structural analysis of an enzyme-inhibitor complex.

Mammalian alpha-amylases catalyze the hydrolysis of alpha-linked glucose polymers according to a complex processive mechanism. We have determined the X-ray structures of porcine pancreatic alpha-amylase complexes with the smallest molecule of the trestatin family (acarviosine-glucose) which inhibits porcine pancreatic alpha-amylase and yet is not hydrolyzed by the enzyme. A structure analysis at 1.38 A resolution of this complex allowed for a clear identification of a genuine single hexasaccharide species composed of two alpha-1,4-linked original molecules bound to the active site of the enzyme. The structural results supported by mass spectrometry experiments provide evidence for an enzymatically catalyzed condensation reaction in the crystal.

Streptomyces matensis laminaripentaose hydrolase is an 'inverting' beta-1,3-glucanase.

The laminaripentaose-producing beta-1,3-glucanase of Streptomyces matensis is a member of the glycoside hydrolase family GH-64. We have constructed and purified a recombinant hexahistidine-tagged form of the enzyme for characterisation. The enzyme, which exists as a monomer in solution, hydrolyses beta-1,3-glucan by a mechanism leading to overall inversion of the anomeric configuration. This is the first determination of the mechanism prevailing in glycoside hydrolase family GH-64 and this is the first characterisation of an 'inverting' beta-1,3-glucanase.

The celA gene, encoding a glycosyl hydrolase family 3 beta-glucosidase in Azospirillum irakense, is required for optimal growth on cellobiosides.

The CelA beta-glucosidase of Azospirillum irakense, belonging to glycosyl hydrolase family 3 (GHF3), preferentially hydrolyzes cellobiose and releases glucose units from the C(3), C(4), and C(5) oligosaccharides. The growth of a DeltacelA mutant on these cellobiosides was affected. In A. irakense, the GHF3 beta-glucosidases appear to be functional alternatives for the GHF1 beta-glucosidases in the assimilation of beta-glucosides by other bacteria.