Chromosome-level genome map provides insights into diverse defense mechanisms in the medicinal fungus Ganoderma sinense

Whole-genome sequencing and comparative genomic analyses of the medicinal fungus Sanguinoderma infundibulare in Ganodermataceae

Sanguinoderma infundibulare is a newly discovered species of Ganodermataceae known to have high medicinal and ecological values. In this study, the whole-genome sequencing and comparative genomic analyses were conducted to further understand Ganodermataceae's genomic structural and functional characteristics. Using the Illumina NovaSeq and PacBio Sequel platforms, 88 scaffolds were assembled to obtain a 48.99-Mb high-quality genome of S. infundibulare. A total of 14,146 protein-coding genes were annotated in the whole genome, with 98.6% of complete benchmarking universal single-copy orthologs (BUSCO) scores. Comparative genomic analyses were conducted among S. infundibulare, Sanguinoderma rugosum, Ganoderma lucidum, and Ganoderma sinense to determine their intergeneric differences. The 4 species were found to share 4,011 orthogroups, and 24 specific gene families were detected in the genus Sanguinoderma. The gene families associated with carbohydrate esterase in S. infundibulare were significantly abundant, which was reported to be involved in hemicellulose degradation. One specific gene family in Sanguinoderma was annotated with siroheme synthase, which may be related to the typical characteristics of fresh pore surface changing to blood red when bruised. This study enriched the available genome data for the genus Sanguinoderma, elucidated the differences between Ganoderma and Sanguinoderma, and provided insights into the characteristics of the genome structure and function of S. infundibulare.

Telomere-to-Telomere Genome Assembly of Tibetan Medicinal Mushroom Ganoderma leucocontextum and the First Copia Centromeric Retrotransposon in Macro-Fungi Genome

A complete telomere-to-telomere (T2T) genome has been a longstanding goal in the field of genomic research. By integrating high-coverage and precise long-read sequencing data using multiple assembly strategies, we present here the first T2T gap-free genome assembly of Ganoderma leucocontextum strain GL72, a Tibetan medicinal mushroom. The T2T genome, with a size of 46.69 Mb, consists 13 complete nuclear chromosomes and typical telomeric repeats (CCCTAA)n were detected at both ends of 13 chromosomes. The high mapping rate, uniform genome coverage, a complete BUSCOs of 99.7%, and base accuracy exceeding 99.999% indicate that this assembly represents the highest level of completeness and quality. Regions characterized by distinct structural attributes, including highest Hi-C interaction intensity, high repeat content, decreased gene density, low GC content, and minimal or no transcription levels across all chromosomes may represent potential centromeres. Sequence analysis revealed the first Copia centromeric retrotransposon in macro-fungi genome. Phylogenomic analysis identified that G. leucocontextum and G. tsugae diverged from the other Ganoderma species approximately 9.8-17.9 MYA. The prediction of secondary metabolic clusters confirmed the capability of this fungus to produce a substantial quantity of metabolites. This T2T gap-free genome will contribute to the genomic 'dark matter' elucidation and server as a great reference for genetics, genomics, and evolutionary studies of G. leucocontextum.

Unraveling the Polysaccharide Biosynthesis Potential of Ganoderma lucidum: A Chromosome-Level Assembly Using Hi-C Sequencing

Ganoderma lucidum exhibits the ability to synthesize a diverse range of biologically active molecules with significant pharmaceutical potential, including xylomannan and fucogalactan, which have demonstrated antitumor activity. However, there exists considerable intra-species variability in the capacity to produce these metabolites at high concentrations, likely reflecting the high genomic diversity observed from a limited number of strains sequenced to date. We employed high-throughput shotgun sequencing to obtain the complete genome sequence of G. lucidum strain 5.1, which is distinguished by its remarkable xylomannan synthesis capabilities. Through the utilization of semi-automatic reordering based on conformation capture (Hi-C) data, we substantially enhanced the assembly process, resulting in the generation of 12 chromosome-level scaffolds with a cumulative length of 39 Mbp. By employing both de novo and homology-based approaches, we performed comprehensive annotation of the genome, thereby identifying a diverse repertoire of genes likely involved in polysaccharide biosynthesis. The genome sequence generated in this study serves as a valuable resource for elucidating the molecular mechanisms underlying the medicinal potential of Ganoderma species, discovering novel pharmaceutically valuable compounds, and elucidating the ecological mechanisms of the species. Furthermore, the chromosome contact map obtained for the first time for this species extends our understanding of 3D fungal genomics and provides insights into the functional and structural organization within the fungal kingdom.

Integrated omic profiling of the medicinal mushroom Inonotus obliquus under submerged conditions

BACKGROUND

The Inonotus obliquus mushroom, a wondrous fungus boasting edible and medicinal qualities, has been widely used as a folk medicine and shown to have many potential pharmacological secondary metabolites. The purpose of this study was to supply a global landscape of genome-based integrated omic analysis of the fungus under lab-growth conditions.

RESULTS

This study presented a genome with high accuracy and completeness using the Pacbio Sequel II third-generation sequencing method. The de novo assembled fungal genome was 36.13 Mb, and contained 8352 predicted protein-coding genes, of which 365 carbohydrate-active enzyme (CAZyme)-coding genes and 19 biosynthetic gene clusters (BCGs) for secondary metabolites were identified. Comparative transcriptomic and proteomic analysis revealed a global view of differential metabolic change between seed and fermentation culture, and demonstrated positive correlations between transcription and expression levels of 157 differentially expressed genes involved in the metabolism of amino acids, fatty acids, secondary metabolites, antioxidant and immune responses. Facilitated by the widely targeted metabolomic approach, a total of 307 secondary substances were identified and quantified, with a significant increase in the production of antioxidant polyphenols.

CONCLUSION

This study provided the comprehensive analysis of the fungus Inonotus obliquus, and supplied fundamental information for further screening of promising target metabolites and exploring the link between the genome and metabolites.

High-quality genome assembly and multi-omics analysis of pigment synthesis pathway in Auricularia cornea

Owing to its great market potential for food and health care, white Auricularia cornea, a rare edible fungus, has received increased attention in recent years. This study presents a high-quality genome assembly of A. cornea and multi-omics analysis of its pigment synthesis pathway. Continuous Long Reads libraries, combined with Hi-C-assisted assembly were used to assemble of white A. cornea. Based on this data, we analyzed the transcriptome and metabolome of purple and white strains during the mycelium, primordium, and fruiting body stages. Finally, we obtained the genome of A.cornea assembled from 13 clusters. Comparative and evolutionary analysis suggests that A.cornea is more closely related to Auricularia subglabra than to Auricularia heimuer. The divergence of white/purple A.cornea occurred approximately 40,000 years ago, and there were numerous inversions and translocations between homologous regions of the two genomes. Purple strain synthesized pigment via the shikimate pathway. The pigment in the fruiting body of A. cornea was γ-glutaminyl-3,4-dihydroxy-benzoate. During pigment synthesis, α-D-glucose-1P, citrate, 2-Oxoglutarate, and glutamate were four important intermediate metabolites, whereas polyphenol oxidase and other 20 enzyme genes were the key enzymes. This study sheds light on the genetic blueprint and evolutionary history of the white A.cornea genome, revealing the mechanism of pigment synthesis in A.cornea. It has important theoretical and practical implications for understanding the evolution of basidiomycetes, molecular breeding of white A.cornea, and deciphering the genetic regulations of edible fungi. Additionally, it provides valuable insights for the study of phenotypic traits in other edible fungi.

Increasing the production of the bioactive compounds in medicinal mushrooms: an omics perspective

Macroscopic fungi, mainly higher basidiomycetes and some ascomycetes, are considered medicinal mushrooms and have long been used in different areas due to their pharmaceutically/nutritionally valuable bioactive compounds. However, the low production of these bioactive metabolites considerably limits the utilization of medicinal mushrooms both in commerce and clinical trials. As a result, many attempts, ranging from conventional methods to novel approaches, have been made to improve their production. The novel strategies include conducting omics investigations, constructing genome-scale metabolic models, and metabolic engineering. So far, genomics and the combined use of different omics studies are the most utilized omics analyses in medicinal mushroom research (both with 31% contribution), while metabolomics (with 4% contribution) is the least. This article is the first attempt for reviewing omics investigations in medicinal mushrooms with the ultimate aim of bioactive compound overproduction. In this regard, the role of these studies and systems biology in elucidating biosynthetic pathways of bioactive compounds and their contribution to metabolic engineering will be highlighted. Also, limitations of omics investigations and strategies for overcoming them will be provided in order to facilitate the overproduction of valuable bioactive metabolites in these valuable organisms.

Comprehensive analysis of codon bias in 13 Ganoderma mitochondrial genomes

Introduction

Codon usage bias is a prevalent phenomenon observed across various species and genes. However, the specific attributes of codon usage in the mitochondrial genome of Ganoderma species remain unknown.

Methods

In this study, we investigated the codon bias of 12 mitochondrial core protein-coding genes (PCGs) in 9 Ganoderma species, including 13 Ganoderma strains.

Results

The codons of all Ganoderma strains showed a preference for ending in A/T. Additionally, correlations between codon base composition and the codon adaptation index (CAI), codon bias index (CBI) and frequency of optimal codons (FOP) were identified, demonstrating the impact of base composition on codon bias. Various base bias indicators were found to vary between or within Ganoderma strains, including GC3s, the CAI, the CBI, and the FOP. The results also revealed that the mitochondrial core PCGs of Ganoderma have an average effective number of codons (ENC) lower than 35, indicating strong bias toward certain codons. Evidence from neutrality plot and PR2-bias plot analysis indicates that natural selection is a major factor affecting codon bias in Ganoderma. Additionally, 11 to 22 optimal codons (ΔRSCU>0.08 and RSCU>1) were identified in 13 Ganoderma strains, with GCA, AUC, and UUC being the most widely used optimal codons in Ganoderma. By analyzing the combined mitochondrial sequences and relative synonymous codon usage (RSCU) values, the genetic relationships between or within Ganoderma strains were determined, indicating variations between them. Nevertheless, RSCU-based analysis illustrated the intra- and interspecies relationships of certain Ganoderma species.

Discussion

This study deepens our insight into the synonymous codon usage characteristics, genetics, and evolution of this important fungal group.

Taxonomy and phylogeny of Sanguinoderma rugosum complex with descriptions of a new species and a new combination

Sanguinoderma is distributed in tropical and subtropical areas as a member of Amauroderma s. lat., and the economic values of Sanguinoderma led to high attention in the taxonomic studies. Previously, 16 species have been developed into Sanguinoderma. In this study, the taxonomic system of Sanguinoderma was reconducted based on morphological and multi-gene phylogenetic analyses, especially making a distinction for Sanguinoderma rugosum complex. Morphological analysis was based on the notes of macro- and micro morphological observations. Multi-gene phylogenetic analyses were used maximum likelihood (ML) and Bayesian inference (BI) analyses inferred from combined dataset of ITS, nLSU, rpb2, tef1, mtSSU, and nSSU. Combined with morphological characters and phylogenetic evidence, the results demonstrated that S. rugosum complex consists of five taxa, in which Sanguinoderma leucomarginatum was described as a new species, and it is characterized by the orbicular pilei with white to buff margin when fresh and clavate apical cells of pileipellis with septa. In addition, Amauroderma preussii was transferred to Sanguinoderma as a new combination due to its blood-red color-changed pore surface; it is characterized by the funnel-shaped, greyish brown, and glabrous pilei with strongly incurved margin. Detailed descriptions and photographs of the two species were provided. With the extension of this study, 18 species were accepted in Sanguinoderma, and 12 species among them were distributed in China. A key to accepted species of Sanguinoderma was also provided.

Massive genome investigations reveal insights of prevalent introgression for environmental adaptation and triterpene biosynthesis in Ganoderma

Genome introgression is one of the driving forces that can increase species and genetic diversity and facilitate the adaptive evolution of organisms and biodiversity conservation. However, the genomic introgression and its contribution to biodiversity of macrofungi are still unclear. The genus Ganoderma is a typical macrofungal group that plays crucial roles in forest ecosystem as saprophytic organisms and plant pathogens, and is also involved in human health as medicinal mushrooms. Most public Ganoderma genomes are fragmented, and reference genomes and whole-genome information of diverse germplasm resources for many Ganoderma species are lacking, thus hindering functional and evolutionary genomic investigations among Ganoderma species. In this study, we provide high-quality genomes of 10 Ganoderma species and whole-genome variants data of 224 individuals from various ecoregions, enabling us to infer the phylogeny of Ganoderma species and their historical population dynamics. Based on whole-genome variants, widespread and genome-wide introgression among Ganoderma species is revealed. Genes with significant introgression signals were related to stress response, digestive absorption, and secondary metabolite synthesis, factors that may contribute to environmental adaptation and important biocomponent metabolism. CYP512U6, an essential functional gene in the CYP450 family related to Ganoderma triterpene synthesis, was detected with significant introgression and selection signals combined with Ganoderma metabolomic analysis, indicating that both ancient gene exchange and recent domestication have contributed to the categories and content of secondary metabolites of Ganoderma. The reference genomes, whole-genome variants, and metabolite profiles could serve as abundant and valuable genetic resources for evolution, ecology, and conservation investigations of Ganoderma species and other macrofungi.

Comparative Mitogenomic Analysis Reveals Intraspecific, Interspecific Variations and Genetic Diversity of Medical Fungus Ganoderma

Ganoderma species are widely distributed in the world with high diversity. Some species are considered to be pathogenic fungi while others are used as traditional medicine in Asia. In this study, we sequenced and assembled four Ganoderma complete mitogenomes, including G. subamboinense s118, G. lucidum s37, G. lingzhi s62, and G. lingzhi s74. The sizes of the four mitogenomes ranged from 50,603 to 73,416 bp. All Ganoderma specimens had a full set of core protein-coding genes (PCGs), and the rps3 gene of Ganoderma species was detected to be under positive or relaxed selection. We found that the non-conserved PCGs, which encode RNA polymerases, DNA polymerases, homing endonucleases, and unknown functional proteins, are dynamic within and between Ganoderma species. Introns were thought to be the main contributing factor in Ganoderma mitogenome size variation (p < 0.01). Frequent intron loss/gain events were detected within and between Ganoderma species. The mitogenome of G. lucidum s26 gained intron P637 in the cox3 gene compared with the other two G. lucidum mitogenomes. In addition, some rare introns in Ganoderma were detected in distinct Basidiomycetes, indicating potential gene transfer events. Comparative mitogenomic analysis revealed that gene arrangements also varied within and between Ganoderma mitogenomes. Using maximum likelihood and Bayesian inference methods with a combined mitochondrial gene dataset, phylogenetic analyses generated identical, well-supported tree topologies for 71 Agaricomycetes species. This study reveals intraspecific and interspecific variations of the Ganoderma mitogenomes, which promotes the understanding of the origin, evolution, and genetic diversity of Ganoderma species.

Whole-Genome Sequencing and Analysis of the White-Rot Fungus Ceriporia lacerata Reveals Its Phylogenetic Status and the Genetic Basis of Lignocellulose Degradation and Terpenoid Synthesis

Ceriporia lacerata is an endophytic white-rot fungus that has lignocellulolytic and terpenoid-biosynthetic abilities. However, little is known about the genomic architecture of this fungus, even at the genus level. In this study, we present the first de novo genome assembly of C. lacerata (CGMCC No. 10485), based on PacBio long-read and Illumina short-read sequencing. The size of the C. lacerata genome is approximately 36 Mb (N50, 3.4 Mb). It encodes a total of 13,243 genes, with further functional analysis revealing that these genes are primarily involved in primary metabolism and host interactions in this strain's saprophytic lifestyle. Phylogenetic analysis based on ITS demonstrated a primary evolutionary position for C. lacerata, while the phylogenetic analysis based on orthogroup inference and average nucleotide identity revealed high-resolution phylogenetic details in which Ceriporia, Phlebia, Phlebiopsis, and Phanerochaete belong to the same evolutionary clade within the order Polyporales. Annotation of carbohydrate-active enzymes across the genome yielded a total of 806 genes encoding enzymes that decompose lignocellulose, particularly ligninolytic enzymes, lytic polysaccharides monooxygenases, and enzymes involved in the biodegradation of aromatic components. These findings illustrate the strain's adaptation to woody habitats, which requires the degradation of lignin and various polycyclic aromatic hydrocarbons. The terpenoid-production potential of C. lacerata was evaluated by comparing the genes of terpenoid biosynthetic pathways across nine Polyporales species. The shared genes highlight the major part of terpenoid synthesis pathways, especially the mevalonic acid pathway, as well as the main pathways of sesquiterpenoid, monoterpenoid, diterpenoid, and triterpenoid synthesis, while the strain-specific genes illustrate the distinct genetic factors determining the synthesis of structurally diverse terpenoids. This is the first genomic analysis of a species from this genus that we are aware of, and it will help advance functional genome research and resource development of this important fungus for applications in renewable energy, pharmaceuticals, and agriculture.

Phylogenomics and Comparative Genomics Highlight Specific Genetic Features in Ganoderma Species

The Ganoderma species in Polyporales are ecologically and economically relevant wood decayers used in traditional medicine, but their genomic traits are still poorly documented. In the present study, we carried out a phylogenomic and comparative genomic analyses to better understand the genetic blueprint of this fungal lineage. We investigated seven Ganoderma genomes, including three new genomes, G. australe, G. leucocontextum, and G. lingzhi. The size of the newly sequenced genomes ranged from 60.34 to 84.27 Mb and they encoded 15,007 to 20,460 genes. A total of 58 species, including 40 white-rot fungi, 11 brown-rot fungi, four ectomycorrhizal fungi, one endophyte fungus, and two pathogens in Basidiomycota, were used for phylogenomic analyses based on 143 single-copy genes. It confirmed that Ganoderma species belong to the core polyporoid clade. Comparing to the other selected species, the genomes of the Ganoderma species encoded a larger set of genes involved in terpene metabolism and coding for secreted proteins (CAZymes, lipases, proteases and SSPs). Of note, G. australe has the largest genome size with no obvious genome wide duplication, but showed transposable elements (TEs) expansion and the largest set of terpene gene clusters, suggesting a high ability to produce terpenoids for medicinal treatment. G. australe also encoded the largest set of proteins containing domains for cytochrome P450s, heterokaryon incompatibility and major facilitator families. Besides, the size of G. australe secretome is the largest, including CAZymes (AA9, GH18, A01A), proteases G01, and lipases GGGX, which may enhance the catabolism of cell wall carbohydrates, proteins, and fats during hosts colonization. The current genomic resource will be used to develop further biotechnology and medicinal applications, together with ecological studies of the Ganoderma species.

Characterization of γ-Cadinene Enzymes in Ganoderma lucidum and Ganoderma sinensis from Basidiomycetes Provides Insight into the Identification of Terpenoid Synthases

Enzymes boost protein engineering, directed evolution, and the biochemical industry and are also the cornerstone of metabolic engineering. Basidiomycetes are known to produce a large variety of terpenoids with unique structures. However, basidiomycetous terpene synthases remain largely untapped. Therefore, we provide a modeling method to obtain specific terpene synthases. Aided by bioinformatics analysis, three γ-cadinene enzymes from Ganoderma lucidum and Ganoderma sinensis were accurately predicted and identified experimentally. Based on the highly conserved amino motifs of the characterized γ-cadinene enzymes, the enzyme was reassembled as model 1. Using this model as a template, 67 homologous sequences of the γ-cadinene enzyme were screened from the National Center for Biotechnology Information (NCBI). According to the 67 sequences, the same gene structure, and similar conserved motifs to model 1, the γ-cadinene enzyme model was further improved by the same construction method and renamed as model 2. The results of bioinformatics analysis show that the conservative regions of models 1 and 2 are highly similar. In addition, five of these sequences were verified, 100% of which were γ-cadinene enzymes. The accuracy of the prediction ability of the γ-cadinene enzyme model was proven. In the same way, we also reanalyzed the identified Δ⁶-protoilludene enzymes in fungi and (-)-α-bisabolol enzymes in plants, all of which have their own unique conserved motifs. Our research method is expected to be used to study other terpenoid synthases with a similar or the same function in basidiomycetes, ascomycetes, bacteria, and plants and to provide rich enzyme resources.

DNA methylation across the tree of life, from micro to macro-organism

DNA methylation is a process in which methyl (CH3) groups are added to the DNA molecule. The DNA segment does not change in the sequence, but DNA methylation could alter the action of DNA. Different enzymes like DNA methyltransferases (DNMTs) take part in methylation of cytosine/adenine nucleosides in DNA. In prokaryotes, DNA methylation is performed to prevent the attack of phage and also plays a role in the chromosome replication and repair. In fungi, DNA methylation is studied to see the transcriptional changes, as in insects, the DNA methylation is not that well-known, it plays a different role like other organisms. In mammals, the DNA methylation is related to different types of cancers and plays the most important role in the placental development and abnormal DNA methylation connected with diseases like cancer, autoimmune diseases, and rheumatoid arthritis.

Feasible strategies for studying the involvement of DNA methylation and histone acetylation in the stress-induced formation of quality-related metabolites in tea (Camellia sinensis)

Tea plants are subjected to multiple stresses during growth, development, and postharvest processing, which affects levels of secondary metabolites in leaves and influences tea functional properties and quality. Most studies on secondary metabolism in tea have focused on gene, protein, and metabolite levels, whereas upstream regulatory mechanisms remain unclear. In this review, we exemplify DNA methylation and histone acetylation, summarize the important regulatory effects that epigenetic modifications have on plant secondary metabolism, and discuss feasible research strategies to elucidate the underlying specific epigenetic mechanisms of secondary metabolism regulation in tea. This information will help researchers investigate the epigenetic regulation of secondary metabolism in tea, providing key epigenetic data that can be used for future tea genetic breeding.

Genome of Ganoderma Species Provides Insights Into the Evolution, Conifers Substrate Utilization, and Terpene Synthesis for Ganoderma tsugae

Ganoderma tsugae is an endemic medicinal mushroom in Northeast China, providing important source of pharmaceutical product. Comparing with other Ganoderma species, wild G. tsugae can utilize coniferous wood. However, functional genes related to medicinal component synthesis and the genetic mechanism of conifer substrate utilization is still obscure. Here, we assembled a high-quality G. tsugae genome with 18 contigs and 98.5% BUSCO genes and performed the comparative genomics with other Ganoderma species. G. tsugae diverged from their common ancestor of G. lingzhi and G. sinense about 21 million years ago. Genes in G. tsugae-specific and G. tsugae-expanded gene families, such as salh, phea, cyp53a1, and cyp102a, and positively selected genes, such as glpk and amie, were functionally enriched in plant-pathogen interaction, benzoate degradation, and fanconi anemia pathway. Those functional genes might contribute to conifer substrate utilization of G. tsugae. Meanwhile, gene families in the terpene synthesis were identified and genome-wide SNP variants were detected in population. Finally, the study provided valuable genomic resources and offered useful hints for the functional gene mapping and investigation of key gene contributing to conifer cultivation substrate utilization and medicinal component biosynthesis.

Evolution and enrichment of CYP5035 in Polyporales: functionality of an understudied P450 family

Abstract

Bioprospecting for innovative basidiomycete cytochrome P450 enzymes (P450s) is highly desirable due to the fungi’s enormous enzymatic repertoire and outstanding ability to degrade lignin and detoxify various xenobiotics. While fungal metagenomics is progressing rapidly, the biocatalytic potential of the majority of these annotated P450 sequences usually remains concealed, although functional profiling identified several P450 families with versatile substrate scopes towards various natural products. Functional knowledge about the CYP5035 family, for example, is largely insufficient. In this study, the families of the putative P450 sequences of the four white-rot fungi Polyporus arcularius, Polyporus brumalis, Polyporus squamosus and Lentinus tigrinus were assigned, and the CYPomes revealed an unusual enrichment of CYP5035, CYP5136 and CYP5150. By computational analysis of the phylogeny of the former two P450 families, the evolution of their enrichment could be traced back to the Ganoderma macrofungus, indicating their evolutionary benefit. In order to address the knowledge gap on CYP5035 functionality, a representative subgroup of this P450 family of P. arcularius was expressed and screened against a test set of substrates. Thereby, the multifunctional enzyme CYP5035S7 converting several plant natural product classes was discovered. Aligning CYP5035S7 to 102,000 putative P450 sequences of 36 fungal species from Joint Genome Institute-provided genomes located hundreds of further CYP5035 family members, which subfamilies were classified if possible. Exemplified by these specific enzyme analyses, this study gives valuable hints for future bioprospecting of such xenobiotic-detoxifying P450s and for the identification of their biocatalytic potential.

Graphical abstract

Key points

• The P450 families CYP5035 and CYP5136 are unusually enriched in P. arcularius.

• Functional screening shows CYP5035 assisting in the fungal detoxification mechanism.

• Some Polyporales encompass an unusually large repertoire of detoxification P450s.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00253-021-11444-2.

Identification of non-ribosomal peptide synthetase in Ganoderma boninense Pat. that was expressed during the interaction with oil palm

Basal stem rot (BSR) of oil palm is a disastrous disease caused by a white-rot fungus Ganoderma boninense Pat. Non-ribosomal peptides (NRPs) synthesized by non-ribosomal peptide synthetases (NRPSs) are a group of secondary metabolites that act as fungal virulent factors during pathogenesis in the host. In this study, we aimed to isolate NRPS gene of G. boninense strain UPMGB001 and investigate the role of this gene during G. boninense-oil palm interaction. The isolated NRPS DNA fragment of 8322 bp was used to predict the putative peptide sequence of different domains and showed similarity with G. sinense (85%) at conserved motifs of three main NRPS domains. Phylogenetic analysis of NRPS peptide sequences demonstrated that NRPS of G. boninense belongs to the type VI siderophore family. The roots of 6-month-old oil palm seedlings were artificially inoculated for studying NRPS gene expression and disease severity in the greenhouse. The correlation between high disease severity (50%) and high expression (67-fold) of G. boninense NRPS gene at 4 months after inoculation and above indicated that this gene played a significant role in the advancement of BSR disease. Overall, these findings increase our knowledge on the gene structure of NRPS in G. boninense and its involvement in BSR pathogenesis as an effector gene.

Identifying a melanogenesis-related candidate gene by a high-quality genome assembly and population diversity analysis in Hypsizygus marmoreus

Hypsizygus marmoreus is one of the most important edible fungi in Basidiomycete division and includes white and gray strains. However, very limited knowledge is known about the genomic structures and the genetic basis for the white/gray diversity of this mushroom. Here, we report the near-complete high-quality H. marmoreus genome at the chromosomal level. Comparative genomics analysis indicates that chromosome structures were relatively conserved, and variations in collinearity and chromosome number were mainly attributed by chromosome split/fusion events in Aragicales, whereas the fungi genome experienced many genomic chromosome fracture, fusion, and genomic replication events after the split of Aragicales from Basidiomycetes. Resequencing of 57 strains allows us to classify the population into four major groups and associate genetic variations with morphological features, indicating that white strains were not originated independently. We further generated genetic populations and identified a cytochrome P450 as the candidate causal gene for the melanogenesis in H. marmoreus based on bulked segregant analysis (BSA) and comparative transcriptome analysis. The high-quality H. marmoreus genome and diversity data compiled in this study provide new knowledge and resources for the molecular breeding of H. marmoreus as well as the evolution of Basidiomycete.

Genome sequencing and annotation and phylogenomic analysis of the medicinal mushroom Amauroderma rugosum, a traditional medicinal species in the family Ganodermataceae

Amauroderma rugosum is one of the traditional Chinese medicinal mushrooms and is used to reduce inflammation, treat diuretic and upset stomach, and prevent cancer. Here, we present a genomic resource of Amauroderma rugosum (ACCC 51706) for further understanding its biology and exploration of the synthesis pathway of bioactive compounds. Genomic DNA was extracted and then subjected to Illumina HiSeq X Ten and PacBio Sequel I sequencing. The final genome is 40.66 Mb in size, with an N₅₀ scaffold size of 36.6 Mb, and encodes 10 181 putative predicted genes. Among them, 6931 genes were functionally annotated. Phylogenomic analysis suggested that A. rugosum and Ganoderma sinense were not clustered together into a group and the latter was grouped with the Polyporaceae. Further, we also identified 377 carbohydrate-active enzymes (CAZymes) and 15 secondary metabolite biosynthetic gene clusters. This is the first genome-scale assembly and annotation for an Amauroderma species. The identification of novel secondary metabolite biosynthetic gene clusters would promote pharmacological research and development of novel bioactive compounds in the future.

Diversity of Fungal DNA Methyltransferases and Their Association With DNA Methylation Patterns

DNA methyltransferases (DNMTs) are a group of proteins that catalyze DNA methylation by transferring a methyl group to DNA. The genetic variation in DNMTs results in differential DNA methylation patterns associated with various biological processes. In fungal species, DNMTs and their DNA methylation profiles were found to be very diverse and have gained many research interests. We reviewed fungal DNMTs in terms of their biological functions, protein domain structures, and their associated epigenetic regulations compared to those known in plant and animal systems. In addition, we summarized recent reports on potential RNA-directed DNA methylation (RdDM) related to DNMT5 in fungi. We surveyed up to 40 fungal species with published genome-wide DNA methylation profiles (methylomes) and presented the associations between the specific patterns of fungal DNA methylation and their DNMTs based on a phylogenetic tree of protein domain structures. For example, the main DNMTs in Basidiomycota, DNMT1 with RFD domain + DNMT5, contributing to CG methylation preference, were distinct from RID + Dim-2 in Ascomycota, resulting in a non-CG methylation preference. Lastly, we revealed that the dynamic methylation involved in fungal life stage changes was particularly low in mycelium and DNA methylation was preferentially located in transposable elements (TEs). This review comprehensively discussed fungal DNMTs and methylomes and their connection with fungal development and taxonomy to present the diverse usages of DNA methylation in fungal genomes.

DNA Methylation Is Responsive to the Environment and Regulates the Expression of Biosynthetic Gene Clusters, Metabolite Production, and Virulence in Fusarium graminearum

Histone modifications play a significant role in the regulation of biosynthetic gene clusters (BGCs) in the phytopathogen Fusarium graminearum, by contrast, epigenetic regulation by DNA methyltransferases (DNMTs) is less documented. In this study, we characterized two DNMTs (FgDIM-2 and FgRID) in F. graminearum, with homologies to "Deficient in methylation" (DIM-2) and "Repeat-induced point (RIP) deficient" (RID) from Neurospora. The loss of DNMTs resulted in not only a decrease in average methylation density in the nutrient-poor, compared to nutrient-rich conditions, but also differences in the genes expressed between the WT and the DNMT mutant strains, implicating the external environment as an important trigger in altering DNA methylation patterns. Consequently, we observed significant changes in the regulation of multiple BGCs and alterations in the pathogenicity of the fungus.

Immunomodulatory Effects of Edible and Medicinal Mushrooms and Their Bioactive Immunoregulatory Products

Mushrooms have been valued as food and health supplements by humans for centuries. They are rich in dietary fiber, essential amino acids, minerals, and many bioactive compounds, especially those related to human immune system functions. Mushrooms contain diverse immunoregulatory compounds such as terpenes and terpenoids, lectins, fungal immunomodulatory proteins (FIPs) and polysaccharides. The distributions of these compounds differ among mushroom species and their potent immune modulation activities vary depending on their core structures and fraction composition chemical modifications. Here we review the current status of clinical studies on immunomodulatory activities of mushrooms and mushroom products. The potential mechanisms for their activities both in vitro and in vivo were summarized. We describe the approaches that have been used in the development and application of bioactive compounds extracted from mushrooms. These developments have led to the commercialization of a large number of mushroom products. Finally, we discuss the problems in pharmacological applications of mushrooms and mushroom products and highlight a few areas that should be improved before immunomodulatory compounds from mushrooms can be widely used as therapeutic agents.

Transcriptome based functional identification and application of regulator AbrB on alkaline protease synthesis in Bacillus licheniformis 2709

Bacillus licheniformis 2709 is the major alkaline protease producer, which has great potential value of industrial application, but how the high-producer can be regulated rationally is still not completely understood. It's meaningful to understand the metabolic processes during alkaline protease production in industrial fermentation medium. Here, we collected the transcription database at various enzyme-producing stages (preliminary stage, stable phase and decline phase) to specifically research the synthesized and regulatory mechanism of alkaline protease in B. licheniformis. The RNA-sequencing analysis showed differential expression of numerous genes related to several processes, among which genes correlated with regulators were concerned, especially the major differential gene abrB on enzyme (AprE) synthesis was investigated. It was further verified that AbrB is a repressor of AprE by plasmid-mediated over-expression due to the severely descending enzyme activity (11,300 U/mL to 2695 U/mL), but interestingly it is indispensable for alkaline protease production because the enzyme activity of the null abrB mutant was just about 2279 U/mL. Thus, we investigated the aprE transcription by eliminating the theoretical binding site (TGGAA) of AbrB protein predicated by computational strategy, which significantly improved the enzyme activity by 1.21-fold and gene transcription level by 1.77-fold in the mid-log phase at a cultivation time of 18 h. Taken together, it is of great significance to improve the production strategy, control the metabolic process and oriented engineering by rational molecular modification of regulatory network based on the high throughput sequencing and computational prediction.

Genomic data mining approaches for the discovery of anticancer peptides from Ganoderma sinense

Drug discovery from traditional Chinese medicine (TCM) typically involves the extraction of active ingredients from natural products with high biological activity and function from the vast repertoire of traditional Chinese medicine. This strategy cannot fully exploit the vast resources of TCM. Known as the longevity mushroom, Ganoderma spp. has been used as medicine for thousands of years. Recent studies have demonstrated its anticancer activity. While most research on Ganoderma spp. has focused on their polysaccharides or small molecules as potential anticancer components, possible anticancer peptides (ACPs) or proteins have been neglected. In this study, genomic data mining approaches were used to discover potential ACPs from Ganoderma sinense. A search against known ACPs identified 477 proteins in the G. sinense proteome that possess putative ACP sequences and that thus may serve as parent proteins. After in silico digestion by trypsin, 34 G. sinense proteins were predicted to release putative ACPs (by the mACPpred program). A subsequent sequence similarity comparison against known ACPs identified 15 trypsin-digested fragments as possible ACPs, of which 3 sequences were identical to known ACPs. The results indicated that ACPs may be involved in the anticancer activity of G. sinense and that genomic mining approaches can be effective strategies for discovering active components in TCM resources. The accumulation of genomic and proteomic data will undoubtedly accelerate drug discovery from TCM resources.

Genome-Wide Characterization and Comparative Analysis of MYB Transcription Factors in Ganoderma Species

Numerous studies in plants have shown the vital roles of MYB transcription factors in signal transduction, developmental regulation, biotic/abiotic stress responses and secondary metabolism regulation. However, less is known about the functions of MYBs in Ganoderma. In this study, five medicinal macrofungi of genus Ganoderma were subjected to a genome-wide comparative analysis of MYB genes. A total of 75 MYB genes were identified and classified into four types: 1R-MYBs (52), 2R-MYBs (19), 3R-MYBs (2) and 4R-MYBs (2). Gene structure analysis revealed varying exon numbers (3-14) and intron lengths (7-1058 bp), and noncanonical GC-AG introns were detected in G. lucidum and G. sinense. In a phylogenetic analysis, 69 out of 75 MYB genes were clustered into 15 subgroups, and both single-copy orthologous genes and duplicated genes were identified. The promoters of the MYB genes harbored multiple cis-elements, and specific genes were co-expressed with the G. lucidum MYB genes, indicating the potential roles of these MYB genes in stress response, development and metabolism. This comprehensive and systematic study of MYB family members provides a reference and solid foundation for further functional analysis of MYB genes in Ganoderma species.

The Pattern and Function of DNA Methylation in Fungal Plant Pathogens

To successfully infect plants and trigger disease, fungal plant pathogens use various strategies that are dependent on characteristics of their biology and genomes. Although pathogenic fungi are different from animals and plants in the genomic heritability, sequence feature, and epigenetic modification, an increasing number of phytopathogenic fungi have been demonstrated to share DNA methyltransferases (MTases) responsible for DNA methylation with animals and plants. Fungal plant pathogens predominantly possess four types of DNA MTase homologs, including DIM-2, DNMT1, DNMT5, and RID. Numerous studies have indicated that DNA methylation in phytopathogenic fungi mainly distributes in transposable elements (TEs), gene promoter regions, and the repetitive DNA sequences. As an important and heritable epigenetic modification, DNA methylation is associated with silencing of gene expression and transposon, and it is responsible for a wide range of biological phenomena in fungi. This review highlights the relevant reports and insights into the important roles of DNA methylation in the modulation of development, pathogenicity, and secondary metabolism of fungal plant pathogens. Recent evidences prove that there are massive links between DNA and histone methylation in fungi, and they commonly regulate fungal development and mycotoxin biosynthesis.

Evaluation of DNA Methylation Changes by CRED-RA Analysis Following Prednisone Treatment of Endophyte, Fusarium oxysporum

Endophytes that represent a sub-set of plant resident microbes are a reservoir of bioactive metabolites. Many of the secondary metabolite biosynthetic gene clusters of endophytes are silent under axenic culture conditions. Epigenetic reprogramming of such cryptic pathways is possible by use of small molecule modulators like prednisone. Methylation changes induced by prednisone, a hypomethylating epigenetic modulator were studied in endophytic Fusarium oxysporum. CRED-RA analysis following exposure to non-cytotoxic dose (300 µM) revealed prednisone as effective in inducing non-methylation and semi-methylation pattern while inhibiting full-methylation of the genome. Effectiveness of prednisone as a DNA methyl transferase inhibitor can be explored in future to study alterations in secondary metabolite gene expression profile in endophytic F. oxysporum.

Draft Genome Sequence of a Phytopathogenic Ganoderma sp. Strain That Causes Basal Stem Rot Disease on Oil Palm in Sabah, Malaysia

Basal stem rot (BSR) disease on Elaeis guineens is known to be caused by members of the pathogenic fungal genus Ganoderma, especially the species Ganoderma boninense. This species affects oil palm plantation in Sabah, Malaysia. The genome sequence (52.28 Mbp) will add to the representation of this genus, especially in regard to BSR disease.

Basal stem rot (BSR) disease on Elaeis guineens is known to be caused by members of the pathogenic fungal genus Ganoderma, especially the species Ganoderma boninense. This species affects oil palm plantation in Sabah, Malaysia. The genome sequence (52.28 Mbp) will add to the representation of this genus, especially in regard to BSR disease.

Genome-wide DNA methylation analysis of Metarhizium anisopliae during tick mimicked infection condition

BACKGROUND

The Metarhizium genus harbors important entomopathogenic fungi. These species have been widely explored as biological control agents, and strategies to improve the fungal virulence are under investigation. Thus, the interaction between Metarhizium species and susceptible hosts have been explored employing different methods in order to characterize putative virulence determinants. However, the impact of epigenetic modulation on the infection cycle of Metarhizium is still an open topic. Among the different epigenetic modifications, DNA methylation of cytosine bases is an important mechanism to control gene expression in several organisms. To better understand if DNA methylation can govern Metarhizium-host interactions, the genome-wide DNA methylation profile of Metarhizium anisopliae was explored in two conditions: tick mimicked infection and a saprophytic-like control.

RESULTS

Using a genome wide DNA methylation profile based on bisulfite sequencing (BS-Seq), approximately 0.60% of the total cytosines were methylated in saprophytic-like condition, which was lower than the DNA methylation level (0.89%) in tick mimicked infection condition. A total of 670 mRNA genes were found to be putatively methylated, with 390 mRNA genes uniquely methylated in the tick mimicked infection condition. GO terms linked to response to stimuli, cell wall morphogenesis, cytoskeleton morphogenesis and secondary metabolism biosynthesis were over-represented in the tick mimicked infection condition, suggesting that energy metabolism is directed towards the regulation of genes associated with infection. However, recognized virulence determinants known to be expressed at distinct infection steps, such as the destruxin backbone gene and the collagen-like protein gene Mcl1, were found methylated, suggesting that a dynamic pattern of methylation could be found during the infectious process. These results were further endorsed employing RT-qPCR from cultures treated or not with the DNA methyltransferase inhibitor 5-Azacytidine.

CONCLUSIONS

The set of genes here analyzed focused on secondary metabolites associated genes, known to be involved in several processes, including virulence. The BS-Seq pipeline and RT-qPCR analysis employing 5-Azacytidine led to identification of methylated virulence genes in M. anisopliae. The results provided evidences that DNA methylation in M. anisopliae comprises another layer of gene expression regulation, suggesting a main role of DNA methylation regulating putative virulence determinants during M. anisopliae infection cycle.

Genomic Analyses Provide Insights Into the Evolutionary History and Genetic Diversity of Auricularia Species

Species in the genus Auricularia play important roles for people’s food and nutrition especially Auricularia cornea and A. heimuer. To understand their evolutionary history, genome structure, and population-level genetic variation, we performed a high-quality genome sequencing of Auricularia cornea and the corresponding comparative genomic analysis. The genome size of A. cornea was similar to Auricularia subglabra, but 1.5 times larger than that of A. heimuer. Several factors were responsible for genome size variation including gene numbers, repetitive elements, and gene lengths. Phylogenomic analysis revealed that the estimated divergence time between A. heimuer and other Auricularia is ∼79.1 million years ago (Mya), while the divergence between A. cornea and A. subglabra occurred in ∼54.8 Mya. Population genomic analysis also provided insight into the demographic history of A. cornea and A. heimuer, indicating that their populations fluctuated over time with global climate change during Marine Isotope Stage 5-2. Moreover, despite the highly similar external morphologies of A. cornea and A. heimuer, their genomic properties were remarkably different. The A. cornea genome only shared 14% homologous syntenic blocks with A. heimuer and possessed more genes encoding carbohydrate-active enzymes and secondary metabolite biosynthesis proteins. The cross-taxa transferability rates of simple sequence repeat (SSR) and insertion or deletion (InDel) markers within the genus Auricularia were also lower than that previously observed for species within the same genus. Taken together, these results indicate a high level of genetic differentiation between these two Auricularia species. Consequently, our study provides new insights into the genomic evolution and genetic differentiation of Auricularia species that will facilitate future genetic breeding.

GSP-2, a polysaccharide extracted from Ganoderma sinense, is a novel toll-like receptor 4 agonist

Ganoderma sinense is a Chinese unique medicinal fungus that has been used in folk medicine for thousands of years. Polysaccharides are considered to be biologically active ingredients due to their immune-modulating functions. Previously we found that GSP-2, a new polysaccharide isolated from Ganoderma sinense, exerts an immunomodulatory effect in human peripheral blood mononuclear cells but the underlying mechanism is unclear. The present study aimed to investigate how GSP-2 triggers immunologic responses and the implicated signaling pathways. GSP-2 effects were investigated both in a macrophagic cell line, RAW264.7, and in primary macrophages. Moreover, the molecular basis of GSP-2 recognition by immune cells, and the consequent activation of signaling cascades, were explored by employing recombinant human HEK293-TLR-Blue clones, individually overexpressing various Toll-like receptors. GSP-2 dose-dependently induced the overexpression of Toll-like receptor 4 (TLR4) but did not affect the expression of other TLRs. Moreover, GSP-2 induced TNFα secretion in primary macrophages from wild-type, but not TLR4-knockout mice. In addition, GSP-2 upregulated TLR4 protein expression and activated the MAPK pathway in RAW246.7 macrophages. Finally, GSP-2 induced the production of the cytokines TNFα, IL1β, and IL6. Our data demonstrated that GSP-2 was specifically recognized by TLR4, promoting cytokine secretion and immune modulation in macrophages.

Transcriptomic Insights into Benzenamine Effects on the Development, Aflatoxin Biosynthesis, and Virulence of Aspergillus flavus

Aspergillus flavus is a soilborne pathogenic fungus that poses a serious public health threat due to it contamination of food with carcinogenic aflatoxins. Our previous studies have demonstrated that benzenamine displayed strong inhibitory effects on the mycelial growth of A. flavus. In this study, we systematically investigated the inhibitory effects of benzenamine on the development, aflatoxin biosynthesis, and virulence in A. flavus, as well as the underlying mechanism. The results indicated that benzenamine exhibited great capacity to combat A. flavus at a concentration of 100 µL/L, leading to significantly decreased aflatoxin accumulation and colonization capacity in maize. The transcriptional profile revealed that 3589 genes show altered mRNA levels in the A. flavus after treatment with benzenamine, including 1890 down-regulated and 1699 up-regulated genes. Most of the differentially expressed genes participated in the biosynthesis and metabolism of amino acid, purine metabolism, and protein processing in endoplasmic reticulum. Additionally, the results brought us to a suggestion that benzenamine affects the development, aflatoxin biosynthesis, and pathogenicity of A. flavus via down-regulating related genes by depressing the expression of the global regulatory factor leaA. Overall, this study indicates that benzenamine have tremendous potential to act as a fumigant against pathogenic A. flavus. Furthermore, this work offers valuable information regarding the underlying antifungal mechanism of benzenamine against A. flavus at the level of transcription, and these potential targets may be conducive in developing new strategies for preventing aflatoxin contamination.

Genome-wide DNA methylation and transcriptomic profiles in the lifestyle strategies and asexual development of the forest fungal pathogen Heterobasidion parviporum

Heterobasidion parviporum is the most devastating fungal pathogen of conifer forests in Northern Europe. The fungus has dual life strategies, necrotrophy on living trees and saprotrophy on dead woods. DNA cytosine methylation is an important epigenetic modification in eukaryotic organisms. Our presumption is that the lifestyle transition and asexual development in H. parviporum could be driven by epigenetic effects. Involvements of DNA methylation in the regulation of aforementioned processes have never been studied thus far. RNA-seq identified lists of highly induced genes enriched in carbohydrate-active enzymes during necrotrophic interaction with host trees and saprotrophic sawdust growth. It also highlighted signaling- and transcription factor-related genes potentially associated with the transition of saprotrophic to necrotrophic lifestyle and groups of primary cellular activities throughout asexual development. Whole-genome bisulfite sequencing revealed that DNA methylation displayed pronounced preference in CpG dinucleotide context across the genome and mostly targeted transposable element (TE)-rich regions. TE methylation level demonstrated a strong negative correlation with TE expression, reinforcing the protective function of DNA methylation in fungal genome stability. Small groups of genes putatively subject to methylation transcriptional regulation in response to saprotrophic and necrotrophic growth in comparison with free-living mycelia were also explored. Our study reported on the first methylome map of a forest pathogen. Analysis of transcriptome and methylome variations associated with asexual development and different lifestyle strategies provided further understanding of basic biological processes in H. parviporum. More importantly, our work raised additional potential roles of DNA methylation in fungi apart from controlling the proliferation of TEs.

Trends in herbgenomics

From Shen Nong's Herbal Classic (Shennong Bencao Jing) to the Compendium of Materia Medica (Bencao Gangmu) and the first scientific Nobel Prize for the mainland of China, each milestone in the historical process of the development of traditional Chinese medicine (TCM) involves screening, testing and integrating. After thousands of years of inheritance and development, herbgenomics (bencaogenomics) has bridged the gap between TCM and international advanced omics studies, promoting the application of frontier technologies in TCM. It is a discipline that uncovers the genetic information and regulatory networks of herbs to clarify their molecular mechanism in the prevention and treatment of human diseases. The main theoretical system includes genomics, functional genomics, proteomics, transcriptomics, metabolomics, epigenomics, metagenomics, synthetic biology, pharmacogenomics of TCM, and bioinformatics, among other fields. Herbgenomics is mainly applicable to the study of medicinal model plants, genomic-assisted breeding, herbal synthetic biology, protection and utilization of gene resources, TCM quality evaluation and control, and TCM drug development. Such studies will accelerate the application of cutting-edge technologies, revitalize herbal research, and strongly promote the development and modernization of TCM.

Whole-Genome Resequencing and Pan-Transcriptome Reconstruction Highlight the Impact of Genomic Structural Variation on Secondary Metabolite Gene Clusters in the Grapevine Esca Pathogen Phaeoacremonium minimum

The Ascomycete fungus Phaeoacremonium minimum is one of the primary causal agents of Esca, a widespread and damaging grapevine trunk disease. Variation in virulence among Pm. minimum isolates has been reported, but the underlying genetic basis of the phenotypic variability remains unknown. The goal of this study was to characterize intraspecific genetic diversity and explore its potential impact on virulence functions associated with secondary metabolism, cellular transport, and cell wall decomposition. We generated a chromosome-scale genome assembly, using single molecule real-time sequencing, and resequenced the genomes and transcriptomes of multiple isolates to identify sequence and structural polymorphisms. Numerous insertion and deletion events were found for a total of about 1 Mbp in each isolate. Structural variation in this extremely gene dense genome frequently caused presence/absence polymorphisms of multiple adjacent genes, mostly belonging to biosynthetic clusters associated with secondary metabolism. Because of the observed intraspecific diversity in gene content due to structural variation we concluded that a transcriptome reference developed from a single isolate is insufficient to represent the virulence factor repertoire of the species. We therefore compiled a pan-transcriptome reference of Pm. minimum comprising a non-redundant set of 15,245 protein-coding sequences. Using naturally infected field samples expressing Esca symptoms, we demonstrated that mapping of meta-transcriptomics data on a multi-species reference that included the Pm. minimum pan-transcriptome allows the profiling of an expanded set of virulence factors, including variable genes associated with secondary metabolism and cellular transport.

Co-extraction of genomic DNA & total RNA from recalcitrant woody tissues for next-generation sequencing studies

The successful implementation of next-generation sequencing techniques in plant and woody tissues depends on the quality of initial starting material. This study demonstrated the use of a modified protocol that enabled the simultaneous extraction of both genomic DNA and total RNA from recalcitrant woody material. The genetic material obtained by this protocol is of high quality and can be directly used in downstream analysis (e.g., next-generation sequencing). This protocol is particularly useful not only when the initial plant material is limited but also when genomic DNA features (e.g., methylation) have to be compared with the total RNA (e.g., gene expression). For such studies, the extraction from the same materials is highly preferred to minimize sample variation.

The advancement of next-generation sequencing techniques has greatly facilitated almost all branches of biological and life-science studies. To guarantee reliable results and conclusions, high-quality starting genetic material (genomic DNA and total RNA) is essential. Our modified protocol, which is based on existing methods, enabled the efficient simultaneous extraction of genomic DNA and total RNA from the same recalcitrant woody sample. Thereby, it greatly minimizes the amount of initial plant materials needed for a specific study.

Comprehensive Characterization for Ginsenosides Biosynthesis in Ginseng Root by Integration Analysis of Chemical and Transcriptome

Herbgenomics provides a global platform to explore the genetics and biology of herbs on the genome level. Panax ginseng C.A. Meyer is an important medicinal plant with numerous pharmaceutical effects. Previous reports mainly discussed the transcriptome of ginseng at the organ level. However, based on mass spectrometry imaging analyses, the ginsenosides varied among different tissues. In this work, ginseng root was separated into three tissues-periderm, cortex and stele-each for five duplicates. The chemical analysis and transcriptome analysis were conducted simultaneously. Gene-encoding enzymes involved in ginsenosides biosynthesis and modification were studied based on gene and molecule data. Eight widely-used ginsenosides were distributed unevenly in ginseng roots. A total of 182,881 unigenes were assembled with an N50 contig size of 1374 bp. About 21,000 of these unigenes were positively correlated with the content of ginsenosides. Additionally, we identified 192 transcripts encoding enzymes involved in two triterpenoid biosynthesis pathways and 290 transcripts encoding UDP-glycosyltransferases (UGTs). Of these UGTs, 195 UGTs (67.2%) were more highly expressed in the periderm, and that seven UGTs and one UGT were specifically expressed in the periderm and stele, respectively. This genetic resource will help to improve the interpretation on complex mechanisms of ginsenosides biosynthesis, accumulation, and transportation.

Transcriptome analyses provide insights into the difference of alkaloids biosynthesis in the Chinese goldthread (Coptis chinensis Franch.) from different biotopes

Coptis chinensis Franch., the Chinese goldthread (‘Weilian’ in Chinese), one of the most important medicinal plants from the family Ranunculaceae, and its rhizome has been widely used in Traditional Chinese Medicine for centuries. Here, we analyzed the chemical components and the transcriptome of the Chinese goldthread from three biotopes, including Zhenping, Zunyi and Shizhu. We built comprehensive, high-quality de novo transcriptome assemblies of the Chinese goldthread from short-read RNA-Sequencing data, obtaining 155,710 transcripts and 56,071 unigenes. More than 98.39% and 95.97% of core eukaryotic genes were found in the transcripts and unigenes respectively, indicating that this unigene set capture the majority of the coding genes. A total of 520,462, 493,718, and 507,247 heterozygous SNPs were identified in the three accessions from Zhenping, Zunyi, and Shizhu respectively, indicating high polymorphism in coding regions of the Chinese goldthread (∼1%). Chemical analyses of the rhizome identified six major components, including berberine, palmatine, coptisine, epiberberine, columbamine, and jatrorrhizine. Berberine has the highest concentrations, followed by coptisine, palmatine, and epiberberine sequentially for all the three accessions. The drug quality of the accession from Shizhu may be the highest among these accessions. Differential analyses of the transcriptome identified four pivotal candidate enzymes, including aspartate aminotransferaseprotein, polyphenol oxidase, primary-amine oxidase, and tyrosine decarboxylase, were significantly differentially expressed and may be responsible for the difference of alkaloids contents in the accessions from different biotopes.

Genomic CG dinucleotide deficiencies associated with transposable element hypermutation in Basidiomycetes, some lower fungi, a moss and a clubmoss

Many Basidiomycete genomes include substantial fractions that are deficient in CG dinucleotides, in extreme cases amounting to 70% of the genome. CG deficiency is variable and correlates with genome size and, more closely, with transposable element (TE) content. Many species have limited CG deficiency; it is therefore likely that there are other mechanisms that can control TE proliferation. Examination of TEs confirms that C-to-T transition mutations in CG dinucleotides may comprise a conspicuous proportion of differences between paired elements, however transition/transversion ratios are never as high as those due to RIP in some Ascomycetes, suggesting that repeat-associated CG mutation is not totally pervasive. This has allowed gene family expansion in Basidiomycetes, although CG transition differences are often prominent in paired gene family members, and are evidently responsible for destruction of some copies. A few lower fungal genomes exhibit similar evidence of repeat-associated CG mutation, as do the genomes of the two lower plants Physcomitrella patens and Selaginella moellendorffii, in both of which mutation parallels published methylation of CHG as well as CG nucleotides. In Basidiomycete DNA methylation has been reported to be largely confined to CG dinucleotides in repetitive DNA, but while methylation and mutation are evidently associated, it is not clear which is cause and which effect.

Biologically Active Secondary Metabolites from the Fungi

Many Fungi have a well-developed secondary metabolism. The diversity of fungal species and the diversification of biosynthetic gene clusters underscores a nearly limitless potential for metabolic variation and an untapped resource for drug discovery and synthetic biology. Much of the ecological success of the filamentous fungi in colonizing the planet is owed to their ability to deploy their secondary metabolites in concert with their penetrative and absorptive mode of life. Fungal secondary metabolites exhibit biological activities that have been developed into life-saving medicines and agrochemicals. Toxic metabolites, known as mycotoxins, contaminate human and livestock food and indoor environments. Secondary metabolites are determinants of fungal diseases of humans, animals, and plants. Secondary metabolites exhibit a staggering variation in chemical structures and biological activities, yet their biosynthetic pathways share a number of key characteristics. The genes encoding cooperative steps of a biosynthetic pathway tend to be located contiguously on the chromosome in coregulated gene clusters. Advances in genome sequencing, computational tools, and analytical chemistry are enabling the rapid connection of gene clusters with their metabolic products. At least three fungal drug precursors, penicillin K and V, mycophenolic acid, and pleuromutilin, have been produced by synthetic reconstruction and expression of respective gene clusters in heterologous hosts. This review summarizes general aspects of fungal secondary metabolism and recent developments in our understanding of how and why fungi make secondary metabolites, how these molecules are produced, and how their biosynthetic genes are distributed across the Fungi. The breadth of fungal secondary metabolite diversity is highlighted by recent information on the biosynthesis of important fungus-derived metabolites that have contributed to human health and agriculture and that have negatively impacted crops, food distribution, and human environments.

Survey of the genome of Pogostemon cablin provides insights into its evolutionary history and sesquiterpenoid biosynthesis

Pogostemon cablin (Blanco) Benth. (Patchouli) is an important traditional Chinese medicinal plant that has both essential oil value and a broad range of therapeutic effects. Here we report the first de novo assembled 1.15-Gb draft genome sequence for P. cablin from next-generation sequencing technology. Our assembly, with a misassembly rate of <4 bp per 100 kb, is ~73% of the predicted genome size (1.57 Gb). Analysis of whole-genome sequences identified 3,147,333 heterozygous single-nucleotide polymorphisms and 490,407 insertions and deletions, giving an estimated heterozygosity rate of 0.274%. A comprehensive annotation pipeline indicated that repetitive sequences make up 58.55% of the assemblies, and that there are estimated 45,020 genes. Comparative genomics analysis showed that the Phrymaceae and Lamiaceae family split ~62.80 Mya, and the divergence between patchouli and sesame occurred ~52.42 Mya, implying a potentially shared recent whole-genome duplication event. Analysis of gene homologs involved in sesquiterpenoid biosynthesis showed that patchouli contains key genes involved in more sesquiterpenoid types and has more copies of genes for each sesquiterpenoid type than several other related plant species. The patchouli genome will facilitate future research on secondary metabolic pathways and their regulation as well as potential selective breeding of patchouli.