Large-scale genomic analyses of in vitro yeast-mycelium dimorphism in human, insect and plant pathogenic fungi: From ESTs to RNAseq experiments

Unraveling the transcriptional features and gene expression networks of pathogenic and saprotrophic Ophiostoma species during the infection of Ulmus americana

American elm (Ulmus americana), highly prized for its ornamental value, has suffered two successive outbreaks of Dutch elm disease (DED) caused by ascomycete fungi belonging to the genus Ophiostoma. To identify the genes linked to the pathogenicity of different species and lineages of Ophiostoma, we inoculated 2-year-old U. americana saplings with six strains representing three species of DED fungi, and one strain of the saprotroph Ophiostoma quercus. Differential expression analyses were performed following RNA sequencing of fungal transcripts recovered at 3- and 10-days post-infection. Based on a total of 8,640 Ophiostoma genes, we observed a difference in fungal gene expression depending on the strain inoculated and the time of incubation in host tissue. Some genes overexpressed in the more virulent strains of Ophiostoma encode hydrolases that possibly act synergistically. A mutant of Ophiostoma novo-ulmi in which the gene encoding the ogf1 transcription factor had been deleted did not produce transcripts for the gene encoding the hydrophobin cerato-ulmin and was less virulent. Weighted gene correlation network analyses identified several candidate pathogenicity genes distributed among 13 modules of interconnected genes.IMPORTANCEOphiostoma is a genus of cosmopolitan fungi that belongs to the family Ophiostomataceae and includes the pathogens responsible for two devastating pandemics of Dutch elm disease (DED). As the mechanisms of action of DED agents remain unclear, we carried out the first comparative transcriptomic study including representative strains of the three Ophiostoma species causing DED, along with the phylogenetically close saprotrophic species Ophiostoma quercus. Statistical analyses of the fungal transcriptomes recovered at 3 and 10 days following infection of Ulmus americana saplings highlighted several candidate genes associated with virulence and host-pathogen interactions wherein each strain showed a distinct transcriptome. The results of this research underscore the importance of investigating the transcriptional behavior of different fungal taxa to understand their pathogenicity and virulence in relation to the timeline of infection.

Comparative transcriptomics reveal different mechanisms for hyphal growth across four plant-associated dimorphic fungi

Fungal dimorphism is a phenomenon by which a fungus can grow both as a yeast form and a hyphal form. It is frequently related to pathogenicity as different growth forms are more suitable for different functions during a life cycle. Among dimorphic plant pathogens, the corn smut fungus Ustilago maydis serves as a model organism to understand fungal dimorphism and its effect on pathogenicity. However, there is a lack of data on whether mechanisms elucidated from model species are broadly applicable to other fungi. In this study, two non-model plant-associated species in the smut fungus subphylum (Ustilaginomycotina), Tilletiopsis washingtonensis and Meira miltonrushii, were selected to compare dimorphic mechanisms in these to those in U. maydis. We sequenced transcriptomic profiles during both yeast and hyphal growth in these two species using Tween40, a lipid mimic, as a trigger for hyphal growth. We then compared our data with previously published data from U. maydis and a fourth but unrelated dimorphic phytopathogen, Ophiostoma novo-ulmi. Comparative transcriptomics was performed to identify common genes upregulated during hyphal growth in all four dimorphic species. Intriguingly, T. washingtonensis shares the least similarities of transcriptomic alteration (hyphal growth versus yeast growth) with the others, although it is closely related to M. miltonrushii and U. maydis. This suggests that phylogenetic relatedness is not correlated with transcriptomic similarity under the same biological phenomenon. Among commonly expressed genes in the four species, genes in cell energy production and conversion, amino acid transport and metabolism and cytoskeleton are significantly enriched. Considering dimorphism genes characterized in U. maydis, as well as hyphal tip-associated genes from the literature, we found only genes encoding the cell end marker Tea4/TeaC and the kinesin motor protein Kin3 concordantly expressed in all four species. This suggests a divergence in species-specific mechanisms for dimorphic transition and hyphal growth.

Quantifying the dominant growth mechanisms of dimorphic yeast using a lattice-based model

A mathematical model is presented for the growth of yeast that incorporates both dimorphic behaviour and nutrient diffusion. The budding patterns observed in the standard and pseudohyphal growth modes are represented by a bias in the direction of cell proliferation. A set of spatial indices is developed to quantify the morphology and compare the relative importance of the directional bias to nutrient concentration and diffusivity on colony shape. It is found that there are three different growth modes: uniform growth, diffusion-limited growth (DLG) and an intermediate region in which the bias determines the morphology. The dimorphic transition due to nutrient limitation is investigated by relating the directional bias to the nutrient concentration, and this is shown to replicate the behaviour observed in vivo Comparisons are made with experimental data, from which it is found that the model captures many of the observed features. Both DLG and pseudohyphal growth are found to be capable of generating observed experimental morphologies.