Curation of characterized glycoside hydrolases of fungal origin

Whole-genome sequencing of Ganoderma boninense, the causal agent of basal stem rot disease in oil palm, via combined short- and long-read sequencing

The hemibiotrophic Basidiomycete pathogen Ganoderma boninense (Gb) is the dominant causal agent of oil palm basal stem rot disease. Here, we report a complete chromosomal genome map of Gb using a combination of short-read Illumina and long-read Pacific Biosciences (PacBio) sequencing platforms combined with chromatin conformation capture data from the Chicago and Hi-C platforms. The genome was 55.87 Mb in length and assembled to a high contiguity (N50: 304.34 kb) of 12 chromosomes built from 112 scaffolds, with a total of only 4.34 Mb (~ 7.77%) remaining unplaced. The final assemblies were evaluated for completeness of the genome by using Benchmarking Universal Single Copy Orthologs (BUSCO) v4.1.4, and based on 4464 total BUSCO polyporales group searches, the assemblies yielded 4264 (95.52%) of the conserved orthologs as complete and only a few fragmented BUSCO of 42 (0.94%) as well as a missing BUSCO of 158 (3.53%). Genome annotation predicted a total of 21,074 coding genes, with a GC content ratio of 59.2%. The genome features were analyzed with different databases, which revealed 2471 Gene Ontology/GO (11.72%), 5418 KEGG (Kyoto Encyclopedia of Genes and Genomes) Orthologous/KO (25.71%), 13,913 Cluster of Orthologous Groups of proteins/COG (66.02%), 60 ABC transporter (0.28%), 1049 Carbohydrate-Active Enzymes/CAZy (4.98%), 4005 pathogen-host interactions/PHI (19%), and 515 fungal transcription factor/FTFD (2.44%) genes. The results obtained in this study provide deep insight for further studies in the future.

Different techniques reveal the difference of community structure and function of fungi from root and rhizosphere of Salvia miltiorrhiza Bunge

Fungi have essential functions in plant health and performance. However, the plant-associated functions of many cultured fungi have not been established in detail. Here, the fungal species diversity in Salvia miltiorrhiza roots and rhizosphere was assessed for the first time using culturomics and high-throughput sequencing. We present a comprehensive functional metagenomic analysis of these fungi and verified activity of cellulase and chitinase predicted in the metagenomic analysis. We first collected and cultured fungi from the root and rhizosphere of S. miltiorrhiza. We found 92 species across 37 families and five phyla, with Ascomycota being dominant. Many rDNA internal transcribed spacer sequences could not be assigned to lower taxonomic levels. There were 19 genera of endophytic fungi and 37 genera of rhizosphere fungi. The culturomics approach had lower taxonomic diversity than high-throughput sequencing, but some fungi were only found in cultures. Structural analyses indicated that the dominant species differed in cultured and non-cultured samples at other levels, apart from the phylum level. Functional analysis mapped 223 carbohydrate enzyme families and 393 pathways in the CAZy and KEGG databases, respectively. The most abundant families were glycoside hydrolases and those involved in carbohydrate metabolism. As predicted by metagenomics, we experimentally verified cellulase and chitinase activity for 29 and 74 fungi, respectively. We provide the first evidence of biomass recycling by fungi that are associated with plants. Culturing is essential to reveal the hidden microbial community and critical functions in plant-microbe interactions.

High-quality RNA extraction and the regulation of genes encoding cellulosomes are correlated with growth stage in anaerobic fungi

Anaerobic fungi produce biomass-degrading enzymes and natural products that are important to harness for several biotechnology applications. Although progress has been made in the development of methods for extracting nucleic acids for genomic and transcriptomic sequencing of these fungi, most studies are limited in that they do not sample multiple fungal growth phases in batch culture. In this study, we establish a method to harvest RNA from fungal monocultures and fungal-methanogen co-cultures, and also determine an optimal time frame for high-quality RNA extraction from anaerobic fungi. Based on RNA quality and quantity targets, the optimal time frame in which to harvest anaerobic fungal monocultures and fungal-methanogen co-cultures for RNA extraction was 2-5 days of growth post-inoculation. When grown on cellulose, the fungal strain Anaeromyces robustus cocultivated with the methanogen Methanobacterium bryantii upregulated genes encoding fungal carbohydrate-active enzymes and other cellulosome components relative to fungal monocultures during this time frame, but expression patterns changed at 24-hour intervals throughout the fungal growth phase. These results demonstrate the importance of establishing methods to extract high-quality RNA from anaerobic fungi at multiple time points during batch cultivation.

Preassembled Cas9 Ribonucleoprotein-Mediated Gene Deletion Identifies the Carbon Catabolite Repressor and Its Target Genes in Coprinopsis cinerea

Cre1 is an important transcription factor that regulates carbon catabolite repression (CCR) and is widely conserved across fungi. The cre1 gene has been extensively studied in several Ascomycota species, whereas its role in gene expression regulation in the Basidiomycota species remains poorly understood. Here, we identified and investigated the role of cre1 in Coprinopsis cinerea, a basidiomycete model mushroom that can efficiently degrade lignocellulosic plant wastes. We used a rapid and efficient gene deletion approach based on PCR-amplified split-marker DNA cassettes together with in vitro assembled Cas9-guide RNA ribonucleoproteins (Cas9 RNPs) to generate C. cinerea cre1 gene deletion strains. Gene expression profiling of two independent C. cinerea cre1 mutants showed significant deregulation of carbohydrate metabolism, plant cell wall degrading enzymes (PCWDEs), plasma membrane transporter-related and several transcription factor-encoding genes, among others. Our results support the notion that, like reports in the ascomycetes, Cre1 of C. cinerea orchestrates CCR through a combined regulation of diverse genes, including PCWDEs, transcription factors that positively regulate PCWDEs, and membrane transporters which could import simple sugars that can induce the expression of PWCDEs. Somewhat paradoxically, though in accordance with other Agaricomycetes, genes related to lignin degradation were mostly downregulated in cre1 mutants, indicating they fall under different regulation than other PCWDEs. The gene deletion approach and the data presented here will expand our knowledge of CCR in the Basidiomycota and provide functional hypotheses on genes related to plant biomass degradation. IMPORTANCE Mushroom-forming fungi include some of the most efficient lignocellulosic plant biomass degraders. They degrade dead plant materials by a battery of lignin-, cellulose-, hemicellulose-, and pectin-degrading enzymes, the encoding genes of which are under tight transcriptional control. One of the highest-level regulations of these metabolic enzymes is known as carbon catabolite repression, which is orchestrated by the transcription factor Cre1, and ensures that costly lignocellulose-degrading enzyme genes are expressed only when simple carbon sources (e.g., glucose) are not available. Here, we identified the Cre1 ortholog in a litter decomposer Agaricomycete, Coprinopsis cinerea, knocked it out, and characterized transcriptional changes in the mutants. We identified several dozen lignocellulolytic enzyme genes as well as membrane transporters and other transcription factors as putative target genes of C. cinerea cre1. These results extend knowledge on carbon catabolite repression to litter decomposer Basidiomycota.

Secretome Analysis of Macrophomina phaseolina Identifies an Array of Putative Virulence Factors Responsible for Charcoal Rot Disease in Plants

Macrophomina phaseolina is a global devastating necrotrophic fungal pathogen. It causes charcoal rot disease in more than 500 host plants including major food crops, pulse crops, fiber crops, and oil crops. Despite having the whole-genome sequence of M. phaseolina, understanding the M. phaseolina genome-based plant–pathogen interactions is limited in the absence of direct experimental proof of secretion. Thus, it is essential to understand the host–microbe interaction and the disease pathogenesis, which can ensure global agricultural crop production and security. An in silico–predicted secretome of M. phaseolina is unable to represent the actual secretome. We could identify 117 proteins present in the secretome of M. phaseolina using liquid chromatography–electrospray ionization–tandem mass spectrometry. Data are available via ProteomeXchange with identifier PXD032749. An array of putative virulence factors of M. phaseolina were identified in the present study using solid-state culture. Similar virulence factors have been reported in other plant pathogenic fungi also. Among the secretory fungal proteins with positive economic impacts, lignocellulolytic enzymes are of prime importance. Further, we validated our results by detecting the cell wall–degrading enzymes xylanase, endoglucanase, and amylase in the secretome of M. phaseolina. The present study may provide a better understanding about the necrotrophic fungi M. phaseolina, which modulate the host plant defense barriers using secretory proteins.

Transcriptional Profiles of a Foliar Fungal Endophyte (Pestalotiopsis, Ascomycota) and Its Bacterial Symbiont (Luteibacter, Gammaproteobacteria) Reveal Sulfur Exchange and Growth Regulation during Early Phases of Symbiotic Interaction

Symbiosis with bacteria is widespread among eukaryotes, including fungi. Bacteria that live within fungal mycelia (endohyphal bacteria) occur in many plant-associated fungi, including diverse Mucoromycota and Dikarya. Pestalotiopsis sp. strain 9143 is a filamentous ascomycete isolated originally as a foliar endophyte of Platycladus orientalis (Cupressaceae). It is infected naturally with the endohyphal bacterium Luteibacter sp. strain 9143, which influences auxin and enzyme production by its fungal host. Previous studies have used transcriptomics to examine similar symbioses between endohyphal bacteria and root-associated fungi such as arbuscular mycorrhizal fungi and plant pathogens. However, currently there are no gene expression studies of endohyphal bacteria of Ascomycota, the most species-rich fungal phylum. To begin to understand such symbioses, we developed methods for assessing gene expression by Pestalotiopsis sp. and Luteibacter sp. when grown in coculture and when each was grown axenically. Our assays showed that the density of Luteibacter sp. in coculture was greater than in axenic culture, but the opposite was true for Pestalotiopsis sp. Dual-transcriptome sequencing (RNA-seq) data demonstrate that growing in coculture modulates developmental and metabolic processes in both the fungus and bacterium, potentially through changes in the balance of organic sulfur via methionine acquisition. Our analyses also suggest an unexpected, potential role of the bacterial type VI secretion system in symbiosis establishment, expanding current understanding of the scope and dynamics of fungal-bacterial symbioses.

IMPORTANCE Interactions between microbes and their hosts have important outcomes for host and environmental health. Foliar fungal endophytes that infect healthy plants can harbor facultative endosymbionts called endohyphal bacteria, which can influence the outcome of plant-fungus interactions. These bacterial-fungal interactions can be influential but are poorly understood, particularly from a transcriptome perspective. Here, we report on a comparative, dual-RNA-seq study examining the gene expression patterns of a foliar fungal endophyte and a facultative endohyphal bacterium when cultured together versus separately. Our findings support a role for the fungus in providing organic sulfur to the bacterium, potentially through methionine acquisition, and the potential involvement of a bacterial type VI secretion system in symbiosis establishment. This work adds to the growing body of literature characterizing endohyphal bacterial-fungal interactions, with a focus on a model facultative bacterial-fungal symbiosis in two species-rich lineages, the Ascomycota and Proteobacteria.

Quantitative proteomic analysis of Neosartorya pseudofischeri ascospores subjected to heat treatment

Neosartorya pseudofischeri, a heat- resistant fungus, was isolated from contaminated blueberry juice in our laboratory. To obtain a better understanding of the cellular process of heat- resistant fungus, a quantitative proteomic investigation employing stable isotope labeling by peptide demethylation was conducted on changes in intracellular proteins of N. pseudofischeri ascospores subjected to heat treatment at 93 °C for 0, 1 or 8 min. In total, 811 proteins were identified and quantified. Using the normalized ratio of protein abundance, proteins that changed more than two- fold after heat treatment were identified as significantly increased or decreased proteins and grouped into four clusters based on their quantitative changes. Decreased proteins were found mainly involved in the central carbon metabolism, heat stress responses, reactive oxygen intermediates elimination and translation events. A group of proteins in relation to toxicant degradation and antibiotic neutralization linking to environmental adaptability and tolerance of heat- resistant, was also identified. These findings provide insights into protein changes of N. pseudofischeri ascospores and lay foundations for further investigations on heat- resistant molds using targeted quantitative approaches. SIGNIFICANCE OF THE RESEARCH: Heat- resistant fungi can survive pasteurization processes, and subsequently germinate and grow to cause spoilage of food products, leading to significant economic losses for the food industry and potential health risk for the consumers. To obtain a better understanding of the cellular process of heat- resistant fungi, a quantitative proteomic approach employing stable isotope labeling by peptide demethylation has been used to investigate the intracellular protein changes of N. pseudofischeri ascospores isolated from blueberry juice and subjected to heat treatment. 150 significantly changed proteins were grouped into four clusters based on their quantitative changes. The significant decrease in protein abundance in response to heat treatment revealed possible mechanism that N. pseudofischeri ascospores could survive the heat treatment. This is the first proteomic profile report for N. pseudofischeri. These findings provide insights into protein changes of N. pseudofischeri and lay foundations for further investigations on heat- resistant fungi using targeted quantitative approaches to evaluate the efficiency of thermal treatment for processed food products.

Glycosyl hydrolases family 5, subfamily 5: Relevance and structural insights for designing improved biomass degrading cocktails

Endoglucanases are carbohydrate-degrading enzymes widely used for bioethanol production as part of the enzymatic cocktail. However, family 5 subfamily 5 (GH5_5) endoglucanases are still poorly explored in depth. The Trichoderma reesei representative is the most studied enzyme, presenting catalytic activity in acidic media and mild temperature conditions. Though biochemically similar, its modular structure and synergy with other components vary greatly compared to other GH5_5 members and there is still a lack of specific studies regarding their interaction with other cellulases and application on novel and better mixtures. In this regard, the threedimensional structure elucidation is a highly valuable tool to both uncover basic catalytic mechanisms and implement engineering techniques, proved by the high success rate GH5_5 endoglucanases show. GH5_5 enzymes must be carefully evaluated to fully uncover their potential in biomass-degrading cocktails: the optimal industrial conditions, synergy with other cellulases, structural studies, and enzyme engineering approaches. We aimed to provide the current understanding of these main topics, collecting all available information about characterized GH5_5 endoglucanases function, structure, and bench experiments, in order to suggest future directions to a better application of these enzymes in the industry.

A survey of substrate specificity among Auxiliary Activity Family 5 copper radical oxidases

There is significant contemporary interest in the application of enzymes to replace or augment chemical reagents toward the development of more environmentally sound and sustainable processes. In particular, copper radical oxidases (CRO) from Auxiliary Activity Family 5 Subfamily 2 (AA5_2) are attractive, organic cofactor-free catalysts for the chemoselective oxidation of alcohols to the corresponding aldehydes. These enzymes were first defined by the archetypal galactose-6-oxidase (GalOx, EC 1.1.3.13) from the fungus Fusarium graminearum. The recent discovery of specific alcohol oxidases (EC 1.1.3.7) and aryl alcohol oxidases (EC 1.1.3.47) within AA5_2 has indicated a potentially broad substrate scope among fungal CROs. However, only relatively few AA5_2 members have been characterized to date. Guided by sequence similarity network and phylogenetic analysis, twelve AA5_2 homologs have been recombinantly produced and biochemically characterized in the present study. As defined by their predominant activities, these comprise four galactose 6-oxidases, two raffinose oxidases, four broad-specificity primary alcohol oxidases, and two non-carbohydrate alcohol oxidases. Of particular relevance to applications in biomass valorization, detailed product analysis revealed that two CROs produce the bioplastics monomer furan-2,5-dicarboxylic acid (FDCA) directly from 5-hydroxymethylfurfural (HMF). Furthermore, several CROs could desymmetrize glycerol (a by-product of the biodiesel industry) to D- or L-glyceraldehyde. This study furthers our understanding of CROs by doubling the number of characterized AA5_2 members, which may find future applications as biocatalysts in diverse processes.

Genomic Studies of White-Rot Fungus Cerrena unicolor SP02 Provide Insights into Food Safety Value-Added Utilization of Non-Food Lignocellulosic Biomass

Cerrena unicolor is an ecologically and biotechnologically important wood-degrading basidiomycete with high lignocellulose degrading ability. Biological and genetic investigations are limited in the Cerrena genus and, thus, hinder genetic modification and commercial use. The aim of the present study was to provide a global understanding through genomic and experimental research about lignocellulosic biomass utilization by Cerrena unicolor. In this study, we reported the genome sequence of C. unicolor SP02 by using the Illumina and PacBio 20 platforms to obtain trustworthy assembly and annotation. This is the combinational 2nd and 3rd genome sequencing and assembly of C. unicolor species. The generated genome was 42.79 Mb in size with an N50 contig size of 2.48 Mb, a G + C content of 47.43%, and encoding of 12,277 predicted genes. The genes encoding various lignocellulolytic enzymes including laccase, lignin peroxidase, manganese peroxidase, cytochromes P450, cellulase, xylanase, α-amylase, and pectinase involved in the degradation of lignin, cellulose, xylan, starch, pectin, and chitin that showed the C. unicolor SP02 potentially have a wide range of applications in lignocellulosic biomass conversion. Genome-scale metabolic analysis opened up a valuable resource for a better understanding of carbohydrate-active enzymes (CAZymes) and oxidoreductases that provide insights into the genetic basis and molecular mechanisms for lignocellulosic degradation. The C. unicolor SP02 model can be used for the development of efficient microbial cell factories in lignocellulosic industries. The understanding of the genetic material of C. unicolor SP02 coding for the lignocellulolytic enzymes will significantly benefit us in genetic manipulation, site-directed mutagenesis, and industrial biotechnology.

The Role of Glycoside Hydrolases in Phytopathogenic Fungi and Oomycetes Virulence

Phytopathogenic fungi need to secrete different hydrolytic enzymes to break down complex polysaccharides in the plant cell wall in order to enter the host and develop the disease. Fungi produce various types of cell wall degrading enzymes (CWDEs) during infection. Most of the characterized CWDEs belong to glycoside hydrolases (GHs). These enzymes hydrolyze glycosidic bonds and have been identified in many fungal species sequenced to date. Many studies have shown that CWDEs belong to several GH families and play significant roles in the invasion and pathogenicity of fungi and oomycetes during infection on the plant host, but their mode of function in virulence is not yet fully understood. Moreover, some of the CWDEs that belong to different GH families act as pathogen-associated molecular patterns (PAMPs), which trigger plant immune responses. In this review, we summarize the most important GHs that have been described in eukaryotic phytopathogens and are involved in the establishment of a successful infection.

Defeated by the nines: nine extracellular strategies to avoid microbe-associated molecular patterns recognition in plants

Recognition of microbe-associated molecular patterns (MAMPs) by cell-surface receptors is pivotal in host-microbe interactions. Both pathogens and symbionts establish plant-microbe interactions using fascinating intricate extracellular strategies to avoid recognition. Here we distinguish nine different extracellular strategies to avoid recognition by the host, acting at three different levels. To avoid the accumulation of MAMP precursors (Level 1), microbes take advantage of polymorphisms in both MAMP proteins and glycans, or downregulate MAMP production. To reduce hydrolytic MAMP release (Level 2), microbes shield MAMP precursors with proteins or glycans and inhibit or degrade host-derived hydrolases. And to prevent MAMP perception directly (Level 3), microbes degrade or sequester MAMPs before they are perceived. We discuss examples of these nine strategies and envisage three additional extracellular strategies to avoid MAMP perception in plants.

Harnessing microbial wealth for lignocellulose biomass valorization through secretomics: a review

The recalcitrance of lignocellulosic biomass is a major constraint to its high-value use at industrial scale. In nature, microbes play a crucial role in biomass degradation, nutrient recycling and ecosystem functioning. Therefore, the use of microbes is an attractive way to transform biomass to produce clean energy and high-value compounds. The microbial degradation of lignocelluloses is a complex process which is dependent upon multiple secreted enzymes and their synergistic activities. The availability of the cutting edge proteomics and highly sensitive mass spectrometry tools make possible for researchers to probe the secretome of microbes and microbial consortia grown on different lignocelluloses for the identification of hydrolytic enzymes of industrial interest and their substrate-dependent expression. This review summarizes the role of secretomics in identifying enzymes involved in lignocelluloses deconstruction, the development of enzyme cocktails and the construction of synthetic microbial consortia for biomass valorization, providing our perspectives to address the current challenges.

The Nutritional Significance of Intestinal Fungi: Alteration of Dietary Carbohydrate Composition Triggers Colonic Fungal Community Shifts in a Pig Model

Carbohydrates represent the most important energy source in the diet of humans and animals. A large number of studies have shown that dietary carbohydrates (DCHO) are related to the bacterial community in the gut, but their relationship with the composition of intestinal fungi is still unknown. Here, we report the response of the colonic fungal community to different compositions of DCHO in a pig model. Three factors, ratio (2:1, 1:1, and 1:2) of amylose to amylopectin (AM/AP), level of nonstarch polysaccharides (NSP; 1%, 2%, and 3%), and mannan-oligosaccharide (MOS; 400, 800, and 1,200 mg/kg body weight), were considered according to an L9 (3⁴) orthogonal design to form nine diets with different carbohydrate compositions. Sequencing based on an Illumina HiSeq 2500 platform targeting the internal transcribed spacer 1 region showed that the fungal community in the colon of the pigs responded to DCHO in the order of MOS, AM/AP, and NSP. A large part of some low-abundance fungal genera correlated with the composition of DCHO, represented by Saccharomycopsis, Mrakia, Wallemia, Cantharellus, Eurotium, Solicoccozyma, and Penicillium, were also associated with the concentration of glucose and fructose, as well as the activity of β-d-glucosidase in the colonic digesta, suggesting a role of these fungi in the degradation of DCHO in the colon of pigs. Our study provides direct evidence for the relationship between the composition of DCHO and the fungal community in the colon of pigs, which is helpful to understand the function of gut microorganisms in pigs.IMPORTANCE Although fungi are a large group of microorganisms along with bacteria and archaea in the gut of monogastric animals, the nutritional significance of fungi has been ignored for a long time. Our previous studies revealed a distinct fungal community in the gut of grazing Tibetan pigs (J. Li, D. Chen, B. Yu, J. He, et al., Microb Biotechnol 13:509-521, 2020, https://doi.org/10.1111/1751-7915.13507) and a close correlation between fungal species and short-chain fatty acids, the main microbial metabolites of carbohydrates in the hindgut of pigs (J. Li, Y. Luo, D. Chen, B. Yu, et al., J Anim Physiol Anim Nutr 104:616-628, 2020, https://doi.org/10.1111/jpn.13300). These groundbreaking findings indicate a potential relationship between intestinal fungi and the utilization of DCHO. However, no evidence directly proves the response of intestinal fungi to changes in DCHO. Here, we show a clear alteration of the colonic fungal community in pigs triggered by different compositions of DCHO simulated by varied concentrations of starch, nonstarch polysaccharides (NSP), and oligosaccharides. Our results highlight the potential involvement of intestinal fungi in the utilization of nutrients in monogastric animals.

Polyketide Derivatives in the Resistance of Gerbera hybrida to Powdery Mildew

Powdery mildew is a common disease affecting the commercial production of gerbera flowers (Gerbera hybrida, Asteraceae). Some varieties show a certain degree of resistance to it. Our objective was to identify biomarkers of resistance to powdery mildew using an ¹H nuclear magnetic resonance spectroscopy and chemometrics approach in a complex, fully factorial experiment to suggest a target for selection and breeding. Resistant varieties were found to differ from those that were susceptible in the metabolites of the polyketide pathway, such as gerberin, parasorboside, and gerberinside. A new compound probably involved in resistance, 5-hydroxyhexanoic acid 3-O-β-D-glucoside, was described for the first time. A decision tree model was built to distinguish resistant varieties, with an accuracy of 57.7%, sensitivity of 72%, and specificity of 44.44% in an independent test. Our results suggest the mechanism of resistance to powdery mildew in gerbera and provide a potential tool for resistance screening in breeding programs.

Molecular Dialogues between Early Divergent Fungi and Bacteria in an Antagonism versus a Mutualism

Fungal-bacterial symbioses range from antagonisms to mutualisms and remain one of the least understood interdomain interactions despite their ubiquity as well as ecological and medical importance. To build a predictive conceptual framework for understanding interactions between fungi and bacteria in different types of symbioses, we surveyed fungal and bacterial transcriptional responses in the mutualism between Rhizopus microsporus (Rm) (ATCC 52813, host) and its Mycetohabitans (formerly Burkholderia) endobacteria versus the antagonism between a nonhost Rm (ATCC 11559) and Mycetohabitans isolated from the host, at two time points, before and after partner physical contact. We found that bacteria and fungi sensed each other before contact and altered gene expression patterns accordingly. Mycetohabitans did not discriminate between the host and nonhost and engaged a common set of genes encoding known as well as novel symbiosis factors. In contrast, responses of the host versus nonhost to endobacteria were dramatically different, converging on the altered expression of genes involved in cell wall biosynthesis and reactive oxygen species (ROS) metabolism. On the basis of the observed patterns, we formulated a set of hypotheses describing fungal-bacterial interactions and tested some of them. By conducting ROS measurements, we confirmed that nonhost fungi increased production of ROS in response to endobacteria, whereas host fungi quenched their ROS output, suggesting that ROS metabolism contributes to the nonhost resistance to bacterial infection and the host ability to form a mutualism. Overall, our study offers a testable framework of predictions describing interactions of early divergent Mucoromycotina fungi with bacteria.IMPORTANCE Animals and plants interact with microbes by engaging specific surveillance systems, regulatory networks, and response modules that allow for accommodation of mutualists and defense against antagonists. Antimicrobial defense responses are mediated in both animals and plants by innate immunity systems that owe their functional similarities to convergent evolution. Like animals and plants, fungi interact with bacteria. However, the principles governing these relations are only now being discovered. In a study system of host and nonhost fungi interacting with a bacterium isolated from the host, we found that bacteria used a common gene repertoire to engage both partners. In contrast, fungal responses to bacteria differed dramatically between the host and nonhost. These findings suggest that as in animals and plants, the genetic makeup of the fungus determines whether bacterial partners are perceived as mutualists or antagonists and what specific regulatory networks and response modules are initiated during each encounter.

Genome-wide study of saprotrophy-related genes in the basal fungus Conidiobolus heterosporus

Conidiobolus spp. are important saprophytic basal fungi. However, to date, no genomic-level data for decaying plant materials in the genus Conidiobolus has been reported. Here, we report that the 33.4-Mb genome of Conidiobolus heterosporus encodes 10,857 predicted genes. Conidiobolus heterosporus harbors 394 CAZyme-encoding genes belonging to 4 major modules but does not encode a polysaccharide lyase (PL). Many carbohydrate esterases (CEs) belonging to the family CE12 play crucial roles as pectin acetylesterases, and 14 genes were upregulated in the IM (fungus grown on inducing medium) among 17 expressed CE12 family genes. In addition, most of the genes in the GH132 CAZyme family showed a greater than 5-fold increase in expression in the IM compared with that in the wild type. Furthermore, 122 P450-encoding genes grouped into 11 families were detected in the fungal genome, most of which belonged to the CYP547 family (36 genes) followed by CYP548 (27 genes) and CYP5856 (25 genes). Interestingly, members of the families CYP5014 and CYP5136 were identified, the first time such enzymes have been described in a fungus. Our findings provide new insights into the genomics and genomic features of the saprophytic basal fungus C. heterosporus.Key Points• Genome of the saprobiotic basal fungus C. heterosporus was sequenced and analyzed.• 394 CAZymes but no PL family genes were found and expression levels were determined.• CE12 and GH132 proteins may play roles in the pectin and plant material degradation.• A large number of P450s but few P450 families existed in the fungus.

In silico Proteomic Analysis Provides Insights Into Phylogenomics and Plant Biomass Deconstruction Potentials of the Tremelalles

Basidiomycetes populate a wide range of ecological niches but unlike ascomycetes, their capabilities to decay plant polymers and their potential for biotechnological approaches receive less attention. Particularly, identification and isolation of CAZymes is of biotechnological relevance and has the potential to improve the cache of currently available commercial enzyme cocktails toward enhanced plant biomass utilization. The order Tremellales comprises phylogenetically diverse fungi living as human pathogens, mycoparasites, saprophytes or associated with insects. Here, we have employed comparative genomics approaches to highlight the phylogenomic relationships among thirty-five Tremellales and to identify putative enzymes of biotechnological interest encoded on their genomes. Evaluation of the predicted proteomes of the thirty-five Tremellales revealed 6,918 putative carbohydrate-active enzymes (CAZYmes) and 7,066 peptidases. Two soil isolates, Saitozyma podzolica DSM 27192 and Cryptococcus sp. JCM 24511, show higher numbers harboring an average of 317 compared to a range of 267-121 CAZYmes for the rest of the strains. Similarly, the proteomes of the two soil isolates along with two plant associated strains contain higher number of peptidases sharing an average of 234 peptidases compared to a range of 226-167 for the rest of the strains. Despite these huge differences and the apparent enrichment of these enzymes among the soil isolates, the data revealed a diversity of the various enzyme families that does not reflect specific habitat type. Growth experiment on various carbohydrates to validate the predictions provides support for this view. Overall, the data indicates that the Tremellales could serve as a rich source of both CAZYmes and peptidases with wide range of potential biotechnological relevance.

Comparative Genomics Uncovers the Genetic Diversity and Synthetic Biology of Secondary Metabolite Production of Trametes

The carbohydrate-active enzyme (CAZyme) genes of Trametes contribute to polysaccharide degradation. However, the comprehensive analysis of the composition of CAZymes and the biosynthetic gene clusters (BGCs) of Trametes remain unclear. Here, we conducted comparative analysis, detected the CAZyme genes, and predicted the BGCs for nine Trametes strains. Among the 82,053 homologous clusters obtained for Trametes, we identified 8518 core genes, 60,441 accessory genes, and 13,094 specific genes. A large proportion of CAZyme genes were cataloged into glycoside hydrolases, glycosyltransferases, and carbohydrate esterases. The predicted BGCs of Trametes were divided into six strategies, and the nine Trametes strains harbored 47.78 BGCs on average. Our study revealed that Trametes exhibits an open pan-genome structure. These findings provide insights into the genetic diversity and explored the synthetic biology of secondary metabolite production for Trametes.

Network analysis exposes core functions in major lifestyles of fungal and oomycete plant pathogens

Background

Genomic studies demonstrate that components of virulence mechanisms in filamentous eukaryotic pathogens (FEPs, fungi and oomycetes) of plants are often highly conserved, or found in gene families that include secreted hydrolytic enzymes (e.g., cellulases and proteases) and secondary metabolites (e.g., toxins), central to the pathogenicity process. However, very few large-scale genomic comparisons have utilized complete proteomes from dozens of FEPs to reveal lifestyle-associated virulence mechanisms. Providing a powerful means for exploration, and the discovery of trends in large-scale datasets, network analysis has been used to identify core functions of the primordial cyanobacteria, and ancient evolutionary signatures in oxidoreductases.

Results

We used a sequence-similarity network to study components of virulence mechanisms of major pathogenic lifestyles (necrotroph (ic), N; biotroph (ic), B; hemibiotroph (ic), H) in complete pan-proteomes of 65 FEPs and 17 saprobes. Our comparative analysis highlights approximately 190 core functions found in 70% of the genomes of these pathogenic lifestyles. Core functions were found mainly in: transport (in H, N, B cores); carbohydrate metabolism, secondary metabolite synthesis, and protease (H and N cores); nucleic acid metabolism and signal transduction (B core); and amino acid metabolism (H core). Taken together, the necrotrophic core contains functions such as cell wall-associated degrading enzymes, toxin metabolism, and transport, which are likely to support their lifestyle of killing prior to feeding. The biotrophic stealth growth on living tissues is potentially controlled by a core of regulatory functions, such as: small G-protein family of GTPases, RNA modification, and cryptochrome-based light sensing. Regulatory mechanisms found in the hemibiotrophic core contain light- and CO₂-sensing functions that could mediate important roles of this group, such as transition between lifestyles.

Conclusions

The selected set of enriched core functions identified in our work can facilitate future studies aimed at controlling FEPs. One interesting example would be to facilitate the identification of the pathogenic potential of samples analyzed by metagenomics. Finally, our analysis offers potential evolutionary scenarios, suggesting that an early-branching saprobe (identified in previous studies) has probably evolved a necrotrophic lifestyle as illustrated by the highest number of shared gene families between saprobes and necrotrophs.

Genomes and secretomes of Ascomycota fungi reveal diverse functions in plant biomass decomposition and pathogenesis

BACKGROUND

The dominant fungi in arid grasslands and shrublands are members of the Ascomycota phylum. Ascomycota fungi are important drivers in carbon and nitrogen cycling in arid ecosystems. These fungi play roles in soil stability, plant biomass decomposition, and endophytic interactions with plants. They may also form symbiotic associations with biocrust components or be latent saprotrophs or pathogens that live on plant tissues. However, their functional potential in arid soils, where organic matter, nutrients and water are very low or only periodically available, is poorly characterized.

RESULTS

Five Ascomycota fungi were isolated from different soil crust microhabitats and rhizosphere soils around the native bunchgrass Pleuraphis jamesii in an arid grassland near Moab, UT, USA. Putative genera were Coniochaeta, isolated from lichen biocrust, Embellisia from cyanobacteria biocrust, Chaetomium from below lichen biocrust, Phoma from a moss microhabitat, and Aspergillus from the soil. The fungi were grown in replicate cultures on different carbon sources (chitin, native bunchgrass or pine wood) relevant to plant biomass and soil carbon sources. Secretomes produced by the fungi on each substrate were characterized. Results demonstrate that these fungi likely interact with primary producers (biocrust or plants) by secreting a wide range of proteins that facilitate symbiotic associations. Each of the fungal isolates secreted enzymes that degrade plant biomass, small secreted effector proteins, and proteins involved in either beneficial plant interactions or virulence. Aspergillus and Phoma expressed more plant biomass degrading enzymes when grown in grass- and pine-containing cultures than in chitin. Coniochaeta and Embellisia expressed similar numbers of these enzymes under all conditions, while Chaetomium secreted more of these enzymes in grass-containing cultures.

CONCLUSIONS

This study of Ascomycota genomes and secretomes provides important insights about the lifestyles and the roles that Ascomycota fungi likely play in arid grassland, ecosystems. However, the exact nature of those interactions, whether any or all of the isolates are true endophytes, latent saprotrophs or opportunistic phytopathogens, will be the topic of future studies.

iTRAQ-based proteomic analysis reveals the mechanisms of Botrytis cinerea controlled with Wuyiencin

Background

Grey mould is an important plant disease worldwide, caused by Botrytis cinerea, resulting in serious economic loss. Wuyiencin, a low toxicity, high efficiency, and broad-spectrum agricultural antibiotic, has been demonstrated effectiveness against B. cinerea.

Results

Wuyiencin treatment inhibited growth and sporulation of B. cinerea, specifically altering hypha morphology and intracellular structures. These changes were accompanied by differential expression (fold change > 2.0) of 316 proteins identified by iTRAQ-labelling LC-MS/MS analysis (P < 0.05). Up-regulation of 14 proteins, including carbohydrate metabolism proteins and cell wall stabilization proteins, was validated by parallel reaction monitoring (PRM). Down-regulation of 13 proteins was validated by PRM, including regulators of energy metabolism, nucleotide/protein synthesis, and the biosynthesis of mediators of plant stress and decay.

Conclusion

Our results confirm the inhibitory biological effects of wuyiencin on B. cinereal and elaborate on the differentially expressed proteins and associated pathways implicated in the capacity of wuyiencin to debilitate the growth and pathogenicity of grey mould. This study provides validated candidates for further targeted exploration with the goal of optimizing wuyiencin as a safe, low-toxicity agent for biological control.

Redesigning N-glycosylation sites in a GH3 β-xylosidase improves the enzymatic efficiency

BACKGROUND

β-Xylosidases are glycoside hydrolases (GHs) that cleave xylooligosaccharides and/or xylobiose into shorter oligosaccharides and xylose. Aspergillus nidulans is an established genetic model and good source of carbohydrate-active enzymes (CAZymes). Most fungal enzymes are N-glycosylated, which influences their secretion, stability, activity, signalization, and protease protection. A greater understanding of the N-glycosylation process would contribute to better address the current bottlenecks in obtaining high secretion yields of fungal proteins for industrial applications.

RESULTS

In this study, BxlB-a highly secreted GH3 β-xylosidase from A. nidulans, presenting high activity and several N-glycosylation sites-was selected for N-glycosylation engineering. Several glycomutants were designed to investigate the influence of N-glycans on BxlB secretion and function. The non-glycosylated mutant (BxlB^non-glyc) showed similar levels of enzyme secretion and activity compared to the wild-type (BxlB^wt), while a partially glycosylated mutant (BxlB^N1;5;7) exhibited increased activity. Additionally, there was no enzyme secretion in the mutant in which the N-glycosylation context was changed by the introduction of four new N-glycosylation sites (BxlB^CC), despite the high transcript levels. BxlB^wt, BxlB^non-glyc, and BxlB^N1;5;7 formed similar secondary structures, though the mutants had lower melting temperatures compared to the wild type. Six additional glycomutants were designed based on BxlB^N1;5;7, to better understand its increased activity. Among them, the two glycomutants which maintained only two N-glycosylation sites each (BxlB^N1;5 and BxlB^N5;7) showed improved catalytic efficiency, whereas the other four mutants' catalytic efficiencies were reduced. The N-glycosylation site N5 is important for improved BxlB catalytic efficiency, but needs to be complemented by N1 and/or N7. Molecular dynamics simulations of BxlB^non-glyc and BxlB^N1;5 reveals that the mobility pattern of structural elements in the vicinity of the catalytic pocket changes upon N1 and N5 N-glycosylation sites, enhancing substrate binding properties which may underlie the observed differences in catalytic efficiency between BxlB^non-glyc and BxlB^N1;5.

CONCLUSIONS

This study demonstrates the influence of N-glycosylation on A. nidulans BxlB production and function, reinforcing that protein glycoengineering is a promising tool for enhancing thermal stability, secretion, and enzymatic activity. Our report may also support biotechnological applications for N-glycosylation modification of other CAZymes.

Lignin degradation potential and draft genome sequence of Trametes trogii S0301

BACKGROUND

Trametes trogii is a member of the white-rot fungi family, which has a unique ability to break down recalcitrant lignin polymers to CO2 and water, and they have enormous potential to biodegrade a wide range of toxic environmental pollutants. Because of its industrial potential, the identification of lignin-degrading enzyme systems in Trametes is an important area of research. Development and utilization of industrial value genes are suffering due to deficiency knowledge of genome available for their manipulation.

RESULTS

In the present study, Homokaryotic strains of T. trogii S0301 were screened and sequencing by PacBio Sequel II platform. The final draft genome is ~ 39.88 Mb, with a contig N50 size of 2.4 Mb, this was the first genome sequencing and assembly of T. trogii species. Further analyses predicted 14,508 protein-coding genes. Results showed that T. trogii S0301 contains 602 genes encoding CAZymes, include 211 glycoside hydrolase and 117 lignin-degrading family genes, nine laccases related genes. Small subunit ribosomal RNA gene (18S rRNA) sequencing confirms its phylogenetic position. Moreover, T. trogii S0301 has the largest number of cytochromes P450 (CYPs) superfamily genes compare to other fungi. All these results are consistent with enzymatic assays and transcriptome analysis results. We also analyzed other genome characteristics in the T. trogii S0301genome.

CONCLUSION

Here, we present a nearly complete genome for T. trogii S0301, which will help elucidate the biosynthetic pathways of the lignin-degrading enzyme, advancing the discovery, characterization, and modification of novel enzymes from this genus. This genome sequence will provide a valuable reference for the investigation of lignin degradation in the Trametes genus.

The Genome of the Fungal Pathogen Verticillium dahliae Reveals Extensive Bacterial to Fungal Gene Transfer

Horizontal gene transfer (HGT) involves the transmission of genetic material between distinct evolutionary lineages and can be an important source of biological innovation. Reports of interkingdom HGT to eukaryotic microbial pathogens have accumulated over recent years. Verticillium dahliae is a notorious plant pathogen that causes vascular wilt disease on hundreds of plant species, resulting in high economic losses every year. Previously, the effector gene Ave1 and a glucosyltransferase-encoding gene were identified as virulence factor-encoding genes that were proposed to be horizontally acquired from a plant and a bacterial donor, respectively. However, to what extent HGT contributed to the overall genome composition of V. dahliae remained elusive. Here, we systematically searched for evidence of interkingdom HGT events in the genome of V. dahliae and provide evidence for extensive horizontal gene acquisition from bacterial origin.

Global transcriptome analysis of the aphelid Paraphelidium tribonemae supports the phagotrophic origin of fungi

Aphelids are little-known phagotrophic parasites of algae whose life cycle and morphology resemble those of the parasitic rozellids (Cryptomycota, Rozellomycota). In previous phylogenetic analyses of RNA polymerase and rRNA genes, aphelids, rozellids and Microsporidia (parasites of animals) formed a clade, named Opisthosporidia, which appeared as the sister group to Fungi. However, the statistical support for the Opisthosporidia was always moderate. Here, we generated full life-cycle transcriptome data for the aphelid species Paraphelidium tribonemae. In-depth multi-gene phylogenomic analyses using several protein datasets place this aphelid as the closest relative of fungi to the exclusion of rozellids and Microsporidia. In contrast with the comparatively reduced Rozella allomycis genome, we infer a rich, free-living-like aphelid proteome, with a metabolism similar to fungi, including cellulases likely involved in algal cell-wall penetration and enzymes involved in chitin biosynthesis. Our results suggest that fungi evolved from complex aphelid-like ancestors that lost phagotrophy and became osmotrophic.

Ultra-low input transcriptomics reveal the spore functional content and phylogenetic affiliations of poorly studied arbuscular mycorrhizal fungi

Arbuscular mycorrhizal fungi (AMF) are a group of soil microorganisms that establish symbioses with the vast majority of land plants. To date, generation of AMF coding information has been limited to model genera that grow well axenically; Rhizoglomus and Gigaspora. Meanwhile, data on the functional gene repertoire of most AMF families is non-existent. Here, we provide primary large-scale transcriptome data from eight poorly studied AMF species (Acaulospora morrowiae, Diversispora versiforme, Scutellospora calospora, Racocetra castanea, Paraglomus brasilianum, Ambispora leptoticha, Claroideoglomus claroideum and Funneliformis mosseae) using ultra-low input ribonucleic acid (RNA)-seq approaches. Our analyses reveals that quiescent spores of many AMF species harbour a diverse functional diversity and solidify known evolutionary relationships within the group. Our findings demonstrate that RNA-seq data obtained from low-input RNA are reliable in comparison to conventional RNA-seq experiments. Thus, our methodology can potentially be used to deepen our understanding of fungal microbial function and phylogeny using minute amounts of RNA material.

Genomic and exoproteomic diversity in plant biomass degradation approaches among Aspergilli

We classified the genes encoding carbohydrate-active enzymes (CAZymes) in 17 sequenced genomes representing 16 evolutionarily diverse Aspergillus species. We performed a phylogenetic analysis of the encoding enzymes, along with experimentally characterized CAZymes, to assign molecular function to the Aspergilli CAZyme families and subfamilies. Genome content analysis revealed that the numbers of CAZy genes per CAZy family related to plant biomass degradation follow closely the taxonomic distance between the species. On the other hand, growth analysis showed almost no correlation between the number of CAZyme genes and the efficiency in polysaccharide utilization. The exception is A. clavatus where a reduced number of pectinolytic enzymes can be correlated with poor growth on pectin. To gain detailed information on the enzymes used by Aspergilli to breakdown complex biomass, we conducted exoproteome analysis by mass spectrometry. These results showed that Aspergilli produce many different enzymes mixtures in the presence of sugar beet pulp and wheat bran. Despite the diverse enzyme mixtures produced, species of section Nigri, A. aculeatus, A. nidulans and A. terreus, produce mixtures of enzymes with activities that are capable of digesting all the major polysaccharides in the available substrates, suggesting that they are capable of degrading all the polysaccharides present simultaneously. For the other Aspergilli, typically the enzymes produced are targeted to a subset of polysaccharides present, suggesting that they can digest only a subset of polysaccharides at a given time.

Scaling up data curation using deep learning: An application to literature triage in genomic variation resources

Manually curating biomedical knowledge from publications is necessary to build a knowledge based service that provides highly precise and organized information to users. The process of retrieving relevant publications for curation, which is also known as document triage, is usually carried out by querying and reading articles in PubMed. However, this query-based method often obtains unsatisfactory precision and recall on the retrieved results, and it is difficult to manually generate optimal queries. To address this, we propose a machine-learning assisted triage method. We collect previously curated publications from two databases UniProtKB/Swiss-Prot and the NHGRI-EBI GWAS Catalog, and used them as a gold-standard dataset for training deep learning models based on convolutional neural networks. We then use the trained models to classify and rank new publications for curation. For evaluation, we apply our method to the real-world manual curation process of UniProtKB/Swiss-Prot and the GWAS Catalog. We demonstrate that our machine-assisted triage method outperforms the current query-based triage methods, improves efficiency, and enriches curated content. Our method achieves a precision 1.81 and 2.99 times higher than that obtained by the current query-based triage methods of UniProtKB/Swiss-Prot and the GWAS Catalog, respectively, without compromising recall. In fact, our method retrieves many additional relevant publications that the query-based method of UniProtKB/Swiss-Prot could not find. As these results show, our machine learning-based method can make the triage process more efficient and is being implemented in production so that human curators can focus on more challenging tasks to improve the quality of knowledge bases.

As the volume of literature on genomic variants continues to grow at an increasing rate, it is becoming more difficult for a curator of a variant knowledge base to keep up with and curate all the published papers. Here, we suggest a deep learning-based literature triage method for genomic variation resources. Our method achieves state-of-the-art performance on the triage task. Moreover, our model does not require any laborious preprocessing or feature engineering steps, which are required for traditional machine learning triage methods. We applied our method to the literature triage process of UniProtKB/Swiss-Prot and the NHGRI-EBI GWAS Catalog for genomic variation by collaborating with the database curators. Both the manual curation teams confirmed that our method achieved higher precision than their previous query-based triage methods without compromising recall. Both results show that our method is more efficient and can replace the traditional query-based triage methods of manually curated databases. Our method can give human curators more time to focus on more challenging tasks such as actual curation as well as the discovery of novel papers/experimental techniques to consider for inclusion.

Application of Transcriptomics to Compare the Carbohydrate Active Enzymes That Are Expressed by Diverse Genera of Anaerobic Fungi to Degrade Plant Cell Wall Carbohydrates

The efficiency with which the anaerobic fungi (phylum Neocallimastigomycota) degrade plant biomass is well-recognized and in recent years has received renewed interest. To further understand the biological mechanisms that are utilized by the rumen anaerobic fungi to break down lignocellulose, we have used a transcriptomic approach to examine carbohydrate digestion by Neocallimastix frontalis, Piromyces rhizinflata, Orpinomyces joyonii, and Anaeromyces mucronatus cultured on several carbon sources. The number of predicted unique transcripts ranged from 6,633 to 12,751. Pfam domains were identified in 62–70% of the fungal proteins and were linked to gene ontology terms to infer the biological function of the transcripts. Most of the predicted functions are consistent across species suggesting a similar overall strategy evolved for successful colonization of the rumen. However, the presence of differential profiles in enzyme classes suggests that there may be also be niche specialization. All fungal species were found to express an extensive array of transcripts encoding carbohydrate active enzymes (CAZymes) ranging from 8.3 to 11.3% of the transcriptome. CAZyme families involved in hemicellulose digestion were the most abundant across all four fungi. This study provides additional insight into how anaerobic fungi have evolved to become specialists at breaking down the plant cell wall in the complex and, strictly anaerobic rumen ecosystem.

De novo assembly and comparative transcriptome analysis of Monilinia fructicola, Monilinia laxa and Monilinia fructigena, the causal agents of brown rot on stone fruits

BACKGROUND

Brown rots are important fungal diseases of stone and pome fruits. They are caused by several Monilinia species but M. fructicola, M. laxa and M. fructigena are the most common all over the world. Although they have been intensively studied, the availability of genomic and transcriptomic data in public databases is still scant. We sequenced, assembled and annotated the transcriptomes of the three pathogens using mRNA from germinating conidia and actively growing mycelia of two isolates of opposite mating types per each species for comparative transcriptome analyses.

RESULTS

Illumina sequencing was used to generate about 70 million of paired-end reads per species, that were de novo assembled in 33,861 contigs for M. fructicola, 31,103 for M. laxa and 28,890 for M. fructigena. Approximately, 50% of the assembled contigs had significant hits when blasted against the NCBI non-redundant protein database and top-hits results were represented by Botrytis cinerea, Sclerotinia sclerotiorum and Sclerotinia borealis proteins. More than 90% of the obtained sequences were complete, the percentage of duplications was always less than 14% and fragmented and missing transcripts less than 5%. Orthologous transcripts were identified by tBLASTn analysis using the B. cinerea proteome as reference. Comparative transcriptome analyses revealed 65 transcripts over-expressed (FC ≥ 8 and FDR ≤ 0.05) or unique in M. fructicola, 30 in M. laxa and 31 in M. fructigena. Transcripts were involved in processes affecting fungal development, diversity and host-pathogen interactions, such as plant cell wall-degrading and detoxifying enzymes, zinc finger transcription factors, MFS transporters, cell surface proteins, key enzymes in biosynthesis and metabolism of antibiotics and toxins, and transposable elements.

CONCLUSIONS

This is the first large-scale reconstruction and annotation of the complete transcriptomes of M. fructicola, M. laxa and M. fructigena and the first comparative transcriptome analysis among the three pathogens revealing differentially expressed genes with potential important roles in metabolic and physiological processes related to fungal morphogenesis and development, diversity and pathogenesis which need further investigations. We believe that the data obtained represent a cornerstone for research aimed at improving knowledge on the population biology, physiology and plant-pathogen interactions of these important phytopathogenic fungi.

Manual Gene Curation and Functional Annotation

No genome sequencing project is complete without structural and functional annotation. Gene models and functional predictions for these models can be obtained relatively easily using computational methods, but they are prone to errors. We describe herein the steps we use to manually curate gene models and functionally annotate them. Our approach is to examine each gene model carefully, and improve its structure if necessary, using a comprehensive set of experimental and computational data as evidence. Then, functional predictions are assigned to the gene models based on conserved protein domains and sequence similarities. We use stringent sequence similarity cutoffs and reviewed sequence-database records as external sources for our annotations. By methodically choosing which evidence to use for each annotation, we minimize the risk of adopting and assigning false predictions to the gene models.

Phylogenetic Analysis of Protein Family

With the number of sequenced genomes increasing rapidly, it is impractical to perform functional and structural analyses on all individual proteins. Phylogenetic analysis employs a combination of molecular and statistical approaches to infer or estimate relationships among individuals. It provides a credible method to explore the relationship between sequence similarity and function of proteins belonging to the same family. This chapter describes a standardized framework of phylogenetic analysis to study large protein families. Bioinformatic approaches and online tools used in phylogenetic analyses are presented.

In vitro investigations on extracellular proteins secreted by Aphanomyces invadans, the causative agent of epizootic ulcerative syndrome

BACKGROUND

Proteases produced by many microorganisms, including oomycetes, are crucial for their growth and development. They may also play a critical role in disease manifestation. Epizootic ulcerative syndrome is one of the most destructive fish diseases known. It is caused by the oomycete Aphanomyces invadans and leads to mass mortalities of cultured and wild fish in many countries. The areas of concern are Australia, China, Japan, South and Southeast Asian countries and the USA. Extracellular proteases produced by this oomycete are believed to trigger EUS pathogenesis in fish. To address this activity, we collected the extracellular products (ECP) of A. invadans and identified the secreted proteins using SDS-PAGE and mass spectrometery. A. invadans was cultivated in liquid Glucose-Peptone-Yeats media. The culture media was ultra-filtered through 10 kDa filters and analysed using SDS-PAGE. Three prominent protein bands from the SDS gel were excised and identified by mass spectrometery. Furthermore, we assessed their proteolytic effect on casein and immunoglobulin M (IgM) of rainbow trout (Oncorhynchus mykiss) and giant gourami (Osphronemus goramy). Antiprotease activity of the fish serum was also investigated.

RESULTS

BLASTp analysis revealed that the prominent secreted proteins were proteases, mainly of the serine and cysteine types. Proteins containing fascin-like domain and bromodomain were also identified. We could demonstrate that the secreted proteases showed proteolytic activity against the casein and the IgM of both fish species. The anti-protease activity experiment showed that the percent inhibition of the common carp serum was 94.2% while that of rainbow trout and giant gourami serum was 7.7 and 12.9%, respectively.

CONCLUSIONS

The identified proteases, especially serine proteases, could be the potential virulence factors in A. invadans and, hence, are candidates for further functional and host-pathogen interaction studies. The role of identified structural proteins in A. invadans also needs to be investigated further.

Carbohydrate-active enzymes in Trichoderma harzianum: a bioinformatic analysis bioprospecting for key enzymes for the biofuels industry

BACKGROUND

Trichoderma harzianum is used in biotechnology applications due to its ability to produce powerful enzymes for the conversion of lignocellulosic substrates into soluble sugars. Active enzymes involved in carbohydrate metabolism are defined as carbohydrate-active enzymes (CAZymes), and the most abundant family in the CAZy database is the glycoside hydrolases. The enzymes of this family play a fundamental role in the decomposition of plant biomass.

RESULTS

In this study, the CAZymes of T. harzianum were identified and classified using bioinformatic approaches after which the expression profiles of all annotated CAZymes were assessed via RNA-Seq, and a phylogenetic analysis was performed. A total of 430 CAZymes (3.7% of the total proteins for this organism) were annotated in T. harzianum, including 259 glycoside hydrolases (GHs), 101 glycosyl transferases (GTs), 6 polysaccharide lyases (PLs), 22 carbohydrate esterases (CEs), 42 auxiliary activities (AAs) and 46 carbohydrate-binding modules (CBMs). Among the identified T. harzianum CAZymes, 47% were predicted to harbor a signal peptide sequence and were therefore classified as secreted proteins. The GH families were the CAZyme class with the greatest number of expressed genes, including GH18 (23 genes), GH3 (17 genes), GH16 (16 genes), GH2 (13 genes) and GH5 (12 genes). A phylogenetic analysis of the proteins in the AA9/GH61, CE5 and GH55 families showed high functional variation among the proteins.

CONCLUSIONS

Identifying the main proteins used by T. harzianum for biomass degradation can ensure new advances in the biofuel production field. Herein, we annotated and characterized the expression levels of all of the CAZymes from T. harzianum, which may contribute to future studies focusing on the functional and structural characterization of the identified proteins.

Comparative transcriptome analysis and identification of candidate effectors in two related rust species (Gymnosporangium yamadae and Gymnosporangium asiaticum)

BACKGROUND

Rust fungi constitute the largest group of plant fungal pathogens. However, a paucity of data, including genomic sequences, transcriptome sequences, and associated molecular markers, hinders the development of inhibitory compounds and prevents their analysis from an evolutionary perspective. Gymnosporangium yamadae and G. asiaticum are two closely related rust fungal species, which are ecologically and economically important pathogens that cause apple rust and pear rust, respectively, proved to be devastating to orchards. In this study, we investigated the transcriptomes of these two Gymnosporangium species during the telial stage of their lifecycles. The aim of this study was to understand the evolutionary patterns of these two related fungi and to identify genes that developed by selection.

RESULTS

The transcriptomes of G. yamadae and G. asiaticum were generated from a mixture of RNA from three biological replicates of each species. We obtained 49,318 and 54,742 transcripts, with N50 values of 1957 and 1664, for G. yamadae and G. asiaticum, respectively. We also identified a repertoire of candidate effectors and other gene families associated with pathogenicity. A total of 4947 pairs of putative orthologues between the two species were identified. Estimation of the non-synonymous/synonymous substitution rate ratios for these orthologues identified 116 pairs with Ka/Ks values greater than1 that are under positive selection and 170 pairs with Ka/Ks values of 1 that are under neutral selection, whereas the remaining 4661 genes are subjected to purifying selection. We estimate that the divergence time between the two species is approximately 5.2 Mya.

CONCLUSION

This study constitutes a de novo assembly and comparative analysis between the transcriptomes of the two rust species G. yamadae and G. asiaticum. The results identified several orthologous genes, and many expressed genes were identified by annotation. Our analysis of Ka/Ks ratios identified orthologous genes subjected to positive or purifying selection. An evolutionary analysis of these two species provided a relatively precise divergence time. Overall, the information obtained in this study increases the genetic resources available for research on the genetic diversity of the Gymnosporangium genus.

Protein Bioinformatics Databases and Resources

Many publicly available data repositories and resources have been developed to support protein-related information management, data-driven hypothesis generation, and biological knowledge discovery. To help researchers quickly find the appropriate protein-related informatics resources, we present a comprehensive review (with categorization and description) of major protein bioinformatics databases in this chapter. We also discuss the challenges and opportunities for developing next-generation protein bioinformatics databases and resources to support data integration and data analytics in the Big Data era.

Phylogenetic analysis of β-xylanase SRXL1 of Sporisorium reilianum and its relationship with families (GH10 and GH11) of Ascomycetes and Basidiomycetes

In this paper, the amino acid sequence of the β-xylanase SRXL1 of Sporisorium reilianum, which is a pathogenic fungus of maize was used as a model protein to find its phylogenetic relationship with other xylanases of Ascomycetes and Basidiomycetes and the information obtained allowed to establish a hypothesis of monophyly and of biological role. 84 amino acid sequences of β-xylanase obtained from the GenBank database was used. Groupings analysis of higher-level in the Pfam database allowed to determine that the proteins under study were classified into the GH10 and GH11 families, based on the regions of highly conserved amino acids, 233-318 and 180-193 respectively, where glutamate residues are responsible for the catalysis.

Transcriptome analysis of the white pine blister rust pathogen Cronartium ribicola: de novo assembly, expression profiling, and identification of candidate effectors

Background

The fungus Cronartium ribicola (Cri) is an economically and ecologically important forest pathogen that causes white pine blister rust (WPBR) disease on five-needle pines. To cause stem cankers and kill white pine trees the fungus elaborates a life cycle with five stages of spore development on five-needle pines and the alternate host Ribes plants. To increase our understanding of molecular WP-BR interactions, here we report genome-wide transcriptional profile analysis of C. ribicola using RNA-seq.

Results

cDNA libraries were constructed from aeciospore, urediniospore, and western white pine (Pinus monticola) tissues post Cri infection. Over 200 million RNA-seq 100-bp paired-end (PE) reads from rust fungal spores were de novo assembled and a reference transcriptome was generated with 17,880 transcripts that were expressed from 13,629 unigenes. A total of 734 unique proteins were predicted as a part of the Cri secretome from complete open reading frames (ORFs), and 41 % of them were Cronartium-specific. This study further identified a repertoire of candidate effectors and other pathogenicity determinants. Differentially expressed genes (DEGs) were identified to gain an understanding of molecular events important during the WPBR fungus life cycle by comparing Cri transcriptomes at different infection stages. Large-scale changes of in planta gene expression profiles were observed, revealing that multiple fungal biosynthetic pathways were enhanced during mycelium growth inside infected pine stem tissues. Conversely, many fungal genes that were up-regulated at the urediniospore stage appeared to be signalling components and transporters. The secreted fungal protein genes that were up-regulated in pine needle tissues during early infection were primarily associated with cell wall modifications, possibly to mask the rust pathogen from plant defenses.

Conclusion

This comprehensive transcriptome profiling substantially improves our current understanding of molecular WP-BR interactions. The repertoire of candidate effectors and other putative pathogenicity determinants identified here are valuable for future functional analysis of Cri virulence and pathogenicity.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1861-1) contains supplementary material, which is available to authorized users.

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.

mycoCLAP, the database for characterized lignocellulose-active proteins of fungal origin: resource and text mining curation support

Enzymes active on components of lignocellulosic biomass are used for industrial applications ranging from food processing to biofuels production. These include a diverse array of glycoside hydrolases, carbohydrate esterases, polysaccharide lyases and oxidoreductases. Fungi are prolific producers of these enzymes, spurring fungal genome sequencing efforts to identify and catalogue the genes that encode them. To facilitate the functional annotation of these genes, biochemical data on over 800 fungal lignocellulose-degrading enzymes have been collected from the literature and organized into the searchable database, mycoCLAP (http://mycoclap.fungalgenomics.ca). First implemented in 2011, and updated as described here, mycoCLAP is capable of ranking search results according to closest biochemically characterized homologues: this improves the quality of the annotation, and significantly decreases the time required to annotate novel sequences. The database is freely available to the scientific community, as are the open source applications based on natural language processing developed to support the manual curation of mycoCLAP. Database URL: http://mycoclap.fungalgenomics.ca.

Genome sequence, comparative analysis, and evolutionary insights into chitinases of entomopathogenic fungus Hirsutella thompsonii

Hirsutella thompsonii (Ht) is a fungal pathogen of acarines and the primary cause of epizootics among mites. The draft genomes of two isolates of Ht (MTCC 3556: Ht3, 34.6 Mb and MTCC 6686: Ht6, 34.7 Mb) are presented and compared with the genomes of Beauveria bassiana (Bb) ARSEF 2860 and Ophiocordyceps sinensis (Os) CO18. Comparative analysis of carbohydrate active enzymes, pathogen–host interaction genes, metabolism-associated genes, and genes involved in biosynthesis of secondary metabolites in the four genomes was carried out. Reduction in gene family sizes in Ht3 and Os as compared with Ht6 and Bb is observed. Analysis of the mating type genes in Ht reveals the presence of MAT idiomorphs which is suggestive of cryptic sexual traits in Ht. We further identify and classify putative chitinases that may function as virulence factors in fungal entomopathogens due to their role in degradation of arthropod cuticle.