A global phylogenomic analysis of the shiitake genus Lentinula.

Lentinula is a broadly distributed group of fungi that contains the cultivated shiitake mushroom, L. edodes. We sequenced 24 genomes representing eight described species and several unnamed lineages of Lentinula from 15 countries on four continents. Lentinula comprises four major clades that arose in the Oligocene, three in the Americas and one in Asia-Australasia. To expand sampling of shiitake mushrooms, we assembled 60 genomes of L. edodes from China that were previously published as raw Illumina reads and added them to our dataset. Lentinula edodes sensu lato (s. lat.) contains three lineages that may warrant recognition as species, one including a single isolate from Nepal that is the sister group to the rest of L. edodes s. lat., a second with 20 cultivars and 12 wild isolates from China, Japan, Korea, and the Russian Far East, and a third with 28 wild isolates from China, Thailand, and Vietnam. Two additional lineages in China have arisen by hybridization among the second and third groups. Genes encoding cysteine sulfoxide lyase (lecsl) and γ-glutamyl transpeptidase (leggt), which are implicated in biosynthesis of the organosulfur flavor compound lenthionine, have diversified in Lentinula. Paralogs of both genes that are unique to Lentinula (lecsl 3 and leggt 5b) are coordinately up-regulated in fruiting bodies of L. edodes. The pangenome of L. edodes s. lat. contains 20,308 groups of orthologous genes, but only 6,438 orthogroups (32%) are shared among all strains, whereas 3,444 orthogroups (17%) are found only in wild populations, which should be targeted for conservation.

Evolutionary innovations through gain and loss of genes in the ectomycorrhizal Boletales.

We aimed to identify genomic traits of transitions to ectomycorrhizal ecology within the Boletales by comparing the genomes of 21 symbiotrophic species with their saprotrophic brown-rot relatives. Gene duplication rate is constant along the backbone of Boletales phylogeny with large loss events in several lineages, while gene family expansion sharply increased in the late Miocene, mostly in the Boletaceae. Ectomycorrhizal Boletales have a reduced set of plant cell-wall-degrading enzymes (PCWDEs) compared with their brown-rot relatives. However, the various lineages retain distinct sets of PCWDEs, suggesting that, over their evolutionary history, symbiotic Boletales have become functionally diverse. A smaller PCWDE repertoire was found in Sclerodermatineae. The gene repertoire of several lignocellulose oxidoreductases (e.g. laccases) is similar in brown-rot and ectomycorrhizal species, suggesting that symbiotic Boletales are capable of mild lignocellulose decomposition. Transposable element (TE) proliferation contributed to the higher evolutionary rate of genes encoding effector-like small secreted proteins, proteases, and lipases. On the other hand, we showed that the loss of secreted CAZymes was not related to TE activity but to DNA decay. This study provides novel insights on our understanding of the mechanisms influencing the evolutionary diversification of symbiotic boletes.

Gene family expansions and transcriptome signatures uncover fungal adaptations to wood decay.

Because they comprise some of the most efficient wood-decayers, Polyporales fungi impact carbon cycling in forest environment. Despite continuous discoveries on the enzymatic machinery involved in wood decomposition, the vision on their evolutionary adaptation to wood decay and genome diversity remains incomplete. We combined the genome sequence information from 50 Polyporales species, including 26 newly sequenced genomes and sought for genomic and functional adaptations to wood decay through the analysis of genome composition and transcriptome responses to different carbon sources. The genomes of Polyporales from different phylogenetic clades showed poor conservation in macrosynteny, indicative of genome rearrangements. We observed different gene family expansion/contraction histories for plant cell wall degrading enzymes in core polyporoids and phlebioids and captured expansions for genes involved in signalling and regulation in the lineages of white rotters. Furthermore, we identified conserved cupredoxins, thaumatin-like proteins and lytic polysaccharide monooxygenases with a yet uncharacterized appended module as new candidate players in wood decomposition. Given the current need for enzymatic toolkits dedicated to the transformation of renewable carbon sources, the observed genomic diversity among Polyporales strengthens the relevance of mining Polyporales biodiversity to understand the molecular mechanisms of wood decay.

Novel Monomeric Fungal Subtilisin Inhibitor from a Plant-Pathogenic Fungus, Choanephora cucurbitarum: Isolation and Molecular Characterization.

The bacterial protease inhibitor domains known as Streptomyces subtilisin inhibitors (SSI) are rarely found in fungi. Genome analysis of a fungal pathogen, Choanephora cucurbitarum KUS-F28377, revealed 11 SSI-like domains that are horizontally transferred and sequentially diverged during evolution. We investigated the molecular function of fungal SSI-like domains of C. cucurbitarum, designated "choanepins." Among the proteins tested, only choanepin9 showed inhibitory activity against subtilisin as the target protease, accounting for 47% of the inhibitory activity of bacterial SSI. However, the binding affinity (expressed as the dissociation constant [Kd ]) of choanepin9 measured via microscale thermophoresis was 21 nM, whereas that for bacterial SSI is 34 nM. The trend of binding and inhibitory activity suggests that the two inhibitors exhibit different inhibitory mechanisms for subtilisin protease. Interestingly, choanepin9 was identified as a monomer in studies in vitro, whereas bacterial SSI is a homodimer. Based on these observations, we constructed a monomeric bacterial SSI protein with decreased binding affinity to abrogate its inhibitory activity. By altering the reactive sites of choanepin9 deduced from the P1 and P4 sites of bacterial SSI, we reestablished that these residues in choanepins are also crucial for modulating inhibitory activity. These findings suggest that the fungal SSI evolved to target specific cognate proteases by altering the residues involved in inhibitory reactivity (reactive sites) and binding affinity (structural integrity). The function of fungal SSI proteins identified in this study provides not only a clue to fungal pathogenesis via protease inhibition but also a template for the design of novel serine protease inhibitors.IMPORTANCE Until recently, Streptomyces subtilisin inhibitors (SSI) were reported and characterized only in bacteria. We found SSI-like domains in a plant-pathogenic fungus, Choanephora cucurbitarum KUS-F28377, which contains 11 sequentially diverged SSI-like domains. None of these fungal SSI-like domains were functionally characterized before. The active form of fungal SSI-like protein is a monomer, in contrast to the homodimeric bacterial SSI. We constructed a synthetic monomer of bacterial SSI to demonstrate the modulation of its activity based on structural integrity and not reactive sites. Our results suggest the duplication and divergence of SSI-like domains of C. cucurbitarum within the genome to inhibit various cognate proteases during evolution by modulating both binding and reactivity. The molecular functional characterization of fungal SSI-like domains will be useful in understanding their biological role and future biotechnological applications.

Comprehensive genomic and transcriptomic analysis of polycyclic aromatic hydrocarbon degradation by a mycoremediation fungus, Dentipellis sp. KUC8613.

The environmental accumulation of polycyclic aromatic hydrocarbons (PAHs) is of great concern due to potential carcinogenic and mutagenic risks, as well as their resistance to remediation. While many fungi have been reported to break down PAHs in environments, the details of gene-based metabolic pathways are not yet comprehensively understood. Specifically, the genome-scale transcriptional responses of fungal PAH degradation have rarely been reported. In this study, we report the genomic and transcriptomic basis of PAH bioremediation by a potent fungal degrader, Dentipellis sp. KUC8613. The genome size of this fungus was 36.71 Mbp long encoding 14,320 putative protein-coding genes. The strain efficiently removed more than 90% of 100 mg/l concentration of PAHs within 10 days. The genomic and transcriptomic analysis of this white rot fungus highlights that the strain primarily utilized non-ligninolytic enzymes to remove various PAHs, rather than typical ligninolytic enzymes known for playing important roles in PAH degradation. PAH removal by non-ligninolytic enzymes was initiated by both different PAH-specific and common upregulation of P450s, followed by downstream PAH-transforming enzymes such as epoxide hydrolases, dehydrogenases, FAD-dependent monooxygenases, dioxygenases, and glycosyl- or glutathione transferases. Among the various PAHs, phenanthrene induced a more dynamic transcriptomic response possibly due to its greater cytotoxicity, leading to highly upregulated genes involved in the translocation of PAHs, a defense system against reactive oxygen species, and ATP synthesis. Our genomic and transcriptomic data provide a foundation of understanding regarding the mycoremediation of PAHs and the application of this strain for polluted environments.

Genomic discovery of the hypsin gene and biosynthetic pathways for terpenoids in Hypsizygus marmoreus.

BACKGROUND: Hypsizygus marmoreus (Beech mushroom) is a popular ingredient in Asian cuisine. The medicinal effects of its bioactive compounds such as hypsin and hypsiziprenol have been reported, but the genetic basis or biosynthesis of these components is unknown. RESULTS: In this study, we sequenced a reference strain of H. marmoreus (Haemi 51,987-8). We evaluated various assembly strategies, and as a result the Allpaths and PBJelly produced the best assembly. The resulting genome was 42.7 Mbp in length and annotated with 16,627 gene models. A putative gene (Hypma_04324) encoding the antifungal and antiproliferative hypsin protein with 75% sequence identity with the previously known N-terminal sequence was identified. Carbohydrate active enzyme analysis displayed the typical feature of white-rot fungi where auxiliary activity and carbohydrate-binding modules were enriched. The genome annotation revealed four terpene synthase genes responsible for terpenoid biosynthesis. From the gene tree analysis, we identified that terpene synthase genes can be classified into six clades. Four terpene synthase genes of H. marmoreus belonged to four different groups that implies they may be involved in the synthesis of different structures of terpenes. A terpene synthase gene cluster was well-conserved in Agaricomycetes genomes, which contained known biosynthesis and regulatory genes. CONCLUSIONS: Genome sequence analysis of this mushroom led to the discovery of the hypsin gene. Comparative genome analysis revealed the conserved gene cluster for terpenoid biosynthesis in the genome. These discoveries will further our understanding of the biosynthesis of medicinal bioactive molecules in this edible mushroom.

FunGAP: Fungal Genome Annotation Pipeline using evidence-based gene model evaluation.

MOTIVATION: Successful genome analysis depends on the quality of gene prediction. Although fungal genome sequencing and assembly have become trivial, its annotation procedure has not been standardized yet. RESULTS: FunGAP predicts protein-coding genes in a fungal genome assembly. To attain high-quality gene models, this program runs multiple gene predictors, evaluates all predicted genes, and assembles gene models that are highly supported by homology to known sequences. To do this, we built a scoring function to estimate the congruency of each gene model based on known protein or domain homology. AVAILABILITY AND IMPLEMENTATION: FunGAP is written in Python script and is available in GitHub ( https://github.com/CompSynBioLab-KoreaUniv/FunGAP ). This software is freely available only for noncommercial users. CONTACT: igchoi@korea.ac.kr. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The obligate alkalophilic soda-lake fungus Sodiomyces alkalinus has shifted to a protein diet.

Sodiomyces alkalinus is one of the very few alkalophilic fungi, adapted to grow optimally at high pH. It is widely distributed at the plant-deprived edges of extremely alkaline lakes and locally abundant. We sequenced the genome of S. alkalinus and reconstructed evolution of catabolic enzymes, using a phylogenomic comparison. We found that the genome of S. alkalinus is larger, but its predicted proteome is smaller and heavily depleted of both plant-degrading enzymes and proteinases, when compared to its closest plant-pathogenic relatives. Interestingly, despite overall losses, S. alkalinus has retained many proteinases families and acquired bacterial cell wall-degrading enzymes, some of them via horizontal gene transfer from bacteria. This fungus has very potent proteolytic activity at high pH values, but slowly induced low activity of cellulases and hemicellulases. Our experimental and in silico data suggest that plant biomass, a common food source for most fungi, is not a preferred substrate for S. alkalinus in its natural environment. We conclude that the fungus has abandoned the ancestral plant-based diet and has become specialized in a more protein-rich food, abundantly available in soda lakes in the form of prokaryotes and small crustaceans.

Genomics of wood-degrading fungi.

Woody plants convert the energy of the sun into lignocellulosic biomass, which is an abundant substrate for bioenergy production. Fungi, especially wood decayers from the class Agaricomycetes, have evolved ways to degrade lignocellulose into its monomeric constituents, and understanding this process may facilitate the development of biofuels. Over the past decade genomics has become a powerful tool to study the Agaricomycetes. In 2004 the first sequenced genome of the white rot fungus Phanerochaete chrysosporium revealed a rich catalog of lignocellulolytic enzymes. In the decade that followed the number of genomes of Agaricomycetes grew to more than 75 and revealed a diversity of wood-decaying strategies. New technologies for high-throughput functional genomics are now needed to further study these organisms.

Transcriptomics of Temporal- versus Substrate-Specific Wood Decay in the Brown-Rot Fungus Fibroporia radiculosa.

Brown-rot fungi lack many enzymes associated with complete wood degradation, such as lignin-attacking peroxidases, and have developed alternative mechanisms for rapid wood breakdown. To identify the effects of culture conditions and wood substrates on gene expression, we grew Fibroporia radiculosa in submerged cultures containing Wiley milled wood (5 days) and solid wood wafers (30 days), using aspen, pine, and spruce as a substrate. The comparative analysis revealed that wood species had a limited effect on the transcriptome: <3% of genes were differentially expressed between different wood species substrates. The comparison between gene expression during growth on milled wood and wood wafer conditions, however, indicated that the genes encoding plant cell wall-degrading enzymes, such as glycoside hydrolases and peptidases, were activated during growth on wood wafers, confirming previous reports. On the other hand, it was shown for the first time that the genes encoding Fenton chemistry enzymes, such as hydroquinone biosynthesis enzymes and oxidoreductases, were activated during submerged growth on ground wood. This illustrates the diversity of wood-decay reactions encoded in fungi and activated at different stages of this process.

Draft Genome Sequence of the Reference Strain of the Korean Medicinal Mushroom Wolfiporia cocos KMCC03342.

Wolfiporia cocos is a wood-decay brown rot fungus belonging to the family Polyporaceae. While the fungus grows, the sclerotium body of the strain, dubbed Bokryeong in Korean, is formed around the roots of conifer trees. The dried sclerotium has been widely used as a key component of many medicinal recipes in East Asia. Wolfiporia cocos strain KMCC03342 is the reference strain registered and maintained by the Korea Seed and Variety Service for commercial uses. Here, we present the first draft genome sequence of W. cocos KMCC03342 using a hybrid assembly technique combining both short- and long-read sequences. The genome has a total length of 55.5 Mb comprised of 343 contigs with N50 of 332 kb and 95.8% BUSCO completeness. The GC ratio was 52.2%. We predicted 14,296 protein-coding gene models based on ab initio gene prediction and evidence-based annotation procedure using RNAseq data. The annotated genome was predicted to have 19 terpene biosynthesis gene clusters, which was the same number as the previously sequenced W. cocos strain MD-104 genome but higher than Chinese W. cocos strains. The genome sequence and the predicted gene clusters allow us to study biosynthetic pathways for the active ingredients of W. cocos.

Unusual genome expansion and transcription suppression in ectomycorrhizal Tricholoma matsutake by insertions of transposable elements.

Genome sequencing of Tricholoma matsutake revealed its unusually large size as 189.0 Mbp, which is a consequence of extraordinarily high transposable element (TE) content. We identified that 702 genes were surrounded by TEs, and 83.2% of these genes were not transcribed at any developmental stage. This observation indicated that the insertion of TEs alters the transcription of the genes neighboring these TEs. Repeat-induced point mutation, such as C to T hypermutation with a bias over "CpG" dinucleotides, was also recognized in this genome, representing a typical defense mechanism against TEs during evolution. Many transcription factor genes were activated in both the primordia and fruiting body stages, which indicates that many regulatory processes are shared during the developmental stages. Small secreted protein genes (<300 aa) were dominantly transcribed in the hyphae, where symbiotic interactions occur with the hosts. Comparative analysis with 37 Agaricomycetes genomes revealed that IstB-like domains (PF01695) were conserved across taxonomically diverse mycorrhizal genomes, where the T. matsutake genome contained four copies of this domain. Three of the IstB-like genes were overexpressed in the hyphae. Similar to other ectomycorrhizal genomes, the CAZyme gene set was reduced in T. matsutake, including losses in the glycoside hydrolase genes. The T. matsutake genome sequence provides insight into the causes and consequences of genome size inflation.

Practical Guide for Fungal Gene Prediction from Genome Assembly and RNA-Seq Reads by FunGAP.

FunGAP is a Python-wrapped fungal genome annotation pipeline running under the Linux/Unix operating system. The annotation procedure used in FunGAP requires two inputs, genome assembly and RNA-seq reads. FunGAP aims to predict the most feasible gene from all plausible gene models obtained from various gene prediction programs using multiple strategies such as ab initio, EST-, and/or homology-based methods. This guide covers how to run the FunGAP from the command line and use various options for practical gene prediction. Users can choose options for quality control of the input sequences, selecting model database, filtration of predicted gene models, and post-process such as checking genome completeness and transposable elements. Using FunGAP, the user will acquire a high-quality fungal gene prediction for post-genome sequencing analysis.

Prokaryotic Contig Annotation Pipeline Server: Web Application for a Prokaryotic Genome Annotation Pipeline Based on the Shiny App Package.

Genome annotation is a primary step in genomic research. To establish a light and portable prokaryotic genome annotation pipeline for use in individual laboratories, we developed a Shiny app package designated as "P-CAPS" (Prokaryotic Contig Annotation Pipeline Server). The package is composed of R and Python scripts that integrate publicly available annotation programs into a server application. P-CAPS is not only a browser-based interactive application but also a distributable Shiny app package that can be installed on any personal computer. The final annotation is provided in various standard formats and is summarized in an R markdown document. Annotation can be visualized and examined with a public genome browser. A benchmark test showed that the annotation quality and completeness of P-CAPS were reliable and compatible with those of currently available public pipelines.

Genome Analysis of a Zygomycete Fungus Choanephora cucurbitarum Elucidates Necrotrophic Features Including Bacterial Genes Related to Plant Colonization.

A zygomycete fungus, Choanephora cucurbitarum is a plant pathogen that causes blossom rot in cucurbits and other plants. Here we report the genome sequence of Choanephora cucurbitarum KUS-F28377 isolated from squash. The assembled genome has a size of 29.1 Mbp and 11,977 protein-coding genes. The genome analysis indicated that C. cucurbitarum may employ a plant pathogenic mechanism similar to that of bacterial plant pathogens. The genome contained 11 genes with a Streptomyces subtilisin inhibitor-like domain, which plays an important role in the defense against plant immunity. This domain has been found only in bacterial genomes. Carbohydrate active enzyme analysis detected 312 CAZymes in this genome where carbohydrate esterase family 6, rarely found in dikaryotic fungal genomes, was comparatively enriched. The comparative genome analysis showed that the genes related to sexual communication such as the biosynthesis of ß-carotene and trisporic acid were conserved and diverged during the evolution of zygomycete genomes. Overall, these findings will help us to understand how zygomycetes are associated with plants.

Genome sequence of a white rot fungus Schizopora paradoxa KUC8140 for wood decay and mycoremediation.

Schizopora paradoxa KUC8140 is a white rot wood degrader commonly found in Korea. Tolerance to heavy metals and polycyclic aromatic hydrocarbons and dye decolorization activity make this strain a potential candidate for mycoremediation. We report the genome sequence of S. paradoxa KUC8140 containing 44.4Mbp. Based on ab initio gene prediction, homology search and RNA-seq, total 17,098 gene models were annotated. We identified 17 lignin-modifying peroxidases and other 377 carbohydrate-active enzymes for modeling lignocellulose deconstruction and mycoremediation.

Metabolically engineered glucose-utilizing Shewanella strains under anaerobic conditions.

Comparative genome analysis of Shewanella strains predicted that the strains metabolize preferably two- and three-carbon carbohydrates as carbon/electron source because many Shewanella genomes are deficient of the key enzymes in glycolysis (e.g., glucokinase). In addition, all Shewanella genomes are known to have only one set of genes associated with the phosphotransferase system required to uptake sugars. To engineer Shewanella strains that can utilize five- and six-carbon carbohydrates, we constructed glucose-utilizing Shewanella oneidensis MR-1 by introducing the glucose facilitator (glf; ZMO0366) and glucokinase (glk; ZMO0369) genes of Zymomonas mobilis. The engineered MR-1 strain was able to grow on glucose as a sole carbon/electron source under anaerobic conditions. The glucose affinity (Ks) and glucokinase activity in the engineered MR-1 strain were 299.46 mM and 0.259 ± 0.034 U/g proteins. The engineered strain was successfully applied to a microbial fuel cell system and exhibited current generation using glucose as the electron source.