RNA-seq transcriptome and pathway analysis of the medicinal mushroom Lignosus tigris (Polyporaceae) offer insights into its bioactive compounds with anticancer and antioxidant potential

Medicinal mushrooms belonging to the Lignosus spp., colloquially known as Tiger Milk mushrooms (TMMs), are used as traditional medicine by communities across various regions of China and Southeast Asia to enhance immunity and to treat various diseases. At present, three Lignosus species have been identified in Malaysia: L. rhinocerus, L. tigris, and L. cameronensis. Similarities in their macroscopic morphologies and the nearly indistinguishable appearance of their sclerotia often lead to interchangeability between them. Hence, substantiation of their traditional applications via identification of their individual bioactive properties is imperative in ensuring that they are safe for consumption. L. tigris was first identified in 2013. Thus far, studies on L. tigris cultivar sclerotia (Ligno TG-K) have shown that it possesses significant antioxidant activities and has greater antiproliferative action against selected cancer cells in vitro compared to its sister species, L. rhinocerus TM02®. Our previous genomics study also revealed significant genetic dissimilarities between them. Further -omics investigations on Ligno TG-K hold immense potential in facilitating the identification of its bioactive compounds and their associated bioactivities. The overall aim of this study was to investigate the gene expression profile of Ligno TG-K via de novo RNA-seq and pathway analysis. We also aimed to identify highly expressed genes encoding compounds that contribute to its cytotoxic and antioxidant properties, as well as perform a comparative transcriptomics analysis between Ligno TG-K and its sister species, L. rhinocerus TM02®. Total RNA from fresh 3-month-old cultivated L. tigris sclerotia (Ligno TG-K) was extracted and analyzed via de novo RNA sequencing. Expressed genes were analyzed using InterPro and NCBI-Nr databases for domain identification and homology search. Functional categorization based on gene functions and pathways was performed using Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Clusters of Orthologous Genes (COG) databases. Selected genes were subsequently subjected to phylogenetic analysis. Our transcriptomics analysis of Ligno TG-K revealed that 68.06% of its genes are expressed in the sclerotium; 80.38% of these were coding transcripts. Our analysis identified highly expressed transcripts encoding proteins with prospective medicinal properties. These included serine proteases (FPKM = 7356.68), deoxyribonucleases (FPKM = 3777.98), lectins (FPKM = 3690.87), and fungal immunomodulatory proteins (FPKM = 2337.84), all of which have known associations with anticancer activities. Transcripts linked to proteins with antioxidant activities, such as superoxide dismutase (FPKM = 1161.69) and catalase (FPKM = 1905.83), were also highly expressed. Results of our sequence alignments revealed that these genes and their orthologs can be found in other mushrooms. They exhibit significant sequence similarities, suggesting possible parallels in their anticancer and antioxidant bioactivities. This study is the first to provide a reference transcriptome profile of genes expressed in the sclerotia of L. tigris. The current study also presents distinct COG profiles of highly expressed genes in Ligno TG-K and L. rhinocerus TM02®, highlighting that any distinctions uncovered may be attributed to their interspecies variations and inherent characteristics that are unique to each species. Our findings suggest that Ligno TG-K contains bioactive compounds with prospective medicinal properties that warrant further investigations. Systems biology and omics