Xinyu Zhang

Trajectory and genomic determinants of fungal-pathogen speciation and host adaptation

Significance Fossil records have provided compelling evidence for evolution, but lack of existing transitional species has hindered our understanding of speciation at the molecular level. Genomic analyses of seven Metarhizium species revealed a directional speciation continuum from specialists with narrow host ranges to transitional species and then to generalists that paralleled insect evolution. This diversification was coupled with a complex interplay between an array of genomic features that worked together to drive fungal speciation at an accelerating rate and provided a roadmap for identifying variation underlying adaptation and speciation. In particular, specialization was associated with retention of sexuality and rapid evolution of existing protein sequences whereas generalization was associated with loss of sexuality and protein-family expansion. Much remains unknown regarding speciation. Host–pathogen interactions are a major driving force for diversification, but the genomic basis for speciation and host shifting remains unclear. The fungal genus Metarhizium contains species ranging from specialists with very narrow host ranges to generalists that attack a wide range of insects. By genomic analyses of seven species, we demonstrated that generalists evolved from specialists via transitional species with intermediate host ranges and that this shift paralleled insect evolution. We found that specialization was associated with retention of sexuality and rapid evolution of existing protein sequences whereas generalization was associated with protein-family expansion, loss of genome-defense mechanisms, genome restructuring, horizontal gene transfer, and positive selection that accelerated after reinforcement of reproductive isolation. These results advance understanding of speciation and genomic signatures that underlie pathogen adaptation to hosts.

Genomic and transcriptomic analysis of the endophytic fungus Pestalotiopsis fici reveals its lifestyle and high potential for synthesis of natural products

BackgroundIn recent years, the genus Pestalotiopsis is receiving increasing attention, not only because of its economic impact as a plant pathogen but also as a commonly isolated endophyte which is an important source of bioactive natural products. Pestalotiopsis fici Steyaert W106-1/CGMCC3.15140 as an endophyte of tea produces numerous novel secondary metabolites, including chloropupukeananin, a derivative of chlorinated pupukeanane that is first discovered in fungi. Some of them might be important as the drug leads for future pharmaceutics.ResultsHere, we report the genome sequence of the endophytic fungus of tea Pestalotiopsis fici W106-1/CGMCC3.15140. The abundant carbohydrate-active enzymes especially significantly expanding pectinases allow the fungus to utilize the limited intercellular nutrients within the host plants, suggesting adaptation of the fungus to endophytic lifestyle. The P. fici genome encodes a rich set of secondary metabolite synthesis genes, including 27 polyketide synthases (PKSs), 12 non-ribosomal peptide synthases (NRPSs), five dimethylallyl tryptophan synthases, four putative PKS-like enzymes, 15 putative NRPS-like enzymes, 15 terpenoid synthases, seven terpenoid cyclases, seven fatty-acid synthases, and five hybrids of PKS-NRPS. The majority of these core enzymes distributed into 74 secondary metabolite clusters. The putative Diels-Alderase genes have undergone expansion.ConclusionThe significant expansion of pectinase encoding genes provides essential insight in the life strategy of endophytes, and richness of gene clusters for secondary metabolites reveals high potential of natural products of endophytic fungi.

Genomics-driven discovery of the pneumocandin biosynthetic gene cluster in the fungus Glarea lozoyensis

BackgroundThe antifungal therapy caspofungin is a semi-synthetic derivative of pneumocandin B0, a lipohexapeptide produced by the fungus Glarea lozoyensis, and was the first member of the echinocandin class approved for human therapy. The nonribosomal peptide synthetase (NRPS)-polyketide synthases (PKS) gene cluster responsible for pneumocandin biosynthesis from G. lozoyensis has not been elucidated to date. In this study, we report the elucidation of the pneumocandin biosynthetic gene cluster by whole genome sequencing of the G. lozoyensis wild-type strain ATCC 20868.ResultsThe pneumocandin biosynthetic gene cluster contains a NRPS (GLNRPS4) and a PKS (GLPKS4) arranged in tandem, two cytochrome P450 monooxygenases, seven other modifying enzymes, and genes for L-homotyrosine biosynthesis, a component of the peptide core. Thus, the pneumocandin biosynthetic gene cluster is significantly more autonomous and organized than that of the recently characterized echinocandin B gene cluster. Disruption mutants of GLNRPS4 and GLPKS4 no longer produced the pneumocandins (A0 and B0), and the Δglnrps4 and Δglpks4 mutants lost antifungal activity against the human pathogenic fungus Candida albicans. In addition to pneumocandins, the G. lozoyensis genome encodes a rich repertoire of natural product-encoding genes including 24 PKSs, six NRPSs, five PKS-NRPS hybrids, two dimethylallyl tryptophan synthases, and 14 terpene synthases.ConclusionsCharacterization of the gene cluster provides a blueprint for engineering new pneumocandin derivatives with improved pharmacological properties. Whole genome estimation of the secondary metabolite-encoding genes from G. lozoyensis provides yet another example of the huge potential for drug discovery from natural products from the fungal kingdom.

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.

Genome characteristics reveal the impact of lichenization on lichen-forming fungus Endocarpon pusillum Hedwig (Verrucariales, Ascomycota)

BackgroundLichen is a classic mutualistic organism and the lichenization is one of the fungal symbioses. The lichen-forming fungus Endocarpon pusillum is living in symbiosis with the green alga Diplosphaera chodatii Bialsuknia as a lichen in the arid regions.Results454 and Illumina technologies were used to sequence the genome of E. pusillum. A total of 9,285 genes were annotated in the 37.5xa0Mb genome of E. pusillum. Analyses of the genes provided direct molecular evidence for certain natural characteristics, such as homothallic reproduction and drought-tolerance. Comparative genomics analysis indicated that the expansion and contraction of some protein families in the E. pusillum genome reflect the specific relationship with its photosynthetic partner (D. chodatii). Co-culture experiments using the lichen-forming fungus E. pusillum and its algal partner allowed the functional identification of genes involved in the nitrogen and carbon transfer between both symbionts, and three lectins without signal peptide domains were found to be essential for the symbiotic recognition in the lichen; interestingly, the ratio of the biomass of both lichen-forming fungus and its photosynthetic partner and their contact time were found to be important for the interaction between these two symbionts.ConclusionsThe present study lays a genomic analysis of the lichen-forming fungus E. pusillum for demonstrating its general biological features and the traits of the interaction between this fungus and its photosynthetic partner D. chodatii, and will provide research basis for investigating the nature of its drought resistance and symbiosis.

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.

Drechslerella stenobrocha genome illustrates the mechanism of constricting rings and the origin of nematode predation in fungi

BackgroundNematode-trapping fungi are a unique group of organisms that can capture nematodes using sophisticated trapping structures. The genome of Drechslerella stenobrocha, a constricting-ring-forming fungus, has been sequenced and reported, and provided new insights into the evolutionary origins of nematode predation in fungi, the trapping mechanisms, and the dual lifestyles of saprophagy and predation.ResultsThe genome of the fungus Drechslerella stenobrocha, which mechanically traps nematodes using a constricting ring, was sequenced. The genome was 29.02xa0Mb in size and was found rare instances of transposons and repeat induced point mutations, than that of Arthrobotrys oligospora. The functional proteins involved in nematode-infection, such as chitinases, subtilisins, and adhesive proteins, underwent a significant expansion in the A. oligospora genome, while there were fewer lectin genes that mediate fungus-nematode recognition in the D. stenobrocha genome. The carbohydrate-degrading enzyme catalogs in both species were similar to those of efficient cellulolytic fungi, suggesting a saprophytic origin of nematode-trapping fungi. In D. stenobrocha, the down-regulation of saprophytic enzyme genes and the up-regulation of infection-related genes during the capture of nematodes indicated a transition between dual life strategies of saprophagy and predation. The transcriptional profiles also indicated that trap formation was related to the protein kinase C (PKC) signal pathway and regulated by Zn(2)–C6 type transcription factors.ConclusionsThe genome of D. stenobrocha provides support for the hypothesis that nematode trapping fungi evolved from saprophytic fungi in a high carbon and low nitrogen environment. It reveals the transition between saprophagy and predation of these fungi and also proves new insights into the mechanisms of mechanical trapping.

Comparative genomic and transcriptomic analyses of the Fuzhuan brick tea-fermentation fungus Aspergillus cristatus.

BackgroundAspergillus cristatus is the dominant fungus involved in the fermentation of Chinese Fuzhuan brick tea. Aspergillus cristatus is a homothallic fungus that undergoes a sexual stage without asexual conidiation when cultured in hypotonic medium. The asexual stage is induced by a high salt concentration, which completely inhibits sexual development. The taxon is therefore appropriate for investigating the mechanisms of asexual and sexual reproduction in fungi. In this study, de novo genome sequencing and analysis of transcriptomes during culture under high- and low-osmolarity conditions were performed. These analyses facilitated investigation of the evolution of mating-type genes, which determine the mode of sexual reproduction, in A. cristatus, the response of the high-osmolarity glycerol (HOG) pathway to osmotic stimulation, and the detection of mycotoxins and evaluation of the relationship with the location of the encoding genes.ResultsThe A. cristatus genome comprised 27.9 Mb and included 68 scaffolds, from which 10,136 protein-coding gene models were predicted. A phylogenetic analysis suggested a considerable phylogenetic distance between A. cristatus and A. nidulans. Comparison of the mating-type gene loci among Aspergillus species indicated that the mode in A. cristatus differs from those in other Aspergillus species. The components of the HOG pathway were conserved in the genome of A. cristatus. Differential gene expression analysis in A. cristatus using RNA-Seq demonstrated that the expression of most genes in the HOG pathway was unaffected by osmotic pressure. No gene clusters associated with the production of carcinogens were detected.ConclusionsA model of the mating-type locus in A. cristatus is reported for the first time. Aspergillus cristatus has evolved various mechanisms to cope with high osmotic stress. As a fungus associated with Fuzhuan tea, it is considered to be safe under low- and high-osmolarity conditions.

Comparative transcriptome analysis of the lichen-forming fungus Endocarpon pusillum elucidates its drought adaptation mechanisms

The lichen-forming fungus was isolated from the desert lichen Endocarpon pusillum that is extremely drought resistant. To understand the molecular mechanisms of drought resistance in the fungus, we employed RNA-seq and quantitative real-time PCR to compare and characterize the differentially expressed genes in pure culture at two different water levels and with that in desiccated lichen. The comparative transcriptome analysis indicated that a total of 1781 genes were differentially expressed between samples cultured under normal and PEG-induced drought stress conditions. Similar to those in drought resistance plants and non-lichenized fungi, the common drought-resistant mechanisms were differentially expressed in E. pusillum. However, the expression change of genes involved in osmotic regulation in E. pusillum is different, which might be the evidence for the feature of drought adaptation. Interestingly, different from other organisms, some genes involved in drought adaption mechanisms showed significantly different expression patterns between the presence and absence of drought stress in E. pusillum. The expression of 23 candidate stress responsive genes was further confirmed by quantitative real-time PCR using dehydrated E. pusillum lichen thalli. This study provides a valuable resource for future research on lichen-forming fungi and shall facilitate future functional studies of the specific genes related to drought resistance.

Functional Operons in Secondary Metabolic Gene Clusters in Glarea lozoyensis (Fungi, Ascomycota, Leotiomycetes)

ABSTRACT Operons are multigene transcriptional units which occur mostly in prokaryotes but rarely in eukaryotes. Protein-coding operons have not been reported in the Fungi even though they represent a very diverse kingdom of organisms. Here, we report a functional operon involved in the secondary metabolism of the fungus Glarea lozoyensis belonging to Leotiomycetes (Ascomycota). Two contiguous genes, glpks3 and glnrps7, encoding polyketide synthase and nonribosomal peptide synthetase, respectively, are cotranscribed into one dicistronic mRNA under the control of the same promoter, and the mRNA is then translated into two individual proteins, GLPKS3 and GLNRPS7. Heterologous expression in Aspergillus nidulans shows that the GLPKS3-GLNRPS7 enzyme complex catalyzes the biosynthesis of a novel pyrrolidinedione-containing compound, xenolozoyenone (compound 1), which indicates the operon is functional. Although it is structurally similar to prokaryotic operons, the glpks3-glnrps7 operon locus has a monophylogenic origin from fungi rather than having been horizontally transferred from prokaryotes. Moreover, two additional operons, glpks28-glnrps8 and glpks29-glnrps9, were verified at the transcriptional level in the same fungus. This is the first report of protein-coding operons in a member of the Fungi. IMPORTANCE Operons are multigene transcriptional units which occur mostly in prokaryotes but rarely in eukaryotes. Three operon-like gene structures for secondary metabolism that were discovered in the filamentous fungus Glarea lozoyensis are the first examples of protein-coding operons identified in a member of the Fungi. Among them, the glpks3-glnrps7 operon is responsible for the biosynthesis of xenolozoyenone, which is a novel tetramic acid-containing compound. Although structurally similar to prokaryotic operons, the glpks3-glnrps7 operon locus did not result from horizontal gene transfer from prokaryotes. In addition, operonlike structures have been predicted in silico to be common in other fungi. The common occurrence and operonlike structure in fungi provide evolutionary insight and essential data for eukaryotic gene transcription. Operons are multigene transcriptional units which occur mostly in prokaryotes but rarely in eukaryotes. Three operon-like gene structures for secondary metabolism that were discovered in the filamentous fungus Glarea lozoyensis are the first examples of protein-coding operons identified in a member of the Fungi. Among them, the glpks3-glnrps7 operon is responsible for the biosynthesis of xenolozoyenone, which is a novel tetramic acid-containing compound. Although structurally similar to prokaryotic operons, the glpks3-glnrps7 operon locus did not result from horizontal gene transfer from prokaryotes. In addition, operonlike structures have been predicted in silico to be common in other fungi. The common occurrence and operonlike structure in fungi provide evolutionary insight and essential data for eukaryotic gene transcription.