Ekaterina Shelest

antiSMASH 4.0-improvements in chemistry prediction and gene cluster boundary identification

Abstract Many antibiotics, chemotherapeutics, crop protection agents and food preservatives originate from molecules produced by bacteria, fungi or plants. In recent years, genome mining methodologies have been widely adopted to identify and characterize the biosynthetic gene clusters encoding the production of such compounds. Since 2011, the ‘antibiotics and secondary metabolite analysis shell—antiSMASH’ has assisted researchers in efficiently performing this, both as a web server and a standalone tool. Here, we present the thoroughly updated antiSMASH version 4, which adds several novel features, including prediction of gene cluster boundaries using the ClusterFinder method or the newly integrated CASSIS algorithm, improved substrate specificity prediction for non-ribosomal peptide synthetase adenylation domains based on the new SANDPUMA algorithm, improved predictions for terpene and ribosomally synthesized and post-translationally modified peptides cluster products, reporting of sequence similarity to proteins encoded in experimentally characterized gene clusters on a per-protein basis and a domain-level alignment tool for comparative analysis of trans-AT polyketide synthase assembly line architectures. Additionally, several usability features have been updated and improved. Together, these improvements make antiSMASH up-to-date with the latest developments in natural product research and will further facilitate computational genome mining for the discovery of novel bioactive molecules.

Comparative and functional genomics provide insights into the pathogenicity of dermatophytic fungi.

BackgroundMillions of humans and animals suffer from superficial infections caused by a group of highly specialized filamentous fungi, the dermatophytes, which exclusively infect keratinized host structures. To provide broad insights into the molecular basis of the pathogenicity-associated traits, we report the first genome sequences of two closely phylogenetically related dermatophytes, Arthroderma benhamiae and Trichophyton verrucosum, both of which induce highly inflammatory infections in humans.Results97% of the 22.5 megabase genome sequences of A. benhamiae and T. verrucosum are unambiguously alignable and collinear. To unravel dermatophyte-specific virulence-associated traits, we compared sets of potentially pathogenicity-associated proteins, such as secreted proteases and enzymes involved in secondary metabolite production, with those of closely related onygenales (Coccidioides species) and the mould Aspergillus fumigatus. The comparisons revealed expansion of several gene families in dermatophytes and disclosed the peculiarities of the dermatophyte secondary metabolite gene sets. Secretion of proteases and other hydrolytic enzymes by A. benhamiae was proven experimentally by a global secretome analysis during keratin degradation. Molecular insights into the interaction of A. benhamiae with human keratinocytes were obtained for the first time by global transcriptome profiling. Given that A. benhamiae is able to undergo mating, a detailed comparison of the genomes further unraveled the genetic basis of sexual reproduction in this species.ConclusionsOur results enlighten the genetic basis of fundamental and putatively virulence-related traits of dermatophytes, advancing future research on these medically important pathogens.

Comparative Genome Analysis of Trichophyton rubrum and Related Dermatophytes Reveals Candidate Genes Involved in Infection

ABSTRACT The major cause of athlete’s foot is Trichophyton rubrum, a dermatophyte or fungal pathogen of human skin. To facilitate molecular analyses of the dermatophytes, we sequenced T. rubrum and four related species, Trichophyton tonsurans, Trichophyton equinum, Microsporum canis, and Microsporum gypseum. These species differ in host range, mating, and disease progression. The dermatophyte genomes are highly colinear yet contain gene family expansions not found in other human-associated fungi. Dermatophyte genomes are enriched for gene families containing the LysM domain, which binds chitin and potentially related carbohydrates. These LysM domains differ in sequence from those in other species in regions of the peptide that could affect substrate binding. The dermatophytes also encode novel sets of fungus-specific kinases with unknown specificity, including nonfunctional pseudokinases, which may inhibit phosphorylation by competing for kinase sites within substrates, acting as allosteric effectors, or acting as scaffolds for signaling. The dermatophytes are also enriched for a large number of enzymes that synthesize secondary metabolites, including dermatophyte-specific genes that could synthesize novel compounds. Finally, dermatophytes are enriched in several classes of proteases that are necessary for fungal growth and nutrient acquisition on keratinized tissues. Despite differences in mating ability, genes involved in mating and meiosis are conserved across species, suggesting the possibility of cryptic mating in species where it has not been previously detected. These genome analyses identify gene families that are important to our understanding of how dermatophytes cause chronic infections, how they interact with epithelial cells, and how they respond to the host immune response. IMPORTANCE Athlete’s foot, jock itch, ringworm, and nail infections are common fungal infections, all caused by fungi known as dermatophytes (fungi that infect skin). This report presents the genome sequences of Trichophyton rubrum, the most frequent cause of athlete’s foot, as well as four other common dermatophytes. Dermatophyte genomes are enriched for four gene classes that may contribute to the ability of these fungi to cause disease. These include (i) proteases secreted to degrade skin; (ii) kinases, including pseudokinases, that are involved in signaling necessary for adapting to skin; (iii) secondary metabolites, compounds that act as toxins or signals in the interactions between fungus and host; and (iv) a class of proteins (LysM) that appear to bind and mask cell wall components and carbohydrates, thus avoiding the host’s immune response to the fungi. These genome sequences provide a strong foundation for future work in understanding how dermatophytes cause disease. Athlete’s foot, jock itch, ringworm, and nail infections are common fungal infections, all caused by fungi known as dermatophytes (fungi that infect skin). This report presents the genome sequences of Trichophyton rubrum, the most frequent cause of athlete’s foot, as well as four other common dermatophytes. Dermatophyte genomes are enriched for four gene classes that may contribute to the ability of these fungi to cause disease. These include (i) proteases secreted to degrade skin; (ii) kinases, including pseudokinases, that are involved in signaling necessary for adapting to skin; (iii) secondary metabolites, compounds that act as toxins or signals in the interactions between fungus and host; and (iv) a class of proteins (LysM) that appear to bind and mask cell wall components and carbohydrates, thus avoiding the host’s immune response to the fungi. These genome sequences provide a strong foundation for future work in understanding how dermatophytes cause disease.

Activation of a silent fungal polyketide biosynthesis pathway through regulatory cross talk with a cryptic nonribosomal peptide synthetase gene cluster

ABSTRACT Filamentous fungi produce numerous natural products that constitute a consistent source of potential drug leads, yet it seems that the majority of natural products are overlooked since most biosynthesis gene clusters are silent under standard cultivation conditions. Screening secondary metabolite genes of the model fungus Aspergillus nidulans, we noted a silent gene cluster on chromosome II comprising two nonribosomal peptide synthetase (NRPS) genes, inpA and inpB, flanked by a regulatory gene that we named scpR for secondary metabolism cross-pathway regulator. The induced expression of the scpR gene using the promoter of the alcohol dehydrogenase AlcA led to the transcriptional activation of both the endogenous scpR gene and the NRPS genes. Surprisingly, metabolic profiling of the supernatant of mycelia overexpressing scpR revealed the production of the polyketide asperfuranone. Through transcriptome analysis we found that another silent secondary metabolite gene cluster located on chromosome VIII coding for asperfuranone biosynthesis was specifically induced. Quantitative reverse transcription-PCR proved the transcription not only of the corresponding polyketide synthase (PKS) biosynthesis genes, afoE and afoG, but also of their activator, afoA, under alcAp-scpR-inducing conditions. To exclude the possibility that the product of the inp cluster induced the asperfuranone gene cluster, a strain carrying a deletion of the NRPS gene inpB and, in addition, the alcAp-scpR overexpression cassette was generated. In this strain, under inducing conditions, transcripts of the biosynthesis genes of both the NRPS-containing gene cluster inp and the asperfuranone gene cluster except gene inpB were detected. Moreover, the existence of the polyketide product asperfuranone indicates that the transcription factor ScpR controls the expression of the asperfuranone biosynthesis gene cluster. This expression as well as the biosynthesis of asperfuranone was abolished after the deletion of the asperfuranone activator gene afoA, indicating that ScpR binds to the afoA promoter. To the best of our knowledge, this is the first report of regulatory cross talk between two biosynthesis gene clusters located on different chromosomes.

Novel traits of Trichoderma predicted through the analysis of its secretome

Abstract Mycotrophic species of Trichoderma are among the most common fungi isolated from free soil, dead wood and as parasites on sporocarps of other fungi (mycoparasites). In addition, they undergo various other biotrophic associations ranging from rhizosphere colonization and endophytism up to facultative pathogenesis on such animals as roundworms and humans. Together with occurrence on a variety of less common substrata (marine invertebrates, artificial materials, indoor habitats), these lifestyles illustrate a wealthy opportunistic potential of the fungus. One tropical species, Trichoderma reesei, has become a prominent producer of cellulases and hemicellulases, whereas several other species are applied in agriculture for the biological control of phytopathogenic fungi. The sequencing of the complete genomes of the three species (T. reesei, T. virens, and T. atroviride) has led to a deepened understanding of Trichoderma lifestyle and its molecular physiology. In this review, we present the in silico predicted secretome of Trichoderma, and – in addition to the unique features of carbohydrate active enzymes – demonstrate the importance of such protein families as proteases, oxidative enzymes, and small cysteine-rich proteins, all of that received little attention in Trichoderma genetics so far. We also discuss the link between Trichoderma secretome and biology of the fungus.

Breaking the Silence: Protein Stabilization Uncovers Silenced Biosynthetic Gene Clusters in the Fungus Aspergillus nidulans

ABSTRACT The genomes of filamentous fungi comprise numerous putative gene clusters coding for the biosynthesis of chemically and structurally diverse secondary metabolites (SMs), which are rarely expressed under laboratory conditions. Previous approaches to activate these genes were based primarily on artificially targeting the cellular protein synthesis apparatus. Here, we applied an alternative approach of genetically impairing the protein degradation apparatus of the model fungus Aspergillus nidulans by deleting the conserved eukaryotic csnE/CSN5 deneddylase subunit of the COP9 signalosome. This defect in protein degradation results in the activation of a previously silenced gene cluster comprising a polyketide synthase gene producing the antibiotic 2,4-dihydroxy-3-methyl-6-(2-oxopropyl)benzaldehyde (DHMBA). The csnE/CSN5 gene is highly conserved in fungi, and therefore, the deletion is a feasible approach for the identification of new SMs.

CASSIS and SMIPS: promoter-based prediction of secondary metabolite gene clusters in eukaryotic genomes

Motivation: Secondary metabolites (SM) are structurally diverse natural products of high pharmaceutical importance. Genes involved in their biosynthesis are often organized in clusters, i.e., are co-localized and co-expressed. In silico cluster prediction in eukaryotic genomes remains problematic mainly due to the high variability of the clusters’ content and lack of other distinguishing sequence features. Results: We present Cluster Assignment by Islands of Sites (CASSIS), a method for SM cluster prediction in eukaryotic genomes, and Secondary Metabolites by InterProScan (SMIPS), a tool for genome-wide detection of SM key enzymes (‘anchor’ genes): polyketide synthases, non-ribosomal peptide synthetases and dimethylallyl tryptophan synthases. Unlike other tools based on protein similarity, CASSIS exploits the idea of co-regulation of the cluster genes, which assumes the existence of common regulatory patterns in the cluster promoters. The method searches for ‘islands’ of enriched cluster-specific motifs in the vicinity of anchor genes. It was validated in a series of cross-validation experiments and showed high sensitivity and specificity. Availability and implementation: CASSIS and SMIPS are freely available at https://sbi.hki-jena.de/cassis. Contact: thomas.wolf@leibniz-hki.de or ekaterina.shelest@leibniz-hki.de Supplementary information: Supplementary data are available at Bioinformatics online.

Gene expansion shapes genome architecture in the human pathogen Lichtheimia corymbifera: an evolutionary genomics analysis in the ancient terrestrial mucorales (Mucoromycotina)

Lichtheimia species are the second most important cause of mucormycosis in Europe. To provide broader insights into the molecular basis of the pathogenicity-associated traits of the basal Mucorales, we report the full genome sequence of L. corymbifera and compared it to the genome of Rhizopus oryzae, the most common cause of mucormycosis worldwide. The genome assembly encompasses 33.6 MB and 12,379 protein-coding genes. This study reveals four major differences of the L. corymbifera genome to R. oryzae: (i) the presence of an highly elevated number of gene duplications which are unlike R. oryzae not due to whole genome duplication (WGD), (ii) despite the relatively high incidence of introns, alternative splicing (AS) is not frequently observed for the generation of paralogs and in response to stress, (iii) the content of repetitive elements is strikingly low (<5%), (iv) L. corymbifera is typically haploid. Novel virulence factors were identified which may be involved in the regulation of the adaptation to iron-limitation, e.g. LCor01340.1 encoding a putative siderophore transporter and LCor00410.1 involved in the siderophore metabolism. Genes encoding the transcription factors LCor08192.1 and LCor01236.1, which are similar to GATA type regulators and to calcineurin regulated CRZ1, respectively, indicating an involvement of the calcineurin pathway in the adaption to iron limitation. Genes encoding MADS-box transcription factors are elevated up to 11 copies compared to the 1–4 copies usually found in other fungi. More findings are: (i) lower content of tRNAs, but unique codons in L. corymbifera, (ii) Over 25% of the proteins are apparently specific for L. corymbifera. (iii) L. corymbifera contains only 2/3 of the proteases (known to be essential virulence factors) in comparision to R. oryzae. On the other hand, the number of secreted proteases, however, is roughly twice as high as in R. oryzae.

Identification of the novel penicillin biosynthesis gene aatB of Aspergillus nidulans and its putative evolutionary relationship to this fungal secondary metabolism gene cluster

The final step of penicillin biosynthesis in the filamentous fungus Aspergillus nidulans is catalysed by isopenicillin N acyltransferase encoded by the aatA gene. Because there is no bacterial homologue, its evolutionary origin remained obscure. As shown here, disruption of aatA still enabled penicillin production. Genome mining led to the discovery of the aatB gene (AN6775.3) which has a similar structure and expression pattern as aatA. Disruption of aatB resulted in a reduced penicillin titre. Surface plasmon resonance analysis and Northern blot analysis indicated that the promoters of both aatA and aatB are bound and regulated by the same transcription factors AnCF and AnBH1f. In contrast to aatA, aatB does not encode a peroxisomal targeting signal (PTS1). Overexpression of a mutated aatBPTS1 gene in an aatA‐disruption strain (leading to peroxisomal localization of AatB) increased the penicillin titre more than overexpression of the wild‐type aatB. Homologues of aatA are exclusively part of the penicillin biosynthesis gene cluster, whereas aatB homologues also exist in non‐producing fungi. Our findings suggest that aatB is a paralogue of aatA. They extend the model of evolution of the penicillin biosynthesis gene cluster by recruitment of a biosynthesis gene and its cis‐regulatory sites upon gene duplication.

The fungal α-aminoadipate pathway for lysine biosynthesis requires two enzymes of the aconitase family for the isomerization of homocitrate to homoisocitrate

Fungi produce α‐aminoadipate, a precursor for penicillin and lysine via the α‐aminoadipate pathway. Despite the biotechnological importance of this pathway, the essential isomerization of homocitrate via homoaconitate to homoisocitrate has hardly been studied. Therefore, we analysed the role of homoaconitases and aconitases in this isomerization. Although we confirmed an essential contribution of homoaconitases from Saccharomyces cerevisiae and Aspergillus fumigatus, these enzymes only catalysed the interconversion between homoaconitate and homoisocitrate. In contrast, aconitases from fungi and the thermophilic bacterium Thermus thermophilus converted homocitrate to homoaconitate. Additionally, a single aconitase appears essential for energy metabolism, glutamate and lysine biosynthesis in respirating filamentous fungi, but not in the fermenting yeast S. cerevisiae that possesses two contributing aconitases. While yeast Aco1p is essential for the citric acid cycle and, thus, for glutamate synthesis, Aco2p specifically and exclusively contributes to lysine biosynthesis. In contrast, Aco2p homologues present in filamentous fungi were transcribed, but enzymatically inactive, revealed no altered phenotype when deleted and did not complement yeast aconitase mutants. From these results we conclude that the essential requirement of filamentous fungi for respiration versus the preference of yeasts for fermentation may have directed the evolution of aconitases contributing to energy metabolism and lysine biosynthesis.