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

Transcriptome analysis highlights the influence of temperature on hydrolase and traps in nematode-trapping fungi

Pine wilt disease caused by Bursaphelenchus xylophilus poses a serious threat to the economic and ecological value of forestry. Nematode trapping fungi trap and kill nematodes using specialized trapping devices, which are highly efficient and non-toxic to the environment, and are very promising for use as biological control agents. In this study, we isolated several nematode-trapping fungi from various regions and screened three for their high nematocidal efficiency. However, the effectiveness of these fungi as nematicides is notably influenced by temperature and exhibits different morphologies in response to temperature fluctuations, which are categorized as "NA," "thin," "dense," and "sparse." The trend of trap formation with temperature was consistent with the trend of nematocidal efficiency with temperature. Both of which initially increased and then decreased with increasing temperature. Among them, Arthrobotrys cladodes exhibited the highest level of nematocidal activity and trap formation among the tested species. Transcriptome data were collected from A. cladodes with various trap morphologies. Hydrolase activity was significantly enriched according to GO and KEGG enrichment analyses. Eight genes related to hydrolases were found to be consistent with the trend of trap morphology with temperature. Weighted gene co-expression analysis and the Cytoscape network revealed that these 8 genes are associated with either mitosis or autophagy. This suggests that they contribute to the formation of "dense" structures in nematode-trapping fungi. One of these genes is the serine protein hydrolase gene involved in autophagy. This study reveals a potentially critical role for hydrolases in trap formation and nematocidal efficiency. And presents a model where temperature affects trap formation and nematocidal efficiency by influencing the serine protease prb1 involved in the autophagy process.

Expression profiles and characterization of microRNAs responding to chitin in Arthrobotrys oligospora

Extracellular proteases, such as chitinases secreted by Arthrobotrys oligospora (A. oligospora), play a crucial role in the process of nematode infection. However, post-transcriptional regulation of gene expression involving microRNAs (miRNAs) in A. oligospora remains scarcely described. Hereto, transcriptome sequencing was carried out to analyze the expression profiles of chitin-responsive miRNAs in A. oligospora. Based on the RNA-seq data, the differential expression of miRNAs (DEmiRNAs) in response to chitin was screened, identified and characterized in A. oligospora. Meanwhile, the potential target genes were predicted by the online tools miRanda and Targetscan, respectively. Furthermore, the interaction of DEmiRNA with it's target gene was validated by a dual-luciferase reporter assay system. Among 85 novel miRNAs identified, 25 miRNAs displayed significant differences in expression in A. oligospora in response to chitin. Gene Ontology (GO) analysis showed that the potential genes targeted by DEmiRNAs were enriched in the biological processes such as bio-degradation, extracellular components and cell cycle. KEGG analysis revealed that the target genes were mainly involved in Hippo, carbon and riboflavin metabolic pathway. Outstandingly, chitinase AOL_s00004g379, which is involved in the hydrolysis metabolic pathway of chitin, was confirmed to be a target gene of differential miR_70. These findings suggest that chitin-responsive miRNAs are involved in the regulation of cell proliferation, predator hyphae growth and chitinase expression through the mechanisms of post-transcriptional regulation, which provides a new perspective to the molecular mechanisms underlying miRNAs-mediated control of gene expression in A. oligospora.

Activity of root-knot nematodes associated with composition of a nematode-attached microbiome and the surrounding soil microbiota

We investigated if activity of the pre-infective juveniles (J2s) of root-knot nematodes is linked to the recruitment of a specific microbiome on the nematode surface and/or to the composition of the surrounding microbiota. For this, we determined the J2 activity (active vs. non-motile, which referred to dead and immobile J2s) upon a 3-day incubation in soil suspensions and studied the composition of bacteria, protists, and fungi present on the nematode surface and in the suspensions using amplicon sequencing of the 16S/18S rRNA genes, and ITS region. We also amended suspensions with Pseudomonas protegens strain CHA0 to study its effects on J2 activity and microbial composition. The J2 activity was suppressed in soil suspensions, but increased when suspensions were amended with P. protegens CHA0. The active and non-motile J2s differed in the composition of surface-attached bacteria, which was altered by the presence of P. protegens CHA0 in the soil suspensions. The bacterial genera Algoriphagus, Pedobacter, and Bdellovibrio were enriched on active J2s and may have protected the J2s against antagonists. The incubation time appeared short for attachment of fungi and protists. Altogether, our study is a step forward in disentangling the complex nematode-microbe interactions in soil for more successful nematode control.

Whole genome sequencing of Clarireedia aff. paspali reveals potential pathogenesis factors in Clarireedia species, causal agents of dollar spot in turfgrass

Dollar spot is one of the most damaging diseases in turfgrass, reducing its quality and playability. Two species, Clarireedia monteithiana and C. jacksonii (formerly Sclerotinia homoeocarpa) have been reported so far in the United States To study the Clarireedia genome, two isolates H2 and H3, sampled from seashore paspalum in Hawaii in 2019 were sequenced via Illumina paired-end sequencing by synthesis technology and PacBio SMRT sequencing. Both isolates were identified as C. aff. paspali, a novel species in the United States Using short and long reads, C. aff. paspali H3 contained 193 contigs with 48.6 Mbp and presented the most completed assembly and annotation among Clarireedia species. Out of the 13,428 protein models from AUGUSTUS, 349 cytoplasmic effectors and 13 apoplastic effectors were identified by EffectorP. To further decipher Clarireedia pathogenicity, C. aff. paspali genomes (H2 and H3), as well as available C. jacksonii (LWC-10 and HRI11), C. monteithiana (DRR09 and RB-19) genomes were screened for fifty-four pathogenesis determinants, previously identified in S. sclerotiorum. Seventeen orthologs of pathogenicity genes have been identified in Clarireedia species involved in oxalic acid production (pac1, nox1), mitogen-activated protein kinase cascade (pka1, smk3, ste12), appressorium formation (caf1, pks13, ams2, rgb1, rhs1) and glycolytic pathway (gpd). Within these genes, 366 species-specific SNPs were recorded between Clarireedia species; twenty-eight were non-synonymous and non-conservative. The predicted protein structure of six of these genes showed superimposition of the models among Clarireedia spp. The genomic variations revealed here could potentially lead to differences in pathogenesis and other physiological functions among Clarireedia species.

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.

New Approaches to an Old Problem: Dollar Spot of Turfgrass

Dollar spot, caused by fungal pathogens Clarireedia spp. (formerly Sclerotinia homoeocarpa), is the most common and widely distributed disease of turfgrass worldwide. It can drastically reduce the quality of turfgrass species and affect their aesthetic value and playability. Management of dollar spot typically includes a costly program of multiple application of fungicides within a growing season. Consequently, there have been reported cases of fungicide resistance in populations of Clarireedia spp. Host resistance could be an important component of dollar spot management; however, this approach has been hampered by the lack of sources of resistance because nearly all known warm- and cool-season turfgrass species are susceptible. With the recent advancement in genome sequencing technologies, studies on pathogen genomics and host-pathogen interactions are emerging with the hope of revealing candidate resistance genes in turfgrass and genes for virulence and pathogenicity in Clarireedia spp. Large-scale screening of turfgrass germplasm and quantitative trait locus (QTL) analysis for dollar spot resistance are important for resistance breeding, but only a handful of such studies have been conducted to date. This review summarizes currently available information on the dollar spot pathosystem, taxonomy, pathogen genomics, host-pathogen interaction, genetics of resistance, and QTL mapping and also provides some thoughts for future research prospects to better manage this disease.

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.

Recent advances in fungal serine protease inhibitors

Four types of antifungal drugs are available that include inhibitors of ergosterol synthesis, of fungal RNA biosynthesis, and of cell wall biosynthesis as well as physiochemical regulators of fungal membrane sterols. Increasing resistance to antifungal drugs can severely limit treatment options of fungal nail infections, vaginal candidiasis, ringworm, blastomycosis, histoplasmosis, and Candida infections of the mouth, throat, and esophagus, among other infections. Development of strategies focused on new fungicides can effectively help tackle troublesome fungal diseases. The virulence and optimal growth of fungi depend on various extracellular secreted factors, among which proteases, such as serine proteases, are of particular interest. A specific extracellular proteolytic system enables fungi to survive and penetrate the tissues. Given the role of fungal proteases in infection, any molecule capable of selectively and specifically inhibiting their activity can lead to the development of potential drugs. Owing to their specific mode of action, fungal protease inhibitors can avoid fungal resistance observed with currently available treatments. Although fungal secreted proteases have been extensively studied as potential virulence factors, our understanding of the substrate specificity of such proteases remains poor. In this review, we summarize the recent advances in the design and development of specific serine protease inhibitors and provide a brief history of the compounds that inhibit fungal serine protease activity.

Subtilisin-like Pr1 proteases marking the evolution of pathogenicity in a wide-spectrum insect-pathogenic fungus

Subtilisin-like Pr1 proteases of insect-pathogenic fungi are a large family of extracellular cuticle-degrading enzymes that presumably determine a capability of hyphal invasion into insect hemocoel through normal cuticle infection, but remain poorly understood although often considered as virulence factors for genetic improvement of fungal potential against pests. Here, we report that not all of 11 Pr1 family members necessarily function in Beauveria bassiana, an ancient wide-spectrum pathogen evolved insect pathogenicity ~200 million years ago. These Pr1 proteases are phylogenetically similar to or distinct from 11 homologues (Pr1A–K) early named in Metarhizium anisopliae complex, a young entomopathogen lineage undergoing molecular evolution toward Pr1 diversification, and hence renamed Pr1A1/A2, Pr1B1–B3, Pr1 C, Pr1F1–F4,4 and Pr1 G, respectively. Multiple analyses of all single gene-deleted and rescued mutants led to the recognition of five conserved members (Pr1C, Pr1G, Pr1A2, Pr1B1, and Pr1B2) contributing significantly to the fungal pathogenicity to insect. The conserved Pr1 proteases were proven to function only in cuticle degradation, individually contribute 19–29% to virulence, but play no role in post-infection cellular events critical for fungal killing action. Six other Pr1 proteases were not functional at all in either cuticle degradation during host infection or virulence-related cellular events post-infection. Therefore, only the five conserved proteases are collectively required for, and hence mark evolution of, insect pathogenicity in B. bassiana. These findings provide the first referable base for insight into the evolution of Pr1 family members in different lineages of fungal insect pathogens.

Metagenomics survey unravels diversity of biogas microbiomes with potential to enhance productivity in Kenya

The obstacle to optimal utilization of biogas technology is poor understanding of biogas microbiomes diversities over a wide geographical coverage. We performed random shotgun sequencing on twelve environmental samples. Randomized complete block design was utilized to assign the twelve treatments to four blocks, within eastern and central regions of Kenya. We obtained 42 million paired-end reads that were annotated against sixteen reference databases using two ENVO ontologies, prior to β-diversity studies. We identified 37 phyla, 65 classes and 132 orders. Bacteria dominated and comprised 28 phyla, 42 classes and 92 orders, conveying substrate's versatility in the treatments. Though, Fungi and Archaea comprised 5 phyla, the Fungi were richer; suggesting the importance of hydrolysis and fermentation in biogas production. High β-diversity within the taxa was largely linked to communities' metabolic capabilities. Clostridiales and Bacteroidales, the most prevalent guilds, metabolize organic macromolecules. The identified Cytophagales, Alteromonadales, Flavobacteriales, Fusobacteriales, Deferribacterales, Elusimicrobiales, Chlamydiales, Synergistales to mention but few, also catabolize macromolecules into smaller substrates to conserve energy. Furthermore, δ-Proteobacteria, Gloeobacteria and Clostridia affiliates syntrophically regulate PH2 and reduce metal to provide reducing equivalents. Methanomicrobiales and other Methanomicrobia species were the most prevalence Archaea, converting formate, CO2(g), acetate and methylated substrates into CH4(g). Thermococci, Thermoplasmata and Thermoprotei were among the sulfur and other metal reducing Archaea that contributed to redox balancing and other metabolism within treatments. Eukaryotes, mainly fungi were the least abundant guild, comprising largely Ascomycota and Basidiomycota species. Chytridiomycetes, Blastocladiomycetes and Mortierellomycetes were among the rare species, suggesting their metabolic and substrates limitations. Generally, we observed that environmental and treatment perturbations influenced communities' abundance, β-diversity and reactor performance largely through stochastic effect. Understanding diversity of biogas microbiomes over wide environmental variables and its' productivity provided insights into better management strategies that ameliorate biochemical limitations to effective biogas production.

Mechanisms of Insecticidal Action of Metarhizium anisopliae on Adult Japanese Pine Sawyer Beetles (Monochamus alternatus)

Pine wilt disease, caused by Bursaphelenchus xylophilus (pine wood nematode), leads to severe environmental and economic damage. Here, we report the results of experiments on the biological control of pine wilt disease through termination of the insect vector of the nematode and the mechanism of the insecticidal action of Metarhizium anisopliae JEF-279 against Monochamus alternatus (Japanese pine sawyer). A combined treatment with a fungal conidia suspension and a fungal protease-containing culture filtrate caused 75.8% mortality of the insect vector. Additionally, the presence of destruxins was confirmed in the dead Japanese pine sawyer adults, and half of the 10 protein spots in proteomic analysis were identified as an actin related to muscle contraction. Based on proteomic and microscopic analyses, the infection cycle of the Japanese pine sawyer by M. anisopliae JEF-279 was inferred to proceed in the following sequence: (1) host adhesion and germination, (2) epicuticle degradation, (3) growth as blastospore, (4) killing by various fungal toxins (insecticidal metabolites), (5) immune response as defense mechanism, and (6) hyphal extrusion and conidiation. Consequently, the combined fungal conidia suspension and protease-containing culture filtrate treatment may be applied as an insecticidal agent, and flaccid paralysis is likely a major mechanism underlying the insecticidal action of M. anisopliae JEF-279 on host insects.

Fungi-Nematode Interactions: Diversity, Ecology, and Biocontrol Prospects in Agriculture

Fungi and nematodes are among the most abundant organisms in soil habitats. They provide essential ecosystem services and play crucial roles for maintaining the stability of food-webs and for facilitating nutrient cycling. As two of the very abundant groups of organisms, fungi and nematodes interact with each other in multiple ways. Here in this review, we provide a broad framework of interactions between fungi and nematodes with an emphasis on those that impact crops and agriculture ecosystems. We describe the diversity and evolution of fungi that closely interact with nematodes, including food fungi for nematodes as well as fungi that feed on nematodes. Among the nematophagous fungi, those that produce specialized nematode-trapping devices are especially interesting, and a great deal is known about their diversity, evolution, and molecular mechanisms of interactions with nematodes. Some of the fungi and nematodes are significant pathogens and pests to crops. We summarize the ecological and molecular mechanisms identified so far that impact, either directly or indirectly, the interactions among phytopathogenic fungi, phytopathogenic nematodes, and crop plants. The potential applications of our understanding to controlling phytophagous nematodes and soilborne fungal pathogens in agricultural fields are discussed.

Floral transcriptomes reveal gene networks in pineapple floral growth and fruit development

Proper flower development is essential for sexual reproductive success and the setting of fruits and seeds. The availability of a high quality genome sequence for pineapple makes it an excellent model for studying fruit and floral organ development. In this study, we sequenced 27 different pineapple floral samples and integrated nine published RNA-seq datasets to generate tissue- and stage-specific transcriptomic profiles. Pairwise comparisons and weighted gene co-expression network analysis successfully identified ovule-, stamen-, petal- and fruit-specific modules as well as hub genes involved in ovule, fruit and petal development. In situ hybridization confirmed the enriched expression of six genes in developing ovules and stamens. Mutant characterization and complementation analysis revealed the important role of the subtilase gene AcSBT1.8 in petal development. This work provides an important genomic resource for functional analysis of pineapple floral organ growth and fruit development and sheds light on molecular networks underlying pineapple reproductive organ growth.

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.

The Subtilisin-Like Protease Bcser2 Affects the Sclerotial Formation, Conidiation and Virulence of Botrytis cinerea

Botrytis cinerea, a ubiquitous necrotrophic plant-pathogenic fungus, is responsible for grey mold and rot disease in a very wide range of plant species. Subtilisin-like proteases (or subtilases) are a very diverse family of serine proteases present in many organisms and are reported to have a broad spectrum of biological functions. Here, we identified two genes encoding subtilisin-like proteases (Bcser1 and Bcser2) in the genome of B. cinerea, both of which contain an inhibitor I9 domain and a peptidase S8 domain. The expression levels of Bcser1 and Bcser2 increased during the sclerotial forming stage, as well as during a later stage of hyphal infection on Arabidopsis thaliana leaves, but the up-regulation of Bcser1 was significantly higher than that of Bcser2. Interestingly, deletion of Bcser1 had no effect on the fungal development or virulence of B. cinerea. However, deletion of Bcser2 or double deletion of Bcser1 and Bcser2 severely impaired the hyphal growth, sclerotial formation and conidiation of B. cinerea. We also found that ∆Bcser2 and ∆Bcser1/2 could not form complete infection cushions and then lost the ability to infect intact plant leaves of Arabidopsis and tomato but could infect wounded plant tissues. Taken together, our results indicate that the subtilisin-like protease Bcser2 is crucial for the sclerotial formation, conidiation, and virulence of B. cinerea.

Metarhizium robertsii ammonium permeases (MepC and Mep2) contribute to rhizoplane colonization and modulates the transfer of insect derived nitrogen to plants

The endophytic insect pathogenic fungi (EIPF) Metarhizium promotes plant growth through symbiotic association and the transfer of insect-derived nitrogen. However, little is known about the genes involved in this association and the transfer of nitrogen. In this study, we assessed the involvement of six Metarhizium robertsii genes in endophytic, rhizoplane and rhizospheric colonization with barley roots. Two ammonium permeases (MepC and Mep2) and a urease, were selected since homologous genes in arbuscular mycorrhizal fungi were reported to play a pivotal role in nitrogen mobilization during plant root colonization. Three other genes were selected on the basis on RNA-Seq data that showed high expression levels on bean roots, and these encoded a hydrophobin (Hyd3), a subtilisin-like serine protease (Pr1A) and a hypothetical protein. The root colonization assays revealed that the deletion of urease, hydrophobin, subtilisin-like serine protease and hypothetical protein genes had no impact on endophytic, rhizoplane and rhizospheric colonization at 10 or 20 days. However, the deletion of MepC resulted in significantly increased rhizoplane colonization at 10 days whereas ΔMep2 showed increased rhizoplane colonization at 20 days. In addition, the nitrogen transporter mutants also showed significantly higher 15N incorporation of insect derived nitrogen in barley leaves in the presence of nutrients. Insect pathogenesis assay revealed that disruption of MepC, Mep2, urease did not reduce virulence toward insects. The enhanced rhizoplane colonization of ΔMep2 and ΔMepC and insect derived nitrogen transfer to plant hosts suggests the role of MepC and Mep2 in Metarhizium-plant symbiosis.

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.

Widespread signatures of selection for secreted peptidases in a fungal plant pathogen

BACKGROUND

Fungal plant pathogens secrete a large arsenal of hydrolytic enzymes during the course of infection, including peptidases. Secreted peptidases have been extensively studied for their role as effectors. In this study, we combined transcriptomics, comparative genomics and evolutionary analyses to investigate all 39 secreted peptidases in the fungal wheat pathogen Zymoseptoria tritici and its close relatives Z. pseudotritici and Z. ardabiliae.

RESULTS

RNA-seq data revealed that a majority of the secreted peptidases displayed differential transcription during the course of Z. tritici infection, indicative of specialization for different stages in the life cycle. Evolutionary analyses detected widespread evidence of adaptive evolution acting on at least 28 of the peptidases. A few peptidases displayed lineage-specific rates of molecular evolution, suggesting altered selection pressure in Z. tritici following host specialization on domesticated wheat. The peptidases belonging to MEROPS families A1 and G1 emerged as a particularly interesting group that may play key roles in host-pathogen co-evolution, host adaptation and pathogenicity. Sister genes in the A1 and G1 families showed accelerated substitution rates after gene duplications.

CONCLUSIONS

These results suggest widespread evolution of secreted peptidases leading to novel gene functions, consistent with predicted models of "escape from adaptive conflict" and "neo-functionalization". Our analyses identified candidate genes worthy of functional analyses that may encode effector functions, for example by suppressing plant defenses during the biotrophic phase of infection.

Comparative genomic analyses reveal the features for adaptation to nematodes in fungi

Nematophagous (NP) fungi are ecologically important components of the soil microbiome in natural ecosystems. Esteya vermicola (Ev) has been reported as a NP fungus with a poorly understood evolutionary history and mechanism of adaptation to parasitism. Furthermore, NP fungal genomic basis of lifestyle was still unclear. We sequenced and annotated the Ev genome (34.2 Mbp) and integrated genetic makeup and evolution of pathogenic genes to investigate NP fungi. The results revealed that NP fungi had some abundant pathogenic genes corresponding to their niche. A number of gene families involved in pathogenicity were expanded, and some pathogenic orthologous genes underwent positive selection. NP fungi with diverse morphological features exhibit similarities of evolutionary convergence in attacking nematodes, but their genetic makeup and microscopic mechanism are different. Endoparasitic NP fungi showed similarity in large number of transporters and secondary metabolite coding genes. Noteworthy, expanded families of transporters and endo-beta-glucanase implied great genetic potential of Ev in quickly perturbing nematode metabolism and parasitic behavior. These results facilitate our understanding of NP fungal genomic features for adaptation to nematodes and lay a solid theoretical foundation for further research and application.

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.

Knowledge-transfer learning for prediction of matrix metalloprotease substrate-cleavage sites

Matrix Metalloproteases (MMPs) are an important family of proteases that play crucial roles in key cellular and disease processes. Therefore, MMPs constitute important targets for drug design, development and delivery. Advanced proteomic technologies have identified type-specific target substrates; however, the complete repertoire of MMP substrates remains uncharacterized. Indeed, computational prediction of substrate-cleavage sites associated with MMPs is a challenging problem. This holds especially true when considering MMPs with few experimentally verified cleavage sites, such as for MMP-2, -3, -7, and -8. To fill this gap, we propose a new knowledge-transfer computational framework which effectively utilizes the hidden shared knowledge from some MMP types to enhance predictions of other, distinct target substrate-cleavage sites. Our computational framework uses support vector machines combined with transfer machine learning and feature selection. To demonstrate the value of the model, we extracted a variety of substrate sequence-derived features and compared the performance of our method using both 5-fold cross-validation and independent tests. The results show that our transfer-learning-based method provides a robust performance, which is at least comparable to traditional feature-selection methods for prediction of MMP-2, -3, -7, -8, -9 and -12 substrate-cleavage sites on independent tests. The results also demonstrate that our proposed computational framework provides a useful alternative for the characterization of sequence-level determinants of MMP-substrate specificity.

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.

Phylogenetic analyses reveal molecular signatures associated with functional divergence among Subtilisin like Serine Proteases are linked to lifestyle transitions in Hypocreales

Background

Subtilisin-like serine proteases or Subtilases in fungi are important for penetration and colonization of host. In Hypocreales, these proteins share several properties with other fungal, bacterial, plant and mammalian homologs. However, adoption of specific roles in entomopathogenesis may be governed by attainment of unique biochemical and structural features during the evolutionary course. Due to such functional shifts Subtilases coded by different family members of Hypocreales acquire distinct features according to respective hosts and lifestyle. We conducted phylogenetic and DIVERGE analyses and identified important protein residues that putatively assign functional specificity to Subtilases in fungal families/species under the order Hypocreales.

Results

A total of 161 Subtilases coded by 10 species from five different families under the fungal order Hypocreales was included in the analysis. Based on the presence of conserved domains, the Subtilase genes were divided into three subfamilies, Subtilisin (S08.005), Proteinase K (S08.054) and Serine-carboxyl peptidases (S53.001). These subfamilies were investigated for phylogenetic associations, protein residues under positive selection and functional divergence among paralogous clades. The observations were co-related with the life-styles of the fungal families/species. Phylogenetic and Divergence analyses of Subtilisin (S08.005) and Proteinase K (S08.054) families of proteins revealed that the paralogous clades were clear-cut representation of familial origin of the protein sequences. We observed divergence between the paralogous clades of plant-pathogenic fungi (Nectriaceae), insect-pathogenic fungi (Cordycipitaceae/Clavicipitaceae) and nematophagous fungi (Ophiocordycipitaceae). In addition, Subtilase genes from the nematode-parasitic fungus Purpureocillium lilacinum made a unique cluster which putatively indicated that the fungus might have developed distinctive mechanisms for nematode-pathogenesis. Our evolutionary genetics analysis revealed evidence of positive selection on the Subtilisin (S08.005) and Proteinase K (S08.054) protein sequences of the entomopathogenic and nematophagous species belonging to Cordycipitaceae, Clavicipitaceae and Ophiocordycipitaceae families of Hypocreales.

Conclusions

Our study provided new insights into the evolution of Subtilisin like serine proteases in Hypocreales, a fungal order largely consisting of biological control species. Subtilisin (S08.005) and Proteinase K (S08.054) proteins seemed to play important roles during life style modifications among different families and species of Hypocreales. Protein residues found significant in functional divergence analysis in the present study may provide support for protein engineering in future.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0793-y) contains supplementary material, which is available to authorized users.

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.

Whole genome annotation and comparative genomic analyses of bio-control fungus Purpureocillium lilacinum

Background

The fungus Purpureocillium lilacinum is widely known as a biological control agent against plant parasitic nematodes. This research article consists of genomic annotation of the first draft of whole genome sequence of P. lilacinum. The study aims to decipher the putative genetic components of the fungus involved in nematode pathogenesis by performing comparative genomic analysis with nine closely related fungal species in Hypocreales.

Results

de novo genomic assembly was done and a total of 301 scaffolds were constructed for P. lilacinum genomic DNA. By employing structural genome prediction models, 13, 266 genes coding for proteins were predicted in the genome. Approximately 73 % of the predicted genes were functionally annotated using Blastp, InterProScan and Gene Ontology. A 14.7 % fraction of the predicted genes shared significant homology with genes in the Pathogen Host Interactions (PHI) database. The phylogenomic analysis carried out using maximum likelihood RAxML algorithm provided insight into the evolutionary relationship of P. lilacinum. In congruence with other closely related species in the Hypocreales namely, Metarhizium spp., Pochonia chlamydosporia, Cordyceps militaris, Trichoderma reesei and Fusarium spp., P. lilacinum has large gene sets coding for G-protein coupled receptors (GPCRs), proteases, glycoside hydrolases and carbohydrate esterases that are required for degradation of nematode-egg shell components. Screening of the genome by Antibiotics & Secondary Metabolite Analysis Shell (AntiSMASH) pipeline indicated that the genome potentially codes for a variety of secondary metabolites, possibly required for adaptation to heterogeneous lifestyles reported for P. lilacinum. Significant up-regulation of subtilisin-like serine protease genes in presence of nematode eggs in quantitative real-time analyses suggested potential role of serine proteases in nematode pathogenesis.

Conclusions

The data offer a better understanding of Purpureocillium lilacinum genome and will enhance our understanding on the molecular mechanism involved in nematophagy.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2229-2) contains supplementary material, which is available to authorized users.

Expression of serine proteinase P186 of Arthrobotrys oligospora and analysis of its nematode-degrading activity

The nematode-trapping fungi possess a unique capability of predating and invading nematodes. As a representative nematode-trapping fungus, Arthrobotrys oligospora has been widely used to study the interactions between nematode-trapping fungi and their hosts. Serine proteinase is one of the important virulence factors during process of invasion of the nematode-trapping fungi into nematodes. In this study, using reverse transcription polymerase chain reaction, we amplified the gene sequence of serine proteinase 186 from A. oligospora, cloned it into pPIC9K vector and expressed it in the yeast Pichia pastoris. The expressed recombinant serine proteinase186 (reP186) was purified via Ni-affinity chromatography. The in vitro nematode-degrading activity of reP186 was analyzed. Sodium dodecyl sulfate polyacrylamide gel electrophoresis and Western blot analysis revealed that reP186 with molecular weight of 33 kDa was successfully obtained. ReP186 was capable of degrading a series of protein substrates including casein, gelatin, bovine serum albumin, denatured collagen and nematode cortical layer. The reP186 exhibited the maximal activity at pH 8.0 and 55 °C and was highly sensitive to the inhibitor, phenylmethanesulfonylfluoride. Treatment of Caenorhabditis elegans and Haemonchus contortus with reP186 for 12, 24 and 36 h, respectively, resulted in 62, 88 and 100 % of killing rates for C. elegans, and 52, 65 and 84 % of killing rates for H. contortus, respectively, indicating a relatively strong nematode-degrading bioactivity of reP186.

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.

Destructin-1 is a collagen-degrading endopeptidase secreted by Pseudogymnoascus destructans, the causative agent of white-nose syndrome

Pseudogymnoascus destructans is the causative agent of white-nose syndrome, a disease that has caused the deaths of millions of bats in North America. This psychrophilic fungus proliferates at low temperatures and targets hibernating bats, resulting in their premature arousal from stupor with catastrophic consequences. Despite the impact of white-nose syndrome, little is known about the fungus itself or how it infects its mammalian host. P. destructans is not amenable to genetic manipulation, and therefore understanding the proteins involved in infection requires alternative approaches. Here, we identify hydrolytic enzymes secreted by P. destructans, and use a novel and unbiased substrate profiling technique to define active peptidases. These experiments revealed that endopeptidases are the major proteolytic activities secreted by P. destructans, and that collagen, the major structural protein in mammals, is actively degraded by the secretome. A serine endopeptidase, hereby-named Destructin-1, was subsequently identified, and a recombinant form overexpressed and purified. Biochemical analysis of Destructin-1 showed that it mediated collagen degradation, and a potent inhibitor of peptidase activity was identified. Treatment of P. destructans-conditioned media with this antagonist blocked collagen degradation and facilitated the detection of additional secreted proteolytic activities, including aminopeptidases and carboxypeptidases. These results provide molecular insights into the secretome of P. destructans, and identify serine endopeptidases that have the clear potential to facilitate tissue invasion and pathogenesis in the mammalian host.

Molecular mechanisms of nematode-nematophagous microbe interactions: basis for biological control of plant-parasitic nematodes

Plant-parasitic nematodes cause significant damage to a broad range of vegetables and agricultural crops throughout the world. As the natural enemies of nematodes, nematophagous microorganisms offer a promising approach to control the nematode pests. Some of these microorganisms produce traps to capture and kill the worms from the outside. Others act as internal parasites to produce toxins and virulence factors to kill the nematodes from within. Understanding the molecular basis of microbe-nematode interactions provides crucial insights for developing effective biological control agents against plant-parasitic nematodes. Here, we review recent advances in our understanding of the interactions between nematodes and nematophagous microorganisms, with a focus on the molecular mechanisms by which nematophagous microorganisms infect nematodes and on the nematode defense against pathogenic attacks. We conclude by discussing several key areas for future research and development, including potential approaches to apply our recent understandings to develop effective biocontrol strategies.

Cloning of a serine protease gene from the nematophagous fungus Esteya vermicola and expressed activity of the recombinant enzyme against Bursaphelenchus xylophilus

A serine protease geneEvspwas cloned from the nematophagous fungusEsteya vermicolawith strong virulence againstBursaphelenchus xylophilus. The full-length cDNA ofEvspcontains 2280 nucleotides with a 1656 bp ORF encoding a protein with 551 amino acids. The genomicEvspincludes two exons (396 bp and 1260 bp) separated by an intron (207 bp). There is only one copy ofEvspgene in the fungal genome. The deduced amino acids sequences ofEvspshowed highly homology with the catalytic domains in subtilisin serine proteases. Phylogenetic analyses based on the protein sequences revealed thatE. vermicolais separated from nematode-trapping fungi but close to other nematophagous and entomopathogenic fungi. The recombinant serine protease rEvsp was induced inEscherichia coliwith expression vector pET28a(+). The tests of protease and nematicidal activities for the purified and refolded rEvsp indicated it is possibly involved in the fungal infection process againstB. xylophilus.

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.

Enzymology of the nematode cuticle: A potential drug target?

All nematodes possess an external structure known as the cuticle, which is crucial for their development and survival. This structure is composed primarily of collagen, which is secreted from the underlying hypodermal cells. Extensive studies using the free-living nematode Caenorhabditis elegans demonstrate that formation of the cuticle requires the activity of an extensive range of enzymes. Enzymes are required both pre-secretion, for synthesis of component proteins such as collagen, and post-secretion, for removal of the previous developmental stage cuticle, in a process known as moulting or exsheathment. The excretion/secretion products of numerous parasitic nematodes contain metallo-, serine and cysteine proteases, and these proteases are conserved across the nematode phylum and many are involved in the moulting/exsheathment process. This review highlights the enzymes required for cuticle formation, with a focus on the post-secretion moulting events. Where orthologues of the C. elegans enzymes have been identified in parasitic nematodes these may represent novel candidate targets for future drug/vaccine development.

Leucoagaricus gongylophorus uses leaf-cutting ants to vector proteolytic enzymes towards new plant substrate

The mutualism between leaf-cutting ants and their fungal symbionts revolves around processing and inoculation of fresh leaf pulp in underground fungus gardens, mediated by ant fecal fluid deposited on the newly added plant substrate. As herbivorous feeding often implies that growth is nitrogen limited, we cloned and sequenced six fungal proteases found in the fecal fluid of the leaf-cutting ant Acromyrmex echinatior and identified them as two metalloendoproteases, two serine proteases and two aspartic proteases. The metalloendoproteases and serine proteases showed significant activity in fecal fluid at pH values of 5-7, but the aspartic proteases were inactive across a pH range of 3-10. Protease activity disappeared when the ants were kept on a sugar water diet without fungus. Relative to normal mycelium, both metalloendoproteases, both serine proteases and one aspartic protease were upregulated in the gongylidia, specialized hyphal tips whose only known function is to provide food to the ants. These combined results indicate that the enzymes are derived from the ingested fungal tissues. We infer that the five proteases are likely to accelerate protein extraction from plant cells in the leaf pulp that the ants add to the fungus garden, but regulatory functions such as activation of proenzymes are also possible, particularly for the aspartic proteases that were present but without showing activity. The proteases had high sequence similarities to proteolytic enzymes of phytopathogenic fungi, consistent with previous indications of convergent evolution of decomposition enzymes in attine ant fungal symbionts and phytopathogenic fungi.

Characterization and functional analyses of the chitinase-encoding genes in the nematode-trapping fungus Arthrobotrys oligospora

Nematode-trapping fungi can secrete many extracellular hydrolytic enzymes such as serine proteases and chitinases to digest and penetrate nematode/egg-cuticles. However, little is known about the structure and function of chitinases in these fungi. In this study, 16 ORFs encoding putative chitinases, which all belong to glycoside hydrolase (GH) family 18, were identified from the Arthrobotrys oligospora genome. Bioinformatics analyses showed that these 16 putative chitinases differ in their functional domains, molecular weights and pI. Phylogenetic analysis grouped these A. oligospora chitinases into four clades: clades I, II, III and IV, respectively, including an A. oligospora-specific subclade (Clade IV-B) that contained high-molecular weight chitinases (≥100 kDa). Transcriptional analysis of A. oligospora chitinases suggested that the expression of most chitinases was repressed by carbon starvation, and all chitinases were up-regulated under nitrogen starvation. However, chitinase AO-190 was up-regulated under carbon and/or nitrogen starvation. Moreover, several chitinases (such as AO-59, AO-190 and AO-801) were up-regulated in the presence of chitinous substrates or a plant pathogenic fungus, indicating that they could play a role in biocontrol applications of A. oligospora. Our results provided a basis for further understanding the functions, diversities and evolutionary relationships between chitinase genes in nematode-trapping fungi.

RNA-Seq analysis of the Sclerotinia homoeocarpa--creeping bentgrass pathosystem

Sclerotinia homoeocarpa causes dollar spot disease, the predominate disease on highly-maintained turfgrass. Currently, there are major gaps in our understanding of the molecular interactions between S. homoeocarpa and creeping bentgrass. In this study, 454 sequencing technology was used in the de novo assembly of S. homoeocarpa and creeping bentgrass transcriptomes. Transcript sequence data obtained using Illumina's first generation sequencing-by-synthesis (SBS) were mapped to the transcriptome assemblies to estimate transcript representation in different SBS libraries. SBS libraries included a S. homoeocarpa culture control, a creeping bentgrass uninoculated control, and a library for creeping bentgrass inoculated with S. homoeocarpa and incubated for 96 h. A Fisher's exact test was performed to determine transcripts that were significantly different during creeping bentgrass infection with S. homoeocarpa. Fungal transcripts of interest included glycosyl hydrolases, proteases, and ABC transporters. Of particular interest were the large number of glycosyl hydrolase transcripts that target a wide range of plant cell wall compounds, corroborating the suggested wide host range and saprophytic abilities of S. homoeocarpa. Several of the multidrug resistance ABC transporters may be important for resistance to both fungicides and plant defense compounds. Creeping bentgrass transcripts of interest included germins, ubiquitin transcripts involved in proteasome degradation, and cinnamoyl reductase, which is involved in lignin production. This analysis provides an extensive overview of the S. homoeocarpa-turfgrass pathosystem and provides a starting point for the characterization of potential virulence factors and host defense responses. In particular, determination of important host defense responses may assist in the development of highly resistant creeping bentgrass varieties.

Origin and evolution of carnivorism in the Ascomycota (fungi)

Carnivorism is one of the basic life strategies of fungi. Carnivorous fungi possess the ability to trap and digest their preys by sophisticated trapping devices. However, the origin and development of fungal carnivorism remains a gap in evolution biology. In this study, five protein-encoding genes were used to construct the phylogeny of the carnivorous fungi in the phylum Ascomycota; these fungi prey on nematodes by means of specialized trapping structures such as constricting rings and adhesive traps. Our analysis revealed a definitive pattern of evolutionary development for these trapping structures. Molecular clock calibration based on two fossil records revealed that fungal carnivorism diverged from saprophytism about 419 Mya, which was after the origin of nematodes about 550-600 Mya. Active carnivorism (fungi with constricting rings) and passive carnivorism (fungi with adhesive traps) diverged from each other around 246 Mya, shortly after the occurrence of the Permian-Triassic extinction event about 251.4 Mya. The major adhesive traps evolved around 198-208 Mya, which was within the time frame of the Triassic-Jurassic extinction event about 201.4 Mya. However, no major carnivorous ascomycetes divergence was correlated to the Cretaceous-Tertiary extinction event, which occurred more recently (about 65.5 Mya). Therefore, a causal relationship between mass extinction events and fungal carnivorism evolution is not validated in this study. More evidence including additional fossil records is needed to establish if fungal carnivorism evolution was a response to mass extinction events.

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.

Genomic and proteomic analyses of the fungus Arthrobotrys oligospora provide insights into nematode-trap formation

Nematode-trapping fungi are "carnivorous" and attack their hosts using specialized trapping devices. The morphological development of these traps is the key indicator of their switch from saprophytic to predacious lifestyles. Here, the genome of the nematode-trapping fungus Arthrobotrys oligospora Fres. (ATCC24927) was reported. The genome contains 40.07 Mb assembled sequence with 11,479 predicted genes. Comparative analysis showed that A. oligospora shared many more genes with pathogenic fungi than with non-pathogenic fungi. Specifically, compared to several sequenced ascomycete fungi, the A. oligospora genome has a larger number of pathogenicity-related genes in the subtilisin, cellulase, cellobiohydrolase, and pectinesterase gene families. Searching against the pathogen-host interaction gene database identified 398 homologous genes involved in pathogenicity in other fungi. The analysis of repetitive sequences provided evidence for repeat-induced point mutations in A. oligospora. Proteomic and quantitative PCR (qPCR) analyses revealed that 90 genes were significantly up-regulated at the early stage of trap-formation by nematode extracts and most of these genes were involved in translation, amino acid metabolism, carbohydrate metabolism, cell wall and membrane biogenesis. Based on the combined genomic, proteomic and qPCR data, a model for the formation of nematode trapping device in this fungus was proposed. In this model, multiple fungal signal transduction pathways are activated by its nematode prey to further regulate downstream genes associated with diverse cellular processes such as energy metabolism, biosynthesis of the cell wall and adhesive proteins, cell division, glycerol accumulation and peroxisome biogenesis. This study will facilitate the identification of pathogenicity-related genes and provide a broad foundation for understanding the molecular and evolutionary mechanisms underlying fungi-nematodes interactions.

Comparison of homology models and crystal structures of cuticle-degrading proteases from nematophagous fungi: structural basis of nematicidal activity

Cuticle-degrading proteases secreted by nematophagous fungi can degrade nematode cuticle during infection. Alkaline proteases from nematode-parasitic fungi show stronger nematicidal activity in vitro than neutral proteases from nematode-trapping fungi. Sequence alignment of these proteases revealed that the active-site residues were much conserved. Disulfide bridges in alkaline proteases not only contribute to the thermal stability of enzyme structure but also increase the flexibility of S1 and S4 pockets located at the substrate-binding site. Molecular electrostatic potential surfaces of these proteases change gradually from negative to positive while arranging in the order from neutral to alkaline proteases, possibly contributing to the distinct extent of substrate (nematode cuticle) attraction by proteases. The differences in flexibility of substrate-binding site and in electrostatic surface potential distribution between neutral and alkaline cuticle-degrading proteases are associated with the changes of their catalytic activities and nematicidal activities with fungal species. Our results indicate that nematode-parasitic and nematode-trapping fungi have evolved for distinct adaptation under selective pressure.

Overexpression of a cuticle-degrading protease Ver112 increases the nematicidal activity of Paecilomyces lilacinus

Due to their ability to degrade the proteins in nematode cuticle, serine proteases play an important role in the pathogenicity of nematophagous fungi against nematodes. The serine protease Ver112 was identified from the nematophagous fungus Lecanicillium psalliotae capable of degrading the nematode cuticle and killing nematodes effectively. In this study, the gene ver112 was introduced into the commercial biocontrol fungal agent Paecilomyces lilacinus by the restriction enzyme-mediated integration transformation. Compared to the wild strain, the transformant P. lilacinus 112 showed significantly greater protease activity, with nematicidal activities increased by 79% and 96% to Panagrellus redivivus and Caenorhabditis elegans at the second day, respectively. The crude protein extract isolated from the culture filtrate of P. lilacinus 112 also showed 20-25% higher nematicidal activity than that of the wild-type strain. Reverse transcription PCR results showed that the expression of gene ver112 in P. lilacinus 112 was correlated to protease activity of the culture filtrate. Our results demonstrated the first successful transfer of a virulence gene from one nematophagous fungus to another nematophagous fungus, and improved the pathogenicity of the recipient fungus against pest nematodes.