Comparative genomics of Sporothrix species and identification of putative pathogenic-gene determinants
Abstract
Aim: To understand the phylogenomics, pathogenic/virulence-associated genes and genomic evolution of pathogenic Sporothrix species. Materials & methods: We performed in silico comparative genome analysis of Sporothrix species using ab initio tools and in-house scripts. We predicted genes and repeats, compared genomes based on synteny, identified orthologous clusters, assessed genes family expansion/contraction, predicted secretory proteins and finally searched for similar sequences from various databases. Results: The phylogenomics revealed that Sporothrix species are closely related to Ophiostoma species. The gene family evolutionary analysis revealed the expansion of genes related to virulence (CFEM domain, iron acquisition genes, lysin motif domain), stress response (Su[var]3-9, Enhancer-of-zeste and Trithorax domain and Domain of unknown function 1996), proteases (aspartic protease, x-pro dipeptidyl-peptidase), cell wall composition associated genes (chitin deacetylase, chitinase) and transporters (major facilitator superfamily transporter, oligo-peptide transporter family) in Sporothrix species. Conclusion: The present study documents the putative pathogenic/virulence-associated genes in the Sporothrix species.
References
- 1. . Global epidemiology of sporotrichosis. Med. Mycol. 53(1), 3–14 (2015).
- 2. . Global ITS diversity in the Sporothrix schenckii complex. Fungal Divers. 66(1), 153–165 (2013).
- 3. . Sporotrichosis caused by Sporothrix mexicana, Portugal. Emerg. Infect. Dis. 17(10), 1975–1976 (2011).
- 4. . The divorce of Sporothrix and Ophiostoma: solution to a problematic relationship. Stud. Mycol. 83, 165–191 (2016).
- 5. An atypical cause of sporotrichosis in a cat. Med. Mycol. Case Rep. 23(November 2018), 72–76 (2019).
- 6. Sporothrix humicola (Ascomycota: Ophiostomatales) – A soil-borne fungus with pathogenic potential in the eastern quoll (Dasyurus viverrinus). Med. Mycol. Case Rep. 25(June), 39–44 (2019).
- 7. Is Sporothrix chilensis circulating outside Chile? PLoS Negl. Trop. Dis. 14(3), e0008151 (2020).
- 8. . First case report of cutaneous sporotrichosis (Sporothrix species) in a cat in the UK. J. Feline Med. Surg. Open Rep. 6(1), 1–5 (2020).
- 9. . An environmental Sporothrix as a cause of corneal ulcer. Med. Mycol. Case Rep. 2(1), 88–90 (2013).
- 10. Sporotrichosis caused by Sporothrix brasiliensis in Argentina: case report, molecular identification and in vitro susceptibility pattern to antifungal drugs. J. Mycol. Med. 30(1), 100908 (2019).
- 11. Proteins potentially involved in immune evasion strategies in Sporothrix brasiliensis elucidated by ultra-high-resolution mass spectrometry. mSphere 3(3), 1–10 (2018).
- 12. . Thermally dimorphic human fungal pathogens—polyphyletic pathogens with a convergent pathogenicity trait. Cold Spring Harb. Perspect. Med. 5(8), a019794 (2015).
- 13. . Genome analysis reveals evolutionary mechanisms of adaptation in systemic dimorphic fungi. Sci. Rep. 8(1), 4473 (2018).
- 14. Comparative genomics of the major fungal agents of human and animal Sporotrichosis: sporothrix schenckii and Sporothrix brasiliensis. BMC Genomics 15(1), 943 (2014).
- 15. Whole-genome sequencing and in silico analysis of two strains of Sporothrix globosa. Genome Biol. Evol. 8(11), 3292–3296 (2016).
- 16. Draft genome sequence of the dimorphic fungus Sporothrix pallida, a nonpathogenic species belonging to sporothrix, a genus containing agents of human and feline sporotrichosis. Genome Announc. 4(2), 1–2 (2016).
- 17. Identification of multiple species and subpopulations among Australian clinical Sporothrix isolates using whole genome sequencing. Med. Mycol. 57(7), 905–908 (2019).
- 18. Draft genome sequences of two Sporothrix schenckii clinical isolates associated with human sporotrichosis in colombia. Genome Announc. 6(24), 1–2 (2018).
- 19. . BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31(19), 3210–3212 (2015).
- 20. BUSCO applications from quality assessments to gene prediction and phylogenomics. Mol. Biol. Evol. 35(3), 543–548 (2018).
- 21. . Using RepeatMasker to identify repetitive elements in genomic sequences. Curr. Protoc. Bioinformatics 25(Chapter 4, Unit 4.10), 4.10.1–4.10.14 (2009).
- 22. . AUGUSTUS: a web server for gene finding in eukaryotes. Nucleic Acids Res. 32(Web Server issue), W309–312 (2004).
- 23. . Eukaryotic gene prediction using GeneMark.hmm-E and GeneMark-ES. Curr. Protoc. Bioinformatics 35(Chapter 4(1), Unit 4.6.1-10), 4.6.1–4.6.10 (2011).
- 24. . Gene finding in novel genomes. BMC Bioinformatics. 5, 59 (2004).
- 25. Automated eukaryotic gene structure annotation using EVidenceModeler and the Program to Assemble Spliced Alignments. Genome Biol. 9(1), R7 (2008).
- 26. . Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25(14), 1754–1760 (2009).
- 27. . Sambamba: fast processing of NGS alignment formats. Bioinformatics 31(12), 2032–2034 (2015).
- 28. . bioalcidae, samjs and vcffilterjs: object-oriented formatters and filters for bioinformatics files. Bioinformatics 34(7), 1224–1225 (2018).
- 29. Rfam 12.0: updates to the RNA families database. Nucleic Acids Res. 43(D1), D130–D137 (2015).
- 30. . Infernal 1.1: 100-fold faster RNA homology searches. Bioinformatics 29(22), 2933–2935 (2013).
- 31. . Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25(17), 3389–3402 (1997).
- 32. . MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 30(14), 3059–3066 (2002).
- 33. . trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 25(15), 1972–1973 (2009).
- 34. . Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol. Biol. Evol. 17(4), 540–552 (2000).
- 35. . IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 32(1), 268–274 (2015).
- 36. . MRBAYES: bayesian inference of phylogenetic trees. Bioinformatics 17(8), 754–755 (2001).
- 37. . EvolView, an online tool for visualizing, annotating and managing phylogenetic trees. Nucleic Acids Res. 40(W1), W569–W572 (2012).
- 38. . MUMmer4: a fast and versatile genome alignment system. PLoS Comput. Biol. 14(1), e1005944 (2018).
- 39. Circos: an information aesthetic for comparative genomics. Genome Res. 19(9), 1639–1645 (2009).
- 40. . OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 13(9), 2178–2189 (2003).
- 41. The Pfam protein families database: towards a more sustainable future. Nucleic Acids Res. 44(D1), D279–D285 (2016).
- 42. . CAFE: a computational tool for the study of gene family evolution. Bioinformatics 22(10), 1269–1271 (2006).
- 43. . edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26(1), 139–140 (2010).
- 44. . PHI-base: a new database for pathogen host interactions. Nucleic Acids Res. 34(90001), D459–D464 (2006).
- 45. . The pathogen-host interactions database (PHI-base): additions and future developments. Nucleic Acids Res. 43(D1), D645–D655 (2015).
- 46. . FungalRV: adhesin prediction and immunoinformatics portal for human fungal pathogens. BMC Genomics 12(1), 192 (2011).
- 47. UniProtKB/Swiss-Prot, the manually annotated section of the UniProt KnowledgeBase: how to use the entry view. Methods Mol Biology. springer, 1374, 23–54 (2016).
- 48. . MEROPS: the peptidase database. Nucleic Acids Res. 27(1), 325–331 (1999).
- 49. . Twenty years of the MEROPS database of proteolytic enzymes, their substrates and inhibitors. Nucleic Acids Res. 44(D1), D343–D350 (2016).
- 50. . Challenges in homology search: HMMER3 and convergent evolution of coiled-coil regions. Nucleic Acids Res. 41(12), (2013).
- 51. . Fast and sensitive protein alignment using DIAMOND. Nat. Methods. 12(1), 59–60 (2014).
- 52. . Homology to peptide pattern for annotation of carbohydrate-active enzymes and prediction of function. BMC Bioinformatics 18(1), 1–9 (2017).
- 53. dbCAN2: a meta server for automated carbohydrate-active enzyme annotation. Nucleic Acids Res. 46(W1), W95–W101 (2018).
- 54. Apophysomyces variabilis: draft genome sequence and comparison of predictive virulence determinants with other medically important Mucorales. BMC Genomics 18(1), 736 (2017).
- 55. . Comparative Analysis of Secretomes from Ectomycorrhizal Fungi with an Emphasis on Small-Secreted Proteins. Front. Microbiol. 6(NOV), 1–15 (2015).
- 56. . The Transporter Classification Database (TCDB): recent advances. Nucleic Acids Res. 44(D1), D372–D379 (2016).
- 57. . Molecular mechanisms of iron uptake in fungi. Mol. Microbiol. 47(5), 1185–1197 (2003).
- 58. . An eight-cysteine-containing CFEM domain unique to a group of fungal membrane proteins. Trends Biochem. Sci. 28(3), 118–121 (2003).
- 59. Systematic analyses reveal uniqueness and origin of the CFEM domain in fungi. Sci. Rep. 5(1), 13032 (2015).
- 60. . LysM effectors: secreted proteins supporting fungal life. PLoS Pathog. 9(12), e1003769 (2013).
- 61. . Fungal LysM effectors: extinguishers of host immunity? Trends Microbiol. 17(4), 151–157 (2009).
- 62. Genomic analysis of the basal lineage fungus Rhizopus oryzae reveals a whole-genome duplication. PLoS Genet. 5(7), e1000549 (2009).
- 63. . Comparative phylogenomics of pathogenic and nonpathogenic species. G3; Genes|Genomes|Genetics 6(2), 235–244 (2016).
- 64. . Fungal dimorphism: the switch from hyphae to yeast is a specialized morphogenetic adaptation allowing colonization of a host. FEMS Microbiol. Rev. 39(6), 797–811 (2015).
- 65. . Two-component histidine kinase DRK1 is required for pathogenesis in Sporothrix schenckii. Mol. Med. Rep. 17(1), 721–728 (2018).
- 66. . Extracellular superoxide dismutase protects Histoplasma yeast cells from host-derived oxidative stress. PLoS Pathog. 8(5), e1002713 (2012).
- 67. . Heat shock proteins in Histoplasma and Paracoccidioides. Clin. Vaccine Immunol. 24(11), 1–8 (2017).
- 68. . Three DUF1996 proteins localize in vacuoles and function in fungal responses to multiple stresses and metal ions. Sci. Rep. 6(1), 20566 (2016).
- 69. . Histone methylation by SET domain proteins in fungi. Annu. Rev. Microbiol. 71(1), 413–439 (2017).
- 70. . SET domains and stress: uncovering new functions for yeast Set4. Curr. Genet. 65(3), 643–648 (2019).
- 71. . Histone H3K4 methylation regulates hyphal growth, secondary metabolism and multiple stress responses in Fusarium graminearum. Environ. Microbiol. 17(11), 4615–4630 (2015).
- 72. Histone H3 lysine 9 methyltransferase FvDim5 regulates fungal development, pathogenicity and osmotic stress responses in Fusarium verticillioides. FEMS Microbiol. Lett. 364(19), 1–8 (2017).
- 73. . Methylation of histone H3 lysine 36 is required for normal development in Neurospora crassa. Eukaryot. Cell. 4(8), 1455–1464 (2005).
- 74. The high affinity iron permease is a key virulence factor required for Rhizopus oryzae pathogenesis. Mol. Microbiol. 77(3), 587–604 (2010).
- 75. Components of a new gene family of ferroxidases involved in virulence are functionally specialized in fungal dimorphism. Sci. Rep. 8(1), 7660 (2018).
- 76. . Iron acquisition in fungal pathogens of humans. Metallomics 9(3), 215–227 (2017).
- 77. . The siderophore iron transporter of Candida albicans (Sit1p/Arn1p) mediates uptake of ferrichrome-type siderophores and is required for epithelial invasion. Infect. Immun. 70(9), 5246–5255 (2002).
- 78. . Host iron withholding demands siderophore utilization for Candida glabrata to survive macrophage killing. PLoS Pathog. 7(3), e1001322 (2011).
- 79. . Iron uptake and virulence in Histoplasma capsulatum. Curr. Opin. Microbiol. 16(6), 700–707 (2013).
- 80. . Adhesins in human fungal pathogens: glue with plenty of stick. Eukaryot. Cell. 12(4), 470–481 (2013).
- 81. . The cell wall-associated glyceraldehyde-3-phosphate dehydrogenase of Candida albicans is also a fibronectin and laminin binding protein. Infect. Immun. 66(5), 2052–2059 (1998).
- 82. Glyceraldehyde-3-phosphate dehydrogenase of Paracoccidioides brasiliensis is a cell surface protein involved in fungal adhesion to extracellular matrix proteins and interaction with cells. Infect. Immun. 74(1), 382–389 (2006).
- 83. . Surface-expressed enolase contributes to the adhesion of Paracoccidioides brasiliensis to host cells. FEMS Yeast Res. 12(5), 557–570 (2012).
- 84. . Biofilm formation by Candida albicans mutants for genes coding fungal proteins exhibiting the eight-cysteine-containing CFEM domain. FEMS Yeast Res. 6(7), 1074–1084 (2006).
- 85. . A systematic analysis reveals an essential role for high-affinity iron uptake system, haemolysin and CFEM domain-containing protein in iron homoeostasis and virulence in Candida glabrata. Biochem. J. 463(1), 103–114 (2014).
- 86. BcCFEM1, a CFEM domain-containing protein with putative GPI-anchored site, is involved in pathogenicity, conidial production and stress tolerance in Botrytis cinerea. Front. Microbiol. 8(SEP), 1–11 (2017).
- 87. Hemoglobin uptake by Paracoccidioides spp. is receptor-mediated. PLoS Negl. Trop. Dis. 8(5), e2856 (2014).
- 88. . Yeast transcriptome and in vivo hypoxia detection reveals Histoplasma capsulatum response to low oxygen tension. Med. Mycol. 54(1), myv073 (2015).
- 89. . Molecular diversity of LysM carbohydrate-binding motifs in fungi. Curr. Genet. 61(2), 103–113 (2015).
- 90. . Divergent LysM effectors contribute to the virulence of Beauveria bassiana by evasion of insect immune defenses. PLoS Pathog. 13(9), e1006604 (2017).
- 91. . Characterisation of novel-cell-wall LysM-domain proteins LdpA and LdpB from the human pathogenic fungus Aspergillus fumigatus. Sci. Rep. 9(1), 3345 (2019).
- 92. . Candida albicans secreted aspartyl proteinases in virulence and pathogenesis. Microbiol. Mol. Biol. Rev. 67(3), 400–28 table of contents (2003).
- 93. . Secreted proteases from dermatophytes. Mycopathologia 166(5–6), 285–294 (2008).
- 94. . Fungal lifestyle reflected in serine protease repertoire. Sci. Rep. 7(1), 9147 (2017).
- 95. . Molecular evolution of the deuterolysin (M35) family genes in Coccidioides. PLoS ONE 7(2), e31536 (2012).
- 96. . Lipolytic enzymes involved in the virulence of human pathogenic fungi. Mycobiology 41(2), 67–72 (2013).
- 97. . Genome sequencing and carbohydrate-active enzyme (CAZyme) repertoire of the white rot fungus Flammulina elastica. Int. J. Mol. Sci. 19(8), 2379 (2018).
- 98. . The fungal cell wall: structure, biosynthesis and function. Microbiol. Spectr. 5(3), 3341–3354 (2017).
- 99. . Immune recognition of fungal polysaccharides. J. Fungi. 3(3), 47 (2017).
- 100. . Histoplasma capsulatum -(1,3)-glucan blocks innate immune recognition by the beta-glucan receptor. Proc. Natl Acad. Sci. USA 104(4), 1366–1370 (2007).
- 101. . RNA interference in Histoplasma capsulatum demonstrates a role for α-(1,3)-glucan in virulence. Mol. Microbiol. 53(1), 153–165 (2004).
- 102. Exploring the genomic diversity of black yeasts and relatives (Chaetothyriales, Ascomycota). Stud. Mycol. 86, 1–28 (2017).
- 103. . Essential letters in the fungal alphabet: ABC and MFS transporters and their roles in survival and pathogenicity. Adv. Genet. 85, 201–253 (2014).
- 104. . CgOpt1, a putative oligopeptide transporter from Colletotrichum gloeosporioides that is involved in responses to auxin and pathogenicity. BMC Microbiol. 9(1), 173 (2009).
- 105. . Fusarium graminearum Tri12p influences virulence to wheat and trichothecene accumulation. Mol. Plant Microbe Interact. 25(11), 1408–1418 (2012).