Priyanka JABB 2024
Priyanka JABB 2024
Priyanka JABB 2024
Priyanka N. Shenoy1, Sneha Bhaskar1, M. Manu1, M. P. Likitha1, N. Geetha1, Shailasree Sekhar2*, K Ramachandra Kini1
1
Department of Studies in Biotechnology, University of Mysore, Mysuru, Karnataka, India.
2
Division of Biochemistry, School of Life Sciences, JSS Academy of Higher Education and Research, Sri Shivarathreeshwara Nagara, Mysuru, Karnataka, India.
Article history:
Medicinal plants and their endophytes are one of the efficient producers of diverse secondary metabolites with
Received on: November 29, 2023
therapeutic importance. In the present study, an endophytic fungus Talaromyces spp. isolated from Syzygium
Accepted on: February 27, 2024
samarangense was subjected to whole-genome sequencing and Antibiotics and Secondary Metabolite Analysis
Available online: April 20, 2024
Shell (AntiSMASH) annotation to identify biosynthetic gene clusters (BGCs) of secondary metabolites and their
biosynthetic pathways. The Funannotate results revealed that the Talaromyces isolate has a total of 30.5Mb genome
Key words:
consisting of 372 contigs, 372 scaffolds and 47.64% GC content. In addition, 114 tRNA, 12722 functional mRNA,
Endophytic fungi,
Secondary metabolites, 12721 CDS transcripts, and 12721 protein coding sequences were predicted and annotated using various BLAST
Antibiotics and secondary metabolite databases. AntiSMASH revealed the presence of 76 BGCs, including 28 T1 Polyketide synthase (T1 PKS), 10
analysis shell, Nonribosomal peptide synthetases (NRPSs), 9 terpene, 1 Indole, 12 NRPs like, 4 T1 PKS and NRPs like, 2 T1
Talaromyces, PKS Indole, 1 NRPs like terpene, 3 NRPs T1 PKS, 1 Indole NRPs, 2 betalactone, 1 phosphonate, 1 fungal-RiPP T1
Biosynthetic gene clusters, PKS and 1 Other type. The analysis also predicted the occurrence of enzymes involved in the biosynthesis of some
Whole genome Sequence. of the important secondary metabolites such as Pyranonigrin, Squalestatin -S, Azanigerone -A, Asperterpenoid -A,
Napthopyrone, Clavaric acid, and Fusarin.
© 2024 Priyanka N. Shenoy, et al. This is an open access article distributed under the terms of the Creative Commons Attribution License -NonCommercial-
ShareAlike Unported License (http://creativecommons.org/licenses/by-nc-sa/3.0/).
230 Shenoy, et al.: Journal of Applied Biology & Biotechnology 2024;12(3):229-235
sterile distilled water wash and treated with 1% sodium hypochlorite 2.5. Genome Mining for Natural Product Discovery
solution for 4 min. They were again washed with sterile distilled AntiSMASH tool: Fungal version 6.1.1 pipeline was used for the
water for 2–3 times, and the leaf bits were placed on solidified potato identification of potent secondary metabolite BGCs [16].
dextrose agar (PDA) media supplemented with ampicillin 100 mg/mL.
The Petri dishes were sealed with parafilmTM and incubated at 25 3. RESULTS AND DISCUSSION
± 2°C for 15 days with alternate light and dark cycles. Endophytic
fungal mycelia which were observed on leaf and stem explants, 3.1. YS-1 Strain Identification
were subcultured by point inoculation on PDA and potato dextrose The cultured fungal organism showed yellow-colored cottony mycelial
broth media. Samples were monitored for further mycelial growth mat with white dots covering the entire plate after being incubated
and were subjected to staining with lactophenol cotton blue in their for 15 days at 25 ± 2°C. Annotation and phylogenetic analysis results
sporulation state. Among the four endophytic fungal strains isolated showed that the organism belongs to Talaromyces spp. and the
from S. samarangense, strain YS-1 was selected for the present study. stereomicroscopic view of strain YS-1 was examined [Figure 1a and
and b].
2.2. Genomic DNA isolation and WGS of YS1 strain 3.2. Genome Assembly and Sequencing
DNA was isolated from the fungal mycelia by CTAB method [8]. The genome sequence of strain YS-1 was deposited in the NCBI GenBank
Isolated DNA was eluted with 50 µL of TE buffer and stored at −20°C. database (Accession No. PRJNA953512), and the Genome assembly
The purity of isolated DNA was determined using NanoDrop™ 2000C of Talaromyces spp. was accessed by assemblers such as Biosynthetic
spectrophotometer (Thermo Fisher Scientific). Genomic DNA samples spades, Spades, and Megahit. These assemblers were compared and
were electrophoresed on 1% agarose gel to check purity and integrity. analyzed using Quast (Quality assessment tool for genome assemblies).
The high-quality DNA of strain YS1 was sent for whole genome Megahit was found to be an effective assembler that provided data on an
sequence analysis to Hi-Gx360® Himedia Laboratories Pvt. Ltd., organism’s complex metagenomic assemblies. The whole genome size
Mumbai, India. The sample analysis was as follows: 250 ng of total of Talaromyces spp. was found to be 30.5Mb, consisting of 605 contigs
DNA was used as input for library preparation using QIASeq FX DNA (≥0 bp), 372 contigs (≥1000 bp), 208 contigs (≥5000 bp), 168 contigs
kit to fragment and obtain adapter ligated and indexed library as per (≥10000 bp), 149 (≥25000 bp), 125 contigs (≥50,000 bp) and the largest
manufacturer’s instructions. The indexed library was sequenced on an contig was found to be 1112360. The GC content was 46.74%, N50
Illumina Miseq by paired-end chemistry of 300 cycles. The parameters was 3.33Mb, N75 was 1.86Mb. L50 value was 30 and the L75 value
used to check the raw data quality are included in a MultiQC report was 65. The Funannotate analysis indicated the presence of 114tRNA,
and are used to understand the quantity of data obtained for each of the 12722 functional mRNA, 12721 CDS transcripts, and 12721 protein-
paired read files for an individual sample [9]. The fastp tool (v0.12.4) coding sequences in Talaromyces spp. YS-1 strain. The completeness
was used to remove adapter contamination [10]. of genome assembly, as measured by BUSCO, gave a score of 92.9%.
2.3. Genome Assembly and Annotations 3.3. Phylogenetic Analysis and Phylogenomic Tree
A de novo de-brujin graph-based assembly was carried out to assemble The phylogenetic analysis of the Talaromyces spp. sequenced in the
the short reads into larger stretches of DNA called contigs. These present study (YS-1 Strain) was compared with the isolates available
contigs are the starting material used for executing a genome annotation in the database. Three groups could be identified in the phylogenetic
that assigns functions to various regions of the genome of the organism
in question. Quality assessment for genome assemblies generated by
Spades, bioSpades, and Megahit assemblers was performed using
the Quast tool [11-13]. It was noticed that the assembly generated by
the megahit assembler had a comparatively higher number of contigs
with sizes >10kb. Hence, assembly generated in bioSpades mode used
for phylogenomics analysis and assembly generated using megahit
assembler was masked and further used for annotations and natural
product discovery.
tree, with the Talaromyces YS-1 strain forming a close relationship applications in pharmaceuticals and drug discovery can be analyzed
with Talaromyces amestolkiae, Talaromyces marneffei, and using this comprehensive analysis tool for the fungal genome.
Talaromyces stipitatus in the first group. The second group consisted
AntiSMASH analysis showed the presence of 76 secondary metabolite
of Talaromyces proteolyticus and Talaromyces rugulosus. The other
BGCs, which include 28 T1Polyketide synthases, 10 Nonribosomal
species Talaromyces atroroseus, formed a separate group [Figure 2].
peptide synthetases (NRPSs), 9 terpene, 1 Indole, 12 NRPs like 4 T1
Polyketide synthase (T1 PKS) and NRPs like, 2 T1 PKS Indole, 1
3.4. Genome Mining for Natural Product Discovery NRPs like terpene, 3 NRPs T1 PKS, 1 Indole NRPs, 2betalactone 1
AntiSMASH provides information about the secondary metabolite phosphonate, 1 fungal-RiPP T1 PKS and 1 Other type [Table 1]. Over
genes by comparing the query gene cluster with the BGC library. 40% of these BGCs indicated gene homologies with known clusters in
Potential gene clusters of certain secondary metabolites such as the MIBiG database. Several clusters responsible for coding bioactive
terpenes, alkaloids, and phenolic compounds having industrial compounds such as Pyranonigrin-A, Alternariol, Asperterenoid-A,
h
Figure 3: Comparison of BGCs in Talaromyces spp. YS-1 (a-h) with identified BGCs for biosynthesis of various bioactive compounds (ⅰ-ⅷ).
232 Shenoy, et al.: Journal of Applied Biology & Biotechnology 2024;12(3):229-235
(Contd...)
Shenoy, et al.: Genome mining of Talaromyces sp. for BGCs 2024;12(3):229-235 233
Table 1: (Continued).
Region From To Type Most similar known cluster
Region 34.1 52,652 149,220 Terpene Unknown
Region 37.1 78,840 267,110 NRPS-like Unknown
Region 38.1 128,355 424,526 NRPS, T1PKS Unknown
Region 39.1 111,254 271,795 NRPS-like Unknown
Region 40.1 105,708 141,371 Terpene Unknown
Region 41.1 1 349,214 Phosphonate Unknown
Region 41.2 90,852 224,543 NRPS-like Unknown
Region 47.1 27,647 171,848 Betalactone Unknown
Region 48.1 2189 74,424 T1PKS Naphthopyrone (100%)
Region 49.1 67,645 122,148 fungal-RiPP, T1PKS Ankaflavin/monascin/rubropunctatine/monascorubrin (20%)
Region 54.1 13,165 179,953 T1PKS Pyranonigrin E (100%)
Region 54.2 116,565 155,050 Terpene Unknown
Region 57.1 19,253 127,159 T1PKS, NRPS-like Unknown
Region 59.1 196,661 35,955 Terpene Unknown
Region 60.1 49,877 133,974 Terpene Sordarin (11%)
Region 63.1 163,422 51,896 Terpene Clavaric acid (100%)
Region 67.1 92,953 48,783 T1PKS Chrysoxanthone A/chrysoxanthone B/chrysoxanthone C (8%)
Region 68.1 132,828 174,717 NRPS-like Unknown
Region 69.1 4601 61,197 T1PKS Melanin (100%)
Region 69.2 105,785 139,510 T1PKS Fusarin (100%)
Region 71.1 502 72,951 NRPS Unknown
Region 71.2 73,873 210,071 NRPS Unknown
Region 73.1 1 71,552 Terpene Unknown
Region 97.1 1 185,810 T1PKS Sorbicillin (71%)
Region 98.1 21,324 147,843 T1PKS, NRPS-like Shearinine D (9%)
Region 107.1 27,684 176,667 T1PKS Chrodrimanin B (100%)
Region 127.1 3370 48,170 NRPS-like Unknown
Region 129.1 3888 154,701 T1PKS Depudecin (33%)
Region 130.1 1 43,527 NRPS Unknown
Region 139.1 1 119,563 T1PKS, NRPS-like Fusarin (100%)
Naphthopyrone, Clavaric acid, Melanin, and Fusarin showed 100% from Hypholomas ublateritium (Gen Bank EU665687.1) for synthesis
similarity. of Clavaric acid, a potential anti-tumor compound [21] [Figure 3e].
Biosynthetic cluster region 69.1 showed homology with Melanin
AntiSMASH analysis of Talaromyces spp. strain YS-1indicated that the BGC [Figure 3f]. Melanin is found predominantly in PKS clusters of
gene clusters in the regions 1.1, 31.1, and 54.1 has significant similarity Fungi and is useful for adjusting to different environmental conditions
with the BGCs of Pyranonigrin E from Aspergillus niger ATCC 1015 [22]. The fungal YS-1 strain in regions 69.2 and 139.1 encodes BGC
(GenBank: ACJE01000019.1.) [Figure 3a]. Pyranonigrin- E, isolated with BLAST hits for Fusarin BGC from Fusarium fujikuroi (Gen
from Aspergillus niger, is a pyranopyrrole having effective anti-oxidant Bank: JX308619.1) [Figure 3g]. Fusarin BGCs were been identified
and radical scavenging activities [17]. BGC region 5.1 noted significant in Fusarium spp., along with different polyketide gene clusters [23].
BLAST hits with Alternariol from Parastagono sporanodorum SN15 BGC region 107.1 of fungal YS-1 strain presented 100% similarity with
(GenBank: KP941080.1) [Figure 3b]. Alternariol is reported to be BGC from Talaromyces verruculosus (Gen Bank: LC422696.1) meant
an effective anti-cancerous and antibacterial metabolite [18]. BGC for biosynthesis of secondary metabolite chrodrimanin B [Figure 3h].
13.1 exhibited high similarity with known clusters implicated in Different classes of secondary metabolites and a new monoterpenoid
the synthesis of Asperterpenoid-A from Talaromyces wortmannii Chrodrimanin-T have been reported from T. amestolkiae [24].
(GenBank: MK140602.1) [Figure 3c]. Asperterpenoid-A has effective Chrodrimanin-B- B isolated from Talaromyces funiculosus showed a
anti-inflammatory and antibacterial activities [19]. BGCs found in broad range of anti-bactericidal activity against Staphylococcus aureus,
region 48.1 of the YS-1 strain showed similarity with Naphthopyrone Mycobacterium smegmatis, Mycobacterium phlei, Micrococcus
BGC from Aspergillus nidulans FGSC A4 (Gen Bank: BN001302.1) tetragenus and Escherichia coli [25]. In addition, several other
[Figure 3d]. Li et al. [20] have reported that naphthopyrone is bioactives and Antibiotics such as Terretonin, Zearalenone, Trypacidin,
a potential antiproliferative compound produced from A. niger. Viriditoxin, Tricholignan -A, Squalestatin -S, Emericellamide -A,B,
A notable BLAST hit for the BGC region 63.1 was similar to BGC Fumigaclavine -C, Azanigerone -A, Curvupallide -B, Griseofulvin,
234 Shenoy, et al.: Journal of Applied Biology & Biotechnology 2024;12(3):229-235
Citrinin, Clapurines, Fujikurin A,B,C,D, Ankaflavin, Sordarin, experiments and collected the data. PSN, SB, SS and KRK analysed
Chrysoxanthone -A,B,C, Sorbicillin, Shearinine -D, and Depudecin the data and finalized the manuscript. All authors have read and
have been identified in the present study. approved the final manuscript.
14. Lau YL, Lee WC, Xia J, Zhang G, Razali R, Anwar A, et al. Draft guided discovery of a new polyketide from Talaromyces funiculosus
genome of Brugiapahangi: High similarity between B. pahangi and Salicorn 58. Appl Microbiol Biotechnol 2016;100:5323-38.
B. malayi. Parasites Vectors 2015;8:1-14. 26. Schafhauser T, Wibberg D, Rückert C, Winkler A, Flor L,
15. Waterhouse RM, Seppey M, Simao FA, Manni M, Ioannidis P, van Pée KH, et al. Draft genome sequence of Talaromyces islandicus
Klioutchnikov G, et al. BUSCO applications from quality (“Penicillium islandicum”) WF-38-12, a neglected mold with
assessments to gene prediction and phylogenomics. Mol Biol significant biotechnological potential. J Biotechnol 2015;211:101-2.
Evol 2018;35:543-8. 27. Chen S, Liu Y, Liu Z, Cai R, Lu Y, Huang X, et al. Isocoumarins
16. Medema MH, Blin K, Cimermancic P, De Jager V, Zakrzewski P, and benzofurans from the mangrove endophytic fungus Talaromyces
Fischbach MA, et al. AntiSMASH: Rapid identification, annotation amestolkiae possess α-glucosidase inhibitory and antibacterial
and analysis of secondary metabolite biosynthesis gene clusters activities. RSC Adv 2016;6:26412-20.
in bacterial and fungal genome sequences. Nucleic Acids Res 28. El-Elimat T, Figueroa M, Raja HA, Alnabulsi S, Oberlies NH.
2011;39(Suppl 2):339-46. Coumarins, dihydroisocoumarins, a dibenzo-α-pyrone, a
17. Riko R, Nakamura H, Shindo K. Studies on pyranonigrins-isolation meroterpenoid, and a merodrimane from Talaromyces amestolkiae.
of pyranonigrin E and biosynthetic studies on pyranonigrin A. Tetrahedron Lett 2021;72:153067.
J Antibiot (Tokyo) 2014;67:179-81. 29. Singh SK, Mondal A, Saha N, Husain SM. Identification and
18. Palanichamy P, Kannan S, Murugan D, Alagusundaram P, characterization of an anthrol reductase from: Talaromyces islandicus
Marudhamuthu M. Purification, crystallization and anticancer activity (Penicillium islandicum) WF-38-12. Green Chem 2019;21:6594-9.
evaluation of the compound alternariol methyl ether from endophytic 30. Zhou C, Cao X, Ge Y, Wu X, Zhang Z, Ma Y, et al. Talaropeptins
fungi Alternaria alternata. J Appl Microbiol 2019;127:1468-78. A and B, tripeptides with an N- trans-cinnamoyl moiety from the
19. Yan C, Liu W, Li J, Deng Y, Chen S, Liu H. Bioactive terpenoids marine-derived fungus Talaromyces purpureogenus CX11. J Nat
from: Santalum album derived endophytic fungus Fusarium sp. YD- Prod 2022; 85:2620-5.
2. RSC Adv 2018;8:14823-8. 31. Li CX, Zhao S, Zhang T, Xian L, Liao LS, Liu JL, et al. Genome
20. Li DH, Han T, Guan LP, Bai J, Zhao N, Li ZL, et al. New sequencing and analysis of Talaromyces pinophilus provide insights
naphthopyrones from marine-derived fungus Aspergillus Niger into biotechnological applications. Sci Rep 2017;7:490.
2HL-M-8 and their in vitro antiproliferative activity. Nat Prod Res 32. Rashied RM, Abdelfattah MA, El-Beshbishy HA, ElShazly AM,
2016;30:1116-22. Mahmoud MF, Sobeh M. Syzygium samarangense leaf extract
21. Godio RP, Martín JF. Modified oxidosqualene cyclases in the exhibits distinct antidiabetic activities: Evidences from in silico and
formation of bioactive secondary metabolites: Biosynthesis of the in vivo studies: Syzygium samarangense leaf extract exhibits distinct
antitumor clavaric acid. Fungal Genet Biol 2009;46:232-42. antidiabetic activities. Arab J Chem 2022;15:103822.
22. Jia SL, Chi Z, Chen L, Liu GL, Hu Z, Chi ZM. Molecular evolution 33. Kaul S, Gupta S, Ahmed M, Dhar MK. Endophytic fungi from
and regulation of DHN melanin-related gene clusters are closely medicinal plants: A treasure hunt for bioactive metabolites.
related to adaptation of different melanin-producing fungi. Genomics Phytochem Rev 2012;11:487-505.
2021;113:1962-75.
23. Desjardins AE, Proctor RH. Molecular biology of Fusarium How to cite this article:
mycotoxins. Int J Food Microbiol 2007;119:47-50. Shenoy PN, Bhaskar S, Manu M, Likitha MP, Geetha N, Shailasree S,
24. Fu Y, Li C, Zhu J, Zhang L, Wang Y, Chen Q, et al. A new Kini KR. Genome mining and AntiSMASH analysis of an Endophytic
meroterpenoid from endophytic fungus Talaromyces amestolkiae Talaromyces sp. reveal biosynthetic pathway gene clusters for novel
CS-O-1. Biochem Syst Ecol 2020;93:104186. bioactive compounds. J App Biol Biotech. 2024;12(3):229-235. DOI: 10.7324/
25. Guo J, Ran H, Zeng J, Liu D, Xin Z. Tafuketide, a phylogeny- JABB.2024.154874