A phylogenomic approach to reconstructing the diversification of serine proteases in fungi

In silico structural exploration of serine protease from a CTG-clade yeast Meyerozyma guilliermondii strain SO

In eukaryotes, serine proteases are cellular localized hydrolases reported to regulate essential biological reactions. Improved industrial applications of proteins are aided by prediction and analysis of their 3-dimensional structures (3D). A serine protease was identified from CTG-clade yeast Meyerozyma guilliermondii strain SO and its 3D structure as well as its catalytic attributes have not been fully understood yet, thus we seek to report on the catalytic mechanism of M. guilliermondii strain SO MgPRB1 using substrate PMSF via in silico docking as well as its stability by way of disulfide bonds formation. Herein, bioinformatics tools and techniques were used to predict, validate and analyze the possible changes of CUG ambiguity (if any) in strain SO using template PDB ID: 3F7O. Structural assessments confirmed the classic catalytic triad Asp305, His337, and Ser499. Superimposition of MgPRB1 and template 3F7O structures revealed the unlinked cysteine residues between Cys341, Cys440, Cys471 and Cys506 of MgPRB1 compared to template 3F7O with two disulfide bonds formation, which confers structural stability. In conclusion, serine protease structure from strain SO was successfully predicted and studies towards understanding at the molecular level may be undertaken for its potential applications in the degradation of peptide bonds.

Cloning and sequence analysis of a serine protease gene from Rhizoctonia solani Kühn AG5

The present study aimed to identify the subtilisin-like proteases (SLPs) of Rhizoctonia solani Kühn potentially involved in the virulence of this phytopathogenic fungus, which has 14 anastomosis groups (AGs) responsible for many crop diseases. Through mycelial microscope observation and strain identification of pathogenic fungus MS-3, it was determined to be R. solani AG-5. Both 5' and 3' rapid amplification of cDNA ends were used to clone the serine protease gene RsSLP from R. solani AG-5. The full-length obtained for RsSLP was 1714 bp with an open reading frame of 1587 bp, encoding a protein of 528 amino acids with a molecular mass of 55.8 kDa. This protein contained a predicted signal peptide for secretion but lacked a transmembrane domain or membrane anchor site. Bioinformatics analysis identified this protein as a serine protease with the Peptidase_S8 and Inhibitor_I9 characteristic domains of SLPs. Phylogenetic analysis suggested that frequent gene duplications of the SLPs occurred in R. solani (RsSLP), and RsSLP shares characteristic sequence features with virulence factors of other phytopathogenic fungi. Because the secretory serine protease RsSLP from R. solani AG5 is similar to the virulence factors of other phytopathogenic fungi, its identification will be helpful in studies considering the roles of these proteases in pathogen virulence.

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.

Structured Framework and Genome Analysis of Magnaporthe grisea Inciting Pearl Millet Blast Disease Reveals Versatile Metabolic Pathways, Protein Families, and Virulence Factors

Magnaporthe grisea (T.T. Herbert) M.E. Barr is a major fungal phytopathogen that causes blast disease in cereals, resulting in economic losses worldwide. An in-depth understanding of the basis of virulence and ecological adaptation of M. grisea is vital for devising effective disease management strategies. Here, we aimed to determine the genomic basis of the pathogenicity and underlying biochemical pathways in Magnaporthe using the genome sequence of a pearl millet-infecting M. grisea PMg_Dl generated by dual NGS techniques, Illumina NextSeq 500 and PacBio RS II. The short and long nucleotide reads could be draft assembled in 341 contigs and showed a genome size of 47.89 Mb with the N50 value of 765.4 Kb. Magnaporthe grisea PMg_Dl showed an average nucleotide identity (ANI) of 86% and 98% with M. oryzae and Pyricularia pennisetigena, respectively. The gene-calling method revealed a total of 10,218 genes and 10,184 protein-coding sequences in the genome of PMg_Dl. InterProScan of predicted protein showed a distinct 3637 protein families and 695 superfamilies in the PMg_Dl genome. In silico virulence analysis revealed the presence of 51VFs and 539 CAZymes in the genome. The genomic regions for the biosynthesis of cellulolytic endo-glucanase and beta-glucosidase, as well as pectinolytic endo-polygalacturonase, pectin-esterase, and pectate-lyases (pectinolytic) were detected. Signaling pathways modulated by MAPK, PI3K-Akt, AMPK, and mTOR were also deciphered. Multicopy sequences suggestive of transposable elements such as Type LTR, LTR/Copia, LTR/Gypsy, DNA/TcMar-Fot1, and Type LINE were recorded. The genomic resource presented here will be of use in the development of molecular marker and diagnosis, population genetics, disease management, and molecular taxonomy, and also provide a genomic reference for ascomycetous genome investigations in the future.

Genome-Wide Analysis and Characterization of the Riemerella anatipestifer Putative T9SS Secretory Proteins with a Conserved C-Terminal Domain

Riemerella anatipestifer is a major pathogenic agent of duck septicemic and exudative diseases. Recent studies have shown that the R. anatipestifer type IX secretion system (T9SS) acts as a crucial virulence factor. We previously identified two T9SS component proteins, GldK and GldM, and one T9SS effector metallophosphoesterase, which play important roles in bacterial virulence. In this study, 19 T9SS-secreted proteins that contained a conserved T9SS C-terminal domain (CTD) were predicted in R. anatipestifer strain Yb2 by searching for CTD-encoding sequences in the whole genome. The proteins were confirmed with a liquid chromatography-tandem mass spectrometry analysis of the bacterial culture supernatant. Nine of them were reported in our previous study. We generated recombinant proteins and mouse antisera for the 19 predicted proteins to confirm their expression in the bacterial culture supernatant and in bacterial cells. Western blotting indicated that the levels of 14 proteins were significantly reduced in the T9SS mutant Yb2ΔgldM culture medium but were increased in the bacterial cells. RT-qPCR indicated that the expression of these genes did not differ between the wild-type strain Yb2 and the T9SS mutant Yb2ΔgldM. Nineteen mutant strains were successfully constructed to determine their virulence and proteolytic activity, which indicated that seven proteins are associated with bacterial virulence, and two proteins, AS87_RS04190 and AS87_RS07295, are protease-activity-associated virulence factors. In summary, we have identified at least 19 genes encoding T9SS-secreted proteins in the R. anatipestifer strain Yb2 genome, which encode multiple functions associated with the bacterium's virulence and proteolytic activity. IMPORTANCE Riemerella anatipestifer T9SS plays an important role in bacterial virulence. We have previously reported nine R. anatipestifer T9SS-secreted proteins and clarified the function of the metallophosphoesterase. In this study, we identified 10 more secreted proteins associated with the R. anatipestifer T9SS, in addition to the nine previously reported. Of these, 14 proteins showed significantly reduced secretion into the bacterial culture medium but increased expression in the bacterial cells of the T9SS mutant Yb2ΔgldM; seven proteins were shown to be associated with bacterial virulence; and two proteins, AS87_RS04190 and AS87_RS07295, were shown to be protease-activity-associated virulence factors. Thus, we have demonstrated that multiple R. anatipestifer T9SS-secreted proteins function in virulence and proteolytic activity.

Protease production and molecular characterization of a protease dipeptidyl-aminopeptidase gene from different strains of Sordaria fimicola

The current research was designed to reach extracellular protease production potential in different strains of Sordaria fimicola which were previously obtained from Dr. Lamb (Imperial College, London) from North Facing Slope and South Facing Slope of Evolution Canyon. After initial and secondary screening, two hyper-producers strains S2 and N6 were selected for submerged fermentation and cultural conditions including temperature, pH, incubation period, inoculum size, substrate concentration, and different carbon and nitrogen sources were optimized for enzyme production. S2 strain showed maximum protease production of 3.291 U/mL after 14 days of incubation at 30 °C with 7 pH, 1% substrate concentration and 1 mL inoculum, While N6 strain showed maximum protease production of 1.929 U/mL under fermentation optimized conditions. Another aim of the present research was to underpin the biodiversity of genetics and post-translational modifications (PTMs) of protease DPAP (peptidyl-aminopeptidase) in Sordaria fimicola. Five polymorphic sites were observed in amino acid sequence of S. fimicola strains with reference to Neurospora crassa. PTMs prediction from bioinformatics tools predicted 38 phosphorylation sites on serine residues for protease peptidyl-aminopeptidase in S1 strain of S. fimicola while 45 phosphorylation sites on serine in N7 strain and 47 serine phosphorylation modifications were predicted in N. crassa. Current research gave an insight that change in genetic makeup effected PTMs which ultimately affected the production of protease enzyme in different strains of same organism (S. fimicola). The production and molecular data of the research revealed that environmental stress has strong effects on the specific genes through mutations which may cause genetic diversity. S. fimicola is non- pathogenic fungus and has a short life cycle. This fungus can be chosen to produce protease enzyme on a commercial scale.

Riemerella anatipestifer T9SS Effector SspA Functions in Bacterial Virulence and Defending Natural Host Immunity

Riemerella anatipestifer is a major pathogenic agent of duck septicemic and exudative diseases. Recent studies have shown that the R. anatipestifer type IX secretion system (T9SS) is a crucial factor in bacterial virulence. The AS87_RS04190 protein was obviously missing from the secreted proteins of the T9SS mutant strain Yb2ΔgldM. A bioinformatic analysis indicated that the AS87_RS04190 protein contains a T9SS C-terminal domain sequence and encodes a putative subtilisin-like serine protease (SspA). To determine the role of the putative SspA protein in R. anatipestifer pathogenesis and proteolysis, we constructed two strains with an sspA mutation and complementation, respectively, and determined their median lethal doses, their bacterial loads in infected duck blood, and their adherence to and invasion of cells. Our results demonstrate that the SspA protein functions in bacterial virulence. It is also associated with the bacterial protease activity and has a conserved catalytic triad structure (Asp126, His158, and Ser410), which is necessary for protein function. The optimal reactive pH and temperature were determined to be 7.0 and 50°C, respectively, and K_m and V_max were determined to be 10.15 mM and 246.96 U/mg, respectively. The enzymatic activity of SspA is activated by Ca²⁺, Mg²⁺, and Mn²⁺ and inhibited by Cu²⁺ and EDTA. SspA degrades gelatin, fibrinogen, and bacitracin LL-37. These results demonstrate that SspA is an effector protein of T9SS and functions in R. anatipestifer virulence and its proteolysis of gelatin, fibrinogen, and bacitracin LL-37. IMPORTANCE In recent years, Riemerella anatipestifer T9SS has been reported to act as a virulence factor. However, the functions of the proteins secreted by R. anatipestifer T9SS are not entirely clear. In this study, a secreted subtilisin-like serine protease SspA was shown to be associated with R. anatipestifer virulence, host complement evasion, and degradation of gelatin, fibrinogen, and LL-37. The enzymatic activity of recombinant SspA was determined, and its K_m and V_max were 10.15 mM and 246.96 U/mg, respectively. Three conserved sites (Asp126, His158, and Ser410) are necessary for the protein's function. The median lethal dose of the sspA-deleted mutant strain was reduced >10,000-fold, indicating that SspA is an important virulence factor. In summary, we demonstrate that the R. anatipestifer AS87_RS04190 gene encodes an important T9SS effector, SspA, which plays an important role in bacterial virulence.

Proteomic Analysis of Trichomonas vaginalis Phagolysosome, Lysosomal Targeting, and Unconventional Secretion of Cysteine Peptidases

The lysosome represents a central degradative compartment of eukaryote cells, yet little is known about the biogenesis and function of this organelle in parasitic protists. Whereas the mannose 6-phosphate (M6P)-dependent system is dominant for lysosomal targeting in metazoans, oligosaccharide-independent sorting has been reported in other eukaryotes. In this study, we investigated the phagolysosomal proteome of the human parasite Trichomonas vaginalis, its protein targeting and the involvement of lysosomes in hydrolase secretion. The organelles were purified using Percoll and OptiPrep gradient centrifugation and a novel purification protocol based on the phagocytosis of lactoferrin-covered magnetic nanoparticles. The analysis resulted in a lysosomal proteome of 462 proteins, which were sorted into 21 classes. Hydrolases represented the largest functional class and included proteases, lipases, phosphatases, and glycosidases. Identification of a large set of proteins involved in vesicular trafficking (80) and turnover of actin cytoskeleton rearrangement (29) indicate a dynamic phagolysosomal compartment. Several cysteine proteases such as TvCP2 were previously shown to be secreted. Our experiments showed that secretion of TvCP2 was strongly inhibited by chloroquine, which increases intralysosomal pH, thus indicating that TvCP2 secretion occurs through lysosomes rather than the classical secretory pathway. Unexpectedly, we identified divergent homologues of the M6P receptor TvMPR in the phagolysosomal proteome, although T. vaginalis lacks enzymes for M6P formation. To test whether oligosaccharides are involved in lysosomal targeting, we selected the lysosome-resident cysteine protease CLCP, which possesses two glycosylation sites. Mutation of any of the sites redirected CLCP to the secretory pathway. Similarly, the introduction of glycosylation sites to secreted β-amylase redirected this protein to lysosomes. Thus, unlike other parasitic protists, T. vaginalis seems to utilize glycosylation as a recognition marker for lysosomal hydrolases. Our findings provide the first insight into the complexity of T. vaginalis phagolysosomes, their biogenesis, and role in the unconventional secretion of cysteine peptidases.

Co-option of an extracellular protease for transcriptional control of nutrient degradation in the fungus Aspergillus nidulans

Nutrient acquisition is essential for all organisms. Fungi regulate their metabolism according to environmental nutrient availability through elaborate transcription regulatory programs. In filamentous fungi, a highly conserved GATA transcription factor AreA and its co-repressor NmrA govern expression of genes involved in extracellular breakdown, uptake, and metabolism of nitrogen nutrients. Here, we show that the Aspergillus nidulans PnmB protease is a moonlighting protein with extracellular and intracellular functions for nitrogen acquisition and metabolism. PnmB serves not only as a secreted protease to degrade extracellular nutrients, but also as an intracellular protease to control the turnover of the co-repressor NmrA, accelerating AreA transcriptional activation upon nitrogen starvation. PnmB expression is controlled by AreA, which activates a positive feedback regulatory loop. Hence, we uncover a regulatory mechanism in the well-established controls determining the response to nitrogen starvation, revealing functional evolution of a protease gene for transcriptional regulation and extracellular nutrient breakdown.

Genomic Analysis of the Mycoparasite Pestalotiopsis sp. PG52

Pestalotiopsis sp. is a mycoparasite of the plant pathogen Aecidium wenshanense. To further understand the mycoparasitism mechanism of Pestalotiopsis sp., we assembled and analyzed its genome. The genome of Pestalotiopsis sp. strain PG52 was assembled into 335 scaffolds and had a size of 58.01 Mb. A total of 20,023 predicted genes and proteins were annotated. This study compared PG52 with the mycoparasites Trichoderma harzianum, Trichoderma atroviride, and Trichoderma virens. This study reveals the entirely different mycoparasitism mechanism of Pestalotiopsis compared to Trichoderma and reveals this mycoparasite’s strong ability to produce secondary metabolites.

Metarhizium: jack of all trades, master of many

The genus Metarhizium and Pochonia chlamydosporia comprise a monophyletic clade of highly abundant globally distributed fungi that can transition between long-term beneficial associations with plants to transitory pathogenic associations with frequently encountered protozoans, nematodes or insects. Some very common 'specialist generalist' species are adapted to particular soil and plant ecologies, but can overpower a wide spectrum of insects with numerous enzymes and toxins that result from extensive gene duplications made possible by loss of meiosis and associated genome defence mechanisms. These species use parasexuality instead of sex to combine beneficial mutations from separate clonal individuals into one genome (Vicar of Bray dynamics). More weakly endophytic species which kill a narrow range of insects retain sexuality to facilitate host-pathogen coevolution (Red Queen dynamics). Metarhizium species can fit into numerous environments because they are very flexible at the genetic, physiological and ecological levels, providing tractable models to address how new mechanisms for econutritional heterogeneity, host switching and virulence are acquired and relate to diverse sexual life histories and speciation. Many new molecules and functions have been discovered that underpin Metarhizium associations, and have furthered our understanding of the crucial ecology of these fungi in multiple habitats.

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 reconstruction of the non-culturable spinach downy mildew Peronospora effusa by metagenome filtering

Peronospora effusa (previously known as P. farinosa f. sp. spinaciae, and here referred to as Pfs) is an obligate biotrophic oomycete that causes downy mildew on spinach (Spinacia oleracea). To combat this destructive many disease resistant cultivars have been bred and used. However, new Pfs races rapidly break the employed resistance genes. To get insight into the gene repertoire of Pfs and identify infection-related genes, the genome of the first reference race, Pfs1, was sequenced, assembled, and annotated. Due to the obligate biotrophic nature of this pathogen, material for DNA isolation can only be collected from infected spinach leaves that, however, also contain many other microorganisms. The obtained sequences can, therefore, be considered a metagenome. To filter and obtain Pfs sequences we utilized the CAT tool to taxonomically annotate ORFs residing on long sequences of a genome pre-assembly. This study is the first to show that CAT filtering performs well on eukaryotic contigs. Based on the taxonomy, determined on multiple ORFs, contaminating long sequences and corresponding reads were removed from the metagenome. Filtered reads were re-assembled to provide a clean and improved Pfs genome sequence of 32.4 Mbp consisting of 8,635 scaffolds. Transcript sequencing of a range of infection time points aided the prediction of a total of 13,277 gene models, including 99 RxLR(-like) effector, and 14 putative Crinkler genes. Comparative analysis identified common features in the predicted secretomes of different obligate biotrophic oomycetes, regardless of their phylogenetic distance. Their secretomes are generally smaller, compared to hemi-biotrophic and necrotrophic oomycete species. We observe a reduction in proteins involved in cell wall degradation, in Nep1-like proteins (NLPs), proteins with PAN/apple domains, and host translocated effectors. The genome of Pfs1 will be instrumental in studying downy mildew virulence and for understanding the molecular adaptations by which new isolates break spinach resistance.

Trichoderma Species: Versatile Plant Symbionts

Because of the need to provide food for the growing population, agricultural activity is faced with the huge challenge of counteracting the negative effects generated by adverse environmental factors and diseases caused by pathogens on crops, while avoiding environmental pollution due to the excessive use of agrochemicals. The exploitation of biological systems that naturally increase plant vigor, preparing them against biotic and abiotic stressors that also promote their growth and productivity represents a useful and viable strategy to help face these challenges. Fungi from the genus Trichoderma have been widely used in agriculture as biocontrol agents because of their mycoparasitic capacity and ability to improve plant health and protection against phytopathogens, which makes it an excellent plant symbiont. The mechanisms employed by Trichoderma include secretion of effector molecules and secondary metabolites that mediate the beneficial interaction of Trichoderma with plants, providing tolerance to biotic and abiotic stresses. Here we discuss the most recent advances in understanding the mechanisms employed by this opportunistic plant symbiont as biocontrol agent and plant growth promoter. In addition, through genome mining we approached a less explored factor that Trichoderma could be using to become successful plant symbionts, the production of phytohormones-auxins, cytokinins, abscisic acid, gibberellins, among others. This approach allowed us to detect sets of genes encoding proteins potentially involved in phytohormone biosynthesis and signaling. We discuss the implications of these findings in the physiology of the fungus and in the establishment of its interaction with plants.

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.

Genomic characterization of Trichoderma atrobrunneum (T. harzianum species complex) ITEM 908: insight into the genetic endowment of a multi-target biocontrol strain

BACKGROUND

So far, biocontrol agent selection has been performed mainly by time consuming in vitro confrontation tests followed by extensive trials in greenhouse and field. An alternative approach is offered by application of high-throughput techniques, which allow extensive screening and comparison among strains for desired genetic traits. In the genus Trichoderma, the past assignments of particular features or strains to one species need to be reconsidered according to the recent taxonomic revisions. Here we present the genome of a biocontrol strain formerly known as Trichoderma harzianum ITEM 908, which exhibits both growth promoting capabilities and antagonism against different fungal pathogens, including Fusarium graminearum, Rhizoctonia solani, and the root-knot nematode Meloidogyne incognita. By genomic analysis of ITEM 908 we investigated the occurrence and the relevance of genes associated to biocontrol and stress tolerance, providing a basis for future investigation aiming to unravel the complex relationships between genomic endowment and exhibited activities of this strain.

RESULTS

The MLST analysis of ITS-TEF1 concatenated datasets reclassified ITEM 908 as T. atrobrunneum, a species recently described within the T. harzianum species complex and phylogenetically close to T. afroharzianum and T. guizhouense. Genomic analysis revealed the presence of a broad range of genes encoding for carbohydrate active enzymes (CAZYmes), proteins involved in secondary metabolites production, peptaboils, epidithiodioxopiperazines and siderophores potentially involved in parasitism, saprophytic degradation as well as in biocontrol and antagonistic activities. This abundance is comparable to other Trichoderma spp. in the T. harzianum species complex, but broader than in other biocontrol species and in the species T. reesei, known for its industrial application in cellulase production. Comparative analysis also demonstrated similar genomic organization of major secondary metabolites clusters, as in other Trichoderma species.

CONCLUSIONS

Reported data provide a contribution to a deeper understanding of the mode of action and identification of activity-specific genetic markers useful for selection and improvement of biocontrol strains. This work will also enlarge the availability of genomic data to perform comparative studies with the aim to correlate phenotypic differences with genetic diversity of Trichoderma species.

Secreted protein extract analyses present the plant pathogen Alternaria alternata as a suitable industrial enzyme toolbox

Lignocellulosic plant biomass is the most abundant carbon source in the planet, which makes it a potential substrate for biorefinery. It consists of polysaccharides and other molecules with applications in pharmaceutical, food and feed, cosmetics, paper and textile industries. The exploitation of these resources requires the hydrolysis of the plant cell wall, which is a complex process. Aiming to discover novel fungal natural isolates with lignocellulolytic capacities, a screening for feruloyl esterase activity was performed in samples taken from different metal surfaces. An extracellular enzyme extract from the most promising candidate, the natural isolate Alternaria alternata PDA1, was analyzed. The feruloyl esterase activity of the enzyme extract was characterized, determining the pH and temperature optima (pH 5.0 and 55-60 °C, respectively), thermal stability and kinetic parameters, among others. Proteomic analyses derived from two-dimensional gels allowed the identification and classification of 97 protein spots from the extracellular proteome. Most of the identified proteins belonged to the carbohydrates metabolism group, particularly plant cell wall degradation. Enzymatic activities of the identified proteins (β-glucosidase, cellobiohydrolase, endoglucanase, β-xylosidase and xylanase) of the extract were also measured. These findings confirm A. alternata PDA1 as a promising lignocellulolytic enzyme producer.

SIGNIFICANCE

Although plant biomass is an abundant material that can be potentially utilized by several industries, the effective hydrolysis of the recalcitrant plant cell wall is not a straightforward process. As this hydrolysis occurs in nature relying almost solely on microbial enzymatic systems, it is reasonable to infer that further studies on lignocellulolytic enzymes will discover new sustainable industrial solutions. The results included in this paper provide a promising fungal candidate for biotechnological processes to obtain added value from plant byproducts and analogous substrates. Moreover, the proteomic analysis of the secretome of a natural isolate of Alternaria sp. grown in the presence of one of the most used vegetal substrates on the biofuels industry (sugar beet pulp) sheds light on the extracellular enzymatic machinery of this fungal plant pathogen, and can be potentially applied to developing new industrial enzymatic tools. This work is, to our knowledge, the first to analyze in depth the secreted enzyme extract of the plant pathogen Alternaria when grown on a lignocellulosic substrate, identifying its proteins by means of MALDI-TOF/TOF mass spectrometry and characterizing its feruloyl esterase, cellulase and xylanolytic activities.

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.

Establishment of Neurospora crassa as a host for heterologous protein production using a human antibody fragment as a model product

Background

Filamentous fungi are commonly used as production hosts for bulk enzymes in biotechnological applications. Their robust and quick growth combined with their ability to secrete large amounts of protein directly into the culture medium makes fungi appealing organisms for the generation of novel production systems. The red bread mold Neurospora crassa has long been established as a model system in basic research. It can be very easily genetically manipulated and a wealth of molecular tools and mutants are available. In addition, N. crassa is very fast growing and non-toxic. All of these features point to a high but so far untapped potential of this fungus for biotechnological applications. In this study, we used genetic engineering and bioprocess development in a design-build-test-cycle process to establish N. crassa as a production host for heterologous proteins.

Results

The human antibody fragment HT186-D11 was fused to a truncated version of the endogenous enzyme glucoamylase (GLA-1), which served as a carrier protein to achieve secretion into the culture medium. A modular expression cassette was constructed and tested under the control of different promoters. Protease activity was identified as a major limitation of the production strain, and the effects of different mutations causing protease deficiencies were compared. Furthermore, a parallel bioreactor system (1 L) was employed to develop and optimize a production process, including the comparison of different culture media and cultivation parameters. After successful optimization of the production strain and the cultivation conditions an exemplary scale up to a 10 L stirred tank reactor was performed.

Conclusions

The data of this study indicate that N. crassa is suited for the production and secretion of heterologous proteins. Controlling expression by the optimized promoter Pccg1nr in a fourfold protease deletion strain resulted in the successful secretion of the heterologous product with estimated yields of 3 mg/L of the fusion protein. The fungus could easily be cultivated in bioreactors and a first scale-up was successful. The system holds therefore much potential, warranting further efforts in optimization.

Electronic supplementary material

The online version of this article (doi:10.1186/s12934-017-0734-5) contains supplementary material, which is available to authorized users.

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.

The Insect Pathogens

Fungi are the most common disease-causing agents of insects; aside from playing a crucial role in natural ecosystems, insect-killing fungi are being used as alternatives to chemical insecticides and as resources for biotechnology and pharmaceuticals. Some common experimentally tractable genera, such as Metarhizium spp., exemplify genetic diversity and dispersal because they contain numerous intraspecific variants with distinct environmental and insect host ranges. The availability of tools for molecular genetics and multiple sequenced genomes has made these fungi ideal experimental models for answering basic questions on the genetic and genomic processes behind adaptive phenotypes. For example, comparative genomics of entomopathogenic fungi has shown they exhibit diverse reproductive modes that often determine rates and patterns of genome evolution and are linked as cause or effect with pathogenic strategies. Fungal-insect pathogens represent lifestyle adaptations that evolved numerous times, and there are significant differences in host range and pathogenic strategies between the major groups. However, typically, spores landing on the cuticle produce appressoria and infection pegs that breach the cuticle using mechanical pressure and cuticle-degrading enzymes. Once inside the insect body cavity, fungal pathogens face a potent and comprehensively studied immune defense by which the host attempts to eliminate or reduce an infection. The Fungal Kingdom stands alone in the range, extent, and complexity of their manipulation of arthropod behavior. In part, this is because most only sporulate on cadavers, so they must ensure the dying host positions itself to allow efficient transmission.

Comparative transcriptomics reveal host-specific nucleotide variation in entomophthoralean fungi

Obligate parasites are under strong selection to increase exploitation of their host to survive while evading detection by host immune defences. This has often led to elaborate pathogen adaptations and extreme host specificity. Specialization on one host, however, often incurs a trade-off influencing the capacity to infect alternate hosts. Here, we investigate host adaptation in two morphologically indistinguishable and closely related obligate specialist insect-pathogenic fungi from the phylum Entomophthoromycota, Entomophthora muscae sensu stricto and E. muscae sensu lato, pathogens of houseflies (Musca domestica) and cabbage flies (Delia radicum), respectively. We compared single nucleotide polymorphisms within and between these two E. muscae species using 12 RNA-seq transcriptomes from five biological samples. All five isolates contained intra-isolate polymorphisms that segregate in 50:50 ratios, indicative of genetic duplication events or functional diploidy. Comparative analysis of dN/dS ratios between the multinucleate E. muscae s.str. and E. muscae s.l. revealed molecular signatures of positive selection in transcripts related to utilization of host lipids and the potential secretion of toxins that interfere with the host immune response. Phylogenetic comparison with the nonobligate generalist insect-pathogenic fungus Conidiobolus coronatus revealed a gene-family expansion of trehalase enzymes in E. muscae. The main sugar in insect haemolymph is trehalose, and efficient sugar utilization was probably important for the evolutionary transition to obligate insect pathogenicity in E. muscae. These results support the hypothesis that genetically based host specialization in specialist pathogens evolves in response to the challenge of using resources and dealing with the immune system of different hosts.

Advances in Genomics of Entomopathogenic Fungi

Fungi are the commonest pathogens of insects and crucial regulators of insect populations. The rapid advance of genome technologies has revolutionized our understanding of entomopathogenic fungi with multiple Metarhizium spp. sequenced, as well as Beauveria bassiana, Cordyceps militaris, and Ophiocordyceps sinensis among others. Phylogenomic analysis suggests that the ancestors of many of these fungi were plant endophytes or pathogens, with entomopathogenicity being an acquired characteristic. These fungi now occupy a wide range of habitats and hosts, and their genomes have provided a wealth of information on the evolution of virulence-related characteristics, as well as the protein families and genomic structure associated with ecological and econutritional heterogeneity, genome evolution, and host range diversification. In particular, their evolutionary transition from plant pathogens or endophytes to insect pathogens provides a novel perspective on how new functional mechanisms important for host switching and virulence are acquired. Importantly, genomic resources have helped make entomopathogenic fungi ideal model systems for answering basic questions in parasitology, entomology, and speciation. At the same time, identifying the selective forces that act upon entomopathogen fitness traits could underpin both the development of new mycoinsecticides and further our understanding of the natural roles of these fungi in nature. These roles frequently include mutualistic relationships with plants. Genomics has also facilitated the rapid identification of genes encoding biologically useful molecules, with implications for the development of pharmaceuticals and the use of these fungi as bioreactors.

Insect Pathogenic Fungi as Endophytes

In this chapter, we explore some of the evolutionary, ecological, molecular genetics, and applied aspects of a subset of insect pathogenic fungi that also have a lifestyle as endophytes and we term endophytic insect pathogenic fungi (EIPF). We focus particularly on Metarhizium spp. and Beauveria bassiana as EIPF. The discussion of the evolution of EIPF challenges a view that these fungi were first and foremost insect pathogens that eventually evolved to colonize plants. Phylogenetic evidence shows that the lineages of EIPF are most closely related to grass endophytes that diverged c. 100MYA. We discuss the relationship between genes involved in "insect pathogenesis" and those involved in "endophytism" and provide examples of genes with potential importance in lifestyle transitions toward insect pathogenicity. That is, some genes for insect pathogenesis may have been coopted from genes involved in endophytic colonization. Other genes may be multifunctional and serve in both lifestyle capacities. The interactions of EIPF with their host plants are discussed in some detail. The genetic basis for rhizospheric competence, plant communication, and nutrient exchange is examined and we highlight, with examples, the benefits of EIPF to plants, and the potential reservoir of secondary metabolites hidden within these beneficial symbioses.

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.

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.

A novel subtilisin-like serine protease of Batrachochytrium dendrobatidis is induced by thyroid hormone and degrades antimicrobial peptides

Batrachochytrium dendrobatidis (B. dendrobatidis), a chytrid fungus, is one of the major contributors to the global amphibian decline. The fungus infects both tadpoles and adult amphibians. Tadpoles are infected in their keratinized mouthparts, and infected adults exhibit hyperkeratosis and loss of righting reflex. Infections of adults may result in death from cardiac arrest in susceptible species. Thyroid hormone plays a key role in amphibian metamorphosis. The occurrence of B. dendrobatidis in tadpoles during metamorphosis may result in exposure of the fungus to host morphogens including TH. This exposure may induce gene expression in the fungus contributing to invasion and colonization of the host. Here, we demonstrate movement of fungal zoospores toward TH. Additionally, expression of a subtilisin-like serine protease is up-regulated in B. dendrobatidis cells exposed to TH. A gene encoding this protease was cloned from B. dendrobatidis and expressed in Escherichia coli. The protein was partially purified and characterized. The similarity between subtilases of human dermatophytes and the B. dendrobatidis subtilisin-like serine protease suggests the importance of this enzyme in B. dendrobatidis pathogenicity. Cleavage of frog skin antimicrobial peptides (AMPs) by this B. dendrobatidis subtilisin-like serine protease suggests a role for this enzyme in fungal survival and colonization.

Proteases in malaria parasites - a phylogenomic perspective

Malaria continues to be one of the most devastating global health problems due to the high morbidity and mortality it causes in endemic regions. The search for new antimalarial targets is of high priority because of the increasing prevalence of drug resistance in malaria parasites. Malarial proteases constitute a class of promising therapeutic targets as they play important roles in the parasite life cycle and it is possible to design and screen for specific protease inhibitors. In this mini-review, we provide a phylogenomic overview of malarial proteases. An evolutionary perspective on the origin and divergence of these proteases will provide insights into the adaptive mechanisms of parasite growth, development, infection, and pathogenesis.B

Gene cloning, molecular modeling, and phylogenetics of serine protease P32 and serine carboxypeptidase SCP1 from nematophagous fungi Pochonia rubescens and Pochonia chlamydosporia

The fungi Pochonia chlamydosporia and Pochonia rubescens are parasites of nematode eggs and thus are biocontrol agents of nematodes. Proteolytic enzymes such as the S8 proteases VCP1 and P32, secreted during the pathogenesis of nematode eggs, are major virulence factors in these fungi. Recently, expression of these enzymes and of SCP1, a new putative S10 carboxypeptidase, was detected during endophytic colonization of barley roots by these fungi. In our study, we cloned the genomic and mRNA sequences encoding P32 from P. rubescens and SCP1 from P. chlamydosporia. P32 showed a high homology with the serine proteases Pr1A from the entomopathogenic fungus Metarhizium anisopliae and VCP1 from P. chlamydosporia (86% and 76% identity, respectively). However, the catalytic pocket of P32 showed differences in the amino acids of the substrate-recognition sites compared with the catalytic pockets of Pr1A and VCP1 proteases. Phylogenetic analysis of P32 suggests a common ancestor with protease Pr1A. SCP1 displays the characteristic features of a member of the S10 family of serine proteases. Phylogenetic comparisons show that SCP1 and other carboxypeptidases from filamentous fungi have an origin different from that of yeast vacuolar serine carboxypeptidases. Understanding protease genes from nematophagous fungi is crucial for enhancing the biocontrol potential of these organisms.

CYP52X1, representing new cytochrome P450 subfamily, displays fatty acid hydroxylase activity and contributes to virulence and growth on insect cuticular substrates in entomopathogenic fungus Beauveria bassiana

Infection of insects by the entomopathogenic fungus Beauveria bassiana proceeds via attachment and penetration of the host cuticle. The outermost epicuticular layer or waxy layer of the insect represents a structure rich in lipids including abundant amounts of hydrocarbons and fatty acids. A member of a novel cytochrome P450 subfamily, CYP52X1, implicated in fatty acid assimilation by B. bassiana was characterized. B. bassiana targeted gene knockouts lacking Bbcyp52x1 displayed reduced virulence when topically applied to Galleria mellonella, but no reduction in virulence was noted when the insect cuticle was bypassed using an intrahemoceol injection assay. No significant growth defects were noted in the mutant as compared with the wild-type parent on any lipids substrates tested including alkanes and fatty acids. Insect epicuticle germination assays, however, showed reduced germination of ΔBbcyp52x1 conidia on grasshopper wings as compared with the wild-type parent. Complementation of the gene-knock with the full-length gene restored virulence and insect epicuticle germination to wild-type levels. Heterologous expression of CYP52X1 in yeast was used to characterize the substrate specificity of the enzyme. CYP52X1 displayed the highest activity against midrange fatty acids (C12:0 and C14:0) and epoxy stearic acid, 4-8-fold lower activity against C16:0, C18:1, and C18:2, and little to no activity against C9:0 and C18:0. Analyses of the products of the C12:0 and C18:1 reactions confirmed NADPH-dependent regioselective addition of a terminal hydroxyl to the substrates (ω-hydroxylase). These data implicate CYP52X1 as contributing to the penetration of the host cuticle via facilitating the assimilation of insect epicuticle lipids.

Structural basis of trypsin inhibition and entomotoxicity of cospin, serine protease inhibitor involved in defense of Coprinopsis cinerea fruiting bodies

Cospin (PIC1) from Coprinopsis cinerea is a serine protease inhibitor with biochemical properties similar to those of the previously characterized fungal serine protease inhibitors, cnispin from Clitocybe nebularis and LeSPI from Lentinus edodes, classified in the family I66 of the MEROPS protease inhibitor classification. In particular, it exhibits a highly specific inhibitory profile as a very strong inhibitor of trypsin with K(i) in the picomolar range. Determination of the crystal structure revealed that the protein has a β-trefoil fold. Site-directed mutagenesis and mass spectrometry results have confirmed Arg-27 as the reactive binding site for trypsin inhibition. The loop containing Arg-27 is positioned between the β2 and β3 strands, distinguishing cospin from other β-trefoil-fold serine protease inhibitors in which β4-β5 or β5-β6 loops are involved in protease inhibition. Biotoxicity assays of cospin on various model organisms revealed a strong and specific entomotoxic activity against Drosophila melanogaster. The inhibitory inactive R27N mutant was not entomotoxic, associating toxicity with inhibitory activity. Along with the abundance of cospin in fruiting bodies of C. cinerea and the lack of trypsin-like proteases in the C. cinerea genome, these results suggest that cospin and its homologs are effectors of a fungal defense mechanism against fungivorous insects that function by specific inhibition of serine proteases in the insect gut.

Independent subtilases expansions in fungi associated with animals

Many socially important fungi encode an elevated number of subtilisin-like serine proteases, which have been shown to be involved in fungal mutualisms with grasses and in parasitism of insects, nematodes, plants, other fungi, and mammalian skin. These proteins have endopeptidase activities and constitute a significant part of fungal secretomes. Here, we use comparative genomics to investigate the relationship between the quality and quantity of serine proteases and the ability of fungi to cause disease in invertebrate and vertebrate animals. Our screen of previously unexamined fungi allowed us to annotate and identify nearly 1000 subtilisin-containing proteins and to describe six new categories of serine proteases. Architectures of predicted proteases reveal novel combinations of subtilisin domains with other, co-occurring domains. Phylogenetic analysis of the most common clade of fungal proteases, proteinase K, showed that gene family size changed independently in fungi, pathogenic to invertebrates (Hypocreales) and vertebrates (Onygenales). Interestingly, simultaneous expansions in the S8 and S53 families of subtilases in a single fungal species are rare. Our analysis finds that closely related systemic human pathogens may not show the same gene family expansions, and that related pathogens and nonpathogens may show the same type of gene family expansion. Therefore, the number of proteases does not appear to relate to pathogenicity. Instead, we hypothesize that the number of fungal serine proteases in a species is related to the use of the animal as a food source, whether it is dead or alive.

Contribution of the gas1 gene of the entomopathogenic fungus Beauveria bassiana, encoding a putative glycosylphosphatidylinositol-anchored beta-1,3-glucanosyltransferase, to conidial thermotolerance and virulence

Beauveria bassiana is a mycoinsecticide alternative to chemicals for use in biological pest control. The fungus-insect interaction is also an emerging model system to examine unique aspects of the development, pathogenesis, and diversity of fungal lifestyles. The glycoside hydrolase 72 (GH72) family includes β-1,3-glucanosyltransferases that are glycosylphosphatidylinositol (GPI)-anchored cell wall-modeling enzymes affecting fungal physiology. A putative B. bassiana GPI-anchored β-1,3-glucanosyltransferase (Bbgas1) was isolated and characterized. B. bassiana targeted gene knockouts lacking Bbgas1 were affected in Congo red and salt sensitivity but displayed minor growth defects in the presence of sorbitol, SDS, or calcofluor white. Lectin and antibody mapping of surface carbohydrates revealed increased exposure of carbohydrate epitopes, including β-1,3-glucans, in the ΔBbgas1 strain. Transmission electron micrographs revealed localized destabilization of the cell wall in ΔBbgas1 conidia, in which fraying of the outer cell wall was apparent. Heat shock temperature sensitivity profiling showed that in contrast to the wild-type parent, ΔBbgas1 conidial spores displayed decreased germination after 1 to 4 h of heat shock at temperatures >40°C, and propidium iodide exclusion assays revealed decreased membrane stability in the knockout strain at temperatures >50°C. The ΔBbgas1 knockout showed reduced virulence in Galleria mellonella insect bioassays in both topical and intrahemocoel-injection assays. B. bassiana ΔBbgas1 strains complemented with the complete Bbgas1 open reading frame were indistinguishable from the wild-type parent in all phenotypes examined. The Bbgas1 gene did not complement the phenotype of a Saccharomyces cerevisiae β-1,3-glucanosyltransferase Δgas1 mutant, indicating that this family of enzymes likely possess discrete cellular functions.

Comparative analyses of exoproteinases produced by three phytopathogenic microorganisms

Proteinases secreted by the oomycete Phytophthora infestans (Mont.) de Bary, Rhizoctonia solani, and Fusarium culmorum belonging to different families of fungi have been studied to determine if the exoenzyme secretion depends on the environmental conditions and the phylogenetic position of the pathogen. The substrate specificity of the extracellular proteinases of F. culmorum, R. solani, and P. infestans and their sensitivity to the action of synthetic and protein inhibitors suggest that they contain trypsin-like and subtilisin-like enzymes regardless of culture medium composition. The relation of trypsin-like and subtilisin-like enzymes is dependent on the culture medium composition, especially on the form of nitrogen nutrition, particularly in the case of the exoenzymes secreted by R. solani. Phylogenetic analyses have shown that the exoproteinase set of ascomycetes and oomycetes has more similarities than basidiomycetes although they are more distant relatives. Our data suggests that the multiple proteinases secreted by pathogenic fungi could play different roles in pathogenesis, increasing the adaptability and host range, or could have different functions in survival in various ecological habitats outside the host.

The potential for pathogenicity was present in the ancestor of the Ascomycete subphylum Pezizomycotina

BACKGROUND

Previous studies in Ascomycetes have shown that the function of gene families of which the size is considerably larger in extant pathogens than in non-pathogens could be related to pathogenicity traits. However, by only comparing gene inventories in extant species, no insights can be gained into the evolutionary process that gave rise to these larger family sizes in pathogens. Moreover, most studies which consider gene families in extant species only tend to explain observed differences in gene family sizes by gains rather than by losses, hereby largely underestimating the impact of gene loss during genome evolution.

RESULTS

In our study we used a selection of recently published genomes of Ascomycetes to analyze how gene family gains, duplications and losses have affected the origin of pathogenic traits. By analyzing the evolutionary history of gene families we found that most gene families with an enlarged size in pathogens were present in an ancestor common to both pathogens and non-pathogens. The majority of these families were selectively maintained in pathogenic lineages, but disappeared in non-pathogens. Non-pathogen-specific losses largely outnumbered pathogen-specific losses.

CONCLUSIONS

We conclude that most of the proteins for pathogenicity were already present in the ancestor of the Ascomycete lineages we used in our study. Species that did not develop pathogenicity seemed to have reduced their genetic complexity compared to their ancestors. We further show that expansion of gained or already existing families in a species-specific way is important to fine-tune the specificities of the pathogenic host-fungus interaction.

New insights into the evolution of subtilisin-like serine protease genes in Pezizomycotina

Background

Subtilisin-like serine proteases play an important role in pathogenic fungi during the penetration and colonization of their hosts. In this study, we perform an evolutionary analysis of the subtilisin-like serine protease genes of subphylum Pezizomycotina to find if there are similar pathogenic mechanisms among the pathogenic fungi with different life styles, which utilize subtilisin-like serine proteases as virulence factors. Within Pezizomycotina, nematode-trapping fungi are unique because they capture soil nematodes using specialized trapping devices. Increasing evidence suggests subtilisin-like serine proteases from nematode-trapping fungi are involved in the penetration and digestion of nematode cuticles. Here we also conduct positive selection analysis on the subtilisin-like serine protease genes from nematode-trapping fungi.

Results

Phylogenetic analysis of 189 subtilisin-like serine protease genes from Pezizomycotina suggests five strongly-supported monophyletic clades. The subtilisin-like serine protease genes previously identified or presumed as endocellular proteases were clustered into one clade and diverged the earliest in the phylogeny. In addition, the cuticle-degrading protease genes from entomopathogenic and nematode-parasitic fungi were clustered together, indicating that they might have overlapping pathogenic mechanisms against insects and nematodes. Our experimental bioassays supported this conclusion. Interestingly, although they both function as cuticle-degrading proteases, the subtilisin-like serine protease genes from nematode-trapping fungi and nematode-parasitic fungi were not grouped together in the phylogenetic tree. Our evolutionary analysis revealed evidence for positive selection on the subtilisin-like serine protease genes of the nematode-trapping fungi.

Conclusions

Our study provides new insights into the evolution of subtilisin-like serine protease genes in Pezizomycotina. Pezizomycotina subtilisins most likely evolved from endocellular to extracellular proteases. The entomopathogenic and nematode-parasitic fungi likely share similar properties in parasitism. In addition, our data provided better understanding about the duplications and subsequent functional divergence of subtilisin-like serine protease genes in Pezizomycotina. The evidence of positive selection detected in the subtilisin-like serine protease genes of nematode-trapping fungi in the present study suggests that the subtilisin-like serine proteases may have played important roles during the evolution of pathogenicity of nematode-trapping fungi against nematodes.

Transcriptomic response of the mycoparasitic fungus Trichoderma atroviride to the presence of a fungal prey

Background

Combating the action of plant pathogenic microorganisms by mycoparasitic fungi has been announced as an attractive biological alternative to the use of chemical fungicides since two decades. The fungal genus Trichoderma includes a high number of taxa which are able to recognize, combat and finally besiege and kill their prey. Only fragments of the biochemical processes related to this ability have been uncovered so far, however.

Results

We analyzed genome-wide gene expression changes during the begin of physical contact between Trichoderma atroviride and two plant pathogens Botrytis cinerea and Rhizoctonia solani, and compared with gene expression patterns of mycelial and conidiating cultures, respectively. About 3000 ESTs, representing about 900 genes, were obtained from each of these three growth conditions. 66 genes, represented by 442 ESTs, were specifically and significantly overexpressed during onset of mycoparasitism, and the expression of a subset thereof was verified by expression analysis. The upregulated genes comprised 18 KOG groups, but were most abundant from the groups representing posttranslational processing, and amino acid metabolism, and included components of the stress response, reaction to nitrogen shortage, signal transduction and lipid catabolism. Metabolic network analysis confirmed the upregulation of the genes for amino acid biosynthesis and of those involved in the catabolism of lipids and aminosugars.

Conclusion

The analysis of the genes overexpressed during the onset of mycoparasitism in T. atroviride has revealed that the fungus reacts to this condition with several previously undetected physiological reactions. These data enable a new and more comprehensive interpretation of the physiology of mycoparasitism, and will aid in the selection of traits for improvement of biocontrol strains by recombinant techniques.

Characterization of the trypsin-like protease (Ha-TLP2) constitutively expressed in the integument of the cotton bollworm, Helicoverpa armigera

Trypsins belong to the serine endoproteases. They are the most important proteases in insects because of their key roles in food digestion and zymogens activation. But there has been little study of the trypsins in the integuments of insects. In this work, we cloned a trypsin-like protease gene from Helicoverpa armigera and named it trypsin-like protease 2 (Ha-TLP2). Semi-quantitative reverse transcription PCR analysis showed that Ha-TLP2 is constitutively expressed in the integument and can be down-regulated by 20-hydroxyecdysone (20E) and up-regulated by the juvenile hormone (JH) analog methoprene. Immunohistochemistry showed that Ha-TLP2 is located not only in the epidermis, but also in new and old cuticles. Immunoblotting and gelatin-SDS-PAGE revealed that Ha-TLP2 is constitutively expressed with activity in the integument during larval feeding, molting, and metamorphosis. This evidence suggests that Ha-TLP2 is involved in the remodeling of the integument.

Evolution of a subtilisin-like protease gene family in the grass endophytic fungus Epichloë festucae

BACKGROUND

Subtilisin-like proteases (SLPs) form a superfamily of enzymes that act to degrade protein substrates. In fungi, SLPs can play either a general nutritive role, or may play specific roles in cell metabolism, or as pathogenicity or virulence factors.

RESULTS

Fifteen different genes encoding SLPs were identified in the genome of the grass endophytic fungus Epichloë festucae. Phylogenetic analysis indicated that these SLPs belong to four different subtilisin families: proteinase K, kexin, pyrolysin and subtilisin. The pattern of intron loss and gain is consistent with this phylogeny. E. festucae is exceptional in that it contains two kexin-like genes. Phylogenetic analysis in Hypocreales fungi revealed an extensive history of gene loss and duplication.

CONCLUSION

This study provides new insights into the evolution of the SLP superfamily in filamentous fungi.

Comparative genomics using microarrays reveals divergence and loss of virulence-associated genes in host-specific strains of the insect pathogen Metarhizium anisopliae

Many strains of Metarhizium anisopliae have broad host ranges, but others are specialists and adapted to particular hosts. Patterns of gene duplication, divergence, and deletion in three generalist and three specialist strains were investigated by heterologous hybridization of genomic DNA to genes from the generalist strain Ma2575. As expected, major life processes are highly conserved, presumably due to purifying selection. However, up to 7% of Ma2575 genes were highly divergent or absent in specialist strains. Many of these sequences are conserved in other fungal species, suggesting that there has been rapid evolution and loss in specialist Metarhizium genomes. Some poorly hybridizing genes in specialists were functionally coordinated, indicative of reductive evolution. These included several involved in toxin biosynthesis and sugar metabolism in root exudates, suggesting that specialists are losing genes required to live in alternative hosts or as saprophytes. Several components of mobile genetic elements were also highly divergent or lost in specialists. Exceptionally, the genome of the specialist cricket pathogen Ma443 contained extra insertion elements that might play a role in generating evolutionary novelty. This study throws light on the abundance of orphans in genomes, as 15% of orphan sequences were found to be rapidly evolving in the Ma2575 lineage.

Intracellular infection of tick cell lines by the entomopathogenic fungus Metarhizium anisopliae

Several fungal pathogens are able to enter and persist within eukaryotic cells as part of their infectious life cycle. Metarhizium anisopliae is a saprophytic entomopathogenic fungus virulent towards numerous tick species, including those within the genera Ixodes and Amblyomma. Infection of the target organism by this fungus proceeds via several steps, including adhesion and penetration of the host cuticle, proliferation within tissues and the haemolymph, and eventual eruption through the host cadaver. To determine whether M. anisopliae could enter and persist within tick cells, we examined the uptake of wild-type and green fluorescent protein (GFP)-expressing fungal strains into two different tick cells lines, IDE12 and AAE2, derived from Ixodes scapularis and Amblyomma americanum, respectively. Uptake by tick cells was monitored by confocal fluorescent microscopy, as well as by scanning and transmission electron microscopy. M. anisopliae-specific antibodies were used to discriminate between extracellular and internalized fungal cells and to quantify the rate of uptake. Both cell lines internalized fungal conidia, and quantitative studies using the AAE2 cell line indicated that almost 70% of the AAE2 cells contained internalized conidia after 6 h incubation. Internalization of conidia by AAE2 cells was time and temperature dependent and could be inhibited by 80% with 1 mM cytochalasin D. Internalized conidia remained viable within the AAE2 cells, where they were able to germinate and grow, eventually erupting from the host cells. These data provide evidence that M. anisopliae conidia can be internalized, survive and grow within phagocytic cells in vitro and indicate that phagocytosis may serve as an alternative invasion route facilitating fungal nutrient acquisition, immune system evasion and dissemination throughout the host.

Comparative genome analysis of filamentous fungi reveals gene family expansions associated with fungal pathogenesis

Fungi and oomycetes are the causal agents of many of the most serious diseases of plants. Here we report a detailed comparative analysis of the genome sequences of thirty-six species of fungi and oomycetes, including seven plant pathogenic species, that aims to explore the common genetic features associated with plant disease-causing species. The predicted translational products of each genome have been clustered into groups of potential orthologues using Markov Chain Clustering and the data integrated into the e-Fungi object-oriented data warehouse (http://www.e-fungi.org.uk/). Analysis of the species distribution of members of these clusters has identified proteins that are specific to filamentous fungal species and a group of proteins found only in plant pathogens. By comparing the gene inventories of filamentous, ascomycetous phytopathogenic and free-living species of fungi, we have identified a set of gene families that appear to have expanded during the evolution of phytopathogens and may therefore serve important roles in plant disease. We have also characterised the predicted set of secreted proteins encoded by each genome and identified a set of protein families which are significantly over-represented in the secretomes of plant pathogenic fungi, including putative effector proteins that might perturb host cell biology during plant infection. The results demonstrate the potential of comparative genome analysis for exploring the evolution of eukaryotic microbial pathogenesis.