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.

The western redcedar genome reveals low genetic diversity in a self-compatible conifer.

We assembled the 9.8-Gbp genome of western redcedar (WRC; Thuja plicata), an ecologically and economically important conifer species of the Cupressaceae. The genome assembly, derived from a uniquely inbred tree produced through five generations of self-fertilization (selfing), was determined to be 86% complete by BUSCO analysis, one of the most complete genome assemblies for a conifer. Population genomic analysis revealed WRC to be one of the most genetically depauperate wild plant species, with an effective population size of approximately 300 and no significant genetic differentiation across its geographic range. Nucleotide diversity, π, is low for a continuous tree species, with many loci showing zero diversity, and the ratio of π at zero- to fourfold degenerate sites is relatively high (approximately 0.33), suggestive of weak purifying selection. Using an array of genetic lines derived from up to five generations of selfing, we explored the relationship between genetic diversity and mating system. Although overall heterozygosity was found to decline faster than expected during selfing, heterozygosity persisted at many loci, and nearly 100 loci were found to deviate from expectations of genetic drift, suggestive of associative overdominance. Nonreference alleles at such loci often harbor deleterious mutations and are rare in natural populations, implying that balanced polymorphisms are maintained by linkage to dominant beneficial alleles. This may account for how WRC remains responsive to natural and artificial selection, despite low genetic diversity.

Genome sequencing of Porostereum spadiceum to study the degradation of levofloxacin.

Despite various plans to rationalize antibiotic use, antibiotic resistance in environmental bacteria is increasing due to the accumulation of antibiotic residues in the environment. This study aimed to test the ability of basidiomycete fungal strains to biotransform the antibiotic levofloxacin, a widely-used third-generation broad-spectrum fluoroquinolone, and to propose enzyme targets potentially involved in this biotransformation. The biotransformation process was performed using fungal strains. Levofloxacin biotransformation reached 100% after 9 days of culture with Porostereum spadiceum BS34. Using genomics and proteomics analyses coupled with activity tests, we showed that P. spadiceum produces several heme-peroxidases together with H2O2-producing enzymes that could be involved in the antibiotic biotransformation process. Using UV and high-resolution mass spectrometry, we were able to detect five levofloxacin degradation products. Their putative identity based on their MS2 fragmentation patterns led to the conclusion that the piperazine moiety was the main target of oxidative modification of levofloxacin by P. spadiceum, leading to a decrease in antibiotic activity.

Genomic Analysis Enlightens Agaricales Lifestyle Evolution and Increasing Peroxidase Diversity.

As actors of global carbon cycle, Agaricomycetes (Basidiomycota) have developed complex enzymatic machineries that allow them to decompose all plant polymers, including lignin. Among them, saprotrophic Agaricales are characterized by an unparalleled diversity of habitats and lifestyles. Comparative analysis of 52 Agaricomycetes genomes (14 of them sequenced de novo) reveals that Agaricales possess a large diversity of hydrolytic and oxidative enzymes for lignocellulose decay. Based on the gene families with the predicted highest evolutionary rates-namely cellulose-binding CBM1, glycoside hydrolase GH43, lytic polysaccharide monooxygenase AA9, class-II peroxidases, glucose-methanol-choline oxidase/dehydrogenases, laccases, and unspecific peroxygenases-we reconstructed the lifestyles of the ancestors that led to the extant lignocellulose-decomposing Agaricomycetes. The changes in the enzymatic toolkit of ancestral Agaricales are correlated with the evolution of their ability to grow not only on wood but also on leaf litter and decayed wood, with grass-litter decomposers as the most recent eco-physiological group. In this context, the above families were analyzed in detail in connection with lifestyle diversity. Peroxidases appear as a central component of the enzymatic toolkit of saprotrophic Agaricomycetes, consistent with their essential role in lignin degradation and high evolutionary rates. This includes not only expansions/losses in peroxidase genes common to other basidiomycetes but also the widespread presence in Agaricales (and Russulales) of new peroxidases types not found in wood-rotting Polyporales, and other Agaricomycetes orders. Therefore, we analyzed the peroxidase evolution in Agaricomycetes by ancestral-sequence reconstruction revealing several major evolutionary pathways and mapped the appearance of the different enzyme types in a time-calibrated species tree.

Diversity of genomic adaptations to the post-fire environment in Pezizales fungi points to crosstalk between charcoal tolerance and sexual development.

Wildfires drastically impact the soil environment, altering the soil organic matter, forming pyrolyzed compounds, and markedly reducing the diversity of microorganisms. Pyrophilous fungi, especially the species from the orders Pezizales and Agaricales, are fire-responsive fungal colonizers of post-fire soil that have historically been found fruiting on burned soil and thus may encode mechanisms of processing these compounds in their genomes. Pyrophilous fungi are diverse. In this work, we explored this diversity and sequenced six new genomes of pyrophilous Pezizales fungi isolated after the 2013 Rim Fire near Yosemite Park in California, USA: Pyronema domesticum, Pyronema omphalodes, Tricharina praecox, Geopyxis carbonaria, Morchella snyderi, and Peziza echinospora. A comparative genomics analysis revealed the enrichment of gene families involved in responses to stress and the degradation of pyrolyzed organic matter. In addition, we found that both protein sequence lengths and G + C content in the third base of codons (GC3) in pyrophilous fungi fall between those in mesophilic/nonpyrophilous and thermophilic fungi. A comparative transcriptome analysis of P. domesticum under two conditions - growing on charcoal, and during sexual development - identified modules of genes that are co-expressed in the charcoal and light-induced sexual development conditions. In addition, environmental sensors such as transcription factors STE12, LreA, LreB, VosA, and EsdC were upregulated in the charcoal condition. Taken together, these results highlight genomic adaptations of pyrophilous fungi and indicate a potential connection between charcoal tolerance and fruiting body formation in P. domesticum.

Evolutionary transition to the ectomycorrhizal habit in the genomes of a hyperdiverse lineage of mushroom-forming fungi.

The ectomycorrhizal (ECM) symbiosis has independently evolved from diverse types of saprotrophic ancestors. In this study, we seek to identify genomic signatures of the transition to the ECM habit within the hyperdiverse Russulaceae. We present comparative analyses of the genomic architecture and the total and secreted gene repertoires of 18 species across the order Russulales, of which 13 are newly sequenced, including a representative of a saprotrophic member of Russulaceae, Gloeopeniophorella convolvens. The genomes of ECM Russulaceae are characterized by a loss of genes for plant cell wall-degrading enzymes (PCWDEs), an expansion of genome size through increased transposable element (TE) content, a reduction in secondary metabolism clusters, and an association of small secreted proteins (SSPs) with TE 'nests', or dense aggregations of TEs. Some PCWDEs have been retained or even expanded, mostly in a species-specific manner. The genome of G. convolvens possesses some characteristics of ECM genomes (e.g. loss of some PCWDEs, TE expansion, reduction in secondary metabolism clusters). Functional specialization in ECM decomposition may drive diversification. Accelerated gene evolution predates the evolution of the ECM habit, indicating that changes in genome architecture and gene content may be necessary to prime the evolutionary switch.

Comparative genomics of pyrophilous fungi reveals a link between fire events and developmental genes.

Forest fires generate a large amount of carbon that remains resident on the site as dead and partially 'pyrolysed' (i.e. burnt) material that has long residency times and constitutes a significant pool in fire-prone ecosystems. In addition, fire-induced hydrophobic soil layers, caused by condensation of pyrolysed waxes and lipids, increase post-fire erosion and can lead to long-term productivity losses. A small set of pyrophilous fungi dominate post-fire soils and are likely to be involved with the degradation of all these compounds, yet almost nothing is currently known about what these fungi do or the metabolic processes they employ. In this study, we sequenced and analysed genomes from fungi isolated after Rim fire near Yosemite National Park in 2013 and showed the enrichment/expansion of CAZymes and families known to be involved in fruiting body initiation when compared to other basidiomycete fungi. We found gene families potentially involved in the degradation of the hydrophobic layer and pyrolysed organic matter, such as hydrophobic surface binding proteins, laccases (AA1_1), xylanases (GH10, GH11), fatty acid desaturases and tannases. Thus, pyrophilous fungi are important actors to restate the soil's functional capabilities.

Intra-species genetic variability drives carbon metabolism and symbiotic host interactions in the ectomycorrhizal fungus Pisolithus microcarpus.

Ectomycorrhizal (ECM) fungi are integral to boreal and temperate forest ecosystem functioning and nutrient cycling. ECM fungi, however, originate from diverse saprotrophic lineages and the impacts of genetic variation across species, and especially within a given ECM species, on function and interactions with the environment is not well understood. Here, we explore the extent of intra-species variation between four isolates of the ECM fungus Pisolithus microcarpus, in terms of gene regulation, carbon metabolism and growth, and interactions with a host, Eucalyptus grandis. We demonstrate that, while a core response to the host is maintained by all of the isolates tested, they have distinct patterns of gene expression and carbon metabolism, resulting in the differential expression of isolate-specific response pathways in the host plant. Together, these results highlight the importance of using a wider range of individuals within a species to understand the broader ecological roles of ECM fungi and their host interactions.

Transcriptome analysis of Aspergillus niger xlnR and xkiA mutants grown on corn Stover and soybean hulls reveals a highly complex regulatory network.

BACKGROUND: Enzymatic plant biomass degradation by fungi is a highly complex process and one of the leading challenges in developing a biobased economy. Some industrial fungi (e.g. Aspergillus niger) have a long history of use with respect to plant biomass degradation and for that reason have become 'model' species for this topic. A. niger is a major industrial enzyme producer that has a broad ability to degrade plant based polysaccharides. A. niger wild-type, the (hemi-)cellulolytic regulator (xlnR) and xylulokinase (xkiA1) mutant strains were grown on a monocot (corn stover, CS) and dicot (soybean hulls, SBH) substrate. The xkiA1 mutant is unable to utilize the pentoses D-xylose and L-arabinose and the polysaccharide xylan, and was previously shown to accumulate inducers for the (hemi-)cellulolytic transcriptional activator XlnR and the arabinanolytic transcriptional activator AraR in the presence of pentoses, resulting in overexpression of their target genes. The xlnR mutant has reduced growth on xylan and down-regulation of its target genes. The mutants therefore have a similar phenotype on xylan, but an opposite transcriptional effect. D-xylose and L-arabinose are the most abundant monosaccharides after D-glucose in nearly all plant-derived biomass materials. In this study we evaluated the effect of the xlnR and xkiA1 mutation during growth on two pentose-rich substrates by transcriptome analysis. RESULTS: Particular attention was given to CAZymes, metabolic pathways and transcription factors related to the plant biomass degradation. Genes coding for the main enzymes involved in plant biomass degradation were down-regulated at the beginning of the growth on CS and SBH. However, at a later time point, significant differences were found in the expression profiles of both mutants on CS compared to SBH. CONCLUSION: This study demonstrates the high complexity of the plant biomass degradation process by fungi, by showing that mutant strains with fairly straightforward phenotypes on pure mono- and polysaccharides, have much less clear-cut phenotypes and transcriptomes on crude plant biomass.

Is Working Memory Necessary for Both Selection and Calculation in Level 1 Visual Perspective Taking? Evidence from Children and Adults.

Previous research has suggested that there are two modular processes of perspective selection and calculation in level 1 perspective taking. Evidence further showed that verbal working memory is associated with both processes in adults. However, research has not tested whether verbal working memory may be associated with working memory in children. Moreover, since perspective taking is associated with spatial working memory, it is necessary to investigate whether spatial working memory links to both processes. By recruiting 9-year-old children and college students in the single-task paradigm and the dual-task paradigm, we conducted two experiments to answer these questions. Results in experiment 1 suggested that verbal working memory correlated with adults' perspective calculation, but with both processes in children's perspective taking. Results in experiment 2 showed that spatial working memory is associated with adults' perspective calculation and children's perspective selection. These findings suggest that different components of working memory play distinct roles in the processing of perspective taking, which is moderated by age.

Disentangling the relationship between state self-compassion and state coping self-efficacy using dynamic structural equation modeling.

The cross-sectional association between self-compassion and coping self-efficacy has been well documented, but little is known about the extent to which self-compassion or coping self-efficacy persists in daily life. This study used dynamic structural equation modeling to explore the temporal relationship between self-compassion and coping self-efficacy through a daily diary study. Participants (N = 240, Mage = 18.98 ± 0.99 years, 44.8% female) completed 14 consecutive daily diaries on self-compassion and coping self-efficacy (for a total of 3219 observations). We found that self-compassion and coping self-efficacy demonstrated stability through autocorrelations. Self-compassion was a significant predictor of subsequent coping self-efficacy, and coping self-efficacy was a significant predictor of subsequent self-compassion. These findings suggest that there has been a virtuous cycle between state coping self-efficacy and state self-compassion.

Neurospora intermedia from a traditional fermented food enables waste-to-food conversion.

Fungal fermentation of food and agricultural by-products holds promise for improving food sustainability and security. However, the molecular basis of fungal waste-to-food upcycling remains poorly understood. Here we use a multi-omics approach to characterize oncom, a fermented food traditionally produced from soymilk by-products in Java, Indonesia. Metagenomic sequencing of samples from small-scale producers in Western Java indicated that the fungus Neurospora intermedia dominates oncom. Further transcriptomic, metabolomic and phylogenomic analysis revealed that oncom-derived N. intermedia utilizes pectin and cellulose degradation during fermentation and belongs to a genetically distinct subpopulation associated with human-generated by-products. Finally, we found that N. intermedia grew on diverse by-products such as fruit and vegetable pomace and plant-based milk waste, did not encode mycotoxins, and could create foods that were positively perceived by consumers outside Indonesia. These results showcase the traditional significance and future potential of fungal fermentation for creating delicious and nutritious foods from readily available by-products.

Seagrass genomes reveal ancient polyploidy and adaptations to the marine environment.

We present chromosome-level genome assemblies from representative species of three independently evolved seagrass lineages: Posidonia oceanica, Cymodocea nodosa, Thalassia testudinum and Zostera marina. We also include a draft genome of Potamogeton acutifolius, belonging to a freshwater sister lineage to Zosteraceae. All seagrass species share an ancient whole-genome triplication, while additional whole-genome duplications were uncovered for C. nodosa, Z. marina and P. acutifolius. Comparative analysis of selected gene families suggests that the transition from submerged-freshwater to submerged-marine environments mainly involved fine-tuning of multiple processes (such as osmoregulation, salinity, light capture, carbon acquisition and temperature) that all had to happen in parallel, probably explaining why adaptation to a marine lifestyle has been exceedingly rare. Major gene losses related to stomata, volatiles, defence and lignification are probably a consequence of the return to the sea rather than the cause of it. These new genomes will accelerate functional studies and solutions, as continuing losses of the 'savannahs of the sea' are of major concern in times of climate change and loss of biodiversity.

Hypoxia activates tumor suppressor p53 by inducing ATR-Chk1 kinase cascade-mediated phosphorylation and consequent 14-3-3γ inactivation of MDMX protein.

It has been known that p53 can be induced and activated by hypoxia, an abnormal condition that often occurs in rapidly growing solid tumors or when normal tissues undergo ischemia. Although the ATR-Chk1 kinase cascade was associated with hypoxia-induced p53 activation, molecules that directly link this hypoxia-ATR-Chk1 pathway to p53 activation have been elusive. Here, we showed that hypoxia could induce phosphorylation of MDMX at Ser-367 and enhance the binding of this phosphorylated MDMX to 14-3-3γ, consequently leading to p53 activation. A Chk1 inhibitor or knockdown of ATR and Chk1 inhibited the phosphorylation of MDMX at Ser-367 and impaired the binding of MDMX to 14-3-3γ in addition to p53 activation in response to hypoxia. In primary mouse embryonic fibroblast cells that harbor a mutant MDMX, including the S367A mutation, hypoxia also failed to induce the binding of this mutant MDMX to 14-3-3γ and to activate p53 and its direct targets. These results demonstrate that hypoxia can activate p53 through inactivation of MDMX by the ATR-Chk1-MDMX-14-3-3γ pathway.

Integrating chromatin conformation information in a self-supervised learning model improves metagenome binning.

Metagenome binning is a key step, downstream of metagenome assembly, to group scaffolds by their genome of origin. Although accurate binning has been achieved on datasets containing multiple samples from the same community, the completeness of binning is often low in datasets with a small number of samples due to a lack of robust species co-abundance information. In this study, we exploited the chromatin conformation information obtained from Hi-C sequencing and developed a new reference-independent algorithm, Metagenome Binning with Abundance and Tetra-nucleotide frequencies-Long Range (metaBAT-LR), to improve the binning completeness of these datasets. This self-supervised algorithm builds a model from a set of high-quality genome bins to predict scaffold pairs that are likely to be derived from the same genome. Then, it applies these predictions to merge incomplete genome bins, as well as recruit unbinned scaffolds. We validated metaBAT-LR's ability to bin-merge and recruit scaffolds on both synthetic and real-world metagenome datasets of varying complexity. Benchmarking against similar software tools suggests that metaBAT-LR uncovers unique bins that were missed by all other methods. MetaBAT-LR is open-source and is available at https://bitbucket.org/project-metabat/metabat-lr.

Population genomics provides insights into the genetic basis of adaptive evolution in the mushroom-forming fungus Lentinula edodes.

Introduction: Mushroom-forming fungi comprise diverse species that develop complex multicellular structures. In cultivated species, both ecological adaptation and artificial selection have driven genome evolution. However, little is known about the connections among genotype, phenotype and adaptation in mushroom-forming fungi. Objectives: This study aimed to (1) uncover the population structure and demographic history of Lentinula edodes, (2) dissect the genetic basis of adaptive evolution in L. edodes, and (3) determine if genes related to fruiting body development are involved in adaptive evolution. Methods: We analyzed genomes and fruiting body-related traits (FBRTs) in 133 L. edodes strains and conducted RNA-seq analysis of fruiting body development in the YS69 strain. Combined methods of genomic scan for divergence, genome-wide association studies (GWAS), and RNA-seq were used to dissect the genetic basis of adaptive evolution. Results: We detected three distinct subgroups of L. edodes via single nucleotide polymorphisms, which showed robust phenotypic and temperature response differentiation and correlation with geographical distribution. Demographic history inference suggests that the subgroups diverged 36,871 generations ago. Moreover, L. edodes cultivars in China may have originated from the vicinity of Northeast China. A total of 942 genes were found to be related to genetic divergence by genomic scan, and 719 genes were identified to be candidates underlying FBRTs by GWAS. Integrating results of genomic scan and GWAS, 80 genes were detected to be related to phenotypic differentiation. A total of 364 genes related to fruiting body development were involved in genetic divergence and phenotypic differentiation. Conclusion: Adaptation to the local environment, especially temperature, triggered genetic divergence and phenotypic differentiation of L. edodes. A general model for genetic divergence and phenotypic differentiation during adaptive evolution in L. edodes, which involves in signal perception and transduction, transcriptional regulation, and fruiting body morphogenesis, was also integrated here.

Dynamic genome evolution in a model fern.

The large size and complexity of most fern genomes have hampered efforts to elucidate fundamental aspects of fern biology and land plant evolution through genome-enabled research. Here we present a chromosomal genome assembly and associated methylome, transcriptome and metabolome analyses for the model fern species Ceratopteris richardii. The assembly reveals a history of remarkably dynamic genome evolution including rapid changes in genome content and structure following the most recent whole-genome duplication approximately 60 million years ago. These changes include massive gene loss, rampant tandem duplications and multiple horizontal gene transfers from bacteria, contributing to the diversification of defence-related gene families. The insertion of transposable elements into introns has led to the large size of the Ceratopteris genome and to exceptionally long genes relative to other plants. Gene family analyses indicate that genes directing seed development were co-opted from those controlling the development of fern sporangia, providing insights into seed plant evolution. Our findings and annotated genome assembly extend the utility of Ceratopteris as a model for investigating and teaching plant biology.

Interplay between MDM2, MDMX, Pirh2 and COP1: the negative regulators of p53.

MDM2, Pirh2 and COP1 are important E3 ubiquitin ligases, which directly interact with p53 and target p53 for proteasome-mediated degradation. MDMX, the MDM2 homologous protein, inhibits p53-mediated transcription activity. The interplay between MDM2, MDMX, Pirh2 and COP1 has not been reported, except the interaction between MDM2 and MDMX. Here, we reported that there were interactions between these four proteins independently of p53. The protein levels of MDM2, MDMX, Pirh2 and COP1 changed when any two of them were co-transfected. Our data also showed that the integrity of MDM2 RING finger domain was crucial for its ability to elevate the protein levels of COP1 and Pirh2. Any two of these four proteins could inhibit p53-mediated transcriptional activity synergistically. Furthermore, COP1 inhibited MDM2 self-ubiquitination and interfered with MDMX ubiquitination by MDM2. Our results suggest that MDM2, MDMX, Pirh2 and COP1 might inhibit p53 activity synergistically in vivo.

YY1 is a novel potential therapeutic target for the treatment of HPV infection-induced cervical cancer by arsenic trioxide.

OBJECTIVE: YY1 is a zinc finger transcription factor involved in the regulation of cell growth, development, and differentiation. Although YY1 can regulate human papillomavirus-type (HPV) viral oncogenes E6 and E7, it remains unknown if YY1 plays a key role in carcinoma progression of HPV-infected cells. Here we sought to determine whether YY1 is upregulated in the cervical cancer tissues and YY1 inhibition contributes to apoptosis of cervical cancer cells, which is at least partly p53 dependent. Therefore, YY1 can be a potential therapeutic target for cervical cancer treatment by arsenic trioxide (As2O3). MATERIALS AND METHODS: The expression level of YY1 was examined and analyzed by Western blot in pathologically confirmed primary cervical cancer samples, in the adjacent normal samples, as well as in normal cervix samples. The effects of YY1 inhibition by specific small interfering RNA in HeLa cells were determined by Western blot analysis of p53 level, cell growth curve, colony formation assay, and apoptosis. The contribution of YY1 to As2O3-induced p53 activation and apoptosis was also examined by Western blot and cell cycle analysis. RESULTS: Here we report that the expression level of YY1 is significantly elevated in the primary cancer tissues. In HPV-positive HeLa cells, small interfering RNA-mediated YY1 inhibition induced apoptosis and increased the expression of p53. Treatment of HeLa cells with As2O3, a known anti-cervical cancer agent, reduced both protein and mRNA levels of YY1 in HeLa cells. YY1 knockdown significantly further enhanced As2O3-induced apoptosis. CONCLUSIONS: These results demonstrated that the expression of YY1 is upregulated in cervical carcinomas and that YY1 plays a critical role in the progression of HPV-positive cervical cancer. In addition, YY1 inhibition induces p53 activation and apoptosis in HPV-infected HeLa cells. Thus, YY1 is an As2O3 target and could serve as a potential drug sensitizer for anti-cervical cancer therapy.