Nudix hydrolase WvNUDX24 is involved in borneol biosynthesis in Wurfbainia villosa.

Borneol, camphor, and bornyl acetate are highly promising monoterpenoids widely used in medicine, flavor, food, and chemical applications. Bornyl diphosphate (BPP) serves as a common precursor for the biosynthesis of these monoterpenoids. Although bornyl diphosphate synthase (BPPS) that catalyzes the cyclization of geranyl diphosphate (GPP) to BPP has been identified in multiple plants, the enzyme responsible for the hydrolysis of BPP to produce borneol has not been reported. Here, we conducted in vitro and in vivo functional characterization to identify the Nudix hydrolase WvNUDX24 from W. villosa, which specifically catalyzes the hydrolysis of BPP to generate bornyl phosphate (BP), and then BP forms borneol under the action of phosphatase. Subcellular localization experiments indicated that the hydrolysis of BPP likely occurs in the cytoplasm. Furthermore, site-directed mutagenesis experiments revealed that four critical residues (R84, S96, P98, and G99) for the hydrolysis activity of WvNUDX24. Additionally, the functional identification of phosphatidic acid phosphatase (PAP) demonstrated that WvPAP5 and WvPAP10 were able to hydrolyze geranylgeranyl diphosphate (GGPP) and farnesyl diphosphate (FPP) to generate geranylgeranyl phosphate (GGP) and farnesyl phosphate (FP), respectively, but could not hydrolyze BPP, GPP, and neryl diphosphate (NPP) to produce corresponding monophosphate products. These findings highlight the essential role of WvNUDX24 in the first step of BPP hydrolysis to produce borneol and provide genetic elements for the production of BPP-related terpenoids through plant metabolic engineering and synthetic biology.

Comparing genomes of Fructus Amomi-producing species reveals genetic basis of volatile terpenoid divergence.

Wurfbainia longiligularis and Wurfbainia villosa are both rich in volatile terpenoids and are 2 primary plant sources of Fructus Amomi used for curing gastrointestinal diseases. Metabolomic profiling has demonstrated that bornyl diphosphate (BPP)-related terpenoids are more abundant in the W. villosa seeds and have a wider tissue distribution in W. longiligularis. To explore the genetic mechanisms underlying the volatile terpenoid divergence, a high-quality chromosome-level genome of W. longiligularis (2.29 Gb, contig N50 of 80.39 Mb) was assembled. Functional characterization of 17 terpene synthases (WlTPSs) revealed that WlBPPS, along with WlTPS 24/26/28 with bornyl diphosphate synthase (BPPS) activity, contributes to the wider tissue distribution of BPP-related terpenoids in W. longiligularis compared to W. villosa. Furthermore, transgenic Nicotiana tabacum showed that the GCN4-motif element positively regulates seed expression of WvBPPS and thus promotes the enrichment of BPP-related terpenoids in W. villosa seeds. Systematic identification and analysis of candidate TPS in 29 monocot plants from 16 families indicated that substantial expansion of TPS-a and TPS-b subfamily genes in Zingiberaceae may have driven increased diversity and production of volatile terpenoids. Evolutionary analysis and functional identification of BPPS genes showed that BPP-related terpenoids may be distributed only in the Zingiberaceae of monocot plants. This research provides valuable genomic resources for breeding and improving Fructus Amomi with medicinal and edible value and sheds light on the evolution of terpenoid biosynthesis in Zingiberaceae.

Bioactive α-Pyrone Analogs from the Endophytic Fungus Diaporthe sp. CB10100: α-Glucosidase Inhibitory Activity, Molecular Docking, and Molecular Dynamics Studies.

Two α-pyrone analogs were isolated from the endophytic fungus Diaporthe sp. CB10100, which is derived from the medicinal plant Sinomenium acutum. These analogs included a new compound, diaporpyrone F (3), and a known compound, diaporpyrone D (4). The structure of 3 was identified by a comprehensive examination of HRESIMS, 1D and 2D NMR spectroscopic data. Bioinformatics analysis revealed that biosynthetic gene clusters for α-pyrone analogs are common in fungi of Diaporthe species. The in vitro α-glucosidase inhibitory activity and antibacterial assay of 4 revealed that it has a 46.40% inhibitory effect on α-glucosidase at 800 µM, while no antibacterial activity against methicillin-resistant Staphylococcus aureus (MRSA), Mycolicibacterium (Mycobacterium) smegmatis or Klebsiella pneumoniae at 64 µg/mL. Molecular docking and molecular dynamics simulations of 4 with α-glucosidase further suggested that the compounds are potential α-glucosidase inhibitors. Therefore, α-pyrone analogs can be used as lead compounds for α-glucosidase inhibitors in more in-depth studies.

The Identification of SQS/SQE/OSC Gene Families in Regulating the Biosynthesis of Triterpenes in Potentilla anserina.

The tuberous roots of Potentilla anserina (Pan) are an edible and medicinal resource in Qinghai-Tibetan Plateau, China. The triterpenoids from tuberous roots have shown promising anti-cancer, hepatoprotective, and anti-inflammatory properties. In this study, we carried out phylogenetic analysis of squalene synthases (SQSs), squalene epoxidases (SQEs), and oxidosqualene cyclases (OSCs) in the pathway of triterpenes. In total, 6, 26, and 20 genes of SQSs, SQEs, and OSCs were retrieved from the genome of Pan, respectively. Moreover, 6 SQSs and 25 SQEs genes expressed in two sub-genomes (A and B) of Pan. SQSs were not expanded after whole-genome duplication (WGD), and the duplicated genes were detected in SQEs. Twenty OSCs were divided into two clades of cycloartenol synthases (CASs) and ß-amyrin synthases (ß-ASs) by a phylogenetic tree, characterized with gene duplication and evolutionary divergence. We speculated that ß-ASs and CASs may participate in triterpenes synthesis. The data presented act as valuable references for future studies on the triterpene synthetic pathway of Pan.

Degradation of mirubactin to multiple siderophores with varying Fe(III) chelation properties.

Three siderophores mirubactins B-D (4-6) were identified as the degradation products of previously isolated mirubactin (1). Their structures were revealed by HR-ESI-MS/MS, NMR analyses, and density functional calculations, among which 4 contains an unusual cyclic amidine functionality. Cyclic voltammetry showed that 5 and 6 have reduced iron complexing capacity. Mirubactin (1) and Fe(III) could also form a stable complex, which may be an ingenious approach to compete for iron acquisition by the producing organisms.

Isolation and Biosynthesis of Phenazine-Polyketide Hybrids from Streptomyces sp. KIB-H483.

A phenazine-polyketide hybrid compound, nexphenazine A (1), was isolated from Streptomyces sp. KIB-H483. The bioinformatic analysis of the draft genome of the producing strain and gene inactivation experiments revealed that the biosynthesis of 1 involves a phenazine-polyketide hybrid gene cluster. The abolished production of 1 as well as the accumulation of shunt metabolites 4-7 in mutant strain ΔnpzI revealed the key role of the npzI gene, which encodes an NAD(P)H-dependent ketoreductase, in nexphenazine biosynthesis. The structures and absolute configurations of the isolated intermediates were established on the basis of spectroscopic data analysis, single-crystal X-ray diffraction, chiral chromatography, and chemical conversion experiments. NpzI exhibited stereochemical selectivity in reducing the carbonyl group of 4. Nexphenazine biosynthesis is proposed to involve a condensation of the carboxyl group of phenazine with one molecule of methylmalonyl-CoA by a type I PKS, followed by a ketone reduction by NpzI and an unknown methylation reaction.

Identification of potential molecular mechanism related to craniofacial dysmorphism caused by FOXI3 deficiency.

BACKGROUND: Hemifacial macrosomia (HFM, OMIM 164210) is a complex and highly heterogeneous disease. FORKHEAD BOX I3 (FOXI3) is a susceptibility gene for HFM, and mice with loss of function of Foxi3 did exhibit a phenotype similar to craniofacial dysmorphism. However, the specific pathogenesis of HFM caused by FOXI3 deficiency remains unclear till now. METHOD: In this study, we first constructed a Foxi3 deficiency (Foxi3-/- ) mouse model to verify the craniofacial phenotype of Foxi3-/- mice, and then used RNAseq data for gene differential expression analysis to screen candidate pathogenic genes, and conducted gene expression verification analysis using quantitative real-time PCR. RESULTS: By observing the phenotype of Foxi3-/- mice, we found that craniofacial dysmorphism was present. The results of comprehensive bioinformatics analysis suggested that the craniofacial dysmorphism caused by Foxi3 deficiency may be involved in the PI3K-Akt signaling pathway. Quantitative real-time PCR results showed that the expression of PI3K-Akt signaling pathway-related gene Akt2 was significantly increased in Foxi3-/- mice. CONCLUSION: The craniofacial dysmorphism caused by the deficiency of Foxi3 may be related to the expression of Akt2 and PI3K-Akt signaling pathway. This study laid a foundation for understanding the function of FOXI3 and the pathogenesis and treatment of related craniofacial dysmorphism caused by FOXI3 dysfunction.

Low-Toxicity and High-Efficiency Streptomyces Genome Editing Tool Based on the Miniature Type V-F CRISPR/Cas Nuclease AsCas12f1.

Genome editing tools based on SpCas9 and FnCpf1 have facilitated strain improvements for natural product production and novel drug discovery in Streptomyces. However, due to high toxicity, their editing requires high DNA transformation efficiency, which is unavailable in most streptomycetes. The transformation efficiency of an all-in-one editing tool based on miniature Cas nuclease AsCas12f1 was significantly higher than those of SpCas9 and FnCpf1 in tested streptomycetes, which is due to its small size and weak DNA cleavage activity. Using this tool, in Streptomyces coelicolor, we achieved 100% efficiency for single gene or gene cluster deletion and 46.7 and 40% efficiency for simultaneous deletion of two genes and two gene clusters, respectively. AsCas12f1 was successfully extended to Streptomyces hygroscopicus SIPI-054 for efficient genome editing, in which SpCas9/FnCpf1 does not work well. Collectively, this work offers a low-toxicity, high-efficiency genome editing tool for streptomycetes, particularly those with low DNA transformation efficiency.

Genome-Wide Mining of CULLIN E3 Ubiquitin Ligase Genes from Uncaria rhynchophylla.

CULLIN (CUL) protein is a subtype of E3 ubiquitin ligase that is involved in a variety of biological processes and responses to stress in plants. In Uncaria rhynchophylla, the CUL gene family has not been identified and its role in plant development, stress response and secondary metabolite synthesis has not been studied. In this study, 12 UrCUL gene members all contained the typical N-terminal domain and C-terminal domain identified from the U. rhynchophylla genome and were classified into four subfamilies based on the phylogenetic relationship with CULs in Arabidopsis thaliana. They were unevenly distributed on eight chromosomes but had a similar structural composition in the same subfamily, indicating that they were relatively conserved and potentially had similar gene functions. An interspecific and intraspecific collinearity analysis showed that fragment duplication played an important role in the evolution of the CUL gene family. The analysis of the cis-acting elements suggests that the UrCULs may play an important role in various biological processes, including the abscisic acid (ABA) response. To investigate this hypothesis, we treated the roots of U. rhynchophylla tissue-cultured seedlings with ABA. The expression pattern analysis showed that all the UrCUL genes were widely expressed in roots with various expression patterns. The co-expression association analysis of the UrCULs and key enzyme genes in the terpenoid indole alkaloid (TIA) synthesis pathway revealed the complex expression patterns of 12 UrCUL genes and some key TIA enzyme genes, especially UrCUL1, UrCUL1-likeA, UrCUL2-likeA and UrCUL2-likeB, which might be involved in the biosynthesis of TIAs. The results showed that the UrCULs were involved in the response to ABA hormones, providing important information for elucidating the function of UrCULs in U. rhynchophylla. The mining of UrCULs in the whole genome of U. rhynchophylla provided new information for understanding the CUL gene and its function in plant secondary metabolites, growth and development.

Genetic analysis of the forkhead transcriptional factor 2 gene in three Chinese families with blepharophimosis syndrome.

PURPOSE: Clinically, blepharophimosis syndrome (BPES) has been divided into two subsets according to the association of ocular malformation with (type I) or without (type II) premature ovarian failure (POF). BPES is ascribed to mutations in the forkhead transcriptional factor 2 (FOXL2) gene. This study aimed at identifying clinical features and mutations within the FOXL2 gene in three Chinese families with BPES. METHODS: A clinical and molecular genetic investigation was performed in affected and unaffected members from three Chinese families with BPES. Genomic DNA was prepared from leucocytes of peripheral venous blood, the entire coding region of FOXL2 were amplified with PCR, and direct DNA sequencing of the PCR products was performed for mutations in FOXL2. RESULTS: Three mutations in FOXL2 were found in three families, including c.672_701dup30, c.663_692dup30, and c.475dupC. Of the three, the c.475dupC (p.His159fs) was novel in family C and resulted in a frameshift mutation to generate a truncated protein owing to a premature stop codon at codon 238. The new duplication mutation was associated with BPES type II. The c.672_701dup30 (p.Ala224_Ala234dup10) and the c.663_692dup30 (p.Ala221_Ala231dup10) were detected in family A and family B, respectively, leading to expansions of the polyalanine (poly-Ala) tract that is frequently the hot spot of mutations within FOXL2. CONCLUSIONS: Our results expand the spectrum of FOXL2 mutations, and further indicate the association of a novel duplication mutation leading to a truncated protein with BPES type II. The other two known mutations may support the previous hypothesis regarding expansions of the polyalanine tract associated with BPES type II as a mutational hot spot in FOXL2.

[Genetic diagnosis for a Chinese Han family with hereditary multiple osteochondromas].

OBJECTIVE: To identify the genetic cause for a Chinese Han family affected with hereditary multiple osteochondromas. METHODS: Two patients, five unaffected relatives of the family and 100 unrelated healthy controls were collected. The coding sequences and intron/exon boundaries of EXT1 gene were amplified with polymerase chain reaction (PCR) and sequenced. RESULTS: A heterozygous c.600G>A (p.Trp200X) mutation in exon 1 of the EXT1 gene was detected in the patients. The same mutation was not found in unaffected family members and 100 healthy controls. CONCLUSION: The hereditary multiple osteochondromas in the family is caused by a nonsense mutation (p.Trp200X) in the EXT1 gene.

Genome-wide identification and co-expression network analysis of Aux/IAA gene family in Salvia miltiorrhiza.

The auxin/indole-3-acetic acid (Aux/IAA) gene family serves as a principal group of genes responsible for modulating plant growth and development through the auxin signaling pathway. Despite the significance of this gene family, the identification and characterization of members within the well-known Chinese medicinal herb Salvia miltiorrhiza (S. miltiorrhiza) have not been thoroughly investigated. In this study, we employed bioinformatics methods to identify 23 Aux/IAA genes within the genome of S. miltiorrhiza. These genes were classified into typical IAA and atypical IAA based on their domain structure. Our analysis of the promoter regions revealed that the expression of these genes is regulated not only by auxins, but also by other hormones and environmental factors. Furthermore, we found that the expression patterns of these genes varied across various tissues of S. miltiorrhiza. While our initial hypothesis suggested that the primary function of these genes was the interaction between SmIAA and ARF, gene co-expression network analysis revealed that they are also influenced by various other transcription factors, such as WRKY and ERF. The findings establish a sturdy basis for future investigations into the function of the Aux/IAA gene family and exhibit promising prospects for enhancing the genetics of this medicinal flora and its associated species.

Chromosome-level genome assembly of the endangered plant Tetraena mongolica.

Tetraena mongolica is an endangered xerophytic shrub with high ecological value for the restoration of desert vegetation because of its high tolerance to drought and heat stress. Here, we generated a high-quality chromosome-level reference genome of T. mongolica by combining PacBio HiFi data and Hi-C sequencing technologies, which was approximately 1.12 Gb (contig N50 of 25.5 Mb) in size and contained 61,888 protein-coding genes; repetitive sequences comprised 44.8% of the genome. This genome of T. mongolica is the first published genome sequence of a member of the order Zygophyllales. Genome analysis showed that T. mongolica has undergone a recent whole genome duplication event, and a recent burst of long terminal repeat insertions afterward, which may be responsible for its genome size expansion and drought adaptation. We also conducted searches for gene homologues and identified terpene synthase (TPS) gene families and candidate genes involved in triacylglycerol biosynthesis. The T. mongolica genome sequence could aid future studies aimed at functional gene identification, germplasm resource management, molecular breeding efforts, as well as evolutionary studies of Fabids and angiosperm taxa.

Chloroplast DNA reveals genetic population structure in Sinomenium acutum in subtropical China.

Objective: The population density and diversity of Sinomenium acutum (Menispermaceae) have been greatly reduced recently by overharvesting for medicinal purposes in China. Therefore, it is urgent that the remaining populations are investigated, and that strategies for the utilization and conservation of this species are developed. This study aimed to find the possible glacial refugia and define the genetic diversity of S. acutum for its proper utilization and conservation. Methods: A total of 77 S. acutum samples were collected from four locations, Qinling Mountains, Daba Mountains, Dalou Mountains, and Xuefeng Mountains, in subtropical China. Genetic diversity among and between these populations were phylogenetically analyzed using four chloroplast DNA molecular markers (atpI-atpH, trnQ-5'rps16, trnH-psbA and trnL-trnF). Results: A total of 14 haplotypes (C1 to C14) were found in collected samples. Haplotypes C1 and C3 were shared among all populations, with C3 as the ancestral haplotype. Haplotypes C11 and C12 diverged the most from C3 and other haplotypes. No obvious phylogeographic structure was found in four locations using the GST/NST test. There is no evidence of rapid demographic expansion in S. acutum based on the mismatch distribution, and the results of Tajima's D test, and Fu's FS test. Our analyses of molecular variance revealed a high level of genetic variation within populations. In contrast, the genetic differentiation among S. acutum populations was low, indicating frequent gene flow. Conclusion: Xuefeng, Dalou, and Daba Mountains were possible glacial refugia for the populations of S. acutum. C1, C3, C11 and C12 haplotypes of S. acutum should be carefully preserved and managed for their genetic value.

Diaporpyrone E, an undescribed α-pyrone from the endophytic fungus Diaporthe sp. CB10100.

An undescribed α-pyrone diaporpyrone E (1), and three known nucleotides, 5'-O-acetyl uridine (2), 5'-O-acetyl thymidine (3), and adenine (4), were identified from Diaporthe sp. CB10100, an endophytic fungus isolated from the medicinal plant Sinomenium acutum. The structure of 1 was determined by extensive analysis of its HRMS, 1D and 2D NMR spectroscopic data, as well as electronic circular dichroism calculations and comparison. The in vitro cytotoxic and antibacterial assays of 1 revealed that it has a 30.2% inhibitory effect on HepG2 cells at 50 µM, while no antibacterial activities against Staphylococcus aureus and Klebsiella pneumoniae at 64 µg/mL.

Characterization of two putative norlaudanosoline methyltransferases from Aristolochia debilis.

In view of the nephrotoxicity, hepatotoxicity, and carcinogenicity of aristolochic acids (AAs), the removal of AAs from plants becomes an urgent priority for ensuring the safety of Aristolochia herbal materials. In this study, based on the root-predominant distribution of aristolochic acid I (AAI) in Aristolochia debilis, transcriptome sequencing, in combination with phylogenetic analyses, and gene expression pattern analysis together provided five candidate genes for investigating AAI biosynthesis. Comprehensive in vitro and in vivo enzymatic assays revealed that Ab6OMT1 (6-O-methyltransferase) and AbNMT1 (N-methyltransferase) exhibit promiscuity in substrate recognition, and they could act in a cooperative fashion to achieve conversion of norlaudanosoline, a predicted intermediate in AAI biosynthetic route, into 3'-hydroxy-N-methylcoclaurine through two different methylation reaction sequences. These results shed light on the molecular basis for AAI biosynthesis in Aristolochia herbs. More importantly, Ab6OMT1 and AbNMT1 may be employed as targets for the metabolic engineering of AAI biosynthesis to produce AAs-free Aristolochia herbal materials.

[Clinical investigation of a Chinese family with hypotrichosis simplex of the scalp and mutational analysis of CDSN gene].

OBJECTIVE: To analyze clinical symptoms and disease-causing mutations of corneodesmosin (CDSN) gene in a Chinese family affected with hypotrichosis simplex of the scalp and to establish a method for prenatal diagnosis. METHODS: Family survey and clinical examinations were carried out to determine the inheritance pattern. Three patients and 7 unaffected relatives from the family, in addition with 100 unrelated healthy controls were recruited. Genomic DNA from peripheral blood leukocytes was extracted. Five pairs of primers were designed based on the CDSN gene sequence. Exons and flanking regions of the CDSN gene were amplified using polymerase chain reaction (PCR). Potential mutations were analyzed through direct sequencing and comparison by BLAST. RESULTS: The type of alopecia of the family was diagnosed as hypotrichosis simplex of the scalp with an autosomal dominant inheritance pattern. A nonsense mutation (C717G) in cDNA sequence of the CDSN gene was identified in all three patients of the family, which resulted in a premature stop codon (Y239X). The same mutation was not found among healthy members of the family and 100 healthy controls. CONCLUSION: A Chinese family was diagnosed with hypotrichosis simplex of the scalp, which was caused by a novel nonsense mutation (Y239X) in the CDSN gene.

Full-length transcriptome and metabolite analysis reveal reticuline epimerase-independent pathways for benzylisoquinoline alkaloids biosynthesis in Sinomenium acutum.

Benzylisoquinoline alkaloids (BIAs) are a large family of plant natural products with important pharmaceutical applications. Sinomenium acutum is a medicinal plant from the Menispermaceae family and has been used to treat rheumatoid arthritis for hundreds of years. Sinomenium acutum contains more than 50 BIAs, and sinomenine is a representative BIA from this plant. Sinomenine was found to have preventive and curative effects on opioid dependence. Despite the broad applications of S. acutum, investigation on the biosynthetic pathways of BIAs from S. acutum is limited. In this study, we comprehensively analyzed the transcriptome data and BIAs in the root, stem, leaf, and seed of S. acutum. Metabolic analysis showed a noticeable difference in BIA contents in different tissues. Based on the study of the full-length transcriptome, differentially expressed genes, and weighted gene co-expression network, we proposed the biosynthetic pathways for a few BIAs from S. acutum, such as sinomenine, magnoflorine, and tetrahydropalmatine, and screened candidate genes involved in these biosynthesis processes. Notably, the reticuline epimerase (REPI/STORR), which converts (S)-reticuline to (R)-reticuline and plays an essential role in morphine and codeine biosynthesis, was not found in the transcriptome data of S. acutum. Our results shed light on the biogenesis of the BIAs in S. acutum and may pave the way for the future development of this important medicinal plant.

[Mutation analysis of FGFR3 gene in a family featuring hereditary dwarfism].

OBJECTIVE: To investigate the clinical symptoms and potential mutation in FGFR3 gene for a family featuring hereditary dwarfism in order to attain diagnosis and provide prenatal diagnosis. METHODS: Five patients and two unaffected relatives from the family, in addition with 100 healthy controls, were recruited. Genome DNA was extracted. Exons 10 and 13 of the FGFR3 gene were amplified using polymerase chain reaction (PCR). PCR products were sequenced in both directions. RESULTS: All patients had similar features including short stature, short limbs, lumbar hyperlordosis but normal craniofacial features. A heterozygous mutation G1620T (N540K) was identified in the cDNA from all patients but not in the unaffected relatives and 100 control subjects. A heterozygous G380R mutation was excluded. CONCLUSION: The hereditary dwarfism featured by this family has been caused by hypochondroplasia (HCH) due to a N540K mutation in the FGFR3 gene.

Biosynthetic Pathway and the Potential Role of Melatonin at Different Abiotic Stressors and Developmental Stages in Tolypocladium guangdongense.

Melatonin, a bioactive compound and an important signaling molecule produced in plants and animals, is involved in many biological processes. However, its function and synthetic pathways in fungi are poorly understood. Here, the samples from Tolypocladium guangdongense, a highly valued edible fungus with functional food properties, were collected under different experimental conditions to quantify the levels of melatonin and its intermediates. The results showed that the intracellular melatonin content was markedly improved by Congo red (CR), cold, and heat stresses; the levels of intracellular melatonin and its intermediates increased at the primordial (P) and fruiting body (FB) stages. However, the levels of most intermediates exhibited a notable decrease under CR stress. Several genes related to melatonin synthesis, excluding AADC (aromatic-L-amino-acid decarboxylase), were markedly upregulated at an early stage of CR stress but downregulated later. Compared to the mycelial stage, those genes were significantly upregulated at the P and FB stages. Additionally, exogenous melatonin promoted resistance to several abiotic stressors and P formation in T. guangdongense. This study is the first to report melatonin biosynthesis pathway in macro-fungi. Our results should help in studying the diversity of melatonin function and melatonin-synthesis pathways and provide a new viewpoint for melatonin applications in the edible-medicinal fungus.