Comparative genome analysis across a kingdom of eukaryotic organisms: specialization and diversification in the fungi

Genome sequencing of Elaeocarpus spp. stem blight pathogen Pseudocryphonectria elaeocarpicola reveals potential adaptations to colonize woody bark

BACKGROUND

Elaeocarpus spp. stem blight, caused by Pseudocryphonectria elaeocarpicola, is a destructive disease, which will significantly reduce the productivity and longevity of Elaeocarpus spp. plants, especially in the Guangdong Province of China. However, few information is available for P. elaeocarpicola. To unravel the potential adaptation mechanism of stem adaptation, the whole genome of P. elaeocarpicola was sequenced by using the DNBSEQ and PacBio platforms.

RESULTS

P. elaeocarpicola harbors 44.49 Mb genome with 10,894 predicted coding genes. Genome analysis revealed that the P. elaeocarpicola genome encodes a plethora of pathogenicity-related genes. Analysis of carbohydrate-active enzymes (CAZymes) revealed a rich variety of enzymes participated in plant cell wall degradation, which could effectively degrade cellulose, hemicellulose and xyloglucans in the plant cell wall and promote the invasion of the host plant. There are 213 CAZyme families found in P. elaeocarpicola, among which glycoside hydrolase (GH) family has the largest number, far exceeding other tested fungi by 53%. Besides, P. elaeocarpicola has twice as many genes encoding chitin and cellulose degradation as Cryphonectria parasitica, which belong to the same family. The predicted typical secreted proteins of P. elaeocarpicola are numerous and functional, including many known virulence effector factors, indicating that P. elaeocarpicola has great potential to secrete virulence effectors to promote pathogenicity on host plants. AntiSMASH revealed that the genome encoded 61 secondary metabolic gene clusters including 86 secondary metabolic core genes which was much higher than C. parasitica (49). Among them, two gene cluster of P. elaeocarpicola, cluster12 and cluster52 showed 100% similarity with the mycotoxins synthesis clusters from Aspergillus steynii and Alternaria alternata, respectively. In addition, we annotated cytochrome P450 related enzymes, transporters, and transcription factors in P. elaeocarpicola, which are important virulence determinants of pathogenic fungi.

CONCLUSIONS

Taken together, our study represents the first genome assembly for P. elaeocarpicola and reveals the key virulence factors in the pathogenic process of P. elaeocarpicola, which will promote our understanding of its pathogenic mechanism. The acquired knowledge lays a foundation for further exploration of molecular interactions with the host and provide target for management strategies in future research.

Secondary metabolites and transcriptomic analysis of novel pulcherrimin producer Metschnikowia persimmonesis KIOM G15050: A potent and safe food biocontrol agent

Metschnikowia persimmonesis, a novel endophytic yeast strain isolated from Diospyros kaki calyx, possesses strong antimicrobial activity. We investigated its potential use as an environmentally safe food biocontrol agent through genomics, transcriptomics, and metabolomics. Secondary metabolites were isolated from M. persimmonesis, followed by chemical structure elucidation, PUL gene cluster identification, and RNA sequencing. Pulcherrimin was isolated using 2 M NaOH, its structure was confirmed, and the yield was quantified. Biocontrol efficacy of M. persimmonesis on persimmon fruits and calyx was evaluated by assessing lesion diameter and disease incidence. Following compounds were isolated from M. persimmonesis co-culture with Botrytis cinerea and Fusarium oxysporum: fusaric acid, benzoic acid, benzeneacetic acid, 4-hydroxybenzeneacetic acid, 4-(-2-hydoxyethyl)-benzoic acid, cyclo (Leu-Leu), benzenemethanol, 4-hydroxy-benzaldehide, 2-hydroxy-4-methoxy-benzoic acid, 4-hydroxy-benzoic acid, lumichrome, heptadecanoic acid, and nonadecanoic acid. Exposing M. persimmonesis to different growth media conditions (with or without sugar) resulted in the isolation of five compounds: Tyrosol, Cyclo (Pro-Val), cyclo(L-Pro-L-Tyr), cyclo(Leu-Leu), and cyclo(l-tyrosilylicine). Differentially expressed gene analysis revealed 3264 genes that were significantly expressed (fold change ≥2 and p-value ≤0.05) during M. persimmonesis growth in different media, of which only 270 (8.27%) showed altered expression in all sample combinations with Luria-Bertani Agar as control. Minimal media with ferric ions and tween-80 triggered the most gene expression changes, with the highest levels of PUL gene expression and pulcherrimin yield (262.166 mg/L) among all media treatments. M. persimmonesis also produced a higher amount of pulcherrimin (209.733 mg/L) than Metschnikowia pulcherrima (152.8 mg/L). M. persimmonesis inhibited the growth of Fusarium oxysporum in persimmon fruit and calyx. Toxicity evaluation of M. persimmonesis extracts showed no harmful effects on the liver and mitochondria of zebrafish, and no potential risk of cardiotoxicity in hERG-HEK293 cell lines. Thus, M. persimmonesis can be commercialized as a potent and safe biocontrol agent for preserving food products.

Genome-wide prediction and transcriptome analysis of sugar transporters in four ascomycete fungi

The import of plant-derived small sugars by sugar transporters (STs) has received increasing interest due to its important biological role and great industrial potential. STs are important targets of genetic engineering to improve fungal plant biomass conversion. Comparatively analysis of the genome-wide prevalence and transcriptomics of STs was performed in four filamentous fungi: Aspergillus niger, Aspergillus nidulans, Penicillium subrubescens and Trichoderma reesei. Using phylogenetic analysis and literature mining, their predicted STs were divided into ten subfamilies with putative sugar specificities assigned. In addition, transcriptome analysis revealed complex expression profiles among different STs subfamilies and fungal species, indicating a sophisticated transcriptome regulation and functional diversity of fungal STs. Several STs showed strong co-expression with other genes involved in sugar utilization, encoding CAZymes and sugar catabolic enzymes. This study provides new insights into the diversity of STs at the genomic/transcriptomic level, facilitating their biochemical characterization and metabolic engineering.

The NIH Comparative Genomics Resource: addressing the promises and challenges of comparative genomics on human health

Comparative genomics is the comparison of genetic information within and across organisms to understand the evolution, structure, and function of genes, proteins, and non-coding regions (Sivashankari and Shanmughavel, Bioinformation 1:376-8, 2007). Advances in sequencing technology and assembly algorithms have resulted in the ability to sequence large genomes and provided a wealth of data that are being used in comparative genomic analyses. Comparative analysis can be leveraged to systematically explore and evaluate the biological relationships and evolution between species, aid in understanding the structure and function of genes, and gain a better understanding of disease and potential drug targets. As our knowledge of genetics expands, comparative genomics can help identify emerging model organisms among a broader span of the tree of life, positively impacting human health. This impact includes, but is not limited to, zoonotic disease research, therapeutics development, microbiome research, xenotransplantation, oncology, and toxicology. Despite advancements in comparative genomics, new challenges have arisen around the quantity, quality assurance, annotation, and interoperability of genomic data and metadata. New tools and approaches are required to meet these challenges and fulfill the needs of researchers. This paper focuses on how the National Institutes of Health (NIH) Comparative Genomics Resource (CGR) can address both the opportunities for comparative genomics to further impact human health and confront an increasingly complex set of challenges facing researchers.

Genome-wide association mapping reveals genes underlying population-level metabolome diversity in a fungal crop pathogen

Background

Fungi produce a wide range of specialized metabolites (SMs) involved in biotic interactions. Pathways for the production of SMs are often encoded in clusters of tightly arranged genes identified as biosynthetic gene clusters. Such gene clusters can undergo horizontal gene transfers between species and rapid evolutionary change within species. The acquisition, rearrangement, and deletion of gene clusters can generate significant metabolome diversity. However, the genetic basis underlying variation in SM production remains poorly understood.

Results

Here, we analyzed the metabolite production of a large population of the fungal pathogen of wheat, Zymoseptoria tritici. The pathogen causes major yield losses and shows variation in gene clusters. We performed untargeted ultra-high performance liquid chromatography-high resolution mass spectrometry to profile the metabolite diversity among 102 isolates of the same species. We found substantial variation in the abundance of the detected metabolites among isolates. Integrating whole-genome sequencing data, we performed metabolite genome-wide association mapping to identify loci underlying variation in metabolite production (i.e., metabolite-GWAS). We found that significantly associated SNPs reside mostly in coding and gene regulatory regions. Associated genes encode mainly transport and catalytic activities. The metabolite-GWAS identified also a polymorphism in the 3′UTR region of a virulence gene related to metabolite production and showing expression variation.

Conclusions

Taken together, our study provides a significant resource to unravel polymorphism underlying metabolome diversity within a species. Integrating metabolome screens should be feasible for a range of different plant pathogens and help prioritize molecular studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01422-z.

Identifying and testing marker-trait associations for growth and phenology in three pine species: Implications for genomic prediction

In tree species, genomic prediction offers the potential to forecast mature trait values in early growth stages, if robust marker-trait associations can be identified. Here we apply a novel multispecies approach using genotypes from a new genotyping array, based on 20,795 single nucleotide polymorphisms (SNPs) from three closely related pine species (Pinus sylvestris, Pinus uncinata and Pinus mugo), to test for associations with growth and phenology data from a common garden study. Predictive models constructed using significantly associated SNPs were then tested and applied to an independent multisite field trial of P. sylvestris and the capability to predict trait values was evaluated. One hundred and eighteen SNPs showed significant associations with the traits in the pine species. Common SNPs (MAF > 0.05) associated with bud set were only found in genes putatively involved in growth and development, whereas those associated with growth and budburst were also located in genes putatively involved in response to environment and, to a lesser extent, reproduction. At one of the two independent sites, the model we developed produced highly significant correlations between predicted values and observed height data (YA, height 2020: r = 0.376, p < 0.001). Predicted values estimated with our budburst model were weakly but positively correlated with duration of budburst at one of the sites (GS, 2015: r = 0.204, p = 0.034; 2018: r = 0.205, p = 0.034-0.037) and negatively associated with budburst timing at the other (YA: r = -0.202, p = 0.046). Genomic prediction resulted in the selection of sets of trees whose mean height was taller than the average for each site. Our results provide tentative support for the capability of prediction models to forecast trait values in trees, while highlighting the need for caution in applying them to trees grown in different environments.

Evolutionary Insights Into Two Widespread Ectomycorrhizal Fungi (Pisolithus) From Comparative Analysis of Mitochondrial Genomes

The genus Pisolithus is a group of global ectomycorrhizal fungi. The characterizations of Pisolithus mitochondrial genomes have still been unknown. In the present study, the complete mitogenomes of two Pisolithus species, Pisolithus microcarpus, and Pisolithus tinctorius, were assembled and compared with other Boletales mitogenomes. Both Pisolithus mitogenomes comprised circular DNA molecules with sizes of 43,990 and 44,054 bp, respectively. Comparative mitogenomic analysis showed that the rps3 gene differentiated greatly between Boletales species, and this gene may be subjected to strong pressure of positive selection between some Boletales species. Several plasmid-derived genes and genes with unknown functions were detected in the two Pisolithus mitogenomes, which needs further analysis. The two Pisolithus species show a high degree of collinearity, which may represent the gene arrangement of the ancestors of ectomycorrhizal Boletales species. Frequent intron loss/gain events were detected in Boletales and basidiomycetes, and intron P717 was only detected in P. tinctorius out of the eight Boletales mitogenomes tested. We reconstructed phylogeny of 79 basidiomycetes based on combined mitochondrial gene dataset, and obtained well-supported phylogenetic topologies. This study served as the first report on the mitogenomes of the family Pisolithaceae, which will promote the understanding of the evolution of Pisolithus species.

The mitogenomes of two saprophytic Boletales species (Coniophora) reveals intron dynamics and accumulation of plasmid-derived and non-conserved genes

The order Boletales is a group of fungi with complex life styles, which include saprophytic and ectomycorrhizal mushroom-forming fungi. In the present study, the complete mitogenomes of two saprophytic Boletales species, Coniophora olivacea, and C. puteana, were assembled and compared with mitogenomes of ectomycorrhizal Boletales. Both mitogenomes comprised circular DNA molecules with sizes of 78,350 bp and 79,655 bp, respectively. Comparative mitogenomic analysis indicated that the two saprophytic Boletales species contained more plasmid-derived (7 on average) and unknown functional genes (12 on average) than the four ectomycorrhizal Boletales species previously reported. In addition, the core protein coding genes, nad2 and rps3, were found to be subjected to positive selection pressure between some Boletales species. Frequent intron gain/loss events were detected in Boletales and Basidiomycetes, and several novel intron classes were found in two Coniophora species. A total of 33 introns were detected in C. olivacea, and most were found to have undergone contraction in the C. olivacea mitogenome. Mitochondrial genes of Coniophora species were found to have undergone large-scale gene rearrangements, and the accumulation of intra-genomic repeats in the mitogenome was considered as one of the main contributing factors. Based on combined mitochondrial gene sets, we obtained a well-supported phylogenetic tree for 76 Basidiomycetes, demonstrating the utility of mitochondrial gene analysis for inferring Basidiomycetes phylogeny. The study served as the first report on the mitogenomes of the family Coniophorineae, which will help to understand the origin and evolution patterns of Boletales species with complex lifestyles.

In silico Proteomic Analysis Provides Insights Into Phylogenomics and Plant Biomass Deconstruction Potentials of the Tremelalles

Basidiomycetes populate a wide range of ecological niches but unlike ascomycetes, their capabilities to decay plant polymers and their potential for biotechnological approaches receive less attention. Particularly, identification and isolation of CAZymes is of biotechnological relevance and has the potential to improve the cache of currently available commercial enzyme cocktails toward enhanced plant biomass utilization. The order Tremellales comprises phylogenetically diverse fungi living as human pathogens, mycoparasites, saprophytes or associated with insects. Here, we have employed comparative genomics approaches to highlight the phylogenomic relationships among thirty-five Tremellales and to identify putative enzymes of biotechnological interest encoded on their genomes. Evaluation of the predicted proteomes of the thirty-five Tremellales revealed 6,918 putative carbohydrate-active enzymes (CAZYmes) and 7,066 peptidases. Two soil isolates, Saitozyma podzolica DSM 27192 and Cryptococcus sp. JCM 24511, show higher numbers harboring an average of 317 compared to a range of 267-121 CAZYmes for the rest of the strains. Similarly, the proteomes of the two soil isolates along with two plant associated strains contain higher number of peptidases sharing an average of 234 peptidases compared to a range of 226-167 for the rest of the strains. Despite these huge differences and the apparent enrichment of these enzymes among the soil isolates, the data revealed a diversity of the various enzyme families that does not reflect specific habitat type. Growth experiment on various carbohydrates to validate the predictions provides support for this view. Overall, the data indicates that the Tremellales could serve as a rich source of both CAZYmes and peptidases with wide range of potential biotechnological relevance.

Greater genetic and regulatory plasticity of retained duplicates in Epichloë endophytic fungi

Gene duplicates can act as a source of genetic material from which new functions arise. Most duplicated genes revert to single copy genes and only a small proportion are retained. However, it remains unclear why some duplicate genes persist in the genome for an extended time. We investigate this question by analysing retained gene duplicates in the fungal genus Epichloë, ascomycete fungi that form close endophytic symbioses with their host grasses. Retained duplicates within this genus have two independent origins, but both long pre-date the origin and diversification of the genus Epichloë. We find that loss of retained duplicates within the genus is frequent and often associated with speciation. Retained duplicates have faster evolutionary rates (Ka) and show relaxed selection (Ka/Ks) compared to single copy genes. Both features are time-dependent. Through comparison of conspecific strains, we find greater evolutionary rates in coding regions and sequence divergence in regulatory regions of retained duplicates than single copy genes, with this pattern more pronounced for strains adapted to different grass host species. Consistent with this sequence divergence in regulatory regions, transcriptome analyses show greater expression variation of retained duplicates than single copy genes. This suggest that cis-regulatory changes make important contributions to the expression patterns of retained duplicates. Coupled with supporting observations from the model yeast Saccharomyces cerevisiae, these data suggest that genetic robustness and regulatory plasticity are common drivers behind the retention of duplicated genes in fungi.

In Silico Analysis of Putative Sugar Transporter Genes in Aspergillus niger Using Phylogeny and Comparative Transcriptomics

Aspergillus niger is one of the most widely used fungi to study the conversion of the lignocellulosic feedstocks into fermentable sugars. Understanding the sugar uptake system of A. niger is essential to improve the efficiency of the process of fungal plant biomass degradation. In this study, we report a comprehensive characterization of the sugar transportome of A. niger by combining phylogenetic and comparative transcriptomic analyses. We identified 86 putative sugar transporter (ST) genes based on a conserved protein domain search. All these candidates were then classified into nine subfamilies and their functional motifs and possible sugar-specificity were annotated according to phylogenetic analysis and literature mining. Furthermore, we comparatively analyzed the ST gene expression on a large set of fungal growth conditions including mono-, di- and polysaccharides, and mutants of transcriptional regulators. This revealed that transporter genes from the same phylogenetic clade displayed very diverse expression patterns and were regulated by different transcriptional factors. The genome-wide study of STs of A. niger provides new insights into the mechanisms underlying an extremely flexible metabolism and high nutritional versatility of A. niger and will facilitate further biochemical characterization and industrial applications of these candidate STs.

Identification of candidate pathogenicity determinants of Rhizoctonia solani AG1-IA, which causes sheath blight disease in rice

Sheath blight disease is one of the predominant diseases of rice and it is caused by the necrotrophic fungal pathogen Rhizoctonia solani. The mechanistic insight about its widespread success as a broad host range pathogen is limited. In this study, we endeavor to identify pathogenicity determinants of R. solani during infection process in rice. Through RNAseq analysis, we identified a total of 65 and 232 R. solani (strain BRS1) genes to be commonly upregulated in three different rice genotypes (PB1, Tetep, and TP309) at establishment and necrotrophic phase, respectively. The induction of genes encoding extracellular protease, ABC transporter, and transcription factors were notable during establishment phase. While during necrotrophic phase, several CAZymes, sugar transporters, cellular metabolism, and protein degradation-related genes were prominently induced. We have also identified few putative secreted effector encoding genes that were upregulated during pathogenesis. The qPCR analysis further validated the phase-specific expression dynamics of some selected putative effectors and pathogenicity-associated genes. Overall, the present study reports identification of key genes and processes that might be crucial for R. solani pathogenesis. The ability to effectively damage host cell wall and survive in hostile plant environment by managing oxidative stress, cytotoxic compounds, etc. is being proposed to be important for pathogenesis of R. solani in rice. The functional characterization of these genes would provide key insights about this important pathosystem and facilitate development of strategies to control this devastating disease.

Evolutionary Analysis of Pectin Lyases of the Genus Colletotrichum

Pectin lyases (PNLs) are important enzymes that are involved in plant cell wall degradation during the infection process. Colletotrichum is a diverse genus of fungi, which allows the study of the evolution of PNLs and their possible role in pathogen-host interactions and lifestyle adaptations. The phylogenetic reconstruction of PNLs from Colletotrichum and analysis of selection pressures showed the formation of protein lineages by groups of species with different selection pressures and specific patterns. The analysis of positive selection at individual sites using different methods allowed for the identification of three codons with evidence of positive selection in the oligosaccharide-binding region and two codons on the antiparallel sheet, which may influence the interaction with the substrate. Seven codons on the surface of the protein, mainly in the peripheral helices of the PNLs, could have an important function in evasion of plant defenses, as has been proposed in other enzymes. According to our results, it is possible that events of genetic duplication occurred in ancestral lines, followed by episodes of genetic diversification and gene loss, probably influenced by differences in the composition of the host cell wall. Additionally, different patterns of evolution in Colletotrichum appear to be molded by a strong purifying selection and positive selection episodes that forged the observed evolutionary patterns, possibly influenced by host interaction or substrate specificity. This work represents a starting point for the study of sites that may be important for evasion of plant defenses and biotechnological purposes.

Immunity and starvation: new opportunities to elevate disease resistance in crops

Plants use multiple mechanisms to defend themselves against invading microbes. Besides using their immune system to surveil and eliminate pathogens, plants actively block the pathogens' access to nutrients as an alternative way to prevent colonization. In this review, we focus on immunity and starvation as major obstacles for pathogens' adaptation. We summarize the key mechanisms employed by pathogens to modulate host immunity and to guarantee sugar uptake. In contrast to genes that deal with the immune system and show high levels of plasticity, pathogen genes involved in sugar acquisition are highly conserved, and may not have adapted to co-evolving interactions with the host. We propose a model to assess the durability of different control strategies based on the ability of pathogens to deal with host immunity or starvation. This analysis opens new opportunities to elevate disease resistance in crops by reducing the likelihood of pathogen adaptation.

A plant pathology perspective of fungal genome sequencing

The majority of plant pathogens are fungi and many of these adversely affect food security. This mini-review aims to provide an analysis of the plant pathogenic fungi for which genome sequences are publically available, to assess their general genome characteristics, and to consider how genomics has impacted plant pathology. A list of sequenced fungal species was assembled, the taxonomy of all species verified, and the potential reason for sequencing each of the species considered. The genomes of 1090 fungal species are currently (October 2016) in the public domain and this number is rapidly rising. Pathogenic species comprised the largest category (35.5 %) and, amongst these, plant pathogens are predominant. Of the 191 plant pathogenic fungal species with available genomes, 61.3 % cause diseases on food crops, more than half of which are staple crops. The genomes of plant pathogens are slightly larger than those of other fungal species sequenced to date and they contain fewer coding sequences in relation to their genome size. Both of these factors can be attributed to the expansion of repeat elements. Sequenced genomes of plant pathogens provide blueprints from which potential virulence factors were identified and from which genes associated with different pathogenic strategies could be predicted. Genome sequences have also made it possible to evaluate adaptability of pathogen genomes and genomic regions that experience selection pressures. Some genomic patterns, however, remain poorly understood and plant pathogen genomes alone are not sufficient to unravel complex pathogen-host interactions. Genomes, therefore, cannot replace experimental studies that can be complex and tedious. Ultimately, the most promising application lies in using fungal plant pathogen genomics to inform disease management and risk assessment strategies. This will ultimately minimize the risks of future disease outbreaks and assist in preparation for emerging pathogen outbreaks.

Multiple isoforms for the catalytic subunit of PKA in the basal fungal lineage Mucor circinelloides

Protein kinase A (PKA) activity is involved in dimorphism of the basal fungal lineage Mucor. From the recently sequenced genome of Mucor circinelloides we could predict ten catalytic subunits of PKA. From sequence alignment and structural prediction we conclude that the catalytic core of the isoforms is conserved, and the difference between them resides in their amino termini. This high number of isoforms is maintained in the subdivision Mucoromycotina. Each paralogue, when compared to the ones form other fungi is more homologous to one of its orthologs than to its paralogs. All of these fungal isoforms cannot be included in the class I or II in which fungal protein kinases have been classified. mRNA levels for each isoform were measured during aerobic and anaerobic growth. The expression of each isoform is differential and associated to a particular growth stage. We reanalyzed the sequence of PKAC (GI 20218944), the only cloned sequence available until now for a catalytic subunit of M. circinelloides. PKAC cannot be classified as a PKA because of its difference in the conserved C-tail; it shares with PKB a conserved C2 domain in the N-terminus. No catalytic activity could be measured for this protein nor predicted bioinformatically. It can thus be classified as a pseudokinase. Its importance can not be underestimated since it is expressed at the mRNA level in different stages of growth, and its deletion is lethal.

Differential gene retention as an evolutionary mechanism to generate biodiversity and adaptation in yeasts

The evolutionary history of the characters underlying the adaptation of microorganisms to food and biotechnological uses is poorly understood. We undertook comparative genomics to investigate evolutionary relationships of the dairy yeast Geotrichum candidum within Saccharomycotina. Surprisingly, a remarkable proportion of genes showed discordant phylogenies, clustering with the filamentous fungus subphylum (Pezizomycotina), rather than the yeast subphylum (Saccharomycotina), of the Ascomycota. These genes appear not to be the result of Horizontal Gene Transfer (HGT), but to have been specifically retained by G. candidum after the filamentous fungi-yeasts split concomitant with the yeasts' genome contraction. We refer to these genes as SRAGs (Specifically Retained Ancestral Genes), having been lost by all or nearly all other yeasts, and thus contributing to the phenotypic specificity of lineages. SRAG functions include lipases consistent with a role in cheese making and novel endoglucanases associated with degradation of plant material. Similar gene retention was observed in three other distantly related yeasts representative of this ecologically diverse subphylum. The phenomenon thus appears to be widespread in the Saccharomycotina and argues that, alongside neo-functionalization following gene duplication and HGT, specific gene retention must be recognized as an important mechanism for generation of biodiversity and adaptation in yeasts.

Distinctive expansion of gene families associated with plant cell wall degradation, secondary metabolism, and nutrient uptake in the genomes of grapevine trunk pathogens

BACKGROUND

Trunk diseases threaten the longevity and productivity of grapevines in all viticulture production systems. They are caused by distantly-related fungi that form chronic wood infections. Variation in wood-decay abilities and production of phytotoxic compounds are thought to contribute to their unique disease symptoms. We recently released the draft sequences of Eutypa lata, Neofusicoccum parvum and Togninia minima, causal agents of Eutypa dieback, Botryosphaeria dieback and Esca, respectively. In this work, we first expanded genomic resources to three important trunk pathogens, Diaporthe ampelina, Diplodia seriata, and Phaeomoniella chlamydospora, causal agents of Phomopsis dieback, Botryosphaeria dieback, and Esca, respectively. Then we integrated all currently-available information into a genome-wide comparative study to identify gene families potentially associated with host colonization and disease development.

RESULTS

The integration of RNA-seq, comparative and ab initio approaches improved the protein-coding gene prediction in T. minima, whereas shotgun sequencing yielded nearly complete genome drafts of Dia. ampelina, Dip. seriata, and P. chlamydospora. The predicted proteomes of all sequenced trunk pathogens were annotated with a focus on functions likely associated with pathogenesis and virulence, namely (i) wood degradation, (ii) nutrient uptake, and (iii) toxin production. Specific patterns of gene family expansion were described using Computational Analysis of gene Family Evolution, which revealed lineage-specific evolution of distinct mechanisms of virulence, such as specific cell wall oxidative functions and secondary metabolic pathways in N. parvum, Dia. ampelina, and E. lata. Phylogenetically-informed principal component analysis revealed more similar repertoires of expanded functions among species that cause similar symptoms, which in some cases did not reflect phylogenetic relationships, thereby suggesting patterns of convergent evolution.

CONCLUSIONS

This study describes the repertoires of putative virulence functions in the genomes of ubiquitous grapevine trunk pathogens. Gene families with significantly faster rates of gene gain can now provide a basis for further studies of in planta gene expression, diversity by genome re-sequencing, and targeted reverse genetic approaches. The functional validation of potential virulence factors will lead to a more comprehensive understanding of the mechanisms of pathogenesis and virulence, which ultimately will enable the development of accurate diagnostic tools and effective disease management.

A Fox2-dependent fatty acid ß-oxidation pathway coexists both in peroxisomes and mitochondria of the ascomycete yeast Candida lusitaniae

It is generally admitted that the ascomycete yeasts of the subphylum Saccharomycotina possess a single fatty acid ß-oxidation pathway located exclusively in peroxisomes, and that they lost mitochondrial ß-oxidation early during evolution. In this work, we showed that mutants of the opportunistic pathogenic yeast Candida lusitaniae which lack the multifunctional enzyme Fox2p, a key enzyme of the ß-oxidation pathway, were still able to grow on fatty acids as the sole carbon source, suggesting that C. lusitaniae harbored an alternative pathway for fatty acid catabolism. By assaying 14Cα-palmitoyl-CoA consumption, we demonstrated that fatty acid catabolism takes place in both peroxisomal and mitochondrial subcellular fractions. We then observed that a fox2Δ null mutant was unable to catabolize fatty acids in the mitochondrial fraction, thus indicating that the mitochondrial pathway was Fox2p-dependent. This finding was confirmed by the immunodetection of Fox2p in protein extracts obtained from purified peroxisomal and mitochondrial fractions. Finally, immunoelectron microscopy provided evidence that Fox2p was localized in both peroxisomes and mitochondria. This work constitutes the first demonstration of the existence of a Fox2p-dependent mitochondrial β-oxidation pathway in an ascomycetous yeast, C. lusitaniae. It also points to the existence of an alternative fatty acid catabolism pathway, probably located in peroxisomes, and functioning in a Fox2p-independent manner.

Functional study of the Hap4-like genes suggests that the key regulators of carbon metabolism HAP4 and oxidative stress response YAP1 in yeast diverged from a common ancestor

The transcriptional regulator HAP4, induced by respiratory substrates, is involved in the balance between fermentation and respiration in S. cerevisiae. We identified putative orthologues of the Hap4 protein in all ascomycetes, based only on a conserved sixteen amino acid-long motif. In addition to this motif, some of these proteins contain a DNA-binding motif of the bZIP type, while being nonetheless globally highly divergent. The genome of the yeast Hansenula polymorpha contains two HAP4-like genes encoding the protein HpHap4-A which, like ScHap4, is devoid of a bZIP motif, and HpHap4-B which contains it. This species has been chosen for a detailed examination of their respective properties. Based mostly on global gene expression studies performed in the S. cerevisiae HAP4 disruption mutant (ScΔhap4), we show here that HpHap4-A is functionally equivalent to ScHap4, whereas HpHap4-B is not. Moreover HpHAP4-B is able to complement the H₂O₂ hypersensitivity of the ScYap1 deletant, YAP1 being, in S. cerevisiae, the main regulator of oxidative stress. Finally, a transcriptomic analysis performed in the ScΔyap1 strain overexpressing HpHAP4-B shows that HpHap4-B acts both on oxidative stress response and carbohydrate metabolism in a manner different from both ScYap1 and ScHap4. Deletion of these two genes in their natural host, H. polymorpha, confirms that HpHAP4-A participates in the control of the fermentation/respiration balance, while HpHAP4-B is involved in oxidative stress since its deletion leads to hypersensitivity to H₂O₂. These data, placed in an evolutionary context, raise new questions concerning the evolution of the HAP4 transcriptional regulation function and suggest that Yap1 and Hap4 have diverged from a unique regulatory protein in the fungal ancestor.

The evolution of fungal metabolic pathways

Fungi contain a remarkable range of metabolic pathways, sometimes encoded by gene clusters, enabling them to digest most organic matter and synthesize an array of potent small molecules. Although metabolism is fundamental to the fungal lifestyle, we still know little about how major evolutionary processes, such as gene duplication (GD) and horizontal gene transfer (HGT), have interacted with clustered and non-clustered fungal metabolic pathways to give rise to this metabolic versatility. We examined the synteny and evolutionary history of 247,202 fungal genes encoding enzymes that catalyze 875 distinct metabolic reactions from 130 pathways in 208 diverse genomes. We found that gene clustering varied greatly with respect to metabolic category and lineage; for example, clustered genes in Saccharomycotina yeasts were overrepresented in nucleotide metabolism, whereas clustered genes in Pezizomycotina were more common in lipid and amino acid metabolism. The effects of both GD and HGT were more pronounced in clustered genes than in their non-clustered counterparts and were differentially distributed across fungal lineages; specifically, GD, which was an order of magnitude more abundant than HGT, was most frequently observed in Agaricomycetes, whereas HGT was much more prevalent in Pezizomycotina. The effect of HGT in some Pezizomycotina was particularly strong; for example, we identified 111 HGT events associated with the 15 Aspergillus genomes, which sharply contrasts with the 60 HGT events detected for the 48 genomes from the entire Saccharomycotina subphylum. Finally, the impact of GD within a metabolic category was typically consistent across all fungal lineages, whereas the impact of HGT was variable. These results indicate that GD is the dominant process underlying fungal metabolic diversity, whereas HGT is episodic and acts in a category- or lineage-specific manner. Both processes have a greater impact on clustered genes, suggesting that metabolic gene clusters represent hotspots for the generation of fungal metabolic diversity.

On the origins of parasite-extended phenotypes

Parasites that adaptively manipulate the behavior of their host are among the most exciting adaptations that we can find in nature. The behavior of the host can become an extended phenotype of the parasites within animals such that the success and failure of the parasite's genome rely on precise change of the host's behavior. Evolutionary biology was borne from the close attention of naturalists such as Wallace and Darwin to phenotypic variation in seeking to understand the origins of new species. In this essay, I argue that we also need to think about the origins of parasite-extended phenotypes. This is a more difficult task than understanding the evolution of textbook examples of novelty such as the eyes of vertebrates or the hooves of horses. However, new tools such as phylogenomics provide an important opportunity to make significant progress in understanding the extended phenotypes of parasites. Knowing the origins of parasite-extended phenotypes is important as a goal all by itself. But the knowledge gained will also help us understand why complex manipulation is so rare and to identify the evolutionary tipping points driving its appearance.

Genomic evolution of 11 type strains within family Planctomycetaceae

The species in family Planctomycetaceae are ideal groups for investigating the origin of eukaryotes. Their cells are divided by a lipidic intracytoplasmic membrane and they share a number of eukaryote-like molecular characteristics. However, their genomic structures, potential abilities, and evolutionary status are still unknown. In this study, we searched for common protein families and a core genome/pan genome based on 11 sequenced species in family Planctomycetaceae. Then, we constructed phylogenetic tree based on their 832 common protein families. We also annotated the 11 genomes using the Clusters of Orthologous Groups database. Moreover, we predicted and reconstructed their core/pan metabolic pathways using the KEGG (Kyoto Encyclopedia of Genes and Genomes) orthology system. Subsequently, we identified genomic islands (GIs) and structural variations (SVs) among the five complete genomes and we specifically investigated the integration of two Planctomycetaceae plasmids in all 11 genomes. The results indicate that Planctomycetaceae species share diverse genomic variations and unique genomic characteristics, as well as have huge potential for human applications.

Genome sequencing of high-penicillin producing industrial strain of Penicillium chrysogenum

BACKGROUND

Due to the importance of Penicillium chrysogenum holding in medicine, the genome of low-penicillin producing laboratorial strain Wisconsin54-1255 had been sequenced and fully annotated. Through classical mutagenesis of Wisconsin54-1255, product titers and productivities of penicillin have dramatically increased, but what underlying genome structural variations is still little known. Therefore, genome sequencing of a high-penicillin producing industrial strain is very meaningful.

RESULTS

To reveal more insights into the genome structural variations of high-penicillin producing strain, we sequenced an industrial strain P. chrysogenum NCPC10086. By whole genome comparative analysis, we observed a large number of mutations, insertions and deletions, and structural variations. There are 69 new genes that not exist in the genome sequence of Wisconsin54-1255 and some of them are involved in energy metabolism, nitrogen metabolism and glutathione metabolism. Most importantly, we discovered a 53.7 Kb "new shift fragment" in a seven copies of determinative penicillin biosynthesis cluster in NCPC10086 and the arrangement type of amplified region is unique. Moreover, we presented two large-scale translocations in NCPC10086, containing genes involved energy, nitrogen metabolism and peroxysome pathway. At last, we found some non-synonymous mutations in the genes participating in homogentisate pathway or working as regulators of penicillin biosynthesis.

CONCLUSIONS

We provided the first high-quality genome sequence of industrial high-penicillin strain of P. chrysogenum and carried out a comparative genome analysis with a low-producing experimental strain. The genomic variations we discovered are related with energy metabolism, nitrogen metabolism and so on. These findings demonstrate the potential information for insights into the high-penicillin yielding mechanism and metabolic engineering in the future.

Fungal genes in context: genome architecture reflects regulatory complexity and function

Gene context determines gene expression, with local chromosomal environment most influential. Comparative genomic analysis is often limited in scope to conserved or divergent gene and protein families, and fungi are well suited to this approach with low functional redundancy and relatively streamlined genomes. We show here that one aspect of gene context, the amount of potential upstream regulatory sequence maintained through evolution, is highly predictive of both molecular function and biological process in diverse fungi. Orthologs with large upstream intergenic regions (UIRs) are strongly enriched in information processing functions, such as signal transduction and sequence-specific DNA binding, and, in the genus Aspergillus, include the majority of experimentally studied, high-level developmental and metabolic transcriptional regulators. Many uncharacterized genes are also present in this class and, by implication, may be of similar importance. Large intergenic regions also share two novel sequence characteristics, currently of unknown significance: they are enriched for plus-strand polypyrimidine tracts and an information-rich, putative regulatory motif that was present in the last common ancestor of the Pezizomycotina. Systematic consideration of gene UIR in comparative genomics, particularly for poorly characterized species, could help reveal organisms’ regulatory priorities.

Molecular evolution and patterns of duplication in the SEP/AGL6-like lineage of the Zingiberales: a proposed mechanism for floral diversification

The diversity of floral forms in the plant order Zingiberales has evolved through alterations in floral organ morphology. One striking alteration is the shift from fertile, filamentous stamens to sterile, laminar (petaloid) organs in the stamen whorls, attributed to specific pollination syndromes. Here, we examine the role of the SEPALLATA (SEP) genes, known to be important in regulatory networks underlying floral development and organ identity, in the evolution of development of the diverse floral organs phenotypes in the Zingiberales. Phylogenetic analyses show that the SEP-like genes have undergone several duplication events giving rise to multiple copies. Selection tests on the SEP-like genes indicate that the two copies of SEP3 have mostly evolved under balancing selection, probably due to strong functional restrictions as a result of their critical role in floral organ specification. In contrast, the two LOFSEP copies have undergone differential positive selection, indicating neofunctionalization. Reverse transcriptase-polymerase chain reaction, gene expression from RNA-seq data, and in situ hybridization analyses show that the recovered genes have differential expression patterns across the various whorls and organ types found in the Zingiberales. Our data also suggest that AGL6, sister to the SEP-like genes, may play an important role in stamen morphology in the Zingiberales. Thus, the SEP-like genes are likely to be involved in some of the unique morphogenetic patterns of floral organ development found among this diverse order of tropical monocots. This work contributes to a growing body of knowledge focused on understanding the role of gene duplications and the evolution of entire gene networks in the evolution of flower development.

Physiological role of Acyl coenzyme A synthetase homologs in lipid metabolism in Neurospora crassa

Acyl coenzyme A (CoA) synthetase (ACS) enzymes catalyze the activation of free fatty acids (FAs) to CoA esters by a two-step thioesterification reaction. Activated FAs participate in a variety of anabolic and catabolic lipid metabolic pathways, including de novo complex lipid biosynthesis, FA β-oxidation, and lipid membrane remodeling. Analysis of the genome sequence of the filamentous fungus Neurospora crassa identified seven putative fatty ACSs (ACS-1 through ACS-7). ACS-3 was found to be the major activator for exogenous FAs for anabolic lipid metabolic pathways, and consistent with this finding, ACS-3 localized to the endoplasmic reticulum, plasma membrane, and septa. Double-mutant analyses confirmed partial functional redundancy of ACS-2 and ACS-3. ACS-5 was determined to function in siderophore biosynthesis, indicating alternative functions for ACS enzymes in addition to fatty acid metabolism. The N. crassa ACSs involved in activation of FAs for catabolism were not specifically defined, presumably due to functional redundancy of several of ACSs for catabolism of exogenous FAs.

Genome analyses of the wheat yellow (stripe) rust pathogen Puccinia striiformis f. sp. tritici reveal polymorphic and haustorial expressed secreted proteins as candidate effectors

BACKGROUND

Wheat yellow (stripe) rust caused by Puccinia striiformis f. sp. tritici (PST) is one of the most devastating diseases of wheat worldwide. To design effective breeding strategies that maximize the potential for durable disease resistance it is important to understand the molecular basis of PST pathogenicity. In particular, the characterisation of the structure, function and evolutionary dynamics of secreted effector proteins that are detected by host immune receptors can help guide and prioritize breeding efforts. However, to date, our knowledge of the effector repertoire of cereal rust pathogens is limited.

RESULTS

We re-sequenced genomes of four PST isolates from the US and UK to identify effector candidates and relate them to their distinct virulence profiles. First, we assessed SNP frequencies between all isolates, with heterokaryotic SNPs being over tenfold more frequent (5.29 ± 2.23 SNPs/kb) than homokaryotic SNPs (0.41 ± 0.28 SNPs/kb). Next, we implemented a bioinformatics pipeline to integrate genomics, transcriptomics, and effector-focused annotations to identify and classify effector candidates in PST. RNAseq analysis highlighted transcripts encoding secreted proteins that were significantly enriched in haustoria compared to infected tissue. The expression of 22 candidate effector genes was characterised using qRT-PCR, revealing distinct temporal expression patterns during infection in wheat. Lastly, we identified proteins that displayed non-synonymous substitutions specifically between the two UK isolates PST-87/7 and PST-08/21, which differ in virulence to two wheat varieties. By focusing on polymorphic variants enriched in haustoria, we identified five polymorphic effector candidates between PST-87/7 and PST-08/21 among 2,999 secreted proteins. These allelic variants are now a priority for functional validation as virulence/avirulence effectors in the corresponding wheat varieties.

CONCLUSIONS

Integration of genomics, transcriptomics, and effector-directed annotation of PST isolates has enabled us to move beyond the single isolate-directed catalogues of effector proteins and develop a framework for mining effector proteins in closely related isolates and relate these back to their defined virulence profiles. This should ultimately lead to more comprehensive understanding of the PST pathogenesis system, an important first step towards developing more effective surveillance and management strategies for one of the most devastating pathogens of wheat.

Development of highly efficient, low-cost lignocellulolytic enzyme systems in the post-genomic era

The current high cost of lignocellulolytic enzymes is a major bottleneck in the economic bioconversion of lignocellulosic biomass to fuels and chemicals. Fungal lignocellulolytic enzyme systems are secreted at high levels, making them the most promising starting points for further development of highly efficient lignocellulolytic enzyme systems. In this paper, recent advances in improvement of fungal lignocellulolytic enzyme systems are reviewed, with an emphasis on the achievements made using genomic approaches. A general strategy for lignocellulolytic enzyme system development is proposed, including the improvement of the hydrolysis efficiencies and productivities of current enzyme systems. The applications of genomic, transcriptomic and proteomic analysis methods in examining the composition of native enzyme systems, discovery of novel enzymes and synergistic proteins from natural sources, and understanding of regulatory mechanisms for lignocellulolytic enzyme biosynthesis are summarized. By combining systems biology and synthetic biology tools, engineered fungal strains are expected to produce high levels of optimized lignocellulolytic enzyme systems.

Knowledge representation in metabolic pathway databases

The accurate representation of all aspects of a metabolic network in a structured format, such that it can be used for a wide variety of computational analyses, is a challenge faced by a growing number of researchers. Analysis of five major metabolic pathway databases reveals that each database has made widely different choices to address this challenge, including how to deal with knowledge that is uncertain or missing. In concise overviews, we show how concepts such as compartments, enzymatic complexes and the direction of reactions are represented in each database. Importantly, also concepts which a database does not represent are described. Which aspects of the metabolic network need to be available in a structured format and to what detail differs per application. For example, for in silico phenotype prediction, a detailed representation of gene-protein-reaction relations and the compartmentalization of the network is essential. Our analysis also shows that current databases are still limited in capturing all details of the biology of the metabolic network, further illustrated with a detailed analysis of three metabolic processes. Finally, we conclude that the conceptual differences between the databases, which make knowledge exchange and integration a challenge, have not been resolved, so far, by the exchange formats in which knowledge representation is standardized.

CFGP 2.0: a versatile web-based platform for supporting comparative and evolutionary genomics of fungi and Oomycetes

In 2007, Comparative Fungal Genomics Platform (CFGP; http://cfgp.snu.ac.kr/) was publicly open with 65 genomes corresponding to 58 fungal and Oomycete species. The CFGP provided six bioinformatics tools, including a novel tool entitled BLASTMatrix that enables search homologous genes to queries in multiple species simultaneously. CFGP also introduced Favorite, a personalized virtual space for data storage and analysis with these six tools. Since 2007, CFGP has grown to archive 283 genomes corresponding to 152 fungal and Oomycete species as well as 201 genomes that correspond to seven bacteria, 39 plants and 105 animals. In addition, the number of tools in Favorite increased to 27. The Taxonomy Browser of CFGP 2.0 allows users to interactively navigate through a large number of genomes according to their taxonomic positions. The user interface of BLASTMatrix was also improved to facilitate subsequent analyses of retrieved data. A newly developed genome browser, Seoul National University Genome Browser (SNUGB), was integrated into CFGP 2.0 to support graphical presentation of diverse genomic contexts. Based on the standardized genome warehouse of CFGP 2.0, several systematic platforms designed to support studies on selected gene families have been developed. Most of them are connected through Favorite to allow of sharing data across the platforms.

Identification and functional analysis of the erh1(+) gene encoding enhancer of rudimentary homolog from the fission yeast Schizosaccharomyces pombe

The ERH gene encodes a highly conserved small nuclear protein with a unique amino acid sequence and three-dimensional structure but unknown function. The gene is present in animals, plants, and protists but to date has only been found in few fungi. Here we report that ERH homologs are also present in all four species from the genus Schizosaccharomyces, S. pombe, S. octosporus, S. cryophilus, and S. japonicus, which, however, are an exception in this respect among Ascomycota and Basidiomycota. The ERH protein sequence is moderately conserved within the genus (58% identity between S. pombe and S.japonicus), but the intron-rich genes have almost identical intron-exon organizations in all four species. In S. pombe, erh1(+) is expressed at a roughly constant level during vegetative growth and adaptation to unfavorable conditions such as nutrient limitation and hyperosmotic stress caused by sorbitol. Erh1p localizes preferentially to the nucleus with the exception of the nucleolus, but is also present in the cytoplasm. Cells lacking erh1(+) have an aberrant cell morphology and a comma-like shape when cultured to the stationary phase, and exhibit a delayed recovery from this phase followed by slower growth. Loss of erh1(+) in an auxotrophic background results in enhanced arrest in the G1 phase following nutritional stress, and also leads to hypersensitivity to agents inducing hyperosmotic stress (sorbitol), inhibiting DNA replication (hydroxyurea), and destabilizing the plasma membrane (SDS); this hypersensitivity can be abolished by expression of S. pombe erh1(+) and, to a lesser extent, S. japonicus erh1(+) or human ERH. Erh1p fails to interact with the human Ciz1 and PDIP46/SKAR proteins, known molecular partners of human ERH. Our data suggest that in Schizosaccharomyces sp. erh1(+) is non-essential for normal growth and Erh1p could play a role in response to adverse environmental conditions and in cell cycle regulation.

Voriconazole hepatotoxicity in severe liver dysfunction

There are no studies regarding to these effects in patients with severe liver dysfunction.

OBJECTIVES

The aims of this study were to characterize voriconazole hepatotoxicity in patients with severe liver dysfunction and to compare it with a matched cohort treated with liposomal amphotericin B.

METHODS

This is an observational study, in which adults patients treated with at least 4 doses of voriconazole were included. Patients treated with liposomal amphotericin B were used as control group.

RESULTS

Sixty nine percent of patients treated with voriconazole showed changes in liver function tests (LFTs) during therapy. They showed elevated transaminases in 35%, cholestasis in 15% or a combination of both in 45%. According to the CTC classification, all patients with hepatotoxicity had a severe reaction. The Roussel Uclaf Causality Assessment Method score in all patients with hepatotoxicity was greater than 8. There was a correlation between initial loading dose greater than 300 mg (4.5 mg/kg) and the risk of hepatotoxicity (p < 0.001). The control group developed alterations in the LFTs in only 10.3% of patients.

CONCLUSION

Voriconazole should be used with caution in patients with severe liver dysfunction and following liver transplantation, with frequent monitoring of LFTs or using liposomal amphotericin B instead.

The fate of gene duplicates in the genomes of fungal pathogens

Understanding how molecular changes underlie phenotypic variation within and between species is one of the main goals of evolutionary biology and comparative genetics. The recent proliferation of sequenced fungal genomes offers a unique opportunity to start elucidating the extreme phenotypic diversity in the Kingdom Fungi.1-4 We attempted to investigate the contribution of gene families to the evolutionary forces shaping the diversity of pathogenic lifestyles among the fungi.5 We studied a family of secreted enzymes which is present and expanded in all genomes of fungal pathogens sequenced to date and absent from the genomes of true yeasts.3,4 This family of cutinases6 predates the division between the two major fungal phyla, Ascomycota and Basidiomycota.5 We discuss our molecular phylogenetic analyses, the number and sequence diversity, and gene gains and losses of cutinase family members between five Ascomycetes: the phytopathogens Magnaporthe oryzae, Fusarium graminearum and Botrytis cinerea; and the model organisms Neurospora crassa and Aspergillus nidulans.5 The functional characterization of three members of the M. oryzae cutinase family,6-10 coupled with the regulatory subfunctionalization and neofunctionalization of most gene pairs5 provide the first justification for the retention of paralogs after duplication and for gene redundancy in the genomes of fungal pathogens.

Defects in mitochondrial and peroxisomal β-oxidation influence virulence in the maize pathogen Ustilago maydis

An understanding of metabolic adaptation during the colonization of plants by phytopathogenic fungi is critical for developing strategies to protect crops. Lipids are abundant in plant tissues, and fungal phytopathogens in the phylum basidiomycota possess both peroxisomal and mitochondrial β-oxidation pathways to utilize this potential carbon source. Previously, we demonstrated a role for the peroxisomal β-oxidation enzyme Mfe2 in the filamentous growth, virulence, and sporulation of the maize pathogen Ustilago maydis. However, mfe2 mutants still caused disease symptoms, thus prompting a more detailed investigation of β-oxidation. We now demonstrate that a defect in the had1 gene encoding hydroxyacyl coenzyme A dehydrogenase for mitochondrial β-oxidation also influences virulence, although its paralog, had2, makes only a minor contribution. Additionally, we identified a gene encoding a polypeptide with similarity to the C terminus of Mfe2 and designated it Mfe2b; this gene makes a contribution to virulence only in the background of an mfe2Δ mutant. We also show that short-chain fatty acids induce cell death in U. maydis and that a block in β-oxidation leads to toxicity, likely because of the accumulation of toxic intermediates. Overall, this study reveals that β-oxidation has a complex influence on the formation of disease symptoms by U. maydis that includes potential metabolic contributions to proliferation in planta and an effect on virulence-related morphogenesis.

Peroxisomal and mitochondrial β-oxidation pathways influence the virulence of the pathogenic fungus Cryptococcus neoformans

An understanding of the connections between metabolism and elaboration of virulence factors during host colonization by the human-pathogenic fungus Cryptococcus neoformans is important for developing antifungal therapies. Lipids are abundant in host tissues, and fungal pathogens in the phylum basidiomycota possess both peroxisomal and mitochondrial β-oxidation pathways to utilize this potential carbon source. In addition, lipids are important signaling molecules in both fungi and mammals. In this report, we demonstrate that defects in the peroxisomal and mitochondrial β-oxidation pathways influence the growth of C. neoformans on fatty acids as well as the virulence of the fungus in a mouse inhalation model of cryptococcosis. Disease attenuation may be due to the cumulative influence of altered carbon source acquisition or processing, interference with secretion, changes in cell wall integrity, and an observed defect in capsule production for the peroxisomal mutant. Altered capsule elaboration in the context of a β-oxidation defect was unexpected but is particularly important because this trait is a major virulence factor for C. neoformans. Additionally, analysis of mutants in the peroxisomal pathway revealed a growth-promoting activity for C. neoformans, and subsequent work identified oleic acid and biotin as candidates for such factors. Overall, this study reveals that β-oxidation influences virulence in C. neoformans by multiple mechanisms that likely include contributions to carbon source acquisition and virulence factor elaboration.

Selection of a DNA barcode for Nectriaceae from fungal whole-genomes

A DNA barcode is a short segment of sequence that is able to distinguish species. A barcode must ideally contain enough variation to distinguish every individual species and be easily obtained. Fungi of Nectriaceae are economically important and show high species diversity. To establish a standard DNA barcode for this group of fungi, the genomes of Neurospora crassa and 30 other filamentous fungi were compared. The expect value was treated as a criterion to recognize homologous sequences. Four candidate markers, Hsp90, AAC, CDC48, and EF3, were tested for their feasibility as barcodes in the identification of 34 well-established species belonging to 13 genera of Nectriaceae. Two hundred and fifteen sequences were analyzed. Intra- and inter-specific variations and the success rate of PCR amplification and sequencing were considered as important criteria for estimation of the candidate markers. Ultimately, the partial EF3 gene met the requirements for a good DNA barcode: No overlap was found between the intra- and inter-specific pairwise distances. The smallest inter-specific distance of EF3 gene was 3.19%, while the largest intra-specific distance was 1.79%. In addition, there was a high success rate in PCR and sequencing for this gene (96.3%). CDC48 showed sufficiently high sequence variation among species, but the PCR and sequencing success rate was 84% using a single pair of primers. Although the Hsp90 and AAC genes had higher PCR and sequencing success rates (96.3% and 97.5%, respectively), overlapping occurred between the intra- and inter-specific variations, which could lead to misidentification. Therefore, we propose the EF3 gene as a possible DNA barcode for the nectriaceous fungi.

Conserved regulators of the cell separation process in Schizosaccharomyces

The fission yeasts (Schizosaccharomyces) representing a highly divergent phylogenetic branch of Fungi evolved from filamentous ancestors by gradual transition from mycelial growth to yeast morphology. For the transition, a mechanism had been developed that separates the sister cells after the completion of cytokinesis. Numerous components of the separation mechanism have been characterised in Schizosaccharomycespombe, including the zinc-finger transcription factor Ace2p and the fork-head transcription factor Sep1p. Here we show that both regulators have regions conserved within the genus. The most conserved parts contain the DNA-binding domains whose amino-acid sequences perfectly reflect the phylogenetic positions of the species. The less conserved parts of the proteins contain sequence blocks specific for the whole genus or only for the species propagating predominantly or exclusively as yeasts. Inactivation of either gene in the dimorphic species Schizosaccharomycesjaponicus abolished cell separation in the yeast phase conferring hypha-like morphology but did not change the growth pattern to unipolar and did not cause extensive polar vacuolation characteristic of the true mycelium. Neither mutation affected the mycelial phase, but both mutations hampered the hyphal fragmentation at the mycelium-to-yeast transition. Ace2p(Sj) acts downstream of Sep1p(Sj) and regulates the orthologues of the Ace2p-dependent S.pombe genes agn1(+) (1,3-alpha-glucanase) and eng1(+) (1,3-beta-glucanase) but does not regulate the orthologue of cfh4(+) (chitin synthase regulatory factor). These results and the complementation of the cell separation defects of the ace2(-) and sep1(-) mutations of S.pombe by heterologously expressed ace2(Sj) and sep1(Sj) indicate that the cell separation mechanism is conserved in the Schizosaccharomyces genus.

Causes and consequences of genome expansion in fungi

Fungi display a large diversity in genome size and complexity, variation that is often considered to be adaptive. But because nonadaptive processes can also have important consequences on the features of genomes, we investigated the relationship of genetic drift and genome size in the phylum Ascomycota using multiple indicators of genetic drift. We detected a complex relationship between genetic drift and genome size in fungi: genetic drift is associated with genome expansion on broad evolutionary timescales, as hypothesized for other eukaryotes; but within subphyla over smaller timescales, the opposite trend is observed. Moreover, fungi and bacteria display similar patterns of genome degradation that are associated with initial effects of genetic drift. We conclude that changes in genome size within Ascomycota have occurred using two different routes: large-scale genome expansions are catalyzed by increasing drift as predicted by the mutation-hazard model of genome evolution and small-scale modifications in genome size are independent of drift.

A consistent phylogenetic backbone for the fungi

The kingdom of fungi provides model organisms for biotechnology, cell biology, genetics, and life sciences in general. Only when their phylogenetic relationships are stably resolved, can individual results from fungal research be integrated into a holistic picture of biology. However, and despite recent progress, many deep relationships within the fungi remain unclear. Here, we present the first phylogenomic study of an entire eukaryotic kingdom that uses a consistency criterion to strengthen phylogenetic conclusions. We reason that branches (splits) recovered with independent data and different tree reconstruction methods are likely to reflect true evolutionary relationships. Two complementary phylogenomic data sets based on 99 fungal genomes and 109 fungal expressed sequence tag (EST) sets analyzed with four different tree reconstruction methods shed light from different angles on the fungal tree of life. Eleven additional data sets address specifically the phylogenetic position of Blastocladiomycota, Ustilaginomycotina, and Dothideomycetes, respectively. The combined evidence from the resulting trees supports the deep-level stability of the fungal groups toward a comprehensive natural system of the fungi. In addition, our analysis reveals methodologically interesting aspects. Enrichment for EST encoded data—a common practice in phylogenomic analyses—introduces a strong bias toward slowly evolving and functionally correlated genes. Consequently, the generalization of phylogenomic data sets as collections of randomly selected genes cannot be taken for granted. A thorough characterization of the data to assess possible influences on the tree reconstruction should therefore become a standard in phylogenomic analyses.

Berberine and itraconazole are not synergistic in vitro against Aspergillus fumigatus isolated from clinical patients

The incidence of Aspergillus fumigatus infections has become more frequent as a consequence of widespread immunosuppression. At present, the number of available antifungal agents in the clinic is limited, and most of them, such as itraconazole (ICZ), are toxic and show resistance. Berberine (BER) is a plant alkaloid used in the clinic mainly for alimentary infections. We have used BER and ICZ to measure in vitro resistance in A. fumigatus isolated from clinical patients. The minimum inhibitory concentration ranges of BER and ICZ were 4–256 and 0.031–0.250 μg/mL, respectively. In addition, against A. fumigatus IFM 40808 strain, the MIC₅₀ values of BER and ICZ were 8 and 0.125 μg/mL. Using this strain, we compared the giant colonies with or without BER, and concluded that BER could restrain A. fumigatus mycelial growth and conidial pigment production. Combinations of the two drugs were also tested by the checkerboard assay to identify any functional interactions between them. Thirty-two out of 42 isolates had FICI values > 4.0, indicating that two drugs were mutually antagonistic. In conclusion, it is not advised that BER and ICZ be used in the clinic at the same time. Our results indicated that BER may inhibit A. fumigatus through the ergosterol biosynthesis pathway, like ICZ.

Fusion of the subunits α and β of succinyl-CoA synthetase as a phylogenetic marker for Pezizomycotina fungi

Gene fusions, yielding the formation of multidomain proteins, are evolutionary events that can be utilized as phylogenetic markers. Here we describe a fusion gene comprising the α and β subunits of succinyl-coA synthetase, an enzyme of the TCA cycle, in Pezizomycotina fungi. This fusion is present in all Pezizomycotina with complete genome sequences and absent from all other organisms. Phylogenetic analysis of the α and β subunits of succinyl-CoA synthetase suggests that both subunits were duplicated and retained in Pezizomycotina while one copy was lost from other fungi. One of the duplicated copies was then fused in Pezizomycotina. Our results suggest that the fusion of the α and β subunits of succinyl-CoA synthetase can be used as a molecular marker for membership in the Pezizomycotina subphylum. If a species has the fusion it can be reliably classified as Pezizomycotina, while the absence of the fusion is suggestive that the species is not a member of this subphylum.

Addressing inter-gene heterogeneity in maximum likelihood phylogenomic analysis: yeasts revisited

Phylogenomic approaches to the resolution of inter-species relationships have become well established in recent years. Often these involve concatenation of many orthologous genes found in the respective genomes followed by analysis using standard phylogenetic models. Genome-scale data promise increased resolution by minimising sampling error, yet are associated with well-known but often inappropriately addressed caveats arising through data heterogeneity and model violation. These can lead to the reconstruction of highly-supported but incorrect topologies. With the aim of obtaining a species tree for 18 species within the ascomycetous yeasts, we have investigated the use of appropriate evolutionary models to address inter-gene heterogeneities and the scalability and validity of supermatrix analysis as the phylogenetic problem becomes more difficult and the number of genes analysed approaches truly phylogenomic dimensions. We have extended a widely-known early phylogenomic study of yeasts by adding additional species to increase diversity and augmenting the number of genes under analysis. We have investigated sophisticated maximum likelihood analyses, considering not only a concatenated version of the data but also partitioned models where each gene constitutes a partition and parameters are free to vary between the different partitions (thereby accounting for variation in the evolutionary processes at different loci). We find considerable increases in likelihood using these complex models, arguing for the need for appropriate models when analyzing phylogenomic data. Using these methods, we were able to reconstruct a well-supported tree for 18 ascomycetous yeasts spanning about 250 million years of evolution.

Elucidation of functional markers from Aspergillus nidulans developmental regulator FlbB and their phylogenetic distribution

Aspergillus nidulans is a filamentous fungus widely used as a model for biotechnological and clinical research. It is also used as a platform for the study of basic eukaryotic developmental processes. Previous studies identified and partially characterized a set of proteins controlling cellular transformations in this ascomycete. Among these proteins, the bZip type transcription factor FlbB is a key regulator of reproduction, stress responses and cell-death. Our aim here was the prediction, through various bioinformatic methods, of key functional residues and motifs within FlbB in order to inform the design of future laboratory experiments and further the understanding of the molecular mechanisms that control fungal development. A dataset of FlbB orthologs and those of its key interaction partner FlbE was assembled from 40 members of the Pezizomycotina. Unique features were identified in each of the three structural domains of FlbB. The N-terminal region encoded a bZip transcription factor domain with a novel histidine-containing DNA binding motif while the dimerization determinants exhibited two distinct profiles that segregated by class. The C-terminal region of FlbB showed high similarity with the AP-1 family of stress response regulators but with variable patterns of conserved cysteines that segregated by class and order. Motif conservation analysis revealed that nine FlbB orthologs belonging to the Eurotiales order contained a motif in the central region that could mediate interaction with FlbE. The key residues and motifs identified here provide a basis for the design of follow-up experimental investigations. Additionally, the presence or absence of these residues and motifs among the FlbB orthologs could help explain the differences in the developmental programs among fungal species as well as define putative complementation groups that could serve to extend known functional characterizations to other species.

Haploinsufficiency and the sex chromosomes from yeasts to humans

BACKGROUND

Haploinsufficient (HI) genes are those for which a reduction in copy number in a diploid from two to one results in significantly reduced fitness. Haploinsufficiency is increasingly implicated in human disease, and so predicting this phenotype could provide insights into the genetic mechanisms behind many human diseases, including some cancers.

RESULTS

In the present work we show that orthologues of Saccharomyces cerevisiae HI genes are preferentially retained across the kingdom Fungi, and that the HI genes of S. cerevisiae can be used to predict haploinsufficiency in humans. Our HI gene predictions confirm known associations between haploinsufficiency and genetic disease, and predict several further disorders in which the phenotype may be relevant. Haploinsufficiency is also clearly relevant to the gene-dosage imbalances inherent in eukaryotic sex-determination systems. In S. cerevisiae, HI genes are over-represented on chromosome III, the chromosome that determines yeast's mating type. This may be a device to select against the loss of one copy of chromosome III from a diploid. We found that orthologues of S. cerevisiae HI genes are also over-represented on the mating-type chromosomes of other yeasts and filamentous fungi. In animals with heterogametic sex determination, accumulation of HI genes on the sex chromosomes would compromise fitness in both sexes, given X chromosome inactivation in females. We found that orthologues of S. cerevisiae HI genes are significantly under-represented on the X chromosomes of mammals and of Caenorhabditis elegans. There is no X inactivation in Drosophila melanogaster (increased expression of X in the male is used instead) and, in this species, we found no depletion of orthologues to yeast HI genes on the sex chromosomes.

CONCLUSION

A special relationship between HI genes and the sex/mating-type chromosome extends from S. cerevisiae to Homo sapiens, with the microbe being a useful model for species throughout the evolutionary range. Furthermore, haploinsufficiency in yeast can predict the phenotype in higher organisms.

Ancient gene duplication provided a key molecular step for anaerobic growth of Baker's yeast

Mitochondria are essential organelles required for a number of key cellular processes. As most mitochondrial proteins are nuclear encoded, their efficient translocation into the organelle is critical. Transport of proteins across the inner membrane is driven by a multicomponent, matrix-localized "import motor," which is based on the activity of the molecular chaperone Hsp70 and a J-protein cochaperone. In Saccharomyces cerevisiae, two paralogous J-proteins, Pam18 and Mdj2, can form the import motor. Both contain transmembrane and matrix domains, with Pam18 having an additional intermembrane space (IMS) domain. Evolutionary analyses revealed that the origin of the IMS domain of S. cerevisiae Pam18 coincides with a gene duplication event that generated the PAM18/MDJ2 gene pair. The duplication event and origin of the Pam18 IMS domain occurred at the relatively ancient divergence of the fungal subphylum Saccharomycotina. The timing of the duplication event also corresponds with a number of additional functional changes related to mitochondrial function and respiration. Physiological and genetic studies revealed that the IMS domain of Pam18 is required for efficient growth under anaerobic conditions, even though it is dispensable when oxygen is present. Thus, the gene duplication was beneficial for growth capacity under particular environmental conditions as well as diversification of the import motor components.