Neurospora from Natural Populations: Population Genomics Insights into the Life History of a Model Microbial Eukaryote

Draft genome sequences of Neurospora crassa clade B, isolated from burned Cytisus sp. plants in France

Neurospora crassa clade A is a model system for genetics, biochemistry, molecular biology, and experimental evolution. Here, we present the draft genome sequences of four isolates of N. crassa clade B. These data represent a valuable resource to investigate the population biology and evolutionary history of N. crassa sensu lato.

Genome-scale phylogeny and comparative genomics of the fungal order Sordariales

The order Sordariales is taxonomically diverse, and harbours many species with different lifestyles and large economic importance. Despite its importance, a robust genome-scale phylogeny, and associated comparative genomic analysis of the order is lacking. In this study, we examined whole-genome data from 99 Sordariales, including 52 newly sequenced genomes, and seven outgroup taxa. We inferred a comprehensive phylogeny that resolved several contentious relationships amongst families in the order, and cleared-up intrafamily relationships within the Podosporaceae. Extensive comparative genomics showed that genomes from the three largest families in the dataset (Chaetomiaceae, Podosporaceae and Sordariaceae) differ greatly in GC content, genome size, gene number, repeat percentage, evolutionary rate, and genome content affected by repeat-induced point mutations (RIP). All genomic traits showed phylogenetic signal, and ancestral state reconstruction revealed that the variation of the properties stems primarily from within-family evolution. Together, the results provide a thorough framework for understanding genome evolution in this important group of fungi.

Lineage-specific genes are clustered with HET-domain genes and respond to environmental and genetic manipulations regulating reproduction in Neurospora

Lineage-specific genes (LSGs) have long been postulated to play roles in the establishment of genetic barriers to intercrossing and speciation. In the genome of Neurospora crassa, most of the 670 Neurospora LSGs that are aggregated adjacent to the telomeres are clustered with 61% of the HET-domain genes, some of which regulate self-recognition and define vegetative incompatibility groups. In contrast, the LSG-encoding proteins possess few to no domains that would help to identify potential functional roles. Possible functional roles of LSGs were further assessed by performing transcriptomic profiling in genetic mutants and in response to environmental alterations, as well as examining gene knockouts for phenotypes. Among the 342 LSGs that are dynamically expressed during both asexual and sexual phases, 64% were detectable on unusual carbon sources such as furfural, a wildfire-produced chemical that is a strong inducer of sexual development, and the structurally-related furan 5-hydroxymethyl furfural (HMF). Expression of a significant portion of the LSGs was sensitive to light and temperature, factors that also regulate the switch from asexual to sexual reproduction. Furthermore, expression of the LSGs was significantly affected in the knockouts of adv-1 and pp-1 that regulate hyphal communication, and expression of more than one quarter of the LSGs was affected by perturbation of the mating locus. These observations encouraged further investigation of the roles of clustered lineage-specific and HET-domain genes in ecology and reproduction regulation in Neurospora, especially the regulation of the switch from the asexual growth to sexual reproduction, in response to dramatic environmental conditions changes.

Origins of lineage-specific elements via gene duplication, relocation, and regional rearrangement in Neurospora crassa

The origin of new genes has long been a central interest of evolutionary biologists. However, their novelty means that they evade reconstruction by the classical tools of evolutionary modelling. This evasion of deep ancestral investigation necessitates intensive study of model species within well-sampled, recently diversified, clades. One such clade is the model genus Neurospora, members of which lack recent gene duplications. Several Neurospora species are comprehensively characterized organisms apt for studying the evolution of lineage-specific genes (LSGs). Using gene synteny, we documented that 78% of Neurospora LSG clusters are located adjacent to the telomeres featuring extensive tracts of non-coding DNA and duplicated genes. Here, we report several instances of LSGs that are likely from regional rearrangements and potentially from gene rebirth. To broadly investigate the functions of LSGs, we assembled transcriptomics data from 68 experimental data points and identified co-regulatory modules using Weighted Gene Correlation Network Analysis, revealing that LSGs are widely but peripherally involved in known regulatory machinery for diverse functions. The ancestral status of the LSG mas-1, a gene with roles in cell-wall integrity and cellular sensitivity to antifungal toxins, was investigated in detail alongside its genomic neighbours, indicating that it arose from an ancient lysophospholipase precursor that is ubiquitous in lineages of the Sordariomycetes. Our discoveries illuminate a "rummage region" in the N. crassa genome that enables the formation of new genes and functions to arise via gene duplication and relocation, followed by fast mutation and recombination facilitated by sequence repeats and unconstrained non-coding sequences.

Genetic structure and triazole resistance among Aspergillus fumigatus populations from remote and undeveloped regions in Eastern Himalaya

Aspergillus fumigatus is a ubiquitous mold and a common human fungal pathogen. Recent molecular population genetic and epidemiological analyses have revealed evidence for long-distance gene flow and high genetic diversity within most local populations of A. fumigatus. However, little is known about the impact of regional landscape factors in shaping the population diversity patterns of this species. Here we sampled extensively and investigated the population structure of A. fumigatus from soils in the Three Parallel Rivers (TPR) region in Eastern Himalaya. This region is remote, undeveloped and sparsely populated, bordered by glaciated peaks more than 6,000 m above sea level, and contained three rivers separated by tall mountains over very short horizontal distances. A total of 358 A. fumigatus strains from 19 sites along the three rivers were isolated and analyzed at nine loci containing short tandem repeats. Our analyses revealed that mountain barriers, elevation differences, and drainage systems all contributed low but statistically significant genetic variations to the total A. fumigatus population in this region. We found abundant novel alleles and genotypes in the TPR population of A. fumigatus and significant genetic differentiation between this population and those from other parts of Yunnan and the globe. Surprisingly, despite limited human presence in this region, about 7% of the A. fumigatus isolates were resistant to at least one of the two medical triazoles commonly used for treating aspergillosis. Our results call for greater surveillance of this and other human fungal pathogens in the environment. IMPORTANCE The extreme habitat fragmentation and substantial environmental heterogeneity in the TPR region have long known to contribute to geographically shaped genetic structure and local adaptation in several plant and animal species. However, there have been limited studies of fungi in this region. Aspergillus fumigatus is a ubiquitous pathogen capable of long-distance dispersal and growth in diverse environments. In this study, using A. fumigatus as a model, we investigated how localized landscape features contribute to genetic variations in fungal populations. Our results revealed that elevation and drainage isolation rather than direct physical distances significantly impacted genetic exchange and diversity among the local A. fumigatus populations. Interestingly, within each local population, we found high allelic and genotypic diversities, and with evidence ~7% of all isolates being resistant to two medical triazoles, itraconazole and voriconazole. Given the high frequency of ARAF found in mostly natural soils of sparsely populated sites in the TPR region, close monitoring of their dynamics in nature and their effects on human health is needed.

Clonality, inbreeding, and hybridization in two extremotolerant black yeasts

BACKGROUND

The great diversity of lifestyles and survival strategies observed in fungi is reflected in the many ways in which they reproduce and recombine. Although a complete absence of recombination is rare, it has been reported for some species, among them 2 extremotolerant black yeasts from Dothideomycetes: Hortaea werneckii and Aureobasidium melanogenum. Therefore, the presence of diploid strains in these species cannot be explained as the product of conventional sexual reproduction.

RESULTS

Genome sequencing revealed that the ratio of diploid to haploid strains in both H. werneckii and A. melanogenum is about 2:1. Linkage disequilibrium between pairs of polymorphic loci and a high degree of concordance between the phylogenies of different genomic regions confirmed that both species are clonal. Heterozygosity of diploid strains is high, with several hybridizing genome pairs reaching the intergenomic distances typically seen between different fungal species. The origin of diploid strains collected worldwide can be traced to a handful of hybridization events that produced diploids, which were stable over long periods of time and distributed over large geographic areas.

CONCLUSIONS

Our results, based on the genomes of over 100 strains of 2 black yeasts, show that although they are clonal, they occasionally form stable and highly heterozygous diploid intraspecific hybrids. The mechanism of these apparently rare hybridization events, which are not followed by meiosis or haploidization, remains unknown. Both extremotolerant yeasts, H. werneckii and even more so A. melanogenum, a close relative of the intensely recombining and biotechnologically relevant Aureobasidium pullulans, provide an attractive model for studying the role of clonality and ploidy in extremotolerant fungi.

On the Mechanistic Basis of Killer Meiotic Drive in Fungi

Spore killers are specific genetic elements in fungi that kill sexual spores that do not contain them. A range of studies in the last few years have provided the long-awaited first insights into the molecular mechanistic aspects of spore killing in different fungal models, including both yeast-forming and filamentous Ascomycota. Here we describe these recent advances, focusing on the wtf system in the fission yeast Schizosaccharomyces pombe; the Sk spore killers of Neurospora species; and two spore-killer systems in Podospora anserina, Spok and [Het-s]. The spore killers appear thus far mechanistically unrelated. They can involve large genomic rearrangements but most often rely on the action of just a single gene. Data gathered so far show that the protein domains involved in the killing and resistance processes differ among the systems and are not homologous. The emerging picture sketched by these studies is thus one of great mechanistic and evolutionary diversity of elements that cheat during meiosis and are thereby preferentially inherited over sexual generations.

From Glaciers to Refrigerators: the Population Genomics and Biocontrol Potential of the Black Yeast Aureobasidium subglaciale

Apples are affected by numerous fungi known as storage rots, which cause significant losses before and after harvest. Concerns about increasing antimicrobial resistance, bans on various fungicides, and changing consumer preferences are motivating the search for safer means to prevent fruit rot. The use of antagonistic microbes has been shown to be an efficient and environmentally friendly alternative to conventional phytopharmaceuticals. Here, we investigate the potential of Aureobasidium subglaciale for postharvest rot control. We tested the antagonistic activity of 9 strains of A. subglaciale and 7 closely related strains against relevant phytopathogenic fungi under conditions simulating low-temperature storage: Botrytis cinerea, Penicillium expansum, and Colletotrichum acutatum. We also investigated a selection of phenotypic traits of all strains and sequenced their whole genomes. The tested strains significantly reduced postharvest rot of apples at low temperatures caused by B. cinerea, C. acutatum (over 60%), and P. expansum (about 40%). Several phenotypic traits were observed that may contribute to this biocontrol capacity: growth at low temperatures, tolerance to high temperatures and elevated solute concentrations, and strong production of several extracellular enzymes and siderophores. Population genomics revealed that 7 of the 15 strains originally identified as A. subglaciale most likely belong to other, possibly undescribed species of the same genus. In addition, the population structure and linkage disequilibrium of the species suggest that A. subglaciale is strictly clonal and therefore particularly well suited for use in biocontrol. Overall, these data suggest substantial biological control potential for A. subglaciale, which represents another promising biological agent for disease control in fresh fruit.

IMPORTANCE After harvest, fruits are often stored at low temperatures to prolong their life. However, despite the low temperatures, much of the fruit is lost to rot caused by a variety of fungi, resulting in major economic losses and food safety risks. An increasingly important environmentally friendly alternative to conventional methods of mitigating the effects of plant diseases is the use of microorganisms that act similarly to probiotics—occupying the available space, producing antimicrobial compounds, and consuming the nutrients needed by the rot-causing species. To find a new microorganism for biological control that is particularly suitable for cold storage of fruit, we tested different isolates of the cold-loving yeast Aureobasidium subglaciale and studied their phenotypic characteristics and genomes. We demonstrated that A. subglaciale can significantly reduce rotting of apples caused by three rot-causing molds at low temperatures and thus has great potential for preventing fruit rot during cold storage.

Secondary Metabolism Gene Clusters Exhibit Increasingly Dynamic and Differential Expression during Asexual Growth, Conidiation, and Sexual Development in Neurospora crassa

Secondary metabolite clusters (SMCs) encode the machinery for fungal toxin production. However, understanding their function and analyzing their products requires investigation of the developmental and environmental conditions in which they are expressed. Gene expression is often restricted to specific and unexamined stages of the life cycle. Therefore, we applied comparative genomics analyses to identify SMCs in Neurospora crassa and analyzed extensive transcriptomic data spanning nine independent experiments from diverse developmental and environmental conditions to reveal their life cycle-specific gene expression patterns. We reported 20 SMCs comprising 177 genes-a manageable set for investigation of the roles of SMCs across the life cycle of the fungal model N. crassa-as well as gene sets coordinately expressed in 18 predicted SMCs during asexual and sexual growth under three nutritional and two temperature conditions. Divergent activity of SMCs between asexual and sexual development was reported. Of 126 SMC genes that we examined for knockout phenotypes, al-2 and al-3 exhibited phenotypes in asexual growth and conidiation, whereas os-5, poi-2, and pmd-1 exhibited phenotypes in sexual development. SMCs with annotated function in mating and crossing were actively regulated during the switch between asexual and sexual growth. Our discoveries call for attention to roles that SMCs may play in the regulatory switches controlling mode of development, as well as the ecological associations of those developmental stages that may influence expression of SMCs. IMPORTANCE Secondary metabolites (SMs) are low-molecular-weight compounds that often mediate interactions between fungi and their environments. Fungi enriched with SMs are of significant research interest to agriculture and medicine, especially from the aspects of pathogen ecology and environmental epidemiology. However, SM clusters (SMCs) that have been predicted by comparative genomics alone have typically been poorly defined and insufficiently functionally annotated. Therefore, we have investigated coordinate expression in SMCs in the model system N. crassa, and our results suggest that SMCs respond to environmental signals and to stress that are associated with development. This study examined SMC regulation at the level of RNA to integrate observations and knowledge of these genes in various growth and development conditions, supporting combining comparative genomics and inclusive transcriptomics to improve computational annotation of SMCs. Our findings call for detailed study of the function of SMCs during the asexual-sexual switch, a key, often-overlooked developmental stage.

Synthetic Soil Aggregates: Bioprinted Habitats for High-Throughput Microbial Metaphenomics

The dynamics of microbial processes are difficult to study in natural soil, owing to the small spatial scales on which microorganisms operate and to the opacity and chemical complexity of the soil habitat. To circumvent these challenges, we have created a 3D-bioprinted habitat that mimics aspects of natural soil aggregates while providing a chemically defined and translucent alternative culturing method for soil microorganisms. Our Synthetic Soil Aggregates (SSAs) retain the porosity, permeability, and patchy resource distribution of natural soil aggregates—parameters that are expected to influence emergent microbial community interactions. We demonstrate the printability and viability of several different microorganisms within SSAs and show how the SSAs can be integrated into a multi-omics workflow for single SSA resolution genomics, metabolomics, proteomics, lipidomics, and biogeochemical assays. We study the impact of the structured habitat on the distribution of a model co-culture microbial community and find that it is significantly different from the spatial organization of the same community in liquid culture, indicating a potential for SSAs to reproduce naturally occurring emergent community phenotypes. The SSAs have the potential as a tool to help researchers quantify microbial scale processes in situ and achieve high-resolution data from the interplay between environmental properties and microbial ecology.

Virulence Traits and Population Genomics of the Black Yeast Aureobasidium melanogenum

The black yeast-like fungus Aureobasidium melanogenum is an opportunistic human pathogen frequently found indoors. Its traits, potentially linked to pathogenesis, have never been systematically studied. Here, we examine 49 A. melanogenum strains for growth at 37 °C, siderophore production, hemolytic activity, and assimilation of hydrocarbons and human neurotransmitters and report within-species variability. All but one strain grew at 37 °C. All strains produced siderophores and showed some hemolytic activity. The largest differences between strains were observed in the assimilation of hydrocarbons and human neurotransmitters. We show for the first time that fungi from the order Dothideales can assimilate aromatic hydrocarbons. To explain the background, we sequenced the genomes of all 49 strains and identified genes putatively involved in siderophore production and hemolysis. Genomic analysis revealed a fairly structured population of A.melanogenum, raising the possibility that some phylogenetic lineages have higher virulence potential than others. Population genomics indicated that the species is strictly clonal, although more than half of the genomes were diploid. The existence of relatively heterozygous diploids in an otherwise clonal species is described for only the second time in fungi. The genomic and phenotypic data from this study should help to resolve the non-trivial taxonomy of the genus Aureobasidium and reduce the medical hazards of exploiting the biotechnological potential of other, non-pathogenic species of this genus.

Revisiting the Neurospora crassa mitochondrial genome

The mitochondrial genome of Neurospora crassa has been less studied than its nuclear counterpart, yet it holds great potential for understanding the diversity and evolution of this important fungus. Here we describe a new mitochondrial DNA (mtDNA) complete sequence of a N. crassa wild type strain. The genome with 64 839 bp revealed 21 protein-coding genes and several hypothetical open reading frames with no significant homology to any described gene. Five large repetitive regions were identified across the genome, including partial or complete genes. The largest repeated region holds a partial nd2 section that was also detected in Neurospora intermedia, suggesting a rearrangement that occurred before the N. crassa speciation. Interestingly, N. crassa has a palindrome adjacent to the partial nd2 repeated region possibly related to the genomic rearrangement, which is absent in N. intermedia. Finally, we compared the sequences of the three available N. crassa complete mtDNAs and found low levels of intraspecific variability. Most differences among strains were due to small indels in noncoding regions. The revisiting of the N. crassa mtDNA forms the basis for future studies on mitochondrial genome organization and variability.