Independent subtilases expansions in fungi associated with animals

Signatures of transposon-mediated genome inflation, host specialization, and photoentrainment in Entomophthora muscae and allied entomophthoralean fungi

Despite over a century of observations, the obligate insect parasites within the order Entomophthorales remain poorly characterized at the genetic level. In this manuscript, we present a genome for a laboratory-tractable Entomophthora muscae isolate that infects fruit flies. Our E. muscae assembly is 1.03 Gb, consists of 7810 contigs and contains 81.3% complete fungal BUSCOs. Using a comparative approach with recent datasets from entomophthoralean fungi, we show that giant genomes are the norm within Entomophthoraceae owing to extensive, but not recent, Ty3 retrotransposon activity. In addition, we find that E. muscae and its closest allies possess genes that are likely homologs to the blue-light sensor white-collar 1, a Neurospora crassa gene that has a well-established role in maintaining circadian rhythms. We uncover evidence that E. muscae diverged from other entomophthoralean fungi by expansion of existing families, rather than loss of particular domains, and possesses a potentially unique suite of secreted catabolic enzymes, consistent with E. muscae's species-specific, biotrophic lifestyle. Finally, we offer a head-to-head comparison of morphological and molecular data for species within the E. muscae species complex that support the need for taxonomic revision within this group. Altogether, we provide a genetic and molecular foundation that we hope will provide a platform for the continued study of the unique biology of entomophthoralean fungi.

Signatures of transposon-mediated genome inflation, host specialization, and photoentrainment in Entomophthora muscae and allied entomophthoralean fungi

Despite over a century of observations, the obligate insect parasites within the order Entomophthorales remain poorly characterized at the genetic level. This is in part due to their large genome sizes and difficulty in obtaining sequenceable material. In this manuscript, we leveraged a recently-isolated, laboratory-tractable Entomophthora muscae isolate and improved long-read sequencing to obtain a largely-complete entomophthoralean genome. Our E. muscae assembly is 1.03 Gb, consists of 7,810 contigs and contains 81.3% complete fungal BUSCOs. Using a comparative approach with other available (transcriptomic and genomic) datasets from entomophthoralean fungi, we provide new insight into the biology of these understudied pathogens. We offer a head-to-head comparison of morphological and molecular data for species within the E. muscae species complex. Our findings suggest that substantial taxonomic revision is needed to define species within this group and we provide recommendations for differentiating strains and species in the context of the existing body of E. muscae scientific literature. We show that giant genomes are the norm within Entomophthoraceae owing to extensive, but not recent, Ty3 retrotransposon activity, despite the presence of machinery to defend against transposable elements(RNAi). In addition, we find that E. muscae and its closest allies are enriched for M16A peptidases and possess genes that are likely homologs to the blue-light sensor white-collar 1, a Neurospora crassa gene that has a well-established role in maintaining circadian rhythms. We find that E. muscae has an expanded group of acid-trehalases, consistent with trehalose being the primary sugar component of fly (and insect) hemolymph. We uncover evidence that E. muscae diverged from other entomophthoralean fungi by expansion of existing families, rather than loss of particular domains, and possesses a potentially unique suite of secreted catabolic enzymes, consistent with E. muscae's species-specific, biotrophic lifestyle. Altogether, we provide a genetic and molecular foundation that we hope will provide a platform for the continued study of the unique biology of entomophthoralean fungi.

Serine protease CrKP43 interacts with MAPK and regulates fungal development and mycoparasitism in Clonostachys chloroleuca

IMPORTANCE

Mycoparasites play important roles in the biocontrol of plant fungal diseases, during which they secret multiple hydrolases such as serine proteases to degrade their fungal hosts. In this study, we demonstrated that the serine protease CrKP43 was involved in C. chloroleuca development and mycoparasitism with the regulation of Crmapk. To the best of our knowledge, it is the first report on the functions and regulatory mechanisms of serine proteases in C. chloroleuca. Our findings will provide new insight into the regulatory mechanisms of serine proteases in mycoparasites and contribute to clarifying the mechanisms underlying mycoparasitism of C. chloroleuca, which will facilitate the development of highly efficient fungal biocontrol agents as well.

Biochemical and molecular characterization of novel keratinolytic protease from Bacillus licheniformis (KRLr1)

The keratin-degrading bacterium Bacillus licheniformis secretes a keratinase with potential industrial interest. Here, the Keratinase gene was intracellularly expressed in Escherichia coli BL21(DE3) using pET-21b (+) vector. Phylogenetic tree analysis showed that KRLr1 is closely related to Bacillus licheniformis keratinase that belongs to the serine peptidase/subtilisin-like S8 family. Recombinant keratinase appeared on the SDS-PAGE gel with a band of about 38 kDa and was confirmed by western blotting. Expressed KRLr1 was purified by Ni-NTA affinity chromatography with a yield of 85.96% and then refolded. It was found that this enzyme has optimum activity at pH 6 and 37°C. PMSF inhibited the KRLr1 activity and Ca²⁺ and Mg²⁺ increased the KRLr1 activity. Using keratin 1% as the substrate, the thermodynamic values were determined as K_m 14.54 mM, k_cat 912.7 × 10^-3 (S^-1), and k_cat/K_m 62.77 (M^-1 S^-1). Feather digestion by recombinant enzyme using HPLC method, showed that the amino acids cysteine, phenylalanine, tyrosine and lysine had the highest amount compared to other amino acids obtained from digestion. Molecular dynamics (MD) simulation of HADDOCK docking results exhibited that KRLr1 enzyme was able to interact strongly with chicken feather keratine 4 (FK4) compared to chicken feather keratine 12 (FK12). These properties make keratinase KRLr1 a potential candidate for various biotechnological applications.

Mycoparasites, Gut Dwellers, and Saprotrophs: Phylogenomic Reconstructions and Comparative Analyses of Kickxellomycotina Fungi

Improved sequencing technologies have profoundly altered global views of fungal diversity and evolution. High-throughput sequencing methods are critical for studying fungi due to the cryptic, symbiotic nature of many species, particularly those that are difficult to culture. However, the low coverage genome sequencing (LCGS) approach to phylogenomic inference has not been widely applied to fungi. Here we analyzed 171 Kickxellomycotina fungi using LCGS methods to obtain hundreds of marker genes for robust phylogenomic reconstruction. Additionally, we mined our LCGS data for a set of nine rDNA and protein coding genes to enable analyses across species for which no LCGS data were obtained. The main goals of this study were to: 1) evaluate the quality and utility of LCGS data for both phylogenetic reconstruction and functional annotation, 2) test relationships among clades of Kickxellomycotina, and 3) perform comparative functional analyses between clades to gain insight into putative trophic modes. In opposition to previous studies, our nine-gene analyses support two clades of arthropod gut dwelling species and suggest a possible single evolutionary event leading to this symbiotic lifestyle. Furthermore, we resolve the mycoparasitic Dimargaritales as the earliest diverging clade in the subphylum and find four major clades of Coemansia species. Finally, functional analyses illustrate clear variation in predicted carbohydrate active enzymes and secondary metabolites (SM) based on ecology, that is biotroph versus saprotroph. Saprotrophic Kickxellales broadly lack many known pectinase families compared with saprotrophic Mucoromycota and are depauperate for SM but have similar numbers of predicted chitinases as mycoparasitic.

Phylogenetic survey of the subtilase family and a data-mining-based search for new subtilisins from Bacillaceae

The subtilase family (S8), a member of the clan SB of serine proteases are ubiquitous in all kingdoms of life and fulfil different physiological functions. Subtilases are divided in several groups and especially subtilisins are of interest as they are used in various industrial sectors. Therefore, we searched for new subtilisin sequences of the family Bacillaceae using a data mining approach. The obtained 1,400 sequences were phylogenetically classified in the context of the subtilase family. This required an updated comprehensive overview of the different groups within this family. To fill this gap, we conducted a phylogenetic survey of the S8 family with characterised holotypes derived from the MEROPS database. The analysis revealed the presence of eight previously uncharacterised groups and 13 subgroups within the S8 family. The sequences that emerged from the data mining with the set filter parameters were mainly assigned to the subtilisin subgroups of true subtilisins, high-alkaline subtilisins, and phylogenetically intermediate subtilisins and represent an excellent source for new subtilisin candidates.

Evolutionary transition to the ectomycorrhizal habit in the genomes of a hyperdiverse lineage of mushroom-forming fungi

The ectomycorrhizal (ECM) symbiosis has independently evolved from diverse types of saprotrophic ancestors. In this study, we seek to identify genomic signatures of the transition to the ECM habit within the hyperdiverse Russulaceae. We present comparative analyses of the genomic architecture and the total and secreted gene repertoires of 18 species across the order Russulales, of which 13 are newly sequenced, including a representative of a saprotrophic member of Russulaceae, Gloeopeniophorella convolvens. The genomes of ECM Russulaceae are characterized by a loss of genes for plant cell wall-degrading enzymes (PCWDEs), an expansion of genome size through increased transposable element (TE) content, a reduction in secondary metabolism clusters, and an association of small secreted proteins (SSPs) with TE 'nests', or dense aggregations of TEs. Some PCWDEs have been retained or even expanded, mostly in a species-specific manner. The genome of G. convolvens possesses some characteristics of ECM genomes (e.g. loss of some PCWDEs, TE expansion, reduction in secondary metabolism clusters). Functional specialization in ECM decomposition may drive diversification. Accelerated gene evolution predates the evolution of the ECM habit, indicating that changes in genome architecture and gene content may be necessary to prime the evolutionary switch.

Comparative Genomics and Molecular Analysis of Epidermophyton floccosum

Epidermophyton floccosum is one of the most common agents of human superficial fungal infections, compared with genus Trichophyton and Microsporum, it possesses uniqueness in ecology traits and rarely causing hair infections. E. floccosum is so far the only representative species of genera Epidermophyton, and it is known as anthropophilic dermatophytes. To further reveal the genome sequences and clues of virulence factors, thus in this study, we sequenced the genome of E. floccosum (CGMCC (F) E1d), and performed comparative genomic analysis with other dermatophytes. It is revealed that E. floccosum owns the largest genome size and similar GC content compared with other dermatophytes. A total of 7565 genes are predicted. By comparing with the closest species N. gypseum, our study reveals that number and structure of adhesion factors, secreted proteases and LysM domain might contribute to the pathogenic and ecological traits of E. floccosum. Mating genes is also detected in genome data. Furthermore, we performed AFLP analysis trying to discuss intraspecific differences of E. floccosum, but no significant relationship is found between genotype and geographical distribution. Upon above, our study provides a deeper understanding and strong foundation for future researches about E. floccosum.

Interaction between the cell walls of microalgal host and fungal carbohydrate-activate enzymes is essential for the pathogenic parasitism process

Fungi can parasitize microalgae, exerting profound impacts on both the aquatic ecosystems and microalgal mass cultures. In this study, the unicellular green alga Haematococcus pluvialis and the blastocladialean fungus Paraphysoderma sedebokerense were used as a model system to address the mechanisms underlying the fungal parasitism on the algal host. High-throughput metabolic assay indicated that P. sedebokerense can utilize several carbon sources with a preference for mannose, glucose and their oligosaccharides, which was compatible with the profile of the host algal cell walls enriched with glucan and mannan. The results of dual transcriptomics analysis suggested that P. sedebokerense can upregulate a large number of putative carbohydrate-activate enzymes (CAZymes) encoding genes, including those coding for the endo-1,4-β-glucanase and endo-1,4-β-mannanase during the infection process. The cell walls of H. pluvialis can be decomposed by both P. sedebokerense and commercial CAZymes (e.g. cellulase and endo-1,4-β-mannanase) to produce mannooligomers, while several putative parasitism-related genes of P. sedebokerense can be in turn upregulated by mannooligomers. In addition, the parasitism can be blocked by interfering the selected CAZymes including glucanase, mannanase and lysozyme with the specific inhibitors, which provided a framework for screening suitable compounds for pathogen mitigation in algal mass culture.

System-wide characterization of subtilases reveals that subtilisin-like protease FgPrb1 of Fusarium graminearum regulates fungal development and virulence

Subtilases represent the second largest subfamily of serine proteases, and are important for various biological processes. However, the biological function of subtilases has not been systematically characterized in plant pathogens. In present study, 32 subtilases were identified in the genome of wheat scab fungus Fusarium graminearum, a devastating cereal plant pathogen. Deletion mutants of each subtilase were obtained and functionally characterized. Among them, the deletion of FgPrb1 resulted in greatly reduced virulence of F. graminearum. The regulatory mechanisms of FgPrb1 in virulence were investigated in details. Our results showed that the loss of FgPrb1 led to defects in deoxynivalenol (DON) production, responses to environmental stimuli, and lipid metabolism. Additionally, we found that FgPrb1 was involved in autophagy regulation. Taken together, the systematic functional characterization of subtilases showed that the FgPrb1 of F. graminearum is critical for plant infection by regulating multiple different cellular processes.

Genome sequence of Isaria javanica and comparative genome analysis insights into family S53 peptidase evolution in fungal entomopathogens

The fungus Isaria javanica is an important entomopathogen that parasitizes various insects and is effective for pest control. In this study, we sequenced and assembled the genomes (IJ1G and IJ2G) of two I. javanica strains isolated from different insects. The genomes were approximately 35 Mb in size with 11,441 and 11,143 protein-coding genes, respectively. Using a phylogenomic approach, we evaluated genome evolution across five entomopathogenic fungi in Cordycipitaceae. By comparative genome analysis, it was found that family S53 serine peptidases were expanded in Cordycipitaceae entomopathogens, particularly in I. javanica. Gene duplication events were identified based on phylogenetic relationships inferred from 82 S53 peptidases within six entomopathogenic fungal genomes. Moreover, we found that carbohydrate-active enzymes and proteinases were the largest secretory protein groups encoded in the I. javanica genome, especially chitinases (GH18), serine and aspartic peptidases (S53, S08, S10, A01). Pathogenesis-related genes and genes for bacterial-like toxins and secondary metabolites were also identified. By comparative transcriptome analysis, differentially expressed genes in response to insect nutrients (in vitro) were identified. Moreover, most S53 peptidases were detected to be significantly upregulated during the initial fungal infection process in insects (in vivo) by RT-qPCR. Our results provide new clues about understanding evolution of pathogenic proteases and may suggest that abundant S53 peptidases in the I. javanica genome may contribute to its effective parasitism on various insects.

Comparative evolutionary histories of fungal proteases reveal gene gains in the mycoparasitic and nematode-parasitic fungus Clonostachys rosea

Background

The ascomycete fungus Clonostachys rosea (order Hypocreales) can control several important plant diseases caused by plant pathogenic fungi and nematodes. Subtilisin-like serine proteases are considered to play an important role in pathogenesis in entomopathogenic, mycoparasitic, and nematophagous fungi used for biological control. In this study, we analysed the evolutionary histories of protease gene families, and investigated sequence divergence and regulation of serine protease genes in C. rosea.

Results

Proteases of selected hypocrealean fungal species were classified into families based on the MEROPS peptidase database. The highest number of protease genes (590) was found in Fusarium solani, followed by C. rosea with 576 genes. Analysis of gene family evolution identified non-random changes in gene copy numbers in the five serine protease gene families S1A, S8A, S9X, S12 and S33. Four families, S1A, S8A, S9X, and S33, displayed gene gains in C. rosea. A gene-tree / species-tree reconciliation analysis of the S8A family revealed that the gene copy number increase in C. rosea was primarily associated with the S08.054 (proteinase K) subgroup. In addition, regulatory and predicted structural differences, including twelve sites evolving under positive selection, among eighteen C. rosea S8A serine protease paralog genes were also observed. The C. rosea S8A serine protease gene prs6 was induced during interaction with the plant pathogenic species F. graminearum.

Conclusions

Non-random increases in S8A, S9X and S33 serine protease gene numbers in the mycoparasitic species C. rosea, Trichoderma atroviride and T. virens suggests an involvement in fungal-fungal interactions. Regulatory and predicted structural differences between C. rosea S8A paralogs indicate that functional diversification is driving the observed increase in gene copy numbers. The induction of prs6 expression in C. rosea during confrontation with F. graminearum suggests an involvement of the corresponding protease in fungal-fungal interactions. The results pinpoint the importance of serine proteases for ecological niche adaptation in C. rosea, including a potential role in the mycoparasitic attack on fungal prey.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1291-1) contains supplementary material, which is available to authorized users.

Leveraging single-cell genomics to expand the fungal tree of life

Environmental DNA surveys reveal that most fungal diversity represents uncultured species. We sequenced the genomes of eight uncultured species across the fungal tree of life using a new single-cell genomics pipeline. We show that, despite a large variation in genome and gene space recovery from each single amplified genome (SAG), ≥90% can be recovered by combining multiple SAGs. SAGs provide robust placement for early-diverging lineages and infer a diploid ancestor of fungi. Early-diverging fungi share metabolic deficiencies and show unique gene expansions correlated with parasitism and unculturability. Single-cell genomics holds great promise in exploring fungal diversity, life cycles and metabolic potential.

Early Diverging Insect-Pathogenic Fungi of the Order Entomophthorales Possess Diverse and Unique Subtilisin-Like Serine Proteases

Insect-pathogenic fungi use subtilisin-like serine proteases (SLSPs) to degrade chitin-associated proteins in the insect procuticle. Most insect-pathogenic fungi in the order Hypocreales (Ascomycota) are generalist species with a broad host-range, and most species possess a high number of SLSPs. The other major clade of insect-pathogenic fungi is part of the subphylum Entomophthoromycotina (Zoopagomycota, formerly Zygomycota) which consists of high host-specificity insect-pathogenic fungi that naturally only infect a single or very few host species. The extent to which insect-pathogenic fungi in the order Entomophthorales rely on SLSPs is unknown. Here we take advantage of recently available transcriptomic and genomic datasets from four genera within Entomophthoromycotina: the saprobic or opportunistic pathogens Basidiobolus meristosporus, Conidiobolus coronatus, C. thromboides, C. incongruus, and the host-specific insect pathogens Entomophthora muscae and Pandora formicae, specific pathogens of house flies (Muscae domestica) and wood ants (Formica polyctena), respectively. In total 154 SLSP from six fungi in the subphylum Entomophthoromycotina were identified: E. muscae (n = 22), P. formicae (n = 6), B. meristosporus (n = 60), C. thromboides (n = 18), C. coronatus (n = 36), and C. incongruus (n = 12). A unique group of 11 SLSPs was discovered in the genomes of the obligate biotrophic fungi E. muscae, P. formicae and the saprobic human pathogen C. incongruus that loosely resembles bacillopeptidase F-like SLSPs. Phylogenetics and protein domain analysis show this class represents a unique group of SLSPs so far only observed among Bacteria, Oomycetes and early diverging fungi such as Cryptomycota, Microsporidia, and Entomophthoromycotina. This group of SLSPs is missing in the sister fungal lineages of Kickxellomycotina and the fungal phyla Mucoromyocta, Ascomycota and Basidiomycota fungi suggesting interesting gene loss patterns.

Fungal Genomes and Insights into the Evolution of the Kingdom

The kingdom Fungi comprises species that inhabit nearly all ecosystems. Fungi exist as both free-living and symbiotic unicellular and multicellular organisms with diverse morphologies. The genomes of fungi encode genes that enable them to thrive in diverse environments, invade plant and animal cells, and participate in nutrient cycling in terrestrial and aquatic ecosystems. The continuously expanding databases of fungal genome sequences have been generated by individual and large-scale efforts such as Génolevures, Broad Institute's Fungal Genome Initiative, and the 1000 Fungal Genomes Project (http://1000.fungalgenomes.org). These efforts have produced a catalog of fungal genes and genomic organization. The genomic datasets can be utilized to better understand how fungi have adapted to their lifestyles and ecological niches. Large datasets of fungal genomic and transcriptomic data have enabled the use of novel methodologies and improved the study of fungal evolution from a molecular sequence perspective. Combined with microscopes, petri dishes, and woodland forays, genome sequencing supports bioinformatics and comparative genomics approaches as important tools in the study of the biology and evolution of fungi.

Proteases and protease inhibitors in infectious diseases

There are numerous proteases of pathogenic organisms that are currently targeted for therapeutic intervention along with many that are seen as potential drug targets. This review discusses the chemical and biological makeup of some key druggable proteases expressed by the five major classes of disease causing agents, namely bacteria, viruses, fungi, eukaryotes, and prions. While a few of these enzymes including HIV protease and HCV NS3-4A protease have been targeted to a clinically useful level, a number are yet to yield any clinical outcomes in terms of antimicrobial therapy. A significant aspect of this review discusses the chemical and pharmacological characteristics of inhibitors of the various proteases discussed. A total of 25 inhibitors have been considered potent and safe enough to be trialed in humans and are at different levels of clinical application. We assess the mechanism of action and clinical performance of the protease inhibitors against infectious agents with their developmental strategies and look to the next frontiers in the use of protease inhibitors as anti-infective agents.

A critical review on serine protease: Key immune manipulator and pathology mediator

Proteolytic activity is fundamental to survival, so it is not surprising that all living organisms have proteases, especially seine protease. This enzyme in its numerous isoforms and homologues, constitutes the quintessential offence and defence factors, in the form of surface proteins, secreted molecules, gut digestive enzymes, venom in specialised glands or plant latex, among other manifestations. Occurring as trypsin, chymotrypsin, elastase, collagenase, thrombin, subtilisin etc., it mediates a diverse array of functions, including pathological roles as inflammatory, coagulatory to haemorrhagic. This review emphasizes that despite the superficial differences in mechanisms, most health issues, be they infectious, allergic, metabolic, or neural have a common conduit. This enzyme, in its various glycosylated forms leads to signal misinterpretations, wreaking havoc. However, organisms are endowed with serine protease inhibitors which might restrain this ubiquitous yet deleterious enzyme. Hence, serine proteases-driven pathogenesis and antagonising role of inhibitors is the focal point of this critical review.

Fungal lifestyle reflected in serine protease repertoire

Fungi are able to switch between different lifestyles in order to adapt to environmental changes. Their ecological strategy is connected to their secretome as fungi obtain nutrients by secreting hydrolytic enzymes to their surrounding and acquiring the digested molecules. We focus on fungal serine proteases (SPs), the phylogenetic distribution of which is barely described so far. In order to collect a complete set of fungal proteases, we searched over 600 fungal proteomes. Obtained results suggest that serine proteases are more ubiquitous than expected. From 54 SP families described in MEROPS Peptidase Database, 21 are present in fungi. Interestingly, 14 of them are also present in Metazoa and Viridiplantae - this suggests that, except one (S64), all fungal SP families evolved before plants and fungi diverged. Most representatives of sequenced eukaryotic lineages encode a set of 13-16 SP families. The number of SPs from each family varies among the analysed taxa. The most abundant are S8 proteases. In order to verify hypotheses linking lifestyle and expansions of particular SP, we performed statistical analyses and revealed previously undescribed associations. Here, we present a comprehensive evolutionary history of fungal SP families in the context of fungal ecology and fungal tree of life.

A Three-Way Transcriptomic Interaction Study of a Biocontrol Agent (Clonostachys rosea), a Fungal Pathogen (Helminthosporium solani), and a Potato Host (Solanum tuberosum)

Helminthosporium solani causes silver scurf, which affects the quality of potato. The biocontrol agent Clonostachys rosea greatly limited the severity of silver scurf symptoms and amount of H. solani genomic DNA in laboratory experiments. Transcriptomic analysis during interaction showed that H. solani gene expression was highly reduced when coinoculated with the biocontrol agent C. rosea, whereas gene expression of C. rosea was clearly boosted as a response to the pathogen. The most notable upregulated C. rosea genes were those encoding proteins involved in cellular response to oxidative stress, proteases, G-protein signaling, and the methyltransferase LaeA. The most notable potato response to both fungi was downregulation of defense-related genes and mitogen-activated protein kinase kinase kinases. At a later stage, this shifted, and most potato defense genes were turned on, especially those involved in terpenoid biosynthesis when H. solani was present. Some biocontrol-activated defense-related genes in potato were upregulated during early interaction with C. rosea alone that were not triggered by H. solani alone. Our results indicate that the reductions of silver scurf using C. rosea are probably due to a combination of mechanisms, including mycoparasitism, biocontrol-activated stimulation of plant defense mechanisms, microbial competition for nutrients, space, and antibiosis.

Pathogenicity-associated protein domains: The fiercely-conserved evolutionary signatures

Proteins have highly conserved domains that determine their functionality. Out of the thousands of domains discovered so far across all living forms, some of the predominant clinically-relevant domains include IENR1, HNHc, HELICc, Pro-kuma_activ, Tryp_SPc, Lactamase_B, PbH1, ChtBD3, CBM49, acidPPc, G3P_acyltransf, RPOL8c, KbaA, HAMP, HisKA, Hr1, Dak2, APC2, Citrate_ly_lig, DALR, VKc, YARHG, WR1, PWI, ZnF_BED, TUDOR, MHC_II_beta, Integrin_B_tail, Excalibur, DISIN, Cadherin, ACTIN, PROF, Robl_LC7, MIT, Kelch, GAS2, B41, Cyclin_C, Connexin_CCC, OmpH, Bac_rhodopsin, AAA, Knot1, NH, Galanin, IB, Elicitin, ACTH, Cache_2, CHASE, AgrB, PRP, IGR, and Antimicrobial21. These domains are distributed in nucleases/helicases, proteases, esterases, lipases, glycosylase, GTPases, phosphatases, methyltransferases, acyltransferase, acetyltransferase, polymerase, kinase, ligase, synthetase, oxidoreductase, protease inhibitors, nucleic acid binding proteins, adhesion and immunity-related proteins, cytoskeletal component-manipulating proteins, lipid biosynthesis and metabolism proteins, membrane-associated proteins, hormone-like and signaling proteins, etc. These domains are ubiquitous stretches or folds of the proteins in pathogens and allergens. Pathogenesis alleviation efforts can benefit enormously if the characteristics of these domains are known. Hence, this review catalogs and discusses the role of such pivotal domains, suggesting hypotheses for better understanding of pathogenesis at molecular level.

•

Proteins have highly conserved regions or domains across pathogens and allergens.

•

Knowledge on these critical domains can facilitate our understanding of pathogenesis mechanisms.

•

Such immune manipulation-related domains include IENR1, HNHc, HELICc, ACTIN, PROF, Robl_LC7, OmpH etc.

•

These domains are presnt in enzyme, transcription regulators, adhesion proteins, and hormones.

•

This review discusses and hypothesizes on these domains.

Phylogenomic evolutionary surveys of subtilase superfamily genes in fungi

Subtilases belong to a superfamily of serine proteases which are ubiquitous in fungi and are suspected to have developed distinct functional properties to help fungi adapt to different ecological niches. In this study, we conducted a large-scale phylogenomic survey of subtilase protease genes in 83 whole genome sequenced fungal species in order to identify the evolutionary patterns and subsequent functional divergences of different subtilase families among the main lineages of the fungal kingdom. Our comparative genomic analyses of the subtilase superfamily indicated that extensive gene duplications, losses and functional diversifications have occurred in fungi, and that the four families of subtilase enzymes in fungi, including proteinase K-like, Pyrolisin, kexin and S53, have distinct evolutionary histories which may have facilitated the adaptation of fungi to a broad array of life strategies. Our study provides new insights into the evolution of the subtilase superfamily in fungi and expands our understanding of the evolution of fungi with different lifestyles.

Exploring the genomic diversity of black yeasts and relatives (Chaetothyriales, Ascomycota)

The order Chaetothyriales (Pezizomycotina, Ascomycetes) harbours obligatorily melanised fungi and includes numerous etiologic agents of chromoblastomycosis, phaeohyphomycosis and other diseases of vertebrate hosts. Diseases range from mild cutaneous to fatal cerebral or disseminated infections and affect humans and cold-blooded animals globally. In addition, Chaetothyriales comprise species with aquatic, rock-inhabiting, ant-associated, and mycoparasitic life-styles, as well as species that tolerate toxic compounds, suggesting a high degree of versatile extremotolerance. To understand their biology and divergent niche occupation, we sequenced and annotated a set of 23 genomes of main the human opportunists within the Chaetothyriales as well as related environmental species. Our analyses included fungi with diverse life-styles, namely opportunistic pathogens and closely related saprobes, to identify genomic adaptations related to pathogenesis. Furthermore, ecological preferences of Chaetothyriales were analysed, in conjuncture with the order-level phylogeny based on conserved ribosomal genes. General characteristics, phylogenomic relationships, transposable elements, sex-related genes, protein family evolution, genes related to protein degradation (MEROPS), carbohydrate-active enzymes (CAZymes), melanin synthesis and secondary metabolism were investigated and compared between species. Genome assemblies varied from 25.81 Mb (Capronia coronata) to 43.03 Mb (Cladophialophora immunda). The bantiana-clade contained the highest number of predicted genes (12 817 on average) as well as larger genomes. We found a low content of mobile elements, with DNA transposons from Tc1/Mariner superfamily being the most abundant across analysed species. Additionally, we identified a reduction of carbohydrate degrading enzymes, specifically many of the Glycosyl Hydrolase (GH) class, while most of the Pectin Lyase (PL) genes were lost in etiological agents of chromoblastomycosis and phaeohyphomycosis. An expansion was found in protein degrading peptidase enzyme families S12 (serine-type D-Ala-D-Ala carboxypeptidases) and M38 (isoaspartyl dipeptidases). Based on genomic information, a wide range of abilities of melanin biosynthesis was revealed; genes related to metabolically distinct DHN, DOPA and pyomelanin pathways were identified. The MAT (MAting Type) locus and other sex-related genes were recognized in all 23 black fungi. Members of the asexual genera Fonsecaea and Cladophialophora appear to be heterothallic with a single copy of either MAT-1-1 or MAT-1-2 in each individual. All Capronia species are homothallic as both MAT1-1 and MAT1-2 genes were found in each single genome. The genomic synteny of the MAT-locus flanking genes (SLA2-APN2-COX13) is not conserved in black fungi as is commonly observed in Eurotiomycetes, indicating a unique genomic context for MAT in those species. The heterokaryon (het) genes expansion associated with the low selective pressure at the MAT-locus suggests that a parasexual cycle may play an important role in generating diversity among those fungi.

Comparative Genomic Analysis of Drechmeria coniospora Reveals Core and Specific Genetic Requirements for Fungal Endoparasitism of Nematodes

Drechmeria coniospora is an obligate fungal pathogen that infects nematodes via the adhesion of specialized spores to the host cuticle. D. coniospora is frequently found associated with Caenorhabditis elegans in environmental samples. It is used in the study of the nematode's response to fungal infection. Full understanding of this bi-partite interaction requires knowledge of the pathogen's genome, analysis of its gene expression program and a capacity for genetic engineering. The acquisition of all three is reported here. A phylogenetic analysis placed D. coniospora close to the truffle parasite Tolypocladium ophioglossoides, and Hirsutella minnesotensis, another nematophagous fungus. Ascomycete nematopathogenicity is polyphyletic; D. coniospora represents a branch that has not been molecularly characterized. A detailed in silico functional analysis, comparing D. coniospora to 11 fungal species, revealed genes and gene families potentially involved in virulence and showed it to be a highly specialized pathogen. A targeted comparison with nematophagous fungi highlighted D. coniospora-specific genes and a core set of genes associated with nematode parasitism. A comparative gene expression analysis of samples from fungal spores and mycelia, and infected C. elegans, gave a molecular view of the different stages of the D. coniospora lifecycle. Transformation of D. coniospora allowed targeted gene knock-out and the production of fungus that expresses fluorescent reporter genes. It also permitted the initial characterisation of a potential fungal counter-defensive strategy, involving interference with a host antimicrobial mechanism. This high-quality annotated genome for D. coniospora gives insights into the evolution and virulence of nematode-destroying fungi. Coupled with genetic transformation, it opens the way for molecular dissection of D. coniospora physiology, and will allow both sides of the interaction between D. coniospora and C. elegans, as well as the evolutionary arms race that exists between pathogen and host, to be studied.

Proteomic response of the phytopathogen Phyllosticta citricarpa to antimicrobial volatile organic compounds from Saccharomyces cerevisiae

Volatile organic compounds (VOCs) released by Saccharomyces cerevisiae inhibit plant pathogens, including the filamentous fungus Phyllosticta citricarpa, causal agent of citrus black spot. VOCs mediate relevant interactions between organisms in nature, and antimicrobial VOCs are promising, environmentally safer fumigants to control phytopathogens. As the mechanisms by which VOCs inhibit microorganisms are not well characterized, we evaluated the proteomic response in P. citricarpa after exposure for 12h to a reconstituted mixture of VOCs (alcohols and esters) originally identified in S. cerevisiae. Total protein was extracted and separated by 2D-PAGE, and differentially expressed proteins were identified by LC-MS/MS. About 600 proteins were detected, of which 29 were downregulated and 11 were upregulated. These proteins are involved in metabolism, genetic information processing, cellular processes, and transport. Enzymes related to energy-generating pathways, particularly glycolysis and the tricarboxylic acid cycle, were the most strongly affected. Thus, the data indicate that antimicrobial VOCs interfere with essential metabolic pathways in P. citricarpa to prevent fungal growth.

Comparative Phylogenomics of Pathogenic and Nonpathogenic Species

The Ascomycete Onygenales order embraces a diverse group of mammalian pathogens, including the yeast-forming dimorphic fungal pathogens Histoplasma capsulatum, Paracoccidioides spp. and Blastomyces dermatitidis, the dermatophytes Microsporum spp. and Trichopyton spp., the spherule-forming dimorphic fungal pathogens in the genus Coccidioides, and many nonpathogens. Although genomes for all of the aforementioned pathogenic species are available, only one nonpathogen had been sequenced. Here, we enhance comparative phylogenomics in Onygenales by adding genomes for Amauroascus mutatus, Amauroascus niger, Byssoonygena ceratinophila, and Chrysosporium queenslandicum—four nonpathogenic Onygenales species, all of which are more closely related to Coccidioides spp. than any other known Onygenales species. Phylogenomic detection of gene family expansion and contraction can provide clues to fungal function but is sensitive to taxon sampling. By adding additional nonpathogens, we show that LysM domain-containing proteins, previously thought to be expanding in some Onygenales, are contracting in the Coccidioides-Uncinocarpus clade, as are the self-nonself recognition Het loci. The denser genome sampling presented here highlights nearly 800 genes unique to Coccidiodes, which have significantly fewer known protein domains and show increased expression in the endosporulating spherule, the parasitic phase unique to Coccidioides spp. These genomes provide insight to gene family expansion/contraction and patterns of individual gene gain/loss in this diverse order—both major drivers of evolutionary change. Our results suggest that gene family expansion/contraction can lead to adaptive radiations that create taxonomic orders, while individual gene gain/loss likely plays a more significant role in branch-specific phenotypic changes that lead to adaptation for species or genera.

Comparative Analysis of Secretomes from Ectomycorrhizal Fungi with an Emphasis on Small-Secreted Proteins

Fungi are major players in the carbon cycle in forest ecosystems due to the wide range of interactions they have with plants either through soil degradation processes by litter decayers or biotrophic interactions with pathogenic and ectomycorrhizal symbionts. Secretion of fungal proteins mediates these interactions by allowing the fungus to interact with its environment and/or host. Ectomycorrhizal (ECM) symbiosis independently appeared several times throughout evolution and involves approximately 80% of trees. Despite extensive physiological studies on ECM symbionts, little is known about the composition and specificities of their secretomes. In this study, we used a bioinformatics pipeline to predict and analyze the secretomes of 49 fungal species, including 11 ECM fungi, wood and soil decayers and pathogenic fungi to tackle the following questions: (1) Are there differences between the secretomes of saprophytic and ECM fungi? (2) Are small-secreted proteins (SSPs) more abundant in biotrophic fungi than in saprophytic fungi? and (3) Are there SSPs shared between ECM, saprotrophic and pathogenic fungi? We showed that the number of predicted secreted proteins is similar in the surveyed species, independently of their lifestyle. The secretome from ECM fungi is characterized by a restricted number of secreted CAZymes, but their repertoires of secreted proteases and lipases are similar to those of saprotrophic fungi. Focusing on SSPs, we showed that the secretome of ECM fungi is enriched in SSPs compared with other species. Most of the SSPs are coded by orphan genes with no known PFAM domain or similarities to known sequences in databases. Finally, based on the clustering analysis, we identified shared- and lifestyle-specific SSPs between saprotrophic and ECM fungi. The presence of SSPs is not limited to fungi interacting with living plants as the genome of saprotrophic fungi also code for numerous SSPs. ECM fungi shared lifestyle-specific SSPs likely involved in symbiosis that are good candidates for further functional analyses.

Genome and secretome analyses provide insights into keratin decomposition by novel proteases from the non-pathogenic fungus Onygena corvina

Poultry processing plants and slaughterhouses produce huge quantities of feathers and hair/bristle waste annually. These keratinaceous wastes are highly resistant to degradation. Onygena corvina, a non-pathogenic fungus, grows specifically on feathers, hooves, horn, and hair in nature. Hence, the proteases secreted by O. corvina are interesting in view of their potential relevance for industrial decomposition of keratinaceous wastes. We sequenced and assembled the genome of O. corvina and used a method called peptide pattern recognition to identify 73 different proteases. Comparative genome analysis of proteases in keratin-degrading and non-keratin-degrading fungi indicated that 18 putative secreted proteases from four protease families (M36, M35, M43, and S8) may be responsible for keratin decomposition. Twelve of the 18 predicted protease genes could be amplified from O. corvina grown on keratinaceous materials and were transformed into Pichia pastoris. One of the recombinant proteases belonging to the S8 family showed high keratin-degrading activity. Furthermore, 29 different proteases were identified by mass spectrometry in the culture broth of O. corvina grown on feathers and bristle. The culture broth was fractionated by ion exchange chromatography to isolate active fractions with five novel proteases belonging to three protease families (S8, M28, and M3). Enzyme blends composed of three of these five proteases, one from each family, showed high degree of degradation of keratin in vitro. A blend of novel proteases, such as those we discovered, could possibly find a use for degrading keratinaceous wastes and provide proteins, peptides, and amino acids as valuable ingredients for animal feed.

Characterization of the Carbohydrate Binding Module 18 gene family in the amphibian pathogen Batrachochytrium dendrobatidis

Batrachochytrium dendrobatidis (Bd) is the causative agent of chytridiomycosis responsible for worldwide decline in amphibian populations. Previous analysis of the Bd genome revealed a unique expansion of the carbohydrate-binding module family 18 (CBM18) predicted to be a sub-class of chitin recognition domains. CBM expansions have been linked to the evolution of pathogenicity in a variety of fungal species by protecting the fungus from the host. Based on phylogenetic analysis and presence of additional protein domains, the gene family can be classified into 3 classes: Tyrosinase-, Deacetylase-, and Lectin-like. Examination of the mRNA expression levels from sporangia and zoospores of nine of the cbm18 genes found that the Lectin-like genes had the highest expression while the Tyrosinase-like genes showed little expression, especially in zoospores. Heterologous expression of GFP-tagged copies of four CBM18 genes in Saccharomyces cerevisiae demonstrated that two copies containing secretion signal peptides are trafficked to the cell boundary. The Lectin-like genes cbm18-ll1 and cbm18-ll2 co-localized with the chitinous cell boundaries visualized by staining with calcofluor white. In vitro assays of the full length and single domain copies from CBM18-LL1 demonstrated chitin binding and no binding to cellulose or xylan. Expressed CBM18 domain proteins were demonstrated to protect the fungus, Trichoderma reeseii, in vitro against hydrolysis from exogenously added chitinase, likely by binding and limiting exposure of fungal chitin. These results demonstrate that cbm18 genes can play a role in fungal defense and expansion of their copy number may be an important pathogenicity factor of this emerging infectious disease of amphibians.

Isolation and identification of an extracellular subtilisin-like serine protease secreted by the bat pathogen Pseudogymnoascus destructans

White nose syndrome (WNS) is a cutaneous fungal disease of bats. WNS is responsible for unprecedented mortalities in North American cave bat populations. There have been few descriptions of enzyme activities that may function in WNS host/pathogen interactions, while no study has isolated and described secreted proteases. To address the hypothesis that Pseudogymnoascus destructans secretes extracellular proteases that function in wing necrosis during WNS infection, the object of this study was to culture P. destructans on various media, then isolate and structurally identify those proteases accumulated stably in the culture medium. We found a single dominant protease activity on minimal nutrient broth enriched with protein substrates, which was strongly inhibited by phenylmethylsulfonyl fluoride. This P. destructans serine protease (PdSP1) was isolated by preparative isoelectric focusing and concanavalin A lectin affinity chromatography. PdSP1 showed a molecular weight 27,900 (estimated by SDS-PAGE), broad pH optimum 6-8, and temperature optimum 60°C. Structural characterization of PdSP1 by MALDI-TOF MS, Orbitrap MS/MS, and Edman amino-terminal peptide sequencing matched it directly to a hypothetical protein accession from the sequenced P. destructans genome that is further identified as a MEROPS family S8A subtilisin-like serine peptidase. Two additional isoforms, PdSP2 and PdSP3, were identified in the P. destructans genome with 90% and 53% homology, respectively. P. destructans S8A serine proteases showed closer sequence conservation to P. pannorum and plant pathogenic fungi than to human pathogenic dermatophytes. Peptide-specific polyclonal antibodies developed from the PdSP1 sequence detected the protein in western blots. These subtilisin-like serine proteases are candidates for further functional studies in WNS host-pathogen interaction.

Characterisation of the secretome of the clam parasite, QPX

Secreted and cell surface-associated molecules play a major role in disease development processes and host-pathogen interactions, and usually determine the virulence of invading organisms. In this study, we investigated proteins secreted by quahog parasite unknown, a thraustochytrid protist that infects the hard clam, Mercenaria mercenaria. In silico analysis of quahog parasite unknown transcripts predicted over 1200 proteins to possess an amino-terminal signal peptide which directs proteins into the classical eukaryotic secretory pathway. Proteomic analysis using LC/MS technology identified 56 proteins present in the extracellular secretion of quahog parasite unknown cells grown in vitro, including six mucin-like molecules, four glycosyl hydrolases and eight peptidases. Transcription levels of 19 quahog parasite unknown extracellular proteins were investigated in clam tissue lesions (in vivo) using quantitative PCR. The overexpression of six of these extracellular proteins in clam tissues compared with in vitro cultures suggests that they are involved in interaction with the clam host.

Genome-wide and molecular evolution analysis of the subtilase gene family in Vitis vinifera

Background

Vitis vinifera (grape) is one of the most economically significant fruit crops in the world. The availability of the recently released grape genome sequence offers an opportunity to identify and analyze some important gene families in this species. Subtilases are a group of subtilisin-like serine proteases that are involved in many biological processes in plants. However, no comprehensive study incorporating phylogeny, chromosomal location and gene duplication, gene organization, functional divergence, selective pressure and expression profiling has been reported so far for the grape.

Results

In the present study, a comprehensive analysis of the subtilase gene family in V. vinifera was performed. Eighty subtilase genes were identified. Phylogenetic analyses indicated that these subtilase genes comprised eight groups. The gene organization is considerably conserved among the groups. Distribution of the subtilase genes is non-random across the chromosomes. A high proportion of these genes are preferentially clustered, indicating that tandem duplications may have contributed significantly to the expansion of the subtilase gene family. Analyses of divergence and adaptive evolution show that while purifying selection may have been the main force driving the evolution of grape subtilases, some of the critical sites responsible for the divergence may have been under positive selection. Further analyses of real-time PCR data suggested that many subtilase genes might be important in the stress response and functional development of plants.

Conclusions

Tandem duplications as well as purifying and positive selections have contributed to the functional divergence of subtilase genes in V. vinifera. The data may contribute to a better understanding of the grape subtilase gene family.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1116) contains supplementary material, which is available to authorized users.

Evolutionary Perspectives on Human Fungal Pathogens

Owing to their small size and paucity of phenotypic characters, progress in the evolutionary biology of microbes in general, and human pathogenic fungi in particular, has been linked to a series of advances in DNA sequencing over the past quarter century. Phylogenetics was the first area to benefit, with the achievement of a basic understanding of fungal phylogeny. Population genetics was the next advance, finding cryptic species everywhere, and recombination in species previously thought to be asexual. Comparative genomics saw the next advance, in which variation in gene content and changes in gene family size were found to be important sources of variation. Fungal population genomics is showing that gene flow among closely related populations and species provides yet another source of adaptive, genetic variation. Now, two means to associate genetic variation with phenotypic variation, "reverse ecology" for adaptive phenotypes, and genome-wide association of any phenotype, are letting evolutionary biology make a profound contribution to molecular developmental biology of pathogenic fungi.

Comparative genomics and transcriptomics analyses reveal divergent lifestyle features of nematode endoparasitic fungus Hirsutella minnesotensis

Hirsutella minnesotensis [Ophiocordycipitaceae (Hypocreales, Ascomycota)] is a dominant endoparasitic fungus by using conidia that adhere to and penetrate the secondary stage juveniles of soybean cyst nematode. Its genome was de novo sequenced and compared with five entomopathogenic fungi in the Hypocreales and three nematode-trapping fungi in the Orbiliales (Ascomycota). The genome of H. minnesotensis is 51.4 Mb and encodes 12,702 genes enriched with transposable elements up to 32%. Phylogenomic analysis revealed that H. minnesotensis was diverged from entomopathogenic fungi in Hypocreales. Genome of H. minnesotensis is similar to those of entomopathogenic fungi to have fewer genes encoding lectins for adhesion and glycoside hydrolases for cellulose degradation, but is different from those of nematode-trapping fungi to possess more genes for protein degradation, signal transduction, and secondary metabolism. Those results indicate that H. minnesotensis has evolved different mechanism for nematode endoparasitism compared with nematode-trapping fungi. Transcriptomics analyses for the time-scale parasitism revealed the upregulations of lectins, secreted proteases and the genes for biosynthesis of secondary metabolites that could be putatively involved in host surface adhesion, cuticle degradation, and host manipulation. Genome and transcriptome analyses provided comprehensive understanding of the evolution and lifestyle of nematode endoparasitism.

Characterization of the transcriptome and temperature-induced differential gene expression in QPX, the thraustochytrid parasite of hard clams

BACKGROUND

The hard clam or northern quahog, Mercenaria mercenaria, is one of the most valuable seafood products in the United States representing the first marine resource in some Northeastern states. Severe episodes of hard clam mortality have been consistently associated with infections caused by a thraustochytrid parasite called Quahog Parasite Unknown (QPX). QPX is considered as a cold/temperate water organism since the disease occurs only in the coastal waters of the northwestern Atlantic Ocean from Maritime Canada to Virginia. High disease development at cold temperatures was also confirmed in laboratory studies and is thought to be caused predominantly by immunosuppression of the clam host even though the effect of temperature on QPX virulence has not been fully investigated. In this study, the QPX transcriptome was sequenced using Roche 454 technology to better characterize this microbe and initiate research on the molecular basis of QPX virulence towards hard clams.

RESULTS

Close to 18,000 transcriptomic sequences were generated and functionally annotated. Results revealed a wide array of QPX putative virulence factors including a variety of peptidases, antioxidant enzymes, and proteins involved in extracellular mucus production and other secretory proteins potentially involved in interactions with the clam host. Furthermore, a 15 K oligonucleotide array was constructed and used to investigate the effect of temperature on QPX fitness and virulence factors. Results identified a set of QPX molecular chaperones that could explain its adaptation to cold temperatures. Finally, several virulence-related factors were up-regulated at low temperature providing molecular targets for further investigations of increased QPX pathogenicity in cold water conditions.

CONCLUSIONS

This is one of the first studies to characterize the transcriptome of a parasitic labyrinthulid, offering new insights into the molecular bases of the pathogenicity of members of this group. Results from the oligoarray study demonstrated the ability of QPX to cope with a wide range of environmental temperatures, including those considered to be suboptimal for clam immunity (low temperature) providing a mechanistic scenario for disease distribution in the field and for high disease prevalence and intensity at low temperature. These results will serve as basis for studies aimed at a better characterization of specific putative virulence factors.

Genomics in Coccidioides: insights into evolution, ecology, and pathogenesis

Coccidioides immitis and C. posadasii, the causative agents of the mammalian disease coccidioidomycosis, are dimorphic fungal pathogens distributed throughout desert-like environments in North and South America. Coccidioides spp. are members of the Onygenales, a diverse group of pathogenic and nonpathogenic fungi. Recently, full genomes have been published for Coccidioides and a number of other Onygenales species. Phylogenomic comparisons and additional studies in Coccidioides population genomics and gene expression have shed light on the ecology and pathogenesis of Coccidioides and the other medically important species in this clade. Observed patterns of gene family expansion/contraction and evidence of gene flow have provided insight to the evolution of Coccidioides and greatly broadened our understanding of the diversity and sources of genetic variation found in fungi. In the future, expansion of the number of sequenced isolates from all populations will allow deeper insight into the evolutionary processes that have shaped this unique human pathogen. In addition, deep sequencing of isolates from a single Coccidioides population and pairing of those data with phenotype information on growth and pathogenicity for genome-wide association analysis will allow researchers to find genes responsible for any phenotype, virulence included, that shows variation in the population.

Horizontal transfer of a subtilisin gene from plants into an ancestor of the plant pathogenic fungal genus Colletotrichum

The genus Colletotrichum contains a large number of phytopathogenic fungi that produce enormous economic losses around the world. The effect of horizontal gene transfer (HGT) has not been studied yet in these organisms. Inter-Kingdom HGT into fungal genomes has been reported in the past but knowledge about the HGT between plants and fungi is particularly limited. We describe a gene in the genome of several species of the genus Colletotrichum with a strong resemblance to subtilisins typically found in plant genomes. Subtilisins are an important group of serine proteases, widely distributed in all of the kingdoms of life. Our hypothesis is that the gene was acquired by Colletotrichum spp. through (HGT) from plants to a Colletotrichum ancestor. We provide evidence to support this hypothesis in the form of phylogenetic analyses as well as a characterization of the similarity of the subtilisin at the primary, secondary and tertiary structural levels. The remarkable level of structural conservation of Colletotrichum plant-like subtilisin (CPLS) with plant subtilisins and the differences with the rest of Colletotrichum subtilisins suggests the possibility of molecular mimicry. Our phylogenetic analysis indicates that the HGT event would have occurred approximately 150–155 million years ago, after the divergence of the Colletotrichum lineage from other fungi. Gene expression analysis shows that the gene is modulated during the infection of maize by C. graminicola suggesting that it has a role in plant disease. Furthermore, the upregulation of the CPLS coincides with the downregulation of several plant genes encoding subtilisins. Based on the known roles of subtilisins in plant pathogenic fungi and the gene expression pattern that we observed, we postulate that the CPLSs have an important role in plant infection.

Comparative genome analysis of Trichophyton rubrum and related dermatophytes reveals candidate genes involved in infection

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.

Genome sequence of an unclassified pleosporales species isolated from human nasopharyngeal aspirate

Pleosporales is the largest order in the fungal class Dothideomycetes. We report the 36,814,818-bp draft genome sequence and gene annotation of UM1110, a Pleosporales isolate associated with unclassified genera that is potentially a new fungal species. Analysis of the genome sequence led to the finding of genes associated with fungal adhesive proteins, secreted proteases, allergens, and pseudohyphal development.